Á¿µÄʵÑéÊý¾ÝÖ®Í⣬²¢²ÎÕÕÁ˾«È·µÄÁ¿×Ó¼ÆËã½á¹û£¬ÓÈÆäÊÊÓÃÓÚÓлú·Ö×Óϵͳ¡£CFF91¡¢CFF95¡¢COMPASSÓëPCFF¶¼ÊôÓÚµÚ¶þ´úÁ¦³¡¡£
Ò»°ãÇé¿öÏ£¬·Ö×ÓµÄ×ÜÄÜÁ¿Îª¶¯ÄÜÓëÊÆÄܵĺͣ¬¶ø·Ö×ÓµÄÊÆÄܿɱíʾΪ¼òµ¥µÄ¼¸ºÎ×ø±êµÄº¯Êý¡£¸´ÔÓ·Ö×ÓµÄ×ÜÊÆÄÜÒ»°ã¿É·ÖΪ¸÷ÀàÐÍÊÆÄܵÄ×ܺͣ¬¼´ÓɳɼüÏUbonded£©¡¢½»²æÏî(Ucross)ºÍ·Ç¼üÏî(Unobonded)Èý²¿·Ö¹¹³É£¬¼´£º
×ÜÊÆÄÜ=³É¼üÏîÊÆÄÜ+½»²æÏîÊÆÄÜ+·Ç¼üÏîÊÆÄÜ Ï°¹ßÉϳ£ÓõÄÊýѧ±í´ïʽΪ£ºU=Ubonded+Ucross+Unobonded
³É¼üÏBonded term£©°üÀ¨¼üÉìËõÏî(Ub)¡¢¼ü½ÇÍäÇúÏî(U¦È)¡¢¶þÃæ½ÇŤÇúÏî(U¦Õ)ºÍÀëÆ½ÃæÕñ¶¯Ïî(U¦Ö)¡£
½»²æÄÜÁ¿Ïî(Cross energy term)Ö÷ÒªÖ¸³É¼ü×÷ÓÃÖ®¼äñîºÏÒýÆðµÄÄÜÁ¿±ä»¯£¬°üÀ¨¼üÉìËõ-¼üÉìËõ(Ubond-bond)¡¢¼üÉìËõ-¼ü½ÇÍäÇú-¼üÉìËõ(Ubond-angle-bond)¡¢¼ü½ÇÍäÇú-¼ü½ÇÍäÇú(Uangel-angle)¡¢¶þÃæ½ÇŤתÏî-¼üÉìËõ(Utorsion-bond)¡¢¶þÃæ½ÇŤתÏî-¼ü½ÇÍäÇú-¼ü½ÇÍäÇú(Utorsion-angle-angle)¡¢¼ü½ÇÍäÇú-¶þÃæ½ÇŤתÏî-¼ü½ÇÍäÇú(Uangle-torsion-angle)¡¢ÀëÆ½ÃæÕñ¶¯-¶þÃæ½ÇŤתÏî-ÀëÆ½ÃæÕñ¶¯(Uoop-torsion-oop)£¬Ò»°ãÇé¿öÏ£¬½»²æÄÜÁ¿Ïî¶ÔÕûÌåÊÆÄÜÓ°Ïì½ÏС,ËùÒÔÔÚÒ»°ãµÄÁ¦³¡Öв¢²»º¬ÓÐÕâЩÏî,Ö»ÓÐÔÚ¼ÆËã½á¹û¾«¶ÈÒªÇó½Ï¸ßʱ²ÅÓн»²æÄÜÁ¿Ïî¡£½»²æÏîͨ³£ÐÎʽ±È½Ï¸´ÔÓ,¼ÆËãÆðÀ´·Ç³£ºÄʱ,ÓÈÆäÊÇÔÚ¼ÆËãÔ×ÓÊýÄ¿½Ï¶àµÄÌåϵʱ¡£
·Ç¼üÏNobonded term£©°üÀ¨·¶µÂÍß¶ú˹×÷Óã¨Uvdw£©¡¢¿âÂ×¾²µç×÷Óã¨Uel£©ºÍÇâ¼ü×÷Óã¨UH£©µÈ¡£
ÒÔÏÂ×ÅÖØÇ¿µ÷¼¸¸öÖØÒªµÄÊÆÄÜÏ
(1) ·Ç¼ü½á·¶µÂÍß¶úË¹ÊÆÄÜ(Uvdw) ÈôA¡¢BÁ½Ô×ÓÊôÓÚͬһ·Ö×Óµ«Æä¼ä¸ô¶àÓÚÁ½¸öÁ¬½ÓµÄ»¯Ñ§¼ü£¨Èç A¡ªC¡ªC¡ªB£©£¬»òÕßÁ½Ô×ÓÊôÓÚ²»Í¬µÄ·Ö×Ó£¬ËäÎÞÖ±½ÓÁ¬½Ó³É¼üµ«ÒòΪ¾àÀë½Ï½ü£¬Ôò´ËÁ½Ô×Ó¼äÓзǼü½á·¶µÂÍß¶ú˹Á¦£¨Nobonded van der waals force£©¡£¼ÆËã·Ç¼ü½á×÷ÓÃʱ£¬Í¨³£½«Ä³¸öÔ×ÓÊÓΪλÓÚÆäÔ×ÓºË×ø±êµÄÒ»µã£¬Á¦³¡ÖÐ×î³£¼ûµÄ·Ç¼ü½áÊÆÄÜÐÎʽΪLennard-Jones(L-J)ÊÆÄÜ£¬ÕâÖÖÊÆÄÜÓÖ³ÆÎª12-6ÊÆÄÜ£¬ÆäÊýѧ±í´ïʽΪ£º
£¨2. ´íÎ󣡽öÖ÷Îĵµ¡££©
ʽÖЦÅÓë¦ÆÎªÊÆÄܲÎÊý£¨ÒòÔ×Ó¶Ô¶øÒ죩£¬rΪÔ×Ó¶Ô¼äµÄ¾àÀ룬r-12ÏîΪ
ÅųâÏr-6ÏîΪÎüÒýÏµ±rºÜ´óÊ±ÊÆÄÜÇ÷½üÓÚÁ㣬Õâ±íÃ÷Ô×Ó¶Ô¾àÀëºÜԶʱÒÑÎ޷Ǽü½á×÷Óã¬Èçͼ 2- ´íÎ󣡽öÖ÷Îĵµ¡£ L-JÊÆÄÜÇúÏߣ¬±í 2. ´íÎ󣡽öÖ÷Îĵµ¡£²¿·Ö³£¼û·Ö×ÓµÄLJÊÆÄܲÎÊý¡£
ͼ 2- ´íÎ󣡽öÖ÷Îĵµ¡£ L-JÊÆÄÜÇúÏß
±í 2. ´íÎ󣡽öÖ÷Îĵµ¡£²¿·Ö³£¼û·Ö×ÓµÄLJÊÆÄܲÎÊý
Ô×Ó»ò·Ö×Ó He Ar O2 CH4 C3H8
×¢: 1=0.1nm ¡£
£¨¦Å/kb£©/ K 10.41 141.6 126.3 161.3 268.5
r/ 2.602 3.350 3.382 3.721 4.992
(2)¼üÉìËõÊÆÄÜ(Ub) ÔÚ·Ö×ÓÖÐÏ໥¼ü½áµÄÔ×ÓÐγɻ¯Ñ§¼ü£¬ÈçC¡ªH£¬O¡ªH¼üµÈ£¬»¯Ñ§¼üµÄ¼ü³¤²¢·Çºã¶¨²»±ä£¬¶øÊÇÔÚÆäÆ½ºâÖµ¸½½ü³ÊС·ù¶ÈµÄÕñ¶¯£¬ÕâÖÖÕñ¶¯Ëù²úÉúµÄÊÆÄܳÆÎª¼üÉìËõÊÆÄÜ»ò¼üÉìËõÊÆÄÜÏBond stretching term£©¡£
¼üÉìËõÊÆÄÜÏîµÄÒ»°ãÐÎʽΪ(½«´ËÐÎʽÊÓΪ¼òгÕñ¶¯)£º
£¨2. ´íÎ󣡽öÖ÷Îĵµ¡££©
ʽÖÐkbΪ¼üÉìËõµÄµ¯Á¦³£Êý£¬µ¯Á¦³£ÊýÔ½´ó£¬Õñ¶¯ÆµÂÊÔ½¸ß£¬Õñ¶¯Ô½¿ì£»riºÍ·Ö±ð±íʾµÚi¸ö¼üµÄ¼ü³¤ÓëµÚi¸ö¼üµÄƽºâ¼ü³¤¡£
±í 2. ´íÎ󣡽öÖ÷Îĵµ¡£ MM2Á¦³¡Ò»Ð©¼üµÄÉìËõÊÆÄܲÎÊý
¼ü½á Csp3¡ªCsp3 Csp3¡ªCsp2 Csp2=Csp2 Csp3¡ªH Csp2=O Csp3¡ªNsp3
ƽºâ¼ü³¤r/ 1.523 1.497 1.337 1.113 1.438 1.438
µ¯Á¦³£Êýkb/[kcal/(mol¡¤)] 317 317 690 662 777 367
×¢£º1kcal=4.187 kJ, 1=0.1nm ¡£
ÓÐЩÁ¦³¡£¨ÈçMM3Á¦³¡£©ÎªÁËÔö¼Ó¼ÆËãµÄ׼ȷÐÔ£¬³ýÁ˼òгÕñ¶¯ÏîÍâÔö¼ÓÁ˷ǼòгÕñ¶¯Ïî(³ýÁ˶þ´ÎµÄ¼òгÕñ¶¯ÏîÍ⣬»¹Ôö¼ÓÁËÈý´ÎÓëËĴεķǼäЪÕñ¶¯Ïî)ÒÔ¸ü׼ȷµÄÃèÊö¼üÉìËõÊÆÄÜ£¬È磺
£¨2. ´íÎ󣡽öÖ÷Îĵµ¡££©
(3) ¼ü½ÇÍäÇúÊÆÄÜ(U¦È) Èô·Ö×ÓÖÐÁ¬Ðø¼ü½áµÄÈýÔ×ÓÐγÉÒ»¶¨µÄ¼ü½Ç¡£Óë¼üµÄÉìËõÒ»Ñù£¬ÕâЩ¼ü½Ç²¢·Çά³Ö²»±ä£¬¶øÊÇÔÚÆäÆ½ºâÖµ¸½½ü³ÊС·ù¶ÈµÄÕñµ´£¨Ò»°ãÊÓΪ¼òгÕñµ´£©£¬Í¨¹ý¼ü½ÇÍäÇúÏAngle bending term£©¿ÉÒÔÃèÊöÕâÏîÊÆÄÜÏ¼ü½ÇÍäÇúÏîµÄÒ»°ãÐÎʽΪ£º
£¨2. ´íÎ󣡽öÖ÷Îĵµ¡££©
ʽÖЦÈi¼°·Ö±ð±íʾµÚi¸ö¼üµÄ¼ü½Ç¼°Æäƽºâ¼ü½Ç£¬k¦ÈΪ¼ü½ÇÍäÇúµÄµ¯Á¦³£Êý¡£ÎªÁËÔö¼Ó¼ÆËãµÄ¾«È·ÐÔ£¬ÎÒÃÇ¿ÉÒÔÏñ¼üÉìËõÊÆÄÜÏîÒ»ÑùÔö¼Ó¸ß´Î·Ç¼òгÏîÒÔÌá¸ß¼ÆËãµÄ¾«È·¶È¡£
(4) ¶þÃæ½ÇŤÇúÊÆÄÜ(U¦Õ) Èô·Ö×ÓÖÐÁ¬Ðø¼ü½áËĸöÔ×Ó²¢ÐγɶþÃæ½Ç£¨Dihedral angle£©£¬Èç±ûÍéCH3CH2CH3ÖеÄH¡ªC¡ªC¡ªC¡¢H¡ªC¡ªC¡ªHµÈ£¬Ò»°ã·Ö×ÓÖеĶþÃæ½ÇÒ×ÓÚŤ¶¯£¬ËùÒÔÐèÒªÒ»¸öÊÆÄÜÏîÀ´ÃèÊöÕâÖÖÊÆÄÜ£¬ÎÒÃdzÆÖ®Îª¶þÃæ½ÇŤÇúÏTorsion angle term£©£¬¶þÃæ½ÇŤÇúÏîµÄÒ»°ãÐÎʽΪ£º
£¨2. ´íÎ󣡽öÖ÷Îĵµ¡££©
ʽÖУ¬¦ÕΪ¶þÃæ½Ç½Ç¶È£¬V1¡¢V2ºÍV3¾ùΪ¶þÃæ½ÇŤÇúÏîµÄµ¯Á¦³£Êý¡£ (5) ÀëÆ½ÃæÕñ¶¯ÊÆÄÜ(U¦Ö) ÔÚ·Ö×ÓÖУ¬ÓÐЩÔ×ÓÓй²Æ½ÃæµÄÇãÏò£¬Èç±ûͪ
CH3C(O)CH3ÖУ¬Ì¼Ô×ÓÓëÑõÔ×ӵį½ºâλÖÃλÓÚ¹²Í¬Æ½Ã棬ÆäËûÈç±½»·ÉϵÄ̼ÇâÔ×Ó¡¢Ï©ÀàÖеÄ̼ÇâÔ×ӵȾùÓй²Æ½ÃæµÄÇãÏò£¬Í¨³£¹²Æ½ÃæµÄËĸöÔ×ÓµÄÖÐÐÄÔ×Ó»áÀë¿ªÆ½ÃæÐ¡·ù¶ÈµÄÕñ¶¯£¬ËùÒÔÎÒÃÇʹÓÃÀëÆ½ÃæÕñ¶¯ÏOut-of-plane bending term£©À´ÃèÊöÕâÖÖÊÆÄÜ£¬ÀëÆ½ÃæÕñ¶¯ÏîµÄÒ»°ãÐÎʽΪ£º
(2. ´íÎ󣡽öÖ÷Îĵµ¡£)
ʽÖУ¬¦ÖΪÀë¿ªÆ½ÃæÕñ¶¯µÄ½Ç¶È£¬k¦ÖΪÀë¿ªÆ½ÃæÕñ¶¯ÏîµÄµ¯Á¦³£Êý¡£ (6) ¿âÂ×¾²µçÊÆÄÜ£¨Uel£© Àë×Ó»ò´øÓв¿·ÖµçºÉµÄ·Ö×Ó£¬ÔòÕâЩ´øµçºÉµÄÁ£×Ó¼äÒòÏà¶Ô¾àÀë´æÔÚמ²µçÎüÒý»òÅųâ×÷Óã¬ÎÒÃÇͨ¹ý¿âÂ××÷ÓÃÏColumbic interaction term£©À´ÃèÊöÕâÖÖ¾²µç×÷Ó㬿âÂ××÷ÓÃÏîµÄÒ»°ãÐÎʽΪ£º
(2. ´íÎ󣡽öÖ÷Îĵµ¡£)
ʽÖУ¬qiºÍqjΪµÚi¸öÁ£×ÓÓëµÚj¸öÁ£×ÓËù´øµÄµçºÉ£¬rijΪËù´øµçºÉÁ£×Ó¼äµÄ¾àÀ룬DΪÓÐЧ½éµç³£Êý£¨effective dielectric constant£©¡£
2.2.3 Lennard-Jones(L-J)×÷ÓÃÊÆÓëTersoff×÷ÓÃÊÆ
¶ÔÓÚµ¥±Ú̼ÄÉÃ׵ķÖ×Ó¶¯Á¦Ñ§Ä£Ä⣬ÓÉÓÚµ¥±Ú̼ÄÉÃ×¹ÜÄ£ÐÍÖмÈÓÐÁ½ÌåÏ໥×÷ÓÃÒ²ÓÐÈýÌåÏ໥×÷Óã¬ËùÒÔÕë¶Ô´ËÄ£ÄâÌåϵ£¬²ÉÓÃÁËÁ½ÖÖ×÷ÓÃÊÆ£¬·Ö±ðÊÇLennard-Jones×÷ÓÃÊÆÓëTersoff×÷ÓÃÊÆ¡£
Lennard-Jones(L-J)×÷ÓÃÊÆ£¨ÓÖ³ÆL-JÊÆÄܺ¯Êý»ò6-12ÊÆÄܺ¯Êý£©ÊǼÆË㻯ѧÖÐÓÃÀ´Ä£ÄâÁ½·Ö×Ó¼ä×÷ÓÃÊÆÄܵÄÒ»¸öº¯Êý£¬ÊÇÒ»ÖÖÁ½Ìå×÷ÓÃÊÆ¡£×îÔçÓÉÂü³¹Ë¹ÌØ´óѧµÄÊýѧ¼ÒJohn Lennard-JonesÓÚ1931ÄêÌá³ö¡£ÓÉÓÚÆä½âÎöÐÎʽ¼òµ¥ËùÒÔ±»¹ã·ºÊ¹Óã¬L-J×÷ÓÃÊÆµÄ±í´ïʽΪ£º
L-J×÷ÓÃÊÆÏàÓ¦µÄÁ½Ìå×÷ÓÃÁ¦Îª£º
£¨2. ´íÎ󣡽öÖ÷Îĵµ¡££©
ʽÖЦÅΪL-J×÷ÓÃÊÆµÄÊÆÄÜÚåµÄ²ÎÊý£¬¦ÆÊÇL-J×÷ÓÃÊÆµÄ»¥Ïà×÷ÓÃÊÆÄÜÕýºÃΪÁãʱµÄÁ½Ìå¾àÀ루Ҳ³ÆÎªÆ½ºâ³£Êý£©,rijΪÁ£×Ói¡¢j¼äµÄ¾àÀë¡£
Tersoff×÷ÓÃÊÆÊÇÒ»ÖÖµäÐ͵ÄÈýÌå×÷ÓÃÊÆ£¬ËüÒѱ»¹ã·ºÓÃÓÚÃèÊö̼×廯ºÏÎïµÄÔ×Ó¼äµÄ×÷ÓÃÁ¦[i]¡£Tersoff×÷ÓÃÊÆµÄ¹«Ê½ÈçÏ£º
ÆäÖУº
£¨2. ´íÎ󣡽öÖ÷Îĵµ¡££©
£¨2. ´íÎ󣡽öÖ÷Îĵµ¡££©