明确了计算需求，想知道用什么量化程序来做比较好，求推荐！

13799725819

我想使用CCSD(T)，B2PLYP等级的方法做IRC-势能面，程序要能给出反应坐标 vs. 电子能的数据；还要关注体系的多参考态特征，比如直接做“后-CASSCF”等级的计算，有什么程序能满足需求吗？如果有免费的程序就更好了。谢谢！

wal

看看ORCA

YueliLan

ORCA

dlldwtz

OOOOOORCA

1264622935

OOOOOORCA

sobereva

小体系用Gaussian跑CCSD(T)、B2PLYP的能量完全没问题，ORCA里的标准的CCSD(T)也并不比Gaussian明显更快。而中等大小的体系，Gaussian跑CCSD(T)根本跑不动（ORCA照样跑不动，但可以用ORCA支持的其近似版DLPNO-CCSD(T)）、Gaussian跑B2PLYP也非常耗时而且很容易内存不足和硬盘写满（ORCA可以开RI近似轻松算得动）。后-CASSCF学名叫多参考方法，Gaussian虽然有OVB-MP2但不是主流，ORCA支持主流的NEVPT2和CASPT2以及更高精度的MRCI（还支持很多其它的，但不算流行）。

顺带一提，最古老的双杂化泛函B2PLYP已经完全过时了，如今没有任何理由用。双杂化范畴里，Gaussian里可以用精度拔尖的revDSD-PBEP86-D3(BJ)，参考《Gaussian中非内置的理论方法和泛函的用法》（http://sobereva.com/344），ORCA里除了能用这个外，还可以用轻微更好的顶级的双杂化泛函wB97M(2)。

在北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）里非常全面系统讲了ORCA的使用，包括上面提到的所有方法，可以容易上手。ORCA对学术界免费。

13799725819

sobereva 发表于 2025-2-25 11:21
小体系用Gaussian跑CCSD(T)、B2PLYP的能量完全没问题，ORCA里的标准的CCSD(T)也并不比Gaussian明显更快。而 ...

谢谢sob老师的讲解！我还有一些疑问：
ORCA 5.0的IRC计算只能给出笛卡尔坐标vs.电子能的数据，那么这个程序算的到底是MEP，还是和Gaussian算的IRC方法一样，只是最后给的是笛卡尔坐标呢？
如果是前者，我想了解一下MEP和IRC路径的区别大吗？
如果是后者，我是不是可以手动计算出ORCA输出路径的质重坐标vs.电子能的数据呢？（方法：计算ORCA输出路径上两帧结构的质重坐标变化，以此类推）

13799725819

wal 发表于 2025-2-25 08:41
看看ORCA

谢谢！

zjxitcc

建议使用MOKIT进行自动的多组态/多参考方法 单点计算，程序会自动构造活性轨道，自动确定活性空间，自动调用相应的量化软件衔接计算。也可以看引用MOKIT的已发表文章，看是否有与自己类似需求/使用场景的https://jeanwsr.gitlab.io/mokit-doc-mdbook/citing.html

sobereva

13799725819 发表于 2025-2-25 13:05
谢谢sob老师的讲解！我还有一些疑问：
ORCA 5.0的IRC计算只能给出笛卡尔坐标vs.电子能的数据，那么这个 ...

ORCA的目前的IRC功能做的较差，不是质量权重坐标做的（所以严格来说算MEP），而且无法给出反应坐标。可以用Gaussian调用ORCA计算：
将Gaussian与ORCA联用搜索过渡态、产生IRC、做振动分析
http://sobereva.com/422（http://bbs.keinsci.com/thread-10141-1-1.html）

13799725819

zjxitcc 发表于 2025-2-25 13:15
建议使用MOKIT进行自动的多组态/多参考方法 单点计算，程序会自动构造活性轨道，自动确定活性空间，自动调 ...

谢谢，我关注一下

13799725819

sobereva 发表于 2025-2-26 12:47
ORCA的目前的IRC功能做的较差，不是质量权重坐标做的（所以严格来说算MEP），而且无法给出反应坐标。可以 ...

明白了，多谢sob老师！

DavyHou

既然跑CCSD(T),为啥不用molpro

tonganlhy

DavyHou 发表于 2025-3-18 13:59
既然跑CCSD(T),为啥不用molpro

不是免费。跑irc也麻烦