计算化学公社

 找回密码 Forget password
 注册 Register
Views: 8097|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] ORCA 优化UDP失败

[复制链接 Copy URL]

156

帖子

0

威望

1637

eV
积分
1793

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
使用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,也都是不收敛,哪位朋友能帮忙看看,谢谢!


6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

2#
发表于 Post on 2020-11-15 10:24:35 | 只看该作者 Only view this author
先拿xtb预优化一下。或者先甭用tightSCF,等结构大致收敛了再上tightSCF

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

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

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

339

帖子

0

威望

5053

eV
积分
5392

Level 6 (一方通行)

3#
发表于 Post on 2020-11-15 10:35:27 来自手机 | 只看该作者 Only view this author
阴离子带这多电荷,确定不加个溶剂压一压

156

帖子

0

威望

1637

eV
积分
1793

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2020-11-15 19:55:04 | 只看该作者 Only view this author
sobereva 发表于 2020-11-15 10:24
先拿xtb预优化一下。或者先甭用tightSCF,等结构大致收敛了再上tightSCF

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

谢谢SOB!

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

好的,下次直接上传输入输出文件了!

156

帖子

0

威望

1637

eV
积分
1793

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2020-11-15 19:56:18 | 只看该作者 Only view this author
chrinide 发表于 2020-11-15 10:35
阴离子带这多电荷,确定不加个溶剂压一压

谢谢建议!

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

在用xtb优化的时候,加了水了!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-25 14:11 , Processed in 0.177881 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list