I am studying Hematite (Fe2O3) Antiferromagnetic configuration using QE 7.0 (EDIT: compiled and tried 7.1) at the PBE+U level (U = 4.5).
I obtained an excellent agreement with experimental data on various properties, but there is one thing I can't wrap my head around.
The traces of the occupancy matrices after diagonalization for the minority spin on iron atoms are roughly close to 1 (0.88). I would expect them to be close to 0.
Iron atoms have a d5 configuration which I expect would lead to 5 electrons with either spin up or down. In my case, I get that the total number of electrons is 5 in the majority spin channel as expected, but roughly 1 in the minority channel. This is also the case in the example08 in the /PW/tutorials folder of QE for a FeO system.
I have been trying unsuccessfully to force the system into having the traces close to 0 for the minority channel.
- Applying starting magnetizations and charges on iron and oxygen atoms.
- Changing U values
- Reducing the Gaussian smearing to almost 0
- A mix of all the previous ones.
- Starting_ns_eigenvalues for the third eigenvalue in the minority channel initialised to 0
- Specified the Hubbard_occ to 5.00
I am now wondering if this is to be expected or maybe there is something I am missing about the theory or the implementation of DFT+U or maybe my understanding of the occupancy matrix is simply wrong.
EDIT: I am now trying the procedure following http://theossrv1.epfl.ch/Main/OxidationStates. I will try to investigate whether this works but also if there is a big difference in the results. I am wondering if, in the end, this occupancy problem is not a problem at all