|
|
各位老师,大家好,我的一个关于丙酸乙酯(EP)光电离解离的工作,对中性分子EP和光电离得到的分子离子EP+进行了LBO分析,但是审稿人提出了这样一个问题:Sec 3.1 and 3.2. It is clear there is a large difference between the adiabatic and vertical ionization energy of EP. This speaks to a large geometry difference between the neutral and ion. Some of this is discussed with the changing bond lengths. The authors could easily do a Franck Condon simulation (in Gaussian) to see the vibrational excitation that occurs upon ionization and relate this to the changing geometry. It may provide more insight.
审稿人建议进行Franck-Condon因子的计算和分析。针对这个问题,我用Gaussian软件进行了Franck-Condon因子的计算,完整的输入文件和输出文件如附件所示(另外调用的输入文件是第二个结构优化得到的00001-.chk)。
我现在的问题是,
一、审稿人让我补算的具体是Franck-Condon因子的哪些分析呢?
二、对于输入文件部分,我的第二个文件的结构写成分子离子的结构是否正确?我的命令语句是否存在问题?
三、我按照当前的输入文件格式,出现了输出的高斯报错,结构中有两个虚频,我重算了好几遍,没找到出现问题的原因,请问出现报错可能是我的哪些操作有问题呢?
非常感谢各位老师,希望能得到各位老师们的指教!
输入文件:
%chk=000001.chk
#p freq=savenm wb97xd/6-311++g(2df,2p)
000001
0 1
XX
--link1--
%chk=000001-.chk
#p freq=savenm wb97xd/6-311++g(2df,2p) guess=read
000001-
1 2
XX
--link1--
%oldchk=000001.chk
%chk=fc.chk
#p freq=(fc,ReadFCHT) geom=check
FC
0 1
initial=source=calc final=source=chk
spectrum=(lower=-500.,upper=+5000.,grain=1.,HWHM=150)
temperature=value=20
reddim=clearimfreq
prescreening=(maxc1=30.,maxc2=30)
print=matrix=JK
000001-.chk
|
|