精华内容
下载资源
问答
  • AMBER分子动力学模拟软件(AMBER18, AMBER19) AMBER是著名的老牌开源分子动力学模拟软件,最新版本是AMBER19(但是19版本的并行模拟模块还没公布,只有AmberTools19)如果需要并行计算可以下载编译18版本 由于...
  • 随着计算机技术的深入发展和不同类型分子动力学模拟软件的开发运用,更多的煤科学科研人员从事煤大分子的动力学模拟计算领域研究,从煤的结构特征、孔径分布特性、热力学特性和基于孔径分布的吸附特性等方面,深入...
  • 基于Python的分子动力学模拟分析脚本| 该资料库包括一些必须在您的工具包中找到的实用脚本 要记住的好资源: :link: Langevin Dynamics: : :link: 联系人(官方实验室帐户): ://compbio.bioe.eng.marmara....
  • MS分子动力学模拟具体实施步骤,应用material studio进行分子动力模拟
  • 基于python的分子动力学模拟程序编写及应用.pdf
  • lammps分子动力学、gromacs分子动力学、与第一性原理vasp
  • 方法:采用分子动力学模拟与实验相结合的研究方法.结果:研究表明,色素溶液的热稳定性与其体系的势能值呈负线性相关,温度对其构象的影响可能通过调节侧链基团的分布来实现;当溶液pH为4.0~10.0时,体系势能值...
  • 采用分子动力学(MD)方法对激光诱导碎裂效应进行了数值模拟。通过分子动力学方法给出了样品受到激光辐照后的物理图像。其中激光能量为60 mJ,计算步长取3 fs。数值模拟结果表明分子动力学方法完全适用于激光诱导碎裂...
  • 抑癌小分子RITA配体与ASK1受体相互作用的分子对接和分子动力学模拟,吴清华,蒋泉,采用分子对接和分子动力学(MD)模拟方法研究了最新抑癌小分子RITA与ASK1蛋白的相互作用机制。先用分子对接程序将RITA对接到ASK1...
  • 分子动力学模拟lammps 粘度表面张力分析in文件分子动力学模拟lammps 粘度、表面张力分析in文件
  • 单晶Cu纳米线单向拉伸和压缩的分子动力学模拟,袁林,张江伟,加工硬化现象普遍存在于金属塑性加工过程中。长久以来,这种现象一直吸引着人们的注意力。当在低温(T<0.3Tm)条件下施加应力时,加工
  • 二维纳米板裂纹扩展的分子动力学模拟,龙述尧,童杰, 本文采用分子动力学研究了纳米板在拉伸的作用下裂纹扩展。采用Lennard-Jones势模型描述原子间作用,同时,基于二维密排矩形晶格常数计
  • 使用结构建模和分子动力学模拟,残基Leu84,Gly113,Asp116,Ala207和Leu286被认为是结构稳定性的关键残基。 然后将预测的残基Leu84突变为ALA。 结果表明,L84A突变可以提高壳聚糖酶EAG1的热稳定性。 与以前的研究...
  • 孔洞对纳米单晶Cu力学性能影响的分子动力学模拟,袁林,单德彬,当金属材料进入到纳米尺度,由于其极小的尺寸,含有的缺陷数量必然很少,单个或者几个孔洞显然会极大影响纳米材料的力学行为。鉴
  • 采用结合双温方程的分子动力学方法,数值模拟脉宽为200 fs,能量密度为30~45 mJ/cm2 的超快激光与CuZr非晶合金的相互作用过程。模拟结果表明:超快激光作用下CuZr 非晶材料中原子加热速度比普通晶态金属慢得多;作用...
  • Cu熔化及凝固过程的分子动力学模拟,黄维,梁工英,采用Embedded-Atom Method (EAM)作用势,利用分子动力学方法模拟Cu的熔化及凝固过程,研究了不同冷却速率对液态Cu凝固过程的影响,分析了升
  • 为了研究点突变(Met108→Leu108)对树胶醛糖结合蛋白(ABP)与配体结合能力的影响, 对ABP、 ABP结合树胶醛糖复合物及ABP结合半乳糖复合物以及它们各自的突变体分别进行60 ns的分子动力学模拟. 模拟结果表明, 108号残基...
  • 分子动力学轨迹的 Matlab 脚本。 当前功能 使用 AIMD_Toolkit,您可以 从 VASP 输出 XDATCAR 读取轨迹 使用 VESTA 可读的 CHG 格式创建核密度 引用 AIMD_Toolkit 如果您使用此代码,请引用以下论文 陆子恒等物理化学...
  • 一套命令行工具,用于分析Amber Molecular Dynamics模拟的化学相关特征。 可以检测两种物质之间的平面界面的位置,并针对任意类型,名称,电荷或残基的原子绘制相对于该界面的密度。 可以计算并绘制三个原子残基...
  • 扩散系数的分子动力学模拟,刘娟芳,,本文选用平衡分子动力学模拟方法,应用不同的势能模型对强极性分子水的扩散系数进行了模拟计算。结果表明,用分子动力学模拟方法
  • 包含liscense的 分子动力学模拟软件MS 8.0 ,很好使用,推荐给大家
  • 对晶层间不同含水量的蒙脱土-水-离子体系以及NaCl水溶液进行分子动力学模拟,并比较了矿物晶层间水分子及离子特性与水溶液中的差别。模拟结果表明:层间补偿阳离子受蒙脱矿物同晶置换位置影响,而层间水分子则与四面体...
  • 运用分子动力学模拟方法研究原子弛豫对于裂尖场的影响,发现弛豫后裂尖离散非线性区半径约为150?;如果原子区大于这一尺寸,连续介质力学的弹性场可以通过边界条件,有效影响原子尺度的裂尖行为。进一步研究表明:...
  • 先生 分子动力学模拟
  • qaa通过使用联合对角化(JADE)来分离信息来分析分子动力学(MD)轨迹。 JADE和QAA代码基于用Matlab编写的原始代码。 卡多佐(JF); Souloumiac,A。“非高斯信号的盲波束成形”。 IEE Proc F雷达信号过程1993,140...
  • 采用全原子分子动力学模拟研究了丙氨酸二肽分子在水溶液中的结构和氢键作用。径向分布函数和氢键统计用于分子动力学模拟的分析。从径向分布函数(RDF)发现不同类型的强氢键显示出不同的形成能力,而体系中的C-H…O弱...
  • 本发明属于纳米材料技术领域,涉及一种基于分子动力学模拟测试碳纳米管力学性能的方法。技术背景碳纳米管作为一维纳米材料,具有高强度,高柔韧性等优异的力学性能,其特殊的导电性能、良好的导热性能且重量轻,使其...

    本发明属于纳米材料技术领域,涉及一种基于分子动力学模拟测试碳纳米管力学性能的方法。

    技术背景

    碳纳米管作为一维纳米材料,具有高强度,高柔韧性等优异的力学性能,其特殊的导电性能、良好的导热性能且重量轻,使其在材料、电子、航天航空等领域得到广泛的应用。因此对碳纳米管相关性能的研究得到了越来越多的重视。

    碳纳米管属于纳米尺度材料,在实验方面,对其力学性能的研究易受到加载实验装置、实验环境和检测分析仪器的限制,并且实验消耗大,成本高。如公开号为CN 105300794A的发明专利公开了一种通过实验测纳米纤维的杨氏模量的方法和系统,其需要多种仪器设备相互配合,因此成本较高。

    分子动力学模拟作为一种“数值试验”的科学研究方法,由于其高效并且有效节约了实验研究成本,因此分子动力学模拟被广泛应用于材料、医药、机械等领域。

    关于分子动力学模拟的方法研究碳纳米管力学性能亦得到了不断的发展。但当前的研究仍存在一些问题:首先在微观尺度下基于维里公式计算应力的方法在相关学术界上一直存在争议,如文献The virial stress is not a measure of mechanical stress(Materials Research,2002,731:59-70)对于基于维里公式计算应力的方法提出了质疑。其次,由于早期的势能函数(如Tersoff势能、Brenner势能、COMPASS势能等)描述碳纳米管原子间相互作用精度相对较低,如公开号为CN 104535859A的发明专利公开了一种测试碳纳米管温度特性的方法,其力场选择COMPASS势能描述,这是一种通用型共价分子体系力场,对于碳原子体系精度稍差。此外,对于分子动力学模拟来说,系综、模拟参数、加载方式、统计方式等的选择很大程度上影响着模拟的准确性,因此正确地使用系综,采用有效加载和统计方式能够提高模拟的准确性。最后,还有一些关于分子动力学模拟的方法采用粗粒化模型,如公开号为CN 106934137A的发明专利公开了一种对石墨烯粗粒化的分子动力学模拟,由于采用粗粒化模型计算,其计算精度很难保证。

    技术实现要素:

    针对以上的问题和研究现状,本发明提供一种基于分子动力学模拟测试碳纳米管力学性能的方法,具体为采用一种适用于分子动力学模拟中的新型加载和统计模型受力的方式,避免了基于维里公式计算应力方法存在的争议性,并通过恰当选择原子间作用势能,设定合理的模拟参数和方法步骤,可以准确仿真计算碳纳米管(包括不同手性的碳纳米管、多壁碳纳米管、缺陷碳纳米管等)的力学性能参数,如应力应变关系、弹性模量及泊松比等。

    为了达到上述目的,本发明的技术方案如下:

    一种基于分子动力学模拟测试碳纳米管力学性能的方法,包括以下步骤:

    本发明通过Lammps中“Spring”(“弹簧”)命令对碳纳米管进行加载,并作为一种力的传感器,用于统计碳纳米管的受力。更具体地说,即在分子动力学模拟中使用“弹簧”作为加载工具,实现对碳纳米管加载;并根据牛顿第三定律,统计“弹簧”上的力即为碳纳米管所受的力,进而方便有效地得到碳纳米管所受的载荷。其中,本发明采用Airebo多体势能描述碳纳米管原子之间相互作用,它是对Tersoff作用势的改进,使其更适用于碳原子系统,因此数值模拟碳纳米管力学性能更加准确。

    所述的分子动力学模拟包含以下步骤:

    (1)由Materials Studio软件构建碳纳米管模型,通过Lammps工具包中msi2lmp程序将建立好的碳纳米管模型的坐标文件转换为Lammps可识别的数据文件。使用Lammps编写程序实现对碳纳米管力学相关性能的分子动力学模拟计算。所述的模型包括不同手性的碳纳米管、多壁碳纳米管以及缺陷碳纳米管。

    (2)设定模拟初始参数,使用三维模拟模式,计量单位选择metal格式,设定碳原子质量,设置模拟空间三个方向边界均为非周期边界。

    (3)定义原子间相互作用力,采用Airebo势能精确的描述C原子系统之间相互作用,即Pair_style类型设定为Airebo,在Lammps程序中将碳纳米管底端固定,顶端使用“弹簧”进行加载,具体为使用spring命令。设置0.25飞秒为模拟的时间步长,这是因为模拟中采用的是Verlet算法,时间步长的减小可以增加模拟的精确性。

    (4)为方便施加边界条件和数据统计,将碳纳米管模型整体分为3组,分别为上部分组,中部分组和下部分组。上部分组包括顶端分组、近顶端分组,下部分组包括近底端分组、底端分组。

    (5)对整体模型施加随机初速度,通过能量最小化的方式对模型初始构型进行优化。

    (6)优化结构后,重置时间步。使用recenter命令保证下部分组的初始位置仅在局部进行结构的调整,而不离开其位置。

    (7)对模型整体采用NVT系综进行驰豫,驰豫中的控温方法采用Nose-Hoover方法,将温度控制在0.01K左右,保证模拟不受温度影响。

    (8)将下部分组三个方向的速度和力置为零,约束碳纳米管下部分组的位移。

    (9)对中部分组采用NVT系综进行驰豫,并继续保持温度在0.01K,对上部分组采用NVE系综进行驰豫,保证系统总动量为零。

    (10)在上部分组通过“弹簧”对碳纳米管进行加载。更具体地说,此“弹簧”的初始位置为碳纳米管上部分组经过上述驰豫后的质心位置,然后通过对“弹簧”施加位移实现对碳纳米管的加载。其中,通过“弹簧”实现加载,采用一步一平衡的模拟机制,即加载一次,系统进行一次驰豫。统计碳纳米管各个分组的三个方向质心坐标,以及弹簧整体的力等,完成对碳纳米管力学相关性能的分子动力学模拟。

    (11)得到相应数据文件及模拟的图形文件等,根据需要对统计的信息等进行后处理,得到如弹性模量、泊松比等力学性能参数。

    本发明主要具有以下优点:

    a)本发明提供的数值模拟方法计算碳纳米管相关力学性能,不但可以有效减少因实验的带来成本与消耗,而且还可以模拟因相关设备和材料的限制等无法准确进行和完成的实验。

    b)通过本发明方法使用“弹簧”作为传感器来统计受力,避免了基于维里公式计算应力方法存在的争议性,而且可以快捷方便的计算应力。

    c)通过上述碳纳米管模拟方法的步骤及参数设置,仅通过改变上述步骤(4)的边界条件,即可实现对碳纳米管拉伸,压缩,弯曲等分子动力学模拟,计算弹性模量、泊松比等力学参数。

    d)此外,通过Materials Studio可以建立不同手性的碳纳米管、多壁碳纳米管以及缺陷碳纳米管模型,通过本发明分子动力学模拟方法亦可得到其力学性能相关参数。

    附图说明

    图1含碳帽的单壁碳纳米管拉伸模拟示意图;

    图2为含碳帽的单壁碳纳米管拉伸分子动力学模拟流程图;

    图中:1含碳帽的单壁碳纳米管;2“弹簧”;3顶端分组;4近顶端分组;5中部分组;6近底端分组;7底端分组;8坐标轴;9“弹簧”位移方向。

    具体实施例

    下面结合附图与实施例对本发明方法作进一步的详细说明:

    本发明具体实施例为对含碳帽单壁碳纳米管的拉伸分子动力学模拟,如图2为含碳帽的单壁碳纳米管拉伸分子动力学模拟流程图所示:

    1)含碳帽的单壁碳纳米管1模型由Materials Studio软件建立,通过Lammps工具包msi2lmp程序将建立好的含碳帽的单壁碳纳米管模型的坐标文件转换为Lammps可识别的数据文件。使用Lammps编写程序实现对含碳帽的单壁碳纳米管的拉伸分子动力学模拟,具体步骤如下:

    2)设定模拟初始参数,使用三维模拟模式,计量单位选择metal格式,设定碳原子质量,设置模拟空间三个方向边界均为非周期边界。

    3)定义原子间相互作用力,本发明采用Airebo势能,可以精确的描述C原子系统之间相互作用。设置0.25飞秒为模拟的时间步长,这是因为模拟中采用的是Verlet算法,时间步长的减小可以增加模拟的精确性。

    4)为了方便施加边界条件,对含碳帽的单壁碳纳米管1模型整体进行分组。具体分组如图1含碳帽的单壁碳纳米管拉伸模拟示意图所示:含碳帽的单壁碳纳米管1自上而下分为5组,分别为顶端3分组、近顶端4分组、中部5分组、近底端6分组、底端7分组。

    5)对整体模型施加随机初速度,并通过能量最小化的方式对模型初始构型进行优化。

    6)优化结构后,重置时间步。使用recenter命令保证底端与近底端分组的初始位置仅在局部进行结构的调整,而不离开其位置。

    7)对模型整体采用NVT系综进行驰豫,其中控温方法采用Nose-Hoover方法,将温度控制在0.01K左右,以此来保证模拟不受温度影响。

    8)通过将近底端、底端分组三个方向的速度和力置为零,从而约束了含碳帽的单壁碳纳米管近底端、底端分组的位移。

    9)然后对中部分组采用NVT系综进行驰豫,并继续保持温度在0.01K,对近顶端和顶端分组采用NVE系综进行驰豫,保证系统总动量为零。

    10)在近顶端分组通过“弹簧2”对含碳帽的单壁碳纳米管进行加载。更具体地说,此“弹簧”的初始位置为含碳帽的单壁碳纳米管近顶端分组经过上述驰豫后的质心位置,然后通过对“弹簧”施加图示坐标轴8中的Z轴正方向,即“弹簧”位移方向9实现对含碳帽的单壁碳纳米管的加载。其中,通过“弹簧”实现加载,采用一步一平衡的模拟机制,即加载一次,系统进行一次驰豫。最后统计含碳帽的单壁碳纳米管各个分组的三个方向质心坐标,以及弹簧整体的力等,完成对含碳帽的单壁碳纳米管拉伸的分子动力学模拟。

    11)最后可得到相应数据文件及模拟的图形文件等,根据需要对统计的数据信息等进行后处理,得到弹性模量、泊松比,应力应变关系等力学性能参数。

    综上所述,仅为本发明的具体实施方式,但发明的保护范围并不局限于此,任何熟悉本技术领域的工程技术人员在本发明的技术范围内,可做的一些变换,都应该作为侵犯本发明的保护范围。因此本发明的保护范围应该以权利要求书的保护范围为准。

    展开全文
  • 分子动力学方法用Tersoff势函数研究了α-SiO2晶体的力学性能。 结果表明,室温下单轴加载过程中,α-SiO2晶体经历了弹性变形,塑性变形和断裂变形。 α-SiO2在塑性变形中从晶相转变为非晶相。 并且通过研究温度对...
  • 采用分子动力学的方法,对无限稀释MBCl2z水溶液进行了模拟。计算了包括3个超临界点在内的7个状态下无限稀释MgCl2水溶液中Mg2+-O,Mg2+-H,Cl--O与Cl--H的径向分布函数,Mg2+-水(质心)和Cl--水(质心)的配位数。采用...

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 4,839
精华内容 1,935
关键字:

分子动力学模拟