计算化学公社

标题: MOLPRO计算出错，RHF不收敛 [打印本页]
作者
Author: ZQX    时间: 2022-11-2 17:32
标题: MOLPRO计算出错，RHF不收敛
molpro计算单点能是报错，Error: RHF not converged. This error exit can be avoided using the IGNORE_ERROR option on the ORBITAL directive
作者
Author: zjxitcc    时间: 2022-11-2 19:15
本帖最后由 zjxitcc 于 2022-11-2 19:17 编辑

简单，可以先用Gaussian算个ROHF，然后用fch2com小程序传ROHF轨道至Molpro，后者将会ROHF 1圈收敛，马上进入UCCSD(T)-F12计算。假设基组是jun-cc-pV(D+d)Z，gjf示例
(, 下载次数 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会被自动读取。






欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3