计算化学公社

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

[辅助/分析程序] MOKIT做CASSCF计算出错:*_uhf_uno.fch未正常产生

[复制链接 Copy URL]

178

帖子

3

威望

1922

eV
积分
2160

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 Kalinite 于 2022-12-21 00:02 编辑

请问大家,MOKIT做CASSCF force计算出现以下报错应该如何解决:
ERROR in subroutine require_file_exist: file does not exist!
Filename=fe3oxy-casscf_uhf_uno.fch
以上缺少的文件未能正常产生。

输入文件:
  1. %mem=96GB
  2. %nprocshared=32
  3. #p CASSCF/6-31G(d)

  4. mokit{force}

  5. -1 3
  6. 坐标略
复制代码

临时文件uno.out的内容如下:
  1. nbf=  518
  2. nif=  518
  3. ON_criteria=  0.0000100
  4. uno_thres=  0.0000100

  5. The orthonormality of initial Alpha MO:
  6. maxv=    0.0000000000
  7. abs_mean=    0.0000000000

  8. The orthonormality of initial Beta MO:
  9. maxv=    0.0000000000
  10. abs_mean=    0.0000000000

  11. ndb  =  107
  12. nact =   24
  13. nact0=   22
  14. idx  =  108  154    2

  15. The orthonormality of final Alpha MO:
  16. maxv=    0.0000000001
  17. abs_mean=    0.0000000000

  18. ndb  =  107
  19. nact =   24
  20. nact0=   22
  21. idx  =  108  154    2
复制代码


评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
zjxitcc + 1 赞!

查看全部评分 View all ratings

178

帖子

3

威望

1922

eV
积分
2160

Level 5 (御坂)

2#
 楼主 Author| 发表于 Post on 2022-12-21 02:09:17 | 只看该作者 Only view this author
已解决。
邹神:是fch_u2r的问题,部分用户出现过。

对于开壳层体系,手动调用fch_u2r:
  1. $MOKIT_ROOT/bin/fch_u2r *_uhf.fch
  2. mv *_uhf_r.fch *_uhf_uno.fch
复制代码
然后续算。

或将$MOKIT_ROOT/src/automr.f90调用fch_u2r的代码修改:
  1. write(fid1,'(A)') "with os.popen('fch_u2r "//TRIM(hf_fch)//' '//TRIM(uno_fch)//"') as run:"
  2. write(fid1,'(A)') '  null=run.read()'
复制代码
然后make automr即可。
原因应该是python以os.system/popen等方式调用子进程时,子进程并未结束,但主进程会继续运行。当源代码中sleep(1)的时间不够完成fch_u2r时就会出错。

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 17:00 , Processed in 0.156872 second(s), 21 queries , Gzip On.

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