计算化学公社

标题: 求助：得到残缺波函数再优化计算频率后依旧不破缺 [打印本页]

作者
Author: 不干正事 时间: 2022-3-11 09:01
标题: 求助：得到残缺波函数再优化计算频率后依旧不破缺
大家好，计算一个双自由基体系，打算通过两步计算，第一步的关键词是 #p 6-311g(d) nosymm upbe1pbe guess=mix stable=opt；计算的结果的确显示不破缺了，自旋密度也处处不为零，然后就对其优化和计算频率，第二布的关键词是 #p 6-311g(d) nosymm upbe1pbe opt=tight freq guess=read；但这一步计算完后，结果显示还是不破缺，自旋密度也是处处为0，不知道哪里出现了问题。附上了两步的输入与输出文件。

作者
Author: sobereva 时间: 2022-3-12 04:58
用词要严谨，对称破缺不叫“残缺”
基于优化后的结构，结合guess=read关键词从优化任务的chk里读取波函数当初猜，做个stable任务，判断当前结构下稳定波函数是否理应是对称破缺的。如果当前波函数不稳定，基于当前结构和stable=opt搞出来的稳定波函数进一步优化

作者
Author: 不干正事 时间: 2022-3-12 08:43

sobereva 发表于 2022-3-12 04:58
用词要严谨，对称破缺不叫“残缺”
基于优化后的结构，结合guess=read关键词从优化任务的chk里读取波函数 ...

好的，谢谢卢老师，我再看看

作者
Author: 不干正事 时间: 2022-3-13 14:53

sobereva 发表于 2022-3-12 04:58
用词要严谨，对称破缺不叫“残缺”
基于优化后的结构，结合guess=read关键词从优化任务的chk里读取波函数 ...

卢老师我又尝试了，不加opt和freq关键词，采用“nosymm guess=mix stable =opt 的关键词时能有对称破缺，而且也显示波函数稳定，然而当加上opt，freq和guess=read读取波函数后的计算结果后发现，又无法产生对称破缺了

作者
Author: sobereva 时间: 2022-3-14 02:06

不干正事发表于 2022-3-13 14:53
卢老师我又尝试了，不加opt和freq关键词，采用“nosymm guess=mix stable =opt 的关键词时能有对称破缺， ...

既可能极小点结构下本来稳定波函数就不是对称破缺的，也可能优化过程中又恰好变成闭壳层状态了，后者还得重新产生对称破缺波函数继续优化

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