计算化学公社

标题: 求助：态特定法激发态计算输入文件格式问题(http://sobereva.com/314) [打印本页]

作者
Author: liaorongbao 时间: 2024-4-2 08:22
标题: 求助：态特定法激发态计算输入文件格式问题(http://sobereva.com/314)
执行卢老师博文链接http://sobereva.com/314的节选部分：

4.5 基于特定态溶剂模型做TDDFT计算

(4)计算溶液中的荧光发射能：

第一步：在(3)得到的结构下，获得平衡溶剂下的激发态能量，并将激发态溶剂信息写入chk：# PBE1PBE/基组 TD scrf(externaliteration,read,solvent=ethanol)，末尾空一行写noneq=write。

对于存在赝势基组的情况，我不知道输入格式怎么写。对于不存在赝势基组的情况并无疑惑。对应Lanl2DZ基组的案例文件liaotestwriteLanl2DZ.gjf的输入格式如下，出错信息见后边。文件见liaotestwriteLanl2DZ.gjf和liaotestwriteLanl2DZ.log。请各位老师给出针对这一步计算的正确的输入文件格式。
%chk=liaotestwriteLanl2DZ.chk
#P MN15/genecp  scrf(externaliteration,read,solvent=dimethylsulfoxide)  TD(nstates=3)  IOp(9/40=4) nosymm

test

0 1
Zn             -14.85184619 4.02136229 2.47722318
H             -13.30184619 4.02136229 2.47722318
H             -16.40184619 4.02136229 2.47722318

H 0
6-311G*
****
Cu 0
Lanl2DZ
****

Cu 0
Lanl2DZ

noneq=write

liaotestwriteLanl2DZ.log末尾的出错信息如下：
(Enter /soft/g16/l124.exe)
SC-PCM: Self-consistent PCM reaction field calculation.
Using the following non-standard input for PCM:
H 0
------------------------------------------------------------------------------
The following wrong keyword(s) have been found:
H
0
------------------------------------------------------------------------------
Error termination via Lnk1e in /soft/g16/l124.exe at Tue Apr  2 15:58:23 2024.

其他赝势基组，比如SDD赝势，输入文件格式这样写也会出错，见附件liaotestwrite.gjf和liaotestwrite.log

作者
Author: liaorongbao 时间: 2024-4-2 16:15
以上出错文件格式修改好了。我有两处错误，一是把分子说明的Zn原子写成了后边的Cu原子。二是Noneq=write的位置不对，应该前置到分子说明后边。附件其他文件也有类似错误，比如赝势中有S原子而分子说明却没有。另外，在群里问问题时也有这个错误，非常sorry，完全梦游。修改如下：
%chk=liaotestwriteLanl2DZ.chk
#P MN15/genecp  scrf(externaliteration,read,solvent=dimethylsulfoxide)  TD(nstates=3)  IOp(9/40=4) nosymm

test

0 1
Zn             -14.85184619 4.02136229 2.47722318
H             -13.30184619 4.02136229 2.47722318
H             -16.40184619 4.02136229 2.47722318

noneq=write

H 0
6-311G*
****
Zn 0
Lanl2DZ
****

Zn 0
Lanl2DZ

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