用mulitiwfn算NTO轨道贡献的时候，所有的激发态算出的NTO的本征值都是0

mly12345

输入文件：
%mem=2000MW
%NprocShared=8
%LindaWorker=cu30-c-ib
%chk=abs-air.chk
%nproc=8
# cam-b3lyp/6-31g(d,p) td(nstates=40,singlets) gfinput iop(6/7=3)

Title Card Required

0 1
==============================================
算出最大激发态是S4，然后用 mulitiwfn计算NTO贡献的时候，所有激发态的所有NTO贡献都是0：

Multiplicity of this excited state is   1
Excitation energy is   4.6678000 eV
There are       3 orbital pairs in this transition mode
The sum of square of excitation coefficients:  0.469354
The negative of the sum of square of de-excitation coefficients:  0.000000
The sum of above two values  0.469354
The highest 10 eigenvalues of occupied NTOs:
    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000
The highest 10 eigenvalues of virtual NTOs:
    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000
0 Return

求问，是使用的fchk文件有问题吗？

mly12345

让你变成回忆 发表于 2018-6-6 18:17
我根据你提供的体系，改了IOp(9/40=4)，分析NTO没有任何问题。下面是过程：

There are    3 transitio ...

谢谢，谢谢。

让你变成回忆

mly12345 发表于 2018-6-6 16:13
你好，具体输入文件我传了附加文件，麻烦帮忙看一下~

我根据你提供的体系，改了IOp(9/40=4)，分析NTO没有任何问题。下面是过程：

There are    3 transition modes, analyze which one?  e.g. 2
1
Multiplicity of this excited state is   1
Excitation energy is   3.6046000 eV
There are    8387 orbital pairs in this transition mode
The sum of square of excitation coefficients:  0.502523
The negative of the sum of square of de-excitation coefficients: -0.002526
The sum of above two values  0.499997

The highest 10 eigenvalues of occupied NTOs:
    0.000150    0.000162    0.000205    0.000244    0.000258
    0.000350    0.000421    0.000745    0.018735    0.981152

The highest 10 eigenvalues of virtual NTOs:
    0.981152    0.018735    0.000745    0.000421    0.000350
    0.000258    0.000244    0.000205    0.000162    0.000150

0 Return
1 Output NTO orbitals to .molden file
2 Output NTO orbitals to .fch file

mly12345

让你变成回忆 发表于 2018-6-6 15:59
上传一个完成的输入文件？

你好，具体输入文件我传了附加文件，麻烦帮忙看一下~

mly12345

让你变成回忆 发表于 2018-6-6 15:59
上传一个完成的输入文件？

%chk=oxd9-cam-abs-air.chk
#p PBE1PBE/6-31G* TD  IOp(9/40=6)

B3LYP/6-31G* opted

让你变成回忆

mly12345 发表于 2018-6-6 15:55
你好，我用

%chk=nto-air.chk

上传一个完成的输入文件？

mly12345

sobereva 发表于 2018-6-6 04:10
明显又没用IOP(9/40=4)，否则The sum of above two values后头不可能偏离0.5那么多
先严格重复一次博 ...

老师，我写了IOp的之后还是依旧这样，为什么呢？

mly12345

让你变成回忆 发表于 2018-6-6 08:08
如果不使用IOp(9/40=4)，Gaussian默认情况下只会输出组态系数比较大的组态（默认是>0.1），Gaussian计算NTO ...

你好，我用

%chk=nto-air.chk
#p PBE1PBE/6-31G* TD IOp(9/40=4)

B3LYP/6-31G* opted

mulitiwfn分析还是出不来，这又是怎么回事。

mly12345

sobereva 发表于 2018-6-6 04:10
明显又没用IOP(9/40=4)，否则The sum of above two values后头不可能偏离0.5那么多
先严格重复一次博 ...

谢谢老师

mly12345

让你变成回忆 发表于 2018-6-6 08:08
如果不使用IOp(9/40=4)，Gaussian默认情况下只会输出组态系数比较大的组态（默认是>0.1），Gaussian计算NTO ...

谢谢指点！

让你变成回忆

如果不使用IOp(9/40=4)，Gaussian默认情况下只会输出组态系数比较大的组态（默认是>0.1），Gaussian计算NTO的时候由于是调用的chk文件，而chk文件存的组态系数较多，因此不存在这个问题。 但是如果利用Multiwfn来计算的话，由于Multiwfn是从log文件中提取组态系数，因此这就需要log文件中保存的组态系数足够多，否则就会出现如你这样的问题“The sum of above two values后头偏离0.5那么远”。
PS: Multiwfn做NTO分析肯定是没有任何问题的（如果严格按照博文的步骤），自从Multiwfn有了这个功能以后，我就没用过Gaussian来做NTO了。

sobereva

mly12345 发表于 2018-6-5 16:08
老师，我的FCHK文件算NTO结果：
S1：

明显又没用IOP(9/40=4)，否则The sum of above two values后头不可能偏离0.5那么多
先严格重复一次博文里尿嘧啶的例子，确保操作过程无误，之后再算你自己的体系

mly12345

sobereva 发表于 2018-6-5 10:34
这个可以

老师，我的FCHK文件算NTO结果：
S1：

Multiplicity of this excited state is   1
Excitation energy is   4.6678000 eV
There are       3 orbital pairs in this transition mode
The sum of square of excitation coefficients:  0.469354
The negative of the sum of square of de-excitation coefficients:  0.000000
The sum of above two values  0.469354

The highest 10 eigenvalues of occupied NTOs:
    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000

The highest 10 eigenvalues of virtual NTOs:
    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000

---------------
最大振子强度是S5：
5
Multiplicity of this excited state is   1
Excitation energy is   5.4171000 eV
There are       4 orbital pairs in this transition mode
The sum of square of excitation coefficients:  0.487983
The negative of the sum of square of de-excitation coefficients:  0.000000
The sum of above two values  0.487983

The highest 10 eigenvalues of occupied NTOs:
    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000

The highest 10 eigenvalues of virtual NTOs:
    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000

这是怎么回事

sobereva

mly12345 发表于 2018-6-5 10:25
好的好的，谢谢老师指导
按照例子写的：

这个可以

mly12345

sobereva 发表于 2018-6-5 10:22
什么新写的？
严格按照我博文里的例子写关键词

好的好的，谢谢老师指导
按照例子写的：

%mem=2GB
%chk=oxd8c-nto-air.chk
#p PBE1PBE/6-31G* TD(nstates=20,singlets)  IOp(9/40=4)

B3LYP/6-31G* opted

0 1