分子模拟的一些基本资料.docx
- 文档编号:6339261
- 上传时间:2023-01-05
- 格式:DOCX
- 页数:19
- 大小:210.21KB
分子模拟的一些基本资料.docx
《分子模拟的一些基本资料.docx》由会员分享,可在线阅读,更多相关《分子模拟的一些基本资料.docx(19页珍藏版)》请在冰豆网上搜索。
分子模拟的一些基本资料
第1章前言
1.1研究背景
1.1.1分子模拟及其发展
分子模拟(MolecularSimulation)为二十世纪初发展起来的一种计算机模拟方法,它泛指用于模拟分子或分子体系性质的方法,主要用于探索研究具有三维结
构的分子结构和分子的性能[i]。
分子模拟是根据物理和化学的基本原理构建一个模型(通常是数学模型,是对某种分子体系或反应过程的理想化描述),建立一种以计算数据(由计算机来执行)来代替实验测量的研究方法,并获取相关的物理和化学信息。
分子模拟在材料科学方面的应用包括模拟材料的结构、计算材料的性质、预测材料的行为、验证实验结果(重现实验过程)、从微观角度认识材料,总之是为了更深层次理解材料的结构,认识材料的各种行为。
分子模拟的主要优势在于可以降低实验成本、具有较高的安全性、实现通常条件下较难或无法进行的实验(例如:
超低温,低于-100C;超高压,高于100Mpa)、研究极快速的反应和变化等。
R.S.Mulliken获得诺贝尔获奖时的感言是:
总之,我愿意强调我的信念:
计
算化学的年代已经到来,成千上百的化学家以计算机代替实验室,来获得众多的
化学信息。
唯一的障碍是你必须偿付机时费。
”阿基米德曾说过:
给我一个支
点,我就能够翘起地球。
”,但是分子模拟告诉我们:
“Givemeanenoughpowerfulcomputer,Icansimulatethewholeworld”。
从1980年开始,每年在EngineeringVillage中关于“MoleculaSimulation的文章数目由37篇递增到最高5209(2008年)篇。
与分子模拟有关的论文,美国(UnitedStates发表的篇数最多,高达16351篇,其次是日本,中国名列第三。
分子模拟作为一种计算机模拟技术,主要可以进行解释工作和预测工作。
前
者为实验奠定理论基础,通过模拟解释实验现象、建立理论模型、探讨过程机理等,后者为实验过程提供可能性和可行性研究,进行方案辅助设计、材料性能预
测、过程优化筛选等。
不同的分子模拟方法可以得到不同的信息。
量子力学模拟
方法可以计算得到分子的大多数性质,如结构、构象、偶极矩、电离能、电子亲和力、电子密度、过渡态和反应途径等;分子力学可以计算分子体系的稳定构象、热力学特性、振动光谱等;能量最小化可以探索相空间(phasespace和势能面
(potentialcurve),可以找出局部(local)与全局(global)的最小点及转化过程的马鞍点(saddlepoint);MonteCarlo可以计算复杂分子体系的结构变化,特别是相变化;分子动力学可以得到复杂分子的热力学性质、结构、力学性质,特别是可以观察体系的动态演变,得到许多与时间有关的热动力学性质;布朗动力学可以研究蛋白质在水溶液中的折叠过程;构象分析可以研究复杂分子稳态和亚稳态结构之间的演变等等。
分子模拟所涉及的研究领域,涵盖了物理、化学、化工、材料、生化等几乎一切可以通过建立理论模型进行研究的体系,多数能够得到与实验结果近似的计算结果,所以分子模拟己经逐渐成为与实验技术并重的强有力的研究手段。
分子模拟实际上并不仅仅局限于计算机模拟,但今天的分子模拟己和计算机模拟密不可分,正是由于计算机的高速计算技术的发展才使得分子模拟能够像今天这样发挥如此重要的作用。
利用计算先行以了解更多的分子特性,已成为合成
化学家和药物设计学家所依赖的重要方法。
化学家们通过这种方法可以设计出最
佳的反应途径,预测合成的可能性,且节省许多时间和避免材料的浪费。
分子模
拟除了在药物设计方面的应用外,已被广泛地应用于研究金属材料、无机材料、高分子聚合物材料、生物材料等复杂庞大的体系。
由于模拟结果随着计算机的发展与模拟方法(算法)的改良而更加精确,分子模拟已成为如今化学各个领域中不可缺少的工具。
1.1.2分子模拟的主要方法
分子模拟的主要模拟方法有量子力学模拟和经典力学模拟,量子力学模拟主要根据从头算方法、半经验方法、DFT(DensityFunctionalTheory,密度泛函数理论)方法,而经典力学模拟的方法主要依据分子力学、分子动力学、蒙地卡罗模拟、布朗动力学。
量子力学是利用波函数来研究微观粒子的运动规律的一个物理学分支学科,它以分子中电子的非定域化(delocalization)为基础,一切电子的行为以其波函数(wavefunction)表示。
根据海森伯(Heisenberg)的测不准原理(uncertainlyprinciple),量子力学可计算区间内电子出现的概率,其概率正比于波函数绝对值
的平方,而欲得到电子的波函数,则需解薛定谔方程式(Schr?
dingerequation),
即:
此式中,为薛定谔算子(是一些数学指令)、
量。
但是,分子中含有许多电子,要解此方程并非易事。
方面致力于改进量子计算的方法及增进其精准度,如
颁给的密度泛函数理论(densityfunctionaltheory,DFT),即为非常精确的量子计算方法;另一方面,科学家们积极研究提升计算机的计算速度,以期计算较多电
子的分子系统。
通过量子力学,我们可以研究离子、原子、凝聚态物质,以及原子核和基本粒子的结构、性质,它与相对论一起构成了现代物理学的理论基础,也正是量子力学的出现,许多现象才得以真正被解释,新的、无法被直觉想象出来的现象被预言,而且基本上所有的物理间的基本相互作用都可以在量子力学的框架内描写。
随着量子力学理论的逐步完善、经验力场的不断开发和更快更准确的算法以及计算机的普及和计算速度、计算量的不断提升,半个世纪以来,分子模拟的理论和方法得到了快速的发展,在化学化工、材料科学、医药科学、生命科学等诸多领域发挥着越来越重要的作用,并形成了一门专门的学科一一分子模拟。
今天,国内外的绝大多数的优秀科研机构和高等院校都已或正在涉足分子模拟领域,每
年在世界范围内发表的相关研究论文数以万计,并且呈不断上升趋势。
分子模拟
所担当的角色也由早期纯粹的解释型逐渐过渡到解释、指导及预测并重型。
通过上文,我们可以知道量子力学的方法适用于简单的分子或电子数目较少的体系。
但是在自然界与工业上的许多系统中,比如聚合物(脂肪、橡胶、安全玻璃等等)、生化分子(核酸、酶、蛋白质、多糖类等)、合金材料等均含有大量的原子及电子。
对于金属材料、聚合物材料、浓稠溶液、固态混合物、纳米材料等系统,我们不需要了解单一分子的性质及分子间的相互作用,我们需要了解的是整个系统的各种集合的性质、动态行为与热力学的性质。
对于这样的复杂体系,因其含有的电子数目较多,迄今为止,仍不可能完全仰赖量子力学计算。
为解决这样复杂庞大的体系,科学家们在二十世纪六十年代就开始着手研究各种可行的非量子力学计算方法,利用这些非量子力学的计算方法可以解决许多复杂系统的问题,并得到相当精确的结果。
而这些非量子力学计算方法通常指的就是经典力学模拟方法,即:
分子力学、分子动力学、蒙地卡罗模拟、布朗动力学。
分子力学方法(molecularmechanics,MM)是依据经典力学(classical
mechanicS)的计算方法。
因为用量子力学的方法去模拟计算大体系,如上千个原子的聚合物,是不实际的,即便是这种模拟能够进行,大多数情况下得到的大部分信息都要被舍弃。
这是因为在模拟大体系的时候,常常是为了得到统计性质,如模量、扩散系数等,而这些量主要依赖于原子核的位置,此时计算量极大的电
子运动信息并不重要。
根据量子力学的波恩-奥本海默近似,分子的能量可以近似看做构成分子的各个原子的空间坐标的函数,简单的讲就是分子的能量随分子
构型、相对位置变化而变化,而描述这种分子能量和分子结构之间关系的就是分子力场函数(forcefield)。
而分子力学就是主要依据分子的力场计算分子的各种特性。
力场函数中含有许多参数,这些参数可由量子力学计算或通过科学实验获得。
利用分子力学的方法可计算体系庞大的分子的稳定构象、热力学特性及振
动光谱等信息。
与量子力学相比较,此方法简便、快速,往往可快速得到分子的各种性质。
在某些情况下,由分子力学方法所得到的结果几乎可以与高阶量子力学方法所得到的结果一致,但其所需的计算时间却远远小于量子力学的计算。
故分子力学方法常被用于药物、团簇体、聚合物大分子的研究。
分子动力学模拟(moleculardynamicssimulation,MDS)是目前最广泛采用的计算庞大复杂体系的方法。
自1966年起,由于分子力学的发展,人们又系统的建立了许多适用于生化分子体系、聚合物、金属及非金属材料的力场,使得计算复杂体系的结构与一些热力学、光谱性质的能力及精准性大为提升。
分子动力学模拟是应用这些力场及根据牛顿运动力学原理所发展的计算方法。
此方法的优
点为精准性高,同时可以获得系统的动态与热力学统计资料,并可广泛地适用于各种系统及各类特性的研究。
与蒙地卡罗计算方法相比较,分子动力学模拟系统中粒子的运动有正确的物理依据,并且计算技巧经过许多修改,现已日趋成熟,由于其计算能力强,能满足各类问题的需求,许多使用方便的分子动力学模拟商业软件也已问世。
在先进国家的学校、工厂、医院等的实验室里,这些商业软件是不可或缺的研究工具。
虽然分子动力学模拟具有很多优点,但是由于分子动力学模拟计算需要引用数理积分方法,因此仅能研究系统在短时间内的运动,无法
模拟一些运动时间长的运动问题(如由很多氨基酸组成的大分子蛋白质的折叠)。
蒙地卡罗计算方法(MonteCarlomethod)为最早针对庞大系统所采用的非量子计算方法,即MC计算法。
蒙地卡罗计算法借由系统中质点(原子或分子)
的随机运动,结合统计力学的概率分配原理来得到体系的统计及热力学资料,即
根据待求问题的变化规律,构造合适的概率模型,然后进行大量统计实验,使模型的某些统计参量正好是待求问题的解。
MonteCarlo模型的建立可以分为三个
具体的操作步骤,首先将所研究的物理问题演变为类似的概率或统计模型;其次
通过数值随机抽样实验对概率模型进行求解,其中包括大量的算术运算和逻辑操作;最后使用统计方法对得到的结果进行分析处理。
这种计算方法可以研究复杂
体系及金属的结构及其相变化性质。
蒙地卡罗计算的弱点在于只能计算统计平均值,无法得到系统的动态信息。
此计算所依据的随机运动并不适于物理学的运动原理,且与其他的非量子计算方法相比并非特别的经济快速,因此自从分子动力
学计算逐渐盛行后,使用蒙地卡罗计算方法的人已逐渐减少。
还有一种与分子动力学模拟类似的计算方法为布朗动力模拟(Brownian
dynamicssimulation,BD)。
布朗动力模拟适用于大分子的溶液系统,计算中,将大分子的运动分为根据力场作用的运动与来自溶剂分子的随机力作用。
通过布
朗运动方程式可得到大分子运动的轨迹及一些统计与热力学的性质。
布朗动力模
拟通常适用于计算生化分子(如DNA、RNA、蛋白质、多肽)的水溶液。
此算法的优点在于能计算大分子在较长时间范围内(纳秒级)的运动,其缺点为将溶剂分子的运动视为布朗运动粒子的假设未必是正确的。
第2章分子动力学模拟法
分子模拟作为一种计算机模拟技术,主要可以进行解释工作和预测工作。
前者为实验奠定理论基础,通过模拟解释实验现象、建立理论模型、探讨过程机理等,后者为实验过程提供可能性和可行性研究,进行方案辅助设计、材料性能预
测、过程优化筛选等。
分子模拟所涉及的研究领域,涵盖了物理、化学、化工、材料、生化等几乎一切可以通过建立理论模型进行研究的体系,多数能够得到与
实验结果近似的计算结果,所以分子模拟己经逐渐成为与实验技术并重的强有力的研究手段。
20世界90年代以来,由于计算机科学技术的飞速发展,分子模拟的地位日显突出。
现在,理论分析、实验测定和模拟分析已成为现代科学的三种主要方法,在化工新产品、新材料和新药物的研究和开发中,采用分子模拟技术,可以从分子的微观性质推算及预测产品及材料的介观、宏观性质,已成为新兴学术方向。
为探求纳米尺度下传热现象的规律和内在机理,我们必须从微观细节着手,研究载热粒子的行为,并依据统计力学原理得到系统的统计性质。
分子模拟技术为某些问题(如统计力学)的精确解答起到非常重要的作用,否则这些问题很难或只能近似求解。
此外,一些在极温极压条件下很难或者根本无法进行的实验,利
用分子模拟却比较容易做到,在这种情况下,它实际上是对理论应用范围的一个较好的检测。
通过对比分子模拟的结果与真实的实验结果,分子模拟的结果可以
检验理论的正确性,相反,真实的实验结果与理论可以检测模型构建的合理性,此两者相互联系相互进步,通过对比,我们可以修正所建立的模型,并从模拟结果中进一步加深对所研究问题的理解,有助于解释出现的新现象。
2.1分子动力学概述
分子动力学模拟(MolecularDynamicsSimulation,MDS)是目前最广泛计算
庞大复杂体系所采用的方法,并且现已广泛应用于计算物理、化学、材料、生物科学等领域内的一种计算机模拟技术,它的思想是将组成系统的微观粒子视为牛顿经典粒子,将所研究的系统视为经典多体系统,通过选择恰当的场函数来描述系统内微观粒子间相互作用的势函数及系统外加约束条件后,利用牛顿运动定律
来求解微观粒子的牛顿运动方程。
在MDS中,势能模型的选取是非常重要的环节。
模拟能否成功取决于能否正确地选择势能模型。
对于简单分子,常采用硬球、软球、Lennard-Jones模型、Kihara模型等等。
对复杂分子,可以采用多中心的位置一位置相互作用模型,但中心间的相互作用仍为简单模型,对带电分子或离子,还要引入库仑作用,常用的多体作用势有Stillinger-Weber(SW)作用势、Tersoff
势等。
自1966年起,由于分子力学的发展,人们又系统的建立了许多适用于生化分子体系、聚合物、金属及非金属材料的力场,使得计算复杂体系的结构与一些热力学、光谱性质的能力及精准性大为提升。
分子动力学模拟是应用这些力场及根据牛顿运动力学原理所发展的计算方法。
此方法的优点为精准性高,同时可
以获得系统的动态与热力学统计资料,并可广泛地适用于各种系统及各类特性的研究。
与蒙地卡罗计算方法相比较,分子动力学模拟系统中粒子的运动有正确的物理依据,并且计算技巧经过许多修改,现已日趋成熟,由于其计算能力强,能满足各类问题的需求,许多使用方便的分子动力学模拟商业软件也已问世。
2.2分子动力学模拟的基本步骤
分子动力学模拟(MolecularDynamicsSimulation,MDS)的出发点是将组成系统的微观粒子视为牛顿经典粒子,将所研究的系统视为经典多体系统,通过选
择恰当的场函数来描述系统内微观粒子间相互作用的势函数及系统外加约束条件后,利用牛顿运动定律来求解微观粒子的牛顿运动方程,的到系统随时间演进
的微观过程,最后基于统计力学计算得到系统的各种参数。
分子动力学模拟一般由下面几个步骤组成:
1.选取要研究的系统及其边界和系统内粒子间的作用势(力场函数);
2.设定系统中粒子的初始位置和动量初始化;
3.建立模拟算法,计算粒子间作用力及各粒子的速度和位置,并求解Newton运动方程,得到下一时刻粒子的空间位置和速度;
4•当体系达到平衡后,依据相关的统计公式,获得各宏观参数和输运性质;
2.2.1选择恰当系统与边界
分子动力学模拟是在一定的系综下进行的,常见的系综有微正则系综、正则系综、等温等压系综和等温等焓系综等。
在统计物理中,系综(Ensembie代表一群类似的体系的集合。
对一类相同性质的体系,其微观状态(比如每个粒子的位置和速度)仍然可以大不相同。
实际上,对于一个宏观体系,所有可能的微观状态数是天文数字。
统计物理的一个基本假设(各态历经假设)是:
对于一个处于平衡的体系,物理量的时间平均等于对应系综里所有体系进行平均的结果。
体系的平衡态的物理性质可以对不同的微观状态求和来得到。
系综的概念是由约西亚・威拉德•吉布斯(J.WillardGibbs)在1878年提出的。
在分子动力学模拟过程中,经常采用的系综有微正则系综(MicrocanonicalEnsemble)、正贝U系综(CanonicalEnsembie。
正则系综(CanonicalEnsemble也称为NVT系综,它是统计力学中系综的一种。
它代表了许多具有相同温度的体系的集合,该系统具有恒定的粒子数N、体
积V和温度T。
正则系综是最普遍应用的系综。
通常,系综内每个体系的粒子数和体积都是相同的。
但每个体系都可以和系综内其他体系交换能量。
同时系综里
所有体系的能量总和,以及所有体系的总个数是固定的。
但是正则系综的恒温条
件是不能自发实现的,因此在分子动力学模拟过程中经常会采用一定的措施来保证系统温度的恒定。
在分子动力学模拟过程中,系统的温度T又Boltzmann能量均分定理得出川,即:
—(2.错误!
仅主文档。
)
剳k呂
231
式中,N代表区域中原子的数量,Vi为粒子i的速度,kB=1.38X10J・K。
此公式成立的前提是系统区域内达到热力学平衡,若系统区域内无法达到热力学平衡,温度的定义就失去了意义。
对于NVT系综,调节温度恒定的有效方法是直接调节系统内微观粒子的速度。
在这些条件下,当系综内所有体系被分配到不同的微观态上,我们发现:
每
个微观态上的体系个数正比于exp(-BE)。
其中B=1/kBT,kB是玻尔兹曼常数,T是绝对温度。
正则系综的配分函数是:
配分函数的对数就是亥姆霍兹自由能(HelmholtzFreeEnergy,符号F),当配分函数计算出以后,平均能量可以直接从lnZ对B一阶导数中求得。
(2.错误!
仅主文档。
)
微正则系综(MicrocanonicalEnsemble)又称为NVE系综,该系综由许多具有相同能量E,粒子数N,体积V的体系的集合。
微正则系综可以通过平衡态分子动力学由系统内微观粒子的牛顿运动方程自动实验,比较简单方便,因此多数平
衡态分子动力学模拟是在微正则系综下进行的。
对于分子动力学模拟除了要选定合适的系综之外,还需要对模拟体系加以适当的边界条件。
这是因为材料的宏观性质是由大量粒子的相互作用所表现出来的,在模拟过程中,如果所建立的模型含有足够多的粒子,那么才能真实的模拟
出宏观材料的性质。
但是计算机的计算能力是有限的,分子动力学模拟时粒子数目一般不超过10000个,加上系综的某些限制性条件(例如NVT系综中,调节
其表面原子与体内原子之比翌很大因此会造成明显的表明效应。
因此,为了减小小体系的影响且使计算中系统的密度维持恒定,在分子动力学模拟过程中通常会采用有效的边界条件。
常用于分子动力学模拟中的边界条件有三种:
周期性
;厶-A一^L
边界条件、对称性边界条件和固壁边界条件川O
在周期性边界条件下,位于图
中央的盒子表示所模拟的系统,其周围的盒子与模拟系统具有相同的排列及运
动,我们称之为周期性镜像(Periodicmirrorimage)系统,当计算系统中任意微粒移出盒子后,则必有一粒子由相对方向移入,这样的限制条件使得系统中的粒子数目维持恒定,密度不变,符合实际要求。
在周期性边界条件下,计算系统中分子间的作用力变得复杂,我们通常采取最近镜像法(Nearestmirrorimage)来简化问题。
如图2-错误!
仅主文档。
周期性边界示意图。
图2-错误!
仅主文档。
周期性边界示意图
2.2.2分子间作用力的描述一一力场
分子力场(Forcefields)根据量子力学的波恩-奥本海默近似可以将力场看作是粒子势能面的经验表达式,即将一个分子的势能近似看作构成分子的各个原子的空间坐标的函数,简单地讲就是分子的势能随分子构型的变化而变化,而描述
这种分子势能和分子结构之间关系的就是分子力场函数。
分子动力学模拟过程中计算系统内粒子间的作用力主要依靠的就是力场给各个粒子间分配的作用力模型,分子动力学模拟的主要依据是经典力学,那么分子动力学模拟结果的准确性
也主要取决于力场函数对粒子间相互作用力描述的正确与否。
由于分子动力学模
拟研究的广泛性,分子间相互作用模型并非完全相同,针对不同的分子类型、分子体系,分别有不同的力场函数可以使用。
分子力学与分子动力学发展至今已有40多年,其计算的系统由原来最初的单原子分子系统延伸至多原子分子、聚合物分子、生化分子系统。
计算时所需的力场也随着系统的复杂度的增加而增加,在执行分子动力学计算时,选择合适的
力场十分重要。
常见的力场有CVFF、CFF91、CFF95、PCFF、AMBER、CHARMM、COMPASS和MM形态力场等。
MM形态力场为Allinger等所发展,此力场适合用于模拟有机化合物、自由基、离子等随着分子模拟的进步MM形态力场也发展出MM2、MM3和MM4形态力场等,MM力场将一些常见的原子细分(例如,将碳原子分为sp3、sp2、
sp、酮基碳、碳自由基和碳阳离子等并且这些不同形态的碳原子具有不同的力场参数),并为这些原子补充更多的力场参数,以更精确详细的来描述相同类型、不同电子结构原子的受力。
CVFF力场其全名为一致性价力场(Consistentvaleneeforcefield,为Dauber
Osguthope等所发展。
此力场以计算系统的结构与结合能最为准确,亦可以提供合理的构型能和振动频率,常用于计算氨基酸、水及含各种官能团的分子体系,后经过不断改进,CVFF力场已适用于计算各种多肽、蛋白质与大量的有机分子。
AMBER力场(Assistedmodelbuildingwithenergyminimization)为美国加
州大学的PeterKollman等所发展,此力场通常可以得到分子的几何结构、构型能、振动频率与溶剂化自由能(Solvationfreeenergy)等数据,所以常用于较小的蛋白质、核酸、多糖等生化分子的模拟计算。
CHARMM力场(ChemistryatHarvardmacromolecularmeehanic$为美国哈佛大学所发展,此力场可研究许多分子系统,例如小的有机分子、溶液、聚合物、生化分子等。
此力场参数除来自计算结果与实验值的比对外,还引用了大量的量
子计算结果。
应用此力场能够得到与实验值相近的结构、作用能、构型能、转动能、振动频率、自由能与许多和时间有关的物理量。
随着力场应用研究的深入,随后出现了第二代力场。
第二代力场的形式远比经典里场复杂,并需要大量的力场参数。
其设计的目的是为了更加精确的计算分子的各种性质、结构、光谱、热力学特性、晶体特性等,其力场参数除了引用大量的实验数据之外,并参照了精确的量子计算结果,尤其适用于有机分子系统。
CFF91、CFF95、COMPASS与PCFF都属于第二代力场。
一般情况下,分子的总能量为动能与势能的和,而分子的势能可表示为简单的几何坐标的函数。
复杂分子的总势能一般可分为各类型势能的总和,即由成键
项(Ubonded)、父叉项(Ucross)和非键项(Unobonded)三部分构成,即:
总势能=成键项势能+交叉项势能+非键项势能
习惯上
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 分子 模拟 一些 基本资料