计算化学公社

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

[Molpro] MOLPRO计算出错,RHF不收敛

[复制链接 Copy URL]

2

帖子

0

威望

73

eV
积分
75

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
molpro计算单点能是报错,Error: RHF not converged. This error exit can be avoided using the IGNORE_ERROR option on the ORBITAL directive

202211021731026116..png (40.89 KB, 下载次数 Times of downloads: 22)

202211021731026116..png

ts5.inp

2.15 KB, 下载次数 Times of downloads: 25

输入文件

4289

帖子

4

威望

9546

eV
积分
13915

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2022-11-2 19:15:56 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2022-11-2 19:17 编辑

简单,可以先用Gaussian算个ROHF,然后用fch2com小程序传ROHF轨道至Molpro,后者将会ROHF 1圈收敛,马上进入UCCSD(T)-F12计算。假设基组是jun-cc-pV(D+d)Z,gjf示例
ts5.gjf (4.42 KB, 下载次数 Times of downloads: 43)
几十秒就能算完,然后执行
  1. fch2com ts5.chk
复制代码

获得ts5.com和ts5.a文件。打开ts5.com文件,在末尾追加辅助基和UCCSD(T)-F12关键词
  1. basis={
  2. set,df
  3. default,avdz/mp2fit
  4. set,jk
  5. default,avdz/jkfit
  6. set,ri
  7. default,avdz/optri
  8. }
  9. {uccsd(t)-f12a,ri_basis=ri,df_basis=df,df_basis_exch=jk}
复制代码

提交ts5.com给molpro,搞定。ts5.a会被自动读取。


自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 00:12 , Processed in 0.518862 second(s), 29 queries , Gzip On.

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