气相下优化单分子激发态时，结构不收敛，且氢原子解离

爱学习的傻子 · 发表于 Post on 2022-2-24 17:11:40

本帖最后由爱学习的傻子于 2022-2-24 19:57 编辑

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

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

我手动稍微改变分子的初始结构，s1会优化为下图

虽然没有虚频，但结构相比s0变化很大，是不是也不对？

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

wzkchem5 · 发表于 Post on 2022-2-24 17:36:34

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

爱学习的傻子 · 发表于 Post on 2022-2-24 20:05:52

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

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

wzkchem5 · 发表于 Post on 2022-2-24 20:48:14

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

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

爱学习的傻子 · 发表于 Post on 2022-2-24 21:00:21

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

文献.pdf (3.18 MB, 下载次数 Times of downloads: 6) 麻烦您看一下谢谢

wzkchem5 · 发表于 Post on 2022-2-24 21:56:45

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

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

爱学习的傻子 · 发表于 Post on 2022-2-24 22:57:38

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

抱歉没有贴出

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

wzkchem5 · 发表于 Post on 2022-2-24 23:06:42

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

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

爱学习的傻子 · 发表于 Post on 2022-2-24 23:57:14

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

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

wzkchem5 · 发表于 Post on 2022-2-25 01:12:39

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

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

爱学习的傻子 · 发表于 Post on 2022-2-25 09:30:30

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

是的

wzkchem5 · 发表于 Post on 2022-2-25 16:01:34

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

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

爱学习的傻子 · 发表于 Post on 2022-2-25 17:14:01

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

好的谢谢您

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 气相下优化单分子激发态时，结构不收敛，且氢原子解离