关于用Gaussian寻找圆锥交叉点的关键词问题求助

陈小北

Sob 老师好


我有一个疑惑，我想知道在找圆锥交叉点的时候，关键词写成如下这种：


#p opt=conical casscf(6,6,nroot=2)/6-31g(d) guess=read nosymm


title card


0 1
c     -0.00080300    1.41440000    0.65507800
...


0.500000000.50000000


我知道写opt=conical，此时CASSCF默认为态平均CASSCF，所以需要指定各态权重，关键词写的nroot=2，所以是两个态，每个是0.5，此时是在找S0态和S1态之间的交叉点，但是我如果想找S1态和S2态之间的交叉点应该怎么写关键词呢？关键词写nroot=3，那结尾赋权重是每个态0.33333333吗？这样是默认在找三个态（S0、S1和S2）之间的交叉点吗？

zjxitcc

你可以试试平均3个态，然后权重写0, 0.5, 0.5（具体格式自己调整）

pwzhou

是的，nroot=3，权重0.33，找的就是S1和S2之间的交叉点

陈小北

pwzhou 发表于 2020-12-16 16:09
是的，nroot=3，权重0.33，找的就是S1和S2之间的交叉点

老师您好！

看过您很多关于光激发计算的文章，对于您文章中提到的LIIC分析方法十分感兴趣，想知道您找圆锥交点用的是高斯吗？还是有其他别的软件呢？

我最近在学这一部分的内容，但是用高斯做CAS发现很不容易收敛，所以想要问您一下

非常感谢！

陈小北

zjxitcc 发表于 2020-12-16 15:46
你可以试试平均3个态，然后权重写0, 0.5, 0.5（具体格式自己调整）

老师好！

我关注了您的公众号！！而且我学CAS就是在里面学的！！受益匪浅！！非常感谢您深入浅出的指导！以及您为后辈指路所付出的努力！！

祝好！

劲仔小鱼干

zjxitcc 发表于 2020-12-16 15:46
你可以试试平均3个态，然后权重写0, 0.5, 0.5（具体格式自己调整）

感谢您的回复! 我在您的公众号里也收获到了很多！

pwzhou

陈小北 发表于 2020-12-16 16:23
老师您好！

看过您很多关于光激发计算的文章，对于您文章中提到的LIIC分析方法十分感兴趣，想知道您找 ...

我用的Firefly、molpro和OpenMolcas，Firefly的计算速度和收敛情况是最好的，但是这个软件不太容易使用。Gaussian没用过，活化空间超过(10,8)以后Gaussian就基本不行了

陈小北

pwzhou 发表于 2020-12-16 18:07
我用的Firefly、molpro和OpenMolcas，Firefly的计算速度和收敛情况是最好的，但是这个软件不太容易使用。 ...

老师好

非常感谢您的指导，我现在在学习molpro

还有问题想要请教您一下，我理解的用高斯寻找圆锥交叉点相关流程是这样的：（感谢二楼的zjxitcc老师）

1.优化基态结构，为CASSCF选取合适的活性空间

2.用CASSCF优化基态结构，在挑选了合适的活性空间后，做结构优化及频率计算

3.用CASSCF计算第一激发态，获得垂直激发能

4.用CASSCF优化激发态的结构并做频率计算。

5.寻找圆锥交叉点。可以使用的方法是利用激发态下柔性扫描的结果作为优化交叉点的初始结构。

以上是寻找S1和S0态之间的圆锥交叉点，如果想要寻找S2和S1态之间的圆锥交叉点，还需要优化S2态的结构，以及第5步的初始结构可以在S1优化之后的结构和S2优化之后的结构的中间

这都是我理解的步骤，我想问您的是，如果这是对的，那么用molpro做交叉点的寻找过程也是同样的步骤吗？

pwzhou

陈小北 发表于 2020-12-16 21:32
老师好

非常感谢您的指导，我现在在学习molpro

一样的，这个流程和软件无关