gamess-us运行时自旋多重度为2老是出现911错误，修改过内存还是出错

nzh19981019

老师，计算的作业时自旋多重度为2，具gamsss说明书中描述，mult=2，方法要切换成CCTYP=CCSD或者CR-CCL。当我更改后一直报错。我在运行gamess-us（GAMESS-US/20230930-R2-gompi-2022a）运算过程老是出现911错误，看咱们相关的贴子说内存不足，我更改了好几次内存都不行，给增大的核数也不行，希望老师给予解答。谢谢您。万分感谢！

zjxitcc

（1）是gamess，不是gamsss；
（2）你需要用文字清晰写出 你的体系电荷、自旋多重度，你想做什么计算，这样读者才能对照你的文字，判断你的输入文件是否合理。如果没有给出计算体系信息和计算目标，变成无头苍蝇，自然无法得到有效回答。
回答示例：我的体系是H2O与NH3的复合物，体系净电荷为0，自旋多重度为1，每个片段的电荷和自旋多重度均是0 1，我现在想在CCSD/aug-cc-pVTZ级别下进行LMO-EDA计算，请问我的输入文件是否合理，如果不合理，应当如何改进？

nzh19981019

老师，我是使用gamess-us进行能量分解，我的体系是Ar与BF2+的复合物，体系净电荷为1，自旋多重度为2，第一和第二片段的电荷和自旋多重度均是0 1，1 2。我现在想在CCSD/aug-cc-pVTZ或者CR-CCL级别下进行LMO-EDA计算，请问我的输入文件是否合理，如果不合理，应当如何改进？无论我怎么更改输入文件，均出现ERROR: ROHF'S CCTYP MUST BE CCSD OR CR-CCL, WITH SERIAL EXECUTION。
但是自旋多重度为1时，在CCSD(T)/aug-cc-pVTZ级别下进行LMO-EDA计算，可以正常运算的。

zjxitcc

nzh19981019 发表于 2026-4-3 19:38
老师，我是使用gamess-us进行能量分解，我的体系是Ar与BF2+的复合物，体系净电荷为1，自旋多重度为2，第一 ...

回复其他人时，要点击对应楼层的“回复 Reply”，否则对方收不到回复提醒。ROHF和ROCCSD 比 RHF和RCCSD 计算难度大很多。后者计算简单，不能推理前者计算也简单。

SERIAL EXECUTION的意思是单核串行计算，无法并行，这是GAMESS程序ROCCSD方法的限制。当前你上传的输入文件没法做ROCCSD的LMO-EDA计算，可参考以下关键词

1. $CONTRL SCFTYP=ROHF RUNTYP=EDA ICHARG=1 MULT=2 NOSYM=1 ICUT=11
   
2.    MAXIT=200 ISPHER=1 DFTTYP=NONE CCTYP=CCSD $END
   
3. $SYSTEM MWORDS=1500 MEMDDI=1125 $END
   
4. $BASIS GBASIS=acct $END
   
5. $SCF DIRSCF=.T. DIIS=.T. SOSCF=.F. SHIFT=.T. RESET=.F. $END
   
6. $LMOEDA MATOM(1)=1,4 MCHARG(1)=0,1 MMULT(1)=1,2 $END
   
7. $CCINP MAXCC=100 $END
   
8. $GUESS GUESS=HCORE $END
   

复制代码

底下$DATA部分就不展示了。注意每个$符号之前必须有且仅有1个空格。

事实上，我们可以采用更简洁、自动化的方式产生LMO-EDA输入，只需要写一个ArBF3.gjf文件

1. %mem=180GB
   
2. %nprocshared=64
   
3. #p ROCCSD/aug-cc-pVTZ guess(fragment=2)
   
4. 
   
5. {lmo}
   
6. 
   
7. 1 2 0 1 1 2
   
8. Ar(fragment=1)   0.00000000    0.01224939   -1.93470092
   
9. F(fragment=2)    0.00000000    0.83616032    0.01421492
   
10. B(fragment=2)    0.00000000   -0.08183149    1.11978199
    
11. F(fragment=2)    0.00000000   -1.32793992    0.80653336
    
12. F(fragment=2)    0.00000000    0.51171775    2.25183303

复制代码

这个文件就是Gaussian的输入文件，大家都会写，不需要经过各种EDA教程培训。片段已经在gjf文件中划分好；如果用GaussView打开，点击主面板上的Tools -> Atom Groups，可以看到各个片段以不同颜色显示，一目了然。提交任务，即运行

1. frag_guess_wfn ArBF3.gjf >ArBF3.out 2>&1

复制代码

几十秒算完，获得文件ArBF3.inp，包含坐标、基组、ROHF轨道数据，ROCCSD和LMO-EDA关键词全都写好，可以确保ROHF收敛，直接进入ROCCSD计算步骤，这里展示inp文件开头

1. $CONTRL SCFTYP=ROHF RUNTYP=EDA ICHARG=1 MULT=2 NOSYM=1 ICUT=11
   
2.    MAXIT=200 ISPHER=1 DFTTYP=NONE CCTYP=CCSD $END
   
3. $SYSTEM MWORDS=1500 MEMDDI=1125 $END
   
4. $SCF DIRSCF=.T. DIIS=.T. SOSCF=.F. SHIFT=.T. RESET=.F. $END
   
5. $LMOEDA MATOM(1)=1,4 MCHARG(1)=0,1 MMULT(1)=1,2
   
6.    EDATYP=NONE RDVECM=.T. $END
   
7. $CCINP MAXCC=100 $END
   
8. $GUESS GUESS=HCORE $END
   

复制代码

我上面说的“可参考以下关键词”就是从这组关键词修改而来的。frag_guess_wfn运行时有输出提示 应当使用单核串行计算。这个自动产生的ArBF3.inp文件可以用XEDA插件版GAMESS 提交LMO-EDA计算，例如运行

1. xeda ArBF3.inp 00 1 >ArBF3.gms 2>&1

复制代码

如果未获取、没安装XEDA插件版GAMESS，想使用官方GAMESS做LMO-EDA计算，则需使用本回答一开始展示的那套关键词，并且无法读取收敛的轨道数据。
ArBF3.gjf (389 Bytes, 下载次数 Times of downloads: 1)

4 day ago 上传 Uploaded

点击下载
Click to download

， ArBF3.out (1.51 KB, 下载次数 Times of downloads: 1)

4 day ago 上传 Uploaded

点击下载
Click to download

nzh19981019

好的，非常感谢老师的详细解决方案。马上尝试。真心的深深感激。