计算化学公社

标题: 气相下优化单分子激发态时，结构不收敛，且氢原子解离 [打印本页]

作者
Author: 爱学习的傻子 时间: 2022-2-24 17:11
标题: 气相下优化单分子激发态时，结构不收敛，且氢原子解离
本帖最后由爱学习的傻子于 2022-2-24 19:57 编辑

大家好，我在s0的基础上优化硼酸的s1，但是结构不能收敛，且氢原子还会跑散，这是为什么？是泛函或者基组选则的问题吗？

(, 下载次数 Times of downloads: 38)

设置的参数：#p opt freq td=(nstates=5) wb97xd/def2svp

我手动稍微改变分子的初始结构，s1会优化为下图
(, 下载次数 Times of downloads: 38)
虽然没有虚频，但结构相比s0变化很大，是不是也不对？

还请各位老师指点，谢谢！！！

作者
Author: wzkchem5 时间: 2022-2-24 17:36
图没有贴成功
首先，并不是所有分子的S1优化都可以收敛，而且并不是所有分子的S1优化都不会把氢原子跑散。你先说一下你为什么认为这个优化应该收敛，为什么你认为氢原子不应该跑散（比如有没有文献优化这个体系收敛了而且没有跑散），然后我们才能看问题出在了哪里

作者
Author: 爱学习的傻子 时间: 2022-2-24 20:05

wzkchem5 发表于 2022-2-24 17:36
图没有贴成功
首先，并不是所有分子的S1优化都可以收敛，而且并不是所有分子的S1优化都不会把氢原子跑散。 ...

感谢您的回复，我把图重新编辑了一下。我们在实验上测得了偏硼酸粉末的荧光和磷光发射，所以想在量子计算上对发光性质进行进一步研究。之前确实没有查到这个结构的理论计算实验，但有一篇B(OH)3结构的硼酸理论计算，我在复现的时候也遇到了类似问题，一个氢原子跑开了，所以不知道问题出在哪里。

作者
Author: wzkchem5 时间: 2022-2-24 20:48

爱学习的傻子发表于 2022-2-24 13:05
感谢您的回复，我把图重新编辑了一下。我们在实验上测得了偏硼酸粉末的荧光和磷光发射，所以想在量子计算 ...

也就是说文章里面氢原子没跑开，但是你的计算跑开了是吗？
把那篇文章贴上来，我们看一下？

作者
Author: 爱学习的傻子 时间: 2022-2-24 21:00

wzkchem5 发表于 2022-2-24 20:48
也就是说文章里面氢原子没跑开，但是你的计算跑开了是吗？
把那篇文章贴上来，我们看一下？

(, 下载次数 Times of downloads: 6) 麻烦您看一下谢谢

作者
Author: wzkchem5 时间: 2022-2-24 21:56

爱学习的傻子发表于 2022-2-24 14:00
麻烦您看一下谢谢

这篇文章好像没涉及优化单独的硼酸分子的S1态结构？

作者
Author: 爱学习的傻子 时间: 2022-2-24 22:57

wzkchem5 发表于 2022-2-24 21:56
这篇文章好像没涉及优化单独的硼酸分子的S1态结构？

(, 下载次数 Times of downloads: 42)
抱歉没有贴出 (, 下载次数 Times of downloads: 2) 。
他在支撑里说一句，只是没有贴s1的结构，但是他的意思也是应该说明s1结构收敛了吧。我不太清楚这种激发态结构优化，氢原子为什么会跑散。

作者
Author: wzkchem5 时间: 2022-2-24 23:06

爱学习的傻子发表于 2022-2-24 15:57
抱歉没有贴出。
他在支撑里说一句，只是没有贴s1的结构，但是他的意思也是应该说明s1结构收敛了吧。我 ...

你用文献的理论级别优化一下看看，你的泛函和文献不一样，结果不一样也正常

作者
Author: 爱学习的傻子 时间: 2022-2-24 23:57

wzkchem5 发表于 2022-2-24 23:06
你用文献的理论级别优化一下看看，你的泛函和文献不一样，结果不一样也正常

(, 下载次数 Times of downloads: 4) 这个我算过了不行

作者
Author: wzkchem5 时间: 2022-2-25 01:12

爱学习的傻子发表于 2022-2-24 16:57
这个我算过了不行

基组也换成和文献一样的话，也重复不出来文献的结果吗？

作者
Author: 爱学习的傻子 时间: 2022-2-25 09:30

wzkchem5 发表于 2022-2-25 01:12
基组也换成和文献一样的话，也重复不出来文献的结果吗？

是的

作者
Author: wzkchem5 时间: 2022-2-25 16:01

爱学习的傻子发表于 2022-2-25 02:30
是的

那要不发邮件问一下作者吧，可能作者用的初始结构和你的不一样，所以优化到了激发态势能面的不同地方

作者
Author: 爱学习的傻子 时间: 2022-2-25 17:14

wzkchem5 发表于 2022-2-25 16:01
那要不发邮件问一下作者吧，可能作者用的初始结构和你的不一样，所以优化到了激发态势能面的不同地方

好的谢谢您

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