·Ö×ÓÄ£ÄâµÄһЩ»ù±¾×ÊÁÏ(3)

2019-04-09 13:29

Á¿µÄʵÑéÊý¾ÝÖ®Í⣬²¢²ÎÕÕÁ˾«È·µÄÁ¿×Ó¼ÆËã½á¹û£¬ÓÈÆäÊÊÓÃÓÚÓлú·Ö×Óϵͳ¡£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. ´íÎ󣡽öÖ÷Îĵµ¡££©


·Ö×ÓÄ£ÄâµÄһЩ»ù±¾×ÊÁÏ(3).doc ½«±¾ÎĵÄWordÎĵµÏÂÔØµ½µçÄÔ ÏÂÔØÊ§°Ü»òÕßÎĵµ²»ÍêÕû£¬ÇëÁªÏµ¿Í·þÈËÔ±½â¾ö£¡

ÏÂһƪ£º¡°Ê®ÈýÎ塱֨µãÏîÄ¿-¾Û°±õ¥ÅÝÄ­ËÜÁÏÖÆÆ·Éú²ú½¨ÉèÏîÄ¿½ÚÄÜÆÀ¹À±¨

Ïà¹ØÔĶÁ
±¾ÀàÅÅÐÐ
¡Á ×¢²á»áÔ±Ãâ·ÑÏÂÔØ£¨ÏÂÔØºó¿ÉÒÔ×ÔÓɸ´ÖƺÍÅŰ棩

ÂíÉÏ×¢²á»áÔ±

×¢£ºÏÂÔØÎĵµÓпÉÄÜ¡°Ö»ÓÐĿ¼»òÕßÄÚÈݲ»È«¡±µÈÇé¿ö£¬ÇëÏÂÔØÖ®Ç°×¢Òâ±æ±ð£¬Èç¹ûÄúÒѸ¶·ÑÇÒÎÞ·¨ÏÂÔØ»òÄÚÈÝÓÐÎÊÌ⣬ÇëÁªÏµÎÒÃÇЭÖúÄã´¦Àí¡£
΢ÐÅ£º QQ£º