data/attachment/portal/201111/06/0930426dg787ixizkmmwd2.gif
分子动力学
分子动力学是一门结合物理，数学和化学的综合技术。分子动力学是一套分子模拟方法，该方法主要是依靠牛顿力学来模拟分子体系的运动，以在由分子体系的不同状态构成的系综中抽取样本，从而计算体系的构型积分，并以构型积分的结果为基础进一步计算体系的热力学量和其他宏观性质。
目录
简史
基本步骤
作用势与动力学计算
时间步长与约束动力学
应用
模拟简述
参见
参考资料
简史
1957年：基于刚球势的分子动力学法（AlderandWainwright）data/attachment/portal/201111/06/093042we00ese1vyfvk1eh.gif
1964年：质点系への拡张(Rahman)
1971年：刚体系への拡张(RahmanandStillinger)
1977年：约束动力学方法（Rychaert等）
1980年：恒压条件下的动力学方法（Andersenの方法、Parrinello-Rahman法）
1983年：非平衡态动力学方法(GillanandDixon)
1984年：恒温条件下的动力学方法（能势━诈`バーの方法）
1985年：第一原理分子动力学法（→カー?パリネロ法）
1991年：巨正则系综的分子动力学方法(CaginandPettit)
基本步骤
确定起始构型
data/attachment/portal/201111/06/0930438w89zfltz99vlvnv.gif
分子动力学
进行分子动力学模拟的第一步是确定起始构型，一个能量较低的起始构型是进行分子模拟的基础，一般分子的起始构型主要来自实验数据或量子化学计算。
在确定起始构型之后要赋予构成分子的各个原子速度，这一速度是根据波尔兹曼分布随机生成的，由于速度的分布符合波尔兹曼统计，因此在这个阶段，体系的温度是恒定的。另外，在随机生成各个原子的运动速度之后须进行调整，使得体系总体在各个方向上的动量之和为零，即保证体系没有平动位移。
进入平衡相
由上一步确定的分子组建平衡相，在构建平衡相的时候会对构型、温度等参数加以监控。
进入生产相
进入生产相之后体系中的分子和分子中的原子开始根据初始速度运动，可以想象其间会发生吸引、排斥乃至碰
data/attachment/portal/201111/06/093045u37tutbocgxu3u3k.jpg
分子动力学
撞，这时就根据牛顿力学和预先给定的粒子间相互作用势来对各个粒子的运动轨迹进行计算，在这个过程中，体系总能量不变，但分子内部势能和动能不断相互转化，从而体系的温度也不断变化，在整个过程中，体系会遍历势能面上的各个点，计算的样本正是在这个过程中抽取的。+
计算结果
用抽样所得体系的各个状态计算当时体系的势能，进而计算构型积分。
作用势与动力学计算
作用势的选择与动力学计算的关系极为密切，选择不同的作用势，体系的势能面会有不同的形状，动力学计算所
data/attachment/portal/201111/06/09311529vs0ixsczlzy717.gif
分子动力学
得的分子运动和分子内部运动的轨迹也会不同，进而影响到抽样的结果和抽样结果的势能计算，在计算宏观体积和微观成分关系的时候主要采用刚球模型的二体势，计算系统能量，熵等关系时早期多采用Lennard-Jones、morse势等双体势模型，对于金属计算，主要采用morse势，但是由于通过实验拟合的对势容易导致柯西关系，与实验不符，因此在后来的模拟中有人提出采用EAM等多体势模型，或者采用第一性原理计算结果通过一定的物理方法来拟合二体势函数。但是相对于二体势模型，多体势往往缺乏明确的表达式，参量很多，模拟收敛速度很慢，给应用带来很大的困难，因此在一般应用中，通过第一性原理计算结果拟合势函数的L-J，morse等势模型的应用仍然非常广泛。
时间步长与约束动力学
分子动力学计算的基本思想是赋予分子体系初始运动状态之后利用分子的自然运动在相空间中抽取样本进行统计
data/attachment/portal/201111/06/0931162cc9dkbd11kzqjdc.gif
分子动力学
计算，时间步长就是抽样的间隔，因而时间步长的选取对动力学模拟非常重要。太长的时间步长会造成分子间的激烈碰撞，体系数据溢出；太短的时间步长会降低模拟过程搜索相空间的能力，因此一般选取的时间步长为体系各个自由度中最短运动周期的十分之一。
但是通常情况下，体系各自由度中运动周期最短的是各个化学键的振动，而这种运动对计算某些宏观性质并不产生影响，因此就产生了屏蔽分子内部振动或其他无关运动的约束动力学，约束动力学可以有效地增长分子动力学模拟的时间步长，提高搜索相空间的能力。
应用
分子动力学可以用于NPT，NVE，NVT等系综的计算，是一种基于牛顿力学确定论的热力学计算方法，与蒙特卡洛
data/attachment/portal/201111/06/09311828i04kcv19vk2xvh.gif
分子动力学
法相比在宏观性质计算上具有更高的准确度和有效性，可以广泛应用于物理，化学，生物，材料，医学等各个领域。
另外，在实际应用中，经常把分子动力学方法和蒙特卡罗法联合使用。
模拟简述
分子动力学计算机模拟是研究复杂的凝聚态系统的有力工具。这一技术既能得到原子的运动轨迹,还能象做实验一样作各种观察。对于平衡系统,可以在一个分子动力学观察时间(observationtime)内作时间平均来计算一个物理量的统计平均值,对于一个非平衡系统过程,只要发生在一个分子
data/attachment/portal/201111/06/093118gogor8gn6jshol9o.jpg
分子动力学
动力学观察时间内(一般为1～10ps)的物理现象也可以用分子动力学计算进行直接模拟。可见数值实验是对理论和实验的有力补充,特别是许多与原子有关的微观细节,在实际实验中无法获得,而在计算机模拟中可以方便地得到。这些优点使分子动力学在材料科学中显得非常有吸引力。因为许多人们感兴趣的领域,如晶格生长、外延生长(epitaxy)、离子移植、缺陷运动、无定型结构(amorphousstructure)、表面与界面的重构等问题,原则上都可以进行计算机模拟研究。分子动力学假定原子的运动是由牛顿运动方程决定的,这意味着原子的运动是与特定的轨道联系在一起的。当核运动的量子效应可以忽略,以及绝热近似严格成立时,分子动力学的这一假定是可行的。绝热近似也就是要求在分子动力学过程中,每一时刻电子均要处在相应原子结构的基态。大多数情形下,这一条件都是满足的。要进行分子动力学模拟需要知道原子间正确的相互作用势,从而必须知道相应的电子基态。电子基态的计算是一个非常困难的量子多体问题。好在密度泛函的引入使这方面的计算有了很大的简化,意味着我们可以把这一多体问题转化为一组自洽的单粒子轨道方程。对交换相关势采用局域密度近似,这组方程就实际可解了。这是目前凝聚态物理电子结构计算中普遍采用的方法。对于非强相关系统,局域密度泛函非常有效。在局域密度泛函基础上的第一原理性计算,在研究原子、分子、和晶体的结构中取得了巨大的成功。然而由于计算上的复杂性,基于局域密度泛函的第一性原理计算长期以来认为是不可能直接用于统计力学模拟,因为上百个原子的模拟就需要对104-106种原子构型作电子结构计算。因此,在分子动力学模拟中一般只能采用经验势来代替原子间实际作用势。一般就用Lennard-Jones势之类的二体作用势来替代。这种替代可能对隋性气体之类的系统进行模拟有效,但许多实际材料科学中感兴趣的问题都是多体效应的结果,不可能光用二体势代替能解决的。尽管人们对金属和共价系统的模型势作了许多尝试,并取得了一些成果,但一般一种模型势只能适用于一种原子,不能满足实际模拟的需要。
data/attachment/portal/201111/06/0931184lskpl6qsrgcjrgb.gif
分子动力学
另外经验势近似有一个本身的局限性,它丢失了局域电子结构之间存在着的强相关作用信息,也就是说,不能得到成键性质(bondingproperties),以及原子动力学过程中的电子性质。尽管在某些特殊情况下,电子性质也能通过经验势近似来得到,例如用经验势的分子动力学来计算原子结构时,选取某几种原子构型来作电子结构计算,这样不仅耗时,而且原子动力学和电子结构计算成了相互独立的过程。1985年,Car和Parrinello在传统的分子动力学中引入了电子的虚拟动力学,把电子和核的自由度作统一的考虑,首次把密度泛函理论与分子动力学有机地结合起来,提出了从头计算分子动力学方法(也称CP方法),使基于局域密度泛函理论的第一原理计算直接用于统计力学模拟成为可能,极大地扩展了计算机模拟实验的广度和深度。
参见
计算化学分子模拟en:moleculardynamicsja:分子動力学法nl:Moleculairedynamica
参考资料
http://wiki.keyin.cn/index.php/%E5%88%86%E5%AD%90%E5%8A%A8%E5%8A%9B%E5%AD%A6