计算化学公社

标题: CP2K计算Ni单点能不收敛 [打印本页]

作者
Author: zhifanfan 时间: 2025-4-15 09:21
标题: CP2K计算Ni单点能不收敛
各位老师好，我在使用CP2K计算周期性Ni体系的单点能时，就是不收敛，现附上输入文件，请各位老师帮忙检查错误，非常感谢

作者
Author: snljty2 时间: 2025-4-15 10:18
金属Ni单质是铁磁性的，先尝试用MAGNETIZATION设置合理的初猜磁矩。

作者
Author: zhifanfan 时间: 2025-4-16 11:15

snljty2 发表于 2025-4-15 10:18
金属Ni单质是铁磁性的，先尝试用MAGNETIZATION设置合理的初猜磁矩。

好的，谢谢您，请问用multiwfn就能直接设置吗，这个不太会

作者
Author: qczcb2 时间: 2025-4-17 09:10
本帖最后由 qczcb2 于 2025-4-17 09:14 编辑

zhifanfan 发表于 2025-4-16 11:15
好的，谢谢您，请问用multiwfn就能直接设置吗，这个不太会

能是能，但是cp2k算Ni金属就是不容易收敛，设置磁矩也没用。建议换qe或者vasp。如果需要设置磁矩，在multiwfn里cp2k菜单主界面选-9 9 进入磁矩设置，选择*Ni 1 设置所有的Ni磁矩为1就行了，这是铁磁性的。设置反铁磁需要设置完磁矩后，选择需要设置为相反磁矩的原子号，在上一步基础上选择原子序号如1-10，再设置-1磁矩。另外，这个Ni晶胞大小为10 A左右，不能用gamma，应该用3*3*3的k点

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