宸牧羊 发表于 2024-12-24 11:29 大约相当于scf=conver=6 算能量凑凑合合，如果是涉及到能量导数的任务就太糙了 |
本帖最后由 宸牧羊 于 2024-12-24 11:36 编辑 sobereva 发表于 2024-12-11 02:24 参照帖子 CP2K中遇到SCF难收敛时的解决方法 http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html） 我尝试增加了K点数量和cutoff后，出现了收敛趋势，优化过程中convergence到达0.00000几后一直波动，一直无法达到1*10-6的收敛限（此时我尝试适当的加大MAX_SCF也不奏效）， 于是我把收敛限改为1*10-5后重新跑，就可以很快收敛了，请问sob老师这样的结果可信度高么。 提出这样问题的原因是，我之前参加量子化学基础班的时候，讲义上好像说几何优化的收敛标准应当适当的高一些；对于cp2k程序，1*10-5的scf收敛限属于什么精度（勉强接受，还是很粗糙）。 |
宸牧羊 发表于 2024-12-10 11:27 有的情况通过调整SCF相关设置，虽然可能收敛慢，但最后还是能收敛。有的情况就是各种办法都弄不收敛。实际试试便知 XY周期性本身没问题。关键差异在于slab还是体相，体相比slab容易收敛得多，因为表面部分的原子电子结构比体相明显更复杂，导致SCF更难收敛 |
sobereva 发表于 2024-12-10 07:57 谢谢sob老师的回答，关于d族金属难收敛的问题，是指需要调整合适的参数后SCF才能够收敛，还是指SCF一开始无法收敛，需要计算较长时间后才会出现收敛？ 另外，关于现在的模型是只考虑xy周期性的，这对体系的收敛会有影响么。（我对Co的原胞计算xyz周期性，很快就可以收敛，但是一旦扩胞，加上周期性限制，就无法收敛） |
|
CP2K中遇到SCF难收敛时的解决方法 http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html） 上文该说的都说了，没有可补充的。CP2K算d族纯金属slab体系往往就是很难收敛，下文也说了，实在搞不定就用SIRIUS模块基于平面波计算或者QE 驳网上流传的对CP2K缺点的不实描述 http://sobereva.com/729（http://bbs.keinsci.com/thread-50094-1-1.html） 另外，磁性金属slab表面和内部的金属的自旋磁矩是不同的，统一用相同的值不是合适做法。北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里专门提了这点： 
|
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-17 03:43 , Processed in 0.196784 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites