用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

后面我重新走了一个，也不知道可不可以算出

%mem=2GB
%NprocShared=8
%LindaWorker=cu30-c-ib
%chk=nto-abs-air.chk
%nproc=8
# cis/6-31+G* pop(saveNTO,NTO) density=transition=1 td(nstates=20,singlets) gfinput iop(6/7=3)

sobereva

仔细看此文，你的输入文件明显不合要求
使用Multiwfn做自然跃迁轨道(NTO)分析
http://sobereva.com/377

最重要的9/40=4都没写，反倒写了完全多余的gfinput iop(6/7=3)

2L的关键词不是用来结合Multiwfn做NTO分析的

mly12345

sobereva 发表于 2018-6-5 10:19
仔细看此文，你的输入文件明显不合要求
使用Multiwfn做自然跃迁轨道(NTO)分析
http://sobereva.com/377

我之前以为旧的fchk文件可以，新写的那个可以吗老师？

mly12345

我看到了实例是这个，觉得和我最开始那个输入文件差别不大，所以就直接拿mulitiwfn试了一下

%chk=C:\uracil.chk
#p PBE1PBE/6-31G* TD IOp(9/40=4)
B3LYP/6-31G* opted

sobereva

mly12345 发表于 2018-6-5 10:20
我之前以为旧的fchk文件可以，新写的那个可以吗老师？

什么新写的？
严格按照我博文里的例子写关键词

TD任务直接产生的fch文件直接拿来用

如果不熟悉流程，先把博文里的例子重复一遍，几分钟的事

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

sobereva

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

这个可以

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 16:08
老师，我的FCHK文件算NTO结果：
S1：

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

让你变成回忆

如果不使用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了。

mly12345

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

谢谢指点！

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 ...

你好，我用

%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那么多
先严格重复一次博 ...

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