计算化学公社
标题: 请教luscus的使用方法,如何用来帮助编写molcas的输入文件。 [打印本页]
作者Author: TQcomputation 时间: 7 day ago
标题: 请教luscus的使用方法,如何用来帮助编写molcas的输入文件。
请教各位老师,如何用luscus编写OpenMolcas的输入文件,查看输出结果,特别是如何设置活性空间参数?
我从官网下载luscus_0.8.6_windows.zip解压到笔记本电脑的C盘,然后双击luscus_0.8.6_windows文件夹里面的luscus.exe后,虽然能打开luscus面板,并通过File-Open也能打开.xyz格式的文件显示出分子结构,但是其它格式的文件都不能正常打开!例如File-Open打开.input文件时,报错提示“文件不能识别”。打开.molden文件时,面板既不弹出报错信息,也不显色分子结构,面板框依然是白板,与此同时文件夹中新生成两个后缀为lus的文件,其中一个是0KB的空白文件,另一个1KB的文件显示如下:
0
File generated by luscus version 0.8.6 beta
<END>
<BREAK>
直接双击luscus_0.8.6_windows里面的其它.exe程序文件,例如molden2lus.exe和gamess2lus.exe等,均闪一下就没有反应了。后面我在luscus_0.8.6_windows文件夹里新建luscus.bat文件及命令,通过双击luscus.bat同样可以打开luscus面板,但是通过File-Open仍然只能打开 xyz 文件,其它格式的文件依然出现相同问题。我所用电脑是Windows 10家庭中文版,64位操作系统。换其它电脑尝试过,也是这样的问题。
其中luscus_0.8.6_windows文件夹里面的plugin.rc文件用记事本打开,内容如下:
libpath=. extension=outdescription="GAMESS output" forward=gamess2lus.exe
libpath=. extension=GvOrbdescription="Molcas orbital data" backward=lus2gvorb.exe
libpath=. extension=moldendescription=molden forward=molden2lus.exe
libpath=. extension=xyzdescription="cartesian coordinates" forward=xyz2gv.exebackward=gv2xyz.exe
(, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 0)
作者Author: zjxitcc 时间: 7 day ago
本帖最后由 zjxitcc 于 2026-2-9 10:42 编辑
不建议使用LUSCUS,完全没必要安装。如果每个量化软件都搭配一个可视化界面看轨道,那5个量化软件岂不要5个可视化软件?
建议直接使用生产力工具MOKIT,无论调用哪个量化软件做计算,算完就有fch文件(程序自动做波函数文件转化),可以用GaussView/Multiwfn打开看轨道。以双原子分子C2分子为例,经典的双键-四键争议分子,但从CASSCF计算角度讲,大家都会使用CASSCF(8,8)。输入文件就用Gaussian输入文件(gjf)
- %mem=8GB
- %nprocshared=4
- #p CASSCF(8,8)/cc-pVDZ guess(fragment=2)
- mokit{CASSCF_prog=OpenMolcas}
- 0 1 0 5 0 -5
- C(fragment=1) 0.0 0.0 0.0
- C(fragment=2) 0.0 0.0 1.5
复制代码 提交任务,即运行
- 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可视化,不需要安装其他
(, 下载次数 Times of downloads: 0)
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.
当前目录下所有文件一应俱全
- C2.gjf
- C2.out
- C2_uhf.fch
- C2_uhf.gjf
- C2_uhf_gvb4_CASSCF.INPORB
- C2_uhf_gvb4_CASSCF.input
- C2_uhf_gvb4_CASSCF_NO.fch
- C2_uhf_gvb4_CASSCF.out
- C2_uhf_gvb4_CASSCF.RasOrb
- C2_uhf_gvb4_CASSCF.RasOrb.1
- C2_uhf_gvb4_CASSCF.rasscf.h5
- C2_uhf_gvb4_CASSCF.rasscf.molden
- C2_uhf_gvb4_CASSCF.SpdOrb.1
- C2_uhf.log
- C2_uhf_uno_asrot2gvb4.dat
- C2_uhf_uno_asrot2gvb4.gms
- C2_uhf_uno_asrot2gvb4.inp
- C2_uhf_uno_asrot2gvb4_s.dat
- C2_uhf_uno_asrot2gvb4_s.fch
- C2_uhf_uno_asrot.fch
- C2_uhf_uno_asrot.out
- C2_uhf_uno_asrot.py
- C2_uhf_uno.fch
- xmldump
复制代码 如果你还想算NEVPT2,只需把gjf文件中的方法从CASSCF(8,8)改成NEVPT2(8,8)。
作者Author: TQcomputation 时间: 7 day ago
感谢邹老师的回复!我进行量化计算的超算平台不支持联网安装,在个人目录里已安装Gaussian和OpenMolcas,是不是还需要安装PySCF和GAMESS才行,然后再安装预编译版的MOKIT?
作者Author: zjxitcc 时间: 7 day ago
本帖最后由 zjxitcc 于 2026-2-9 18:45 编辑
是的。MOKIT可以用预编译版,也可以自己本地编译。相应安装教程全都有https://gitlab.com/jxzou/qcinstall
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |