Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法

wq8484 · 发表于 Post on 2017-4-28 18:49:50

请问如果想计算中性有机分子dimer的磷光，应该如何用UDFT优化T1结构呢？自旋多重度是设为3还是5呢？

赵云跳槽 · 发表于 Post on 2017-4-30 19:32:14

wq8484 发表于 2017-4-28 18:49
请问如果想计算中性有机分子dimer的磷光，应该如何用UDFT优化T1结构呢？自旋多重度是设为3还是5呢？

可以分别使用自选多重读3, 5. 7 等试算一把，哪个能量最低哪个就是T1态

wq8484 · 发表于 Post on 2017-5-8 09:47:09

赵云跳槽发表于 2017-4-30 19:32
可以分别使用自选多重读3, 5. 7 等试算一把，哪个能量最低哪个就是T1态

多谢多谢！

daodaoliu · 发表于 Post on 2017-6-12 11:54:32

sob老师，按照你的方法，做荧光的计算，前面的计算结果都与实验结果相符了，但在，4.5中（4）中所述，计算溶液中的荧光发射能的第二步：做基态计算时，在gview中打开上一步的log文件，生成基态计算的gif文件，提交以后出现错误：
Z-Matrix unexpectedly found in input stream.
Error termination via Lnk1e in /public/software/gauss_09/g09/l101.exe at Sun Jun 11 23:17:54 2017.
请教老师

如下：
%chk=A_C_NH_5.chk
%nprocshared=2
%nprocl=8
%mem=8GB
#P Pbe1pbe/6-31g*
scrf=(solvent=generic,read) geom=allcheck guess=read

A_C_NH EM

0 1
C                4.67564200 2.01879500 0.05511500
C                3.40736200 2.36323900 0.48885400
C                2.36743400 1.43511300 0.51850000
N                2.64431000 0.16549200 -0.02040100
C                3.92859200 -0.24044300 -0.31831200
C                4.97004800 0.65255700 -0.30945800
C                1.83358000 -0.96874300 -0.00461400
C                2.70947800 -2.03304300 -0.33441200
N                3.93944000 -1.60716400 -0.54987400
C                0.42270900 -0.95679000 0.13146100
C                -0.53207800 0.03287400 -0.12561200
C                -1.77473300 -0.57482100 0.12446400
C                -3.12299700 0.00638100 -0.00410800
N                -1.60767000 -1.84603100 0.51169300
C                -4.25835900 -0.81032900 0.04218100
C                -5.52897100 -0.25933100 -0.07885100
C                -5.68576500 1.11461600 -0.25194400
C                -4.56069100 1.93374800 -0.30054500
C                -3.28860400 1.38424500 -0.17614800
N                -0.29389100 -2.06327600 0.49721000
H                5.46744300 2.75881500 0.04041500
H                3.20171300 3.36930400 0.83784600
H                1.36739400 1.61912900 0.87655000
H                5.97049100 0.31855200 -0.54769000
H                2.42744700 -3.06879700 -0.47229100
H                -0.34029100 1.02063900 -0.51246700
H                -4.13836200 -1.88028200 0.17004500
H                -6.39997300 -0.90581100 -0.04189900
H                -6.67796000 1.54348900 -0.34808000
H                -4.67158200 3.00528000 -0.43184300
H                -2.42249700 2.03719300 -0.20494900
H                0.06497800 -2.94906500 0.82454800

eps=48.52765
epsinf=1.8378245

noneq=read

sobereva · 发表于 Post on 2017-6-12 13:54:27

daodaoliu 发表于 2017-6-12 11:54
sob老师，按照你的方法，做荧光的计算，前面的计算结果都与实验结果相符了，但在，4.5中（4）中所述，计算 ...

既然写了allcheck就别自己写坐标
为了避免麻烦，甭用allcheck，直接在输入文件里写明坐标
epsinf后面别有空行

daodaoliu · 发表于 Post on 2017-6-12 15:20:43

本帖最后由 daodaoliu 于 2017-6-12 15:25 编辑

sobereva 发表于 2017-6-12 13:54
既然写了allcheck就别自己写坐标
为了避免麻烦，甭用allcheck，直接在输入文件里写明坐标
epsinf后面别 ...

输入写成下面这样，还是出现了同样的错误，是不是上一步chk文件被修改了导致。没有读取到相应的信息？另外，上一步的chk文件名为A_C_NH_4.chk 。
%chk=A_C_NH_5.chk
%nprocshared=2
%nprocl=8
%mem=8GB
#P Pbe1pbe/6-31g*
scrf=(solvent=generic,read) geom=check guess=read

A_C_NH EM
0 1
C                4.67564200 2.01879500 0.05511500
C                3.40736200 2.36323900 0.48885400
C                2.36743400 1.43511300 0.51850000
N                2.64431000 0.16549200 -0.02040100
C                3.92859200 -0.24044300 -0.31831200
C                4.97004800 0.65255700 -0.30945800
C                1.83358000 -0.96874300 -0.00461400
C                2.70947800 -2.03304300 -0.33441200
N                3.93944000 -1.60716400 -0.54987400
C                0.42270900 -0.95679000 0.13146100
C                -0.53207800 0.03287400 -0.12561200
C                -1.77473300 -0.57482100 0.12446400
C                -3.12299700 0.00638100 -0.00410800
N                -1.60767000 -1.84603100 0.51169300
C                -4.25835900 -0.81032900 0.04218100
C                -5.52897100 -0.25933100 -0.07885100
C                -5.68576500 1.11461600 -0.25194400
C                -4.56069100 1.93374800 -0.30054500
C                -3.28860400 1.38424500 -0.17614800
N                -0.29389100 -2.06327600 0.49721000
H                5.46744300 2.75881500 0.04041500
H                3.20171300 3.36930400 0.83784600
H                1.36739400 1.61912900 0.87655000
H                5.97049100 0.31855200 -0.54769000
H                2.42744700 -3.06879700 -0.47229100
H                -0.34029100 1.02063900 -0.51246700
H                -4.13836200 -1.88028200 0.17004500
H                -6.39997300 -0.90581100 -0.04189900
H                -6.67796000 1.54348900 -0.34808000
H                -4.67158200 3.00528000 -0.43184300
H                -2.42249700 2.03719300 -0.20494900
H                0.06497800 -2.94906500 0.82454800

log文件最后的几行如下：
A_C_NH EM
---------
Z-Matrix unexpectedly found in input stream.
Error termination via Lnk1e in /public/software/gauss_09/g09/l101.exe at Mon Jun 12 15:10:20 2017.
Job cpu time:  0 days  0 hours  0 minutes  0.7 seconds.
File lengths (MBytes):  RWF=    5 Int=    0 D2E=    0 Chk=    1 Scr=    1

sobereva · 发表于 Post on 2017-6-12 15:25:06

daodaoliu 发表于 2017-6-12 15:20
输入写成下面这样，还是出现了同样的错误，是不是上一步chk文件被修改了导致。没有读取到相应的信息？
% ...

你在末尾都没写自定义的与溶剂相关的信息，干嘛写solvent=generic,read

daodaoliu · 发表于 Post on 2017-6-12 15:29:13

sobereva 发表于 2017-6-12 15:25
你在末尾都没写自定义的与溶剂相关的信息，干嘛写solvent=generic,read

例子里面用的是乙醇，我的溶剂是一个混合溶剂，所以输入了eps和epsinf，是这一步不需要输入这两个数据吗？直接由上一步的chk文件读取？

sobereva · 发表于 Post on 2017-6-12 15:34:43

daodaoliu 发表于 2017-6-12 15:29
例子里面用的是乙醇，我的溶剂是一个混合溶剂，所以输入了eps和epsinf，是这一步不需要输入这两个数据吗 ...

得输入啊，你66L的输入文件里又没写

daodaoliu · 发表于 Post on 2017-6-12 15:41:18

sobereva 发表于 2017-6-12 15:34
得输入啊，你66L的输入文件里又没写

抱歉sob老师，最后写写了的，贴的时候漏掉了。

daodaoliu · 发表于 Post on 2017-6-12 15:57:21

daodaoliu 发表于 2017-6-12 15:41
抱歉sob老师，最后写写了的，贴的时候漏掉了。

我重新贴一下输入和输出文件，请sob老师指点迷津
输入文件：

%chk=A_C_NH_4.chk
%nprocshared=2
%nprocl=8
%mem=8GB
#P Pbe1pbe/6-31g*
scrf=(solvent=generic,read) geom=check guess=read

A_C_NH EM

0 1
C                4.67564200 2.01879500 0.05511500
C                3.40736200 2.36323900 0.48885400
C                2.36743400 1.43511300 0.51850000
N                2.64431000 0.16549200 -0.02040100
C                3.92859200 -0.24044300 -0.31831200
C                4.97004800 0.65255700 -0.30945800
C                1.83358000 -0.96874300 -0.00461400
C                2.70947800 -2.03304300 -0.33441200
N                3.93944000 -1.60716400 -0.54987400
C                0.42270900 -0.95679000 0.13146100
C                -0.53207800 0.03287400 -0.12561200
C                -1.77473300 -0.57482100 0.12446400
C                -3.12299700 0.00638100 -0.00410800
N                -1.60767000 -1.84603100 0.51169300
C                -4.25835900 -0.81032900 0.04218100
C                -5.52897100 -0.25933100 -0.07885100
C                -5.68576500 1.11461600 -0.25194400
C                -4.56069100 1.93374800 -0.30054500
C                -3.28860400 1.38424500 -0.17614800
N                -0.29389100 -2.06327600 0.49721000
H                5.46744300 2.75881500 0.04041500
H                3.20171300 3.36930400 0.83784600
H                1.36739400 1.61912900 0.87655000
H                5.97049100 0.31855200 -0.54769000
H                2.42744700 -3.06879700 -0.47229100
H                -0.34029100 1.02063900 -0.51246700
H                -4.13836200 -1.88028200 0.17004500
H                -6.39997300 -0.90581100 -0.04189900
H                -6.67796000 1.54348900 -0.34808000
H                -4.67158200 3.00528000 -0.43184300
H                -2.42249700 2.03719300 -0.20494900
H                0.06497800 -2.94906500 0.82454800

eps=48.52765
epsinf=1.8378245
noneq=read

输出文件最后几行：
A_C_NH EM
---------
Z-Matrix unexpectedly found in input stream.
Error termination via Lnk1e in /public/software/gauss_09/g09/l101.exe at Mon Jun 12 15:10:20 2017.
Job cpu time:  0 days  0 hours  0 minutes  0.7 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1

daodaoliu · 发表于 Post on 2017-6-12 17:11:50

本帖最后由 daodaoliu 于 2017-6-12 17:13 编辑

sob老师，只做到 (4)计算溶液中的荧光发射能中的第一步，查找log文件中的
Excited state 1：
最后一条对应发射。
不做第二步的计算可以吗？

sobereva · 发表于 Post on 2017-6-12 17:43:52

daodaoliu 发表于 2017-6-12 15:57
我重新贴一下输入和输出文件，请sob老师指点迷津
输入文件：

你既然输入文件里已经有了坐标，就不要再用geom=check

sobereva · 发表于 Post on 2017-6-12 17:45:58

daodaoliu 发表于 2017-6-12 17:11
sob老师，只做到 (4)计算溶液中的荧光发射能中的第一步，查找log文件中的
Excited state 1：
最后一条 ...

你是指就做“计算溶液中的荧光发射能：直接取(3)输出文件当中最后一次的S0->S1激发能即可”这一步？如果你只需要得到荧光发射能，这就够了

daodaoliu · 发表于 Post on 2017-6-12 18:38:03

sobereva 发表于 2017-6-12 17:45
你是指就做“计算溶液中的荧光发射能：直接取(3)输出文件当中最后一次的S0->S1激发能即可”这一步？如果 ...

(3)中得到的发射波长误差接近50nm，因此，又继续考虑溶剂的影响，继续往下做了(4)中的第一步，log文件中最后一条excited state 1中的波长跟实验值只相差9nm，第二步基态计算出错。我想就用(4)的第一步得到的波长作为发射波长可以吗？

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法