计算化学公社

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

[综合交流] 关于计算threshold photoelectron (TPE) spectrum的相关问题

[复制链接 Copy URL]

122

帖子

0

威望

2336

eV
积分
2458

Level 5 (御坂)

本帖最后由 my_aman 于 2019-10-22 15:32 编辑

各位老师同学好,下图是通过光电子谱来研究CN自由基的电离能,请教几个相关的问题,谢谢:
1.图中通过理论拟合出来的谱图,在高斯中是加入什么关键词实现的呢?分别计算什么呢?freq吗?
2.自由基被光子打掉电子处在不同的振动,所以有不同的峰?图中峰的强度代表什么呢?3.假如有若干个异构体都对谱图有贡献,此时怎么拟合谱图呢?



2)2[3%L_`H1}9L8CHRX5L$K.png (112.74 KB, 下载次数 Times of downloads: 12)

2)2[3%L_`H1}9L8CHRX5L$K.png

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

2#
发表于 Post on 2019-10-23 01:54:24 | 只看该作者 Only view this author
那是展宽,不叫拟合

需要对电离态和中性状态都做振动分析,参考下文。和文中的区别是要把下文中的激发态计算改成电离态
振动分辨的电子光谱的计算
http://sobereva.com/223

目前文中说的是G09的情况,G16里用法发生了不小变化,见手册freq关键词的说明

看了文章就知道为什么峰强度不同。

应当根据波尔兹曼分布做权重平均,相关知识参考
使用Multiwfn绘制构象权重平均的光谱
http://sobereva.com/383http://bbs.keinsci.com/thread-6186-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

122

帖子

0

威望

2336

eV
积分
2458

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2019-10-23 21:34:21 | 只看该作者 Only view this author
sobereva 发表于 2019-10-23 01:54
那是展宽,不叫拟合

需要对电离态和中性状态都做振动分析,参考下文。和文中的区别是要把下文中的激发态 ...

谢谢老师,有个地方请教一下。
博文里苯甲醚的优化S1激发态和基态用的多重度都是1呢?光子打掉一个电子不应该是2吗?

~3LDE7_}RR0YPA0AJ~P4IRI.png (18.7 KB, 下载次数 Times of downloads: 19)

~3LDE7_}RR0YPA0AJ~P4IRI.png

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

4#
发表于 Post on 2019-10-24 07:09:54 | 只看该作者 Only view this author
my_aman 发表于 2019-10-23 21:34
谢谢老师,有个地方请教一下。
博文里苯甲醚的优化S1激发态和基态用的多重度都是1呢?光子 ...

S1是被激发到第一单重态激发态,不是被电离
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

122

帖子

0

威望

2336

eV
积分
2458

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2019-10-24 11:15:02 | 只看该作者 Only view this author
sobereva 发表于 2019-10-24 07:09
S1是被激发到第一单重态激发态,不是被电离

感谢耐心解答,还有些疑问。我现在重复文献上CN自由基的电离
问题1:是不是需要按着博文上操作,把第一步的优化S1激发态改成CN电离后的自旋是1或者3,第二步优化基态CN自由基自旋仍然是2。问题2:CN自由基电离后自旋有1或者3,这个时候要按着上面操作两次?得到的数据要怎么弄呢?

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

6#
发表于 Post on 2019-10-26 04:39:49 | 只看该作者 Only view this author
my_aman 发表于 2019-10-24 11:15
感谢耐心解答,还有些疑问。我现在重复文献上CN自由基的电离
问题1:是不是需要按着博文上操作,把第一 ...

1 是
2 电子电离掉之后,体系可以处于阳离子基态,也可能处于阳离子激发态,两种情况在实验谱上应当都能观测到(我不很了解这方面实验,具体不清楚),想考察哪个情况就算哪种
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

122

帖子

0

威望

2336

eV
积分
2458

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2019-12-6 07:16:42 | 只看该作者 Only view this author
本帖最后由 my_aman 于 2019-12-6 11:02 编辑
sobereva 发表于 2019-10-26 04:39
1 是
2 电子电离掉之后,体系可以处于阳离子基态,也可能处于阳离子激发态,两种情况在实验谱上应当都能 ...

因为我实验上的光电子谱已经出来。现在急需重复文献的工作,进行下面的研究。希望各位老师同学热心指导,谢谢。
我按着博客的操作,计算由CN自由基2重态基态电离到单重态,的电离能操作如下:第一步操作:

第二步操作:

结果在第二步报错l718
请问,按着博客的操作有问题吗?若没错,报的718该如何解决呢?谢谢老师


$@58JM{}9TL5C}5_W21M7HM.png (14.02 KB, 下载次数 Times of downloads: 12)

$@58JM{}9TL5C}5_W21M7HM.png

J}G[6KT`BDU0KU%]7DM1{ZV.png (15.31 KB, 下载次数 Times of downloads: 14)

J}G[6KT`BDU0KU%]7DM1{ZV.png

2.log

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

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

8#
发表于 Post on 2019-12-8 19:13:03 | 只看该作者 Only view this author
my_aman 发表于 2019-12-6 07:16
因为我实验上的光电子谱已经出来。现在急需重复文献的工作,进行下面的研究。希望各位老师同学热心指导, ...

换个泛函或者基组再试
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

122

帖子

0

威望

2336

eV
积分
2458

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2019-12-10 15:00:16 | 只看该作者 Only view this author
sobereva 发表于 2019-12-8 19:13
换个泛函或者基组再试

谢谢社长,问题已解决。论坛的好友有需要这方面的可以联系我,一起学习。

再请教个问题“normal mode vectors for each electronic state“里的normal mode vectors信息,我是加个freq关键词就能算出来吗?

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

10#
发表于 Post on 2019-12-15 10:50:31 | 只看该作者 Only view this author
my_aman 发表于 2019-12-10 15:00
谢谢社长,问题已解决。论坛的好友有需要这方面的可以联系我,一起学习。

再请教个问题“normal mode  ...

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

122

帖子

0

威望

2336

eV
积分
2458

Level 5 (御坂)

11#
 楼主 Author| 发表于 Post on 2019-12-16 10:54:01 | 只看该作者 Only view this author

请问社长我用gaussian结合ezSpectrum程序算FC因子,但是ezSpectrum程序需要这两个值:maximum number of the vibrational excitations in the initial and target electronic states (the total number of quanta in all normal modes),这两个值我怎么得到呢?没有明白这句话的意思

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

12#
发表于 Post on 2019-12-16 11:52:44 | 只看该作者 Only view this author
my_aman 发表于 2019-12-16 10:54
请问社长我用gaussian结合ezSpectrum程序算FC因子,但是ezSpectrum程序需要这两个值:maximum number of  ...

我没用过,应该是用户自己设的考虑振动态的上限值,根据需要设置。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1

帖子

0

威望

9

eV
积分
10

Level 1 能力者

13#
发表于 Post on 2025-1-6 13:21:29 | 只看该作者 Only view this author
my_aman 发表于 2019-12-10 15:00
谢谢社长,问题已解决。论坛的好友有需要这方面的可以联系我,一起学习。

再请教个问题“normal mode  ...

您好,我利用G09计算TPES时,同样遇到了l718报错的问题,特向您请教如何解决这个问题呢?百忙之中打扰您,期待您的回复!

本版积分规则 Credits rule

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

GMT+8, 2026-2-18 18:47 , Processed in 0.192820 second(s), 29 queries , Gzip On.

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