计算化学公社

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

[新手求助] 几何优化收敛与SCF收敛

[复制链接 Copy URL]

40

帖子

0

威望

373

eV
积分
413

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
大家好,有个问题想请教大家,怎么区分是几何优化收敛和SCF收敛啊?优化收敛是不是在构型优化过程中遇到的啊即优化输出文件中的那四个收敛条件?而SCF收敛是指在单点能计算过程中遇到的啊?谢谢大家

1102

帖子

18

威望

6693

eV
积分
8155

Level 6 (一方通行)

計算化学の社畜

2#
发表于 Post on 2016-5-10 19:25:11 | 只看该作者 Only view this author
简单的来说,SCF收敛是几何优化过程中【每一步】几何优化的一次循环(以RMSDP、MaxDP和DE为收敛限);而几何优化其实相当于大循环(以Maximum Force 、RMS  Force 、Maximum Displacement 和RMS  Displacement 四个几何收敛判据为收敛限)

评分 Rate

参与人数
Participants 3
eV +7 收起 理由
Reason
gezq + 1
Henge + 5
sobereva + 1

查看全部评分 View all ratings

Stand on the shoulders of giants

40

帖子

0

威望

373

eV
积分
413

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2016-5-16 09:07:01 | 只看该作者 Only view this author
冰释之川 发表于 2016-5-10 19:25
简单的来说,SCF收敛是几何优化过程中【每一步】几何优化的一次循环(以RMSDP、MaxDP和DE为收敛限);而几 ...

好的 谢谢

40

帖子

0

威望

373

eV
积分
413

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2016-5-16 10:44:34 | 只看该作者 Only view this author
冰释之川 发表于 2016-5-10 19:25
简单的来说,SCF收敛是几何优化过程中【每一步】几何优化的一次循环(以RMSDP、MaxDP和DE为收敛限);而几 ...

那我接着请教一下  那我该如何从我的输出文件中判断它是scf不收敛还是优化不收敛呢 这是我的输出文件中的最后 他只告诉我收敛失败
Cycle 128  Pass 1  IDiag  1:
RMSU=  3.79D-06    CP:  1.00D+00  3.00D+00  1.34D+00  1.53D+00  1.55D+00
E= -3413.82812695882     Delta-E=       -0.000009770827 Rises=F Damp=F
DIIS: error= 1.41D-04 at cycle 125 NSaved=  20.
NSaved=20 IEnMin=20 EnMin= -3413.82812695882     IErMin=20 ErrMin= 1.41D-04
ErrMax= 1.41D-04 EMaxC= 1.00D-01 BMatC= 4.98D-06 BMatP= 4.55D-06
IDIUse=3 WtCom= 9.99D-01 WtEn= 1.41D-03
Matrix for removal  1 Erem= -3413.82695539264     Crem= 0.000D+00
Coeff-Com:  0.152D+00-0.116D+00 0.150D+00 0.688D+00-0.113D+00-0.289D+00
Coeff-Com: -0.350D-01-0.114D+01 0.360D+00-0.291D+00 0.294D+00-0.821D-01
Coeff-Com: -0.135D+00 0.390D+00 0.632D+00 0.181D+00 0.124D+00-0.499D-02
Coeff-Com:  0.809D-01 0.156D+00
Coeff-En:   0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En:   0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En:   0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En:   0.000D+00 0.100D+01
Coeff:      0.152D+00-0.116D+00 0.150D+00 0.687D+00-0.113D+00-0.289D+00
Coeff:     -0.350D-01-0.114D+01 0.359D+00-0.291D+00 0.294D+00-0.820D-01
Coeff:     -0.135D+00 0.390D+00 0.631D+00 0.181D+00 0.123D+00-0.498D-02
Coeff:      0.808D-01 0.157D+00
Gap=     0.007 Goal=   None    Shift=    0.000
RMSDP=1.22D-05 MaxDP=9.60D-04 DE=-9.77D-06 OVMax= 8.08D-03

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RPBE1PBE) =  -3413.82812696     A.U. after  129 cycles
             Convg  =    0.1223D-04             -V/T =  2.0183
KE= 3.352616297113D+03 PE=-2.545790568027D+04 EE= 9.907571306294D+03
Convergence failure -- run terminated.
Error termination via Lnk1e in /share/apps/gauss/g09/l502.exe at Mon May 16 02:49:09 2016.
Job cpu time:  2 days  6 hours 19 minutes 48.7 seconds.
File lengths (MBytes):  RWF=    524 Int=      0 D2E=      0 Chk=     57 Scr=      1

1102

帖子

18

威望

6693

eV
积分
8155

Level 6 (一方通行)

計算化学の社畜

5#
发表于 Post on 2016-5-16 10:50:59 | 只看该作者 Only view this author
最瘦二师兄 发表于 2016-5-16 10:44
那我接着请教一下  那我该如何从我的输出文件中判断它是scf不收敛还是优化不收敛呢 这是我的输出文件中的 ...

SCF不收敛

解决办法参考Sob老师的博文:http://www.sobereva.com/61
Stand on the shoulders of giants

40

帖子

0

威望

373

eV
积分
413

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2016-5-16 14:04:05 | 只看该作者 Only view this author
本帖最后由 最瘦二师兄 于 2016-5-16 14:10 编辑
冰释之川 发表于 2016-5-16 10:50
SCF不收敛

解决办法参考Sob老师的博文:http://www.sobereva.com/61

请问是根据这个RMSDP=1.22D-05 MaxDP=9.60D-04 DE=-9.77D-06 OVMax= 8.08D-03

>>>>>>>>>> Convergence criterion not met.因为RMSDP,MaxDP,DE是scf的收敛限而 Convergence criterion not met,从而判断出来它是scf不收敛吗

6万

帖子

99

威望

5万

eV
积分
120141

管理员

公社社长

7#
发表于 Post on 2016-5-16 16:11:11 | 只看该作者 Only view this author
最瘦二师兄 发表于 2016-5-16 14:04
请问是根据这个RMSDP=1.22D-05 MaxDP=9.60D-04 DE=-9.77D-06 OVMax= 8.08D-03

>>>>>>>>>> ...

这提示一看就是SCF步骤的输出,稍有使用经验就能判断。你说的那种判断也对。
北京科音自然科学研究中心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

40

帖子

0

威望

373

eV
积分
413

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2016-5-16 19:44:53 | 只看该作者 Only view this author
sobereva 发表于 2016-5-16 16:11
这提示一看就是SCF步骤的输出,稍有使用经验就能判断。你说的那种判断也对。

明白了 谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 22:12 , Processed in 0.147601 second(s), 21 queries , Gzip On.

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