计算化学公社

标题: gaussian做142个原子的体系用16核机子优化需要多久呢？ [打印本页]

作者
Author: wad 时间: 2023-4-2 14:17
标题: gaussian做142个原子的体系用16核机子优化需要多久呢？
本帖最后由 wad 于 2023-4-2 14:20 编辑

各位大佬、社长好，
我目前希望用gaussian优化一个142个原子的分子，希望得到.chk文件从而计算resp电荷。
这是我的输入文件
%nprocshared=16
%mem=16GB
%chk=\home\zby\g16\aCD.chk
# opt hf/6-31g(d) nosymm scf=qc pop=mk iop(6/33=2,6/42=1,6/50=1) geom=nocrowd

因为考虑到体系不算小，参数都设置得不精细。
我是16核的机子，3080Ti，想问问这样一个体系大概多久能够跑完，目前已经优化了十几个小时了。如果需要太长时间的话，有没有更好的方法能够解决呢？
谢谢！

作者
Author: sobereva 时间: 2023-4-2 14:32
用scf=qc自取其辱。若无碰到SCF不收敛，绝对不要写这个。就算碰到了，也别优先考虑用scf=qc，仔细看
解决SCF不收敛问题的方法
http://sobereva.com/61

绝对甭用HF优化，该用什么好好看下文。通常优先考虑用B3LYP-D3(BJ)
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）
谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557（http://bbs.keinsci.com/thread-17899-1-1.html）

写nosymm完全多余

自行监控优化过程，别傻等着。对于按下文用Multiwfn计算RESP电荷的目的，通常opt=loose用较宽松的收敛限就够了。写6/50=1完全多余
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）

如果缺乏量子化学计算常识，远不如直接用下文的脚本，省事极了，还避免犯低级错误
计算RESP原子电荷的超级懒人脚本（一行命令就算出结果）
http://sobereva.com/476（http://bbs.keinsci.com/thread-12858-1-1.html）

如果初始结构比较烂，先用PM6D3或PM7做预优化。

作者
Author: wzkchem5 时间: 2023-4-2 16:12
按楼上的要求修改理论级别以后，可以按以下方法估计计算时间：先跑到SCF计算了至少一步，看看第一步需要多久，然后乘以10~20（一般纯有机体系大概20步左右收敛，但后面的迭代比前面的迭代快，所以乘以的数要略小于预期的迭代收敛步数），如果分子含有金属或者自由基，可以预期收敛困难，那么可以乘以30~50。这个只是数量级的估计，会很不准，但是很难做到更准了。

作者
Author: wad 时间: 2023-4-2 16:27

sobereva 发表于 2023-4-2 14:32
用scf=qc自取其辱。若无碰到SCF不收敛，绝对不要写这个。就算碰到了，也别优先考虑用scf=qc，仔细看
解决S ...

明白了，谢谢sob老师

作者
Author: wad 时间: 2023-4-2 16:27

wzkchem5 发表于 2023-4-2 16:12
按楼上的要求修改理论级别以后，可以按以下方法估计计算时间：先跑到SCF计算了至少一步，看看第一步需要多 ...

明白了，谢谢

作者
Author: mfdsrax2 时间: 2023-4-3 10:10
高斯只会用CPU计算，所以你的3080Ti一点作用也没有

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