Page 1 of 1

Qualitative Difference from VASP for Perovskites

Posted: Tue Jul 30, 2019 7:02 pm
by jleehand
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.

Re: Qualitative Difference from VASP for Perovskites

Posted: Wed Jul 31, 2019 11:13 am
by amadon
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

Re: Qualitative Difference from VASP for Perovskites

Posted: Wed Jul 31, 2019 5:45 pm
by ebousquet
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