计算化学公社

标题: 求助gamess计算开壳层单重态总是收敛到闭壳层 [打印本页]

作者
Author: 498746012 时间: 2021-4-1 09:45
标题: 求助gamess计算开壳层单重态总是收敛到闭壳层
各位老师打了，用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

作者
Author: zjxitcc 时间: 2021-4-1 09:59
本帖最后由 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参数了么？

作者
Author: 498746012 时间: 2021-4-1 10:10

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

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

作者
Author: zjxitcc 时间: 2021-4-1 10:21

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

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

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

作者
Author: 498746012 时间: 2021-4-1 10:41

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

谢谢老师，太感谢了

作者
Author: 498746012 时间: 2021-4-1 12:36
本帖最后由 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

作者
Author: paramecium86 时间: 2021-4-1 13:35

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各种组合都试试

作者
Author: 498746012 时间: 2021-4-1 13:54

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

谢谢老师，太感谢了

作者
Author: zjxitcc 时间: 2021-4-1 14:01

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

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

作者
Author: 498746012 时间: 2021-4-1 14:11

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

谢谢老师，感恩

作者
Author: 498746012 时间: 2021-4-1 14:41

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对吗

作者
Author: 498746012 时间: 2021-4-1 15:07

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

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

作者
Author: paramecium86 时间: 2021-4-1 15:19

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 好用。基本上这是我自己拿来做计算时候的经验之谈，自己琢磨着试出来的一些经验。

作者
Author: zjxitcc 时间: 2021-4-1 15:22

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

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

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

作者
Author: 498746012 时间: 2021-4-1 16:07

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

谢谢老师，感恩

作者
Author: 498746012 时间: 2021-4-1 16:14

zjxitcc 发表于 2021-4-1 15:22
有，但是功能有限，只能基于RHF/RDFT 检验RHF->UHF稳定性，不能基于UHF/UDFT 检验UHF的内部不稳定性（而 ...

老师打扰了，刚才我没有说清楚，我就是想基于RHF/RDFT 检验RHF->UHF稳定性，从而得到开壳层单重态

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