计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › eda结果出错

eda结果出错

查看数: 15997 | 评论数: 16 | 收藏 Add to favorites 2

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

zpk

发布时间: 2017-2-17 16:27

正文摘要:

请问一下，为什么EDA的结果中各能量值都特别大，正常情况下一般不都是几到十几吗？

回复 Reply

liusy 发表于 Post on 2019-4-2 21:48:17

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

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

zjxitcc 发表于 Post on 2019-4-2 19:29:11

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

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

liusy 发表于 Post on 2019-4-1 11:00:05

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

周子曦 发表于 Post on 2018-6-1 19:23:35

好的，谢谢老师

sobereva 发表于 Post on 2018-5-30 18:15:05

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

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

周子曦 发表于 Post on 2018-5-30 09:05:15

老师，您好，我也出现了跟楼主一样的情况，做能量分解单体一没有收敛，极化能特别高，我已经把迭代步骤加到了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

sobereva 发表于 Post on 2017-2-20 19:06:13

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

不用

zpk 发表于 Post on 2017-2-20 10:27:25

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

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

sobereva 发表于 Post on 2017-2-20 09:38:37

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

没有

zpk 发表于 Post on 2017-2-20 09:38:21

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

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

sobereva 发表于 Post on 2017-2-19 18:14:44

加上$DFT SWITCH=1.0 $END
如果还不行把没什么意义的氢上的弥散去掉

paramecium86 发表于 Post on 2017-2-19 15:02:25

GAMESS 的SCF默认迭代次数太少了才30 我一般默认全改成200步也就是程序上限.具体做法就是在$contrl里规定 maxit=200

zpk 发表于 Post on 2017-2-19 11:34:37

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

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

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