计算化学公社

标题: gamess-us运行时自旋多重度为2老是出现911错误，修改过内存还是出错 [打印本页]

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

作者
Author: zjxitcc 时间: 2026-4-3 19:25
本帖最后由 zjxitcc 于 2026-4-3 19:27 编辑

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

作者
Author: nzh19981019 时间: 2026-4-3 19:38
老师，我是使用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计算，可以正常运算的。

作者
Author: zjxitcc 时间: 4 day ago
本帖最后由 zjxitcc 于 2026-4-10 22:19 编辑

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

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

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

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

复制代码

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

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

%mem=180GB
%nprocshared=64
#p ROCCSD/aug-cc-pVTZ guess(fragment=2)
{lmo}
1 2 0 1 1 2
Ar(fragment=1) 0.00000000 0.01224939 -1.93470092
F(fragment=2) 0.00000000 0.83616032 0.01421492
B(fragment=2) 0.00000000 -0.08183149 1.11978199
F(fragment=2) 0.00000000 -1.32793992 0.80653336
F(fragment=2) 0.00000000 0.51171775 2.25183303

复制代码

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

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

复制代码

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

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

复制代码

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

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

复制代码

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

作者
Author: nzh19981019 时间: 3 day ago
好的，非常感谢老师的详细解决方案。马上尝试。真心的深深感激。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)