准确模拟纳米螺旋结构对于推动新型纳米技术应用的预测性仿真研究至关重要。 目前,软件包和代码在构建原子级螺旋模型方面存在局限性。 我们提出了两种用于构建原子级纳米螺旋模型的模拟方法,并提供了一个图形界面,以通过可视化手段增强研究能力。
方法文章
准确模拟纳米螺旋结构对于推动新型纳米技术应用的预测性仿真研究至关重要。 目前,软件包和代码在构建原子级螺旋模型方面存在局限性。 我们提出了两种用于构建原子级纳米螺旋模型的模拟方法,并提供了一个图形界面,以通过可视化手段增强研究能力。
弹簧状材料在自然界中普遍存在,并在纳米技术中具有重要应用价值,可用于能量收集、氢气存储和生物传感等领域。 为了实现可预测的模拟,准确建模纳米螺旋结构变得日益重要。 为了研究局部结构对这些复杂几何形貌性能的影响,必须建立真实可靠的模型。 迄今为止,现有的软件包在构建原子级螺旋模型方面仍相当有限。 本研究聚焦于为分子动力学(MD)模拟生成原子级的二氧化硅玻璃(SiO2)纳米带和纳米弹簧模型。基于“体相”二氧化硅玻璃的MD模型,本文提出了两种计算方法,用于精确构建纳米带和纳米弹簧的形状。 第一种方法采用AWK编程语言和开源软件,通过设定目标尺寸和螺旋参数方程,从初始体相模型中高效地切割出各种形状的二氧化硅纳米带。 利用该方法,可针对不同螺距和尺寸参数生成精确的原子级二氧化硅纳米带。 第二种方法涉及一个更为强大的代码系统,可在建模纳米螺旋结构时提供更高的灵活性。 该方法采用专门编写的C++代码,集成了预筛选算法以及螺旋的数学方程,从而在构建纳米弹簧模型时实现更高的精度和效率。 利用这些代码,可高效地生成结构明确且可扩展的纳米带和纳米弹簧模型,适用于原子级模拟。 这两个开源代码的附加价值在于,它们可被修改以复现不同类型的螺旋结构,且不受材料种类限制。 此外,还使用MATLAB图形用户界面(GUI)增强可视化与交互功能,帮助普通用户更好地理解原子级螺旋结构,促进学习与应用。 这些方法的一个应用实例是近期通过分子动力学模拟研究纳米螺旋结构在机械能量收集方面的性能。
螺旋纳米结构通常在实验室中采用化学气相沉积技术制备1-2,同时文献中也报道了新的制备方法3。特别是纳米弹簧和纳米带因其独特的性质以及在传感器、光学、机电和流体器件中的广阔应用前景而受到广泛研究4-7。已有报道的合成方法可制备二氧化硅(SiO2)纳米带,使这些结构成为构建多级体系的潜在基本单元。新型三维二氧化硅纳米弹簧的合成进一步拓展了其应用范围,例如涂覆氧化锌(ZnO)8后可用于化学电阻型传感器,或结合纳米颗粒用于诊断应用9-10。
关于二氧化硅纳米弹簧和纳米带力学性能的实验研究较为稀少,这主要归因于当前在操控和测试方法及设备方面的局限性。已有文献通过理论和模拟方法报道了对纳米结构和纳米弹簧纳米力学的研究11-14。一些模拟研究13聚焦于非晶态纳米弹簧的纳米力学行为,因为模拟能够探索实验手段尚难以完全触及的物理机制。文献中已有关于金属纳米弹簧的原子级研究,用于探究弹性性能的尺寸依赖性15,以及近期对螺旋状晶态二氧化硅纳米结构纳米力学的研究14。 针对纳米弹簧结构的实验测试也在其他材料体系中开展,例如螺旋碳纳米结构和碳纳米线圈16-17。尽管目前已积累了部分知识,但为了推动未来纳米器件的制备工作,仍需对这些新型纳米结构的力学性能获得更全面的理解。
由于关于二氧化硅玻璃(非晶态二氧化硅)纳米螺旋的分子动力学(MD)研究仍然十分有限,此类结构的原子尺度建模需要开发定制化的程序代码。根据近期文献检索,目前尚未发现其他可用于构建二氧化硅玻璃螺旋MD模型的方法。本研究采用自下而上的方法,对包括纳米弹簧和纳米带在内的螺旋状二氧化硅玻璃纳米结构进行原子级建模,以用于未来的大型分子动力学纳米力学模拟。该方法首先构建一个此前已报道过的MD“体相”二氧化硅玻璃模型18,然后通过为此专门开发的两个高效且可调的计算机程序,从该“体相”样品中切割出各种螺旋纳米结构。这两种计算流程各自提供了一种独特的方式,能够高效且精细地构建纳米带和纳米弹簧模型,所得结构适用于大规模原子级模拟。此外,还使用了一个定制的图形用户界面,以辅助螺旋结构的构建与可视化。
“块体”二氧化硅玻璃模型的结构最初在室温下构建。为此目的,采用与先前研究类似的Garofalini原子间势函数进行大规模分子动力学(MD)模拟18,该势函数在计算上相对高效,适用于大尺度体系。初始的“块体”二氧化硅玻璃结构包含一个立方模型(14.3 × 14.3 × 14.3 nm3),共含有192,000个原子。该“块体”二氧化硅玻璃模型在300 K下使用周期性边界条件弛豫0.5纳秒,以获得初始平衡状态。
设计并采用了两种计算方法来构建原子级二氧化硅纳米带和纳米弹簧模型。第一种方法涉及使用定义螺旋线的参数方程及其几何参数(螺距、螺旋半径和线丝半径),从“体相”结构中切割出二氧化硅纳米带。该过程包括使用AWK编程语言、LINUX操作系统以及开源可视化软件19。构建纳米带原子模型的一般迭代步骤包括:(1)在“体相”二氧化硅玻璃模型中选择一个原子;(2)计算所选原子与预定义螺旋函数上某空间点之间的距离;(3)将该距离与目标纳米带的半径进行比较;(4)决定是否在输出数据模型中保留或舍弃该原子。该方法的详细逐步说明包含在可扩展开源代码补充材料中。利用此方法,通过不同的螺距、螺旋半径和纳米带半径值构建了多个二氧化硅纳米带,并随后使用分子分析与可视化软件19-20将其实际尺寸与目标尺寸进行比对以验证准确性。生成了具有功能性几何特征(高螺距值和低纳米带半径值)的二氧化硅纳米带原子模型。在纳米带半径值极高而螺距值极低的情况下,观察到一些由原子误删导致的伪影,表现为纳米带表面不够平滑。类似的方法此前已被用于构建二氧化硅纳米线21-23。
本文介绍的第二种方法包括通过预筛选方法提高效率,并结合螺旋线的数学方程,从“整体”二氧化硅结构中切割出二氧化硅纳米弹簧。该方法需要编写更强大的 C++ 代码,以在建模这些螺旋纳米结构时提供更大的灵活性。构建纳米弹簧原子级模型的迭代方法包括:(1)剔除所有确定位于螺旋路径之外的原子,(2)确定性地选择螺旋路径上的一个点,(3)比较该选定点特定距离范围内的所有原子,以及(4)将每个原子判断为舍弃或存入输出数据模型中。该方法的逐步说明也包含在可扩展开源代码补充材料中。利用该方法,获得了多个具有不同尺寸(线半径、螺旋半径和纳米弹簧螺距)的二氧化硅纳米弹簧模型,如图1所示。该方法高效地获得了高精度的二氧化硅纳米弹簧模型,在纳米弹簧的极端(低值和高值)螺距条件下均未发现伪影。该方法的图形用户界面的创建与使用在实验方案部分中进行了描述。

图 1: 通用螺旋结构示意图,显示其特征尺寸,其中 r、R 和 p 分别表示导线半径、螺旋半径和螺距。H 表示螺旋结构的总高度23。
本方案介绍了如何准备NanospringCarver文件,在LINUX24 PC上运行MATLAB25,并使用图形用户界面构建原子级纳米弹簧模型。这些此前无法获得的模型可作为新型分子动力学(MD)模拟23的基础,推动材料创新研究。
构建原子级纳米弹簧模型的通用逐步操作流程包括使用以下要素:(a)NanospringCarver(0.5 测试版)代码(以 C++ 语言编写的开源程序),(b)体相二氧化硅玻璃模型(输入文件),(c)MATLAB 图形用户界面(GUI)及相关文件,以及(d)在 LINUX 计算机上使用本地许可证的 MATLAB 软件(版本 7)。上述项目(a)至(c)(即 NanospringCarver 代码、二氧化硅玻璃模型、MATLAB GUI 文件)均可免费在线下载26。MATLAB(Matrix Laboratory)是由 MathWorks 公司开发的一种用于数值计算、可视化和应用程序开发的高级语言24,广泛应用于数据可视化与分析、图像处理以及计算生物学领域。
访问受限。请登录或开始试用以查看此内容。
1. 准备NanospringCarver文件并在LINUX电脑上启动MATLAB
以下步骤旨在帮助普通用户利用在线提供的文件26。
2. 修改和使用 NanospringCarver 程序的图形用户界面(GUI)
使用在线提供的文件,按照以下步骤操作26。

图 2:MATLAB 用户界面,展示如何打开 MATLAB GUIDE。

图3:MATLAB GUIDE 接口初始化。

图 4: MATLAB GUIDE 界面,展示如何打开现有的 GUI 图形文件。

图 5: MATLAB GUIDE 界面,显示用于修改现有 GUI 图形的工具。

图6:使用图形用户界面创建二氧化硅纳米弹簧模型示例的截图。

图 7: 基于图形用户界面运行纳米弹簧模拟时 MATLAB 命令窗口的反馈信息。
注意:在上述示例中,文件“model”包含组成目标弹簧的 5,176 个原子,每个原子占一行,第一行为文件中原子的总数。每一行定义一个原子,包括原子 ID、原子类型以及该原子的 x、y、z 坐标。
3. 在开源可视化工具中验证 NanospringCarver 结果19
以下步骤旨在帮助普通用户可视化并验证由 NanospringCarver 生成的弹簧模型输出。
4. 使用 NanospringCarver 进行纳米弹簧的分子动力学拉伸模拟
以下步骤概述了普通用户如何将 NanospringCarver 创建的弹簧模型作为输入,用于传统的开源分子动力学(MD)代码32。

图8: 二氧化硅纳米弹簧在拉伸模拟过程中的截图(另见动画图1)。
访问受限。请登录或开始试用以查看此内容。
通过第一种计算程序(纳米带代码)构建的原子级纳米带模型及其相应尺寸如图9所示。利用第二种计算程序(纳米弹簧代码)生成的纳米弹簧模型及其相应尺寸如图10所示。

图9. 具有所需尺寸的二氧化硅纳米带的原子级模型:r(纳米带半径)= 1.07 nm,R(螺旋半径)= 5.37 nm,以及 p(螺距)= 7.16 nm。快照展示了该纳米结构的不同视图:(a) 顶视图,(b) 侧视图,(c) 附加旋转后的侧视图,以及 (d)-(...
访问受限。请登录或开始试用以查看此内容。
对原始方法进行改进以构建纳米螺旋结构,由此开发出两种不同的代码,可从初始的体相二氧化硅玻璃分子动力学模型出发,生成纳米带和纳米弹簧结构。 通过使用不同的软件包19-20对二氧化硅纳米带和纳米弹簧模型进行了验证,确认其尺寸精度在程序的测量能力范围内。 还通过从不同侧面和角度叠加模型,对纳米弹簧与纳米带进行了比较,进一步实现了几何结构的验证。 本项目开发的两种计算方法以不同的方式生成螺旋纳米结构,其附加价值在于可扩展用于任意尺寸的体相材料模型,并具备用于模拟其他材料来源的纳米螺旋结构的潜力。 本文展示的结果表明,使用任一方法均未产生可检测到的伪影(即目标纳米螺旋结构中缺失的原子)。 此外,本研究开发的计算方法在生成右手螺旋或左手螺旋纳米结构方面具有灵活性,仅需反转定义螺旋的正弦和余弦函数的顺序即可实现。 该方法未来的应用将包括扩展至更大的螺旋结构,以实现更广泛的参数变化,并探索其在不同初始材料中的应用。
该方法的局限性包括所构建的纳米螺旋结构的尺寸受限于所采用的初始块体二氧化硅模型,且随着模型...
访问受限。请登录或开始试用以查看此内容。
作者声明不存在任何竞争性经济利益。
作者感谢加州大学默塞德分校的 Tim Allis 在本项目中的协助。 加州大学默塞德分校的 NSF-COINS 项目在本研究早期阶段为 (KAM) 提供了支持。 一项 NSF-BRIGE 奖项支持了合著者 (BND 和 KAM),为本研究以及参加学术会议的差旅费用提供了资金。
研究团队首要感谢美国国家科学基金会通过BRIGE奖项对本工作的资助。 本材料基于美国国家科学基金会在资助号1032653下支持的工作。
访问受限。请登录或开始试用以查看此内容。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| MATLAB 数值计算软件 | Mathworks | http://www.mathworks.com/products/matlab/description1.html | 参见实验方案引言部分及参考文献[24] |
| NanospringCarver 程序代码与文件 | 加州大学默塞德分校 - 开源 | http://tinyurl.com/qame8dj | 参见实验方案第1节(步骤1.2)及参考文献[26] |
| MATLAB GUI 文件 | 加州大学默塞德分校 - 开源 | http://tinyurl.com/qame8dj | 参见实验方案第1节(步骤1.2)及参考文献[26] |
| 原子级块体玻璃输入文件 | 加州大学默塞德分校 - 开源 | http://tinyurl.com/qame8dj | 参见实验方案第1节(步骤1.2)及参考文献[26] |
| IFrIT 可视化软件 | 开源软件 | http://sites.google.com/site/ifrithome/ | 参见实验方案第3节及参考文献[19] |
| LAMMPS 分子动力学软件 | 开源软件 | http://lammps.sandia.gov/ | 参见实验方案第4节及参考文献[32] |
访问受限。请登录或开始试用以查看此内容。
申请许可以重复使用本 JoVE 文章的文本或图表
申请许可A correction was made to Scalable Nanohelices for Predictive Studies and Enhanced 3D Visualization. Online article is missing its 29th author reference.
The reference section has added:
29. MathWorks. Use and create MATLAB MEX-files [Internet]. [Massachusetts]: MathWorks; [cited 2013 May 30]. Available from: http:// www.mathworks.com/help/matlab/call-mex-files-1.html (2014).