计算化学公社

标题: eda结果出错 [打印本页]

作者
Author: zpk 时间: 2017-2-17 16:27
标题: eda结果出错
(, 下载次数 Times of downloads: 62)
请问一下，为什么EDA的结果中各能量值都特别大，正常情况下一般不都是几到十几吗？

作者
Author: sobereva 时间: 2017-2-17 17:12
肯定输入文件有问题，或者计算过程出现问题（检查整个计算中各个能量项输出看问题从哪里开始）

作者
Author: paramecium86 时间: 2017-2-17 21:55
看起来是某一个或者几个monomer 或者 complex 没有converge。 GAMESS在计算解离能的时候没有converge的那一项就认为是0，少了这部分能量所以算完解离能非常的负如果这个能量和你以往算某一个monomer的能量很接近就很容易知道哪个monomer没有converge 回到ouput文件里面翻翻就能确认了。

作者
Author: zpk 时间: 2017-2-19 11:34

paramecium86 发表于 2017-2-17 21:55
看起来是某一个或者几个monomer 或者 complex 没有converge。 GAMESS在计算解离能的时候没有converge的那 ...

谢谢。我在out文件中查到是monomer1没收敛。
输入文件：
$contrl scftyp=rhf runtyp=eda dfttyp=m06-2x nprint=-5 $end
$contrl icharg=0 mult=1 $end
$basis gbasis=n311 ngauss=6 npfunc=1 ndfunc=1
  diffsp=.t. diffs=.t. $end
$lmoeda matom(1)=13,3 mcharg(1)=0,0 mmult(1)=1,1 supbas=.t.
dftxc=.t. usevec=.t.  $end
$system mwords=3000 memddi=6000 $end

SCF迭代部分：
ITER EX DEM    TOTAL ENERGY       E CHANGE  DENSITY CHANGE DIIS ERROR
1  0  0    -304.1693850393  -304.1693850393 7.014814364 0.880549243
      * * * INITIATING DIIS PROCEDURE * * *
2  1  0    -300.5306954648    3.6386895745  51.535722464 0.485255491
3  2  0    -276.0853503417 24.4453451231  50.918211084 0.995947077
4  3  0    -300.6428371077 -24.5574867659 3.869036329 0.481414656
5  4  0    -303.8310207099 -3.1881836022 1.105190294 0.381490392
6  5  0    -304.8395126981 -1.0084919882 1.097061795 0.197149910
7  6  0    -305.1444809047 -0.3049682066 0.699773819 0.093085327
8  7  0    -305.2543134334 -0.1098325287 1.553695733 0.049831020
9  8  0    -305.2859683361 -0.0316549027 0.543785856 0.012347135
  10  9  0    -305.2874425542 -0.0014742181 0.076968418 0.003952934
  11 10  0    -305.2880087786 -0.0005662244 0.062400458 0.004262990
  12 11  0    -305.2887675244 -0.0007587457 0.035984657 0.004511638
  13 12  0    -305.2891876053 -0.0004200809 0.120205153 0.004660689
  14 13  0    -305.2880793417    0.0011082636 0.267093253 0.003616527
  15 14  0    -305.2859388286    0.0021405131 0.054903263 0.007577236
  16 15  0    -305.2861268932 -0.0001880646 0.103227825 0.006985524
  17 16  0    -305.2864196977 -0.0002928045 0.200592569 0.007556871
  18 17  0    -305.2862507311    0.0001689666 0.248641449 0.006327123
  19 18  0    -305.2870110993 -0.0007603682 0.121889054 0.002015137
  20 19  0    -305.2871563808 -0.0001452815 0.012823916 0.001196449
  21 20  0    -305.2871196449    0.0000367359 0.034746596 0.001782930
  22 21  0    -305.2870438454    0.0000757995 0.059881067 0.002798897
  23 22  0    -305.2871680797 -0.0001242344 0.047512482 0.001181841
  24 23  0    -305.2872273642 -0.0000592845 0.006117008 0.001351479
  25 24  0    -305.2872334639 -0.0000060996 0.005208208 0.001330580
  26 25  0    -305.2872429365 -0.0000094726 0.021413785 0.001200927
  27 26  0    -305.2872990087 -0.0000560722 0.008825007 0.001188627
  28 27  0    -305.2872465204    0.0000524883 0.002448868 0.001139170
CONVERGED TO SWOFF, SO DFT CALCULATION IS NOW SWITCHED ON.
      * * * INITIATING DIIS PROCEDURE * * *
  29 28  0    -307.0500975709 -1.7628510505 3.094103061 0.095739440
  30 29  0    -306.1962673558    0.8538302151 1.571677233 0.220837628

SCF IS UNCONVERGED, TOO MANY ITERATIONS

请问输入文件应该如何修改

作者
Author: paramecium86 时间: 2017-2-19 15:02
GAMESS 的SCF默认迭代次数太少了才30 我一般默认全改成200步也就是程序上限.具体做法就是在$contrl里规定 maxit=200

作者
Author: zpk 时间: 2017-2-19 16:31

paramecium86 发表于 2017-2-19 15:02
GAMESS 的SCF默认迭代次数太少了才30 我一般默认全改成200步也就是程序上限.具体做法就是在$contrl里规 ...

我把SCF迭代次数调到200，结果还是没有收敛

作者
Author: sobereva 时间: 2017-2-19 18:14
加上$DFT SWITCH=1.0 $END
如果还不行把没什么意义的氢上的弥散去掉

作者
Author: zpk 时间: 2017-2-20 09:38

sobereva 发表于 2017-2-19 18:14
加上$DFT SWITCH=1.0 $END
如果还不行把没什么意义的氢上的弥散去掉

老师你好，我在你的博文中看到”对于使用了弥散函数的情况出现不收敛时建议先尝试SCF=NoVarAcc来避免因自动降低积分精度而导致的这个问题”，请问在gamess有相应的选项吗？

作者
Author: sobereva 时间: 2017-2-20 09:38

zpk 发表于 2017-2-20 09:38
老师你好，我在你的博文中看到”对于使用了弥散函数的情况出现不收敛时建议先尝试SCF=NoVarAcc来避免因自 ...

没有

作者
Author: zpk 时间: 2017-2-20 10:27

sobereva 发表于 2017-2-20 09:38
没有

gamess手册上说：It is reasonable to pick SWOFF > SWITCH > CONV in $SCF，如果把switch调成1， SWOFF 要调吗？

作者
Author: sobereva 时间: 2017-2-20 19:06

zpk 发表于 2017-2-20 10:27
gamess手册上说：It is reasonable to pick SWOFF > SWITCH > CONV in $SCF，如果把switch调成1， SWOFF ...

不用

作者
Author: 周子曦 时间: 2018-5-30 09:05
老师，您好，我也出现了跟楼主一样的情况，做能量分解单体一没有收敛，极化能特别高，我已经把迭代步骤加到了200，也试了$DFT SWITCH=1.0 $END还是没能收敛，以下是我的输入文件，请问老师这种情况该怎么处理
$contrl scftyp=uhf runtyp=eda dfttyp=B3LYP coord=UNIQUE icharg=0 mult=2
  nosym=1 ispher=1 maxit=200 $end
$guess  guess=huckel $end
$system mwords=125 memddi=2000 $end
$basis gbasis=ACCT $end
$SCF DIIS=.t. SOSCF=.F. DIRSCF=.T. FDIFF=.F. NPUNCH=1 CONV=1e-7 SHIFT=.T. $END
$DFT SWITCH=1.0 $END
$lmoeda matom(1)=5,2 mcharg(1)=0,0 mmult(1)=2,1 supbas=.true. $end
$data

作者
Author: sobereva 时间: 2018-5-30 18:15

周子曦发表于 2018-5-30 09:05
老师，您好，我也出现了跟楼主一样的情况，做能量分解单体一没有收敛，极化能特别高，我已经把迭代步骤加到 ...

基组改小点，改用月份基组，比如may-cc-pVTZ，砍掉氢的弥散函数和高角动量弥散，弥散越多越难收敛
谈谈弥散函数和“月份”基组
http://sobereva.com/119

作者
Author: 周子曦 时间: 2018-6-1 19:23
好的，谢谢老师

作者
Author: liusy 时间: 2019-4-1 11:00

sobereva 发表于 2018-5-30 18:15
基组改小点，改用月份基组，比如may-cc-pVTZ，砍掉氢的弥散函数和高角动量弥散，弥散越多越难收敛
谈 ...

老师，您好，我也是出现这种无法收敛的情况，是否aug-cc-PVDZ的基组也不适用，以下是我的输入文件，请教一下如何处理
$contrl dfttyp=b3lyp scftyp=rhf runtyp=eda nosym=1 maxit=100 $end
$basis basnam(1)=ligh,ligo,ligh,ligh,lign,ligh,ligh $end
$lmoeda matom(1)=3,4 mcharg(1)=0,0 mmult(1)=1,1 $end
$system mwords=1000 $end
$SCF DIIS=.T. SOSCF=.F. DIRSCF=.T. FDIFF=.F. NPUNCH=1 $END
$DFT SWITCH=1.0 $END
$DATA
nh3_h2o_1 by b3lyp/aug-cc-PVDZ

作者
Author: zjxitcc 时间: 2019-4-2 19:29

liusy 发表于 2019-4-1 11:00
老师，您好，我也是出现这种无法收敛的情况，是否aug-cc-PVDZ的基组也不适用，以下是我的输入文件，请教 ...

没什么百试百灵的办法。对泛函不讲究的话，可以试试wB97X之类的，总有一个能收敛。

作者
Author: liusy 时间: 2019-4-2 21:48

zjxitcc 发表于 2019-4-2 19:29
没什么百试百灵的办法。对泛函不讲究的话，可以试试wB97X之类的，总有一个能收敛。

好的，我试试，非常感谢您，总是您在帮我解答

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