|
|
请教各位有QE使用经验的老师，最近学生刚刚接触第一性原理的计算，需要计算一个金属有机化合物的反铁磁态，晶胞中含有二茂钴和一个有机物（标记为A），A作电子受体（也许可以理解为是A和二茂钴的共晶？）。
由于对第一性原理的计算没有经验，所以我首先抠出了晶胞中的分子团簇，这一团簇包含四个二茂钴和四个A，用高斯16优化了其中氢原子的位置，而其他重原子的位置没有优化，优化结束后得到了体系的对称破缺波函数解，从自旋密度等值面图看出，自旋密度在各二茂钴的钴中心上均有分布，且符号均为正，而在四个A上也有自旋密度的分布，且符号为负，但是自旋密度是离域地分布在A上的，并没有局域于A的某个原子上，电子结构是反铁磁态。基于这一结果，我开始使用QE进行第一性原理计算。
我根据丁老师的帖子（http://bbs.keinsci.com/thread-31439-1-1.html）建立了输入文件，赝势直接使用了Standard solid-state pseudopotentials的SSSP Efficiency，由于目前处于截断能的收敛测试，所以设置的ecutwfc为12Ry。设定的starting_magnetization都在四个钴上，其中两个为0.5，另外两个为-0.5。我本来以为，随着SCF的进行，最终的收敛结果将会得到类似分子团簇的自旋密度等值面图那样。然而最终的SCF结果是total magnetization为0.00 Bohr mag/cell，absolute magnetization则为0.19 Bohr mag/cell，从QE输出文件的Magnetic moment per site段落中看出只有钴上有磁矩，而A上的原子磁矩都为0。
学生有以下的问题，希望各位老师抽空指导：
（1）查看最终SCF结果的磁矩信息是否在我上述描述的那个输出文件的段落呢？还是说我还要进一步做其他的分析或者计算才能得知有机物A上有没有自旋密度呢？
（2）如果上述所得的SCF结果不是能量最低的反铁磁电子态，那么像这种单电子离域分布的体系，应该如何设置初始磁矩呢？
|
|
发表于 Post on 2024-2-18 15:32:53
收藏 Add to favorites
楼主 Author