求助：cp2k处理磁性体系收敛极差的解决办法

Penson

输入文件见附件
如题，用cp2k做含有磁性原子LiCoO2的(001)切面计算时（已加z方向表面校正），使用ot做结构优化的能量收敛极差，请教各位大佬有什么比较快速高效的参数调节方法使收敛更快更方便？目前的输入文件并没有带入初始磁性，因为之前试过加了磁性收敛会差更多。

ChemG

cp2k对于有些磁性体系就是很难收敛，只能是用QE、vasp这些纯平面波的程序来算了。这种情况我出现了好机会，只能被迫忍受龟速的计算速度。

Penson

ChemG 发表于 2024-11-11 16:42
cp2k对于有些磁性体系就是很难收敛，只能是用QE、vasp这些纯平面波的程序来算了。这种情况我出现了好机会， ...

这样子吗，那真是头疼

卡开发发

Penson 发表于 2024-11-11 19:29
这样子吗，那真是头疼

1、你的coord.inc我下载后将你的晶胞参数加入改为GView的gjf格式，然后看了下似乎不是很正常，我不确定是你上传的结构没有对应还是这个结构建模的时候就是存在问题，截图如下。

2、LiCoO2我在Materials Project检索得到的磁性是无磁性，但仔细去看properties的部分给出来的能带是有alpha-和beta-band，很难说这个体系一定是有磁性的，或许需要更多的文献和计算上的尝试来确定这个问题。
3、这个体系在不考虑DFT+U与切面后可能存在表面态，或呈现金属性，OT按道理应当谨慎使用。

就收敛难易而言，只能说我个人经验上CP2K的数值稳定性不是很好，虽然经常吐槽，但这种问题很难进行量化讨论，不过在此之前或许我们应该先排除自身的一些不当操作。

Penson

卡开发发 发表于 2024-11-11 21:59
1、你的coord.inc我下载后将你的晶胞参数加入改为GView的gjf格式，然后看了下似乎不是很正常，我不确定是 ...

感谢大佬。
1.结构我是拿晶体结构切(001)之后加真空层再扩胞获得的。用MS导出.car文件再转的xyz。原子配比应该是合适，但其分布情况我还需要确认一下，之前确实没有太多考量这一块，以后一定会注意。
2.除了OT，对角化我也尝试过，收敛性会更差。
具体情况，我后面有时间再仔细验证和跟帖说明。

卡开发发

Penson 发表于 2024-11-12 12:17
感谢大佬。
1.结构我是拿晶体结构切(001)之后加真空层再扩胞获得的。用MS导出.car文件再转的xyz。原子配 ...

1、配比我检查过确实没问题，但位置我看上去存在问题。
2、标准对角化一般是收敛性更差，但OT的适用性需要做预判。

Penson

卡开发发 发表于 2024-11-12 13:14
1、配比我检查过确实没问题，但位置我看上去存在问题。
2、标准对角化一般是收敛性更差，但OT的适用性需 ...

好的，多谢