氧气激发能的计算与实验差距较大

zjxitcc · 发表于 Post on 2021-3-15 15:28:40

paramecium86 发表于 2021-3-15 15:25
嗯我是一次性算的三重态和单重态。不是分别跑的MRPT2

Thanks!

zjxitcc · 发表于 Post on 2021-3-15 18:29:15

paramecium86 发表于 2021-3-15 15:25
嗯我是一次性算的三重态和单重态。不是分别跑的MRPT2

再请教一下哈，你是先跑SA-CASSCF平均2个根，然后再跑MC-QDPT呢；还是只在1个电子态下优化轨道，然后写nstate=2让GAMESS算两个根，然后再跑MC-QDPT？

我虽然知道怎么算，但是不知道一般大家的用户习惯。

paramecium86 · 发表于 Post on 2021-3-15 19:31:57

本帖最后由 paramecium86 于 2021-3-15 19:33 编辑

zjxitcc 发表于 2021-3-15 18:29
再请教一下哈，你是先跑SA-CASSCF平均2个根，然后再跑MC-QDPT呢；还是只在1个电子态下优化轨道，然后写ns ...

您客气，我用的是先SA-MCSCF做态平均之后再MCQDPT的方法。

量化小菜鸡 · 发表于 Post on 2021-4-27 14:28:02

zjxitcc 发表于 2021-3-15 11:13
具体实现方式不一样，二者不等价。不过这差距不算大（看起来一个偏高、另一个偏低），二阶微扰就这精度。 ...

ORCA好像默认IP-EA shift是个0。

zjxitcc · 发表于 Post on 2021-4-27 15:02:50

量化小菜鸡发表于 2021-4-27 14:28
ORCA好像默认IP-EA shift是个0。

嗯，是的。不过表中没列ORCA算的结果。如果用MOKIT调ORCA算CASPT2的话，会自动加上关键词CASPT2_IPEAshift 0.25

EmperorXiang · 发表于 Post on 2024-4-1 19:48:55

把计算水平提高到UCCSD(T,Full)/CBS//UCCSD/cc-pVTZ，分别优化结构，最后算出来也是1.27eV，居然反过来高了不少。结合UCCSD(T,Full)/CBS//UCCSD/cc-pVTZ可以精确计算O2的键能这件事情，它确实刻画不好单线态氧。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 氧气激发能的计算与实验差距较大