计算化学公社

标题: 简谈量子化学计算中DFT泛函的选择 [打印本页]

作者
Author: sobereva 时间: 2014-12-29 12:27
标题: 简谈量子化学计算中DFT泛函的选择

简谈量子化学计算中DFT泛函的选择

A brief discussion on the choice of DFT functional in quantum chemistry calculations

文/Sobereva @北京科音

First release: 2014-Dec-29 Last update: 2024-Oct-23

1 前言

经常有人问怎么选择DFT泛函，或要求推荐一些专门讲各种泛函特点、适用场合的文章。这比较为难，要想详细剖析各类泛函的特点，那得长篇大论一番，几句话根本说不清楚。至于推荐文献，也真没什么和适合推荐的。虽然泛函的综述很多，专著也很多，但是极少有全面评述的，而且写这种文章的人多半是研究泛函且提出过泛函的，观点、内容取材总是很有偏向性，有失客观，一个人一个说法。而且新泛函不断出现，稍老一点的文章就会和主流形势相脱节，甚至对读者产生误导。所以，真想对泛函的选择搞得十分灵通的话，必须认真读读DFT专著打好基础，并且多关注主流文献，逐渐形成自己的认识。光是看几篇综述根本没用，看10篇也不够。这就像测试数据，测一次或几次可能会有很大误差，但测试得多了，平均值就会很接近实际值。

初学者切勿盲目效仿文献里对泛函的选用。已发表的文章里有大量是使用不合理的泛函的，也有很多文献过于老旧，完全脱离时下潮流。效仿这些文献容易吃大亏，被审稿人猛批。初学者更不要瞎做所谓的泛函的“标定”，这只会白浪费时间，结论还可能是荒谬的，注意看《谈谈量子化学研究中什么样的benchmark才有意义》（http://sobereva.com/554）。

下面罗列一下对于各种情况，结合笔者计算经验和大量测试文章给出的泛函选择建议。本文会不断更新以与时俱进、体现笔者最新的看法。这里只涉及量子化学计算，第一性原理计算不涉及。本文说的“有机体系”确切来说是指的由主族原子构成的体系，因此也包括含有碱、碱土金属的情况。而“过渡金属体系”在本文里也包含镧系、锕系的情况。本文为了打字省事，所有ω都写成了w。本文里管双杂化泛函以外的泛函统称为普通泛函。本文里凡是带DFT-D3校正的情况，如果你的程序支持DFT-D4，而且对相应的泛函定义了参数，那么都可以改用DFT-D4，对于有机类体系结果相对于DFT-D3不会有太大差别，但是对于计算过渡金属配合物的金属有机反应、弱相互作用的情况可能精度会有明显的改进。下文推荐的很多泛函不一定在你当前用的程序里能用，能用哪些仔细看手册的DFT部分，有些在Gaussian中不直接支持的泛函通过特殊方式可能可以用，参见《Gaussian中非内置的理论方法和泛函的用法》（http://sobereva.com/344）。

由于很多初学者缺乏一些常识，这里先明确一下，计算能量相关问题的精度的顺序通常是：CCSD(T)>双杂化泛函>恰当选用的普通泛函>MP2>HF>半经验方法。
这些方法的理论方面的知识介绍不在本文范畴内。在北京科音初级量子化学培训班（内容见http://www.keinsci.com/workshop/KEQC_content.html）里我会对本文涉及方法的背景知识进行系统性的介绍使读者更充分地了解他们。在北京科音中级量子化学培训班（内容见http://www.keinsci.com/workshop/KBQC_content.html）里我会进一步展开做更深入全面透彻的讲解，给出更丰富、更具体的推荐，非常欢迎参加。

2 双杂化泛函的选择

能用得起双杂化泛函而且要求计算精度较高时通常建议用双杂化泛函。这类泛函精度大多比普通泛函能提升一个档，更是甩MP2一条街，而耗时比MP2仅略高一丁点。不过，目前一些较新的普通泛函（如wB97M-V）的精度在一些情况下的精度能赶上甚至超越双杂化泛函中比较老、比较中庸的那些，所以不要简单认为双杂化泛函的精度就一定比普通泛函更好。各种双杂化泛函的精度好坏没有定论，和体系、计算的问题有极大关系，不同测试文章的结论出入很大。即便不算弱相互作用也应当加DFT-D3(BJ)以改善综合性能。双杂化泛函做几何优化，特别是振动分析，耗时远高于普通泛函，而在这方面比合适的普通泛函的精度优势并不大，所以若无特殊必要不要用双杂化泛函做这种任务，而适合在普通泛函做优化、振动分析后算单点。

值得注意的双杂化泛函：
• PWPB95-D3(BJ)：稳健性高，整体精度也较好。对于计算金属-有机反应尤为适合优先考虑。算分子间弱相互作用在双杂化里属于排名较前的。ORCA等支持，Gaussian没法用
• B2PLYP-D3(BJ)：历史最老，精度在双杂化泛函里平庸，甚至有时已经落后于后来发展的较好的普通泛函如wB97M-V。但被程序支持最广泛，稳健性还不错。如果是Gaussian用户的话仍可以考虑，但如果是ORCA用户的话建议则不建议用，因为PWPB95-D3(BJ)明显更好。
• DSD-PBEP86-D3(BJ)和DSD-BLYP-D3(BJ)：在一些测试文章中表现很好，但也有不少令人大跌眼镜的时候，故使用风险大。G16已内置了DSD-PBEP86-D3(BJ)，直接写DSDPBEP86即可
• revDSD-PBEP86-D3(BJ)：2019提出，是DSD-PBEP86-D3(BJ)的改良，精度号称有不小提升，但是否对各种情况都表现得稳健还不好说。对于计算弱相互作用能、主族体系反应能来说是很推荐使用的，ORCA 6.0开始直接内置了。在Gaussian 16里截止到C.01版尚没有内置，但可以通过自定义泛函方式使用，见《Gaussian中非内置的理论方法和泛函的用法》（http://sobereva.com/344）。这是Gaussian用户的首选
• wB97X-2-D3(BJ)：整体精度是双杂化泛函里非常好的，也很稳健，不过算分子间弱相互作用很一般（分子内弱相互作用算得则很好）。wB97X-2仅有ORCA、GAMESS-US、Q-Chem支持，至少ORCA 5.0和GAMESS-US里没直接内置其DFT-D3(BJ)色散校正参数，得自己去PCCP, 20, 23175 (2018)的补充材料里获得并在程序中定义才能用，ORCA用户可以按照http://sobereva.com/490说的直接用Multiwfn创建此泛函的输入文件
• wB97M(2)：2018年提出，是双杂化泛函中性能最优秀的，ORCA 6.0、从Q-Chem 5.2后期版本开始支持

关于用双杂化泛函算激发态：局域激发和分子内电荷转移激发用DSD-PBEP86做TDDFT计算最理想，ORCA直接支持。分子间电荷转移激发用RS-PBEP86双杂化泛函基于SOS-CIS(D)形式计算最佳，这对于局域激发和分子内电荷转移激发描述得也较好，但目前只有MRCC程序支持。ORCA里算分子间电荷转移激发建议基于RSX-QIDH双杂化泛函做TDDFT，虽然略逊于RS-PBEP86/SOS-CIS(D)但远强于DSD-PBEP86做TDDFT。这些结论来自JCTC (2022) DOI: 10.1021/acs.jctc.1c01307的测试。

3 普通泛函的选择

对于普通泛函的情况，一般情况笔者推荐的选择如下：
• 计算垂直局域价层激发能：PBE0（指结合TDDFT，后同。平均误差约在0.3 eV，平均好于B3LYP一丝）
• 计算大共轭体系激发能、电荷转移/里德堡垂直激发：CAM-B3LYP、wB97XD。大共轭体系也可以考虑PBE38
• 计算激发态是三重态的激发能：M06-2X
• 计算各类激发能，只用一个泛函的情况，尝试：PBE0、CAM-B3LYP、PBE38、wB97XD
注：特别推荐、非常稳、普适性极好的是用w调控后的LC-wPBE，是对每个体系的每个结构专门优化最佳的w参数，见《优化长程校正泛函w参数的简便工具optDFTw》（http://sobereva.com/346），在《北京科音高级量子化学培训班》（http://www.keinsci.com/workshop/KAQC_content.html）里“交换-相关泛函的调控”一节做了特别全面深入的讲解并给了例子，但代价是需要额外花一些时间优化w参数。此外，用ORCA做基于DSD-PBEP86等双杂化泛函的TDDFT的精度和普适性也比用前述那些普通泛函整体好得多，值得尝试，但明显更耗时、更费内存。
• 计算极化率：PBE0或BHandHLYP（大共轭体系建议用后者）。双杂化里DSD-PBEPBE算极化率最佳，仅次于CCSD。
• 计算超极化率：wB97或LC-τHCTH。BHandHLYP和CAM-B3LYP也可以考虑。
注：计算极化率和第一超极化率几乎最理想的是用前述的w调控后的LC-wPBE。
• 计算NMR：对碳谱、氢谱，推荐用B3LYP结合标度法计算，见《谈谈如何又好又快地计算NMR化学位移》（http://sobereva.com/354）。如果不用标度法而想用普通方法的话，用revTPSS结合pcSseg-1，看《revTPSS泛函结合pcSseg-1基组是计算NMR很好的选择》（http://sobereva.com/623）
• 计算分子振动（谐振、非谐振都包括）：B3LYP。对于色散作用严重影响体系构象的情况，用B3LYP-D3(BJ)
• 主族元素构成的体系结构优化：B3LYP。对于色散作用严重影响体系构象的情况，必须加DFT-D3(BJ)校正。如果是大范围显著共轭体系，用wB97XD更好。此文务必阅读：《谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的》（http://sobereva.com/557）
• 计算HOMO-LUMO gap：B3LYP。更多信息见《正确地认识分子的能隙(gap)、HOMO和LUMO》（http://sobereva.com/543）
• 计算有机体系热力学数据（反应能、构型能量差、能垒等）：首选wB97M-V，其次M06-2X
• 计算各种弱相互作用：wB97M-V>=wB97X-V>=PW6B95-D3(BJ)>=M06-2X（加D3后更稳健）≈wB97X-D3>=B3LYP-D3(BJ)。在ORCA中追求速度用r2SCAN-3c（可利用RIJ显著加速计算）
• 计算卤键：M06-2X
• 主族元素形成的原子团簇：PBE0或TPSSh。碳团簇用B3LYP也合适
• Cu、Ag、Au团簇：TPSS或TPSSh。用revDSD-PBEP86双杂化泛函明显更好
• 计算多参考特征很强的体系：r2SCAN、TPSS、MN15L
• 计算有机体系垂直电子亲和能、Fundamental gap：wB97X或wB97XD（用w调控的泛函更佳）
• 计算有机体系垂直电离能：M06-2X
• 计算电子密度分布：PBE0、TPSSh
• 计算偶极矩、多极矩、静电势：wB97M-V。TPSSh、B3LYP也还行。revDSD-PBEP86-D3(BJ)双杂化泛函表现明显更好

涉及过渡金属、镧系锕系配合物有关的情况：
• 闭壳层过渡金属-有机反应能和势垒（涉及到配位键的形成/断裂、多参考特征不很强的情况）：首选wB97X-V或wB97M-V，其次PBE0-D3(BJ)和MN15。TPSSh-D3(BJ)也不赖。必须用纯泛函的话，用B97M-rV或差一点的B97M-D3(BJ)，再其次是SCAN-D3(BJ)或TPSS-D3(BJ)。用r2SCAN-3c组合方法的话不仅便宜且精度大概率强于任何纯泛函+色散校正且结合大基组
• 开壳层过渡金属-有机反应能和势垒（涉及到配位键的形成/断裂、多参考特征不很强的情况）：首选TPSS0-D4，其次PBE0-D4（PBE0-D3(BJ)稍逊一丝）。纯泛函当中B97M-V表现最好，但用r2SCAN-3c组合方法的话不仅便宜且精度还强于B97M-V/def2-QZVPP
• 过渡金属-配位键键长优化：
3d金属：适合纯泛函或低HF成份泛函。首选TPSS或TPSSh。r2SCAN也不错
4d、5d金属：适合中等HF成份的泛函。PBE0最优
• 过渡金属/镧系锕系配合物的杂问题：PBE0（镧系优先考虑这个）、TPSSh（锕系优先考虑这个）、MN15、MN15L、r2SCAN都可以考虑。如果没有实验参照的话，建议同时用其中几个泛函都算算，尽量都加DFT-D3(BJ)或DFT-D4校正。算单点用wB97X-V或wB97M-V可能结果明显更好。算激发能不宜用纯泛函，应当用前文提到的其它的。B3LYP也不是不能用于算这类体系，但非常不推荐，因为表现很平庸
• 计算过渡金属之间通过d轨道形成的化学键（包括过渡金属团簇体系）：由于静态相关非常强，必须用纯泛函，可考虑r2SCAN、TPSS、PBE。MN15L也可考虑
• 主族原子桥联过渡金属情况（如Fe-S-Fe）：TPSSh
• 计算配合物不同自旋多重度能量差：MN15L（见JCTC,16,4416）、TPSSh

上面没列举的其它情况或模棱两可的时候首先尝试B3LYP-D3(BJ)，若发现其结果不理想，可尝试wB97M-V、MN15。当然，如果能接受更大计算量，最好用双杂化泛函。

关于激发态、弱相互作用、NMR的计算笔者专门写过文章，上面只是十分简略地一说，强烈建议看看以下文章了解详情。

关于激发态计算：
乱谈激发态的计算方法
http://sobereva.com/265
几篇横测电子激发能计算精度的文章
http://bbs.keinsci.com/thread-470-1-1.html

关于弱相互作用计算：
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413
DFT-D色散校正的使用
http://sobereva.com/210
DFT-D4色散校正的简介与使用
http://sobereva.com/464
乱谈DFT-D
http://sobereva.com/83
大体系弱相互作用计算的解决之道
http://sobereva.com/214

关于NMR计算：
谈谈如何又好又快地计算NMR化学位移
http://sobereva.com/354
使用Multiwfn绘制NMR谱
http://sobereva.com/565

4 对一些泛函的评论

顺便说几句B3LYP。这把老骨头，从1994年到现在已过了30多年，依然还是量化界用得最多的泛函，也是多数人默认用的泛函，在《2021年计算化学公社论坛“你最常用的计算化学程序和DFT泛函”投票结果统计》（http://sobereva.com/599）的统计中就可看出来。尽管有很大程度上是因为大家用习惯了，但关键还是B3LYP整体来说很皮实，鲁棒性好。虽然目前有几百种泛函，不少都能在自己的专长的问题上秒杀B3LYP，但是很少有泛函在普适性上显著超过B3LYP，因此整体流行度、知名度和B3LYP也根本没法比。B3LYP主要有两个关键性软肋 1：范德华作用中的色散吸引部分完全没法描述，但加上色散校正如DFT-D3后可完全解决，这个校正也使B3LYP在热力学数据计算精度上有少量改进 2：电荷转移(CT)激发、里德堡激发计算很烂，但变体CAM-B3LYP也大体解决了这个问题（CAM-B3LYP用在其它场合往往明显比B3LYP差）。另外，B3LYP还有些其它相对次要的问题，比如静态相关强的体系描述不理想（杂化泛函的通病）、对某些过渡金属配合物高估了高自旋态稳定性（将HF成份从20%调低至15%可解决，见JCP,117,4729、TCA,107,48）、某些化学反应计算不好（如F+H2->HF+H没有预测出势垒）、对于大范围共轭体系有把体系描述成平面而且键长均等的倾向（详见http://sobereva.com/557）。但是，泛函无完泛函（除非是传说中的圣杯泛函，就算找到了估计也根本算不动），不能苛求B3LYP尽善尽美，或者说，B3LYP在一定程度上粗略地达到了一个平衡点，即如果再专为了改进A性质的精度而优化参数，那么B性质的计算精度就可能会下降。总之，除了色散作用主导的弱相互作用和CT/里德堡激发这两个绝对不能用B3LYP的问题以外，碰到一个新问题，若不知道用什么泛函合适同时又懒得去查文献找最适合的泛函，那么大可先用B3LYP算算试试。

作为B3LYP这种准万能泛函的潜在替代品有二，一个是M06-2X，一个是wB97XD，都是2008年诞生的，主流的量化程序如Gaussian也都支持。它俩整体都还不错，且不仅没有B3LYP那两根关键性软肋，这软肋反倒还成了它们的长处。M06-2X是通过拟合参数时考虑弱相互作用体系使得弱相互作用描述颇不错，而且由于HF成份高达54%，使得算CT、里德堡激发算得不错。wB97XD使用了DFT-D2色散校正使得弱相互作用计算精度不错，又引入了长程校正使得CT、里德堡激发也算得不错。这两个泛函都有个通病，就是计算速度慢，比B3LYP慢不少（不过，B3LYP算得动的，这两个泛函肯定也能算得动）。而且，M06-2X等明尼苏达系列泛函对DFT积分格点精度要求比普通泛函高很多，有时候因为格点不够精细导致几何优化和SCF不容易收敛、优化之后出现小虚频。提高积分格点精度可解决（如Gaussian 09里用int=ultrafine。极个别时候，诸如某些弱相互作用体系甚至需要int=superfine），但会明显增加耗时。另外，M06-2X是针对主族进行参数化的，对于含过渡金属、镧系锕系的体系非常不合适（几乎等于找死。除非研究的问题只涉及配体部分），这种问题非要用明尼苏达系列泛函的话可改用后来的MN15或MN15L。wB97XD和M06-2X一样也是给主族参数化的，但计算含过渡金属的体系却不至于像M06-2X总是那么差。由于wB97XD和M06-2X的HF成分都比B3LYP高得多，故对于静态相关的描述也糟得多，而且会非常显著高估HOMO-LUMO gap，对单重态局域价层激发能也产生了偏高的倾向。另外，根据JCTC,6,2115的测试，M06-2X计算非谐振频率极烂，wB97X也不咋地（因此其姊妹wB97XD这方面也应该不怎么样）。由于这些原因，特别是计算速度偏慢，导致它们还没法真正替代B3LYP成为默认的泛函，只能是在有确切需要的时候才用。M06-2X的一个鲜明优点是在有机体系反应能、异构化能、势垒等热力学量方面的计算精度很好，有时候接近双杂化泛函的精度（往往比一些垫底的双杂化泛函表现还好），这也跟拟合参数时引入了这类问题直接相关。总的来说，算有机体系热化学问题，用M06-2X的话可以不必担心被审稿人质疑。

M06系列明尼苏达系列泛函提出后，虽然又提出了很多新的明尼苏达系列泛函，但对于分子体系，精度都没有实质性的进步（泛函构建思想不行，就知道拟合来拟合去，显然很快就会陷入瓶颈），有的还明显不如M06系列，特别是M11系列搞得极差。2016年提出的MN15和MN15L可以稍微留意下，从Gaussian 16开始都支持了。MN15提出的目的是对过渡金属和主族可以描述得比较均衡，但是算主族比M06-2X没优势，算含过渡金属的体系有时不错有时一般（和体系关系很大），实际用处较有限。MN15L的提出意图是取代M06L的地位，但也没太实质性的进步。

wB97XD提出后，其中作者之一蔡政达又提出了wB97X-D3，把色散校正从DFT-D2改成了DFT-D3，也改了些其它参数。wB97X-D3的精度比wB97XD整体有所改善，包括热化学、弱相互作用等方面，只要程序能支持wB97X-D3（ORCA支持，但Gaussian不支持），那么前面说的wB97XD都建议改为wB97X-D3。

2016年提出的wB97M-V基本上是所有非双杂化泛函里算基态能量方面问题的平均精度最高的，比双杂化泛函里排名靠后的那些如B2PLYP精度还更高。虽然wB97M-V在参数化过程中主要考虑的是主族，但很多文章测试发现它对于含过渡金属的体系的表现也很好。但此泛函只有ORCA、Q-Chem等极少数程序支持，而最主流的Gaussian（截止到G16 C.02）尚不支持。ORCA挺好用，如果你肯用的话，那么在普通泛函的范畴里，很建议你计算能量的时候用wB97M-V，尽管它的耗时比B3LYP明显要高。不过如果你在ORCA里利用RIJCOSX技术加速的话，算较大体系时ORCA里用wB97M-V的耗时比Gaussian里计算B3LYP的耗时都低。注意ORCA从5.0开始虽然已支持了wB97M-V的解析梯度，但不支持解析Hessian，故振动分析极其昂贵。由于做优化之后通常需要做振动分析检验虚频，因此opt freq不建议用wB97M-V而建议用B3LYP-D3(BJ)，速度又快结果通常又不错，容易收敛，还不容易有虚频（需要注意B3LYP局限性，看《谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的》http://sobereva.com/557，很多情况用wB97XD会好得多）。wB97M-V至少在ORCA里有个问题就是SCF收敛难度往往较大，有时我发现改用比wB97M-V更早提出的wB97X-V更容易SCF收敛。wB97X-V各方面特征和wB97M-V相近，虽然整体精度稍逊一点（不过仍远好于更老的wB97XD），但在计算过渡金属体系上根据有的文章的测试还更好一点。PS：虽然wB97M带DFT-D3(BJ)或DFT-D4校正的版本可以算解析Hessian，但我还发现ORCA（至少是5.0）里wB97M-D3(BJ)或wB97M-D4做振动分析时候CPSCF往往很难收敛，而wB97X带色散校正的版本都没这个问题。故总的来说wB97M开头的泛函我建议只在算单点时候用，本来这类泛函做opt freq任务比起wB97XD或wB97X-D3也没什么明显优势。

Gaussian官方的Exploring那本书的第3版里从头到尾几乎都用APFD泛函，令初学者以为APFD万能，这真是超级坑爹！在诸多泛函的测试文章里，APFD无一例外表现都极其平庸，比B3LYP-D3(BJ)没有什么优势，目前文献里也极少有人用APFD，这是个绝对冷门泛函。如果你也用这个泛函，到时候审稿人问你为什么要用这个，你就傻眼了，根本找不到什么像样的文章能支持你用这个泛函。所以大家切勿被那本书误导！（那本书里错误太多了，初学者看了那本书可真是会被害惨。量子化学的学习正途见《谈谈学量子化学如何正确地入门》http://sobereva.com/355）。

5 DFT计算的几个经验和相关问题

下面我提及一些其它和方法选择有关的问题和实际计算经验方面的内容。

经常有人问算A问题和B问题用的泛函不相同能不能行。关于这点，我专门写过一篇博文，没看过者强烈建议一看：《量子化学研究中切换泛函应当注意的问题》（http://sobereva.com/415）。

如果是比较大的有机体系，特别是体系很柔的情况，虽然用算有机体系常用的M06-2X做优化的结果是挺好的，但是由于其本身泛函形式复杂而耗时比B3LYP高，再加上有时候需要用很高质量积分格点才能解决其几何优化难收敛、有细小虚频之类麻烦事，所以从实际角度，我明显更推荐用B3LYP-D3(BJ)做优化和振动分析，这样又便宜又好。而能量计算精度明显更好的M06-2X一般我建议只在计算能量的时候再用。关于什么时候很适合用B3LYP优化结构而什么时候不能用，我专门写过一篇很详细的文章进行讨论，见《谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的》（http://sobereva.com/557），非常建议一看。

一般来说，除非研究的体系很小，否则一般都是建议用普通泛函做几何优化和振动分析，而算能量的时候才用双杂化泛函。双杂化泛函做优化，尤其是振动分析非常昂贵。优化和振动分析对计算级别的敏感度远低于能量计算，见《浅谈为什么优化和振动分析不需要用大基组》（http://sobereva.com/387），因此双杂化泛函比普通泛函在优化和振动分析方面的整体优势远低于算能量方面的优势。

如果你是Gaussian用户，希望改用ORCA来以显著更低的耗时达到相似的精度，或者以相仿佛的耗时达到明显更高精度，强烈建议读者看一下此文里面对ORCA中推荐使用的计算级别的详细介绍：《详谈Multiwfn产生ORCA量子化学程序的输入文件的功能》（http://sobereva.com/490），用上ORCA后你会打开新世界的大门。ORCA的安装见《量子化学程序ORCA的安装方法》（http://sobereva.com/451）。更多与ORCA有关的文章见http://sobereva.com右侧的ORCA分类。笔者在北京科音高级量子化学培训班（http://www.keinsci.com/workshop/KAQC_content.html）里对ORCA有全面详细的讲解、给了大量例子，可以快速上手并用得游刃有余，很推荐参加。

做双杂化泛函单点计算或优化任务时，尤其强烈建议用ORCA，因为ORCA利用RI技术，算双杂化泛函能量的耗时和内存/硬盘使用量远远低于Gaussian。有个性能较好的双路服务器的话，比如笔者的Intel 36核机子，ORCA开RI时用双杂化泛函结合def2-TZVPP算六十多个原子的有机体系六分钟左右就能完成，临时文件也就占3GB，而使用Gaussian则需要三个小时左右，rwf文件最多能占1TB硬盘（但可以用maxdisk关键词事先指定硬盘空余空间来避免实际占用那么多）。这样的级别用ORCA算两百个以内原子的有机体系都能应付得了，结合def2-QZVPP的话算一百个以内原子的体系也没问题，而用Gaussian是绝对不可能算得动的。更多信息看此文的对比测试：《大体系弱相互作用计算的解决之道》（http://sobereva.com/214）。

如果你要计算的体系含有几百个非氢原子，而且体系较为接近球形，最典型的就是很大的原子团簇、从晶体结构里抠出来的簇模型等，用CP2K是效率最高的选择，比ORCA还要快得多得多，普通双路服务器上对上千原子算单点能也完全不是难事（对于纯泛函甚至能算得动几千原子）。因为CP2K利用平面波作为辅助基，算这样体系的库仑作用耗时远低于量子化学程序。CP2K完全免费、功能强大，结合Multiwfn使用很容易上手，通过北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）可以迅速掌握。

6 关于CCSD(T)

如果你想得到比双杂化泛函再高一个档次的结果，通常就是使用CCSD(T)。注意CCSD(T)不是什么情况都是非常理想的，碰到静态相关很强的情况结果并不令人满意，需要用多参考方法，诸如NEVPT2乃至更高级的MRCI。用CCSD(T)时通常是cc-pVTZ起步（起码不能低于def2-TZVP），此级别已经很昂贵了，对于没有对称性的有机体系来说，哪怕是用很好的服务器，一般最多也就能算得动不超过30个原子的体系。想算到更大的体系需要用低标度耦合簇方法，主流的是ORCA支持的DLPNO-CCSD(T)，结合tightPNO选项用cc-pVTZ的话对于>70个原子也照样算得动（比如笔者36核机子算66个原子有机体系花了8个多小时，最多占耗硬盘117GB），而且精度和原版CCSD(T)差异甚微。在北京科音高级量子化学培训班（http://www.keinsci.com/workshop/KAQC_content.html）里专门有一节“低标度耦合簇方法”对此类方法做了全面讲解和示例。

对于CCSD(T)这样很高档次方法而言，从cc-pVTZ/def2-TZVPP进一步提升到cc-pVQZ/def2-QZVPP也对结果有不可忽视的改进，但CCSD(T)结合QZ档次基组耗时更是巨高，有两个降低耗时的策略：
(1)用显式相关F12计算。用CCSD(T)-F12结合cc-pVDZ-F12基组就能达到CCSD(T)/cc-pVQZ的精度而耗时低得多得多。北京科音高级量子化学培训班（http://www.keinsci.com/workshop/KAQC_content.html）里专门有一节“显式相关计算”对这种方法做了全面深入的讲解并给出了在ORCA中使用的示例。
(2)借助MP2来降低耗时。即CCSD(T)/cc-pVDZ（或TZ）加上MP2计算的cc-pVQZ相对于cc-pVDZ（或TZ）对电子相关能的改进量。这种做法很常用，但没有上一种做法优雅、稳健，《各种后HF方法精度简单横测》（http://sobereva.com/378）的第5节给了个例子。

7 不要用HF和MP2！

我不得不提一下Hartree-Fock (HF)。居然这年头，在网上答疑时，还时常看到初学者居然使用HF来做计算！用HF在现在来看完全就是在搞笑！这些初学者多半是被一些特别过时，或者作者水平很低劣的文章或书籍给严重误导。除非泛函选用得超级烂，否则DFT的结果绝对比HF强得多！而在Gaussian里，用普通泛函做DFT只比HF贵一点点而已，显然如今没有任何理由用HF。审稿人只要稍微有一丁点最基本的理论计算常识，几乎都不可能让使用HF计算的文章发表。自从DFT兴起之后，特别是90年代初B3LYP提出并立刻流行后，HF根本就无任何用武之地了。而且利用RIJ技术，纯泛函的耗时比HF还能更低，从计算耗时层面上也完全没有用HF的理由。

还有些初学者居然还用MP2，并且自我感觉良好，在这里也猛烈批评一下。MP2计算耗时显著高于普通泛函，尤其是对于大体系、较大的基组来说，而对于需要计算能量导数的情况耗时更是贼高，还特别耗内存、硬盘。随便选个稍微靠谱的泛函，结果都不比MP2差，如果泛函选好了，结果精度更是吊打MP2，这在《各种后HF方法精度简单横测》（http://sobereva.com/378）就已经清楚地体现出来。而且，双杂化泛函的耗时和MP2是同一档次，精度更是远远甩开MP2几条街。因此，这年头用MP2完全是是费力不讨好，耗时又高精度又烂，很容易糟有经验的审稿人和同行吐槽、批评、鄙视，被认为缺乏计算化学常识。对于算开壳层体系、过渡金属体系，MP2表现得更是垃圾。MP2倒也不是一点使用价值也没有，但少得可怜，也就这些：(1)计算氢键键能挺好（但必须结合足够大的基组如jul-cc-pVTZ，而且比起revDSD-PBEP86-D3(BJ)等很好的双杂化泛函并没额外优势） (2)计算NMR、超极化率比现有的各种普通泛函通常都要好（但双杂化泛函里有能超过MP2的） (3)可以借助MP2的数据以比较便宜的方式估算CCSD(T)/大基组的能量（此做法在前面已提到了）。

作者
Author: gengle 时间: 2014-12-30 12:01
这总结太好了！

作者
Author: alexwpch 时间: 2014-12-30 13:54
拜读！
要是再有基组的对比更好了，像我们新手，除了泛函把握不准，基组也会把握不准。
例如，6-311+g(d)对应月份基组哪个级别，对应 def2-XXX哪个水平，有时还有些相对冷门的；做配合物时，过渡金属的LanL2DZ又与哪个级别差不多，能一起用。
搞张类似于天梯的那种图，在需要对比不同的方法时，能很方便的找出差不多水平的基组。
另外，我是不是要的太多了。。。。

作者
Author: sobereva 时间: 2014-12-30 14:10

alexwpch 发表于 2014-12-30 13:54
拜读！
要是再有基组的对比更好了，像我们新手，除了泛函把握不准，基组也会把握不准。
例如，6-311+g(d) ...

基组以后会专门再写帖子的。
基组的选择很难简单地轨道天梯那种图，不同基组有不同特点，适合不同计算任务、计算方法、体系。
不同基组在构建思想上差异较大，不好比。比如6-311+G*就没有类似的月份基组，因为月份基组是基于dunning相关一致性基组，这种基组专用在后HF计算，极化函数很多，到了3-zeta时极化函数就很多了。

ahlrichs基组有大致这种关系：
def-SV(P)/SVP≈6-31G*/6-31G**。用于DFT定性计算
def-TZVP>=6-311G(d,p)。用于DFT精确计算。
def-TZVPP≈6-311G(2df,2pd)。用于DFT精确计算，MP2较高精度计算。
def2-SV(P)/SVP≈6-31G*/6-31G**
def2-TZVP>=6-311G(2df,p)
def2-TZVPP>=6-311G(2df,2pd)，特别是有些元素多一层更高角动量的极化函数

值得一提的是，ahlrichs基组质量比pople的好很多，可能表面看起来都是同一档次的，但是实际结果会更好（但由于收缩度高，而且没有利用sp壳层，也更耗时）

作者
Author: 小范范1989 时间: 2015-1-1 19:25
敢问sob，为啥这么牛，懂得这么多.......膜拜啊，大神就是大神。
元旦快乐。

作者
Author: fcfsdu 时间: 2015-1-3 02:41
干了好几年的计算了，自己最近才对各种泛函有了点浅显的认识。目前主要用的还是wB97X-D（主要做几十个原子的中性有机分子的光化学性质），觉得这个最可靠，激发态的处理尤其不错；基态的结果也很好，不过有时候还需要和B3LYP来对照一下。要是能早点看到Sob原来的帖子就能省很多力气了。

作者
Author: brothers 时间: 2015-1-24 16:47
受教~~

作者
Author: nkallwar 时间: 2015-3-11 23:11
计算双核金属用PBE、BP86，勿用杂化（see JCTC,8,908）

对于这一句，想请教一下，双核金属是指在有M-M金属键的情况下勿用杂化，还是 M-OO-M这种桥联的双核金属也不要用杂化呢？因为后面又看到一句杂化泛函对双核金属键描述不好....

另外一个，M06-2x本身处理弱相互作用比较好，（加上DFT-D3校正更好）这句指的是在M06-2x基础上再加上EmpiricalDispersion=GD3关键词来进行D3校正会有更好的效果吗？

还有这一句“另外M06-2X是针对主族进行优化的，对于过渡金属，改用M06会更合适”，那对于诸如金属配合物这种大量主族元素，只有一两个过渡金属原子的体系，哪个更合适一些？因为金属配合物更多的是关注金属和配位原子之间的作用，是不是M06更合适些？

作者
Author: sobereva 时间: 2015-3-12 02:23

nkallwar 发表于 2015-3-11 23:11
计算双核金属用PBE、BP86，勿用杂化（see JCTC,8,908）

对于这一句，想请教一下，双核金属是指在有M-M金 ...

指的是M-M这种情况

是的

一般用B3LYP就行了，比较快，研究配合物大多也比较可靠。如果要考虑弱相互作用，可以再加上D3校正。用M06也可以。

作者
Author: 止水 时间: 2016-3-29 20:12
赞，，真好，谢谢

作者
Author: 498746012 时间: 2016-4-3 23:24
谢谢老师

作者
Author: 10yaoxue 时间: 2016-4-4 10:33
不错。

作者
Author: Frank 时间: 2017-2-15 09:03
请问社长，计算局域价层激发时，PBE0表现要比B3LYP要好，应该怎么来解释呢？是因为两种泛函的HF%不同吗？

作者
Author: yangyinuo823 时间: 2017-2-15 14:54
请教一下，如果用m06的话，可以加D3矫正吗？

作者
Author: sobereva 时间: 2017-2-15 18:01

yangyinuo823 发表于 2017-2-15 14:54
请教一下，如果用m06的话，可以加D3矫正吗？

可

作者
Author: sobereva 时间: 2017-2-15 18:02

Frank 发表于 2017-2-15 09:03
请问社长，计算局域价层激发时，PBE0表现要比B3LYP要好，应该怎么来解释呢？是因为两种泛函的HF%不同吗？

只是好一点而已，整体差不多。
HF成份是主要原因

作者
Author: sobereva 时间: 2018-4-30 13:25
今日对此文做了较大update

作者
Author: 研道 时间: 2018-5-2 15:58
写的真好

作者
Author: SongChao 时间: 2018-5-11 13:54

有机分子结构优化：PBE0（弱相互作用明显时应加DFT-D3校正）、M06-2X
计算分子振动（谐振、非谐振都算）：B3LYP
计算有机体系热力学数据（反应能、构型能量差、能垒等）：M06-2X（加上DFT-D3校正更好）

Gaussian计算中分子的振动分析和热力学数据在使用Freq关键词后会一同给出，请问为什么B3lYP和M06-2X一个对频率更加准确一个对热力学数据更加准确？
如果我用M06-2X计算得到接结构来用B3LYP计算分子振动数据，是不是要比同时用B3LYP计算结构和能量更加准确，而参考热力学数据时使用M06-2X计算的结果？

作者
Author: sobereva 时间: 2018-5-12 04:13

SongChao 发表于 2018-5-11 13:54
Gaussian计算中分子的振动分析和热力学数据在使用Freq关键词后会一同给出，请问为什么B3lYP和M06-2X一 ...

热力学数据来自于电子能量以及热力学校正量
B3LYP和M06-2X算热力学校正量这一部分的时候半斤八两，但由于M06-2X算的电子能量比B3LYP准，所以算热力学数据好。
当然你也可以用B3LYP的热力学校正量加到M06-2X的电子能量上，结果也一样好

作者
Author: GoldenBaby 时间: 2018-10-16 11:44
请问sob，如果想用纯泛函计算有机配体的配合物，用什么泛函比较好呢？

作者
Author: sobereva 时间: 2018-10-16 13:21

GoldenBaby 发表于 2018-10-16 11:44
请问sob，如果想用纯泛函计算有机配体的配合物，用什么泛函比较好呢？

仔细看帖子，写明了

作者
Author: sobereva 时间: 2018-12-29 12:12
鉴于ORCA 4.1已经支持了wB97M-V，对帖子相应部分内容做了些修改

作者
Author: sobereva 时间: 2019-4-30 15:55
今日对本文做了极大的更新

作者
Author: gauss98 时间: 2019-6-4 13:53
感谢老师介绍！

您对def2基组比较推崇，但从数量上看，def2svp碳原子16个基， 6-31G*15个，可以说是同一档次，但 6-311G**才18个基，def2tzvp高达31个，计算量大了几倍啊，如果想在def2svp上增加一点，tzvp又算不动，有没有折中的方案？

另外，我原来一直用b3lyp-d3,现在看似乎PBE0也受您推荐，有没有比较优劣的介绍文章？我的体系时过渡金属催化剂比如 Ti， Rh催化有机反应

谢谢！

作者
Author: sobereva 时间: 2019-6-5 10:35

gauss98 发表于 2019-6-4 13:53
感谢老师介绍！

您对def2基组比较推崇，但从数量上看，def2svp碳原子16个基， 6-31G*15个，可以说是同 ...

def2-TZVP算不动就def-TZVP，只比6-311G**大一丁点
def 1代里就只有def-TZVP还有存在意义，即弥补def2-SVP和def2-TZVP的空白
还有种做法是用def2-TZVP(-f)，即去掉f极化，这对DFT计算影响不算特别大，但耗时低得多得多
有机体系算弱相互作用B3LYP-D3(BJ)完全没问题，也比PBE0-D3(BJ)好。但近些年有不少可靠的文章指出对于算过渡金属反应相关问题，PBE0-D3(BJ)更佳，因为本身PBE0算这类问题就比B3LYP更强点。

作者
Author: gauss98 时间: 2019-6-5 14:06

sobereva 发表于 2019-6-5 10:35
def2-TZVP算不动就def-TZVP，只比6-311G**大一丁点
def 1代里就只有def-TZVP还有存在意义，即弥补def2-S ...

谢谢！
不过看到您以前的帖子说 deftzvp对 Mg 有缺陷？我的体系恰好有Mg，另外，def2tzvp（-f）在gaussian里边是怎么输入的？

作者
Author: sobereva 时间: 2019-6-6 20:28

gauss98 发表于 2019-6-5 14:06
谢谢！
不过看到您以前的帖子说 deftzvp对 Mg 有缺陷？我的体系恰好有Mg，另外，def2tzvp（-f）在gaus ...

def-对碱土金属有问题
自行去BSE上拷定义，把里面f给删了即可

作者
Author: gray 时间: 2019-10-25 13:57
感谢大师

作者
Author: 乐嘻嘻嘻 时间: 2020-3-18 11:18
老师好，最近我师妹做了个体系，优化有机大分子，在考虑了B3LYP,PBE0,CAM-B3LYP泛函之后，选择了B3PW91/6-31G(d) 方法进行了结构优化，总觉得不太合适，又没有比较正式的理由~麻烦问一下老师，能否选用这种方法进行优化~谢谢~

作者
Author: snljty 时间: 2020-3-18 11:30

乐嘻嘻嘻发表于 2020-3-18 11:18
老师好，最近我师妹做了个体系，优化有机大分子，在考虑了B3LYP,PBE0,CAM-B3LYP泛函之后，选择了B3PW91/6-3 ...

建议先花十分钟把这篇博文读一遍......

作者
Author: sobereva 时间: 2020-3-19 21:19

乐嘻嘻嘻发表于 2020-3-18 11:18
老师好，最近我师妹做了个体系，优化有机大分子，在考虑了B3LYP,PBE0,CAM-B3LYP泛函之后，选择了B3PW91/6-3 ...

不要用B3PW91，非主流，也没有什么理由用这个（无论是精度还是流行程度上）
有机体系就用B3LYP-D3(BJ)优化就很适合，非常流行，基本不会被找碴（除了个例，诸如18轮烯、18碳环之类涉及长程大共轭体系，此类情况B3LYP表现不如wB97XD等合理，参考诸如《一篇最全面、系统的研究新颖独特的18碳环的理论文章》http://sobereva.com/524）

作者
Author: Q-Chembio-llg 时间: 2020-3-25 22:14

小范范1989 发表于 2015-1-1 19:25
敢问sob，为啥这么牛，懂得这么多.......膜拜啊，大神就是大神。
元旦快乐。

Sob老师好像是北科的博士，666

作者
Author: CCCC 时间: 2020-7-11 18:43
请问像Fe4S4这种团簇，也是用纯泛函比较好吗？之所以选择纯泛函是因为描述键长不合理吗？

作者
Author: sobereva 时间: 2020-7-12 09:04

CCCC 发表于 2020-7-11 18:43
请问像Fe4S4这种团簇，也是用纯泛函比较好吗？之所以选择纯泛函是因为描述键长不合理吗？

“纯泛函比较好吗？之所以选择纯泛函是因为描述键长不合理吗？”你的问题明显前后矛盾

用纯泛函，或者HF成分低的，诸如TPSSh

作者
Author: CCCC 时间: 2020-7-12 16:43

sobereva 发表于 2020-7-12 09:04
“纯泛函比较好吗？之所以选择纯泛函是因为描述键长不合理吗？”你的问题明显前后矛盾

用纯泛函，或者 ...

嗯嗯，我后面这句话掉了主语，谢谢回复~

作者
Author: 乐嘻嘻嘻 时间: 2020-7-23 20:57
老师好，我模拟了有机大分子的吸收光谱，如图所示，之后想考虑分子内的电荷转移等特征，我的问题是，从图来看，B3LYP泛函和实验值吻合的稍微好一点，但看这篇博文，B3LYP不太利于描述电荷转移特征。但是如果要用wB97XD泛函，光谱来看和实验值差的较大，想问一下这种情况应该如何处理？麻烦老师。

作者
Author: sobereva 时间: 2020-7-24 07:45

乐嘻嘻嘻发表于 2020-7-23 20:57
老师好，我模拟了有机大分子的吸收光谱，如图所示，之后想考虑分子内的电荷转移等特征，我的问题是，从图来 ...

能否与实验较好吻合是由多方面因素共同影响的，泛函仅仅是其中一个方面
在当前情况下，如果你只是单纯想通过选择泛函来让模拟结果与实验较好吻合而不管其它因素带来的问题，就用B3LYP就完了。

作者
Author: 乐嘻嘻嘻 时间: 2020-7-24 09:40

sobereva 发表于 2020-7-24 07:45
能否与实验较好吻合是由多方面因素共同影响的，泛函仅仅是其中一个方面
在当前情况下，如果你只是单纯想 ...

谢谢老师，主要是之后想考察激发特征，怕B3LYP会引起假的激发态。

作者
Author: sobereva 时间: 2020-12-26 17:21
今日对本文进行了大幅改写和扩充

作者
Author: snljty 时间: 2020-12-26 18:02

sobereva 发表于 2020-12-26 17:21
今日对本文进行了大幅改写和扩充

卢老师是不是可以考虑抽时间再开一次“最常用的DFT泛函”和“最常用的量子化学软件”投票了，感觉这几年会有不少变化了。另外“最常用的分子动力学软件”等等是不是也可以考虑开一个

作者
Author: 喵星大佬 时间: 2020-12-27 00:58

ts58 发表于 2020-12-26 21:39
感谢站长，希望卢老师也能顺带更新周期性体系计算中的泛函选择（虽然大部分人都是PBE/HSE/SCAN或者+U）。在 ...

VASP之类的在周期性里用杂化泛函计算量会暴增而几乎算不动，内存消耗也很大，这事在cp2k里通过ADMM可以一定程度上缓解。大部分情况优化结构，晶格常数之类的都是纯泛函，只有在要精确计算能带结构之类的时候会去用杂化泛函。纯泛函的话几个常用的区别不算很大，PBE在大多数情况下比较合理，+U这事情就另说了。

作者
Author: sobereva 时间: 2020-12-27 06:34

snljty 发表于 2020-12-26 18:02
卢老师是不是可以考虑抽时间再开一次“最常用的DFT泛函”和“最常用的量子化学软件”投票了，感觉这几年 ...

打算21年弄
每三年一次

作者
Author: sobereva 时间: 2021-1-12 05:00
对本文进行了更新，wB2GP-PLYP不再是首要推荐

关于用双杂化泛函基于TDDFT算激发态：JCTC (2021) DOI: 10.1021/acs.jctc.0c01135的测试表明一般激发态用DSD-PBEP86计算最理想，其次是此文提出的范围分离双杂化泛函RS-PBEP86，另一个范围分离双杂化泛函wB2GP-PLYP表现差一些，但仍比用普通泛函强。CT特征显著的激发态则应当用RS-PBEP86或wB2GP-PLYP。各种类型激发态都涉及的话最好用RS-PBEP86。ORCA支持wB2GP-PLYP和DSD-PBEP86算TDDFT。

作者
Author: 佛祖家的娃儿 时间: 2021-4-1 20:19
文中说“首选wB97X-V或wB97M-V，其次PBE0-D3(BJ)和MN15”，那请问在仅含有一个Rh原子的配合物中（总共70-100个原子），PBE0-D3(BJ)在优化和计算中跟wB97X-V有哪些差距呢？
另外博文中提到过当弱相互作用明显时B3LYP-D3(BJ)也是一个不错的选择，该体系是不是也可以选用B3LYP-D3(BJ)呢？

作者
Author: wzkchem5 时间: 2021-4-1 20:30

佛祖家的娃儿发表于 2021-4-1 20:19
文中说“首选wB97X-V或wB97M-V，其次PBE0-D3(BJ)和MN15”，那请问在仅含有一个Rh原子的配合物中（总共70-10 ...

wB97X-V、wB97M-V准但是慢，而且支持它的梯度的程序少。不过这种过渡金属体系可能HF成分太高了不好，统一用PBE0-D3我觉得就行了。但这个也看看其他人的意见。
有过渡金属的体系，一般认为PBE0-D3总比B3LYP-D3好

作者
Author: 佛祖家的娃儿 时间: 2021-4-2 10:53

wzkchem5 发表于 2021-4-1 20:30
wB97X-V、wB97M-V准但是慢，而且支持它的梯度的程序少。不过这种过渡金属体系可能HF成分太高了不好，统一 ...

体系当中存在弱相互作用，PBE0是否还适用呢？
还是说应该换用别的方法，比如说文献中提到的M06-2X之类的？

作者
Author: 喵星大佬 时间: 2021-4-2 11:20

佛祖家的娃儿发表于 2021-4-2 10:53
体系当中存在弱相互作用，PBE0是否还适用呢？
还是说应该换用别的方法，比如说文献中提到的M06-2X之类的 ...

加了D3(BJ)显然问题不大

作者
Author: 佛祖家的娃儿 时间: 2021-4-2 11:23

喵星大佬发表于 2021-4-2 11:20
加了D3(BJ)显然问题不大

谢谢~

作者
Author: snljty 时间: 2021-5-6 15:48
各位老师，求推荐一些说明PWPB95-D3(BJ)稳健性和精度都很好的Benchmark文献，谢谢！

作者
Author: wzkchem5 时间: 2021-5-6 22:11

snljty 发表于 2021-5-6 08:48
各位老师，求推荐一些说明PWPB95-D3(BJ)稳健性和精度都很好的Benchmark文献，谢谢！

搜一下Grimme做的benchmark吧

作者
Author: 喵星大佬 时间: 2021-5-7 01:08

snljty 发表于 2021-5-6 15:48
各位老师，求推荐一些说明PWPB95-D3(BJ)稳健性和精度都很好的Benchmark文献，谢谢！

Phys. Chem. Chem. Phys. 2018, 20, 23175

作者
Author: sobereva 时间: 2021-5-7 03:14

snljty 发表于 2021-5-6 15:48
各位老师，求推荐一些说明PWPB95-D3(BJ)稳健性和精度都很好的Benchmark文献，谢谢！

J. Chem. Theory Comput. 2018, 14, 5, 2596-2608

作者
Author: 北大-陶豫 时间: 2021-6-15 14:09
请问，涉及过渡金属、镧系锕系配合物有关的情况，有没有一些推荐的综述文献呢？
我最近看了一篇超人气的综述：Christopher J. Cramer and Donald G. Truhlar. Density functional theory for transition metals and transition metal chemistry. Phys. Chem. Chem. Phys., 2009, 11, 10757–10816 但是里面感觉都是一些评测文献的杂乱罗列（也可能我没get到它的点）？有没有综述文献能够告诉我们对于过渡金属的各种问题（能量和势垒计算，结构优化，等）分别用什么泛函比较好呢？

作者
Author: sobereva 时间: 2021-6-15 18:17

北大-陶豫发表于 2021-6-15 14:09
请问，涉及过渡金属、镧系锕系配合物有关的情况，有没有一些推荐的综述文献呢？
我最近看了一篇超人气的综 ...

专门的综述文章没有，benchmark文章倒是很多

你提到的那篇文章是流水账，记得以前还有知名学者吐槽过此文

作者
Author: 北大-陶豫 时间: 2021-6-15 19:13

sobereva 发表于 2021-6-15 18:17
专门的综述文章没有，benchmark文章倒是很多

你提到的那篇文章是流水账，记得以前还有知名学者吐槽过 ...

那看来不是我文献搜索能力的问题，谢谢sob大大

作者
Author: 15029457177 时间: 2021-6-29 15:43
【小白受教

作者
Author: sobereva 时间: 2021-7-2 06:08
结合ORCA 5.0的改变，更新了本文的一些说法

作者
Author: 杜坤 时间: 2021-7-19 17:23
总结的很好，非常感谢

作者
Author: sobereva 时间: 2021-9-25 21:19
今日新加了一部分

• 闭壳层过渡金属-有机反应能和势垒（涉及到配位键的形成/断裂、多参考特征不很强的情况）：首选TPSS0-D4，其次PBE0-D4（PBE0-D3(BJ)稍逊一丝）。纯泛函当中B97M-V表现最好，但用r2SCAN-3c组合方法的话不仅便宜且精度还强于B97M-V/def2-QZVPP

作者
Author: 喵星大佬 时间: 2021-9-28 05:03

sobereva 发表于 2021-9-25 21:19
今日新加了一部分

应该是开壳层

作者
Author: sobereva 时间: 2021-9-28 08:30

喵星大佬发表于 2021-9-28 05:03
应该是开壳层

笔误改了

作者
Author: sobereva 时间: 2022-3-6 17:53
今日更新了以下内容

关于用双杂化泛函算激发态：局域激发和分子内电荷转移激发用DSD-PBEP86做TDDFT计算最理想，ORCA直接支持。分子间电荷转移激发用RS-PBEP86双杂化泛函基于SOS-CIS(D)形式计算最佳，这对于局域激发和分子内电荷转移激发描述得也较好，但目前只有MRCC程序支持。ORCA里算分子间电荷转移激发建议基于RSX-QIDH双杂化泛函做TDDFT，虽然略逊于RS-PBEP86/SOS-CIS(D)但远强于DSD-PBEP86做TDDFT。这些结论来自JCTC (2022) DOI: 10.1021/acs.jctc.1c01307的测试。

作者
Author: danni511 时间: 2022-4-28 10:00
sob老师，您好，我计算的体系既涉及过渡金属与过渡金属成键，又涉及到过渡金属与惰性气体或主族元素原子的相互作用，请问我用PBE0，TPSSh还是哪个纯泛函比较好呢，谢谢老师

作者
Author: sobereva 时间: 2022-4-29 08:15

danni511 发表于 2022-4-28 10:00
sob老师，您好，我计算的体系既涉及过渡金属与过渡金属成键，又涉及到过渡金属与惰性气体或主族元素原子的 ...

ORCA支持的r2SCAN-3c比较理想

作者
Author: danni511 时间: 2022-4-29 15:14

sobereva 发表于 2022-4-29 08:15
ORCA支持的r2SCAN-3c比较理想

老师您好，我用的是g09，可以用这个泛函吗

作者
Author: wzkchem5 时间: 2022-4-29 15:26

danni511 发表于 2022-4-29 08:14
老师您好，我用的是g09，可以用这个泛函吗

不能，必须用orca。
反正orca对学术界免费，而且只跑r2SCAN-3c的话输入文件也好写，没啥门槛

作者
Author: sobereva 时间: 2022-4-30 07:33

danni511 发表于 2022-4-29 15:14
老师您好，我用的是g09，可以用这个泛函吗

如果非要用G09的话，看你研究的具体体系是什么，研究的问题重点是什么，要算什么任务
r2SCAN-3c对于过渡金属之间成键，以及弱相互作用，描述得都比较不错。用G09的话，就得考虑要以什么为优先考虑。如果不具体区分的话，TPSS-D3(BJ)可以优先考虑。

作者
Author: danni511 时间: 2022-4-30 08:45

wzkchem5 发表于 2022-4-29 15:26
不能，必须用orca。
反正orca对学术界免费，而且只跑r2SCAN-3c的话输入文件也好写，没啥门槛

谢谢老师的回答，我没用过orca呢，以后也学习一下。

作者
Author: danni511 时间: 2022-4-30 09:13

sobereva 发表于 2022-4-30 07:33
如果非要用G09的话，看你研究的具体体系是什么，研究的问题重点是什么，要算什么任务
r2SCAN-3c对于过渡 ...

老师，我想用g09研究过渡金属团簇和CO2、惰性气体元素等主族元素构成的体系，主要用泛函优化基态结构，振动分析，轨道分析，单点能计算想用CCSD（T）方法，就不用泛函了，谢谢老师

作者
Author: wzkchem5 时间: 2022-4-30 15:08

danni511 发表于 2022-4-30 02:13
老师，我想用g09研究过渡金属团簇和CO2、惰性气体元素等主族元素构成的体系，主要用泛函优化基态结构，振 ...

稍微大一点的团簇用高斯跑CCSD(T)是跑不动的，必须用orca做DLPNO-CCSD(T)计算

作者
Author: danni511 时间: 2022-4-30 17:24

wzkchem5 发表于 2022-4-30 15:08
稍微大一点的团簇用高斯跑CCSD(T)是跑不动的，必须用orca做DLPNO-CCSD(T)计算

跑了一个确实很慢

作者
Author: tiandikuoyuan 时间: 2022-5-20 15:51
老师，您好，我看了J. Chem. Theory Comput. 2020, 16, 7, 4416–4428关于配合物不同多重度能量差计算的泛函对比，发现M06-2X在计算3d金属配合物能量差时同样比较优秀，那这个泛函能不能使用呢，因为您也说M06-2X是针对主族元素进行参数优化的。我用MN15L计算配合物能量时，有一个多重度的计算试了一些方法SCF始终不收敛，所以想换一个泛函。

作者
Author: sobereva 时间: 2022-5-21 04:22

tiandikuoyuan 发表于 2022-5-20 15:51
老师，您好，我看了J. Chem. Theory Comput. 2020, 16, 7, 4416–4428关于配合物不同多重度能量差计算的泛 ...

不具备普遍性，被批的概率极大
如果你用的Gaussian，可以尝试ORCA 5来算，默认会用的trah对于一些难收敛的过渡金属配合物可能最终能收敛。到时候如果想和之前Gaussian算的其它的对比，可以把ORCA收敛的波函数按下文的做法“gbw到chk的批量转换脚本” 作为Gaussian的初猜，应当能顺利收敛
将Gaussian等程序收敛的波函数作为ORCA的初猜波函数的方法
http://sobereva.com/517（http://bbs.keinsci.com/thread-14905-1-1.html）

作者
Author: tiandikuoyuan 时间: 2022-5-21 09:08

sobereva 发表于 2022-5-21 04:22
不具备普遍性，被批的概率极大
如果你用的Gaussian，可以尝试ORCA 5来算，默认会用的trah对于一些难收敛 ...

好的，非常感谢

作者
Author: mfdsrax2 时间: 2022-10-22 18:59
这些计算任务都太专业了，如果只是最基础的有机小分子的结构优化和频率计算，是不是用什么泛函都可以？

作者
Author: sobereva 时间: 2022-10-23 04:09

mfdsrax2 发表于 2022-10-22 18:59
这些计算任务都太专业了，如果只是最基础的有机小分子的结构优化和频率计算，是不是用什么泛函都可以？

当然不是
诸如用LDA泛函计算直接被审稿人打回去，无论什么分子

作者
Author: chenf 时间: 2023-11-30 20:59
很棒

作者
Author: JCenter 时间: 2024-11-15 12:02
老师，您这个帖子中，有写到，在计算电子密度分布时，普通泛函推荐PBE0、TPSSh。但我另外看了老师您写的有关于电子密度差、原子电荷计算等与电子密度、电子转移相关的博文中，基态情况下，计算波函数信息所用的级别基本都用的B3LYP泛函，级别差不多在B3LYP/6-311g**,或者B3LYP/def2tzvp，外加iefpcm隐式溶剂模型表现溶剂对溶质电子的极化。这方面，B3Lyp泛函与PBE0比如何呢。

作者
Author: sobereva 时间: 2024-11-15 12:14

JCenter 发表于 2024-11-15 12:02
老师，您这个帖子中，有写到，在计算电子密度分布时，普通泛函推荐PBE0、TPSSh。但我另外看了老师您写的有 ...

没有太大差别

作者
Author: JCenter 时间: 2024-11-15 13:15

sobereva 发表于 2024-11-15 12:14
没有太大差别

好嘞，感谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)