按http://sobereva.com/314测试带有赝势的特定态法激发态计算出错

liaorongbao

我的第一步获得基态平衡溶剂信息输入格式如下：
%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
然后将第一步的能量减去第二步的能量就是荧光发射能。

sobereva

博文内容没必要复制到这里

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

liaorongbao

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