计算化学公社

 找回密码 Forget password
 注册 Register
Views: 15996|回复 Reply: 16
打印 Print 上一主题 Last thread 下一主题 Next thread

[GAMESS-US] eda结果出错

[复制链接 Copy URL]

16

帖子

0

威望

3951

eV
积分
3967

Level 5 (御坂)


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

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

2#
发表于 Post on 2017-2-17 17:12:42 | 只看该作者 Only view this author
肯定输入文件有问题,或者计算过程出现问题(检查整个计算中各个能量项输出看问题从哪里开始)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

328

帖子

1

威望

5841

eV
积分
6189

Level 6 (一方通行)

3#
发表于 Post on 2017-2-17 21:55:25 | 只看该作者 Only view this author
看起来是某一个或者几个monomer 或者 complex 没有converge。 GAMESS在计算解离能的时候 没有converge的那一项就认为是0, 少了这部分能量 所以算完解离能非常的负 如果这个能量和你以往算某一个monomer的能量很接近 就很容易知道哪个monomer没有converge 回到ouput文件里面翻翻就能确认了。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

16

帖子

0

威望

3951

eV
积分
3967

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2017-2-19 11:34:37 | 只看该作者 Only view this author
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

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

328

帖子

1

威望

5841

eV
积分
6189

Level 6 (一方通行)

5#
发表于 Post on 2017-2-19 15:02:25 | 只看该作者 Only view this author
GAMESS 的SCF默认迭代次数太少了 才30 我一般默认全改成200步 也就是程序上限.具体做法就是在$contrl里 规定 maxit=200

16

帖子

0

威望

3951

eV
积分
3967

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2017-2-19 16:31:23 | 只看该作者 Only view this author
paramecium86 发表于 2017-2-19 15:02
GAMESS 的SCF默认迭代次数太少了 才30 我一般默认全改成200步 也就是程序上限.具体做法就是在$contrl里 规 ...

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

360截图20170219163006691.jpg (40.21 KB, 下载次数 Times of downloads: 65)

360截图20170219163006691.jpg

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

7#
发表于 Post on 2017-2-19 18:14:44 | 只看该作者 Only view this author
加上$DFT SWITCH=1.0 $END
如果还不行把没什么意义的氢上的弥散去掉
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

16

帖子

0

威望

3951

eV
积分
3967

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2017-2-20 09:38:21 | 只看该作者 Only view this author
sobereva 发表于 2017-2-19 18:14
加上$DFT SWITCH=1.0 $END
如果还不行把没什么意义的氢上的弥散去掉

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

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

9#
发表于 Post on 2017-2-20 09:38:37 | 只看该作者 Only view this author
zpk 发表于 2017-2-20 09:38
老师你好,我在你的博文中看到”对于使用了弥散函数的情况出现不收敛时建议先尝试SCF=NoVarAcc来避免因自 ...

没有
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

16

帖子

0

威望

3951

eV
积分
3967

Level 5 (御坂)

10#
 楼主 Author| 发表于 Post on 2017-2-20 10:27:25 | 只看该作者 Only view this author

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

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

11#
发表于 Post on 2017-2-20 19:06:13 | 只看该作者 Only view this author
zpk 发表于 2017-2-20 10:27
gamess手册上说:It is reasonable to pick SWOFF > SWITCH > CONV in $SCF,如果把switch调成1, SWOFF  ...

不用
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

5

帖子

0

威望

789

eV
积分
794

Level 4 (黑子)

12#
发表于 Post on 2018-5-30 09:05:15 | 只看该作者 Only view this author
老师,您好,我也出现了跟楼主一样的情况,做能量分解单体一没有收敛,极化能特别高,我已经把迭代步骤加到了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   

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

13#
发表于 Post on 2018-5-30 18:15:05 | 只看该作者 Only view this author
周子曦 发表于 2018-5-30 09:05
老师,您好,我也出现了跟楼主一样的情况,做能量分解单体一没有收敛,极化能特别高,我已经把迭代步骤加到 ...


基组改小点,改用月份基组,比如may-cc-pVTZ,砍掉氢的弥散函数和高角动量弥散,弥散越多越难收敛
谈谈弥散函数和“月份”基组
http://sobereva.com/119
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

5

帖子

0

威望

789

eV
积分
794

Level 4 (黑子)

14#
发表于 Post on 2018-6-1 19:23:35 | 只看该作者 Only view this author
好的,谢谢老师

26

帖子

0

威望

121

eV
积分
147

Level 2 能力者

15#
发表于 Post on 2019-4-1 11:00:05 | 只看该作者 Only view this author
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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-16 17:11 , Processed in 0.296918 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list