计算化学公社

标题: 用Columbus的MRCI算这个体系，算得动吗？ [打印本页]

作者
Author: Freeman 时间: 2026-4-23 18:43
标题: 用Columbus的MRCI算这个体系，算得动吗？
本帖最后由 Freeman 于 2026-4-23 18:59 编辑

大家好，我需要算一些分子（信息如下）MRCI级别的能量、SOC和NAC。请问Columbus算得动吗？

原子数：40-60
活性空间：有些是（10o，10e），有些是（10o，16e）
基组：最好用3-zeta，实在不行也可以降级到2-zeta
激发态数：4个左右
核数：管够（不过听说这些高级波函数方法的并行效率都一般般）
时间：2周

作者
Author: kkkjjjlll 时间: 2026-5-26 08:38
不是说核数管够就行，你的显存也很重要，你需要自己先看下总的CSFs个数有多少，原子数在四五十的话，应该需要很大很大显存了，综上：你需要先看下你的机器情况

作者
Author: kkkjjjlll 时间: 2026-5-26 08:39
不是说核数管够就行，你的显存也很重要，你需要自己先看下总的CSFs个数有多少，原子数在四五十的话，应该需要很大很大显存了，综上：你需要先看下你的机器情况

作者
Author: biogon 时间: 2026-5-26 11:48
不可能，这么大的体系MRCI基本没有程序能算的动，能算动NEVPT3就不错了

作者
Author: beefly 时间: 2026-5-26 19:31
本帖最后由 beefly 于 2026-5-26 19:47 编辑

Columbus的mrci没做任何近似，大约只能算到CAS(8o,8e)的活性空间。所以CAS(10o,16e)原则上可以，而CAS(10o,10e)应该算不动。Columbus作者过去算过一些大活性空间的体系，要么组态做了限制，要么没有考虑双激发（MRCIS）。

CAS(10o,10e)活性空间的MRCI可以用内收缩MRCI（Molpro、Xian-CI），以及更快的SDS-MRCI（Xian-CI；单根大约比Molpro快5-10倍，多根更快）。

但是以上所有MRCI方法仅能用传统的积分，上限大约是300-350个基函数，40-60原子还是算不动。应该考虑用密度拟合积分，不过目前未见有主流程序实现了密度拟合MRCI的报道。仅有NEVPT2和CASPT2实现了密度拟合（Mopro，OpenMolcas，Orca），其中OpenMolcas的Cholesky分解+密度拟合CASPT2要比Molpro同类方法快很多。我用它算过40原子有机分子的CASPT2(10o,10e)，3个根几小时就算完了。

但是CASPT2有入侵态的问题，最后的参考态权重可能很小，最终的结果不可靠。反复试活性空间、调能级移动要浪费不少时间。NEVPT2原则上没有入侵态问题，即便活性空间不合理，也不会像CASPT2那么离谱。

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