求助Gaussian圆锥交叉优化过程

zhanfei

如题Gaussian
手册里说“预测圆锥交叉点。在计算路径中使用选项 Opt=Conical ，把使用 CASSCF 对指定态进行
优化的计算变成寻找这个态的圆锥交点或避免交叉点的计算。优化的结构将是在圆锥交叉点
或避免交叉点处的结构。”
然后找了一个小的体系当nroot=1时候报一个FILEIO错误
nroot=2时候
EIGENVALUES AND  EIGENVECTORS OF CI MATRIX
( 1)     EIGENVALUE    -154.0507499371
( 2)     EIGENVALUE    -154.0502715500
nroot=3时候
( 1)     EIGENVALUE    -154.1715170635
( 2)     EIGENVALUE    -153.8778632741
( 3)     EIGENVALUE    -153.8777496141
可以理解成是找nroot=n的态和前一个激发态的交叉点吗，
那么比如找基态和第二激发态的圆锥交叉该怎么处理？
另外程序是如何搜索圆锥交叉的希望能大致说下，给出参考文献。
谢谢大家
sincerely，

zhanfei

再追问一个问题，我找到的参考文献都是普通有机物的，金属配合物一定是避免交叉的吗，为什么找不到金属配合物圆锥交叉的文献？

sobereva

nroot=1不行，只算一个态，没法算交叉
你的nroot=2的情况奏效了，确实是交叉点（不过也建议确认一下计算过程，看看是不是两个态之间的能量随着优化逐渐减小）

nroot=3的情况，可以把各个态在优化过程中的能量变化画个曲线看看。

搜索交叉点，具体来说是搜索极小能量交叉点(MECP)，算法思路很简单，比如搜索i,j两个态的，就是让(Ej-Ei)^2和Ei同时最小化（平时做单个态的几何优化只是让一个态的能量最小化）。细节可以看看之前提到过的A Computational Strategy for Organic Photochemistry那篇，Robb写的，他也是高斯的MCSCF代码的作者。另外可以看看Harvey的Theor Chem Acc (1998) 99:95-99，交叉点的搜索是用Harvey自己的MECP程序。

关于过渡金属配合物，不必纠结圆锥交叉和避免交叉，只是一墙之隔，我也没注意到过相关研究，大抵是因为其结构又大又比较呈刚性，不像有机小分子稍微动动结构就有出现交叉的机会。

zhanfei

sobereva 发表于 2014-11-5 20:23
nroot=1不行，只算一个态，没法算交叉
你的nroot=2的情况奏效了，确实是交叉点（不过也建议确认一下计算过 ...

又大又成刚性原来如此。thanks sob sennsei

378024572

学习了。。。。。。。。。。。。。。