计算化学公社

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

[CP2K] cp2k-2026.2MIXING参数对收敛性的影响

[复制链接 Copy URL]

80

帖子

0

威望

759

eV
积分
839

Level 4 (黑子)

      学习CP2K已有一段时间,期间也同时接触了VASP。个人感觉,CP2K的上手难度确实高于VASP,此前也曾有过不少抱怨。所幸在Sob老师的指导下,我对CP2K的理解有了显著提升。
      以往,用CP2K计算磁性体系,尤其是施加k点的对角化算法时,收敛极其困难。CP2K 2026.2版本引入了HARRIS方法,允许读取Gamma点计算生成的波函数作为对角化计算的初始猜测,从而有效加速磁性体系的收敛。之所以先用**OT算法**提供初猜,是因为OT算法减少了待计算波函数的数量,在处理金属氧化物时,相比对角化算法更容易收敛。此外,该版本还新增了波函数混合方法及相关参数。
     近期,我利用CP2K‑2026.2的HARRIS方法,对我研究体系(掺杂型氢氧化镍晶胞,结构不一定正确)进行了系列SCF测试。具体流程为:先用OT算法获得波函数初猜,再读取该波函数,采用 MODIFIED_BROYDEN_MIXING 混合方法,并系统测量不同关键词的影响。测试分为两个部分:
       第一部分:依次测试了 ALPHA、ALPHA_MAG、BETA、BETA_MAG,每次均使用上一轮获得的最佳参数作为固定值。  
       第二部分:由于时间较晚,我采用控制变量法,保持其他参数为默认值(默认参数:NBUFFER=16,MAX_GVEC_EXP=3,BROY_W0=0.01,BROY_WMAX=100000.0,BROY_WREF=0.025,N_SIMPLE_MIX=16),依次测试上述变量。
      图1和图2展示了测试结果,横坐标为电子步数(每次SCF固定300步),纵坐标为收敛值(取常用对数)。  
      图1a 显示,ALPHA=0.2 时收敛值明显更低,表明该值最为合适。固定 ALPHA=0.2 后,测试 ALPHA_MAG(图1b),发现ALPHA_MAG=0.8时收敛更佳,更有利于加速收敛。随后固定这两项,测试 BETA,结果表明 BETA=1.5 效果更好;再测试 BETA_MAG,发现BETA_MAG=0 时最优(该部分结果因时间原因未等全部结束即记录)。

        第二部分测试(图2)主要结论如下:  
       1. NBUFFER 增大有助于收敛,其中 NBUFFER=28 时表现最佳。
       2. MAX_GVEC_EXP、BROY_W0和 BROY_WREF 对收敛影响不显著。
       3. BROY_WMAX 几乎无影响。
       4. N_SIMPLE_MIX 影响最大,它代表使用 KERKER_MIXING 进行初始混合的步数。当该值过大时,KERKER_MIXING本身即可直接收敛,无需再启用 MODIFIED_BROYDEN_MIXING。
       另外,若不采用HARRIS方法提供的波函数初猜,则无论何种混合方法,都很难在300步内收敛。
       最终总结:
      1. 对于磁性金属氧化物,可借助HARRIS方法,先用OT算法获得初始波函数,再读取该波函数进行后续计算。
      2. 在各种混合方法中,KERKER_MIXING的收敛性能最佳。
   
      测试结果图及原始数据文件已附上,供大家参考。本输入文件由Sob老师编写的Multiwfn生成,在此衷心感谢Sob老师的帮助!测试脚本由deepseek辅助完成,cp2k采用自己编译的cp2k2026.2。如果结果不对的,请各位老师提出批评,以及修改意见。最后一个问题:为什么现在论坛要将帖子上传的照片大小限制得那么小呢?图片经过较大程度的压缩,都不太清楚了。

第一步测试2.jpg (124.1 KB, 下载次数 Times of downloads: 0)

第一部分测试

第一部分测试

第2步测试2-1.jpg (110 KB, 下载次数 Times of downloads: 0)

第二部分测试结果1

第二部分测试结果1

第2步测试-3.jpg (60.14 KB, 下载次数 Times of downloads: 0)

第二部分测试结果2

第二部分测试结果2

cp2k.inp

12.57 KB, 下载次数 Times of downloads: 0

测试MIXING的参数

cp2k_harris.inp

13.62 KB, 下载次数 Times of downloads: 0

利用ENERGY_CORRECTION,采用OT算法提供对角化算法的初猜

测试文件.zip

1.63 MB, 下载次数 Times of downloads: 1

测试结果

6万

帖子

99

威望

6万

eV
积分
128792

管理员

公社社长

2#
发表于 Post on 1 hour ago | 只看该作者 Only view this author
为什么现在论坛要将帖子上传的照片大小限制得那么小呢?

论坛对jpg格式的图片上传限制是600KB、对png格式的限制是1000KB,无论如何对于所有情况都绝对完全够用了。只有那些没有图片尺寸、图片格式、图片压缩常识的人才会因为这个限制而导致上传的图片不够清楚
强烈建议有这种疑惑的人仔细把此文看了:
谈谈与计算化学有关的作图的图像清晰问题
http://sobereva.com/607http://bbs.keinsci.com/thread-24177-1-1.html

诸如你上传的图,用对于这种图最合适的无损压缩的png格式,图像尺寸保存得足够大比如800*600,能展示得十分清楚而且完全不会导致图片文件大小超标。你对图片处理仍缺乏必要的常识。本论坛对图片尺寸要求再合理不过,不仅完全够用,还杜绝了没有图片处理常识的人上传完全不必要的巨大图片文件(比如上传手机直接拍照的jpg原图一个好几兆、直接上传没有任何压缩的bmp格式,等等)导致严重浪费论坛服务器资源、带宽,而且使得访问者加载图片非常慢。


例如此图,png格式上传的,原图95KB,上传到论坛后自动加水印后也仅仅108KB,就已经十分清楚

北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-8-8 04:14 , Processed in 0.231846 second(s), 24 queries , Gzip On.

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