计算化学公社

标题: 求助MOLPRO CASSCF 打印每个电子态的轨道信息及保存密度矩阵 [打印本页]

作者
Author: leogrimm 时间: 2022-4-19 03:51
标题: 求助MOLPRO CASSCF 打印每个电子态的轨道信息及保存密度矩阵
请教大家，在 casscf 计算中，应怎样设置，才可以输出每个电子态的轨道信息呢？
下面是计算 2 个 A1 态的结果，权重相等，且只设置了一个未配对的电子，因此6.1 和 7.1 占据数都是 0.5，可推知每个态各自占据了 6.1 和 7.1。但实际计算往往有更多态，很难反推结果，所以希望打印每个电子态的轨道信息。

NATURAL ORBITALS (state averaged)
=================================

  Orbital  Occupation Energy    Coefficients
1.1    2.00000 -20.81004    2 1s 1.00066
2.1    2.00000 -15.80207    1 1s 1.00036
3.1    2.00000 -1.78265    1 1s 0.55802    1 2pz 0.33771    2 1s 0.76946    2 2pz  -0.29206
4.1    2.00000 -0.98283    1 1s -0.66166    2 1s 0.62276    2 2pz 0.60806
5.1    2.00000 -0.74471    1 1s -0.53189    1 1s -0.26158    1 2pz 0.71302    2 2pz  -0.49109
6.1    0.50000    0.86091    1 1s -0.28897    1 1s -1.39774    1 1s 2.13618    1 2pz 1.27994 ...
7.1    0.50000    0.44996    1 1s 2.45415    1 2pz 0.25961    1 2pz 1.96600    2 1s -0.39833 ...

一个相关的问题是，怎样设置能够保存每个电子态的密度矩阵，输入文件中尝试使用下面设置保存的密度矩阵，也是态平均意义的结果。

{MATROP
LOAD, DMAT, DEN, 2140.2
WRITE, DMAT, DMAT
}

谢谢！

作者
Author: jiangning198511 时间: 2022-4-19 11:08
我的看法：从CASSCF计算原理来看，你做的是态平均计算，所以程序优化了分子轨道系数(活性空间)+不同态的CI系数，所以不同态对应的分子轨道系数是一样的，不同的是CI系数，所以你说的不同态对应的应该是CI系数不同而不是轨道不同

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