I am trying to determine what the lowest energy spin configuration is for a variety of crystal structures of Perovskite Oxides. My research supervisor is also doing some of the ground state calculations in VASP so that we can compare results. We have found, that when using LDA+U, abinit predicts an antiferromagnetic structure whereas VASP predicts a ferromagnetic structure. I think I have all of the settings the same, so we are surprised by the qualitative difference.
I have attached the input files I have used for a particular calculation in abinit and (if anyone has any experience which might help diagnose this) a sample VASP input file for the same structure. I am probably just making a stupid mistake somewhere, so any help is much appreciated.
Qualitative Difference from VASP for Perovskites
Moderator: bguster
Qualitative Difference from VASP for Perovskites
- Attachments
-
- fmpol.in
- Structure optimization
- (646 Bytes) Downloaded 228 times
-
- fmpol_band.in
- Bandstructure calculation
- (825 Bytes) Downloaded 246 times
-
- INCAR.in
- VASP INCAR
- (378 Bytes) Downloaded 273 times
Re: Qualitative Difference from VASP for Perovskites
Dear user,
Could you try LDAUTYPE=1 with VASP ? LDAUTYPE=2 is a simplified approach (it should correspond to use U=U-J and J=0 in ABINIT).
LDAUTYPE=1 should compare well to the ABINIT result.
Best regards
Could you try LDAUTYPE=1 with VASP ? LDAUTYPE=2 is a simplified approach (it should correspond to use U=U-J and J=0 in ABINIT).
LDAUTYPE=1 should compare well to the ABINIT result.
Best regards
Bernard Amadon
CEA
France
CEA
France
Re: Qualitative Difference from VASP for Perovskites
Dear jleehand,
On the top of Bernard's comment I would also add that it can depend on the pseudopotential file you are using (the most common are JTH and GBRV tables).
Another point is that you defined a smearing value in Abinit (tsmear=0.005Ha=0.13605eV) while not in Vasp such that in your Vasp calculation you have the default smearing which is not 0.13605 eV. This smearing could be dramatic on the magnetic energy differences in such metallic case.
Best wishes,
Eric
On the top of Bernard's comment I would also add that it can depend on the pseudopotential file you are using (the most common are JTH and GBRV tables).
Another point is that you defined a smearing value in Abinit (tsmear=0.005Ha=0.13605eV) while not in Vasp such that in your Vasp calculation you have the default smearing which is not 0.13605 eV. This smearing could be dramatic on the magnetic energy differences in such metallic case.
Best wishes,
Eric