计算化学公社

标题: 用mulitiwfn算NTO轨道贡献的时候，所有的激发态算出的NTO的本征值都是0 [打印本页]

作者
Author: mly12345 时间: 2018-6-5 10:17
标题: 用mulitiwfn算NTO轨道贡献的时候，所有的激发态算出的NTO的本征值都是0
本帖最后由 mly12345 于 2018-6-6 16:13 编辑

输入文件：
%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文件有问题吗？

作者
Author: mly12345 时间: 2018-6-5 10:18
后面我重新走了一个，也不知道可不可以算出

%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)

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

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

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

作者
Author: mly12345 时间: 2018-6-5 10:20

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

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

作者
Author: mly12345 时间: 2018-6-5 10:21
我看到了实例是这个，觉得和我最开始那个输入文件差别不大，所以就直接拿mulitiwfn试了一下

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

作者
Author: sobereva 时间: 2018-6-5 10:22

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

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

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

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

作者
Author: mly12345 时间: 2018-6-5 10:25

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

作者
Author: sobereva 时间: 2018-6-5 10:34

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

这个可以

作者
Author: mly12345 时间: 2018-6-5 16:08

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

这是怎么回事

作者
Author: sobereva 时间: 2018-6-6 04:10

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

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

作者
Author: 让你变成回忆 时间: 2018-6-6 08:08
本帖最后由让你变成回忆于 2018-6-6 08:12 编辑

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

作者
Author: mly12345 时间: 2018-6-6 08:10

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

谢谢指点！

作者
Author: mly12345 时间: 2018-6-6 08:11

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

谢谢老师

作者
Author: mly12345 时间: 2018-6-6 15:55

让你变成回忆发表于 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分析还是出不来，这又是怎么回事。

作者
Author: mly12345 时间: 2018-6-6 15:56

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

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

作者
Author: 让你变成回忆 时间: 2018-6-6 15:59

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

%chk=nto-air.chk

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

作者
Author: mly12345 时间: 2018-6-6 16:11

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

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

B3LYP/6-31G* opted

作者
Author: mly12345 时间: 2018-6-6 16:13

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

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

作者
Author: 让你变成回忆 时间: 2018-6-6 18:17

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

作者
Author: mly12345 时间: 2018-6-7 20:55

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

There are 3 transitio ...

谢谢，谢谢。

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