计算化学公社

标题: 遇到 “NBO计算让原本稳定的波函数变得不稳定” 的情况 [打印本页]

作者
Author: Finale 时间: 2023-3-1 21:26
标题: 遇到 “NBO计算让原本稳定的波函数变得不稳定” 的情况
本帖最后由 Finale 于 2023-3-1 21:27 编辑

背景：最近在计算一个双稀土金属（Gd）配合物的体系，两个Gd 都是 +3 价，使用PBE0-D3（BJ）泛函，配体部分使用 def2tzvp，金属使用SDD小核赝势和对应的SDD基组。

通过片段波函数的方法，成功找到稳定的自旋极化波函数，并且自旋密度符合预期，但是为了进行 AdNDP 分析，需要进行 NAO 计算，使用的是稳定的自旋极化波函数，但是计算总是失败，显示
" The wavefunction has an RHF -> UHF instability. " 以及 ” Failed to find a lower point in StbLin: “ 这两个提示。使用的是内嵌在 Gaussian 中的 NBO 7.0 程序。

查看了输出文件，发现前面的 SCF 成功收敛，并且能量也和自旋极化态能量一致，NAO 计算也成功完成，但是后面在分析波函数稳定的时候出现错误，如图是输入文件的字段

这里有两个疑问想请教各位老师：
1、在找到稳定波函数后，再进行后续其他计算（比如 NAO 计算）的时候是否还需要检查波函数？
因为这种计算不涉及到波函数的优化（虽然经过了一两步 SCF 迭代），所以是否可以去掉 stable=opt 呢？

2、这种情况大家有遇到过吗，如果在得到稳定的波函数后，如果进行其他计算（使用初猜），波函数竟然变得不稳定，那么前一步 stable=opt 是否可信呢？
我目前去掉了 NAO 计算，单纯做一个单点（使用 Stable =opt）看看是否出现相同情况（结果来没出来）

谢谢各位老师！

作者
Author: hebrewsnabla 时间: 2023-3-1 22:11
guess=read时不要用vshift

作者
Author: sobereva 时间: 2023-3-2 06:18
先做个单点任务照常得到稳定波函数。之后再单独做pop=nboread任务，届时直接guess=read，不用加其它的

作者
Author: Finale 时间: 2023-3-2 15:29

hebrewsnabla 发表于 2023-3-1 22:11
guess=read时不要用vshift

谢谢您的提醒！

作者
Author: Finale 时间: 2023-3-2 15:29

sobereva 发表于 2023-3-2 06:18
先做个单点任务照常得到稳定波函数。之后再单独做pop=nboread任务，届时直接guess=read，不用加其它的

谢谢社长的建议！

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