本帖最后由 糖糖糖豆9988 于 2022-9-20 20:23 编辑 zjxitcc 发表于 2020-9-9 20:14 老师您好，看了您的回复我不明白①为什么先算单点？以“谈谈片段组合波函数与自旋极化单重态”一文的计算单重态双自由基·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方法计算？ |
斯卡姆 发表于 2020-9-11 09:58 没问题 |
zjxitcc 发表于 2020-9-10 09:39  可能是我的表述有些问题下面以丁烷双自由基为例，我的理解是： 第一步做单点的目的是获得稳定的对称破缺波函数，为下一步结构优化提供初猜，但是这个时候得到的能量可能并不是在最合理的结构下的能量。 第2步做优化，这时候是在上一步自旋极化波函数的基础上优化结构，得到体系的合理结构 第3步，利用上面的这个结构，用更好的方法基组重新再算一次单点；或者直接读取上边的能量也可以；这时候的能量应该会比第一步的能量低些 |
斯卡姆 发表于 2020-9-9 21:03 没看懂你在说啥 |
zjxitcc 发表于 2020-9-9 20:14 啊，有点明白了，还以为是先优化出波函数，再算单点呢； 那这个2步走完成之后，是不是还得用之前的波函数再算一次单点？ 啊，不好意思，认错人啦；我以为量子化学公众号是倪老师的呢  。 |
斯卡姆 发表于 2020-9-9 20:01 先算单点，确保得到自旋极化波函数，且确保波函数稳定性，然后读取波函数做结构优化。分两步走。 PS：你想打出的字是倪？他是我师兄  |
zjxitcc 发表于 2020-9-9 13:06 啊，谢谢ni老师的回复。我的目的是优化丁烷双自由基的结构，获得含有对称破缺的函数。啊，这个是在说如何得到对称破缺的初猜波函数。 |
|
问题就是你的“优化丁烷双自由基的结构”。首先你得确定这个练习是让你学会算个单点，还是有让你做结构优化，不要自己添加多余的步骤。这双自由基你还优化，那就是要么断掉要么成键，无论哪种都不是双自由基。 另外，真正的双自由基结构也不是你这么优化的，缺乏UHF/UDFT的基本常识，看sob老师博文中介绍的自旋极化单重态《谈谈片段组合波函数与自旋极化单重态》http://sobereva.com/82 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-28 22:18 , Processed in 0.211831 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites