计算化学公社

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

[CP2K] OT开OUTER SCF在即将收敛时收敛缓慢

[复制链接 Copy URL]

361

帖子

0

威望

1987

eV
积分
2348

Level 5 (御坂)

所念皆星河

本帖最后由 含光君 于 2024-3-3 16:43 编辑

各位老师好,我最近在计算Pd纳米颗粒负载于石墨烯的体系时,发现开了OT和OUTER SCF后,在Convergence收敛至0.00000200左右时开始收敛得非常缓慢:
  1.   Step     Update method      Time    Convergence         Total energy    Change
  2.   ------------------------------------------------------------------------------
  3.      1 OT DIIS     0.15E+00   13.3     0.00000190    -16401.8884893442 -9.46E-05
  4.      2 OT DIIS     0.15E+00    5.9     0.00000190    -16401.8884951057 -5.76E-06
  5.      3 OT DIIS     0.15E+00    5.9     0.00000188    -16401.8884978986 -2.79E-06
  6.      4 OT DIIS     0.15E+00    5.9     0.00000188    -16401.8885196240 -2.17E-05
  7.      5 OT DIIS     0.15E+00    5.9     0.00000186    -16401.8885421686 -2.25E-05
  8.      6 OT DIIS     0.15E+00    5.9     0.00000186    -16401.8885582783 -1.61E-05
  9.      7 OT DIIS     0.15E+00    6.0     0.00000185    -16401.8886189631 -6.07E-05
  10.      8 OT DIIS     0.15E+00    6.6     0.00000185    -16401.8886542629 -3.53E-05
  11.      9 OT DIIS     0.15E+00    6.0     0.00000185    -16401.8886870421 -3.28E-05
  12.     10 OT DIIS     0.15E+00    6.0     0.00000184    -16401.8887392789 -5.22E-05
  13.     11 OT DIIS     0.15E+00    5.9     0.00000183    -16401.8888117160 -7.24E-05
复制代码
如果按这样的速度收敛的话,极有可能要上千圈才能达到收敛标准。不知道是我输入文件的参数哪边有问题还是OT本来就会这样?

输入文件、结构文件和输出文件已上传至附件,恳请各位老师赐教!
g-Pd116.pdb (40.16 KB, 下载次数 Times of downloads: 2) g-Pd116.inp (4.01 KB, 下载次数 Times of downloads: 10) cp2k.out (94.59 KB, 下载次数 Times of downloads: 3)




心之所向,日复一日,必有精进

361

帖子

0

威望

1987

eV
积分
2348

Level 5 (御坂)

所念皆星河

2#
 楼主 Author| 发表于 Post on 2024-3-3 16:53:34 | 只看该作者 Only view this author
似乎第一圈结构优化没收敛后,第二圈结构优化开始变得可以收敛。在第二圈结构优化中,虽然即将收敛时Convergence变化依然很慢,但会在开始几步快速收敛至0.00000100附近,因此好像不用管第一步了...
心之所向,日复一日,必有精进

40

帖子

1

威望

732

eV
积分
792

Level 4 (黑子)

3#
发表于 Post on 2024-3-4 06:45:38 | 只看该作者 Only view this author
这个体系感觉不是很适合用OT,换对角化试试

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

4#
发表于 Post on 2024-3-4 07:23:56 | 只看该作者 Only view this author
无gap或gap极小的情况应当默认用对角化结合smearing
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

361

帖子

0

威望

1987

eV
积分
2348

Level 5 (御坂)

所念皆星河

5#
 楼主 Author| 发表于 Post on 2024-3-4 15:11:51 | 只看该作者 Only view this author
ultramanm87 发表于 2024-3-4 06:45
这个体系感觉不是很适合用OT,换对角化试试

感谢老哥提醒!学习啦!
心之所向,日复一日,必有精进

361

帖子

0

威望

1987

eV
积分
2348

Level 5 (御坂)

所念皆星河

6#
 楼主 Author| 发表于 Post on 2024-3-4 15:13:13 | 只看该作者 Only view this author
sobereva 发表于 2024-3-4 07:23
无gap或gap极小的情况应当默认用对角化结合smearing

明白啦!非常感谢sob老师!
心之所向,日复一日,必有精进

48

帖子

0

威望

662

eV
积分
710

Level 4 (黑子)

7#
发表于 Post on 2024-5-25 10:15:49 | 只看该作者 Only view this author
sobereva 发表于 2024-3-4 07:23
无gap或gap极小的情况应当默认用对角化结合smearing

sob老师,想问下无gap或gap极小的体系在结构优化时,如果OT能比开smearing的对角化方法更快收敛,是不是可以用OT

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

8#
发表于 Post on 2024-5-25 10:51:31 | 只看该作者 Only view this author
gejuqing 发表于 2024-5-25 10:15
sob老师,想问下无gap或gap极小的体系在结构优化时,如果OT能比开smearing的对角化方法更快收敛,是不是 ...

对于不考虑k点时可以,OT不支持k点
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 15:00 , Processed in 0.168123 second(s), 24 queries , Gzip On.

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