计算化学公社

标题: ORCA 优化UDP失败 [打印本页]

---

作者
Author: beyond    时间: 2020-11-15 09:06
标题: ORCA 优化UDP失败
使用ORCA优化一个小分子UDP （PDB ID 3U2U中）， 第一步就失败了，
! B3LYP D3 RIJCOSX def2-SVP def2/J TIGHTSCF Grid4 OPT
%pal nprocs 32
end
%maxcore 2048

*xyz -3 1
N         40.20300        8.90800       49.91500
C         40.03900        8.22200       51.07500
N         40.29700        6.89600       51.12400
C         40.70600        6.23000       50.02100
C         40.89200        6.92100       48.81800
C         40.60500        8.27200       48.78200
O         39.65800        8.83100       52.10900
O         40.92400        4.99200       50.10800
C         39.84800       10.33700       49.91800
C         40.98200       11.25500       50.32600
O         40.44900       12.47600       50.83100
C         41.64300       11.52800       49.02100
C         40.42500       11.72900       48.12500
O         39.50500       10.74200       48.56300
O         42.55000       12.61800       49.03400
C         40.68900       11.50700       46.62600
O         41.55400       10.36000       46.37500
P         42.53600       10.28500       45.11000
O         43.61100        9.29700       45.35900
O         42.89200       11.73900       44.65300
O         41.55400        9.68800       43.98600
P         41.77800        9.92500       42.42200
O         40.68900        9.09700       41.76700
O         43.11800        9.37100       42.08300
O         41.56400       11.41100       42.14200
H         40.18300        6.39700       51.99500
H         41.25200        6.40900       47.93800
H         40.70000        8.82100       47.85700
H         38.99300       10.49900       50.57500
H         41.65900       10.77900       51.03600
H         42.17300       10.63200       48.70000
H         40.00700       12.72300       48.28500
H         39.73600       11.34200       46.12400
H         41.15700       12.40100       46.21300
H         40.01900       12.31400       51.67400
H         42.93200       12.72800       48.16000
*

还试过其它的functional，还有半经验的PM3，也都是第一步SCF不收敛，能量一直振荡。
也调整过初始结构，改过二面角以让小分子更加extended，也都是不收敛，哪位朋友能帮忙看看，谢谢！

---

作者
Author: sobereva    时间: 2020-11-15 10:24
先拿xtb预优化一下。或者先甭用tightSCF，等结构大致收敛了再上tightSCF

另外，都写了Grid4，没道理不把COSX积分格点也提升

以后问这类问题尽量上传输出文件，至少把收敛过程的信息贴出来，否则他人难以判断怎么处理

---

作者
Author: chrinide    时间: 2020-11-15 10:35
阴离子带这多电荷，确定不加个溶剂压一压

---

作者
Author: beyond    时间: 2020-11-15 19:55

sobereva 发表于 2020-11-15 10:24
先拿xtb预优化一下。或者先甭用tightSCF，等结构大致收敛了再上tightSCF

另外，都写了Grid4，没道理不把 ...

谢谢SOB！

用xtb预先优化一下，再用DFT，scf很快就收敛了！

好的，下次直接上传输入输出文件了！

---

作者
Author: beyond    时间: 2020-11-15 19:56

chrinide 发表于 2020-11-15 10:35
阴离子带这多电荷，确定不加个溶剂压一压

谢谢建议！

之前也有试过，只是也是没收敛，没有在这里提了！

在用xtb优化的时候，加了水了！

---

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

Powered by Discuz! X3.3