计算化学公社

标题: 按http://sobereva.com/314测试带有赝势的特定态法激发态计算出错 [打印本页]

作者
Author: liaorongbao 时间: 2019-4-12 12:17
标题: 按http://sobereva.com/314测试带有赝势的特定态法激发态计算出错
我的第一步获得基态平衡溶剂信息输入格式如下：
%chk=test1.chk
#P CAM-B3LYP/genecp  scrf=(read,solvent=n,n-DiMethylAcetamide)

the following coordinate is already optimizaed

0 1
C                1.30099500 1.07021700 0.00000000
H                2.05294200 0.27505800 0.00000000
H                1.65145900 2.10937300 0.00000000
C                0.00000000 0.83140900 0.00000000
H                -0.83298100 1.54076700 0.00000000
Cl             -0.62808200 -0.90205600 0.00000000

Cl 0
lanl2dz
****
C H 0
STO-3G
****

Cl 0
lanl2dz

noneq=write

第二步获得激发态在非平衡溶剂(基态溶剂)下的能量的输入格式如下：
%chk=test2.chk    !!!!!!!!!!!!!这个test2.chk是直接cp上面那个test1.chk获得的
#P CAM-B3LYP/genecp td=(singlets,nstates=3,root=1)
scrf(externaliteration,read,solvent=n,n-DiMethylAcetamide) geom=check guess=read

the test2.chk is the cp of test1.chk

0 1

C H 0
STO-3G
****
Cl 0
lanl2dz
****

Cl 0
lanl2dz

noneq=read

最后的输出信息，似乎也没明确给出在什么地方出错
99/5=1,9=1/99;
Leave Link 1 at Fri Apr 12 12:11:25 2019, MaxMem=          0 cpu:       0.0
(Enter E:\Program Files (x86)\G09W\l101.exe)
------------------------------------
the test2.chk is the cp of test1.chk
------------------------------------
Structure from the checkpoint file:  "test2.chk"
Charge =  0 Multiplicity = 1
Z-Matrix found in file.
C,0,1.300995,1.070217,0.
H,0,2.052942,0.275058,0.
H,0,1.651459,2.109373,0.
C,0,0.,0.831409,0.
H,0,-0.832981,1.540767,0.
Cl,0,-0.628082,-0.902056,0.
Recover connectivity data from disk.
NAtoms=    6 NQM=       6 NQMF=    0 NMMI=    0 NMMIF=    0
            NMic=    0 NMicF=    0.
                  Isotopes and Nuclear Properties:
(Nuclear quadrupole moments (NQMom) in fm**2, nuclear magnetic moments (NMagM)
  in nuclear magnetons)

  Atom       1          2          3          4          5          6
IAtWgt=       12          1          1       12          1       35
AtmWgt=  12.0000000 1.0078250 1.0078250  12.0000000 1.0078250  34.9688527
NucSpn=          0          1          1          0          1          3
AtZEff=  -3.6000000  -1.0000000  -1.0000000  -3.6000000  -1.0000000-204.0000000
NQMom= 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000  -8.1650000
NMagM= 0.0000000 2.7928460 2.7928460 0.0000000 2.7928460 0.8218740
AtZNuc= 6.0000000 1.0000000 1.0000000 6.0000000 1.0000000  17.0000000
Leave Link  101 at Fri Apr 12 12:11:25 2019, MaxMem= 33554432 cpu:       0.0

附http://sobereva.com/314链接包含如下内容：
4.5 基于特定态溶剂模型做TDDFT计算

这个例子包含4个任务，你需要研究什么问题就follow哪个，没有顺序关系。此例是通过特定态方法，靠PCM溶剂模型表现乙醇环境。

(1)溶剂下基态几何结构优化：和4.4节的(1)等同。

(2)计算溶液中基态到S1态的激发能：
特定态方法计算激发能时，默认情况下激发态处于平衡溶剂下，但实际中应当处于非平衡溶剂下（慢部分对基态是平衡的），因此需要做两步才能正确表现非平衡溶剂效应
  第一步：在(1)得到的结构下，先在溶剂下做基态单点计算，把基态的溶剂信息写到chk里，同时得到基态在平衡溶剂下的能量（SCF Done后面的值）：
# PBE1PBE/6-311G* scrf(read,solvent=ethanol)，末尾空一行写noneq=write
  第二步：做电子激发计算，并读取上一步在chk里写入的基态溶剂信息，由此得到激发态在非平衡溶剂下的能量：
# PBE1PBE/6-311G* TD scrf(externaliteration,read,solvent=ethanol) geom=check guess=read，末尾空一行写noneq=read
这一步计算耗时会很长，会反复计算许多次激发态，这是为了让溶剂场与S1态电子密度相互自洽（即外迭代过程）。输出文件末尾的下面这条信息告诉了你特定态方法+非平衡溶剂下的激发态的能量：
After PCM corrections, the energy is  -153.527526510    a.u.
将这个能量减去第一步的能量就是所需的激发能了。

(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
然后将第一步的能量减去第二步的能量就是荧光发射能。

作者
Author: sobereva 时间: 2019-4-12 14:20
博文内容没必要复制到这里

如果末尾就是你贴的那样，程序莫名其妙就断了，换个版本或平台再试

作者
Author: liaorongbao 时间: 2019-4-12 21:42
Sob老师好，尾部noneq=read和尾部基组调换位置似乎可以执行。即为noneq=read放在前而基组放在后，也许仅仅是个顺序问题。G09的D01的win和linux都测试了，应该是这样的。实际上noneq=read应该后于基组写出才感觉合理

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