计算化学公社

 找回密码 Forget password
 注册 Register

请教luscus的使用方法,如何用来帮助编写molcas的输入文件。

查看数: 142 | 评论数: 3 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2026-2-9 09:55

正文摘要:

   请教各位老师,如何用luscus编写OpenMolcas的输入文件,查看输出结果,特别是如何设置活性空间参数?      我从官网下载luscus_0.8.6_windows.zip解压到笔记本电脑的C盘,然后双击 ...

回复 Reply

zjxitcc 发表于 Post on 7 day ago
本帖最后由 zjxitcc 于 2026-2-9 18:45 编辑
TQcomputation 发表于 2026-2-9 18:23
感谢邹老师的回复!我进行量化计算的超算平台不支持联网安装,在个人目录里已安装Gaussian和OpenMolcas,是 ...


是的。MOKIT可以用预编译版,也可以自己本地编译。相应安装教程全都有https://gitlab.com/jxzou/qcinstall
TQcomputation 发表于 Post on 7 day ago
感谢邹老师的回复!我进行量化计算的超算平台不支持联网安装,在个人目录里已安装Gaussian和OpenMolcas,是不是还需要安装PySCF和GAMESS才行,然后再安装预编译版的MOKIT?
zjxitcc 发表于 Post on 7 day ago
本帖最后由 zjxitcc 于 2026-2-9 10:42 编辑

不建议使用LUSCUS,完全没必要安装。如果每个量化软件都搭配一个可视化界面看轨道,那5个量化软件岂不要5个可视化软件?

建议直接使用生产力工具MOKIT,无论调用哪个量化软件做计算,算完就有fch文件(程序自动做波函数文件转化),可以用GaussView/Multiwfn打开看轨道。以双原子分子C2分子为例,经典的双键-四键争议分子,但从CASSCF计算角度讲,大家都会使用CASSCF(8,8)。输入文件就用Gaussian输入文件(gjf)
  1. %mem=8GB
  2. %nprocshared=4
  3. #p CASSCF(8,8)/cc-pVDZ guess(fragment=2)

  4. mokit{CASSCF_prog=OpenMolcas}

  5. 0 1 0 5 0 -5
  6. C(fragment=1)    0.0   0.0   0.0
  7. C(fragment=2)    0.0   0.0   1.5
复制代码
提交任务,即运行
  1. automr C2.gjf >C2.out 2>&1
复制代码
十几秒完成,包括HF计算,产生UNO轨道、局域配对UNO轨道,GVB计算,CASSCF计算,活性轨道是自动构建的。这里特地写了CASSCF_prog=OpenMolcas,调用你想用的OpenMolcas做计算。注意,我们在没学OpenMolcas的前提下,程序就自动传轨道给它 令其完成了CASSCF计算。算完之后有CASSCF自然轨道 及其 占据数,保存在C2_uhf_gvb4_CASSCF_NO.fch文件中,任何量化软件的轨道文件只需用GaussView/Multiwfn可视化,不需要安装其他
C2.out文件中有各种方法的能量
E(UHF) =       -75.46417879 a.u., <S**2>=  2.305
E(GVB) =       -75.47185919 a.u.
E(CASCI)  =       -75.57068570 a.u.
E(CASSCF) =       -75.57685721 a.u.
当前目录下所有文件一应俱全
  1. C2.gjf
  2. C2.out
  3. C2_uhf.fch
  4. C2_uhf.gjf
  5. C2_uhf_gvb4_CASSCF.INPORB
  6. C2_uhf_gvb4_CASSCF.input
  7. C2_uhf_gvb4_CASSCF_NO.fch
  8. C2_uhf_gvb4_CASSCF.out
  9. C2_uhf_gvb4_CASSCF.RasOrb
  10. C2_uhf_gvb4_CASSCF.RasOrb.1
  11. C2_uhf_gvb4_CASSCF.rasscf.h5
  12. C2_uhf_gvb4_CASSCF.rasscf.molden
  13. C2_uhf_gvb4_CASSCF.SpdOrb.1
  14. C2_uhf.log
  15. C2_uhf_uno_asrot2gvb4.dat
  16. C2_uhf_uno_asrot2gvb4.gms
  17. C2_uhf_uno_asrot2gvb4.inp
  18. C2_uhf_uno_asrot2gvb4_s.dat
  19. C2_uhf_uno_asrot2gvb4_s.fch
  20. C2_uhf_uno_asrot.fch
  21. C2_uhf_uno_asrot.out
  22. C2_uhf_uno_asrot.py
  23. C2_uhf_uno.fch
  24. xmldump
复制代码
如果你还想算NEVPT2,只需把gjf文件中的方法从CASSCF(8,8)改成NEVPT2(8,8)。

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

GMT+8, 2026-2-16 20:22 , Processed in 0.173451 second(s), 25 queries , Gzip On.

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