关于计算开壳层单重态的问题求助

Jameslzy

冰释之川 发表于 2018-5-19 11:50
必须保留chk文件，然后在新任务中guess=read 读取之前的开壳层单重态的波函数。如果光给出坐标直接算单点 ...

恩恩好的，谢谢老师解答，那我还是注意保留chk文件好了

Jameslzy

Accelerator 发表于 2018-5-19 11:14
每一个单点都写stable=opt啊

恩恩好的，谢谢老师回复，我也是好奇地问一下，可能还是读取chk文件比较好

FANYI

老师好，我现在也有一个关于开壳层单重态计算的问题。我的反应基本过程是金属络合物催化剂和反应物A氧化加成，再与反应物B氧化加成，然后整个体系还原消除得到产物。 我最初只考虑的是催化剂的闭壳层单重态的催化反应历程，现在突然发现其实催化剂在三重态的能量比单重态的低，于是在重新做三重态位能面上的反应； 有师兄建议说，再做下开壳层单重态的路径，这个路径上的优化直接去读取三重态的chk 文件就可以。 现在，我有个疑惑：
对于这个催化剂是开壳层单重态的反应路径，除了催化剂以外，另外两个反应物A、B 的单独能量优化频率分析 以及 反应过程中经历的中间体和过渡态，我的输入文件中都得设置 关键词 ub3lyp 吗？ 期待老师的指点迷津，谢谢老师！

sobereva

FANYI 发表于 2019-5-31 09:27
老师好，我现在也有一个关于开壳层单重态计算的问题。我的反应基本过程是金属络合物催化剂和反应物A氧化加 ...

是否计算对称破缺单重态跟写不写U毫无必然关系，仔细看
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

FANYI

sobereva 发表于 2019-6-1 09:21
是否计算对称破缺单重态跟写不写U毫无必然关系，仔细看
谈谈片段组合波函数与自旋极化单重态
http://so ...

非常感谢老师的指点~

lycheeho

sob老师您好！我在计算稠环芳烃双自由基体系时，经常是先按b3lyp/6-31g* opt freq算闭壳层单重态，ub3lyp/6-31g* opt freq算开壳层单重态和三重态，然后比较三者能量，往往闭壳层单重态和开壳层单重态能量相等，对闭壳层单重态stable=opt后显示有向开壳层的趋势，但得到的能量还是跟原来一样，请问这样正常吗？

Jameslzy

lycheeho 发表于 2023-4-20 22:15
sob老师您好！我在计算稠环芳烃双自由基体系时，经常是先按b3lyp/6-31g* opt freq算闭壳层单重态，ub3lyp/6 ...

直接在我这个帖子里发帖的话我估计Sob老师看不到呀，可以去简单量化问题答疑帖里面问Sob老师呀。关于你这个问题，我感觉应该是没有顺利找到开壳层单重态，可以直接在输出文件里面搜Mulliken，看看马利肯分析里面不同原子上的电子自旋是有数值的还是全为零，如果全为零的话，就证明没找到开壳层单重态，就可以看看Sob老师在思想家公社门口帖子里介绍的几种方法，比如用guess=mix和stable=opt之类的方法再找一下。用stable=opt之后要读得出的chk文件，再去优化才能得出开壳层单重态。

lycheeho

Jameslzy 发表于 2023-6-8 21:10
直接在我这个帖子里发帖的话我估计Sob老师看不到呀，可以去简单量化问题答疑帖里面问Sob老师呀。关于你这 ...

谢谢~

火星哥

冰释之川 发表于 2018-5-19 11:50
必须保留chk文件，然后在新任务中guess=read 读取之前的开壳层单重态的波函数。如果光给出坐标直接算单点 ...

冰哥，这篇帖子是不是说明对称破缺的态的稳定点包含过渡态的单点能需要校正啊？还是说使用guess=mix stable=opt算出来的就可以用？

火星哥

您好，可以分享一下原文献吗