µ±Ç°Î»Ö㺠Ê×Ò³ > Á¿×Ó»¯Ñ§ >´óÌåϵÈõÏ໥×÷ÓüÆËãµÄ½â¾öÖ®µÀ

´óÌåϵÈõÏ໥×÷ÓüÆËãµÄ½â¾öÖ®µÀ

×÷Õß sobereva
À´Ô´: Сľ³æ 2100 42 ¾Ù±¨Ìû×Ó
+¹Ø×¢

´óÌåϵÈõÏ໥×÷ÓüÆËãµÄ½â¾öÖ®µÀ
ÎÄ/Sobereva   2013-Dec-18



1 Ç°ÑÔ

ÒÔDFT-DΪ´ú±íµÄÒ»´óÀàɫɢУÕý·½·¨Ê¹µÃ´óÌåϵÈõÏ໥×÷ÓõļÆËã³ÉΪÁË¿ÉÄÜ£¬½éÉܺÍʹÓüû¡¶DFT-DɫɢУÕýµÄʹÓá·£¨http://hi.baidu.com/sobereva/item/772d38c6e13165d497445243£©¡£±ÊÕßÔøÔÚ¡¶ÂÒ̸DFT-D¡·£¨http://hi.baidu.com/sobereva/item/4e18b2da4f023317e1f46fde£©Ò»ÎĵÄ×îºóÖ¸³öÁ˲»Í¬¼ÆËãÄÜÁ¦ÏÂ×îÍƼöµÄ¼ÆËã·½·¨£¬µ«ÕâЩÍƼöÖ»ÊÇÀíÂÛÉϵġ£¶øʵ¼Ê¼ÆËãÖУ¬Ó°ÏìºÄʱµÄÒòËؼ«¶à£¬³ýÁËÀíÂÛÉϵļÆËãÁ¿Ö®Í⣬»¹Ö±½Ó¾ö¶¨ÓÚ¼ÆËã³ÌÐòµÄÑ¡Ôñ£¬ÒÔ¼°ÊýÖµÉϵļÓËÙ·½·¨¡£ÔÚ±¾ÎÄÖУ¬±ÊÕß½«´ÓÍêȫʵÓÃÐԵĽǶȳö·¢£¬ÌÖÂÛÔÚÆÕͨ¼ÆËãÌõ¼þÏ£¬ÈçºÎÑ¡Ôñ×î¼ÑµÄ¼ÆË㼶±ð׼ȷÑо¿´óÌåϵµÄÈõÏ໥×÷ÓÃÎÊÌâ¡£

Ïë׼ȷ¼ÆËã´óÌåϵÈõÏ໥×÷Óã¬ÓÖ²»ÏëäĿµØɵƴ¼ÆËãÁ¿£¬±ØÐëÐèÒª¡°·¨±¦¡±¡£ÏÂÃæÒ»½Ú½«»á¼òÒª½éÉÜÕâÑù7¸ö·¨±¦£¬·Ö±ðÊÇ£ºPM7¡¢MOZYME¡¢MOPAC³ÌÐò¡¢DFT-D¡¢gCP¡¢ÃܶÈÄâºÏ½üËÆ&COSX£¬ÒÔ¼°ORCA³ÌÐò¡£È»ºóÔÚµÚ3½Úͨ¹ý×ÛºÏÔËÓÃÕâЩ·¨±¦£¬¸ø³ö¶ÔÓÚ²»Í¬ÊýÄ¿Ô­×ÓµÄÌåϵ×îÊʺÏʹÓõļÆËã·½·¨¡£ÔÚµÚ4½Ú£¬½«»áͨ¹ý¸÷ÖÖʵÀýչʾʵ¼Ê¼ÆËãºÄʱ²¢½øÐÐÌÖÂÛ¡£×îºóÒ»½Ú×ܽ᱾ÎÄ¡£

Ê×ÏÈÌرðÖµµÃÒ»ÌáµÄÊÇ£¬¿¿×îÖ÷Á÷×î³£ÓõÄGaussianÕâÑùµÄ³ÌÐò£¬¼´±ãÓÃÉÏÁËDFT-D£¬Ïë׼ȷ´¦Àí´óÌåϵ£¬ÄÇÒ²Ö»ÄÜһζµØ»¨·ÑºÜ¸ßµÄ´ú¼ÛɵË㣬¾ÍËã½µµÍµ½6-31+G**ÕâÑù¶ÔÓÚÈõÏ໥×÷ÓüÆËãµÍµÃ²»ÄÜÔٵ͵Ļù×飬ÈôÎÞ¶Ô³ÆÐÔ£¨´óµÄÈõÏ໥×÷ÓÃÌåϵ»ù±¾¶¼Ã»ÓжԳÆÐÔ£©ÔÚÆÕͨ¼ÆËãÌõ¼þÏ´¦ÀíÉÏ°Ù¸öÔ­×ÓµÄÌåϵ¶¼Í¦À§ÄÑ£¬´ËʱÈõÏ໥×÷ÓüÆËã׼ȷ¶ÈÒ²±È½Ï²î¡£È»¶ø£¬Èç¹û¸ÄÓÃORCAµÄ»°£¬½«»á¸Ð¾õORCA¼òÖ±¿ªÁ˹ң¬·¢ÏÖÔø¾­½¨Á¢ÔÚÄÔÖеġ°¼ÆË㾫¶È-¼ÆËãºÄʱ¡±·½ÃæµÄ¾­Ñé±»ÍêÈ«µß¸²£¬·Â·ð´ò¿ªÁËÒ»ÉÈͨÍùÐÂÊÀ½çµÄ´óÃÅ£¡


2 Æ߸ö·¨±¦

2.1 PM7°ë¾­Ñé·½·¨

¶ÔÓÚ¾Þ´óµÄÌåϵ£¬DFTÎÞÂÛ½áºÏʲô·¨×Ó£¬ÊÇÎÞÂÛÈçºÎÒ²´¦Àí²»¶¯µÄ£¨µ±È»£¬SCC-DFTBÕâÑùµÄ°ë¾­Ñé²ãÃæÉϵÄDFTÊÇÁíÍâÒ»»ØÊ£©¡£Èç¹ûÓÃÁ¦³¡£¬×¼È·¶ÈÍùÍùºÜÓÐÏÞ£¬ÆÕÊÊÐԲ¶øÇÒ¾­³£È±·¦ºÏÊʵIJÎÊý¡£°ë¾­Ñé·½·¨ÍùÍùÊÇΨһµÄÑ¡Ôñ£¬¾¡¹ÜÆäÖØÒªÐÔÒѾ­ºÜ´ó³Ì¶ÈÉϱ»ÈËÃǺöÊÓºÍÒÅÍü¡£

´«Í³µÄ°ë¾­Ñé·½·¨£¬ÖîÈçÖ÷Á÷µÄAM1¡¢PM3£¬ÒÔ¼°ºóÀ´ÖØвÎÊý»¯¶ø³ÉµÄRM1ºÍPM6£¬¶ÔÓÚÈõÏ໥×÷Óö¼±È½Ï²î¾¢£¬¾ÍÁ¬Çâ¼üÕâÑù´øÓÐÖ÷Òª¾²µç³É·ÝµÄÈõÏ໥×÷ÓÃÉв»ºÃ£¬¾Í¸ü±ðÌá·¶µÂ»ª×÷ÓÃÁË¡£×Ô´ÓɫɢУÕýÔÚ´ÓÍ·Ëã¡¢DFTÁìÓòÁ÷ÐпªÀ´ºó£¬ÈËÃÇ·×·×°Ñ°ë¾­Ñé·½·¨Ò²½áºÏÉÏÉ«É¢ÒÔ¼°ÆäËü·Ö×Ó¼ä×÷ÓõÄУÕýÏʹµÃ°ë¾­Ñé·½·¨¶ÔÈõÏ໥×÷ÓõıíÏÖÄÜÁ¦ÓÐÁËÖʵķÉÔ¾¡£°üÀ¨AM1-D¡¢PM3-D¡¢OM1/2/3-D¡¢PM6-DH1¡¢PM6-DH2¡¢PM6-DH2X¡¢PM6-DH+¡¢PM6-D3H4ºÍPM7¡£ÆäÖÐ2007ÄêÓÉMcNamaraºÍHillierÌá³öµÄAM1-DºÍPM3-DÓõÄÈ˱ȽÏÉÙ¡£2008ÄêÌá³öµÄOMx-DϵÁÐÖ»ÓÐÆä×÷ÕßWalter ThielµÄMNDO2005³ÌÐò²ÅÄÜË㣬ÕⶫÎ÷²»¹«¿ª¸ø´ó¼ÒÓ㨱ÊÕß±ÉÒÄÕâÖÖÐÐΪ£©£¬ËùÒÔ×ÔȻҲÁ÷Ðв»ÆðÀ´£¬¾Ý˵OM2-DµÄÈõÏ໥×÷ÓüÆËã׼ȷ¶ÈÄÜ´ïµ½DFT-DµÄË®×¼¡£2009ÄêÌá³öµÄPM6-DH1ÊÇÔÚPM6»ù´¡ÉϼÓÉÏÁËÉ«É¢¡¢Çâ¼üУÕýÏºóÀ´ºÜ¿ì±»2010ÄêÌá³öµÄ¸Ä½ø°æPM6-DH2ËùÈ¡´ú¡£2011ÄêµÄPM6-DH2XÔòÊÇÔÚPM6-DH2»ù´¡ÉÏÓÖÔö¼ÓÁ˶Ô±-ÑõºÍ±-µªÐ͵ı¼üУÕý¡£2012ÄêÓÖÌá³öÁËPM6-D3H4£¬³ýÁËÇâ¼ü¡¢É«É¢Ð£ÕýÏîÍ⻹ÒýÈëÁËPauli»¥³âУÕýÏ¿ÉÄÜÕâÊÇÄ¿Ç°×îºÃµÄ»ùÓÚPM6¶ÔÈõÏ໥×÷ÓýøÐÐУÕýµÄ·½·¨¡£Ç°Êö¼¸ÖÖÐÞÕýPM6µÄ·½·¨¶¼ÊÇHobzaµÈÈËÌá³öµÄ£¬2010ÄêÁíÒ»¸öÈËKorthÒ²Ìá³öÁËÒ»ÖÖPM6¼ÓÉÏÇâ¼üºÍɫɢУÕýµÄÐÎʽ£¬ÃûΪPM6-DH+£¬±ÈPM6-DH2ºÃ£¬µ«Î´±Ø±ÈPM6-D3H4ºÃ£¬ÁíÍâ´Ë·½·¨»¹´æÔÚÊÆÄÜÃæÔÚһЩ¹ý³ÌÖбíÏֵò»Æ½»¬µÄÎÊÌâ¡£PM7ÊÇPM3¡¢PM6ºÍMOPACµÄ×÷ÕßJ.J.P. StewartÓÚ2013ÄêÌá³öµÄ×îеİ뾭Ñé·½·¨£¬¸Ä½øÁËÆäÉÏÒ»´úPM6´æÔÚµÄһЩȱÏÝ£¬´Óͳ¼ÆÊý¾ÝÉÏ¿´ÈÈÁ¦Ñ§ÐÔÖʼÆË㲢ûʲôÏÔÖøÌá¸ß£¬µ«ÊÇÕâ´ÎÔÚѵÁ·²ÎÊýµÄ¹ý³ÌÖÐ×ÅÖصؿ¼ÂÇÁËÈõÏ໥×÷ÓÃÌåϵ£¬²¢ÇÒÒýÈëÁËÈõÏ໥×÷ÓõÄУÕýÏʹµÃÈõÏ໥×÷ÓüÆËã±ÈÆðPM6ÓÐÁ˳¤×ã½ø²½¡£PM7ÓÉÓڱȽÏУ¬Ä¿Ç°À´ËµºÍPM6-D3H4¡¢PM6-DH+Ïà±È½ÏË­¶ÔÓÚÈõÏ໥×÷Óøü׼ȷ»¹²»ºÃ˵£¬µ«Ó¦¸ÃÔÚ²®ÖÙÖ®¼ä¡£Ä¿Ç°À´Ëµ£¬PM7ÊÇÕûÌåÉÏ×îºÃµÄ°ë¾­Ñé·½·¨£¬¾Þ´óÌåϵÈõÏ໥×÷ÓüÆËã¸öÈË×îÍƼöPM7¡£²»¹ý¼´±ã´Óͳ¼ÆÊý¾ÝÉÏ¿´£¬PM7ÒÔ¼°¸÷ÖÖÐÞÕýPM6µÄ·½·¨ÔÚÈõÏ໥×÷ÓüÆËãÎó²îÉ϶¼ºÜС£¬ºÍDFT-DϵÁÐÏà±È»¹ÊÇÓвî¾à£¬¿É¿¿ÐÔ²î¾à¸ü´ó£¬ËùÒÔ¶ÔÓÚDFT-DÄÜËãµÃ¶¯µÄÌåϵ×îºÃÏÈÓÃDFT-D£¬µ«ÊÇ°ë¾­Ñé·½·¨¿ÉÒÔ×÷Ϊ³õË㣬±ÈÈç¹À¼ÆһϴóÖµÄ×÷ÓÃÄÜ¡¢³õ²½ÓÅ»¯Ò»Ï½ṹ£¬»òÕß×ö¶¯Á¦Ñ§¡¢¹¹ÏóËÑË÷Ö®ÀàµÄºÜºÄʱµÄÈÎÎñ¡£


2.2 MOZYME¼¼Êõ

MOZYMEÊÇJ.J.P.StewartÌá³öµÄÒ»ÖÖ»ùÓÚ¶¨Óò»¯¹ìµÀ¼Ó¿ì°ë¾­Ñé¼ÆËãµÄ¼¼Êõ£¬¾ßÌåϸ½ÚÔÚÕâÀﲻ̸£¬¿ÉÒÔ²ÎÔÄÆäÔ­ÎÄIJQC,58,133(1996)¡£²»ÏñÐí¶à»ùÓÚ¶¨Óò»¯¹ìµÀµÄ·½·¨£¬Ò²¾ÍÊǽöÏÞÓÚ×÷Õß×Ô¼º·¢paperɹɹÊý¾Ý¶øÎÞ·¨ÕæÕýͶÈë¹ã·ºÓ¦Óã¬MOZYMEÊÇÓо޴óʵÓÃÒâÒåµÄ£¬¼ÓËÙЧ¹ûÊǼ«ÆäÏÔÖøµÄ£¬Õ⽫ÔÚºóÃæµÄ²âÊÔÖп´µ½¡£¶øÇÒMOZYMEÊǸöºÚÏä·½·¨£¬²¢²»ÐèÒªÓû§¶îÍâÈ¥¿¼ÂÇʲô¡£MOZYMEµÄ´æÔÚʹµÃ°ë¾­ÑéÄÜ·½·¨´¦ÀíµÄÌåϵ³ß¶È´óÁËÁ½¸öÊýÁ¿¼¶ÒÔÉÏ£¬ÄÜ´ïµ½¼¸Íò¸öÔ­×ÓµÄÌåϵ¡£¶øÇÒMOZYME²»¹âÄÜÓÃÓÚµ¥µã¼ÆË㣬»¹ÄÜÓÃÓÚ¼¸ºÎÓÅ»¯¡¢ËÑË÷¹ý¶É̬¡¢¶¯Á¦Ñ§µÈÐèÒªµ¼ÊýµÄ³¡ºÏ£¬²»¹ýÖ»ÄÜÓÃÓڱտDzã»ù̬¡£¶ÔÓÚ´óÌåϵÈõÏ໥×÷ÓÃÎÊÌ⣬MOZYME¼¼ÊõÒýÈëµÄÎó²îͨ³£ÊǺÜСµÄ£¨ÆðÂë¶ÔÓÚÌåϵµÄ³ß¶È¶øÑÔ£©¡£µ«¶ÔÓÚ¼«ÉÙÊýµç×ӽṹÎÞ·¨Í¨¹ýLewisʽ±íÏÖµÄÌåϵ£¬ÓпÉÄÜSCF¹ý³Ì»áÊÕÁ²µ½¼¤·¢Ì¬É϶ø²»ÕýÈ·£¬²»¹ýÕâ¸öÎÊÌâ¶ÔÓÚÉúÎï·Ö×ÓÕâÀà³É¼ü·½Ê½µäÐ͵ÄÌåϵ¿Ï¶¨²»»á³öÏÖ¡£MOPACµÄºóÆÚµÄÉÌÒµ°æ±¾ÒÔ¼°Ãâ·Ñ°æ±¾ÖдÓ2009°æ¿ªÊ¼¶¼Ö§³ÖMOZYME¡£MOZYMEÒ²½öÔÚMOPACÖÐʵÏÖ£¬ÊÇMOPACµÄÒ»¸öÖØÒªÁÁµã¡£


2.3 MOPAC

MOPACÊÇ×îÖøÃûµÄרÃÅ×ö°ë¾­Ñé¼ÆËãµÄ³ÌÐò¡£µ½7.0°æ±¾Ò»Ö±ÊÇ¿ªÔ´Ãâ·ÑµÄ£¬ÔÚ90Äê´úMOPAC³ÉΪFujitsu¹«Ë¾µÄÉÌÒµÈí¼þ¡£ºóÀ´MOPACµÄ×÷ÕßStewartÀ뿪ÁËFujitsu×Ô¼º³ÉÁ¢Á˹«Ë¾²¢·¢²¼ÁËMOPAC2007£¬È»ºóÓÖ³öÁËMOPAC2009£¬×îеÄÊÇMOPAC2012£¬Ëü¶ÔѧÊõÓû§ÊÇÃâ·ÑµÄµ«²»¿ªÔ´¡£MOPAC2012Ö»ÊÇÒ»¸ö´óµÄ°æ±¾ºÅ£¬¾ßÌå°æ±¾ºÅÔÚ²»¶ÏƵ·±Éý¼¶¡£×îа汾ÏÂÔصØÖ·¼ûhttp://openmopac.net£¬windows°æÖ»ÓÐ5Õ׶࣬license¿ÉÒÔÃâ·ÑË÷È¡¡£MOPAC2012ÒѾ­ÊǸö·Ç³£³ÉÊìÎȶ¨¶øÇÒ¹¦ÄÜÏ൱ȫÃæµÄ³ÌÐò£¬ÊÖ²á׫дҲºÜÏêϸ£¬¼ûhttp://openmopac.net/manual/index.html¡£×îаæµÄMOPAC2012Ö§³ÖµÄÄܽϺô¦ÀíÈõÏ໥×÷Óõķ½·¨°üÀ¨PM6-DH2¡¢PM6-DH2X¡¢PM6-DH+¡¢PM6-D3£¨PM6½áºÏDFT-D3УÕý£©¡¢PM7¡£ËäÈ»ÓкܶàÁ¿×Ó»¯Ñ§³ÌÐò¶¼ÄÜ×ö°ë¾­Ñé¼ÆË㣬°üÀ¨GaussianÔÚÄÚ£¬µ«ÊÇÄÜÓð뾭Ñé·½·¨×öµÄÈÎÎñÀàÐÍÔ¶²»¼°MOPAC£¨ÖîÈçÖÜÆÚÐÔ¼ÆË㣩£¬¶øÇÒ¶Ôеİ뾭Ñé·½·¨µÄÖ§³Ö³Ì¶ÈÔ¶²»ÈçMOPAC£¬ÖîÈçGaussianÖÐÄÜÓõÄ×îºÃµÄ°ë¾­Ñé·½·¨Ö»ÓÐPM6£¬Éв»Ö§³ÖÈκο¼ÂÇÈõÏ໥×÷ÓÃУÕýµÄ°ë¾­Ñé·½·¨¡£MOPACµÄMOZYME¼¼Êõ¸üÊÇÆäËü³ÌÐòËùûÓеġ£MOPAC2012ÖаÑPM7ºÍMOZYME×éºÏʹÓÿÉνÊǶԸ¶¾Þ´óÈõÏ໥×÷ÓÃÌåϵµÄÒ»¸öÀûÆ÷¡£

Windows°æMOPAC2012°²×°·½·¨ÊÇ£º°ÑÏÂÔص½µÄ³ÌÐòѹËõ°ü½âѹµ½C:\Program FilesÀ²¢ÇÒ°ÑĿ¼Ãû¸ÄΪMOPAC¡£½øÈë´ËĿ¼£¬½¨Á¢Ò»¸öÎı¾Îļþ£¬Ãû×ÖΪÄãÏòStewartË÷È¡µÄÃÜÂ룬ºó׺Ϊ.mop£¬ÀýÈç1591535a16518252.mop¡£È»ºóÓÒ¼üµãËüÑ¡¡°´ò¿ª¡±£¬²¢ÓÃMOPAC2012.exeÀ´´ò¿ª¡£ÔÚMOPAC2012µÄ³ÌÐò´°¿ÚÀD»îGraphic1´°¿Ú£¬Çûسµ£¬ÊäÈëYes£¬È»ºóÑ¡¡°ÊÇ¡±£¬¼¤»î¼´Íê³É¡£

MOPAC2012µÄʹÓñȽϼòµ¥£¬µäÐ͵ÄÊäÈëÎļþÀý×ÓÈçÏ£º
PM7 CHARGE=-10 MOZYME
molecule
All coordinates are Cartesian
O       4.71271832 0   9.26971072 0  -7.55066164 0
C       4.13662696 0   9.48826045 0  -8.82968461 0
C       2.75025629 0   8.84701245 0  -8.99475234 0
...£¨Ê£ÓàµÄ×ø±ê£©
µÚÒ»ÐÐÊäÈëµÄÊǹؼü´Ê¡£ÕâÀïÓÃPM7£¬µçºÉΪ-10£¬²¢ÇÒÆôÓÃMOZYME¼¼Êõ£¨Ä¬Èϲ»ÆôÓã©¡£µÚ¶þ¡¢ÈýÐÐÊÇÓû§ËæÒⶨÒåµÄÎÄ×Ö¡£×ø±ê²¿·Öµ¥Î»Ä¬ÈÏΪ°£¡£XYZ×ø±êºóÃæµÄ0´ú±íÕâ¸öÔ­×Ó×ø±ê²»×öÓÅ»¯£¬Èç¹ûËùÓÐ×ø±êºóÃ涼ÊÇ0¾Í´ú±íÕâÊǸöµ¥µãÈÎÎñ£¬Èç¹ûÓеÄ×ø±êºóÃæÊÇ1£¬¾Í´ú±íµ±Ç°ÈÎÎñÊÇÓÅ»¯£¬²¢ÇÒÓÅ»¯ÕâЩд1µÄ×ø±ê¡£Èç¹ûд1SCF¹Ø¼ü´ÊÔòÎÞÂÛÊÇ·ñÓеÄ×ø±êºóÃæÊÇ1£¬¶¼Ö»×öµ¥µã¼ÆËã¡£

ÔÚWindowsÏ£¬½«ÊäÈëÎļþ±£´æΪ.mopΪºó׺£¬È»ºóË«»÷Ö®¾Í¿ªÊ¼ÔËË㣬²¢ÔÚµ±Ç°Ä¿Â¼Ïµõ½.outÊä³öÎļþºÍ.arcµµ°¸Îļþ£¨Ö»¼Ç¼¹Ø¼üÐÔÐÅÏ¢£©¡£ÏÖÓеÄMOPAC½á¹ûµÄ¿ÉÊÓ»¯¹¤¾ß¶¼²»ºÜºÃÓ㬶ÔÓÚ´óÌåϵµÄÓÅ»¯£¬ÍùÍùÏë¿´ÓÅ»¯ÈÎÎñµÄ¹ì¼£À´¿¼²ìÓÅ»¯ÈÎÎñÊÇ·ñºÏÀí£¬Õâ¿ÉÒÔÓñÊÕß¿ª·¢µÄmopac2xyz¹¤¾ß£¨http://hi.baidu.com/sobereva/item/02416ec7f5676010505058e2£©½áºÏVMDʵÏÖ¡£

×¢Òâ°ë¾­Ñé·½·¨¸ø³öµÄÄÜÁ¿²»ÊǾø¶ÔÄÜÁ¿£¬¶øÊÇÉú³ÉìÊ£¬ËùÒÔÊýÁ¿¼¶ÉϺʹÓÍ·Ëã³ÌÐò¸ø³öµÄÌåϵÄÜÁ¿Ïà²îÉõÔ¶¡£µ«ÊÇÇóÏà¶ÔÄÜÁ¿Ê±ºÍƽ³£Ò»Ñù¡£±ÈÈçÇó·Ö×Ó¼äÏ໥×÷ÓÃÄÜ£¬Ö±½Ó°Ñ¸´ºÏÎïµÄÉú³ÉìÊÓëÿ¸öƬ¶ÎµÄÉú³ÉìÊÏà¼õ¼´¿É¡£


2.4 DFT-D

DFT-DÒѾ­ÔÚÉÏÎÄÌáµ½µÄ¡¶DFT-DɫɢУÕýµÄʹÓ᷺͡¶ÂÒ̸DFT-D¡·Öз´¸´½éÉÜÁË£¬ÕâÀï¾Í²»ÔÙÌáÁË¡£


2.5 gCP

gCPÊÇÌá³öDFT-DµÄGrimmeÔÚJCP,136,154101(2012)Ìá³öµÄÒ»ÖֺܼòÒ׵ĻñµÃBSSEУÕýÄܵķ½·¨¡£±ÊÕßµ±³õ¾õµÃÕâ¸ö·½·¨Ã»Ê²Ã´Òâ˼£¬Óеãºú¸ã£¬µ«ºóÀ´·¢ÏÖ£¬gCP¶ÔÓÚ°ïÖú¼¸°Ù¸öÔ­×ӵĴóÌåϵÈõÏ໥×÷ÓüÆËãʵÔÚÊÇÓÐÖØÒªÒâÒå¡£ÕâÀïÏÈ̸һÏÂBSSEÎÊÌâ¡£

BSSE£¨»ù×éÖصþÎó²î£©ÎÊÌâÔÚ¡¶Ì¸Ì¸BSSEУÕýÓëGaussian¶ÔËüµÄ´¦Àí¡·£¨http://hi.baidu.com/sobereva/item/e1c1874c603580adde2a9fdc£©µ±ÖÐÓÐÏêϸÌÖÂÛ¡£BSSE±¾ÖÊÉÏÓÉÓÚ»ù×é²»Í걸ËùÒýÆ𣬻ᵼÖÂÈõÏ໥×÷ÓüÆËã½á¹ûÏÔÖø²»×¼£¨°ë¾­Ñé¡¢SAPT·½Ê½¼ÆËãÎÞBSSEÎÊÌ⣩¡£ºÜ´ó³Ì¶ÈÉϽâ¾öBSSEÎÊÌâµÄ·½·¨Ò»ÖÖÊÇʹÓú¬Óн϶àÃÖÉ¢º¯ÊýµÄ»ù×飬»òÕßÖÁÉÙÊÇ·ÖÁѼÛÊý½Ï¶àµÄ»ù×飬ÀýÈç4-zeta»ù×éµÄBSSE¾ÍºÜСÁË£¬ÒòΪ×îÍâ²ã»ùº¯ÊýÒѾ­ÌåÏÖ³öÁË°ëÃÖÉ¢ÌØÕ÷£»ÁíÒ»ÖÖÊÇʹÓÃCounterpoise (CP)·½·¨½øÐÐУÕý¡£µ«ÕâÁ½ÖÖ×ö·¨¶¼Ê¹¼ÆËãÁ¿¾çÔö¡£ÁíÍâÃÖÉ¢º¯Êý½Ï¶àʱ»¹»áµ¼ÖÂSCFÊÕÁ²À§ÄÑ£¬¶øʹÓÃCPʱÔòû·¨»ñµÃ½âÎöÌݶȣ¬Ê¹µÃ¼¸ºÎÓÅ»¯ºÄʱ¼«¸ß£¬¶øÇÒCPͨ³£Ö»Äܽâ¾ö·Ö×Ó¼äBSSEÎÊÌ⣬·Ö×ÓÄÚBSSEÎÊÌâºÜÄÑ¿¿CP½â¾ö¡£

Ö®ËùÒÔ²»º¬ÃÖÉ¢º¯ÊýµÄС»ù×éÔÚ¼ÆËãÈõÏ໥×÷ÓÃÉϽá¹û²î£¬Æäʵ¹Ø¼üÒòËؾÍÔÚÓÚBSSE¡£ÆÕ±éÌᳫÔÚ¼ÆËãÈõÏ໥×÷ÓÃʱ¼ÓÉÏÃÖÉ¢º¯Êý£¬Æäʵ´ÓʵЧÉϽ²¾ÍÊÇÔÚÓÚ¼õÇáBSSE¡£»»¾ä»°Ëµ£¬¶ÔÓÚ·Ö×Ó¼äÏ໥×÷ÓõļÆË㣬ÃÖÉ¢º¯ÊýÆðµ½µÄ×÷ÓúÍCPʵ¼ÊÉÏÊÇÏà·Â·ðµÄ¡£Õâµã´ÓJCP,134,084107(2011)µÄͼ4ÉÏÌåÏֵ÷dz£Çå³þ£¬cc-pVDZ+CPºÍaug-cc-pVTZµÄƽ¾ùÎó²îºÜÏà½ü¡£ËùÒÔ˵£¬Èç¹ûÓÃCPÀ´½â¾öµôBSSEÎÊÌ⣬ÄÇô¶ÔÓÚ·Ö×Ó¼äÏ໥×÷ÓüÆËã¾Í²»Ì«ÐèÒªÃÖÉ¢º¯ÊýÁË£¨ËäÈ»ÔÙ¼ÓÉÏÃÖÉ¢µÄ»°½á¹û»¹ÄÜÔÙ¸ÄÉÆÒ»µãµã£©¡£

gCPµÄÈ«³ÆÊÇgeometrical counterpoise£¬ÕâÊÇ»ùÓÚÔ­×Ó×ø±êºÍÔ¤ÏÈÄâºÏµÄ²ÎÊý½øÐмòµ¥¾­ÑéУÕýÀ´»ñµÃBSSEУÕýÄܵķ½·¨£¬Ð£Õýº¯ÊýÊ®·Ö¼òµ¥£¬¾ÍÏñDFT-DУÕýÒ»Ñù²»ÐèÒª¶îÍâºÄʱ¡£¶øÇÒgCP¿ÉÒÔÖ±½Ó½â¾ö·Ö×ÓÄÚBSSEÎÊÌ⣬Òò´Ë±ÈÆðCP»¹¸üÓÐÓÅÊÆ¡£²»¹ýgCP±Ï¾¹Ö»ÊǽüËÆ·½·¨£¬ËüµÃµ½µÄBSSEУÕýÄܵÄ׼ȷ¶ÈºÍÖ±½Ó×öCPÏà±È»¹ÊÇ»áÓÐÒ»¶¨Æ«²îµÄ¡£gCPµÄ²ÎÊýÒ»·½Ãæ¾ö¶¨ÓÚ»ù×飬һ·½Ãæ¾ö¶¨ÓÚÀíÂÛ·½·¨¡£gCPÒѾ­Ïò6-31G*¡¢def2-SVP¡¢def2-TZVPÕâÑùÖ÷Á÷µÄË«¡¢Èýzeta»ù×éºÍ¼«Ð¡»ùMINIS£¨scaled MINI¼«Ð¡»ù£©µÈ»ù×éÄâºÏÁ˲ÎÊý¡£ÀíÂÛ·½·¨Ö»·ÖΪÁ½ÖÖ£¬HFºÍDFT£¬²»Í¬µÄ·ºº¯¶ÔÓ¦µÄgCP²ÎÊýÏà²îºÜС£¬ËùÒÔgCP¶Ô¸÷ÖÖ·ºº¯Ê¹ÓÃÁËͳһµÄ²ÎÊý¡£

gCPµÄʵ¼ÊЧ¹ûºÜÃ÷ÏÔ¡£±ÈÈçÓÃdef2-SVPÕâÖÖ²»º¬ÃÖÉ¢º¯ÊýµÄË«·ÖÁѼۻù×é¼ÆËãÈõÏ໥×÷ÓÃÔ­±¾ÊDz»¿ÉÄܵģ¬¾«¶È²»¿É½ÓÊÜ£¬È»¶ø¼ÓÉÏgCPºóÈõÏ໥×÷ÓüÆËã½á¹û¾ÍºÍ²Î¿¼±ê×¼def2-QZVP½ÏΪ½Ó½üÁË¡£¶ødef2-TZVP£¬ÒѾ­±Èdef2-SVPÇ¿²»ÉÙ£¬¼ÓÉÏgCPºó½á¹û»¹ÄܽøÒ»²½¸ÄÉÆ£¬¶øÇÒ±ÈÆðdef2-SVP½áºÏgCPʱ¸ü½Ó½üdef2-QZVP¡£

gCP¿ÉÒÔ½áºÏDFT-D3УÕýÒ»ÆðʹÓã¬ÍêȫûÓÐì¶Ü£¬ÀýÈ綼¼Óµ½BLYPÉÏ¿ÉÒÔ¹¹³ÉBLYP-gCP-D3¡£×¢ÒâÓÐЩÈËÎóÒÔΪDFT-D²ÎÊýµ±ÖоÍÒѾ­°ÑBSSEÎÊÌâ¸øµÈЧ½â¾öÁË£¬Êµ¼ÊÉϸù±¾²»ÊÇÕâô»ØÊ£¡DFT-D3ÔÚÄâºÏ²ÎÊýµÄʱºòÊÇÔÚdef2-QZVPϽøÐеģ¬´Ë»ù×éËäÈ»»¹ÁôÓÐÒ»µãµãBSSEÎÊÌ⣬µ«ÎÊÌâÒѾ­ºÜСÁË£¬ËùÒÔDFT-D3µÄ²ÎÊý¼¸ºõûÓеÈЧ½â¾öBSSEÎÊÌ⣬¶øÇÒDFT-D3µÄº¯ÊýÐÎʽ±¾ÉíҲû·¨ÏñgCPÕâÑù½â¾öBSSEÎÊÌâ¡£Òò´Ë£¬DFT-D3ÓÃÀ´½â¾ö´«Í³·ºº¯ÈõÏ໥×÷ÓÃÃèÊö²î£¬gCP½â¾ö»ù×é²ãÃæÉϵÄÎÊÌ⣬gCPºÍDFT-D3Ò»ÆðÓþͰÑÀíÂۺͻù×éµÄÎÊÌⶼ½â¾öÁË£¬Ê¹µÃ±ãÒË·½·¨+±ãÒË»ù×éϵÄÈõÏ໥×÷ÓüÆË㾫¶ÈÒÀÈ»ºÜ²»´í£¬ºÜÊʺ϶Ը¶´óÌåϵ¡£

HFÒÔ¼°²»Í¬DFT·½·¨½áºÏgCPºÍDFT-D3ºóÈõÏ໥×÷ÓüÆË㾫¶È¸÷²»Ïàͬ¡£BLYP-gCP-D3²»ÊÇ×îºÃµ«Ò²Í¦²»´í£¬Ç¿ÁÒÍƼöʹÓá£B3LYP-gCP-D3±ÈËüûÃ÷ÏÔÓÅÊÆ£¬¶øÇÒ¼ÆËãÁ¿¸ü´ó£¬²»½¨ÒéÓá£ËäÈ»M06-2X±¾Éí¼ÆËãÈõÏ໥×÷Óúܲ»´í£¬µ«ÊÇgCPºÍËü½áºÏЧ¹û²¢²»ºÃ¡£ÁîÈËÒâÍâµÄÊÇ£¬HF-gCP-D3½áºÏMINIS¼«Ð¡»ù¾ÓÈ»¼ÆËã½á¹ûÏ൱²»´í£¬¹À¼ÆÊÇÇ¡ºÃÎó²î±È½ÏºÃµØµÖÏûÁË¡£´óµÖÊÇÓÉÓÚÕâ¸ö·¢ÏÖ£¬GrimmeÔÚJCC,34,1672(2013)Ìá³öÒ»¸öËùνµÄHF-3c·½·¨£¬Ï൱ÓÚHF-gCP-D3-SRB/MINIX£¬SRBÊǶîÍâÒýÈëµÄ½â¾ö¶Ì³Ì»ù×é²»Í걸ÐÔµÄУÕýÏMINIXÊÇÒ»ÖÖ½ÏÇáÔ­×ÓÓü«Ð¡»ù¶ø½ÏÖØÔ­×ÓÓÃË«·ÖÁѼۻùµÄרΪHF-3cÌØÖƵĻù×é¡£HF-3c½á¹û²»´í£¬ÒѾ­±»ORCAÖ§³Ö¡£µ«ÊDZÊÕ߱ȽϱÉÒÄHF-3c·½·¨£¬ÎïÀíÒâÒå²î£¬½á¹ûºÃÖ÷ÒªÊÇÎó²îµÖÏûµÃºÃ£¬¶ÔÓÚÌØÊâÇé¿öû׼ºÜÀ㬲»ËãºÜ¿¿Æ×£¬¶øÇÒHF±ÈÆð´¿·ºº¯¼ÆËãÁ¿´óµÃ¶à£¨ÌرðÊÇʹÓÃÃܶÈÄâºÏ½üËÆʱ£¬¼ûºóÎÄ£©¡£

gCPµÄ¼ÆËã³ÌÐò¿ÉÒÔ´ÓÕâÀïÏÂÔØ£ºhttp://www.thch.uni-bonn.de/tc/downloads/gcp/index.html¡£
±àÒ뷽ʽ£º¼ÙÉèÓõÄÊÇifort£¬¾Í°ÑMakefileÀïµÄ# FC = ifortÈ¥µô×¢ÊÍ£¬²¢°Ñ-staticÈ¥µô£¬È»ºómake¼´¿É¡£
µäÐ͵ÄÔËÐз½Ê½±ÈÈ磺./gcp water.xyz -level dft/svp£¬Ë²¼ä¼´µÃµ½gCPУÕýÄÜ¡£water.xyzÊÇÌåϵ½á¹¹µÄxyzÎļþ£¬-levelºóÃæ½ÓÀíÂÛ·½·¨ºÍ»ù×飬ÕâÀï¼´¼ÙÉè¼ÆË㼶±ðΪDFT·½·¨½áºÏdef2-SVP£¬Ïê¼ûÊֲᡣ¼ÓÉÏ-grad²ÎÊý¿ÉÒÔÊä³öУÕýÄܵÄÌݶȣ¬¼ÓÉÏ-hess¿ÉÒÔÊä³öËüµÄHessian¡£


2.6 ÃܶÈÄâºÏ½üËƺÍCOSX

ÃܶÈÄâºÏ½üËÆ(Density-fitting approximation)ÊÇÒ»ÖÖÊýÖµÉϵĽüËƼ¼Êõ£¬¿ÉÒÔ´ó´ó½µµÍ´¦ÀíË«µç×Ó»ý·ÖÉϵļÆËãÏûºÄ£¬Í¬Ê±ÒýÈëµÄÎó²î¿ÉÒÔ»ù±¾ºöÂÔ²»¼Æ¡£ÃܶÈÄâºÏ½üËÆϵÄË«µç×Ó»ý·ÖµÄ±í´ïʽ´ÓÐÎʽÉÏÀàËÆÓÚ£¨µ«²»Ïàͬ£©²åÈëÁËRI(Resolution-of-the-identity£¬¼´1=¡Æ|m><m|)£¬ËùÒÔÃܶÈÄâºÏ½üËÆÔÚһЩÎÄÏ׺ͳÌÐòÖб»³ÆΪRI½üËÆ¡£

ÃܶÈÄâºÏ½üËƵĹؼüÊÇ°ÑÔ­ÏȵĻùº¯Êý¦Õ£¨Ò²½Ð¹ìµÀ»ùº¯Êý£©µÄ³Ë»ýͨ¹ý¸¨Öú»ùº¯Êý¦ÖµÄÏßÐÔ×éºÏÀ´½üËƱí´ï£¬¼´
¦Õ_p*¦Õ_q¡Ö¡Æ[k]c_pq,k*¦Ö_k
×éºÏϵÊýcÊÇͨ¹ýÄâºÏµÃµ½µÄ£¬ÕâÊÇÃܶÈÄâºÏ½üËÆÃû×ÖµÄÓÉÀ´¡£×¢ÒâÕâ²¢·ÇÊǶÔÓÚµç×ÓÃܶÈÀ´ÄâºÏ£¬¶øÊǶÔÓÚ¹ãÒåµÄÃܶȣ¬¼´»ùº¯ÊýµÄ³Ë»ýÀ´ÄâºÏ¡£ÄâºÏ·½Ê½µÄ·½·¨ºÜ¶à£¬Ò»°ãÓõÄÊÇͨ¹ý×îС»¯¿âÂØ»ý·ÖÎó²îÀ´ÊµÏÖ£¬¼´×îС»¯Õâ¸öË«µç×Ó»ý·Ö
(¦Ñ_pq-¦Ñ'_pq|1/r12|¦Ñ_pq-¦Ñ'_pq)£¬ÆäÖЦÑ_pq=¦Õ_p*¦Õ_qÊǾ«È·µÄ³Ë»ý£¬¦Ñ'_pq=¡Æ[k]c_pq,k*¦Ö_k´ú±íÄâºÏµÄ³Ë»ý¡£

²ÉÓÃÉÏÃæÕâÖÖÄâºÏ·½Ê½Ê±£¬(pq|1/r12|rs)ÕâÑùµÄË«µç×Ó-ËÄÖ¸±ê»ý·Ö¾Í±»Ö±½Óת»¯ÎªÁËͨ¹ý(pq|1/r12|k)ÕâÖÖË«µç×Ó-ÈýÖ¸±êµÄ»ý·ÖµÄ×éºÏÀ´±í´ï£¬»ý·Ö±¾ÉíµÄ¼ÆËãÒò´Ë¸ü¼òµ¥ÁË¡£¶ÔÓÚHF¡¢DFTÖеĿâÂز¿·ÖµÄ¼ÆË㣬ÕâÖÖ¼¼ÊõÒ²±ÜÃâÁË´¢´æºÍ´¦ÀíËÄÖ¸±ê»ý·Ö£¬¶øÖ»ÐèÒª´¢´æºÍ´¦ÀíË«¡¢ÈýÖ¸±êµÄÊý¾Ý£¬¹ÊÔÚ¹¹½¨HF/KSËã·ûµÄ¹ý³ÌÉÏҲʡÁË´óÁ¿Ê±¼ä£¬Ï¸½Ú¿ÉÒԲο´ORCA3.0ÊÖ²á6.2.2.3½Ú¡£¾¡¹ÜÈýÖ¸±ê»ý·Ö->ËÄÖ¸±ê»ý·Ö±ä»»¹ý³Ì±¾ÉíÐèÒªÒ»µãʱ¼ä£¬µ«ÕâºÍÃܶÈÄâºÏ¼¼Êõ´øÀ´µÄÀûÒæÏà±ÈÊÇ΢²»×ãµÀµÄ¡£ÓÉÓÚ´¿·ºº¯£¨¼´·ÇÔÓ»¯·ºº¯£©²»º¬HF½»»»²¿·Ö£¬Ö»º¬¿âÂز¿·Ö(J)£¬Òò´ËÃܶÈÄâºÏ½üËÆÁî´¿·ºº¯µÄ¼ÆËãËÙ¶ÈÓÐÁËÖÊ·ÉÔ¾£¬¿ìÁËÒ»¸öÊýÁ¿¼¶ÒÔÉÏ£¬¶øÇÒ´¢´æ¿Õ¼äÏûºÄÒ²´óΪ½µµÍ£¬Õâ±»³ÆΪRI-J¼¼Êõ¡£RI-J¼¼ÊõÒѾ­¹ã·ºÆÕ¼°£¬ÖîÈçGaussian¡¢Turbmole¡¢ORCA¡¢Molpro¡¢PSI¡¢Q-ChemµÈÖ÷Á÷³ÌÐò¶¼ÒѾ­Ö§³ÖËü¡£

ÃܶÈÄâºÏ¼¼ÊõÒ²¿ÉÒÔÓÃÔÚ¼ÓËÙHF½»»»²¿·Ö(K)µÄ¼ÆË㣬³ÆΪRI-K£¬µ«ÓÉÓÚ»ùÓÚÈýÖ¸±ê»ý·Ö¹¹½¨KµÄ¹ý³Ì±È¹¹½¨JÂé·³¡¢ºÄʱµÃ¶à£¨¿É²Î¼ûPCCP,4,4285(2002)£©£¬ËùÒÔ´øÀ´µÄ¼ÓËÙЧ¹ûÃ÷ÏÔ²»ÈçÓÃÔÚ¿âÂز¿·ÖÄÇô´ó£¬RI-KҲûÓÐRI-JÄÇôÁ÷ÐУ¬ÏñGaussian¾Í²»Ö§³Ö¡£Èç¹ûÃܶÈÄâºÏ¼ÈÓÃÔÚJÒ²ÓÃÔÚKÉÏÃ棬¾Í½Ð×öRI-JK¡£¶ÔÓÚÔÓ»¯·ºº¯»òHF£¬ËäÈ»RI-JK»á´øÀ´ÏÔÖøµÄËÙ¶ÈÌáÉý£¬µ«ÊǼÆËãËÙ¶ÈÔ¶²»¼°´¿·ºº¯½áºÏRI-J¡£ÓÉÓÚRI-K·½Ê½¼ÆËãKµÄËٶȺÍÕ¼¾Ý¹ìµÀÊýÄ¿´æÔÚÕý±È¹Øϵ£¬ËùÒÔËùÓûù×éÔ½´ó£¬RI-K´øÀ´µÄ¼ÓËÙ¾ÍÔ½ÏÔÖø¡£¶ø¶ÔÓÚ´óÌåϵ£¨Õ¼¾Ý¹ìµÀÊýºÜ¶à£©Ê¹ÓÃС»ù×éµÄʱºòRI-K¾ÍÆð²»µ½¶àÉÙ¼ÓËÙ×÷Óá£

ÃܶÈÄâºÏ¼¼ÊõµÄÁíÒ»¸öÖØÒªµÄÓ¦ÓÃÊǼÓËÙMP2£¬Õâ³ÆΪRI-MP2¡£MP2ÈÎÎñÔÚ΢ÈżÆË㲿·ÖÖ÷ÒªµÄʱ¼ä¶¼»¨ÔÚÁËË«µç×Ó»ý·Ö±ä»»¹ý³Ì£¬ÆäÖÐ×îºÄʱµÄÒ»²½±ä»»µÄ¸´ÔÓ¶ÈÊÇO(N_occ*N_AO^4)£¬ÆäÖÐN_occÊÇÕ¼¾Ý¹ìµÀÊýÄ¿£¬N_AOÊÇ»ùº¯ÊýÊýÄ¿¡£¶øʹÓÃÃܶÈÄâºÏ½üËƺó£¬×îºÄʱµÄ±ä»»³ÉΪÁËO(N_occ^2*N_vir^2*N_aux)£¬ÆäÖÐN_auxÊǸ¨Öú»ùº¯ÊýµÄÊýÄ¿£¬Õâ±ÈÆðÔ­ÏÈO(N_occ*N_AO^4)¼ÆËãÁ¿´ó´ó½µµÍÁË£¬ÒòΪN_AO=N_occ+N_vir£¬N_auxËäÈ»´óÓÚN_AOµ«»¹ÊÇÔÚÒ»¸öÊýÁ¿¼¶Éϵġ£Õû¸öMP2µÄ¼ÆËã¹ý³ÌÖУ¬¶ÔÓÚ´ó»ù×éÀ´Ëµ£¬Êµ¼ÊÉÏÍùÍù¶àÊýʱ¼ä²»ÊÇ»¨ÔÚÁË΢ÈżÆËãÉÏ£¬¶øÊÇ»¨ÔÚÁËHFµü´ú¹ý³ÌÉÏ¡£ÎªÁËÕûÌå¼ÆËã¸ü¿ì£¬RI-MP2ÈÎÎñÖпÉÒÔÔÙ½áºÏRI-JKÀ´¼ÓËÙHF²¿·ÖµÄ¼ÆËã¡£ÏÖÔÚºÜÁ÷ÐеÄË«ÔÓ»¯·ºº¯ÊÇÔÓ»¯·ºº¯µÄSCFÄÜÁ¿Óë»ùÓÚÏàÓ¦DFT¹ìµÀµÄMP2ÄÜÁ¿µÄÏßÐÔ×éºÏ£¬Ò²¿ÉÒÔͬÑùµØͬʱÓÃÃܶÈÄâºÏ¼¼Êõ¼ÓËÙDFTºÍ΢ÈÅÕâÁ½²¿·ÖµÄ¼ÆËã¡£ÃܶÈÄâºÏ¼¼ÊõÔÚÁ¿×Ó»¯Ñ§ÉÏ»¹ÓкܶàÖØÒªÓ¦Óã¬ÓÉÓںͱ¾ÎÄÖ÷ÌâÎÞ¹ØÕâÀï¾Í²»ÔÙ¶àÌáÁË¡£

ÔÚORCAÖÐÓиöRIJONX£¬Ò²½ÐRIJDX£¬Êµ¼ÊÉÏÕâ¾ÍÊÇÖ¸¶ÔÓÚÔÓ»¯·ºº¯/HFÖ»ÓÃRI-J¼ÓËÙ¿âÂز¿·Ö£¬¶ø½»»»²¿·ÖÕÕ³£¼ÆËã¡£ONXµÄÈ«³ÆÊÇHFX calculated with O(N)£¬Æ京ÒåÊÇHF½»»»²¿·ÖÓÉÓÚËæ¾àÀëË¥¼õµÃ±È½Ï¿ì£¬Òò¶øÉæ¼°½ÏÔ¶¾àÀëʱ¿ÉÒÔ¶Ô»ý·Ö½øÐÐÆÁ±Î£¬¶ÔÓÚ´óÌåϵKµÄ¼ÆËãÁ¿ÊÇO(N)£¬¼´Ëæ×Å»ùº¯ÊýÔö¼Ó¼ÆËãÁ¿ÊÇÏßÐÔÔö¼ÓµÄ£¬Õâʵ¼ÊÉÏÊDZȽÏÀíÏ뻯µÄÇé¿ö¡£¶ÔÓÚ»ù×é½Ï´óµÄʱºò£¬¹âÊÇÖС¢Ð¡Ìåϵ¾ÍÒѾ­Ëã²»¶¯ÁË£¬¼´±ã¶ÔÓںܴóÌåϵÓÐÏßÐÔ±ê¶ÈÒâÒåÒ²²»´ó¡£

ÓëRI-K²»Í¬µÄ¼Ó¿ì½»»»²¿·Ö¼ÆËãµÄ½üËÆ·½·¨ÊÇCOSX£¬ÕâÊÇORCAµÄÍ·Ä¿Frank NeeseÔÚChem. Phys.,356,98(2009)ÖÐÌá³öµÄ£¬Ò²ÊÇORCAµÄÒ»¸öÖØÒªÁÁµã¡£COSXÊôÓÚ°ëÊýÖµ»ý·Ö·½·¨£¬ËüÔÚ¼ÆËãKµÄʱºòÓÖ×ö½âÎö»ý·ÖÓÖ×ö»ùÓÚÊýÖµ¸ñµãµÄ»ý·Ö¡£Ïà±ÈRI-K£¬COSXÉó¤´¦Àí´óÌåϵ£¬ÌåϵԽ´ó±ÈRI-KµÄÓÅÊƾÍÓúÃ÷ÏÔ¡£Í¨³£À´Ëµ£¬COSXÔÚСÌåϵ£¨Ô¼30¸öÔ­×ÓÒÔÏ£©Ã»RI-K¿ì£¬40×óÓÒ¸öÔ­×Óʱ¶þÕßËÙ¶ÈÏà·Â·ð£¬100¸öÔ­×ÓʱCOSXµÄºÄʱ¾ÍÍêÈ«ºÍRI-K²»ÔÚÒ»¸ö¾³½çÁË£¬¼¸°Ù¸öÔ­×ÓʱÄǾÍÖ»ÄÜÓÃCOSX²ÅËãµÃ¶¯ÁË¡£COSX±ÈRI-KµÄºÃ´¦»¹ÔÚÓÚCOSXµÄÄÚ´æÏûºÄҪСºÜ¶à£¬ÁíÍâ¶ÔÓÚ¿ª¿Ç²ãÌåϵRI-KµÄËٶȱȱտDzãʱºÄʱ¼Ó±¶£¬¶øCOSXÔòûÕâ¸öÎÊÌâ¡£Ò»°ã´óСµÄ¸ñµãÏÂCOSXÒýÈëµÄÎó²îºÍRI-KÏà·Â·ð£¬¶¼ºÜС¶ø»ù±¾¿ÉÒÔºöÂÔ¡£¶ÔÓÚÔÓ»¯·ºº¯µÄ¼ÆËã±¾À´HFÔÓ»¯²¿·Ö¾ÍÒýÈëµÃ±È½ÏÓÐÏÞ£¬Òò´ËRI-KºÍCOSX´øÀ´µÄÎó²î¾Í¸üСÁË¡£ÈôÏë½øÒ»²½¼õСCOSXµÄÎó²î£¬Ö»ÒªÔö¼Ó»ý·Ö¸ñµãµÄÊýÄ¿¾Í¿ÉÒÔ£¨ÔÚORCAÀïÓÃGridxÉ裬±ÈÈçÓÃGridx4ÉõÖÁGridx5£©£¬µ±È»ºÄʱҲ»áÏàÓ¦Ôö¼Ó¡£RI-K±ÈCOSXÓÐÒ»¸öºÃ´¦ÊÇÎó²î×ÜΪÕý£¬¶øCOSXµÄÎó²îʱÕýʱ¸º£¬Òò´ËÇóÏà¶ÔÄÜÁ¿Ê±RI-KµÄÎó²îµÖÏûµÃ¿ÉÄÜ»áÉÔ΢±ÈCOSXºÃÒ»µã¡£COSXºÍRI-J½áºÏʹÓ÷ֱð¼ÓËÙ½»»»ºÍ¿âÂز¿·Ö¾Í½Ð×öRIJCOSX¡£¶ÔÓÚСÌåϵ´ó»ù×é¼ÆËãÒ»°ã½¨ÒéÓÃRIJKÀ´¼ÓËÙ£¬¶ÔÓÚ´óÌåϵÔò½¨ÒéÓÃRIJCOSX¡£¶Ô¼¸Ê®¸öÔ­×ÓµÄÌåϵͨ³£RIJKºÍRIJCOSX¶Ô¾ø¶ÔÄÜÁ¿ÒýÈëµÄÎó²îÃ÷ÏÔСÓÚ1kcal/mol£¬Èç¹ûÊÇÇóÏà¶ÔÄÜÁ¿£¬Îó²î»¹ÄܵÖÏûµô²»ÉÙ¡£

ÃܶÈÄâºÏ¼ÆËã±ØÐëµÃÓи¨Öú»ù×é(Auxiliary basis sets£¬Ò²³ÆÃܶÈÄâºÏ»ù×飩²ÅÐС£¸¨Öú»ù×éÈçͬһ°ãÒâÒåÉϵĻù×飬Ëü¶¨ÒåÁËÿÖÖÔªËØÔ­×ÓËùЯ´øµÄ¸¨Öú»ùº¯ÊýµÄÊýÄ¿ºÍÌØÕ÷¡£Ò»°ãÓõĸ¨Öú»ù×éÓõĶ¼ÊÇÆÕͨµÄÔ­×ÓÖÐÐĵĸß˹º¯Êý£¬¼È°üº¬·ÇÊÕËõµÄÒ²°üº¬ÊÕËõµÄ£¬ËùÒÔºÍÒ»°ãµÄ»ù×éÔÚÐÎʽÉÏÎÞÒì¡£ÏÔÈ»Èç¹û¸¨Öú»ù×éÊÇÍ걸¼¯£¬ÄÇôÃܶÈÄâºÏ½üËƾͲ»ÔÙÊǽüËÆ£¬¶øÊÇÍêÈ«¾«È·µÄÁË£¬µ±È»Õâ²»¿ÉÄÜ×öµ½¡£Êµ¼ÊËùÓõĸ¨Öú»ù×éµÄ³ß´ç±ØÈ»ÊÇÓÐÏ޵ġ£¼È²»ÄÜ̫С£¬Ì«Ð¡ÁËû·¨½Ï׼ȷ»¹Ô­³ö¹ìµÀ»ùº¯ÊýµÄ³Ë»ý£¬½á¹û²»×¼È·£»¶øÌ«´óµÄ»°¼ÆËãÁ¿ÓÖÌ«¸ß¡£ËùÒÔʵ¼ÊËùÓõÄÃܶÈÄâºÏ»ù×é±ØÐë´óСºÏÊÊ£¬²¢ÇÒ×îºÃÊÇרÃÅΪÌض¨»ù×éÓÅ»¯¹ýµÄ¡£AhlrichsϵÁлù×éºÍDunningÏà¹ØÒ»ÖÂÐÔ»ù×éÒÔ¼°ÆäËü¸ö±ð»ù×鶼רÃÅÓÐÈËÓÅ»¯³öÁËÆ¥ÅäµÄ¸¨Öú»ù×é¡£ÕâЩ¸¨Öú»ù×é¸ù¾ÝÓô¦Ö÷Òª·ÖΪÈýÀ࣬һÖÖÊÇ´øÓÐ/Jºó׺µÄ(Coulomb fitting)£¬×¨ÃÅÓÃÓÚRI-J»òÕßRIJCOSX£»Ò»ÖÖÊÇ/JKºó׺µÄ(Coulomb+Exchange fitting)£¬×¨ÃÅÓÃÓÚRI-JK£»Ò»ÖÖÊÇ/Cºó׺µÄ(Correlation fitting)£¬×¨ÃÅÓÃÓÚRI-MP2¡£ÀýÈ磬def2-TZVPP¾ÍÓÐdef2-TZVPP/J¡¢def2-TZVPP/JKºÍdef2-TZVPP/C¡£/J¡¢/JKºÍ/CÐ͸¨Öú»ù×éÔÚÓÅ»¯GTFÖ¸ÊýºÍÈ·¶¨ÊÕËõ·½Ê½Ê±µÄÄ¿±ê·Ö±ðÊÇÁîÃܶÈÄâºÏ½üËÆÔÚÌåϵµÄ¿âÂØÄÜ¡¢½»»»ÄܺÍMP2Ïà¹ØÄÜÉÏÒýÈëµÄÎó²î¾¡Á¿Ð¡¡£ÎÞÂÛÊÇÄÄÀàÃܶÈÄâºÏ»ù×飬Æä³ß´ç¶¼Òª±È±êÅäµÄ¹ìµÀ»ù×éÒª´ó£¬²¢ÇÒº¬Óиü¸ß½Ç¶¯Á¿º¯Êý¡£¶ÔÓÚͬһÖÖ¹ìµÀ»ù×飬ÕâÈýÀàÃܶÈÄâºÏ»ù×éµÄ´óС¹Øϵ×ÜÊÇ/C > /JK > /J£¬ÕâÊÇÒòΪ½ÏÉÙÊýÁ¿µÄ¸¨Öú»ùº¯Êý¾ÍÄÜÈÿâÂØÄܵļÆËãÎó²îÔڽϵÍˮƽ£¬¶øÈôÏë²»¸øÏà¹ØÄÜ´øÀ´Ã÷ÏÔÎó²î£¬¾Í²»µÃ²»Óøü¶àÒÔ¼°º¬Óиü¸ß½Ç¶¯Á¿µÄ¸¨Öú»ùº¯Êý¡£¸ß¼¶±ðÃܶÈÄâºÏ»ù×é¿ÉÒÔÓÃÔڵͼ¶±ð¹ìµÀ»ù×éÉÏ£¬¶øÇÒ±ÈʹÓñêÅäµÄÃܶÈÄâºÏ»ù×é½á¹û¸üºÃ£¬ÀýÈçcc-pVTZ/JKºÍcc-pVTZ/C¶¼¿ÉÒÔÓÃÓÚÔÚcc-pVTZÉÏ×öRI-J£¬¶øÇұȽáºÏcc-pVTZ/JʱÉÔ΢ºÃÒ»µãµã¡£ÔÙ±ÈÈç6-31G**¿ÉÒÔÓÃcc-pVTZ/J×öRI-J»òRIJCOSX£¬ËäÈ»cc-pVTZ/J²»ÊǸø6-31G**ÓÅ»¯µÄ£¬µ«ÒòΪcc-pVTZ±È6-31G**´óµÃ¶à£¬ËùÒÔËüµÄ±êÅäµÄ¸¨Öú»ù×éÓÃÓÚ6-31G**ºÁÎÞÎÊÌ⣬¾¡¹ÜÓÐЩÀË·Ñ¡£ÓеÄʱºòÎÒÃǻῳµôһЩÒâÒå²»´óµÄ¸ß½Ç¶¯Á¿º¯ÊýÀ´½µµÍ¼ÆËãʱ¼ä£¬±ÈÈçdef2-QZVPP(-g,-f)´ú±í´ÓÖØÔ­×ÓºÍÇâÔ­×ÓÉÏ·Ö±ðÈ¥µô×î¸ß½Ç¶¯Á¿º¯ÊýgºÍf£¬ÕâʱûÓбêÅäµÄ¸¨Öú»ù×飬µ«ÊÇÎÒÃÇ¿ÉÒÔ·ÅÐĵØÖ±½ÓÓÃdef2-QZVPP/J¡¢/JK¡¢/C¡£·´Ö®£¬¸ß¼¶±ð¹ìµÀ»ù×é½áºÏµÍ¼¶±ðµÄ¸¨Öú»ù×飬±ÈÈçdef2-TZVPÓÃdef2-SVP/J½á¹û¾Í»á½Ï²î£¬»òÕßÓÃ/J´úÌæ/JK¸¨Öú»ù×é×öRI-JK¼ÆËãÒ²»á´øÀ´½ÏÃ÷ÏÔµÄÎó²î¡£¶ÔÓÚûÓбêÅ丨Öú»ù×éµÄÇé¿ö£¨ÌرðÊǶÔÓÚPopleϵÁлù×飬Ŀǰ¶¼Ã»ÓбêÅäµÄ£©£¬Ò²ÓÐ×Ô¶¯µÄ·½·¨¿ÉÒÔÉú³É¸¨Öú»ù×飬¶øÇÒ½á¹ûÒ²¿É¿¿£¬µ«ÕâÑùµÃµ½µÄ¸¨Öú»ù×é±ÈÆðרÃÅÓÅ»¯µÄ¸¨Öú»ù×éÒª´ó²»ÉÙ£¬¹Ê¼ÆËãÃ÷ÏÔ¸üºÄʱ£¬ËùÒÔ½¨ÒéÔÚʹÓÃÃܶÈÄâºÏ¼ÆËãʱ»¹ÊÇÓ¦¾¡Á¿ÓÃÄÇЩÓкÏÊʵĸ¨Öú»ù×éµÄ¹ìµÀ»ù×飬def2-ϵÁлù×éµÄ¸¨Öú»ù×éºÜÈ«£¬Òò´ËÍùÍùÊÇÊ×Ñ¡¡£×îºóÖµµÃÒ»ÌáµÄÊÇÔÚORCA³ÌÐòÖУ¬¶ÔÓÚºóHFÈÎÎñ£¬¿ÉÒÔͬʱÓÃÁ½ÀàÃܶÈÄâºÏ»ù×飬¼´ÔÚSCF¹ý³ÌÖÐʹÓÃ/J»ò/JK£¬¶øÔÚÏà¹Ø²¿·Ö¼ÆËãÖÐʹÓøü´óµÄ/CÀ´±£Ö¤¾«¶ÈËðʧ½ÏС¡£


2.7 ORCA

ÒªÎÊÄĸöÁ¿×Ó»¯Ñ§³ÌÐò×îÓÐÇ°¾°£¬´ð°¸Ò»¶¨ÊÇORCA¡£ORCA´Ó90Äê´úºóÆÚ¿ªÊ¼¿ª·¢£¬ÊǸö·Ç³£ÓÅÐ㣬¹¦ÄÜÈ«Ã棬´úÂëЧÂʺܸߣ¬²»¶ÏѸËÙ·¢Õ¹ÍêÉƵijÌÐò£¬¶øÇÒ»¹ÊÇÃâ·ÑµÄ£¨²»¹ý²»¿ªÔ´£©¡£´Ó»ù̬µ½¼¤·¢Ì¬£¬´Óµ¥µãµ½¼¸ºÎÓÅ»¯£¬´Ó´óÌåϵDFT¼ÆË㵽СÌåϵ¸ß¾«¶ÈñîºÏ´Ø¼ÆË㣬ÒÔ¼°NEVPT2¡¢MRCIµÈ¶à²Î¿¼Ì¬¼ÆËãÈ«¶¼×öµÃͦºÃ¡£ÈܼÁÄ£ÐÍ¡¢±êÁ¿Ïà¶ÔÂÛ¼ÆËã¡¢NMR¼ÆËãµÈORCAÒ²¶¼Ö§³Ö¡£¶øÇÒORCA»¹Öð½¥Ö§³ÖÆäËü¶àÊý³ÌÐòûÓеŦÄÜ£¬ÖîÈçÃܶȾØÕóÖØÕû»¯Èº(DMRG)¡¢F12ÏÔʽÏà¹Ø¼ÆËã¡¢DLPNO-CCSD£¨ºÅ³ÆÄÜ°ÑñîºÏ´Ø¼ÆËãÓõ½ÉÏ°Ù¸öÔ­×Ó£©¡£ÌرðÊÇORCA³ä·ÖÀûÓÃÁËÃܶÈÄâºÏ¼¼Êõ£¬²¢ÇÒ½áºÏ¶À¼ÒµÄCOSX·½·¨£¬¼ÆËãËٶȱÈÆðÆÕͨÁ¿»¯³ÌÐòÓÐÁ˾ªÈ˵ķÉÔ¾£¬¸üÊǼ«´ó³Ì¶ÈµØµß¸²ÁË»ùÓÚDFTµÄÈõÏ໥×÷ÓüÆËãµÄ¸ñ¾Ö¡£ÔÚRI-JµÄʱºòORCA»¹Ä¬ÈÏʹÓÃRI-split-J·½·¨½øÒ»²½¼Ó¿ì¼ÆË㣬Õâ²»Ó°Ï쾫¶È£¬¶Ô´ó»ù×éʱ¼ÓËÙЧ¹ûÃ÷ÏÔ¡£ÁíÍâORCAÉÏÊÖÒ²±È½ÏÈÝÒ×£¬ÊÖ²áдµÃ±È½ÏÏêϸ£¬ÓÐÔ­Àí½éÉܺÍʵÀý¡£Ïë±ØÈç¹ûORCAÓиögviewË®×¼µÄͼÐνçÃ棬²¢ÇÒÔÙÓиöExploring Chemistry With Electronic StructureÕâÑùµÄ½Ì³ÌÊéµÄ»°£¬½«»á´ò¿ª³õѧÕßÊг¡£¬´Ó¶ø¶ÔÁ¿×Ó»¯Ñ§Ñо¿ÁìÓò²úÉú½«ÄÑÒÔ¹ÀÁ¿µÄÓ°Ï죬GaussianµÄÊг¡¿Ï¶¨»áÑÏÖØËõË®£¬¶øÏñQ-Chem¡¢TurbomoleµÈÆäËüÊշѳÌÐòµÄÉú´æ¿Õ¼ä¾Í¸üСÁË¡£ORCAÔâÈËÚ¸²¡µÄÊDz»Ö§³ÖÀûÓöԳÆÐÔÀ´¼ÓËÙ¼ÆË㣨µ«¿ÉÒÔÅж϶ԳÆÐÔ£©£¬Êµ¼ÊÉÏÕâ²»Ëã´óÎÊÌ⣬ÃܶÈÄâºÏ¼¼ÊõÒѾ­Ê¹ORCA±ÈºÜ¶à³ÌÐò¿ª¶Ô³ÆÐÔʱ¿ìµÃ¶àÁË¡£¶øÇÒÌرðÊǶÔÓÚ´óÌåϵÈõÏ໥×÷ÓõļÆË㣬Ìåϵһ°ã¶¼ÊÇûÓжԳÆÐԵģ¬ËùÒÔ²»Ö§³Ö¶Ô³ÆÐÔ²¢²»ÊÇÎÊÌâ¡£

ORCAºÍGrimmeÓнôÃܵĺÏ×÷£¬ËùÒÔDFT-D3¡¢gCPÒÔ¼°HF-3cµÈGrimmeÌá³öµÄÓëÈõÏ໥×÷ÓÃÏà¹ØµÄ·½·¨È«¶¼Ö±½Ó±»ORCAËùÖ§³Ö¡£

ORCAµÄÖ÷Ò³ÊÇhttp://www.cec.mpg.de/forum/£¬×«ÎÄʱ×îа汾ÊÇ3.0.1¡£ËäÈ»ÔÚÉÏÃæÃâ·Ñ×¢²áһϺó¾Í¿ÉÒÔÖ±½ÓÏÂÔØ£¬µ«ÏÂÔØËٶȿÉÄܱȽÏÂý£¬Òò´ËÎÒÒ²´«ÁËÒ»·Ýµ½°Ù¶ÈÍøÅÌÉÏ£ºhttp://pan.baidu.com/s/1eQkF2¡£

ORCA½âѹºó¾ÍÄÜÖ±½ÓÓã¬ÔËÐз½Ê½Èçorca Mio.txt > akiyama.out£¬ÕâÀï¼ÙÉèorcaĿ¼ÒѾ­ÄÉÈëÁËϵͳµÄPATH»·¾³±äÁ¿¡£ORCAÊÇÒÔMPI·½Ê½²¢Ðеģ¬Èç¹ûÒª²¢ÐÐÔËÐУ¬Ïȵ½http://www.open-mpi.orgÏÂÔØ×îеÄOpenMPI²¢°²×°µ½ÏµÍ³ÉÏ£¬È»ºóÊäÈëÎļþÀïдÉÏpal2¡¢pal3...pal8¾Í´ú±í2ºË¡¢3ºË...8ºË²¢ÐУ¬È»ºó´øמø¶Ô·¾¶È¥µ÷ÓÃorca¼´¿É£¨²»ÐèÒªmpirunÖ®Àࣩ£¬ÀýÈçd:\study\orca301\orca Love_live.txt > niconiconi.out¡£Èç¹ûÒªÓøü¶àºË¾ÍµÃÔÚÊäÈëÎļþÀïд³É%pal nprocs x end£¬xÊǵ÷ÓõĺËÊý£¬ÎÞÉÏÏÞ¡£

ORCAµÄÊäÈëÎļþдÆðÀ´±È½ÏÈÝÒ×£¬ÀýÈç
!BLYP D3 GCP(DFT/SVP) def2-SVP def2-SVP/J pal4 nopop
%maxcore 1000
* xyz 0 1
O                  4.71271832    9.26971072   -7.55066164
C                  4.13662696    9.48826045   -8.82968461
...ÂÔ
C                  2.75025629    8.84701245   -8.99475234
C                  7.66851269    1.97163671    1.43096468
*
!ºóÃæµÄ¹Ø¼ü´Ê´ú±íÓÃBLYP¼ÆË㣬´øDFT-D3УÕýÒÔ¼°DFT/SVP¼¶±ðµÄgCPУÕý£¬¹ìµÀ»ù×éΪdef2-SVP£¬def2-SVP/J×÷ΪRI-JµÄ¸¨Öú»ù×é¡£¶ÔÓÚ´¿·ºº¯ORCAÊÇĬÈÏ¿ªÆôRI-JµÄ£¬ËùÒÔ²»ÓÃд¶îÍâ¹Ø¼ü´Ê¿ªÆôËü¡£nopopÊDZÜÃâ¼ÆËãÍêÖ®ºóÊä³öÒ»´ó¶ÑÌرðÈß³¤Ò»°ãÓÖûÓõIJ¼¾Ó·ÖÎöÊý¾Ý£¨ÕâÀà·ÖÎö×îºÃÓÃMultiwfn³ÌÐò£©¡£maxcoreÉ趨ÿ¸öºË×î´óµÄÄÚ´æʹÓÃÁ¿(MB)£¬Êµ¼ÊʹÓÃÁ¿¿ÉÄܳ¬¹ýÖ®£¬±ÈÈç´ËÀý¿ÉÄÜ×ܹ²»áÓõ½5G¶à£¬ËùÒÔmaxcore±ðÉèµÃÌ«´ó£¬·ñÔòÒ»µ©ÎïÀíÄÚ´æ²»¹»¶øµ÷ÓÃÐéÄâÄÚ´æµÄ»°ËٶȽ«¾ÞÂý¡£xyz´ú±íÊäÈëµÄÊǵѿ¨¶û×ø±ê£¬Ä¬ÈÏΪ°£¡£0 1´ú±í¾»µçºÉΪ0£¬µ¥ÖØ̬¡£

ÔÙ±ÈÈç¹Ø¼ü´ÊÒ»ÐÐд³ÉÕâÑù
!RI-PWPB95 D3 def2-QZVPP def2-QZVPP/JK def2-QZVPP/C RIJK grid4 tightscf nopop pal4
¾Í´ú±íʹÓÃPWPB95Ë«ÔÓ»¯·ºº¯£¬SCF²½ÖèÓÃdef2-QZVPP/JK¸¨Öú»ù×é×öRI-JKÀ´¼ÓËÙ£¬·½·¨Ãû×ÖǰдÉÏÁËRI-˵Ã÷ÔÚ΢ÈżÆËãʱҲʹÓÃÃܶÈÄâºÏ½üËÆÀ´¼ÓËÙ£¬²¢ÇÒ»ùÓÚdef2-QZVPP/C¸¨Öú»ù×é¡£ÓÉÓÚ´ËÀý¼ÆË㾫¶ÈÒªÇó¸ß£¬ÓÃÁËtightscfÀ´ÉèÑÏÊÕÁ²ÏÞ¡£

ORCAÖн»»»Ïà¹Ø·ºº¯µÄ»ý·Ö¸ñµã°üº¬Æ߸ö¼¶±ð£¬´Ógrid1µ½grid7¾«¶ÈÒÀ´ÎÔö¼Ó¡£ORCAĬÈÏÔÚDFTµÄSCF¹ý³ÌÖÐʹÓýϵÍÖÊÁ¿µÄ»ý·Ö¸ñµã(Grid2)£¬¶øÔÚ×îºó¼ÆËã×îÖÕÄÜÁ¿Ê±²ÅÇл»µ½½Ï¸ßÖÊÁ¿¸ñµã(Grid4)¼ÆËã½»»»Ïà¹ØÄÜ£¬Õâ³ÆΪmultigrid·½·¨£¬Ê¹µÃ¼ÆËãÁ¿ÓÖС¾«¶ÈÓÖÓб£Ö¤¡£È»¶øSCF¹ý³ÌÖÐʹÓýϵ͸ñµãÖÊÁ¿»¹ÊÇ»á¸ø½á¹û´øÀ´¿ÉÒÔ²ì¾õµ½µÄÎó²îµÄ¡£ÎªÁ˲»ÒòΪ»ý·Ö¸ñµã´øÀ´µÄÎó²îÎÛȾÁËË«ÔÓ»¯·ºº¯+def2-QZVPPÕâô¸ß¼¶±ðµÄ¼ÆËãË®×¼´øÀ´µÄÓÅÐã½á¹û£¬ÉÏÀýʹÓÃÁËgrid4¹Ø¼ü´ÊÀ´ÈÃSCF¹ý³ÌҲʹÓýϺÃÖÊÁ¿µÄ»ý·Ö¸ñµã£¬´ËʱÔÚ¼ÆËã×îÖÕÄÜÁ¿Ê±³ÌÐò»á×Ô¶¯Çл»µ½¸ü¸ßµÄGrid5¡£¶ÔÓÚÒ»°ã·ºº¯ÏµļÆËãÒ»°ã²»ÐèÒªÌØÒâÖ¸¶¨¸ñµã¾«¶ÈºÍµ÷ÕûÊÕÁ²ÏÞ£¬ÓÃĬÈϵľÍÐÐÁË£¬µ«¶ÔÓÚƵÂʼÆËãµÈÈÎÎñ½¨Òé×ÜÊǼÓÉÏgrid4 tightscf¡£Èç¹û»¹Ïë½øÒ»²½Ìá¸ß¸ñµã¾«¶È£¨ËäÈ»ÒÑÍêȫûÓбØÒª£©£¬¿ÉÒÔдgrid5£¬²¢ÇÒ»¹¿ÉÒÔͨ¹ýдÉÏ%method FinalGrid 6 endÀ´ÈÃ×îºóÒ»´Î¼ÆËãʱ×Ô¶¯Çл»µ½¸ü¸ßgrid6¡£

˳´øÒ»ÌᣬORCA¼ÆËãºó»á²úÉú.gbwÎļþ£¬°üº¬»ù×é¡¢¹ìµÀµÈÐÅÏ¢¡£Í¨¹ý×Ô´øµÄ¹¤¾ßorca_2mkl xx -molden¿ÉÒÔ½«xx.gbwת»»ÎªmoldenÊäÈëÎļþ£¨ºó׺ÃûÓ¦¸ÄΪ.molden£©£¬Í¨¹ýorca_2aim xx¿ÉÒÔ½«xx.gbwת»»Îª.wfn²¨º¯ÊýÎļþ¡£²¨º¯Êý·ÖÎö¹¤¾ßMultiwfn£¨http://multiwfn.codeplex.com£©Ä¿Ç°×îа汾3.2.1(dev)¿ÉÒÔÔØÈëÕâÁ½¸öÎļþÀ´×ö²¨º¯Êý·ÖÎö£¬°üÀ¨¸÷ÖÖ¸÷ÑùµÄÈõÏ໥×÷ÓõķÖÎö£¬ÈçRDG·ÖÎö¡¢AIM·ÖÎö¡¢¶¨Á¿·Ö×Ó±íÃæ¾²µçÊÆ·ÖÎöµÈµÈ£¬²Î¼û¡¶Multiwfn 3.0²¨º¯Êý·ÖÎö³ÌÐòµÄÒâÒå¡¢¹¦ÄÜÓëÓÃ;¡·£¨http://hi.baidu.com/sobereva/item/b35581418db96cea1f19bc09£©¡£Ê¹ÓÃMultiwfnµÄÖ÷¹¦ÄÜ0»¹¿ÉÒԺܷ½±ãµØ²é¿´·Ö×Ó¹ìµÀͼÐΣ¬µ«Ö»ÓÐ.molden²Åº¬ÓпչìµÀÐÅÏ¢£¬.wfnÖ»¼Ç¼ռ¾Ý¹ìµÀÐÅÏ¢¡£ÁíÍâ×¢Ò⣬MultiwfnµÄһЩ¹¦ÄÜÒÀÀµÓÚ»ùº¯ÊýÐÅÏ¢£¬ÖîÈçMayer¼ü¼¶·ÖÎö¡¢SCPA¹ìµÀ³É·Ö·ÖÎöµÈ£¬Ö»ÓÐ.molden²ÅÌṩÕâÑùµÄÐÅÏ¢£¬¶ø.wfn²»Ìṩ¡£Ò²¼´ÊÇ˵£¬ÓÃ.molden×÷ΪMultiwfnÊäÈëÎļþÄÜ×öµÄʱÈÓÃ.wfn×÷ΪÊäÈëʱ¸ü¶à¡£
×¢£ºÓÉÓÚORCA²úÉúµÄ.molden¸ñʽºÜ²»±ê×¼£¬MultiwfnÄ¿Ç°²»Ö§³Öº¬ÓÐg½Ç¶¯Á¿µÄORCAµÄ.moldenÎļþ£¬Êµ¼ÊÉÏÓú¬Óиߴïg½Ç¶¯Á¿µÄÄÇô´ó»ù×é×ö²¨º¯Êý·ÖÎö¸ù±¾Ã»±ØÒª£¬²¨º¯Êý·ÖÎö½á¹ûԶûÓÐÄÜÁ¿¶Ô»ù×éÃô¸Ð¡£µ«Èç¹û±ØÐëÒªÉæ¼°g£¬¿ÉÒÔÓÃmolden2aim³ÌÐò£¨http://people.smu.edu/wzou/program/index.html£©ÔØÈ뺬gµÄ.moldenÎļþ£¬²»½ö»áת»»³ö.wfnÎļþ£¬»¹»á²úÉú±ê×¼»¯ºóµÄ.moldenÎļþ£¬Õâʱ¾Í¿ÉÒÔÕý³£ÔØÈëMultiwfn½øÐзÖÎöÁË¡£


2.8 ÔÓ̸ÆäËüµÄ¼ÓËÙ·½·¨

Ö±½ÓÓôÓÍ·ËãËã²»¶¯µÄÌåϵÈËÃÇÓÐʱ¿¼ÂÇÓÃÔÓ»¯·½·¨£¬±ÈÈçͨ¹ýONIOM¼¼Êõ½«°ë¾­Ñé¡¢·Ö×ÓÁ¦Ñ§ºÍ´ÓÍ·Ëã½áºÏÈýÕßÏ໥×éºÏ¡£ONIOM£¬»òÕ߸üÒ»°ãµÄQM/MM¼¼Êõ¶ÔÓÚÑо¿´óÌåϵȷʵºÜÓÐÓ㬵«ÊǶÔÓÚ´óÌåϵµÄÈõÏ໥×÷ÓüÆË㣬±ÊÕßÈÏΪÕâЩ·½·¨Ã»ÓÐÌ«´óÓÃÎäÖ®µØ£¬Ô­ÒòºÜ¶à¡£Ô­ÒòÖ®Ò»ÊÇPM7¡¢PM6-DH+µÈ·½·¨½áºÏMOZYME¼¼ÊõÒѾ­ÄÜ´¦ÀíºÜ´óÌåϵÁË£¬ÓÃÓÚÉúÎï´ó·Ö×Ó¶¼Ã»ÎÊÌ⣬½á¹ûÒ²²»´í£¬ËùÒÔÓöÉÏ´ÓÍ·ËãËã²»¶¯µÄÇé¿öÖ±½ÓÓð뾭Ñé¾ÍÐÐÁË£¬Óò»×Å°Ñ·½·¨½øÐÐÔÓ»¯¡£ÓÃÁËÔÓ»¯·½·¨£¬¿ÉÄÜÓÉ´Ë´øÀ´µÄÎó²î¶¼±È°ë¾­Ñé·½·¨±¾ÉíµÄÎó²îÒª´óºÜ¶à¡£Ôõôȥ»®·Ö²»Í¬ÀíÂÛ¼¶±ð´¦ÀíµÄÇøÓòÒ²ÊǼþÂ鷳ʡ£¶øÇÒÈôÊÇ¿¿±ÈÈçONIOM½áºÏ´ÓÍ·ËãÓë°ë¾­Ñ飬GaussianÏÖÔÚÖ§³ÖµÄ°ë¾­ÑéÓÖûÓжÔÈõÏ໥×÷ÓÃËãµÃºÃµÄ£¬µÍ¼¶±ðÇøÓòÈõÏ໥×÷Óÿ¼ÂǵÃÌ«²îÒ²²»ÐС£Èç¹ûÊÇ´ÓÍ·Ëã»ò°ë¾­ÑéÓëMM»ìºÏµ½Ò»ÆðÓã¬ÓÖÒªÉæ¼°Ö¸¶¨Á¦³¡²ÎÊýµÄÎÊÌ⣬ÍùÍùÔÚ²Ù×÷ÉϸãµÃºÜÂé·³¡£×ܶøÑÔÖ®£¬ONIOM¡¢QM/MM¶ÔÓÚ´óÌåϵÈõÏ໥×÷ÓõÄÑо¿ºÜÄÑÅÉÉÏÓó¡¡£

»¹ÓÐÒ»ÖÖ¼¼ÇÉÊÇʹÓûìºÏ»ù×é¼õС¼ÆËãÁ¿£¬±ÈÈçÒ»¸ö´ó·Ö×ÓºÍÒ»¸öС·Ö×Ó½áºÏÒªËã½áºÏÄÜ£¬Ö»°Ñ½áºÏµÄÇøÓò¸½½üµÄÔ­×ÓÓôøÃÖÉ¢µÄ´ó»ù×飬ÕâÑù×öûÓÐÎÊÌ⣬ʵÏÖÒ²ÈÝÒס£²»¹ýÒ»°ãµ¹Ò²Ã»±ØÒªÓã¬ÒòΪÈç¹ûʹÓÃÏÂÒ»½ÚÍƼöµÄ·½·¨£¬ÄÄÅÂÊÇ´óÌåϵ£¬Ò²ÄÜËãµÃͦ¿ìͦ׼¡£

ÁíÒ»ÖÖ´¦Àí´óÌåϵ³£ÓõÄ×ö·¨¾ÍÊÇ°ÑÌåϵ¼ò»¯³ÉÄ£ÐÍ£¬Ö»°Ñ¹Ø¼üÇøÓò¿Ù³öÀ´£¨±ÈÈçµ°°×ÖÊ»îÐÔλµã£©£¬°Ñ±ßÔµÔ­×ÓÓñÈÈç¼×»ù¸ø±¥ºÍÉÏ£¬ÆäÓàÔ­×ÓҪôºöÂÔҪôÓñ³¾°µçºÉ·½Ê½À´±íÏÖ³ö¾²µç×÷Óá£Õâô×öÍêÈ«¿ÉÒÔ£¬ËäÈ»±íÃæÉÏûQM/MMÑϸñ£¬µ«ÊÇʡʲ»ÉÙ¡£²»¹ýÈçºÎ¼ò»¯»¹ÊÇÐèÒªÒ»¶¨¼ÆË㻯ѧ¾­ÑéºÍ»¯Ñ§Ö±¾õµÄ£¬¼ò»¯µÃ¹ý¶È½á¹û²»×¼£¬¼ò»¯²»¹»Ôò¼ÆËãÁ¿»¹ÊÇÌ«´ó¡£

»¹ÓкܶàÏßÐÔ±ê¶ÈµÄ·½·¨¶¼ÄÜ°ïÖú´¦Àí´óÌåϵ£¬MOZYMEÇ°ÃæÒѾ­Ìá¹ýÁË£¬ÕâÀïÖ÷Ҫ˵µÄÊÇ´ÓÍ·Ëã/DFT²ãÃæÉϵÄÏßÐÔ±ê¶È·½·¨¡£ÕâÖÖ·½·¨·Ç³£¶à£¬Ô­ÀíÒ²²»Í¬£¬´óÖÂ˼ÏëҪôÊÇ»ùÓÚ¶¨Óò»¯¹ìµÀ£¬ÒªÃ´ÊÇ°ÑÌåϵ·Ö³ÉƬ¶ÎÒÀ´Î¼ÆËãÔÙ×éºÏÉÏ¡£ÎÞÊýÈËÔÚÕâ·½Ãæ×öÁ˼¸Ê®ÄêµÄÑо¿£¬Ïà¹ØpaperÎÞÊý£¬Ã»ÍêûÁ˵Øɹ¼ÓËٱȣ¬µ«Ê¼ÖÕû²úÉú¶àÉÙʵ¼ÊÓ¦ÓüÛÖµ¡£¹Ø¼üµÄÔ­Òò¾ÍÊÇ·½·¨²»¹»ºÚÏ仯£¬ÄÜÓ¦ÓûòÕßÊÊÓõÄÌåϵ·Ç³£ÏÁÕ­£¬±ÈÈçÖ»¶ÔÖ±Á´Ìåϵ¼ÓËÙЧ¹ûºÃ¶øÅöÉϸ´ÔÓµãµÄ½á¹¹¾ÍÑÆ»ð£¬Óеķ½·¨´¦Àí²»Á˺¬Óдó¹²éîµÄÌåϵ£¬Óеķ½·¨ÔòÖ»ÏÞÓÚ±ÈÈçºÜÀíÏ뻯µÄµ°°×Ìåϵ£¬µÈµÈ¡£×ÜÖ®´ÓÔ­ÎÄÊý¾ÝÉÏÀ´¿´½á¹ûÉõºÃ£¬µ«ÊÇʵ¼ÊÓ¦ÓÃÖе½´¦Åö±Ú£¬²Ù×÷ÆðÀ´ÉõÊÇÂé·³¡£Ã»ÓÐÁ÷ÐÐÆðÀ´µÄÁíÒ»¸öÖØÒªÔ­ÒòÊDZ»Ö÷Á÷³ÌÐòÖ§³Ö²»¹»£¬Ö÷Á÷µÄÁ¿»¯³ÌÐò´ó¶à¶¼Ã»Ö§³ÖÄÇЩ×÷Õß×Ô´µÊÇÏßÐÔ±ê¶ÈµÄ·½·¨£¬ºÜ´ó³Ì¶ÈÉÏÒ²ÊÇÒòΪÕâЩ·½·¨ÊµÓÃÐÔ¡¢Êµ¼ÊЧ¹û²î£¬Ïà¶Ô¶øÑԽϺõÄÒ²¾ÍÊÇFMO (Fragment MO)·½·¨£¬ÒÑÔÚGAMESS-USÀï±»³ä·ÖÖ§³Ö¡£×ÜÖ®Ïֽ׶ÎÀ´¿´ÄÇЩÏßÐÔ±ê¶È·½·¨±íÃæÉÏÏԵöÔÓÚÑо¿´óÌåϵÈõÏ໥×÷ÓúÜÓÐÓô¦£¬µ«¾ø´ó¶àÊý¶¼Í£ÁôÔÚÀíÂÛ²ãÃæÉÏ£¬Êµ¼ÊÉÏÓ¦ÓüÛÖµ²»´ó£¬³ý·Ç×Ô¼ºÒѾ­°ÑÕâ·½ÃæÀíÂۺͲÙ×÷Ñо¿µÃ¯»ð´¿ÇàÁË£¬ÄÜÍæµÃµÃÐÄÓ¦ÊÖ¡£

ORCA´Ó3.0¼ÓÈëµÄDLPNO-CCSD(T)¼¼ÊõÖµµÃÒ»Ìá¡£CCSD(T)ÖÚËùÖÜÖªºÄʱ¼«¸ß£¬ÇÒËæ×ÅÌåϵ³ß¶ÈµÄÔö¼Ó¼ÆËãÁ¿ì­Éý£¬¶þÊ®¶à¸öÒÔÉÏÔ­×ÓµÄÌåϵ¶¼¼¸ºõËã²»¶¯¡£DLPNO-CCSD(T)»ùÓÚ¾ÖÓòÏà¹Ø˼Ï룬ÊÇÇ°ÉíLPNO-CCSD(T)µÄ¸Ä½ø°æ£¬ÁîÎÞÂÛ¶à´óÌåϵµÄCCSD(T)¼ÆË㶼ÓÐÏßÐÔ±ê¶ÈµÄЧ¹û£¬ºÅ³ÆÄ¿Ç°ÒѾ­ÄÜ°ÑCCSD(T)×öµ½650¸öÔ­×ӵĵ°°×£¬¶øÇÒ´Ë·½·¨»ù±¾ÊǺÚÏ䣬²»ÐèÒªÓû§µ÷ʲô²ÎÊý£¬ºÜÓÐÇ°¾°¡£ÀíÂÛÉÏ˵DLPNO-CCSD(T)¶ÔÓÚ¾«È·Ñо¿´óÌåϵÈõÏ໥×÷ÓúÜÓÐÓᣲ»¹ý±¾ÎIJ¢²»Éæ¼°DLPNO-CCSD(T)£¬ÒòΪËü±ÈÆðORCA×öË«ÔÓ»¯·ºº¯ºÄʱԶԶҪ¸ßµÃ¶à£¬¶øÇÒΪÁËÈÃCCSD(T)³ä·Ö·¢»ÓЧ¹û±ØÐëµÃÓÃÆÄ´óµÄ»ù×飬¸üÊǸø¼ÆËãÑ©ÉϼÓ˪£¬¶øÇÒDLPNO¼¼Êõ±¾À´Ò²»á´øÀ´Ò»Ð©Îó²î£¬ËùÒÔ´ÓʵЧÉϽ²Ô¶²»ÈçË«ÔÓ»¯·ºº¯ºÏËã¡£



3 ²»Í¬³ß¶ÈÈõÏ໥×÷ÓÃÌåϵµÄ×î¼Ñ¼ÆËã·½·¨

»ùÓÚÉÏÃæ½éÉܵÄ7¸ö·¨±¦£¬¶ÔÓÚ²»Í¬³ß¶ÈµÄÌåϵµÄÈõÏ໥×÷ÓÃÑо¿£¨ÄÜÁ¿¼ÆËã¶ø·Ç¼¸ºÎÓÅ»¯£©±ÊÕßÍƼöµÄ¼ÆËã·½·¨ÈçÏ£¬´Ó´óÌåϵµ½Ð¡ÌåϵÅÅÐò£¬¹²·ÖΪ8¸öµµ¡£ÕâЩÌåϵ¼Ù¶¨ÊÇÒÔC¡¢H¡¢N¡¢OΪÖ÷µÄÓлú·Ö×Ó¡£µÚ1ÏîÊÇÓÃMOPAC2012À´×ö¼ÆË㣬2~6ÏîÓÃORCAÀ´×ö¼ÆË㣬ÁгöµÄ¾ÍÊÇÔÚORCAÖÐʵ¼ÊҪдµÄ¹Ø¼ü´Ê¡£¶ÔÓÚ7¡¢8£¬ÓÃʲô³ÌÐòÀ´×ö¶¼¿ÉÒÔ£¬ÓÉÓÚÕâÊôÓÚСÌåϵ·¶³ëËùÒÔÔÚ±¾ÎÄÖв»×öÌÖÂÛ£¬ÖµµÃÒ»ÌáµÄÊÇORCAµÄCCSD(T)¼ÆËãЧÂÊÊǺܸߵġ£

1 >500¸öÔ­×Ó£ºPM7»òPM6-DH+½áºÏMOZYME
2 >300¸öÔ­×Ó£ºBLYP D3 GCP(DFT/SVP) def2-SVP def2-SVP/J
3 >150¸öÔ­×Ó£ºBLYP D3 GCP(DFT/TZ) def2-TZVP def2-TZVP/J
4 >60¸öÔ­×Ó£ºBLYP D3 def2-QZVPP(-g,-f) def2-QZVPP/J
5 >40¸öÔ­×Ó£ºRI-PWPB95 D3 RIJCOSX def2-QZVPP(-g,-f) def2-QZVPP/J def2-QZVPP/C grid4 tightSCF£¨/J¿É¸ÄÓÃ/JKÒÔ¼õСÃܶÈÄâºÏµÄÎó²î£©
6 ¶þ¡¢ÈýÊ®¸öÔ­×Ó£ºRI-PWPB95 D3 RIJK def2-QZVPP def2-QZVPP/JK def2-QZVPP/C grid4 tightSCF
7 Ê®¼¸¸öÔ­×Ó£ºCCSD(T)/jul-cc-pVTZ»òmay-cc-pVQZ
8 <10¸öÔ­×Ó£ºCCSD(T)/CBS £¨aug-cc-pVTZ->QZÍâÍÆ¡£Ïë¸ü¾«È·ÔÙ¼ÓÉÏcounterpoiseУÕý£¬ºÄʱ½«Ôö¼Ó½ü2±¶£©

ÓÃÉÏÃæÕâЩ¼ÆËã·½·¨¼ÆËãÏàÓ¦³ß¶ÈµÄÌåϵ£¬ÔÚÄ¿Ç°Ö÷Á÷µÄIntel 4ºËCPUµÄ»ú×ÓÉÏ£¬Ò»°ã¶¼ÊÇ¿ÉÒÔÔÚ¼¸¸öСʱÄÚ¼ÆËãÍê±ÏµÄ¡£ÕâЩ¼ÆËã·½·¨µÄÈõÏ໥×÷ÓüÆËãÎó²îºÍ¶ÔÓ¦µÄÌåϵµÄ³ß´çÊÇÕýÏà¹ØµÄ¡£ÌåϵԽ´ó£¬¿ÉÒÔÈÌÊܵľø¶ÔÎó²îµ±È»Ò²¾ÍÔ½´ó¡£ÎÞÂÛÊÇÉÏÃæËùÁеÄÄĸö¼ÆË㼶±ð£¬¶ÔÏàÓ¦³ß¶ÈµÄÌåϵÀ´ËµÈõÏ໥×÷Óþ«¶È¶¼ÊÇ×ãÒÔÁîÈËÂúÒâµÄ¡£

Èç¹û¼ÆËãÄÜÁ¦±È½ÏºÃ£¬±ÈÈçÓÐ16ºË»ú×Ó£¬»òÕß²»Å»¨Ê±¼ä£¬±ÈÈçÒ»¸öÌåϵË㼸ÌìÒ²ÎÞ·Á£¬ÄÇôËùÓõļÆË㼶±ð¾Í¿ÉÒÔÌá¸ßÒ»Á½¸öµµ¡£·´Ö®£¬Èç¹û´òËã¼ÆËã´óÅúÁ¿Ìåϵ£¨±ÈÈçÓÃÓÚ·Ö×ÓɸѡµÄÄ¿µÄ£©£¬»òÕßÓÃÓÚ³õ²½¹ÀËã½á¹û£¬»òÕßÓÃÓÚ¼¸ºÎÓÅ»¯µÈ£¬ÄÇô¾Í¿ÉÒ԰ѼÆËãË®×¼½µµÍÁ½¡¢Èý¸öµµ¡£ÁíÍâ±ÊÕßÌرðÍƼöʹÓÃPM7»òPM6-DH+¶Ô¼¸Ê®¸ö»ò¸ü¶àÔ­×ÓÊýÄ¿µÄÌåϵ½øÐÐÓÅ»¯£¬ËäÈ»¼¸Ê®¸öÔ­×ÓµÄÌåϵÓÃDFTÒ²ÕÕÑùÄÜÓÅ»¯µÃ¶¯£¬µ«ÓÉÓÚÔÚ°ë¾­ÑéÏÂÓÅ»¯Ëٶȼ«¿ì£¬¼¸ºõ²»·Ñʱ¼ä¾ÍÄܵõ½½ÏºÃµÄ³õʼ¹¹ÐÍ£¬ÔÙ½Ó×ÅÓÃDFTÀ´ÓÅ»¯¿ÉÒÔÈÃÊÕÁ²¿ìµÃ¶à¡£

ÕâЩÍƼöµÄ¼ÆËã·½·¨ÀïÓÐһЩֵµÃ˵Ã÷ºÍÌÖÂ۵ĵط½£¬Ç°ÎÄÖÐÉÐδÌá¼°£¬½«ÔÚÏÂÃæ˵˵¡£

µÚ8¸öµµµ±ÖÐÓÃÁË»ù×éÍâÍƼ¼ÊõÒÔ´ïµ½Í걸¼¯£¬ÍêÈ«Ïû³ýÒòÓÐÏÞ»ù×é¶øµ¼ÖµļÆËãÎó²î£¬ÕâÔÚ¡¶Ì¸Ì¸ÄÜÁ¿µÄ»ù×éÍâÍÆ¡·£¨http://hi.baidu.com/sobereva/item/66cff4055a99e718ebfe38df£©µ±ÖÐÓÐÏêϸÌÖÂÛ¡£ÔÚµÚ7¸öµµµ±ÖÐÓÃÁËÔ·ݻù×飬ÕâÖÖ»ù×é°Ñaug-cc-pVnZϵÁлù×éµÄ¶Ô½á¹ûÓ°Ïì²»´óµÄ¸ß½Ç¶¯Á¿ÃÖÉ¢º¯ÊýÊʵ±µØ³ýÈ¥£¬´ó´ó½ÚÔ¼Á˼ÆËãÁ¿£¬ÕâÔÚ¡¶Ì¸Ì¸ÃÖÉ¢º¯ÊýºÍ¡°Ô·ݡ±»ù×é¡·£¨http://hi.baidu.com/sobereva/item/8201de10d92ddd5af0090ec5£©µ±ÖÐÓÐÏêϸÌÖÂÛ¡£ÔÚµÚ5¡¢6¸öµµµ±ÖÐÓÃÁËPWPB95-D3£¬Õ⼸ºõÊǶÔÓÚÈõÏ໥×÷ÓüÆËã×îºÃµÄË«ÔÓ»¯·ºº¯£¬ºÃÓÚË«ÔÓ»¯·ºº¯ÖÐ×îÁ÷ÐеÄB2PLYP-D3¡£

½áºÏDFT-D3УÕýºóÃܶȷºº¯·½·¨¶ÔÈõÏ໥×÷ÓüÆË㾫¶ÈÕûÌåÀ´ËµÊÇ£ºË«ÔÓ»¯·ºº¯>ÔÓ»¯·ºº¯>GGA·ºº¯£¬¼ÆËãºÄʱҲÊÇÕâ¸ö˳Ðò¡£Ö®ËùÒÔÉÏÃæÍƼöµÄ·½·¨ÖÐÖ±½Ó´ÓË«ÔÓ»¯·ºº¯(PWPB95)Ö±½ÓÇл»µ½GGA·ºº¯(BLYP)£¬¶øûÓо­¹ýÆÕͨÔÓ»¯·ºº¯µÄ¹ý¶É£¬ÊÇÒòΪ˫ÔÓ»¯·ºº¯µÄ¼ÆËãÁ¿²»±ÈÔÓ»¯·ºº¯¶à¶àÉÙ£¬¶îÍâ¶à³öµÄ΢ÈŲ¿·ÖµÄ¼ÆËãºÄʱÍùÍù»¹ÉÙÓÚÔÓ»¯·ºº¯ÔÚSCFµü´ú¹ý³ÌÉϵĺÄʱ¡£ËùÒÔÈç¹ûÓõÃÆðÔÓ»¯·ºº¯µÄ»°»¹²»Èç¸É´à¾ÍÓÃË«ÔÓ»¯·ºº¯À´µÃµ½¸üºÃµÄ¾«¶È¡£¶øGGA·ºº¯²»ÐèÒª¼ÆËãHF½»»»²¿·Ö£¬¿¿RI-J¼¼Êõ¿ÉÒÔÈÃGGA·ºº¯µÄ¼ÆËãËٶȱÈÔÓ»¯·ºº¯£¨ÄÄÅÂÒѾ­ÓÃÁËRI-JK»òRIJCOSXÀ´¼ÓËÙ£©¸ßºÃ¼¸±¶ÉõÖÁÉõÖÁÒ»¸öÊýÁ¿¼¶£¬ËùÒÔÔÓ»¯/Ë«ÔÓ»¯·ºº¯ÄÑÒÔ¶Ô¸¶µÄ´óÌåϵ¿ÉÒÔ¿¿GGA·ºº¯ÇáËɵضԸ¶¡£ËäÈ»ÔÓ»¯·ºº¯ÖÚËùÖÜÖª¶ÔÓÚÒ»°ãÎÊÌâÖîÈçÈÈÁ¦Ñ§Á¿µÄ¼ÆËã±ÈÆðGGA·ºº¯Ç¿µÃ¶à£¬µ«ÊǶÔÓÚÈõÏ໥×÷ÓõļÆË㣬½áºÏDFT-D3ºóËüÃǵÄÐÔÄܲîÒìÆäʵºÜС£¬ÀýÈçB3LYP-D3Ö»ÊDZÈBLYP-D3ÔÚ¶¨Á¿ÉÏÂÔÓеã¶ùÓÅÊƶøÒÑ£¬¸ü²»»áÓÐʲô¶¨ÐÔ²îÒì¡£ËùÒÔ´Óʵ¼Ê½Ç¶ÈÀ´¿´ÔÓ»¯·ºº¯+D3µÄÐԼ۱Ȳ»¸ß£¬¹ÊûÓÐÄÉÈëÍƼö¡£

def2ϵÁлù×éÊÇAhlrichsµÈÈËÔÚ2005ÄêÌá³öµÄ£¬ÊÇTurbomole³ÌÐòµÄÓùÓûù×飬def2-QZVPPÊÇÆäÖÐ×î´óµÄ£¬ºÍcc-pVQZÔÚÒ»¸ö¼¶±ð¡£def2-QZVPPÏà¶ÔÓÚdef2-QZVPÖ»ÊǶÔÓÚijЩԪËذѸ߽Ƕ¯Á¿»ùº¯ÊýÊýÄ¿ÉÔ΢Ôö¼ÓÁËÒ»µã¶øÒÑ£¬¶ÔÓÚC¡¢H¡¢N¡¢OûÓÐÇø±ð¡£def2-QZVPP(-g,-f)ÊÇdef2-QZVPPµÄ¸ß½Ç¶¯Á¿Ñ˸î°æ£¬ÕâÔÚÇ°ÃæÒѾ­Ìáµ½ÁË£¬Ñ˸Öµľ«¶ÈËðʧ¼¸ºõ¿ÉÒÔºöÂÔ¡£Ö®ËùÒÔÍƼöµÄ¼ÆË㼶±ðÖÐʹÓõÄÊÇdef2-QZVPP(-g,-f)»òdef2-QZVPPȴûÓÐʹÓöÔÓÚÈõÏ໥×÷ÓüÆËãºÜ³£ÓõÄaug-cc-pVTZ£¬ÕâÊÇÓм¸¸öÖØÒªÔ­ÒòµÄ£º
(1)def2-QZVPP(-g,-f)±¾ÉíÖ»±Èaug-cc-pVTZ´óÒ»µãµã¶øÒÑ£¬µ«aug-cc-pVTZĿǰûÓÐ/J¶øÖ»ÓÐ/JKºÍ/C¸¨Öú»ù×飬aug-cc-pVTZ/JK±ÈÆðdef2-QZVPP/J´óºÜ¶à¡£ÓÃaug-cc-pVTZ/JKµ±¸¨Öú»ù×é×öRI-JºÍRIJCOSX¼ÆËã±ÈÆðÓÃdef2-QZVPP/J¸üºÄʱ£¬ËùÒÔ´Óʵ¼ÊºÄʱÉÏaug-cc-pVTZ³Ô¿÷¡£
(2)aug-cc-pVTZµÄÃÖÉ¢º¯ÊýÈÝÒ×Òý»ùº¯ÊýÏßÐÔÒÀÀµÎÊÌâµ¼ÖÂÊÕÁ²À§ÄÑ¡£¶ødef2-QZVPP/(-g,-f)µÄ×îÍâ²ã»ùº¯ÊýÊÇ°ëÃÖÉ¢µÄ£¬¼ÈÊʺÏÈõÏ໥×÷ÓüÆË㣬Ҳ²»ÖÁÓÚ´øÀ´ÊÕÁ²ÉϵÄÎÊÌâ¡£
(3)DFT-D3²ÎÊý±¾ÉíÊÇÔÚdef2-QZVPÏÂÄâºÏµÄ£¬ÓÃÔÚdef2-QZVPP/(-g,-f)ÉÏÒ²¿ÉÒÔ˵ÊÇÍêÃÀµÄ£¬¶øaug-cc-pVTZËäÈ»±¾ÉíºÜÊʺÏÈõÏ໥×÷ÓõļÆË㣬µ«¶ÔÓÚʹÓÃDFT-D3УÕýµÄÇé¿ö±ÈÆðdef2-QZVPP/(-g,-f)ÊÇûÓÐÓÅÊƵġ£¾¡¹Üaug-cc-pVTZµÄÃÖÉ¢º¯ÊýʹµÃBSSEÎÊÌâ±Èdef2-QZVPP¸üС£¬µ«ºóÕßµÄBSSEÎÊÌâÆäʵÒѾ­ºÜСÁË£¬¶øÇÒÕâÒ»¶¡µãÎÊÌâÒѾ­±»µÈЧµØÎü½øDFT-D3µÄ²ÎÊýÀïÁË¡£
µÚ(3)µãÌÖÂÛÒ²±íÃ÷û±ØÒª£¬Ò²²»Ó¦¸ÃÔÚDFT-D3УÕý½áºÏdef2-QZVPP/(-g,-f)ʹÓÃʱʹÓÃCPУÕý£¬²»¹âÊÇÒòΪУÕýÓàµØÒѾ­ºÜС£¬¸üÒòΪDFT-D3²ÎÊýÒѾ­µÈЧµØ¿¼ÂÇÁË´ËÎÊÌ⣬ÔÙÈ¥×öCPУÕý·´µ¹»áÒýÈëÎó²î¡£

ËäÈ»±¾ÎÄÌÖÂÛµÄÊÇÈõÏ໥×÷ÓüÆË㣬µ«ÊÇÉÏÃæÍƼöµÄ8ÖÖ¼ÆË㼶±ð¾ÍËãÊǶÔÓÚÒ»°ãÀàÐÍÎÊÌâµÄÑо¿Ò²ÊÇÊ®·ÖÀíÏëµÄ¡£µ«ÊÇ´Ëʱdef2-QZVPP¾Íû±ØÒªÓÃÁË£¬½¨Òé½µµ½def2-TZVPPÀ´½Úʡʱ¼ä¡£

MP2ºÍË«ÔÓ»¯·ºº¯µÄ»ý·Ö±ä»»¹ý³ÌÊǺܺĴ¢´æ×ÊÔ´µÄ£¬¶ÔÓÚ´óÌåϵ£¬¹âÊÇ´¢´æÔÚÄÚ´æÍùÍù²»¹»£¬»¹ÐèÒªÔÚÓ²ÅÌÉÏ´óÁ¿¶ÁдÊý¾Ý¡£ÔÚORCAÀï¿ÉÒÔÖ¸¶¨»ý·Ö´¢´æ·½Ê½À´½µµÍ´¢´æÏûºÄ£¬ÒÔ¾¡Á¿¼õÉÙÉõÖÁ±ÜÃâÒòÓ²Å̵ĶÁд¶øÍÏÂýËٶȡ£´¢´æ·½Ê½°üÀ¨UCDOUBLE¡¢CDOUBLE¡¢UCFLOATºÍCFLOATËÄÖÖ¡£UC´ú±í²»Ñ¹ËõÊý¾Ý£¬C´ú±íѹËõ¡£Ñ¹²»Ñ¹Ëõ²»Ó°Ïì½á¹û£¬Ñ¹Ëõºó¿ÉÒ»¶¨³Ì¶È¼õÉÙ´¢´æ×ÊÔ´µÄʹÓÃÁ¿£¨¼õÉٵñȽÏÓÐÏÞ£©¡£DOUBLE´ú±íË«¾«¶È·½Ê½´¢´æ»ý·Ö£¬FLOAT´ú±íµ¥¾«¶È·½Ê½´¢´æ»ý·Ö£¬±ÈË«¾«¶ÈËùÐè´¢´æ×ÊÔ´¼õСÁËÒ»±¶¡£¶ÔÓÚµ¥µã¼ÆË㣬µ¥¾«¶È´¢´æ»ý·ÖÍêȫûÓÐÎÊÌ⣬¶Ô½á¹ûµÄ´øÀ´µÄÎó²î±ÊÕß·¢ÏÖÔÚ1E-8a.u.¼¶±ð£¬±ÈÃܶÈÄâºÏ½üËÆ´øÀ´µÄÎó²îСµÃ¶àµÃ¶à¡£²»¹ý¾Ý˵¸ö±ðÇé¿öÏÂÓõ¥¾«¶È´¢´æʱMP2ÌݶȺͳÚÔ¥µÄÃܶȶԸö±ðÌåϵÓÐÎó¡£ÏÔÈ»CFLOATÊÇ×îÊ¡ÄÚ´æµÄģʽ£¬¶øĬÈϵÄUCDOUBLEÊÇ×îºÄÄÚ´æµÄ¡£¶ÔÓÚСÌåϵÓÃĬÈϵÄUCDOUBLE¾ÍÐУ¬¶ÔÓÚ´óÌåϵ£¬Èç¹û´¢´æ×ÊÔ´½ôÕÅ»òÕß³¤Ê±¼ä¿ñ¶ÁдӲÅÌ£¬½¨ÒéÓÃCFLOAT£¨Ö±½ÓдÉÏÕâ¸ö¹Ø¼ü´Ê¾ÍÐУ©£¬¾¡¹ÜѹËõ/½âѹ¹ý³Ì»áºÄ·Ñ¼ÆËãÁ¿£¬µ«ÊÇÓÉÓÚ½µµÍÁËÓ²Å̶ÁдÁ¿£¬¶øÓ²ÅÌI/OËÙ¶ÈÓÖÊǼÆËã»úµÄÆ¿¾±£¬ËùÒÔ¶ÔÓÚÕû¸ö¼ÆËãÊÇÓмÓËÙЧ¹ûµÄ¡£

ͨ³£ÈÏΪ6-31+G*ÊǼÆËãÈõÏ໥×÷Óõĵ×ÏÞ»ù×飬Ôٵ;Íû¿ÉÐŶÈÁË£¬ÏñÊÇM062X/6-31+G*ÕâÖÖ×éºÏͨ³£À´ËµÊÇÇîµÃ²»ÄÜÔÙÇîµÄÈõÏ໥×÷ÓüÆË㼶±ð¡£¶øʹÓÃORCAµÄ»°£¬ÔÚÏÂÃ潫»á¿´µ½Ô­±¾ÕâÖÖ¼¶±ðËãÆðÀ´¶¼³ÔÁ¦µÄÇé¿ö¶¼¿ÉÒÔÏíÊÜdef2-QZVPP£¨-g,-f)¡£¶ø¸ü´óµÄÌåϵ£¬Ò²²»±Ø¾½ÆȵؾÀ½áÓÚÈçºÎÔÚÓÐÏ޵ļÆËãÄÜÁ¦ÏÂÒÔ×îºÃµÄ·½Ê½É趨ÃÖÉ¢º¯Êý£¬gCPÓÅÑŵذÚÍÑÁËÈõÏ໥×÷ÓüÆËã¶ÔÃÖÉ¢º¯ÊýµÄÒÀÀµÐÔ¡£



4 ʵÀý

ÕâÒ»½Ú½«Í¨¹ý¶à¸öʵÀýÕ¹ÏÖ²»Í¬¼ÆË㼶±ðϵĺÄʱ£¬ÌåϵÓÉ´óµ½Ð¡¡£¼ÆËãÌõ¼þ¾ÍÊÇÆÕͨµÄ¼ÒÓûúµÄÌõ¼þ£¬CPUΪIntel i7-2630QM(4ºË£¬×î¸ßƵÂÊ2.5G)£¬Ë«Í¨DDR3-1333 16GBÄڴ档ϵͳΪWin7-64bit¡£ORCAΪ3.0.1°æ£¬GaussianΪ09W A.02°æ£¬MOPAC2012Ϊ13.331W°æ¡£Gaussian¼ÆËãʱ¶¼ÓÃÁËSCF=conver=6¹Ø¼ü´Ê£¬ÒòΪG09ĬÈϵÄTightÊÕÁ²ÏÞ¹ýÑÏ£¬ËùÒÔÊʵ±½µµÍµãÒÔ½Úʡʱ¼ä£¨ÒÀÈ»×ã¹»¾«È·£©¡£¼ÆËãʱ¼äͳ¼ÆµÄÊÇʵ¼Ê×ܺÄʱ¶ø·ÇCPUʱ¼ä£¬ÈôδעÃ÷¶¼ÊÇ¿ªÆôÁËËĺ˲¢ÐС£¼ÆËãÈÎÎñ¶¼ÊÇËãÌåϵµÄµ¥µãÄÜ¡£

ÏÂÃæʵÀýÖÐÈç¹û²âÊÔÁ˶à¸ö¼ÆËã·½·¨£¬Ôò·½·¨Ç°Ãæ´òÐǺŵÄÊDZÊÕßÍƼöµÄÓÃÓÚÕâ¸öÌåϵµÄ¼ÆË㼶±ð¡£

4.1 Ë®ºÐ×Ó

Á½¸öË®ºÐ×ӵĽṹÈçÏÂ
´óÌåϵÈõÏ໥×÷ÓüÆËãµÄ½â¾öÖ®µÀ

Ë®ºÐ×Ó1£º12426atoms
PM7 MOZYME 4threads 14.4min£¬-282616.10434 KCAL/MOL£¬ÄÚ´æÓÃÁË1.2GB
Ë®ºÐ×Ó2£º25095atoms
PM7 MOZYME 4threads 98.6min£¬-576587.48001 KCAL/MOL£¬ÄÚ´æÓÃÁË4.5GB

¿É¼û¼´±ãÊÇÒ»Á½Íò¸öÔ­×ÓµÄÌåϵ£¬ÔÚÆÕͨ¼ÆËãÌõ¼þÏ£¬ÀûÓÃMOPAC2012½áºÏËüÌØÓеÄMOZYME¼¼ÊõÒ²ÕÕÑùÄÜÓ¦¸¶µÃÁË¡£Èç¹ûÄÚ´æ×ã¹»´ó£¬×öµ½¼¸Íò¸öÔ­×ÓµÄÌåϵҲÊÇÓпÉÄܵģ¬ÍêÈ«ÄܼÆËãÄÉÃ×¼¶±ðÌåϵµÄÈõÏ໥×÷Óá£


4.2 ETKµ°°×£¬3944¸öÔ­×Ó

´óÌåϵÈõÏ໥×÷ÓüÆËãµÄ½â¾öÖ®µÀ-1

PM7 MOZYME 1/4threads 323.8/325.5s -16765.64836 KCAL/MOL£¬ÓÃÁË300MBÄÚ´æ
¿É¼û¶ÔÓÚ¼¸Ç§¸öÔ­×ÓµÄÆÕͨ´óСµÄµ°°×£¬ÓÃMOPAC2012À´¼ÆËãÒ×Èç·´ÕÆ¡£Í¬Ê±Ò²¿´µ½ÔÚÄ¿Ç°µÄMOPAC2012 13.331W°æÀ¶ÔÓÚMOZYME¼ÆË㣬Óõ¥Ï̺߳Í4Ïß³Ìʵ¼ÊÉÏËٶȻù±¾Ã»ÓвîÒì¡£ÓÉÓÚ²¢Ðл¯Ä¿Ç°»¹Ã»´øÀ´Ê²Ã´Ð§¹û£¬ÎªÁË×î´ó³Ì¶ÈÀûÓüÆËã×ÊÔ´£¬Èç¹ûÄãÓÐn¸öºË£¬¿ÉÒÔͬʱËãn¸öÈÎÎñ¡£MOPACÔ­±¾ÊDz»Ö§³Ö²¢Ðеģ¬ÔÚ׫д±¾ÎÄʱMOPAC2012²Å¸Õ¸ÕÖ§³Ö²¢Ðв»¾Ã£¬¹À¼ÆÊdzÌÐòÏֽ׶λ¹Ã»°ÑMOZYMEµÄ²¢Ðл¯¸ãºÃ¡£

ÆäËü°ë¾­Ñé·½·¨µÄºÄʱ
PM6-DH+ MOZYME 4threads 349.6s -16405.40799 KCAL/MOL
PM6 MOZYME 4threads 342s -14168.08555 KCAL/MOL
PM3 MOZYME 4threads 330s -10619.64239 KCAL/MOL
RM1 MOZYME 4threads 295s -13358.60313 KCAL/MOL
AM1 MOZYME 4threads 271s -10033.36276 KCAL/MOL
¿É¼û£¬²»Í¬°ë¾­Ñé·½·¨ºÄʱ´óÖ¶¼ÔÚÒ»¸öˮƽÉÏ¡£´ÓÉú³ÉìÊÀ´¿´£¬PM6-DH+ºÍPM7µÄ½á¹ûºÜ½Ó½ü£¬È»¶øAM1ºÍPM3Ôò´ó´óµÍ¹ÀÁËÉú³ÉìÊ£¬Ôںܴó³Ì¶ÈÉϺÍËüÃDZíÏÖÈõÏ໥×÷ÓúܲîÓйء£¶øRM1ºÍPM6µÄ½á¹ûÔò½éÓÚAM1/PM3ºÍPM7Ö®¼ä£¬ÔÚ²àÃæÉÏ˵Ã÷RM1¡¢PM6¶ÔÈõÏ໥×÷ÓõıíÏÖÄÜÁ¦¿ÉÄÜÔÚAM1/PM3ºÍPM7Ö®¼ä¡£


4.3 DNA

ÕâÊǸöDNAƬ¶Î£¬×óͼ¹²380¸öÔ­×Ó£¬º¬6¸ö¼î»ù¶Ô¶ù£»ÓÒͼ¹²754¸öÔ­×Ó£¬º¬12¸ö¼î»ù¶Ô¶ù
´óÌåϵÈõÏ໥×÷ÓüÆËãµÄ½â¾öÖ®µÀ-2

754¸öÔ­×Ó
PM7 4threads ³¬¹ý7Сʱ¼ÆËãÉÐδÍê³É
PM7 MOZYME 1/4threads 58.88/57.8s -1161.85336 KCAL/MOL
¿É¼ûMOZYME¼«ÆäÖØÒª£¬²»ÓÃMOZYME¸ù±¾Ã»·¨Ëã¡£

389¸öÔ­×Ó
PM7 1/4threads 56.2s/34.6s -2463.90001 KCAL/MOL
PM7 MOZYME 4threads 10.9s -2463.88277 KCAL/MOL
*BLYP D3 GCP(DFT/SVP) def2-SVP def2-SVP/J  124min
ÕâÑù¼¸°Ù¸öÔ­×ÓµÄÌåϵÔÚMOPACÖв»ÓÃMOZYMEÖÕÓÚÒ²ÄÜËãÁË£¬²¢ÇÒ²¢Ðл¯ÖÕÓÚ¶Ô¼ÆËãËÙ¶ÈÓÐÌáÉý£¬µ«ÊDz¢ÐÐЧÂÊ»¹ÊǺܵ͡£¿ªÁËMOZYMEºó¼ÆËãʱ¼äÄÜÊ¡ºÃ¼¸±¶¡£MOZYME´øÀ´µÄÎó²î¿É¼û·Ç³£Ö®Ð¡£¬½öÓÐ0.117kcal/mol£¬¶ÔÓÚ¼¸°Ù¸öÔ­×ÓµÄÌåϵÍêÈ«¿ÉÒÔºöÂÔ²»¼Æ¡£
µÃÒæÓÚRI-J¼¼Êõ£¬ÔÚORCAÖУ¬BLYPÔÚdef-SVP»ù×飨´óСÓë6-31G**Ïà·Â·ð£©ÏÂÒ²ÄÜÈÝÒ׵ش¦ÀíÕâÑùºÃ¼¸°Ù¸öÔ­×ӵĴóÌåϵ£¡²¢ÇÒ½áºÏDFT-D3ºÍgCPУÕý£¬ÈõÏ໥×÷ÓüÆËã׼ȷ¶ÈºÜ²»´í¡£


4.4 DNAСƬ¶Î-С·Ö×Ó¸´ºÏÎ196¸öÔ­×Ó

´óÌåϵÈõÏ໥×÷ÓüÆËãµÄ½â¾öÖ®µÀ-3

PM7 MOZYME 4threads 2.5s
*BLYP D3 GCP(DFT/TZ) def2-TZVP def2-TZVP/J   91.4min
BLYP D3 GCP(DFT/SVP) def2-SVP def2-SVP/J   21.7min
¶ÔÕâÑùµÄÔ¼200¸öÔ­×ÓµÄÌåϵ£¬°ë¾­Ñé¼ÆËãÖ»ÐèÒªÒ»Õ£Ñ۵Ŧ·ò£¬×ö¼¸ºÎÓÅ»¯Ò²·Ç³£ÈÝÒס£BLYPÔÚdef2-TZVP¼¶±ðÒ²ÄܺÜÈÝÒ׵ؼÆË㣬½ö»¨ÁËÒ»¸ö°ëСʱ£¬ËäÈ»±ÈÆðÓÃdef2-SVP³¤Á˼¸±¶£¬µ«»¹ÊÇÖµµÃµÄ¡£


4.5 °ü¼Ð¸´ºÏÎ92¸öÔ­×Ó

Õâʵ¼ÊÉÏÊÇGrimmeµÄS12LÈõÏ໥²âÊÔ¼¯£¨Chem.Eur.J.,18,9955(2012)£©ÖеÄÌåϵ2a¡£
´óÌåϵÈõÏ໥×÷ÓüÆËãµÄ½â¾öÖ®µÀ-4

*BLYP D3 def2-QZVPP(-g,-f) def2-QZVPP/J  58.3min
BLYP D3 def2-QZVPP def2-QZVPP/J  116min
BLYP D3 GCP(DFT/TZ) def2-TZVP def2-TZVP/J  20.3min
B3LYP D3 GCP(DFT/TZ) def2-TZVP def2-TZVP/J RIJCOSX 126min
ÔÚBLYPÏ£¬¿É¼ûdef2-QZVPP(-g,-f)±Èdef2-QZVPP½ÚÊ¡ÁËÒ»°ëʱ¼ä£¬Ñ˸îµôµÄ¸ß½Ç¶¯Á¿º¯Êý¶ÔÓÚDFT¼ÆËã׼ȷ¶ÈµÄÓ°ÏìºÜС£¬ËùÒÔ½¨ÒéÓÃdef2-QZVPP(-g,-f)½áºÏBLYP¡£Èç¹û½µµÍÒ»¸ö¼ÆËãµµ´ÎÓÃTZVP£¬ÄǾ͸üÈÝÒ×ÁË£¬¿´ÍêÒ»»°Ð·¬£¨¶øÇÒ»¹ÌøOP¡¢ED£©£¬ÕýºÃ¾ÍËãÍêÁË¡£Õâô´óµÄÌåϵ£¬¿¿ORCAÔÚQZ¼¶±ðÏÂ×ö¼ÆËã½öÓÃ1¸öСʱ£¬¶øÒªÈÃGaussianË㻹²»ÖªµÀµÃËãµ½ºïÄêÂíÔÂÈ¥¡£
ÔÓ»¯·ºº¯B3LYP£¬¼´±ã¿¿RIJCOSX¼ÓËÙ£¬¼ÆËãºÄʱ»¹ÊDZÈÏàͬ»ù×éϵĿªÁËRI-JµÄ´¿·ºº¯¶àÁ˺ü¸±¶£¬TZ϶¼±ÈBLYPÔÚQZÏÂÂý£¬ËùÒÔ»®²»À´£¬ºÎ¿öÔÚD3+gCPУÕýÏÂB3LYP¼¸ºõ²¢²»±ÈBLYPÇ¿¡£


4.6 C60
B3LYP def2-TZVP(-f) def2-TZVP/JK RIJK 159.8min 11cyc -2285.47162975a.u.
*B3LYP def2-TZVP(-f) def2-TZVP/J RIJCOSX 47.6min 10cyc -2285.47245742a.u.
´ËÀýÖ÷ҪչʾһÏÂRIJCOSXºÍRI-JKµÄ²îÒì¡£¿É¼û¶ÔÓÚÕâÖÖ´óСµÄÌåϵ£¬RIJCOSXµÄ¼ÆËãËÙ¶ÈÒѾ­ÍêʤRI-JK¡£¶þÕß¼äÄÜÁ¿²îÓÐ0.00083£¨0.52kcal/mol£©£¬Ë­¸ü½Ó½ü׼ȷֵ²»ºÃ˵£¬²»¹ý¿Ï¶¨¶¼²îµÃ²»´ó¡£


4.7 GCGCËļî»ù¶Ñµþ£¬57¸öÔ­×Ó

´óÌåϵÈõÏ໥×÷ÓüÆËãµÄ½â¾öÖ®µÀ-5
6.png

RI-PWPB95 D3 RIJCOSX def2-QZVPP(-g,-f) def2-QZVPP/J def2-QZVPP/C grid4 tightSCF 134.1/405.8min
б¸ÜÇ°ÊÇSCF²½ÖèµÄºÄʱ£¬Ð±¸ÜºóÊÇ×ܺÄʱ¡£½«½ü60¸öÔ­×ÓµÄÌåϵ£¬Ë¯Ò»¾õµÄ¹¦·ò£¬ÔÚQZ¼¶±ðÏÂË«ÔÓ»¯·ºº¯¼ÆËã¾ÍÍê³ÉÁË£¬Õâ¸ö¼¶±ðϵÄÈõÏ໥×÷ÓüÆË㾫¶ÈÒѾ­ºÜÍêÃÀÁË¡£


4.8 ¶à»··¼Ìþ£¬37¸öÔ­×Ó

´óÌåϵÈõÏ໥×÷ÓüÆËãµÄ½â¾öÖ®µÀ-6

BLYP D3 def2-QZVPP(-g,-f) def2-QZVPP/J grid5 £¨Í¬Ê±Ð´ÉÏ%method FinalGrid 6 end£© 10.0min -996.4176595a.u.
BLYP D3 def2-QZVPP(-g,-f) def2-QZVPP/J grid4 8min  -996.4176418a.u. Error=0.01kcal/mol
BLYP D3 def2-QZVPP(-g,-f) def2-QZVPP/J 5.8min -996.4173819a.u.  Error=0.17kcal/mol
B3LYP D3 def2-QZVPP(-g,-f) def2-QZVPP/J RIJCOSX grid4 gridx4 56.4min -996.1652457a.u.
B3LYP D3 def2-QZVPP(-g,-f) def2-QZVPP/J RIJCOSX grid4 44min -996.1651966a.u. Error=0.03kcal/mol
B3LYP D3 def2-QZVPP(-g,-f) def2-QZVPP/J RIJCOSX 42.3min -996.1650320a.u. Error=0.13kcal/mol
B3LYP D3 def2-QZVPP(-g,-f) def2-QZVPP/JK RIJK grid4 35.2min -996.1646612a.u.
RI-PWPB95 D3 def2-QZVPP(-g,-f) def2-QZVPP/JK def2-QZVPP/C RIJK grid4 tightscf 39.5/67.3min -996.1346210a.u.
RI-PWPB95 D3 def2-QZVPP(-g,-f) def2-QZVPP/J def2-QZVPP/C RIJCOSX grid4 tightscf 57.2/87.6min -996.1351977a.u.
*RI-PWPB95 D3 def2-QZVPP def2-QZVPP/JK def2-QZVPP/C RIJK grid4 tightscf 57.5/99.3min -996.1694133a.u.

¿É¼ûBLYPÔÚQZ¼¶±ðÏÂËãÕâÖÖ´óСµÄÌåϵ·Ç³£ÇáËÉ£¬Ò²¾ÍÊÇBÕ¾¿´Á½¸öMADµÄ¹¦·ò¡£

²âÊÔÖбȽÏÁË»ý·Ö¸ñµã¶Ô½á¹ûºÍºÄʱµÄÓ°Ïì¡£BLYPµÄ¼ÆËãÒÔGrid5 + FinalGrid 6×÷Ϊ²Î¿¼Öµ£¬¼ìÑéÁËÆäËü¸ñµãÉ趨ϵÄÎó²î¡£¿É¼ûÔÚgrid4µÄʱºò½á¹ûÒѾ­·Ç³£¾«È·ÁË£¬Ã»±ØÒªÔÙÌá¸ßµ½grid5¡£ÓÃĬÈϵĻý·Ö¸ñµã±Ègrid4ʱ¿ì30%£¬µ«»áÒò´Ë´øÀ´Îó²î(0.16kcal)£¬µ¹ÊÇÍêÈ«ÊôÓÚ¿ÉÒÔ½ÓÊܵķ¶Î§¡£

¿ÉÒÔ¿´µ½B3LYPµÄºÄʱ±ÈBLYP¸ßÁËÒ»¸öÊýÁ¿¼¶¡£RI-JK¶Ô´óÌåϵµÄËÙ¶ÈÔ¶²»ÈçRIJCOSX£¬µ«ÊǶÔÓÚÔ­×Ó½ÏÉٵĴËÀý£¬RI-JK·´µ¹»¹Õ¼Á˵ãÓÅÊÆ¡£½¨ÒéÔÚ¼¸Ê®¸öÔ­×ӵķ¶Î§ÄÚʹÓÃRI-JK¡£ÓÉÓÚ´ËÀý±ÈC60ҪС£¬ËùÒÔRIJCOSXºÍRI-JK¼äµÄ½á¹û²îÒì¸üС£¬½öΪ0.23kcal/mol¡£±ê×¢µÄErrorÊÇÏà¶ÔÓÚgrid4 gridx4ʱµÄÎó²î¡£¿É¼ûRIJCOSXµÄ¸ñ×ÓÉ趨¶ÔÓÚ½á¹ûÓ°ÏìÉõ΢£¬ÓÃgridx4±ÈÓÃĬÈϵĸñ×ÓºÄʱÔö¼ÓÁ˲»ÉÙ£¬µ«Êǽá¹û²¢Ã»Ê²Ã´²îÒì¡£¶ÔÓÚÔÓ»¯·ºº¯µÄ¼ÆË㣬ÓÉÓÚÓÃgrid4¶îÍâÔö¼ÓµÄºÄʱÏà¶ÔÓÚ×ܺÄʱÉõ΢£¬ËùÒÔ½¨Òé×ÜÊÇдÉÏgrid4£¬ÌرðÊÇ´òËã×öË«ÔÓ»¯·ºº¯¼ÆËãʱ¡£

¶ÔÓÚË«ÔÓ»¯·ºº¯£¬ÓÃRI-JKÕÕÑù±ÈRIJCOSX¿ì¡£RI-JKʱ£¬Ê¹ÓÃdef2-QZVPP(-g,-f)»ù×飬SCFºÍ΢ÈŲ¿·Ö·Ö±ðÕ¼ÁË×ܺÄʱµÄ58.7%ºÍ41.0%£¬Î¢ÈŲ¿·ÖºÄʱ±ÈSCF²¿·Ö»¹µÍ¡£Õû¸öPWPB95¼ÆËãʱ¼ä±ÈB3LYP½ö¶àÁ˲»µ½Ò»±¶£¬¶ø¾«¶ÈÌáÉý²»ÉÙ£¬ËùÒÔÓÃORCA×öË«ÔÓ»¯·ºº¯ÕæÊdz¬»®Ë㣡ÏñÕâÑù¼¸Ê®¸öÔ­×ÓµÄÌåϵǿÁÒ½¨ÒéÔÚPWPB95/def2-QZVPP(-g,-f)¼¶±ðÏÂ×ö¡£°Ñdef2-QZVPP(-g,-f)ÌáÉýµ½def2-QZVPP¸ø×ܺÄʱ´øÀ´Á˽üÒ»°ëµÄÔö¼Ó£¬ËäÈ»¿´ËƶԾø¶ÔÄÜÁ¿Ó°ÏìÃ÷ÏÔ£¬Êµ¼ÊÉ϶ÔÏà¶ÔÄÜÁ¿Ó°ÏìÉõ΢¡£Èç¹ûÄãÊÇÍêÃÀÖ÷ÒåÕß²»Ï²»¶Ñ˸²»·ÁÓÃdef2-QZVPP¡£

ÓÃÕâ¸öÌåϵ£¬ÎÒÃÇÒ²±È½ÏÒ»ÏÂG09ºÍORCAµÄËٶȲîÒ졣б¸ÜÇ°ÊÇƽ¾ùÿÂÖµü´úµÄºÄʱ¡£ÏÂÃæÌÖÂÛµÄÒ»ÂÉÊÇƽ¾ùÿÂÖµü´úµÄʱ¼ä£¬ÒòΪӰÏìµü´ú´ÎÊýµÄÒòËظ´ÔÓ£¬×Üʱ¼ä²»Ò×¹«Æ½±È½Ï¡£
  Gaussian09 A.02£º
BLYP/def2-SVP 0.29/3.8min 13cyc -995.123921070
B3LYP/def2-SVP 0.48/5.3min 11cyc -995.505684
M062X/def2-SVP 0.59/6.5min 11cyc -995.110736706
BLYP/def2-TZVP 3.0/42.4min 14cyc -996.230052831
B3LYP/def2-TZVP 9.3/102min 11cyc -996.593814788
-----£¨À¨ºÅÄÚ´ú±íÃܶÈÄâºÏ´øÀ´µÄ¼ÓËٱȣ©
BLYP/TZVP 0.68/13.7min 20cyc -996.1872818
BLYP/TZVP/Auto 0.34/4.8min(2.0x) 14cyc -996.1887121 Error=0.90kcal/mol
BLYP/TZVP/TZVPfit 0.25/3.5min(2.7x) 14cyc -996.1908830 Error=2.26kcal/mol

  ORCA£¨noRI´ú±í±ÜÃâĬÈÏʹÓÃÃܶÈÄâºÏ£©£º
BLYP def2-SVP noRI 0.47/5.6min 12cyc -995.12384177
B3LYP def2-SVP noRI 0.5/5min 10cyc -994.92024998
M062X def2-SVP noRI 0.6/4.2min 7cyc -995.10510837
M062X def2-SVP noRI grid4 0.67/4.7min 7cyc -995.11062471
BLYP def2-TZVP noRI 6.5/84.7min 13cyc -996.22983135
B3LYP def2-TZVP noRI 7.7/77.3min 10cyc -996.00838771
-----£¨À¨ºÅÄÚ´ú±íÃܶÈÄâºÏ´øÀ´µÄ¼ÓËٱȣ©
BLYP def2-SVP def2-SVP/J 0.06/0.72min 12cyc(7.8x) -995.12468812 Error=0.53kcal/mol
BLYP def2-TZVP def2-TZVP/J 0.18/2.2min(36.1x) 12cyc -996.23055140 Error=0.31kcal/mol
BLYP def2-SVP def2-SVP/JK 0.06/0.73min 12cyc(7.8x) -995.12417598 Error=0.21kcal/mol
BLYP def2-TZVP def2-TZVP/JK 0.2/2.4min(32.5x) 12cyc -996.23023270 Error=0.25kcal/mol
B3LYP def2-SVP def2-SVP/J RIJCOSX 0.25/2.5min(2.0x) 10cyc -994.9211858 Error=0.59kcal/mol
B3LYP def2-SVP def2-SVP/JK RIJK 0.34/4.1min(1.5x) 12cyc -994.9203942 Error=0.09kcal/mol
B3LYP def2-TZVP def2-TZVP/JK RIJK 0.98/11.8min(7.8x) 12cyc -996.0085860 Error=0.12kcal/mol

ÎÒÃÇÏÈ¿¼²ì²»Ê¹ÓÃÃܶÈÄâºÏʱBLYP¡¢B3LYP¡¢M06-2XÔÚdef2-SVPºÍdef2-TZVPϵÄƽ¾ùÿÂÖµü´úʱ¼ä£¬ORCA¶ÔÓÚ´¿·ºº¯ËÙ¶ÈÂýÓÚGaussian²»ÉÙ£¬¶øÔÓ»¯·ºº¯ºÍmetaÔÓ»¯·ºº¯ÏÂÁ½¸ö³ÌÐòËÙ¶ÈÔÚ²®ÖÙÖ®¼ä¡£¼ÆËãºÄʱÏÔÈ»¶¼ÊÇM062X>B3LYP>BLYP¡£ORCAÏÂʹÓÃmetaÔÓ»¯·ºº¯¼ÆËãʱÎñ±Ø¼ÓÉÏgrid4£¬·ñÔòÓÉÊý¾Ý¿É¼ûÎó²îÄÜ´ïµ½3kcal/molµÄ³Ì¶È£¬ÕâÊÇÒòΪÖÚËùÖÜÖªmetaÔÓ»¯·ºº¯¶Ô¸ñµãÃô¸ÐÐԽϸߣ¨¶øG09ËùÓõĻý·Ö¸ñµã¾«¶È±ÈORCA¸ß²»ÉÙ£¬ËùÒÔ²»ÓÃ×öÌØÊ⿼ÂÇ£©¡£ORCAºÍGaussianµÄBLYP¼°M062X½á¹ûÏà·û³Ì¶È½ÏºÃ£¬µ«B3LYPÏàÈ¥ÉõÔ¶£¬ÕâÊÇÒòΪORCAµÄB3LYPÓõÄÊÇ25%µÄHF³É·Ý¶ø·ÇÔ­ÎĶ¨ÒåÒÔ¼°GaussianÖÐËùÓõÄ20%³É·Ý¡£ORCAÀïдB3LYP/G²ÅÊÇÓÃB3LYPµÄԭʼ¶¨Ò壬´Ëʱdef2-SVPϵĽá¹ûΪ-995.5056201£¬ºÍGaussianµÄ½á¹ûºÜÏà·û¡£

±È½Ï¿É¼û£¬ÔÚGaussianÖÐÓÃB3LYP/def2-TZVPËãÕâ¸öÌåϵ£¬Ã¿ÂÖµü´úµÄʱ¼ä¶¼×ã¹»ORCAÆôÓÃRI-Jʱ×öÒ»´ÎBLYP/def2-QZVPP(-g,-f)ÈÎÎñ£¡Õû¸öB3LYP/def2-TZVPµÄºÄʱÔÚORCAÖÐ×öÒ»´ÎPWPB95/def2-QZVPP¶¼´Â´ÂÓÐÓ࣡ͬʱ¿É¼û£¬Ëæ×Å»ù×é´Ódef2-SVPÔö´óµ½def2-TZVP£¬GaussianµÄ×ܺÄʱÔö¼ÓÒ»¸ö¶àÊýÁ¿¼¶£¬ÈôÈÃGaussian×öB3LYP/def2-QZVPP£¬ºÄʱ¿ÉÄܶ¼¹»ORCAÅܼ¸Ê®´ÎPWPB95/def2-QZVPPÁË£¡

ÏÂÃæÔÙÀ´¿¼²ìÃܶÈÄâºÏ¸ø¼ÆËãËٶȺͽá¹û´øÀ´µÄÓ°Ï졣ͨ³£´ó¼Ò²»ÔõôÌáGaussianµÄÃܶÈÄâºÏ£¬Ö÷ÒªÔ­ÒòÊÇGaussianÔÚÕâ·½Ãæ±È½Ï²Ð¼²£¬Ö§³ÖµÄÃܶÈÄâºÏ»ù×éÌ«ÉÙ£¬¶øÇÒ¼ÓËÙ²»Ã÷ÏÔ£¬ÁíÍâÒ²Ö»Ö§³ÖRI-J£¬¶ø²»Ö§³Ö¼ÓËÙ½»»»²¿·ÖºÍMP2²¿·Ö¡£GaussianÖÐ×Ô´øµÄTZVPÃܶÈÄâºÏ»ù×éÊǸøÀÏ°æ±¾TZVPÓõģ¬Ð´³ÉBLYP/TZVP/TZVPfitÐÎʽÀ´Ê¹Ó㬱ÈÆðBLYP/TZVP¼ÓËٱȽö½öÊÇ2.7x£¬ËäÈ»ÓÐЧ¹û£¬µ«ÊDZȽÏÓÐÏÞ£¬¶øÇÒÒýÈëµÄÎó²î´ïµ½2.26kcal/mol£¬ÊDz»¿ÉºöÂԵġ£GaussianÒ²Ö§³Ö×Ô¶¯Éú³ÉÃܶÈÄâºÏ»ù×飬дÉÏAuto¼´¿É£¬Ç°ÎÄÌáµ½¹ýÕâÑù×Ô¶¯Éú³ÉµÄÃܶÈÄâºÏ»ù×é½Ï´ó£¬ËùÒÔ¼ÓËٱȽµµ½2.0x£¬µ«ÊǾ«¶È±È±êÅäµÄTZVPfitºÃ£¬ÒýÈëµÄÎó²îΪ0.9kcal/mol¡£

ÔÙÀ´¿´ORCAµÄÃܶÈÄâºÏ¡£RI-JʹµÃʹÓÃ/J¸¨Öú»ù×éÏÂBLYP/def2-SVPµÄ¼ÆËãËٶȴﵽnoRIʱµÄ7.8±¶£¬¶ÔÓÚdef2-TZVP£¬¼ÓËÙ±ÈÔò´ïµ½ÁË36±¶£¡ºÍÔ­ÏÈÒѾ­ÍêÈ«²»ÔÚÒ»¸ö¾³½çÁË¡£¶ÔÓÚB3LYP£¬ÔÚʹÓÃdef2-SVP»ù×éʱ£¬RI-JKºÍRIJCOSX¶¼Ã»¸ø¼ÆËãËٶȴøÀ´ÖʵÄÌáÉý£¬È»¶øÔÚdef2-TZVP»ù×éʱ£¬ÌáÉýÏÔÖø£¬ËÙ¶ÈÊÇÔ­ÏȵÄ7.8±¶£¬µ«±ÈÆðRI-JµÄÌáÉý·ù¶È»¹ÊDzîµÃºÜÔ¶¡£Ô­ÒòÔÚ֮ǰÌÖÂÛ¹ý£¬ÃܶÈÄâºÏ¼¼Êõ¼ÓËÙ½»»»²¿·Ö±ÈÆð¼ÓËÙ¿âÂز¿·ÖÄѵöࡣÁíÍâ²âÊÔ½á¹ûÒ²Çå³þÌåÏÖ³ö»ù×éÔ½´óÃܶÈÄâºÏ´øÀ´µÄÊÕÒæÔ½ÏÔÖø£¬ÕâÔÚÇ°ÃæÒ²ÒѾ­ÌÖÂÛ¹ý¡£

ÃܶÈÄâºÏÒýÈëµÄÎó²îºÍÌåϵ´óСÓÐÒ»¶¨¹Øϵ£¬¶ÔÕâ¸öÌåϵORCAµÄRIÒýÈëµÄÎó²î<0.6kcal/mol£¬ÊÇÍêÈ«¿ÉÒÔ½ÓÊܵģ¨±ðÍüÁËʵ¼Ê¹ØÐĵÄÊÇÏà¶ÔÄÜÁ¿£¬ÄܵÖÏûºÜ¶à£©£¬Õâ±ÈÆðGaussianµÄÃܶÈÄâºÏµÄÎó²îСµÃ¶à¡£´ËÀýRI-JKµÄÎó²î¿´ÆðÀ´ÉõС£¬¼ÈÓÐÇɺϳɷݣ¬Ò»¶¨³Ì¶ÈÉÏÒ²µÃÒæÓÚ¼ÆËã¿âÂز¿·ÖÓÃÁ˱È/J¸ü´óµÄ/JK¸¨Öú»ù×飬ÖÁÉÙÖ¤Ã÷ÁËRI-JKÊǺܿɿ¿µÄ£¬Ôڸ߾«¶ÈÈõÏ໥×÷ÓüÆËãÖÐʹÓÃRI-JKÊÇ¿ÉÒÔ·ÅÐĵġ£×¢Òâ¶ÔÓÚBLYPµÄ¼ÆË㣬ÈçÊý¾ÝËùʾ£¬Èç¹û°Ñ/J¸¨Öú»ù×é¸Ä³É¸ü´óµÄ/JK£¬RI-JËäÈ»¶àºÄʱһµãµã£¬µ«ÊÇÎó²î¼õСÁË¡£Ò²¼´ÊÇ˵£¬RI-JºÍRIJCOSX¼ÆËãÖÐÈôÏë¼õÉÙÎó²î£¬Ê¹ÓñÈ/J¸ü´óµÄ/JKÊǸö´ÏÃ÷µÄ°ì·¨£¬¾¡¹Ü/JKÔ­±¾ÊÇΪÁË×öRI-JK²ÅÌá³öµÄ¡£


5 ×ܽá

´óµÖ´Ó2005ÄêÆð£¬ÈõÏ໥×÷ÓüÆËã¼¼ÊõµÄ·¢Õ¹¾ÍÔÚ·ÉËÙ·¢Õ¹ÖУ¬¶øÇÒ»¹ÔÚ²»¶Ï¼ÓËÙ£¬benchmarkÎÄÕÂÒѾ­ÂúÌì·É£¨ÔÚJCTCÉÏÓÈÉõ£¬ÇÒ¶à°ëÊǸߵµË®ÎÄ£©£¬ÌرðÊÇDFT-DÊܹØ×¢¶È¼«¸ß¡£µ«ÊÇ£¬³ýÁËDFT-D£¬ÕýÈç±¾ÎÄËù½éÉÜÌÖÂ۵ģ¬»¹ÓкܶàÆäËüÀíÂÛ·½·¨¡¢ÊýÖµ¼¼ÊõºÍ³ÌÐòÉϵĽøÕ¹¶¼¸øÈõÏ໥×÷ÓüÆËãµÄÁìÓò´øÀ´Éî¿Ì±ä¸ï£¬ÖØÒª³Ì¶ÈË¿ºÁ²»ÑÇÓÚDFT-D¡£±¾ÎÄÔÚµÚ3½Ú¸ø³öµÄÍƼöÊDZÊÕßÈÏΪÔÚ׫ÎÄʱ×î¼ÑµÄÓÃÓÚ¸÷Öֳ߶ÈÌåϵµÄÈõÏ໥×÷ÓüÆËã·½·¨¡£Ïë±ØÔÚ¼¸ÄêÖ®ºó£¬ÓÖÓиü¶à¸üºÃµÄ·½·¨¡¢¼¼Êõ±»Ìá³ö²¢ÔÚÖ÷Á÷³ÌÐòÖÐʵÏÖ£¬±¾ÎÄÍƼöµÄ·½·¨½ìʱ¿Ï¶¨»¹ÐèÒª¸ù¾Ýʵ¼ÊÐÎÊƽøÐÐÒ»¶¨µ÷Õû¡£

[ Last edited by sobereva on 2013-12-18 at 10:15 ] ·µ»ØСľ³æ²é¿´¸ü¶à

½ñÈÕÈÈÌû
  • ¾«»ªÆÀÂÛ
  • wzwanglab

    ·ÖÎöµ½Î»£¬Ð´µÄºÜºÃ¡£
    Ñо¿·Ö×Ó¼äÈõÏ໥×÷ÓÃÌåϵÀë²»¿ªSAPT£¬×î½ü·¢Õ¹µÄSAPT0ºÅ³Æ¿ÉÒÔ¼ÆËã200¸öÔ­×Ó¡£
    PSI×÷Õߣ¬SherrillÊÇJCP±à¼­£¬T. Daniel Crawford ÊÇJPCA±à¼­£¬¶¼ÊÇÅ£ÈË¡£

  • hlq3175

    Ì«Å£ÁË£¬Ñ§Ï°Ð¶«Î÷ÁË£¬Ð»Ð»ÀÏʦ

  • valenhou001

    ·Ç³£Ð»Ð»ÄãµÄºÃÎÄ£¡

  • lihb734

    ÒýÓûØÌû:
    4Â¥: Originally posted by valenhou001 at 2013-12-18 08:21:47
    ·Ç³£Ð»Ð»ÄãµÄºÃÎÄ£¡

    ºîÀÏʦ¶¼À´¿´ÁË£¬ÄãÁ©ÊôÓÚÇ¿Ç¿¶Ô»°°
    £¬

  • lifei2008

    O(¡É_¡É)OллsobÀÏʦ£¬ÓÖѧµ½Ð¶«Î÷ÁË£¡

  • ÂþÌìƮѩ

    Â¥Ö÷Ò»³öÂí£¬¾ÍÊÇÉñÌù£¡ÍƼöÖ¸Êý5ÐÇ¡£

  • fish.yfyh

    ÇëÎÊSob£¬DFT-DÒÔ¼°DFT-D3ÕæµÄ¾Í×öºÃÈõÏ໥×÷ÓõÄУÕýÁËÂð£¿ÔÚ¸´ÔÓÌåϵÖУ¬Ö§Á´½Ï¶àµÄÇé¿öÏ£¬Ò»Î¶Ö»Ç¿µ÷dispersionµÄµØ룬»á²»»á³öÏÖÈԸ㲻¶¨µÄÇéÐÎÄØ£¿

²ÂÄãϲ»¶
ÏÂÔØСľ³æAPP
Óë700Íò¿ÆÑдïÈËËæʱ½»Á÷
  • ¶þάÂë
  • IOS
  • °²×¿