计算化学公社

 找回密码 Forget password
 注册 Register
Views: 139|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 求助审稿意见中关于Franck-Condon因子的问题

[复制链接 Copy URL]

83

帖子

0

威望

1263

eV
积分
1346

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师,大家好,我的一个关于丙酸乙酯(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








000001-FC.gjf

2.52 KB, 下载次数 Times of downloads: 2

000001-FC.log

332 KB, 下载次数 Times of downloads: 4

338

帖子

0

威望

1107

eV
积分
1445

Level 4 (黑子)

2#
发表于 Post on 2024-11-7 09:27:00 | 只看该作者 Only view this author
不同振动能级之间的跃迁概率和对应的F-C因子有关,也许对你寻找光电离反应中间态的几何结构选择和自由能变化上可以得到一些分析结果。
在Gaussian里如何计算FC因子的问题得靠其他高手来解释了。

689

帖子

21

威望

5019

eV
积分
6128

Level 6 (一方通行)

3#
发表于 Post on 2024-11-7 10:07:40 | 只看该作者 Only view this author
FC不是必须算的。可以通过垂直激发做粗糙的估计,也就是中性分子基态平衡结构到阳离子的激发/电离过程,不考虑结构变化。实验条件下,可能会在垂直激发到阳离子某个振动能级附近的+/-3个振动能级都有分布

83

帖子

0

威望

1263

eV
积分
1346

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2024-11-8 20:45:19 | 只看该作者 Only view this author
cokie 发表于 2024-11-7 09:27
不同振动能级之间的跃迁概率和对应的F-C因子有关,也许对你寻找光电离反应中间态的几何结构选择和自由能变 ...

好的,非常感谢您的回复!

83

帖子

0

威望

1263

eV
积分
1346

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2024-11-8 20:49:01 | 只看该作者 Only view this author
beefly 发表于 2024-11-7 10:07
FC不是必须算的。可以通过垂直激发做粗糙的估计,也就是中性分子基态平衡结构到阳离子的激发/电离过程,不 ...

老师,感谢您的回复,请问您提到的垂直激发,也就是第二个结构和第一个结构保持不变,然后电荷和自旋多重度对应 1 2,加上激发态的命令语句td=(nstates=6,root=1),我试算了一下,还是出现了第三个结构虚频的问题,想知道为什么会出现这样的报错呢?

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

6#
发表于 Post on 2024-11-9 07:20:42 | 只看该作者 Only view this author
人工检查1 2那个任务振动分析后到底有没有虚频,倘若有虚频,必然此处用的结构没有经过合理优化
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

83

帖子

0

威望

1263

eV
积分
1346

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2024-11-10 14:12:16 | 只看该作者 Only view this author
sobereva 发表于 2024-11-9 07:20
人工检查1 2那个任务振动分析后到底有没有虚频,倘若有虚频,必然此处用的结构没有经过合理优化

卢老师,检查了1 2任务的振动频率是没有虚频的,也优化了好几次,确认了任务3调用的确实是任务2的.chk,但是3还是出现了两个虚频,所以不知道哪里出现了问题

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 05:48 , Processed in 0.176122 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list