求助：Fe2O3吸附小分子机构优化波函数不收敛

alohomora100

求助大家，最近在用cp2k计算Fe2O3吸附小分子DMI的吸附能，目前优化Fe2O3 2*2的晶胞在z方向构建真空层构象的单点能，波函数是可以收敛的，输入输出文件在附件中（Fe2O3_22），但是在Fe2O3表面放一个DMI分子做结构优化时，波函数出现不收敛的现象，忽略不收敛做结构优化，优化几步之后也未见波函数出现收敛的情况，目前尝试了设置不同的k点，不同的CUTOFF和REL_CUTOFF效果都不理想，想求助大家有什么方法解决吗？输入输出文件也在附件中（Fe2O3_22_DMI_0421）

sobereva

盒子Z方向尺寸不合理
你用Multiwfn产生Fe2O3的二维周期性计算结合PSOLVER MT的输入文件时Multiwfn对当前体系恰到好处地定义盒子，满足MT的要求，然而你加了被吸附分子却不改盒子，盒子Z方向显然就不够大了。必须重新用Multiwfn创建输入文件。
CP2K用户一定要充分了解各种PSOLVER的特点。北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里“能量的计算及相关问题”部分讲“低维体系的计算”的一节非常详细讲了这些知识。