分子模拟简介
- 格式:pdf
- 大小:4.67 MB
- 文档页数:58
化学分子动力学模拟的原理和应用随着计算机技术的不断发展和进步,分子模拟技术在化学、物理、生物等学科中得到了广泛的应用,其中分子动力学模拟是其中比较重要的一种方法。
分子动力学模拟是一种数值模拟技术,利用分子动力学方程模拟分子之间的相互作用和运动规律,从而揭示分子的结构、性质、运动和相互作用等,能够对活性物质的设计与评价起到重要的作用。
一、分子动力学模拟的原理分子动力学模拟是一种基于牛顿力学的方法,它使用运动方程来描述在各种外部场下,分子的运动轨迹。
既反映了分子中各个原子之间的相互作用,也体现了整个系统的运动规律。
简单来说,分子动力学模拟是在已知原子间作用势和运动方程的条件下,以数值方法计算分子的运动和结构的方法。
分子动力学模拟的基本步骤分为以下几部分:1、布朗运动模拟模拟分子在溶液中的布朗运动,通过计算分子的位置和速度之间的关系,可以得出分子受到的作用力。
2、势函数计算计算分子所受到的各个势函数,如位能、马德隆势等。
3、运动方程求解根据分子所受到的力以及它们相互之间的运动规律,求解运动方程,对数值解得出各点的位置和速度。
4、相互作用计算对于每两个相互作用的粒子,根据其位置和速度计算出与一点位置的距离,再代入相互作用的势函数,最后计算出所有相互作用的和。
5、轨迹预测根据初始条件以及数学模型,预测出分子的轨迹和状态,最后得出分子的结构、动力学和热力学等性质。
二、分子动力学模拟的应用分子动力学模拟的应用十分广泛,不同领域有所不同的应用。
下面列举出几个典型的应用场景。
1、药物发现在新药研发过程中,研究分子相互作用和分子构象改变等问题十分重要。
使用分子动力学模拟,可以得到分子的能量、熵、电荷分布等信息,为药物设计和评价提供依据。
2、材料开发分子动力学模拟可以用于模拟材料的力学性能、热导性能和光学性能等。
例如,可以用此模拟在不同应力下的金属疲劳,探究其疲劳机理。
3、化学反应机理在化学反应中,可以使用分子动力学模拟来研究各个物种之间的反应,从而探讨反应的机理。
利用分子动力学模拟研究材料力学性能引言:材料力学性能是指材料在外力作用下的变形和破坏行为,对于材料设计和应用具有重要意义。
随着计算机技术的发展,分子动力学模拟成为一种有效研究材料力学性能的方法。
本文将以分子动力学模拟为基础,探讨材料力学性能的研究方法和应用。
1. 分子动力学模拟简介分子动力学模拟是一种以牛顿力学为基础,通过模拟材料中原子和分子的运动来研究材料性质的方法。
它可以实现对材料的结构、能量、力学性能等多个方面进行研究。
分子动力学模拟基于原子的尺度,通过求解牛顿运动方程,模拟原子间的相互作用和运动轨迹。
2. 分子动力学模拟在材料力学性能研究中的应用2.1 材料的力学性质研究分子动力学模拟可以通过计算应力应变曲线、弹性模量、屈服强度和断裂韧性等力学性质,来评估材料的力学性能。
在研究中,我们可以改变原子的初始位置、应用外力或者改变温度等条件,模拟材料的力学响应过程,并得到与实验结果相符合的数据。
这些数据可以为材料设计和工程应用提供重要参考。
2.2 界面和缺陷研究材料中的界面和缺陷对力学性能有重要影响。
利用分子动力学模拟可以揭示界面的形成、稳定性和力学行为等。
例如,通过模拟晶界和颗粒界面的生成和迁移过程,可以研究材料的晶界塑性和断裂行为。
此外,分子动力学模拟还可以研究缺陷如位错、空位和夹杂物等对材料力学性能的影响。
3. 分子动力学模拟的优势和限制3.1 优势分子动力学模拟具有以下优势:①可以模拟材料在原子级别的行为,可以提供精确的局部信息;②能够模拟多尺度问题,从纳米到宏观材料性能;③可通过改变条件来预测材料在不同环境下的性能;④节省成本和时间,避免传统实验中的大量试错。
3.2 限制分子动力学模拟也存在一些限制:①需要建立合理的原子模型和相互作用势函数,这对于复杂体系来说是挑战;②受到计算资源等约束,模拟时间和空间尺度有限;③需要验证和修正模拟结果,因为理论模型和实验结果之间可能存在差异。
4. 分子动力学模拟研究材料力学性能的前景随着计算能力的提高和模拟方法的不断发展,分子动力学模拟将在材料力学性能研究中发挥更重要的作用。
分子动力学模拟的原理和计算方法分子动力学模拟是一种用于研究分子、原子以及离子等微观粒子在时间和空间上的运动行为的计算方法。
它可以帮助科学家们更好地理解物质的性质和行为,对材料科学、化学、生物学等学科的研究起到了重要的推动作用。
分子动力学模拟的基本原理是基于牛顿力学和统计物理学的原理。
牛顿力学描述了物体的运动规律,而统计物理学则研究了大量微观粒子的整体行为。
分子动力学模拟将这两者结合起来,通过经典力学的运动方程对微观粒子的运动进行模拟与计算。
在分子动力学模拟中,首先需要确定系统的边界条件和初始状态。
边界条件包括系统的尺寸、形状以及宏观环境的温度和压力等。
初始状态则是指系统中各个微观粒子的初始位置和动量。
接下来,通过数值积分方法求解牛顿运动方程。
分子动力学模拟中最常用的数值积分算法是Verlet算法和Leapfrog算法。
这些算法根据粒子的当前位置、速度和加速度等信息,经过一段时间步长的迭代计算,更新粒子的位置和速度。
通过不断迭代计算,分子动力学模拟可以模拟微观粒子在时间上的演化过程。
在每个时间步长内,模拟中的粒子会受到相互作用力的影响,从而改变其位置和动量。
这些相互作用力包括分子间相互作用力、静电相互作用力以及外界外力等。
分子动力学模拟还可以通过引入一些其他的技术和手段来增加计算的准确性和效率。
其中一项常用的技术是周期边界条件,通过在系统的边界上连接系统的各个边界,模拟无限大系统。
另外,还可以利用Monte Carlo方法和多尺度模拟等技术来处理一些特殊的系统和问题。
分子动力学模拟不仅仅是一种计算方法,更是一种对物质和自然现象深入理解的工具。
通过分子动力学模拟,科学家可以观察到一些实验无法观察到的细节,揭示了物质的微观行为和特性。
例如,可以通过模拟水分子的运动来研究水的溶解性和扩散性质,可以模拟蛋白质的折叠过程来研究生物分子的结构和功能等。
分子动力学模拟虽然具有很强的理论基础,但同时也面临着一些挑战和限制。
化学软件——HyperChem(分子模拟)介绍HyperChem是一款以高质量,灵活易操作而闻名的分子模拟软件。
通过利用3D对量子化学计算,分子力学及动力学进行模拟动画,HyperChem 为您提供比其它 Windows 软件更多的模拟工具。
图形界面:图形界面,有半经验方法( AM1 , PM3 等), UHF , RHF 和 CI 和7.0 版新增加的密度泛函。
可进行单点能,几何优化,分子轨道分析,预测可见- 紫外光谱,蒙特卡罗和分子力学计算。
主要功能:1. 结构输入和对分子操作。
2. 显示分子。
3. 化学计算。
用量子化学或经典势能曲面方法,进行单点、几何优化和过渡态寻找计算。
可以进行的计算类型有:单点能,几何优化,计算振动频率得到简正模式,过渡态寻找,分子动力学模拟 Langevin 动力学模拟, Metropolis Monte Carlo 模拟。
支持的计算方法有:从头计算,半经验方法,分子力学,混合计算。
4. 可以用来研究的分子特性有:同位素的相对稳定性;生成热;活化能;原子电荷; HOMO-LUMO能量间隔;电离势;电子亲和力;偶极矩;电子能级; MP2 电子相关能; CI 激发态能量;过渡态结构和能量;非键相互作用能; UV-VIS 吸收谱; IR 吸收谱;同位素对振动的影响;对结构特性的碰撞影响;团簇的稳定性。
5. 支持用户定制的外部程序。
6. 其它模块: RAYTRACE 模块, RMS Fit , SEQUENCE 编辑器,晶体构造器;糖类构造器,构像搜寻,QSAR 特性,脚本编辑器。
7. 新的力场方法: Amber 2 , Amber 3 ,用于糖类的 Amber , Amber 94 ,Amber 96 。
8. ESR 谱。
9. 电极化率。
10. 二维和三维势能图。
11. 蛋白质设计。
12. 电场。
13. 梯度的图形显示。
14. 新增功能:密度泛函理论 (DFT) 计算; NMR 模拟;数据库; Charmm 蛋白质模拟;半经验方法TNDO ;磁场中分子计算;激发态几何优化; MP2 相关结构优化;新的芳香环图;交互式参数控制;增强的聚合物构造功能;新增基组。
分子动力学模拟的原理简介授课人:杨俊升博士
内容
1 2 4
分子模拟的应用
分子动力学计算的原理分子动力学模拟实例
3体系模型构建
一、分子模拟的应用
1.分子模拟概述
定义:
计算机辅助试验技术,以原子水平的分子模型来模拟分子的结构与行为、体系的各种物理化学性质。
2.分子模拟的作用
模拟材料的结构
计算材料的性质
预测材料的行为
验证试验结果
重现试验过程
从微观角度认识材料
总之,是为了更深层次理解结构,认识各种行为。
介观动力学分子力学、动力学
量子力学
密度泛函理论Walter Kohn E ρ[]=T o ρ[]+U ρ[]+E xc ρ[]
[])
()()]([,,,2r r r n v k i k i k i eff ϕεϕ=+∇)]
([)],...([1r n E r r E N =ψr
d r n N r r f r n i i i i ⎰∑Ω=⋅=3*)()
()()( ϕϕHohenberg-Kohn
theorem Kohn-Sham
equations •Exact only for ground state •Needs approximation to E xc
荣获1998年的诺贝尔化学奖
这三位科学家结合经典和量子物理学,设计出多尺度复杂化学系统模型,将传统的化学实验搬到了网络世界。
第一原理研究领域包括:
✓晶体材料结构优化及性质研究(半导体、陶瓷、金属、分子筛等)
✓表面和表面重构的性质、表面化学
✓电子结构(能带、态密度、声子谱、电荷密度、差分电荷密度及轨道波函分析等)✓晶体光学性质(包括EELS, XANS, XES)
✓材料热力学参数计算
✓点缺陷性质(如空位、间隙或取代掺杂)、扩展缺陷(晶体晶界、位错)
✓磁性材料研究
✓材料力学性质研究
✓材料逸出功及电离能计算
✓STM图像模拟
✓红外,拉曼光谱模拟
✓反应过渡态计算
✓动力学方法研究扩散路径
A b s o r p t i o n (c m -1)图1 (a )本征LN 晶体; (b )Mn 替代Li 位
LN 晶体; (c )Mn 替代Nb 位LN 晶体;
System Volume/Å3E total /eV LN
1232.98-73221.751Mn@LN-1
1240.78-73678.119Mn@LN-21225.65-72314.594
从上个世纪九十年代初期以来,计算机模拟技
术得到了飞速发展,主要基于三个方面的发展: 分子力场的发展(基石)
(Amber,OPLS、Compass)
原子间的键长、键角、分子间的内聚能等
模拟算法(途径)
计算机硬件(工具)
HPCx
二、分子动力学计算的原理
力场的概念:
分子力场是原子分子尺度上的一种势能场,它描述决定着分子中原子的拓扑结构和运动行为。
在量子力学上说,描述决定分子中原子的拓扑结构与运动行为的是分子的电子结构的本征性质。
也就是说,分子的基态波函数确定了该体系的分子力场。
从经典力学上说,分子力场是由一套是函数与一套力常数构成的。
这种简单的数学表达形式表示的势能函数称为力场。
PET
范德瓦尔斯项
键伸缩项双面角扭曲项
键角弯曲项离平面振动项库伦作用项对于复杂的分子体系:
nb b el V V V V V V V θφχ=+++++
什么是分子动力学模拟?
分子动力学模拟(molecular dynamics simulation)计算,简称MD计算。
分子动力学方法也可以叫做分子动态法。
由于分子力场所描述的是静态性质的势能,而实际的分子体系除了受势能影响以外,还得受外部环境,如温度,压强等的影响。
所以分子动力学方法是根据力场及牛顿运动学原理所发展的计算方法。
S r ()
三、体系模型构建
Menu
Toolbar Property
View
Project
Job status
晶胞,表面和团簇模型搭建
团簇建模工具界面
Fe
2
O
3
(001) Au晶胞
Au团簇
+
PET
SiO 2
PET/SiO 2
聚合物纳米复合材料模型搭建
聚合物单链及聚合物和水构型搭建
聚合物和水质量比10:90
搭建的初始模型
聚合物单链初始模型
聚合物建模工具界面
Amorphous Cell 建模界面
四、分子动力学模拟在高分子体系中的应用实例
无定型聚合物体系在拉伸场下的断裂机理
交联环氧树脂(2006)
Polymer
选择COMPASS力
场
Mechanical and Thermal Properties of Amorphous Polysulfone
•无定形聚砜
•MM 和MD 模拟得到力学和热学性质;•预测纯气体和混合物的吸附等温线和吸附位.
•聚砜密度1.2g/cm³. NPT 系统MD 计算平均体积
计算得到的力学和热学性质:
Young's Modulus = 3.88 +/-1.51 GPa Poisson's Ratio = 0.24 +/-0.15 GPa Bulk Modulus = 3.02 +/-1.48 Gpa Thermal Expansion Coefficient = 2.16 x 104 (K-1) Fan C.F ., Hsu S.L.,
Macromolecules, 25, 266, 1992
The picture shows a model of the bulk amorphous polymer polysulfone
DKI Ltd. Polyethersulfone
Amorphous Cell Mechanical Properties
Young's Modulus
Poisson's Ratio
Bulk Modulus
Thermal Expansion
Coefficient
Molekulares Modellieren in Polymeresystemen, 1 (2000) 78
无定形和半结晶聚合物的力学性质
交联环氧树脂的结构与性质关系
BHP Steel、RMIT University 交联环氧树脂可以作为钢铁上的底漆,BHP钢铁公司希望能够
深入了解环氧树脂的结构/性能之间的关系,从而设计出具有更好阻
隔性能,同时与钢铁基底结合更紧密的交联材料。
RMIT大学和BHP钢铁公司合作,使用分子动力学方法考察了环
氧树脂在不同固化剂作用下的交联,并得到了相关的交联密度以及
发生交联反应的位点数目等信息。
此外,树脂在固化过程中的收缩情况也得到了重现,并估计了
交联树脂的阻隔性能以及树脂与基底的相互作用强弱。
Polymer 43, (2002),963-969
CYMEL 1158CYMEL 1172
环氧树脂P0 (n=0), P1(n=1)
交联后的体积收缩
Before After
交联树脂的阻隔性能
>-<=∑=n
i n n r t r n t MSD 1
2)0()(1)()(61MSD t
D ∂∂=
与基底的结合能力
结论
通过模拟,能够得到以下关键信息:
交联密度和未反应交联位点数目。
对于性能设计非常重要。
交联树脂的收缩率。
与潜在涂层材料的合成和环境性能关
系密切。
阻隔性能。
对于考察防腐涂层的性能非常关键。
粘附能力。
涂层设计的重要指标。
Polymer 43, (2002),963-969
无定形海藻糖-水混合物的玻璃化转变温度模拟
剑桥大学化学系无定形碳水化合物因其能够形成稳定胶囊,性能稳定并持续释放药物活性组分,因此常被用在药品与食品工业中。
此类材料的优势在于高纯度,低毒性,良好的玻璃体形并具有较高的玻璃化转变温度(Tg)。
具有较高玻璃化转变温度(Tg)的碳水化合物,在高温下其物理性能也更稳定。
注入高Tg无定形碳水化合物中的药物可以在室温下进行储藏,不需要冷藏就可以进行运输。
例如海藻堂就具有较高的玻璃化转变温度。
海藻糖分子
75个海藻糖分子与64个水分子
的共混体系密度vs.模拟时间图。
其中370K为该浓度下混合物的Tg温度
搭建任意配比的海藻糖混合体系,并计算得到相应的Tg温度。
2.9%(w/w)H2O混合体Tg温度(实验) 2.9%(w/w)H
O混合体Tg温度(模拟)
2。