计算化学公社

标题: 求助：IRC跑完后波函数不稳定该如何解决 [打印本页]

作者
Author: Xx_DerEnTh 时间: 2025-1-22 05:54
标题: 求助：IRC跑完后波函数不稳定该如何解决
本帖最后由 Xx_DerEnTh 于 2025-1-22 05:56 编辑

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

作者
Author: AndrewX 时间: 2025-1-22 09:05
过渡态的结构是否用stable检查过，波函数是否稳定？

作者
Author: Xx_DerEnTh 时间: 2025-1-22 09:45

AndrewX 发表于 2025-1-22 09:05
过渡态的结构是否用stable检查过，波函数是否稳定？

过渡态波函数是稳定的

作者
Author: AndrewX 时间: 2025-1-22 10:27

Xx_DerEnTh 发表于 2025-1-22 09:45
过渡态波函数是稳定的

该反应是否存在锕系元素化合价变化？可能存在two-state reactivity的情况。

作者
Author: Xx_DerEnTh 时间: 2025-1-22 12:43

AndrewX 发表于 2025-1-22 10:27
该反应是否存在锕系元素化合价变化？可能存在two-state reactivity的情况。

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

作者
Author: sobereva 时间: 2025-1-22 13:46

Xx_DerEnTh 发表于 2025-1-22 12:43
确实存在化合价变化，想请教老师，如何判断是否存在two-state reactivity的情况？这方面接触较少，不是特 ...

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

(, 下载次数 Times of downloads: 7)
(, 下载次数 Times of downloads: 8)

这种情况如果要把反应过程了解全面透彻，应当在不同透热态势能面下都找驻点、跑IRC，并且用sobMECP（http://sobereva.com/286）之类程序把MECP搜索出来。

作者
Author: Xx_DerEnTh 时间: 2025-1-22 14:15

sobereva 发表于 2025-1-22 13:46
参考北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）以下幻灯片以及引用的文献

好的，感谢sob老师

作者
Author: Xx_DerEnTh 时间: 2025-2-18 16:59

sobereva 发表于 2025-1-22 13:46
参考北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）以下幻灯片以及引用的文献

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

作者
Author: sobereva 时间: 2025-2-19 08:51

Xx_DerEnTh 发表于 2025-2-18 16:59
Sob老师您好！请教一下:波函数不稳定是说在当前结构当前计算级别下存在比此时能量更低的结构。为什么IRC ...

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

诸如当前情况，两个自旋态势能面交叉。倘若你对红线的势能面跑的IRC，此势能面的极小点结构下的稳定波函数显然不是红线对应的态，而是黑线对应的态

作者
Author: Xx_DerEnTh 时间: 2025-2-19 11:50

sobereva 发表于 2025-2-19 08:51
诸如当前情况，两个自旋态势能面交叉。倘若你对红线的势能面跑的IRC，此势能面的极小点结构下的稳定 ...

明白了，感谢sob老师！

作者
Author: Xx_DerEnTh 时间: 2025-3-3 12:53
本帖最后由 Xx_DerEnTh 于 2025-3-4 14:01 编辑

sobereva 发表于 2025-2-19 08:51
诸如当前情况，两个自旋态势能面交叉。倘若你对红线的势能面跑的IRC，此势能面的极小点结构下的稳定 ...

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

作者
Author: sobereva 时间: 2025-3-4 17:03

Xx_DerEnTh 发表于 2025-3-3 12:53
sob老师您好！请教一下：如果通过IRC曲线判断出两个态的势能面完全孤立，明显不会相交，但是其IRC的最后 ...

可能
不能。稳定和不稳定波函数对应的势能面的极小点坐标一般是不同的

作者
Author: Xx_DerEnTh 时间: 2025-3-4 18:12

sobereva 发表于 2025-3-4 17:03
可能
不能。稳定和不稳定波函数对应的势能面的极小点坐标一般是不同的

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

作者
Author: sobereva 时间: 2025-3-4 19:51

Xx_DerEnTh 发表于 2025-3-4 18:12
感谢sob老师回复，想请教sob老师，这种情况该如何进一步分析？目前进度停滞不前了，已想不到其他解决办法

对比稳定和不稳定波函数的特征，弄清楚差异在哪，比如单电子分布情况等

作者
Author: Xx_DerEnTh 时间: 2025-3-4 21:24

sobereva 发表于 2025-3-4 19:51
对比稳定和不稳定波函数的特征，弄清楚差异在哪，比如单电子分布情况等

好的，感谢sob老师！

作者
Author: Xx_DerEnTh 时间: 2025-3-5 12:21

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老师了！

作者
Author: sobereva 时间: 2025-3-5 12:44

Xx_DerEnTh 发表于 2025-3-5 12:21
sob老师您好！自身理论基础实在薄弱，请教一下稳定和不稳定波函数的特征具体是什么对比的？从哪里可以得 ...

诸如按照下文查看自旋密度和自旋布居
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353（http://bbs.keinsci.com/thread-4422-1-1.html）

作者
Author: Xx_DerEnTh 时间: 2025-3-5 12:45

sobereva 发表于 2025-3-5 12:44
诸如按照下文查看自旋密度和自旋布居
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sob ...

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

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