求助gamess计算开壳层单重态总是收敛到闭壳层

498746012 · 发表于 Post on 2021-4-1 09:45:14

各位老师打了，用gamess计算开壳层单重态总是收敛到闭壳层，是不是写的关键词有问题呀？下边是输入文件
  $CONTRL SCFTYP=UHF RUNTYP=ENERGY ICHARG=0 MULT=1
  DFTTYP=B3LYPV1R MAXIT=200 ISPHER=1 $END
$SYSTEM MWORDS=500 $END
$basis gbasis=n31 ngauss=6 ndfunc=1  $end
$GUESS  norder=1 IORDER(16)=17,16 JORDER(16)=17,16 $END
$DATA
GAMESS inp format file produced by MOKIT
C1
C  6.0 -0.74742092 1.77656753 0.00000000
H  1.0 -0.62438907 2.32965189 0.92333358
H  1.0 -0.62438907 2.32965189 -0.92333358
C  6.0 -0.74742092 0.30933780 0.00000000
H  1.0 -1.24978876 -0.09122321 0.88294562
H  1.0 -1.24978876 -0.09122321 -0.88294562
C  6.0 0.74742092 -0.30933780 0.00000000
H  1.0 1.24978876 0.09122321 -0.88294562
H  1.0 1.24978876 0.09122321 0.88294562
C  6.0 0.74742092 -1.77656753 0.00000000
H  1.0 0.62438907 -2.32965189 -0.92333358
H  1.0 0.62438907 -2.32965189 0.92333358
$end

zjxitcc · 发表于 Post on 2021-4-1 09:59:16

本帖最后由 zjxitcc 于 2021-4-1 10:01 编辑

我问个问题哈。。。你都用MOKIT的fch2inp小程序产生GAMESS输入文件了，那算出开壳层单重态就是几秒钟的事，为啥还自己删掉基组和轨道呢？

我举个例子，高斯关键词
#p UB3LYP/6-31G(d) 5D 7F nosymm guess=mix stable=opt int=nobasistransform

算完后执行
formchk xxx.chk xxx.fch
fch2inp xxx.fch -uhf

生成GAMESS .inp文件，加个B3LYP关键词，交上去算，几秒结束，迅速收敛（不能1-2圈收敛，因为GAMESS默认格点设置与高斯不同），结果也是开壳层单重态，自旋污染与高斯一模一样。

你是忘记加-uhf参数了么？

498746012 · 发表于 Post on 2021-4-1 10:10:24

zjxitcc 发表于 2021-4-1 09:59
我问个问题哈。。。你都用MOKIT的fch2inp小程序产生GAMESS输入文件了，那算出开壳层单重态就是几秒钟的事， ...

老师您好，我想算开壳层单重态的EDA，无法读取分子轨道，我在调关键词看看能不能算出开壳层单重态

zjxitcc · 发表于 Post on 2021-4-1 10:21:24

498746012 发表于 2021-4-1 10:10
老师您好，我想算开壳层单重态的EDA，无法读取分子轨道，我在调关键词看看能不能算出开壳层单重态

那很简单，GAMESS有提供与高斯guess=mix等价的关键词
$GUESS MIX=.TRUE. $END
用不着IORDER，JORDER。

另外，这只能保证你这个单点能算对，至于EDA里效果如何，还得具体再看。

498746012 · 发表于 Post on 2021-4-1 10:41:26

zjxitcc 发表于 2021-4-1 10:21
那很简单，GAMESS有提供与高斯guess=mix等价的关键词
$GUESS MIX=.TRUE. $END
用不着IORDER，JORDER。
...

谢谢老师，太感谢了

498746012 · 发表于 Post on 2021-4-1 12:36:46

本帖最后由 498746012 于 2021-4-1 12:37 编辑

zjxitcc 发表于 2021-4-1 10:21
那很简单，GAMESS有提供与高斯guess=mix等价的关键词
$GUESS MIX=.TRUE. $END
用不着IORDER，JORDER。
...

老师我的开壳层单重态不收敛，刚开始能量一直下降到-2060.792，后来能量升高的很厉害，有没有关键词可以解决这个问题呀？输入文件如下 $CONTRL SCFTYP=UHF RUNTYP=ENERGY ICHARG=0 MULT=1 NOSYM=1 ICUT=11 PP=READ
DFTTYP=B3LYPV1R MAXIT=200 ISPHER=1 $END
$SYSTEM MWORDS=500 $END
$DFT DC=.T. NRAD=99 NLEB=590 $END
$SCF DIRSCF=.TRUE. DIIS=.TRUE. $END
$GUESS MIX=.TRUE. $END
$DATA

paramecium86 · 发表于 Post on 2021-4-1 13:35:21

498746012 发表于 2021-4-1 12:36
老师我的开壳层单重态不收敛，刚开始能量一直下降到-2060.792，后来能量升高的很厉害，有没有关键词可以 ...

SCF不容易收敛 GAMESS的老毛病了。
如果DIIS不行就换成SOSCF 如果soscf不行，就先DIIS再SOSCF。或者干脆二者都别用，先用noconv算一个。然后再读
比如可以先试试soscf
$scf dirscf=.true. fdiff=.f. noconv=.f. soscf=.t. diis=.f. ethrsh=0.002
shift=.t. damp=.f. $end

遇到不好收敛的 noconv soscf diis damp 这几个关键词true false各种组合都试试

498746012 · 发表于 Post on 2021-4-1 13:54:49

paramecium86 发表于 2021-4-1 13:35
SCF不容易收敛 GAMESS的老毛病了。
如果DIIS不行就换成SOSCF 如果soscf不行，就先DIIS再SOSCF。或者干 ...

谢谢老师，太感谢了

zjxitcc · 发表于 Post on 2021-4-1 14:01:51

498746012 发表于 2021-4-1 12:36
老师我的开壳层单重态不收敛，刚开始能量一直下降到-2060.792，后来能量升高的很厉害，有没有关键词可以 ...

（1）你这个例子与帖子一开始说的例子不是同一个体系，应该讲清楚。乍一眼看上去咋有机体系还能不收敛。。
（2）贴输入文件时最好直接上传文件，或用论坛插入代码功能贴出，这样空格之类的格式不会乱。你这样直接复制粘贴，有些格式消失了，也不知道你到底写对了没有。
（3）假设你写的输入文件是对的，就按#7楼给你的建议。

498746012 · 发表于 Post on 2021-4-1 14:11:47

zjxitcc 发表于 2021-4-1 14:01
（1）你这个例子与帖子一开始说的例子不是同一个体系，应该讲清楚。乍一眼看上去咋有机体系还能不收敛 ...

谢谢老师，感恩

498746012 · 发表于 Post on 2021-4-1 14:41:09

paramecium86 发表于 2021-4-1 13:35
SCF不容易收敛 GAMESS的老毛病了。
如果DIIS不行就换成SOSCF 如果soscf不行，就先DIIS再SOSCF。或者干 ...

学生是第一次接触gamess，老师您说的先DIIS再SOSCF, 是算完DIIS,读check文件接着算SOSCF,还是如下一个任务就搞定呀？$SCF DIIS=.t. SOSCF=.t. $END对吗

498746012 · 发表于 Post on 2021-4-1 15:07:30

zjxitcc 发表于 2021-4-1 14:01
（1）你这个例子与帖子一开始说的例子不是同一个体系，应该讲清楚。乍一眼看上去咋有机体系还能不收敛 ...

老师您好，学生想请教一下gamess可以做stable=opt的计算吗,有没有相关的关键词呀

paramecium86 · 发表于 Post on 2021-4-1 15:19:23

498746012 发表于 2021-4-1 14:41
学生是第一次接触gamess，老师您说的先DIIS再SOSCF, 是算完DIIS,读check文件接着算SOSCF,还是如下一个任 ...

不是。
如果想先做DIIS 在做 SOSCF大概写法如下
$scf dirscf=.true. fdiff=.f. soscf=.f.
      diis=.t. ethrsh=2.0 shift=.t.
      swdiis=0.001 damp=.t. $end

其中重要的是 swdiis 这个阈值。达到这个阈值之后就停止DIIS 开始soscf 而在设置上还是diis=.t.  soscf=.f.

gamess的scf里面里面选项很多。你最好每个都看看什么意思。不收敛的时候可能哪个就有用。不过我并不知道为什么有时候diis好用有时候 soscf好用。有时候就是先diis后soscf 好用。基本上这是我自己拿来做计算时候的经验之谈，自己琢磨着试出来的一些经验。

zjxitcc · 发表于 Post on 2021-4-1 15:22:29

498746012 发表于 2021-4-1 15:07
老师您好，学生想请教一下gamess可以做stable=opt的计算吗,有没有相关的关键词呀

有，但是功能有限，只能基于RHF/RDFT 检验RHF->UHF稳定性，不能基于UHF/UDFT 检验UHF的内部不稳定性（而高斯都支持）。看GAMESS手册里的$SCF章节，UHFCHK关键词。

你这一开始算的就是UB3LYP，在GAMESS里不支持检验。

498746012 · 发表于 Post on 2021-4-1 16:07:30

paramecium86 发表于 2021-4-1 15:19
不是。
如果想先做DIIS 在做 SOSCF大概写法如下
$scf dirscf=.true. fdiff=.f. soscf=.f.

谢谢老师，感恩

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GAMESS-US] 求助gamess计算开壳层单重态总是收敛到闭壳层