请问如何从组态系数和分子轨道得到自然轨道

zhanfei

请问如何从组态系数和分子轨道得到自然轨道？
自己的理解是
Phi^b_j是CSF





激发态计算得到本征值为激发态的能量，本征向量是组态系数也就是激发态波函数用CSF展开的系数.
找到CSF用MO展开的展开系数加上组态系数就得到这个激发态的波函数，然后再出自然轨道。
不知道自己理解是否有问题。
组态系数是程序的常规输出,但是CSF如何输出？


谢谢
此致

sobereva

输出组态系数的时候程序直接就会给出对应的CSF定义
但有的程序给出的往往是行列式的系数

有了组态系数和CSF表达式，就去计算密度矩阵（利用slater-condon规则），然后做对角化就得到了自然轨道

zhanfei

sobereva 发表于 2016-10-2 19:05
输出组态系数的时候程序直接就会给出对应的CSF定义
但有的程序给出的往往是行列式的系数

sobereva老师你好.请问G09 ORCA能输出CSF和全部组合系数吗？
G09我查了rwf 635R是组态系数但是发现数量不对比如B3LYP/6-31G TD(nstates=5,root=1)设置计算单个氧原子0 1的结果rwf 635R中数据1836个，但是Gaussian 09 Programmer's Reference上没明确的说明，而且log中显示的数值较大的组态系数值rwf文件中也找不到。log文件中没有CSF的明确信息 谢谢
输入文件 附件为rwf 635R
%rwf=c:\ex.rwf
#p B3LYP/6-31G TD(nstates=5,root=1) out=wfn

ex

0 1
O 0.0 0.0 0.0

c:\ex.wfn

sobereva

zhanfei 发表于 2016-10-2 20:58
sobereva老师你好.请问G09 ORCA能输出CSF和全部组合系数吗？
G09我查了rwf 635R是组态系数但是发现数量 ...

你用不着去折腾rwf文件
默认把组态系数绝对值大于0.1的输出，这里提到了
电子激发任务中轨道跃迁贡献的计算
http://sobereva.com/230
比如15->16 0.2就代表15->16的单激发CSF的组态系数是0.2

用IOp(9/40=x)可将组态系数绝对值大于10^-x的输出出来。

PS：对于单重态激发态计算，实际上Gaussian输出的是行列式系数不是组态系数。