计算化学公社

标题: 晶胞优化收敛后再通过MS设置真空层后opt不收敛 [打印本页]

作者
Author: JianZ 时间: 2024-4-23 16:47
标题: 晶胞优化收敛后再通过MS设置真空层后opt不收敛
刚开始学使用cp2k这个软件，采用Fe的fcc结构替换两个Cr原子后进行晶胞优化结果收敛。目的是作为一个吸附的基底，所以再通过MS设置真空层再进行优化结果再也收敛不了，排除问题后去除吸附原子只留下优化后的基底后发现不收敛，实在理解不了这个情况，向各位大佬求解一下。
fefcc是进行晶胞优化的文件，fedzvp是设置真空层后的文件

作者
Author: sobereva 时间: 2024-4-24 07:05
显然得开smearing。当前当前都没设初猜磁矩这些，输入文件也太简单了

参考下文解决
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

d族金属与真空的界面体系本身就是属于较难收敛的一类，表面原子的电子结构比体相更复杂。

作者
Author: JianZ 时间: 2024-5-7 10:03
还有一个用的GFN1-xtb方法收敛，但是结果Cs原子跑到我设的盒子边缘，正常应该是在基底中央上方的，是不是因为GFN1算小体系不准的原因。

作者
Author: JianZ 时间: 2024-5-7 14:54

sobereva 发表于 2024-4-24 07:05
显然得开smearing。当前当前都没设初猜磁矩这些，输入文件也太简单了

参考下文解决

sobereva老师根据您提到的设置磁矩重新设置inp文件，将fe和cr的初始磁矩设为3，通过smearing来使得迭代过程中每个原子的磁矩会慢慢趋于合理。但优化结果还是有很大问题，出现大量原子跑偏。想请教下您输入文件是否有其他问题。

作者
Author: sobereva 时间: 2024-5-8 20:40

JianZ 发表于 2024-5-7 14:54
sobereva老师根据您提到的设置磁矩重新设置inp文件，将fe和cr的初始磁矩设为3，通过smearing来使得迭代过 ...

除了z方向不该设k点外，没有明显问题
注意监控优化过程中SCF收敛情况，严重没收敛可能导致受力离谱

作者
Author: JianZ 时间: 2024-5-9 18:43

sobereva 发表于 2024-5-8 20:40
除了z方向不该设k点外，没有明显问题
注意监控优化过程中SCF收敛情况，严重没收敛可能导致受力离谱

谢谢老师解答，确实出现了严重不收敛受力离谱的情况
我推测是我在第一步晶胞收敛时没有加磁性设置导致后来加入真空层后不收敛，所以先进行加磁矩使得晶胞收敛再进行下一步。
对于您提出的k点设置是否可以理解为因为加入真空层后z方向太大，所以z方向上k点应该设置为1。
但是在晶胞优化时出现

*** WARNING in qs_scf_post_gpw.F:1848 :: Spin contamination estimate not ***
*** implemented for k-points. ***

这个错误是否也和k点设置有关

作者
Author: sobereva 时间: 2024-5-10 07:49

JianZ 发表于 2024-5-9 18:43
谢谢老师解答，确实出现了严重不收敛受力离谱的情况
我推测是我在第一步晶胞收敛时没有加磁性设置导致后 ...

非周期性方向上设k点毫无意义

根本不是报错。参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）的ppt：
(, 下载次数 Times of downloads: 9)

欢迎光临计算化学公社 (http://bbs.keinsci.com/)