计算化学公社

标题: 关于“基于特定态溶剂模型做TDDFT计算”的几个问题 [打印本页]

作者
Author: Luadamcy 时间: 2017-4-12 22:52
标题: 关于“基于特定态溶剂模型做TDDFT计算”的几个问题
求助各位大神，在sob老师写的“Gaussian中用TDDFT计算激发态的吸收、荧光和磷光光谱“一文中，介绍了两种加入溶剂模型的计算荧光（磷光）的方法，我先用线性响应溶剂模型计算后发现与气相结果相差不大，无法解释实验现象，故尝试用基于特定态的溶剂模型来计算，但是文中相应段落的第四步不是很明白
(3)优化溶液中的S1结构：和4.4节的(3)等同。特定态方法下虽然也能做几何优化，但巨慢，所以优化激发态还是应当用线性响应模型。

(4)计算溶液中的荧光发射能：
第一步：在(3)得到的结构下，获得平衡溶剂下的激发态能量，并将激发态溶剂信息写入chk：
# PBE1PBE/6-311G* TD scrf(externaliteration,read,solvent=ethanol)，末尾空一行写noneq=write
激发态能量应当取末尾输出的After PCM corrections后面的能量值。
第二步：做基态计算，读取第上一步在chk里写入的激发态溶剂信息，得到非平衡溶剂下的基态能量（SCF Done后面的值）：
#P PBE1PBE/6-311G* scrf(read,solvent=ethanol) geom=check guess=read，末尾空一行写noneq=read
然后将第一步的能量减去第二步的能量就是荧光发射能。

如上，我的问题是第一，如何把激发态溶剂信息写入chk，这里指的chk是（3）步中优化结果后得到的chk文件吗？是否也要把改好的chk文件与gjf文件放在一起才能完成这步计算；

第二，第一步与第二步是写在一个gjf文件里，有先后之分，还是分成两个gjf文件分别计算？具体的关键词应该如何写？

计算小白，请大神详细解释一下

作者
Author: sobereva 时间: 2017-4-13 02:14
1 noneq=write就说明写入到chk里，chk是指当前这一步产生的chk。
只有这段文字里的“第二步”才需要“第一步”产生的chk。“第一步”不需要之前就有chk，毕竟这一步又没从chk里读取任何信息
2 分别做，写两个gjf，做完第一个做第二个。关键词文中不都给了，就是#后面的

作者
Author: Luadamcy 时间: 2017-4-13 09:37
本帖最后由 Luadamcy 于 2017-4-13 09:38 编辑

sobereva 发表于 2017-4-13 02:14
1 noneq=write就说明写入到chk里，chk是指当前这一步产生的chk。

那读取上一步在chk里写入的激发态溶剂信息具体是什么意思呢？

作者
Author: 冰释之川 时间: 2017-4-13 10:13
本帖最后由冰释之川于 2017-4-13 10:21 编辑

Luadamcy 发表于 2017-4-13 09:37
那读取上一步在chk里写入的激发态溶剂信息具体是什么意思呢？

“上一步”是指 “# PBE1PBE/6-311G* TD scrf(externaliteration,read,solvent=ethanol)，末尾空一行写noneq=write”下计算完的chk文件，此时的溶剂分子对于处于平衡激发态的溶质分子来说是完全平衡的（溶剂的快慢部分都对激发态响应）

作者
Author: Luadamcy 时间: 2017-4-13 11:10

冰释之川发表于 2017-4-13 10:13
“上一步”是指 “# PBE1PBE/6-311G* TD scrf(externaliteration,read,solvent=ethanol)，末尾空一行写no ...

谢谢您的回复！
我刚试着按照步骤计算了一下（4）中的第一步，出现了“There are no atoms!
Error termination via Lnk1e in /usr/local/g09/l101.exe at Thu Apr 13 12:02:24 2017.
”的错误，我的计算文件如下：
%nprocshared=8
%mem=20GB
%chk=LpBr-100-sol-equ.chk
# wB97XD/gen pseudo=read TD(triplets) scrf(externaliteration,read,solvent=dichloromethane)

exp_6257; Picture 1

0 1
原子坐标

NonEq=write

C H O N 0
6-311G**
****
Br 0
lanl2dz
****

Br 0
lanl2dz

请问是什么问题使得没有读取原子信息呢？

作者
Author: sobereva 时间: 2017-4-13 11:42

Luadamcy 发表于 2017-4-13 11:10
谢谢您的回复！
我刚试着按照步骤计算了一下（4）中的第一步，出现了“There are no atoms!
Error ter ...

你没写guess=read当然不会读取原子坐标

作者
Author: Luadamcy 时间: 2017-4-13 14:54

sobereva 发表于 2017-4-13 11:42
你没写guess=read当然不会读取原子坐标

继续救助，已经加了guess=read，还是报错，没有读取原子坐标，计算文件如下，是原子坐标的问题吗？

%nprocshared=8
%mem=20GB
%chk=LpBr-100-sol-equ.chk
# wB97XD/gen pseudo=read Guess=read TD=(triplets) SCRF=(solvent=dichloromethane,externaliteration,read)

exp_6257; Picture 1

0 1

C -2.5985411 -1.0359586 0.1683559
C -3.2788443 -0.0475923 -0.6674372
C -2.8504297 1.1911669 -0.3058980
C -1.8734475 1.0421234 0.7708235
N -1.7654271 -0.3292148 1.0603280
C -2.7851348 -2.3673642 0.0244952
C -3.6932912 -2.8035736 -0.9995446
C -4.3836602 -1.8808469 -1.8108722
C -4.2130039 -0.5266633 -1.6727791
......

C H O N 0
6-311G**
****
Br 0
lanl2dz
****

Br 0
lanl2dz

NonEq=Write

作者
Author: Luadamcy 时间: 2017-4-13 16:36

sobereva 发表于 2017-4-13 11:42
你没写guess=read当然不会读取原子坐标

以及为什么我输入的原子坐标不被读取呢？

作者
Author: sobereva 时间: 2017-4-13 18:38

Luadamcy 发表于 2017-4-13 16:36
以及为什么我输入的原子坐标不被读取呢？

你要写guess=read就别写坐标，要写坐标就别写guess=read

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