简单量子化学问题答疑专帖

糖糖糖豆9988

wzkchem5 发表于 2023-6-27 17:22
检查有没有收敛到正确的自旋对称性破缺波函数。如果波函数是正确的，而且扫描的键长范围够宽，那就是没有 ...

%chk=sfor-ts3_stable.chk
%nprocs=8
#p um062x 6-311g(d,p) nosymm guess=mix stable=opt

坐标

--Link1--  
%oldchk=sfor-ts3_stable.chk        
%chk=sfor-ts3.chk
%nprocs=8
#p opt=modredundant um062x 6-311g(d,p) nosymm guess=read geom=allcheck

B 6 8 S 30 0.050000

老师，如上是我输入文件设置，①我在输入文件中查找<S**2>，发现基本所有的<S**2>= 0.0000，这表示没有收敛到正确波函数吗？②我换成# UM062X/6-31G* guess(mix,always) opt=modredundant nosymm方法寻找是不是可行呢？谢谢老师

wzkchem5

糖糖糖豆9988 发表于 2023-6-27 10:26
%chk=sfor-ts3_stable.chk
%nprocs=8
#p um062x 6-311g(d,p) nosymm guess=mix stable=opt

两个自由基接近完全解离的时候，应该得收敛到自旋破缺波函数，但没有完全解离时可能波函数本来就不该对称性破缺。
如果确实发现完全解离时<S**2>也是0的话，那么要么键长从长往短扫（只需要在一开始准备对称性破缺初猜），要么guess=(mix,always)，其中前者应该更快，因为节省SCF迭代步数

糖糖糖豆9988

wzkchem5 发表于 2023-6-27 17:51
两个自由基接近完全解离的时候，应该得收敛到自旋破缺波函数，但没有完全解离时可能波函数本来就不该对称 ...

好的，谢谢老师我这就试一试~

tsl666

sobereva 发表于 2023-6-26 23:11
没定义赝势。低级错误
仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入

谢谢sobereva老师的解答，已经算出来啦，感恩

sai77

请问用Multiwfn作弱相互作用分析，ORCA生成的gbw文件和Gaussian生成的chk文件可以混用吗？
想比较两个过渡态各自的弱相互作用，其中一个过渡态的chk文件被误删了，可以用gbw代替，与另一个过渡态的chk比较吗？
谢谢！

wzkchem5

sai77 发表于 2023-6-28 05:17
请问用Multiwfn作弱相互作用分析，ORCA生成的gbw文件和Gaussian生成的chk文件可以混用吗？
想比较两个过渡 ...

不太建议这样做。即使结果可比，你想想你文章的计算部分怎么写，如果写两个计算都是用orca或者都是用高斯算的，那是撒谎；但是说一个用高斯算一个用orca算，审稿人必然觉得可疑，要你证明两者的结果可比，到时候还是要把其中一个计算换软件重算。
Multiwfn可以把chk转成gbw格式给orca做初猜，也可以把gbw转成chk格式给高斯做初猜，你把其中一个波函数转一下格式、再做SCF收敛一下，这样两个过渡态就都是同一个软件算出来的了，不会被质疑了，但是计算又比从头跑一个单点能计算要快得多

yjk

各位老师好，最近在算一个分子的激发三重态基于优化好的s0用TD-DFT方法分别计算了T1,T2,T3然后发现T3的能量比T1的能量还低，我想问一下这时候优化的T2态是不是就是重新排布成T1了？

sai77

wzkchem5 发表于 2023-6-28 16:24
不太建议这样做。即使结果可比，你想想你文章的计算部分怎么写，如果写两个计算都是用orca或者都是用高斯 ...

明白了，谢谢！！

zhangxiaolu

snljty2 发表于 2023-6-21 23:45
SSH客户端那是你登录服务器用的软件，跟你cmd里直接用OpenSSH的ssh命令没啥本质区别，这不是计算软件。计算 ...

输入文件贴上了，可以麻烦再看下吗？

zhangxiaolu

sobereva 发表于 2023-6-22 00:37
没输入输出文件没法说

输入文件贴上了，麻烦老师再看下呢

sobereva

zhangxiaolu 发表于 2023-6-28 21:23
输入文件贴上了，麻烦老师再看下呢

输入文件本身没明显问题，可能你用的高斯版本的bug。尝试G16，或者尝试加上int=ultrafine，或者不带溶剂模型，或者换别的溶剂模型如SMD能否避免

sobereva

yjk 发表于 2023-6-28 17:16
各位老师好，最近在算一个分子的激发三重态基于优化好的s0用TD-DFT方法分别计算了T1,T2,T3然后发现T3的能量 ...

你的图里根本没体现出T3，哪显示那是T3了？

reader

老师您好，我使用B2PLYP/ma-qzvp计算阴离子的结合，一直显示Convergence failure -- run terminated.我试过很多方法都没法使SCF收敛，这是为什么呢，但是有的结构只需要输入# b2plyp/gen em=gd3bj就可以收敛，我这边有几个类似的结构改了很多次# b2plyp/gen em=gd3bj SCF(novaracc,noincfock) int=acc2e=12 IOp(5/37=15)都没法收敛，想问一下老师是为什么

sobereva

reader 发表于 2023-6-29 08:58
老师您好，我使用B2PLYP/ma-qzvp计算阴离子的结合，一直显示Convergence failure -- run terminated.我试过 ...

发帖的时候看清楚板块再发，这在首页的公告栏里已经明确提醒了，求助帖别发到资源分享子版块、量子化学问题不得发到分子模拟版块，这次给你移动了，以后发错地方直接删帖+扣分处理。

解决SCF不收敛问题的方法
http://sobereva.com/61


甭写IOp(5/37=15)

yjk

sobereva 发表于 2023-6-28 23:58
输入文件本身没明显问题，可能你用的高斯版本的bug。尝试G16，或者尝试加上int=ultrafine，或者不带溶剂 ...

抱歉老师，我刚才上传了。我是分别优化的,T2目前还没结束。