计算化学公社

标题: 计算三甲基硅基乙炔分子，虚频消不掉（倒也不是消不掉的问题） [打印本页]

作者
Author: 169ccjc 时间: 2022-9-17 19:53
标题: 计算三甲基硅基乙炔分子，虚频消不掉（倒也不是消不掉的问题）
用的是PBE0-D3泛函，def2svp基组，scrf=iefpcm，高斯版本是g09 d.01。
主要任务是计算一个涉及含过渡金属配合物的反应能垒，里面有有一个双分子反应（过渡金属配合物+三甲基硅基乙炔），对三甲基硅基乙炔的孤立分子做opt+freq的时候出现一个-12.88的虚频，直接加opt=tight或verytight以及maxstep=n的方法都消不掉，最后直接加上opt=calcall算，发现结构一直震荡无法收敛。
用上opt=tight和int=ultrafine，虚频就能消掉了，但是因为看到PBE0对积分格点不敏感，之前一直用的是int=fine，而且一直没算孤立分子，算反应复合物的时候一直都挺正常的，算了挺多的了，这个用int=ultrafine的话之前的需要全部加上int=ultrafine重算么？有啥不重算的解决方法么？

顺便问一下，对于4d、5d的过渡金属配合物，计算单点能双杂化泛函B2PLYP-D3(BJ)比PBE0-D3强么？或者只考虑用高斯的前提下，对于这类体系还有啥好的双杂化泛函么？

作者
Author: sobereva 时间: 2022-9-18 04:28
对于DFT消虚频，用opt=tight的时候一定要用ultrafine或者更高质量格点。积分格点质量不够高往往难以达到tight收敛限而震荡，亦有可能由于积分格点质量太次导致出现虚频。
verytight完全没必要尝试
对于对比能量，积分格点必须统一。对于opt freq，可以不统一

算双杂化，除非原子数不多，否则不建议用Gaussian。用ORCA算双杂化快得多，而且有明显更适合过渡金属配合物的PWPB95-D3(BJ)或D4。对于金属有机反应，B2PLYP-D3(BJ)未必比PBE0-D3(BJ)强，在Gaussian里还多花大量时间。

作者
Author: mt13 时间: 2022-9-18 08:12
gaussian双杂化可以用revDSD-PBEP86-D3(BJ)，见sob老师博文http://sobereva.com/272

作者
Author: 169ccjc 时间: 2022-9-18 09:03

sobereva 发表于 2022-9-18 04:28
对于DFT消虚频，用opt=tight的时候一定要用ultrafine或者更高质量格点。积分格点质量不够高往往难以达到tig ...

好的，十分感谢，那我再去看看ORCA怎么用

作者
Author: ionexchangeC 时间: 2022-9-18 21:27
你用int=fine的格点在int=ultrafine格点优化完毕之后的结构下算个单点就行

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