求助MOLPRO CASSCF 打印每个电子态的轨道信息 及 保存密度矩阵

leogrimm

请教大家，在 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
}

谢谢！

jiangning198511

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