计算化学公社

标题: 关于计算threshold photoelectron (TPE) spectrum的相关问题 [打印本页]

作者
Author: my_aman 时间: 2019-10-22 15:07
标题: 关于计算threshold photoelectron (TPE) spectrum的相关问题
本帖最后由 my_aman 于 2019-10-22 15:32 编辑

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

作者
Author: sobereva 时间: 2019-10-23 01:54
那是展宽，不叫拟合

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

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

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

应当根据波尔兹曼分布做权重平均，相关知识参考
使用Multiwfn绘制构象权重平均的光谱
http://sobereva.com/383（http://bbs.keinsci.com/thread-6186-1-1.html）

作者
Author: my_aman 时间: 2019-10-23 21:34

sobereva 发表于 2019-10-23 01:54
那是展宽，不叫拟合

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

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

作者
Author: sobereva 时间: 2019-10-24 07:09

my_aman 发表于 2019-10-23 21:34
谢谢老师，有个地方请教一下。
博文里苯甲醚的优化S1激发态和基态用的多重度都是1呢？光子 ...

S1是被激发到第一单重态激发态，不是被电离

作者
Author: my_aman 时间: 2019-10-24 11:15

sobereva 发表于 2019-10-24 07:09
S1是被激发到第一单重态激发态，不是被电离

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

作者
Author: sobereva 时间: 2019-10-26 04:39

my_aman 发表于 2019-10-24 11:15
感谢耐心解答，还有些疑问。我现在重复文献上CN自由基的电离
问题1：是不是需要按着博文上操作，把第一 ...

1 是
2 电子电离掉之后，体系可以处于阳离子基态，也可能处于阳离子激发态，两种情况在实验谱上应当都能观测到（我不很了解这方面实验，具体不清楚），想考察哪个情况就算哪种

作者
Author: my_aman 时间: 2019-12-6 07:16
本帖最后由 my_aman 于 2019-12-6 11:02 编辑

sobereva 发表于 2019-10-26 04:39
1 是
2 电子电离掉之后，体系可以处于阳离子基态，也可能处于阳离子激发态，两种情况在实验谱上应当都能 ...

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

第二步操作：[attach]22805[/attach]

结果在第二步报错l718
请问，按着博客的操作有问题吗？若没错，报的718该如何解决呢？谢谢老师
(, 下载次数 Times of downloads: 28)

作者
Author: sobereva 时间: 2019-12-8 19:13

my_aman 发表于 2019-12-6 07:16
因为我实验上的光电子谱已经出来。现在急需重复文献的工作，进行下面的研究。希望各位老师同学热心指导， ...

换个泛函或者基组再试

作者
Author: my_aman 时间: 2019-12-10 15:00

sobereva 发表于 2019-12-8 19:13
换个泛函或者基组再试

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

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

作者
Author: sobereva 时间: 2019-12-15 10:50

my_aman 发表于 2019-12-10 15:00
谢谢社长，问题已解决。论坛的好友有需要这方面的可以联系我，一起学习。

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

是

作者
Author: my_aman 时间: 2019-12-16 10:54

sobereva 发表于 2019-12-15 10:50
是

请问社长我用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)，这两个值我怎么得到呢？没有明白这句话的意思

作者
Author: sobereva 时间: 2019-12-16 11:52

my_aman 发表于 2019-12-16 10:54
请问社长我用gaussian结合ezSpectrum程序算FC因子，但是ezSpectrum程序需要这两个值：maximum number of ...

我没用过，应该是用户自己设的考虑振动态的上限值，根据需要设置。

作者
Author: ZoeWu 时间: 2025-1-6 13:21

my_aman 发表于 2019-12-10 15:00
谢谢社长，问题已解决。论坛的好友有需要这方面的可以联系我，一起学习。

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

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

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