计算化学公社

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

[新手求助] 求助:IRC跑完后波函数不稳定该如何解决

[复制链接 Copy URL]

138

帖子

0

威望

662

eV
积分
800

Level 4 (黑子)

本帖最后由 Xx_DerEnTh 于 2025-1-22 05:56 编辑

请教各位老师,我计算带锕系元素的体系时,在保证过渡态波函数稳定的情况下通过guess=read读取过渡态chk文件跑了IRC,跑完后我验证IRC最后一个点的波函数稳定性时,显示不稳定。
问:1.这种情况下这个不稳定的IRC得到的反应物/产物的近似结构是否可以代表该过渡态对应的反应物/产物结构?
      2.如果问题1可以的话,后续是否应该先用stable=opt让IRC的反应物/产物结构波函数稳定后,再进行几何优化?如果问题1不可以,这条IRC该如何处理来解决此波函数不稳定的问题?

138

帖子

0

威望

662

eV
积分
800

Level 4 (黑子)

18#
 楼主 Author| 发表于 Post on 2025-3-5 12:45:28 | 只看该作者 Only view this author
sobereva 发表于 2025-3-5 12:44
诸如按照下文查看自旋密度和自旋布居
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sob ...

感谢sob老师回复!这就去学习

6万

帖子

99

威望

6万

eV
积分
125504

管理员

公社社长

17#
发表于 Post on 2025-3-5 12:44:10 | 只看该作者 Only view this author
Xx_DerEnTh 发表于 2025-3-5 12:21
sob老师您好!自身理论基础实在薄弱,请教一下稳定和不稳定波函数的特征具体是什么对比的?从哪里可以得 ...

诸如按照下文查看自旋密度和自旋布居
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353http://bbs.keinsci.com/thread-4422-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

138

帖子

0

威望

662

eV
积分
800

Level 4 (黑子)

16#
 楼主 Author| 发表于 Post on 2025-3-5 12:21:31 | 只看该作者 Only view this author
sobereva 发表于 2025-3-4 19:51
对比稳定和不稳定波函数的特征,弄清楚差异在哪,比如单电子分布情况等

sob老师您好!自身理论基础实在薄弱,请教一下稳定和不稳定波函数的特征具体是什么对比的?从哪里可以得知数据?请求sob老师指路,我去学习一下。
还想请教sob老师:这是stable=opt优化前后的结果,从这个数据中能够看出什么吗?
Eigenvector   1:  6.021-A    Eigenvalue=-0.0070214  <S**2>=8.814
     21A -> 28A       -0.17945
     22A -> 28A        0.36725
     23A -> 29A        0.28867
     25A -> 28A       -0.15995
     27A -> 28A        0.84143
The wavefunction has an internal instability.

Eigenvector   1:  6.016-A    Eigenvalue= 0.0000881  <S**2>=8.797
     22A -> 29A       -0.22311
     24A -> 28A        0.53952
     25A -> 29A        0.63595
     26A -> 28A       -0.12260
     27A -> 28A       -0.46858
The wavefunction is stable under the perturbations considered.
先谢谢sob老师了!

138

帖子

0

威望

662

eV
积分
800

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2025-3-4 21:24:18 | 只看该作者 Only view this author
sobereva 发表于 2025-3-4 19:51
对比稳定和不稳定波函数的特征,弄清楚差异在哪,比如单电子分布情况等

好的,感谢sob老师!

6万

帖子

99

威望

6万

eV
积分
125504

管理员

公社社长

14#
发表于 Post on 2025-3-4 19:51:03 | 只看该作者 Only view this author
Xx_DerEnTh 发表于 2025-3-4 18:12
感谢sob老师回复,想请教sob老师,这种情况该如何进一步分析?目前进度停滞不前了,已想不到其他解决办法

对比稳定和不稳定波函数的特征,弄清楚差异在哪,比如单电子分布情况等
北京科音自然科学研究中心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

138

帖子

0

威望

662

eV
积分
800

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2025-3-4 18:12:46 | 只看该作者 Only view this author
sobereva 发表于 2025-3-4 17:03
可能
不能。稳定和不稳定波函数对应的势能面的极小点坐标一般是不同的

感谢sob老师回复,想请教sob老师,这种情况该如何进一步分析?目前进度停滞不前了,已想不到其他解决办法

6万

帖子

99

威望

6万

eV
积分
125504

管理员

公社社长

12#
发表于 Post on 2025-3-4 17:03:33 | 只看该作者 Only view this author
Xx_DerEnTh 发表于 2025-3-3 12:53
sob老师您好!请教一下:如果通过IRC曲线判断出两个态的势能面完全孤立,明显不会相交,但是其IRC的最后 ...

可能
不能。稳定和不稳定波函数对应的势能面的极小点坐标一般是不同的
北京科音自然科学研究中心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

138

帖子

0

威望

662

eV
积分
800

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2025-3-3 12:53:37 | 只看该作者 Only view this author
本帖最后由 Xx_DerEnTh 于 2025-3-4 14:01 编辑
sobereva 发表于 2025-2-19 08:51
诸如当前情况,两个自旋态势能面交叉。倘若你对红线的势能面跑的IRC,此势能面的极小点结构下的稳定 ...

sob老师您好!请教一下:如果通过IRC曲线判断出两个态的势能面完全孤立,明显不会相交,但是其IRC的最后一个点波函数仍然显示不稳定,这种情况还可能是因为什么原因造成的?有没有可能过渡态的波函数稳定但跑IRC中途导致波函数不稳定,从而导致IRC最后一个点波函数不稳定?这种情况能否直接对IRC端点不稳定的结构进行stable=opt,然后直接使用优化出来的能量?谢谢sob老师!

138

帖子

0

威望

662

eV
积分
800

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2025-2-19 11:50:35 | 只看该作者 Only view this author
sobereva 发表于 2025-2-19 08:51
诸如当前情况,两个自旋态势能面交叉。倘若你对红线的势能面跑的IRC,此势能面的极小点结构下的稳定 ...

明白了,感谢sob老师!

6万

帖子

99

威望

6万

eV
积分
125504

管理员

公社社长

9#
发表于 Post on 2025-2-19 08:51:49 | 只看该作者 Only view this author
Xx_DerEnTh 发表于 2025-2-18 16:59
Sob老师您好!请教一下:波函数不稳定是说在当前结构当前计算级别下存在比此时能量更低的结构。为什么IRC ...



诸如当前情况,两个自旋态势能面交叉。倘若你对红线的势能面跑的IRC,此势能面的极小点结构下的稳定波函数显然不是红线对应的态,而是黑线对应的态
北京科音自然科学研究中心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

138

帖子

0

威望

662

eV
积分
800

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2025-2-18 16:59:19 | 只看该作者 Only view this author
sobereva 发表于 2025-1-22 13:46
参考北京科音中级量子化学培训班(http://www.keinsci.com/KBQC)以下幻灯片以及引用的文献

Sob老师您好!请教一下:波函数不稳定是说在当前结构当前计算级别下存在比此时能量更低的结构。为什么IRC端点处波函数不稳定会联想到两态反应呢?是因为此时波函数不稳定代表着当前结构下有另一个不同的自旋态能量更低吗?如果IRC一端波函数稳定,另一端不稳定,是说明稳定的那一端不存在自旋翻转,而不稳定的那一端存在自旋翻转吗?

138

帖子

0

威望

662

eV
积分
800

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2025-1-22 14:15:39 | 只看该作者 Only view this author
sobereva 发表于 2025-1-22 13:46
参考北京科音中级量子化学培训班(http://www.keinsci.com/KBQC)以下幻灯片以及引用的文献

好的,感谢sob老师

6万

帖子

99

威望

6万

eV
积分
125504

管理员

公社社长

6#
发表于 Post on 2025-1-22 13:46:25 | 只看该作者 Only view this author
Xx_DerEnTh 发表于 2025-1-22 12:43
确实存在化合价变化,想请教老师,如何判断是否存在two-state reactivity的情况?这方面接触较少,不是特 ...

参考北京科音中级量子化学培训班(http://www.keinsci.com/KBQC)以下幻灯片以及引用的文献





这种情况如果要把反应过程了解全面透彻,应当在不同透热态势能面下都找驻点、跑IRC,并且用sobMECP(http://sobereva.com/286)之类程序把MECP搜索出来。
北京科音自然科学研究中心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

138

帖子

0

威望

662

eV
积分
800

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-1-22 12:43:38 | 只看该作者 Only view this author
AndrewX 发表于 2025-1-22 10:27
该反应是否存在锕系元素化合价变化?可能存在two-state reactivity的情况。

确实存在化合价变化,想请教老师,如何判断是否存在two-state reactivity的情况?这方面接触较少,不是特别懂

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

GMT+8, 2026-3-13 12:20 , Processed in 0.179404 second(s), 24 queries , Gzip On.

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