LMO-EDA分子内弱作用的计算问题

123qwertybobo

想用LMO-EDA计算分子内的弱相互作用，参考了Chemical Physics Letters 635 (2015) 250–256这篇文章，试图重复他们的test分子的Table2(a)的数据。
学着写了一下input，不知道对不对？附件上传，希望大家能够帮我指点迷津，谢谢。

以下是input：

$CONTRL SCFTYP=ROHF DFTTYP=wB97-D RUNTYP=EDA ICHARG=0 MULT=1 MAXIT=200 $END
$BASIS GBASIS=N311 NGAUSS=6 NDFUNC=3 NFFUNC=1 NPFUNC=3 DIFFSP=.T. DIFFS=.T. $END
$LMOEDA MATOM(1)=20,7 MCHARG(1)=0,0 MMULT(1)=-2,2 $END
$GUESS GUESS=HCORE $END
$SYSTEM MWORDS=300 MEMDDI=300 $END
$DATA
LMD-EDA wb97-D/6-311++G**
分子坐标略
$END

wei

作者是基于GKS-EDA做的，LMO-EDA做的结果不可信。你还是先拿GKS-EDA代码吧

123qwertybobo

wei 发表于 2017-9-7 15:35
软件拆开不了，要人为判断合适断键位置。CPL那篇文章的通讯作者当面告诉我的，你信了吧？

ok，明白了。对了，想必您也研究过这篇文章，能给个inp让我学习学习吗？谢谢

wei

软件拆开不了，要人为判断合适断键位置。CPL那篇文章的通讯作者当面告诉我的，你信了吧？

123qwertybobo

wei 发表于 2017-9-6 20:35
不对。
1）他是用gks-eda做的
2）要断键后形成三个输入文件，然后再加加减减

直觉没有这么繁琐，估计是软件能拆开的，如果片段多，难道还手动一个个么

wei

不对。
1）他是用gks-eda做的
2）要断键后形成三个输入文件，然后再加加减减

123qwertybobo

$LMOEDA MATOM(1)=20,7 MCHARG(1)=0,0 MMULT(1)=-2,2 $END     这里一提交，就发现写反了，应该是7,20

123qwertybobo

附件在这