Benchmarking density functional theory on the prediction of antiferromagnetic transition temperatures

Abstract

This study investigates the predictive capabilities of common DFT methods (GGA, GGA+$U$, and GGA+$U$+$V$) for determining the transition temperature of antiferromagnetic insulators. We utilize a dataset of 29 compounds and derive Heisenberg exchanges based on DFT total energies of different magnetic configurations. To obtain exchange parameters within a supercell, we have devised an innovative method that utilizes null space analysis to identify and address the limitations imposed by the supercell on these exchange parameters. With obtained exchanges, we construct Heisenberg Hamiltonian to compute Transition temperatures using classical Monte Carlo simulations. To refine the calculations, we apply linear response theory to compute on-site ($U$) and intersite ($V$) corrections through a self-consistent process. Our findings reveal that GGA significantly overestimates the transition temperature (by ~113%), while GGA+$U$ underestimates it (by ~53%). To improve GGA+$U$ results, we propose adjusting the DFT results with the $(S+1)/S$ coefficient to compensate for quantum effects in Monte Carlo simulation, resulting in a reduced error of 44%. Additionally, we discover a high Pearson correlation coefficient of approximately 0.92 between the transition temperatures calculated using the GGA+$U$ method and the experimentally determined transition temperatures. Furthermore, we explore the impact of geometry optimization on a subset of samples. Using consistent structures with GGA+U and GGA+U+V theories reduced the error.

References in corpus (8)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Accurate and efficient algorithm for Bader charge integration
• Continuum variational and diffusion quantum Monte Carlo calculations
• Self-consistent Hubbard parameters from density-functional perturbation theory in the ultrasoft and projector-augmented wave formulations
• Calculation of the exchange constants of the Heisenberg model in the plane-wave based methods using the Green's function approach
• First-principles calculation of the electronic and optical properties of Gd${}_2$FeCrO${}_6$ double perovskite: Effect of Hubbard U parameter
• Lattice and spin dynamics in a low-symmetry antiferromagnet NiWO${}_4$
• A deep investigation of NiO and MnO through the first principle calculations and Monte Carlo simulations

Cited by in corpus (4)

• Origin of $A$-type antiferromagnetism and chiral split magnons in altermagnetic $\alpha$-MnTe
• Evaluating SCAN and r${}^2$SCAN meta-GGA functionals for predicting transition temperatures in antiferromagnetic materials
• Optimizing Supercell Structures for Heisenberg Exchange Interaction Calculations
• 1D Transition Metal Oxide Chains as a Challenging Model for Ab Initio Calculations