《分子力场简介精选PPT.ppt》由会员分享，可在线阅读，更多相关《分子力场简介精选PPT.ppt（66页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、关于分子力场简介第1页，讲稿共66张，创作于星期日Interatomic interactions原子与分子层次上的物质第2页，讲稿共66张，创作于星期日2003年諾貝爾化學獎年諾貝爾化學獎:-細胞膜通道之謎細胞膜通道之謎2003年的諾貝爾化學獎，頒給了兩位美國科學家：年的諾貝爾化學獎，頒給了兩位美國科學家：約翰霍普金斯大學醫學院的阿格雷（約翰霍普金斯大學醫學院的阿格雷（Peter Agre）與）與洛克斐勒大學霍華休斯醫學研究中心的麥金農洛克斐勒大學霍華休斯醫學研究中心的麥金農（Roderick MacKinnon）。他們獲獎的研究都與細）。他們獲獎的研究都與細胞膜上的通道有關，瑞典皇家科學院
2、在胞膜上的通道有關，瑞典皇家科學院在10月月8日發佈日發佈的新聞稿中指出，阿格雷是因為發現水通道與的新聞稿中指出，阿格雷是因為發現水通道與麥金農在離子通道的結構與機制上的研究，而麥金農在離子通道的結構與機制上的研究，而共享共享20032003的諾貝爾化學獎。的諾貝爾化學獎。第3页，讲稿共66张，创作于星期日The 2003 Nobel prize for Chemistry2003諾貝爾化學獎頒給了兩位醫生，或許有些奇怪，然而仔細的去瞭解他們所做的工作就會發現，其實他們的研究已經深入到了原子的層次，那不是化學又會是什麼呢！這其實也告訴了眾多對生命科學有興趣的年輕人，其實真正精采的還是在這個化學
3、的層次。第4页，讲稿共66张，创作于星期日The 2003 Nobel prize for Chemistry水分子通過 aquaporin AQP2。由於通道中心的正電荷，如H3O+般帶正電荷的離子會被驅逐，這可避免質子的滲漏。第5页，讲稿共66张，创作于星期日Water channels in cell membraneswater permeation movie for the 2003 Nobel prize in Chemistry第6页，讲稿共66张，创作于星期日分子力场分子力场分子力场分子力场:根据量子力学的波恩-奥本海默近似，一个分子的能量可以近似看作构成分子的各个原子的空间
4、坐标的函数，简单地讲就是分子的能量随分子构型的变化而变化，而描述这种分子能量和分子结构之间关系的就是分子力场函数。第7页，讲稿共66张，创作于星期日分子力学和分子力场 很多情况下，对大分子体系的处理可以完全避免使用量子化学计算。分子力学模拟使用古典力学模型（例如谐振子）描述化合物的能量。分子力学模型的所有常数均通过实验数据或第一原理计算结果得到。参数和函数的优化结果称为分子力场。进行参数化的化合物库对分子力学方法的计算成功与否至关重要。针对某类分子优化的力场只有在应用于同类分子时才可保证得到可信的结果。第8页，讲稿共66张，创作于星期日分子力场分子力场 定义:A classical force
5、field consists of:1.analytical functional forms describing the interactions 2.parameters in these functional forms The forcefield is fit to experimental data from a small set of molecules,with the goal that the forcefield will accurately interpolate and extrapolate to other systems.The classical for
6、cefields are fit to reproduce observables such as structures and relative energies,so the significant quantum mechanical effects(i.e.covalent bonds)are included empirically 第9页，讲稿共66张，创作于星期日分子力场Force field(also called a forcefield)refers to the functional form and parameter sets used to describe the
7、 potential energy of a system of particles(typically but not necessarily atoms).第10页，讲稿共66张，创作于星期日分子力场分子力场分子力场函数为来自实验结果的经验公式，可以讲对分子能量的模拟比较粗糙，但是相比于精确的量子力学从头计算方法，分子力场方法的计算量要小数十倍，而且在适当的范围内，分子力场方法的计算精度与量子化学计算相差无几，因此对大分子复杂体系而言，分子力场方法是一套行之有效的方法。第11页，讲稿共66张，创作于星期日分子力场分子力场分子动力学蒙特卡罗方法分子对接等分子模拟方法中有着广泛的应用.第12页
8、，讲稿共66张，创作于星期日分子力场分子力场A force field is used to minimize the bond stretching energy of this ethane molecule.第13页，讲稿共66张，创作于星期日分子力学Molecular mechanics refers to the use of Newtonian mechanics to model molecular systems.The potential energy of all systems in molecular mechanics is calculated using forc
9、e fields.Molecular mechanics can be used to study small molecules as well as large biological systems or material assemblies with many thousands to millions of atoms.第14页，讲稿共66张，创作于星期日第15页，讲稿共66张，创作于星期日分子力学v分分子子力力学学，又又叫叫力力场场方方法法（force field method），目目前前广广泛泛地地用用于于计计算算分分子子的的构象和能量。构象和能量。第16页，讲稿共66张，创作于
10、星期日分子力学v分子力学的基本思想分子力学的基本思想-1930，D.H.Andrews 在在分分子子内内部部，化化学学键键都都有有“自自然然”的的键键长长值值和和键键角角值值。分分子子要要调调整整它它的的几几何何形形状状（构构象象），以以使使其其键键长长值值和和键键角角值值尽尽可可能能接接近近自自然然值值，同同时时也也使使非非键键作作用用（van der Waals力力）处处于于最最小小的的状状态态，给给出出原原子子核核位位置置的的最最佳排布。佳排布。第17页，讲稿共66张，创作于星期日分子力学v分子的经典力学模型分子的经典力学模型-1946，T.L.Hill T.L.Hill提出用提出用va
11、n derWaals作用能和键长、键角的变形能来作用能和键长、键角的变形能来计算分子的能量，以优化分子的空间构型。计算分子的能量，以优化分子的空间构型。Hill指出：指出：“分子内部的空间作用是众所周知的，（分子内部的空间作用是众所周知的，（1）基团或）基团或原子之间靠近时则相互排斥；（原子之间靠近时则相互排斥；（2）为了减少这种作用，基团）为了减少这种作用，基团或原子就趋于相互离开，但是这将使键长伸长或键角发生弯曲，或原子就趋于相互离开，但是这将使键长伸长或键角发生弯曲，又引起了相应的能量升高。最后的构型将是这两种力折衷的结又引起了相应的能量升高。最后的构型将是这两种力折衷的结果，并且是能量
12、最低的构型果，并且是能量最低的构型”。第18页，讲稿共66张，创作于星期日分子力学v分子力学的发展分子力学的发展 虽然分子力学的思想和方法在虽然分子力学的思想和方法在40年代就建立起来了，年代就建立起来了，但是直到但是直到50年代以后，随着电子计算机的发展，用年代以后，随着电子计算机的发展，用分子力学来确定和理解分子的结构和性质的研究才分子力学来确定和理解分子的结构和性质的研究才越来越多。直到这时，才可以说分子力学已成为结越来越多。直到这时，才可以说分子力学已成为结构化学研究的重要方法之一。构化学研究的重要方法之一。第19页，讲稿共66张，创作于星期日分子力学v分子力学的发展分子力学的发展 近
13、几年来，随着现代技术的发展和应用，特别是计算近几年来，随着现代技术的发展和应用，特别是计算机技术的发展，分子力学方法已不仅能处理一般的中机技术的发展，分子力学方法已不仅能处理一般的中小分子，也不仅主要应用于有机化学领域，而且能处小分子，也不仅主要应用于有机化学领域，而且能处理大分子体系。在其他的一些领域，如生物化学、药理大分子体系。在其他的一些领域，如生物化学、药物设计、配位化学中，都有了广泛的应用。物设计、配位化学中，都有了广泛的应用。第20页，讲稿共66张，创作于星期日分子力学v分子力学的基本假设The Born-Oppenheimer Approximation-原原子子核核的的运运动动
14、与电子的运动可以看成是独立的；与电子的运动可以看成是独立的；分分子子是是一一组组靠靠各各种种作作用用力力维维系系在在一一起起的的原原子子集集合合。这这些些原原子子在在空空间间上上若若过过于于靠靠近近，便便相相互互排排斥斥；但但又又不不能能远远离离，否否则则连连接接它它们们的的化化学学键键以以及及由由这这些些键键构构成成的的键键角角等等会会发发生生变变化化，即即出出现现键键的的拉拉伸伸或或压压缩缩、键键角角的的扭扭变变等等，会会引引起起分分子子内内部部应应力力的的增增加加。每每个个真真实实的的分分子子结结构构，都都是是在在上上述述几几种种作用达到平衡状态的表现。作用达到平衡状态的表现。第21页，
15、讲稿共66张，创作于星期日分子力学v分子的空间能分分子子力力学学从从几几个个主主要要的的典典型型结结构构参参数数和和作作用用力力出出发发来来讨讨论论分分子子结结构构，即即用用位位能能函函数数来来表表示示当当键键长长、键键角角、二二面面角角等等结结构构参参数数以以及及非非键键作作用用等等偏偏离离“理理想想”值值时时分分子子能能量量（称称为为空空间间能能，space energy）的的变变化化。采采用用优优化化的的方方法法，寻寻找找分分子子空空间间能能处处于于极极小小值值状状态态时时分分子子的的构型。构型。第22页，讲稿共66张，创作于星期日分子力学v分子的空间能分子的空间能分子的空间能Es可表示
16、为：可表示为：Es=Ec+Eb+Et+Enb+其其中中Ec是是键键的的伸伸缩缩能能，Eb是是键键角角弯弯曲曲能能，Et是是键键的的二二面面角角扭扭转转能能，Enb是是非非键键作作用用能能，它它包包括括van der Waals作用能，偶极（电荷）作用能、氢键作用能等等。作用能，偶极（电荷）作用能、氢键作用能等等。第23页，讲稿共66张，创作于星期日分子力学v分子的空间能位位能能函函数数描描述述了了各各种种形形式式的的相相互互作作用用力力对对分分子子位位能能的的影影响响，它它的的有有关关参参数数、常常数数和和表表达达式式通通常常称称为为力场。力场。对对于于某某个个分分子子来来说说，空空间间能能是
17、是分分子子构构象象的的函函数数。由由于于在在分分子子内内部部的的作作用用力力比比较较复复杂杂，作作用用类类型型也也较较多；对于不同类型的体系作用力的情况也有差别。多；对于不同类型的体系作用力的情况也有差别。第24页，讲稿共66张，创作于星期日分子力学v分子的空间能第25页，讲稿共66张，创作于星期日分子力学v分子的力场形式-键能项：键长伸缩能键长伸缩能 Morse曲线：曲线：V=De 1-exp-a(r-r0)2谐振势函数：谐振势函数：V=(k/2)(r-r0)2第26页，讲稿共66张，创作于星期日分子力学v分子的力场形式-键能项：键长伸缩能键长伸缩能 含非谐项的函数：含非谐项的函数：V=(k
18、/2)(r-r0)21-k1(r-r0)-k2(r-r0)2-k3(r-r0)3第27页，讲稿共66张，创作于星期日分子力学v分子的力场形式-键能项：键角弯曲能键角弯曲能 V=(k/2)(-0)2V=(k/2)(-0)21-k1(-0)-k2(-0)2-k3(-0)3键长及键角交叉项：键长及键角交叉项：Vb/=(1/2)kr(r-r0)(-0)第28页，讲稿共66张，创作于星期日分子力学v分子的力场形式-键能项：扭转能扭转能 分子内部基团绕化学键的旋转会改变分子的构象，同分子内部基团绕化学键的旋转会改变分子的构象，同时也改变了分子的能量状态。用二面角时也改变了分子的能量状态。用二面角(i)来描
19、述这种来描述这种运动形式运动形式,采用采用Fourier级数形式表示为级数形式表示为：第29页，讲稿共66张，创作于星期日分子力学v分子的力场形式-键能项：扭转能扭转能 一一般般情情况况下下，由由于于二二面面角角的的扭扭转转对对总总能能量量的的贡贡献献小小于于键键长长和和键键角角的的贡贡献献，一一般般情情况况下下二二面面角角的的改改变变要要比比键键长长和和键键角角的的变变化化自自由由得得多多。因因此此在在一一些些处处理理大大分分子子的的力力场场中中常常保保持持键键长长、键键角角不不变变，只只考考虑虑二二面面角角及及其其他他的的作作用用而而优优化整个分子的构象和能量。化整个分子的构象和能量。第3
20、0页，讲稿共66张，创作于星期日分子力学v分子的力场形式-非键能项：van der Waals势势 V=Ar-12-Br-6 V=Aexp（B/r）-Cr-6 一般地说前者节省机时而后者却能更好地描述原子间一般地说前者节省机时而后者却能更好地描述原子间的非键作用。的非键作用。第31页，讲稿共66张，创作于星期日van der Waals L-J势势第32页，讲稿共66张，创作于星期日分子力学v分子的力场形式-非键能项：静电相互作用点电荷，部份电荷点电荷，部份电荷第33页，讲稿共66张，创作于星期日分子力学v分子的力场形式-非键能项：静电相互作用对于原子间的静电作用，作为一级近似，仅考虑永久对于
21、原子间的静电作用，作为一级近似，仅考虑永久偶极（偶极（i，j）间的作用）间的作用Vdipl：第34页，讲稿共66张，创作于星期日分子力学v分子的力场形式-氢键V(r)=A/r12-C/r10实例：实例：YETI力场力场VHB =(A/r12-C/r10)cos2 cos4 r第35页，讲稿共66张，创作于星期日Improper torsion第36页，讲稿共66张，创作于星期日分子力学v分子的力场形式-氢键:Charmm力场“第37页，讲稿共66张，创作于星期日分子力学v力场的参数化力场的参数化分子力学力场的性能即它的计算结果的准确性和可分子力学力场的性能即它的计算结果的准确性和可靠性主要取决
22、于势能函数和结构参数。靠性主要取决于势能函数和结构参数。这些有关力常数，结构参数的这些有关力常数，结构参数的“本征值本征值”的置定过的置定过程称为力场的参数化。程称为力场的参数化。参数化的过程要在大量的热力学、光谱学实验数据的参数化的过程要在大量的热力学、光谱学实验数据的基础上进行，有时也需要由量子化学计算的结果提供基础上进行，有时也需要由量子化学计算的结果提供数据。数据。第38页，讲稿共66张，创作于星期日分子力学v力场的参数化力场的参数化各各类类键键长长、键键角角的的“本本征征值值”一一般般取取自自晶晶体体学学、电电子子衍衍射射或或其其他他的的谱谱学学数数据据，键键伸伸缩缩和和角角变变力力
23、常常数数主主要要由由振振动动光光谱谱数数据据确确定定，扭转力常数经常要从分子内旋转位垒来推算。扭转力常数经常要从分子内旋转位垒来推算。对对于于不不同同的的力力场场不不仅仅力力场场参参数数不不同同，函函数数形形式式也也可可能能不不同同。因因此此，在将一个力场中的参数应用于另一个力场时应十分小心。在将一个力场中的参数应用于另一个力场时应十分小心。一一个个好好的的力力场场不不仅仅能能重重现现已已被被研研究究过过的的实实验验观观察察结结果果，而而且且能能有有一定的广泛性，能用于解决未被实验测定过的分子的结构和性质。一定的广泛性，能用于解决未被实验测定过的分子的结构和性质。第39页，讲稿共66张，创作于
24、星期日分子力学v力场的参数化力场的参数化-常见的力场及程序 QCFF/PI A Warshel&M levittMMI/MMPI NL AllingerMM2/MMP2 NL AllingerMM3NL AllingerECEPPHA ScheragaAMBERP KollmanCHARMMM KarplusGROMOSvan GunsterenSYBYLTripos Inc.DISCOVER MSI Inc.第40页，讲稿共66张，创作于星期日分子力学v力场所存在的问题两个相互作用原子间的诱导偶极的作用会受到其它原子的影响；两个相互作用原子间的诱导偶极的作用会受到其它原子的影响；非非键键作作
25、用用势势中中假假定定原原子子为为球球形形，实实际际上上非非键键作作用用受受原原子子形形状状影影响响，还需考虑孤对电子；还需考虑孤对电子；谐振势函数不能精确拟合实验数据谐振势函数不能精确拟合实验数据对于静电作用的处理过于简化。对于静电作用的处理过于简化。第41页，讲稿共66张，创作于星期日分子力学v力场的发展趋势考虑原子极化率考虑原子极化率取用高次项取用高次项发展含金属的力场发展含金属的力场第42页，讲稿共66张，创作于星期日分子力学v力场举例：Bond Stretching:Eb=k2b(b-b0)2+k3b(b-b0)3+k4b(b-b0)4Angle Bending:E =k2(-0)2+
26、k3(-0)3+k4(-0)4Torsion:E =k1(1-cos )+k2(1-cos2 )+k3(1-cos3 )Non-bond:En=qiqj/r+2(r*/r)9-3(r*/r)6第43页，讲稿共66张，创作于星期日分子力学v力场举例：Bond/Bond Coupling:Ebb=kbb(b-b0)(b-b0)Bond/Angle Coupling:Eb =kb(b-b0)(-0)Angle/Angle Coupling:E =k(-0)(-0)第44页，讲稿共66张，创作于星期日分子力学v力场举例：Angle/Angle/Torsion coupling:E =k (-0)(-0
27、)cos Bond/Torsion and Angle/Torsion coupling:E b=(b-b0)(k1 b cos +k2 b cos2 +k3 b cos3 )E =(-0)k1 cos +k 2 cos2 +k 3 cos3 第45页，讲稿共66张，创作于星期日分子力学v分子结构的优化分子结构的优化首先，给出所计算分子的试探结构。不一定是分子的稳首先，给出所计算分子的试探结构。不一定是分子的稳定构象，而且往往不是稳定构象。定构象，而且往往不是稳定构象。然后，将总空间能然后，将总空间能Es对所有描述分子构象的变量即分子各原子的对所有描述分子构象的变量即分子各原子的三维坐标在一定
28、的范围内求极小值。三维坐标在一定的范围内求极小值。由于数学上只能保证求得局部极小值，即实现局部优化，而不一由于数学上只能保证求得局部极小值，即实现局部优化，而不一定能求得全局最小值。所以得到的是在这一构象附近的一相对稳定能求得全局最小值。所以得到的是在这一构象附近的一相对稳定的构象。定的构象。分子力学常用的优化方法有使用一阶导数的最速下降法和分子力学常用的优化方法有使用一阶导数的最速下降法和使用二阶导数的使用二阶导数的Newton-Raphson法。法。第46页，讲稿共66张，创作于星期日分子力学v分子结构的优化分子结构的优化粗结构粗结构能量极小构象能量极小构象分子几何优化分子几何优化E局部极
29、小值问题；鞍点局部极小值问题；鞍点第47页，讲稿共66张，创作于星期日分子力学v分子结构的优化分子结构的优化 用用于于描描述述分分子子初初始始结结构构的的原原子子坐坐标标可可以以使使用用分分子子内内坐坐标标、直直角角坐坐标标或或晶晶体体坐坐标标。从从晶晶体体数数据据得得到到初初始始结结构构数数据据往往往往是是比比较较方方便便的的，若若没没有有晶晶体体数数据据，则可用模型来估计。则可用模型来估计。输入坐标及连接关系输入坐标及连接关系力场选择、作用项、参数力场选择、作用项、参数能量极小化能量极小化最终结构与能量最终结构与能量其它信息其它信息第48页，讲稿共66张，创作于星期日分子力学v分子结构的优
30、化分子结构的优化 除除了了初初始始坐坐标标外外，还还要要提提供供分分子子中中所所有有原原子子的的联联接接关关系系，以以便便自自动动搜搜索索任任何何两两个个原原子子之之间间的的作作用用，按按不不同同的的联联接接关关系系以以不不同同的的能能量量函函数数形形式式计计算算对对总总能能量量的的贡贡献献。计计算算中中所所用用的的能能量量参参数数大大部部分分已已在在程程序序中中准准备备好好，有有时时，要对某些参数进行修改或增补。要对某些参数进行修改或增补。输入坐标及连接关系输入坐标及连接关系力场选择、作用项、参数力场选择、作用项、参数能量极小化能量极小化最终结构与能量最终结构与能量其它信息其它信息第49页，
31、讲稿共66张，创作于星期日分子力学v分子结构的优化分子结构的优化 分子总能量是原子三维坐标的分子总能量是原子三维坐标的函数，在计算完初始构象的分子能函数，在计算完初始构象的分子能量后，要进行能量极小化的迭代，量后，要进行能量极小化的迭代，直到达到收敛标准为止。最终给出直到达到收敛标准为止。最终给出分子体系优化的原子坐标，总空间分子体系优化的原子坐标，总空间能及各能量项的贡献。能及各能量项的贡献。输入坐标及连接关系输入坐标及连接关系力场选择、作用项、参数力场选择、作用项、参数能量极小化能量极小化最终结构与能量最终结构与能量其它信息其它信息第50页，讲稿共66张，创作于星期日分子力学v分子结构的优
32、化分子结构的优化由于一般只是局部优化，这样的计算只能找到所用的初始构象附由于一般只是局部优化，这样的计算只能找到所用的初始构象附近的近的“最优构象最优构象”。所以，选择初始构象是非常关键的。所以，选择初始构象是非常关键的。若为了找到全局能量最低构象，须将所有可能的初始构若为了找到全局能量最低构象，须将所有可能的初始构象分别进行优化，最后进行比较，从而确定分子体系的象分别进行优化，最后进行比较，从而确定分子体系的最优构象。最优构象。对于较大的分子，可能的初始构象的数目会随原子数目的增加而对于较大的分子，可能的初始构象的数目会随原子数目的增加而急剧增加。在选择初始构象时，应把从基本的化学知识方面考
33、虑急剧增加。在选择初始构象时，应把从基本的化学知识方面考虑是不可能的构象略去。是不可能的构象略去。第51页，讲稿共66张，创作于星期日分子力学v能量极小化算法能量极小化算法一级微商算法一级微商算法最陡下降算法最陡下降算法 Steepest Descents-SD共轭梯度算法共轭梯度算法 Conjugate Gradients-CONJ二级微商算法二级微商算法Newton-Raphson Method 第52页，讲稿共66张，创作于星期日分子力学v能量极小化算法能量极小化算法-最陡下降法（SD）f(x i+1)=f(x i)+f(x i)x第53页，讲稿共66张，创作于星期日分子力学v能量极小化
34、算法比较能量极小化算法比较最陡下降法：方向变化大，收敛慢，优化辐度大方向变化大，收敛慢，优化辐度大共轭梯度法共轭梯度法收敛快，易陷入局部势阱，对初始结构偏离不大收敛快，易陷入局部势阱，对初始结构偏离不大 Newton-Raphson法法计算量较大，当微商小时收敛快计算量较大，当微商小时收敛快第54页，讲稿共66张，创作于星期日分子力学v分子力学的特点分子力学的特点概念清楚，便于理解及应用概念清楚，便于理解及应用概念简明易于接受。分子力学中的总概念简明易于接受。分子力学中的总“能量能量”被被分解成键的伸缩、键角弯曲、键的扭曲和非键分解成键的伸缩、键角弯曲、键的扭曲和非键作用等，比起量子化学计算中
35、的作用等，比起量子化学计算中的Fock矩阵等矩阵等概念来要直观易懂。概念来要直观易懂。第55页，讲稿共66张，创作于星期日分子力学v分子力学的特点分子力学的特点计算速度快计算速度快量子化学从头算的计算量随原子轨道数目的增量子化学从头算的计算量随原子轨道数目的增加，按加，按4次方的速度上升，而分子力学的计算量次方的速度上升，而分子力学的计算量仅与原子数目的平方成正比。仅与原子数目的平方成正比。计算时间计算时间 -MM正比于原子数正比于原子数m的平方的平方m2 QM正比于轨道数正比于轨道数n的的n4或或n3第56页，讲稿共66张，创作于星期日分子力学v分子力学的特点分子力学的特点与量子化学计算相辅
36、相成与量子化学计算相辅相成 分分子子力力学学是是一一种种经经验验方方法法，其其力力场场是是在在大大量量的的实实验验数数据据的的基基础础上上产产生生的的。分分子子力力学学宜宜用用于于对对大大分分子子进进行行构构象象分分析析、研研究究与与空空间间效效应应密密切切相相关关的的有有机机反反应应机机理理、反反应应活活性性、有有机机物物的的稳稳定定性性及及生生物物活活性性分分子子的的构构象象与与活活性性的的关关系系；但但是是，当当研研究究对对象象与与所所用用的的分分子子力力学学力力场场参参数数化化基基于于的的分分子子集集合合相相差差甚甚远远时时不不宜宜使使用用，当当然然也也不不能能用用于于人人们感兴趣但没
37、有足够多的实验数据的新类型的分子。们感兴趣但没有足够多的实验数据的新类型的分子。第57页，讲稿共66张，创作于星期日分子力学v分子力学的特点分子力学的特点与量子化学计算相辅相成与量子化学计算相辅相成 对对于于化化合合物物的的电电子子结结构构、光光谱谱性性质质、反反应应能能力力等等涉涉及及电电子子运运动动的的研研究究，则则应应使使用用量量子子化化学学计计算算的的方方法法。然然而而，在在许许多多情情况况下下，将将量量子子化化学学计计算算和和分分子子力力学学计计算算结结合合使使用用能能取取得得较较好好的的效效果果。分分子子力力学学计计算算结结果果可可提提供供量量子子化化学学计计算算所所需需的的分分子
38、子构构象象坐坐标标，而而量量子子化化学学计计算算结结果果又又给给出出了了分分子子力力学学所所不不能能给给出出的的分分子子的的电电子性质。子性质。第58页，讲稿共66张，创作于星期日分子力学v分子力学与量子化学计算比较分子力学与量子化学计算比较 分子力学与量子化学计算的区别分子力学与量子化学计算的区别分子力学是经典模型，以原子为分子力学是经典模型，以原子为“粒子粒子”，按经典力学运动，而，按经典力学运动，而量子化学则主要处理对象为电子，其运动服从量子力学规量子化学则主要处理对象为电子，其运动服从量子力学规律；律；量子化学中，电子或原子核间的相互作用服从库仑定律，而分量子化学中，电子或原子核间的相
39、互作用服从库仑定律，而分子力学中每对原子之间有一特定的作用势函数，原子不同或者子力学中每对原子之间有一特定的作用势函数，原子不同或者原子虽然相同但所处环境不同，则势函数不同，即使对同一对原子虽然相同但所处环境不同，则势函数不同，即使对同一对原子，也无法给出准确的普适势函数。原子，也无法给出准确的普适势函数。第59页，讲稿共66张，创作于星期日ApplicationslNMR or X-ray structure refinementlProtein structure predictionlProtein folding kinetics and mechanicslConformationa
40、l dynamicslGlobal optimizationlDNA/RNA simulationslMembrane proteins/lipid layers simulationslComputer aided drug design第60页，讲稿共66张，创作于星期日MD Simulations第61页，讲稿共66张，创作于星期日Ion permeation through alpha-hemolysin This movie shows induced by a 240mV bias currents of K+and Cl-ions through the transmembran
41、e pore of alpha-hemolysin第62页，讲稿共66张，创作于星期日Water channels in cell membraneswater permeation movie for the 2003 Nobel prize in Chemistry第63页，讲稿共66张，创作于星期日A potassium ion interacts with a 16 Angstroms carbon nanotube.第64页，讲稿共66张，创作于星期日 Argon cluster MD第65页，讲稿共66张，创作于星期日感谢大家观看12.10.2022第66页，讲稿共66张，创作于星期日