求助：丁烷双自由基结构优化求助

斯卡姆

老师，您好，我今天算了一下基础班中丁烷双自由基的对称破缺函数，先是优化丁烷双自由基的结构，采用的是UB3LYP/6-31G* opt，但是优化到最后，结构从中间断开成为了2个片段，软件用的是G16 A.03版本，请问这是哪里出了问题？

zjxitcc

问题就是你的“优化丁烷双自由基的结构”。首先你得确定这个练习是让你学会算个单点，还是有让你做结构优化，不要自己添加多余的步骤。这双自由基你还优化，那就是要么断掉要么成键，无论哪种都不是双自由基。

另外，真正的双自由基结构也不是你这么优化的，缺乏UHF/UDFT的基本常识，看sob老师博文中介绍的自旋极化单重态《谈谈片段组合波函数与自旋极化单重态》http://sobereva.com/82

斯卡姆

zjxitcc 发表于 2020-9-9 13:06
问题就是你的“优化丁烷双自由基的结构”。首先你得确定这个练习是让你学会算个单点，还是有让你做结构优化 ...

啊，谢谢ni老师的回复。我的目的是优化丁烷双自由基的结构，获得含有对称破缺的函数。啊，这个是在说如何得到对称破缺的初猜波函数。

zjxitcc

斯卡姆 发表于 2020-9-9 20:01
啊，谢谢ni老师的回复。我的目的是优化丁烷双自由基的结构，获得含有对称破缺的函数。啊，这个是在说如何 ...

先算单点，确保得到自旋极化波函数，且确保波函数稳定性，然后读取波函数做结构优化。分两步走。

PS：你想打出的字是倪？他是我师兄

斯卡姆

zjxitcc 发表于 2020-9-9 20:14
先算单点，确保得到自旋极化波函数，且确保波函数稳定性，然后读取波函数做结构优化。分两步走。

PS： ...

啊，有点明白了，还以为是先优化出波函数，再算单点呢；

那这个2步走完成之后，是不是还得用之前的波函数再算一次单点？

啊，不好意思，认错人啦；我以为量子化学公众号是倪老师的呢。

zjxitcc

斯卡姆 发表于 2020-9-9 21:03
啊，有点明白了，还以为是先优化出波函数，再算单点呢；

那这个2步走完成之后，是不是还得用之前的波 ...

没看懂你在说啥

斯卡姆

zjxitcc 发表于 2020-9-10 09:39
没看懂你在说啥

可能是我的表述有些问题
下面以丁烷双自由基为例，我的理解是：
第一步做单点的目的是获得稳定的对称破缺波函数，为下一步结构优化提供初猜，但是这个时候得到的能量可能并不是在最合理的结构下的能量。

第2步做优化，这时候是在上一步自旋极化波函数的基础上优化结构，得到体系的合理结构

第3步，利用上面的这个结构，用更好的方法基组重新再算一次单点；或者直接读取上边的能量也可以；这时候的能量应该会比第一步的能量低些

zjxitcc

斯卡姆 发表于 2020-9-11 09:58
可能是我的表述有些问题
下面以丁烷双自由基为例，我的理解是：
第一步做单点的目的是获得稳定的对 ...

没问题

糖糖糖豆9988

zjxitcc 发表于 2020-9-9 20:14
先算单点，确保得到自旋极化波函数，且确保波函数稳定性，然后读取波函数做结构优化。分两步走。

PS： ...

老师您好，看了您的回复我不明白①为什么先算单点？以“谈谈片段组合波函数与自旋极化单重态”一文的计算单重态双自由基·CH2-CH2-CH2-CH2·为例，我的理解是第一步：# ub3lyp/6-31g(d) guess(fragment=2)得到自旋极化初猜波函数；第二步：# ub3lyp/6-31g(d) guess=read nosymm，读取上述初猜波函数计算单点能，同时可以在multiwfn查看自旋密度确认其双自由基特征；我还有有一个问题②如果我想计算·CH2-CH2-CH2-CH2·精度较高的能量，是不是可以第一步：# ub3lyp/6-31g(d) guess(fragment=2)得到自旋极化初猜波函数；第二步：# ub3lyp/6-31g(d) opt guess=read nosymm优化；第三步CCSD(T)/cc-pVTZ算能量呢？不过双自由基好像还是要用CASSCF方法计算？