DFT+$U$ within the framework of linear combination of numerical atomic orbitals

Abstract

We present a formulation and implementation of the DFT+\textit{U} method within the framework of linear combination of numerical atomic orbitals (NAO). Our implementation not only enables single-point total energy and electronic-structure calculations but also provides access to atomic forces and stresses, hence allowing for full structure relaxations of periodic systems. Furthermore, our implementation allows one to deal with non-collinear spin texture, with the spin-orbit coupling (SOC) effect treated self-consistently. The key aspect behind our implementation is a suitable definition of the correlated subspace when multiple atomic orbitals with the same angular momentum are used, and this is addressed via the "Mulliken charge projector" constructed in terms of the first (most localized) atomic orbital within the $d/f$ angular momentum channel. The important Hubbard $U$ and Hund $J$ parameters can be estimated from a screened Coulomb potential of the Yukawa type, with the screening parameter either chosen semi-empirically or determined from the Thomas-Fermi screening model. Benchmark calculations are performed for four late transition metal monoxide bulk systems, i.e., MnO, FeO, CoO, and NiO, and for the 5$d$-electron compounds IrO${}_2$. For the former type of systems, we check the performance of our DFT+$U$ implementation for calculating band gaps, magnetic moments, electronic band structures, as well as forces and stresses; for the latter, the efficacy of our DFT+$U$+SOC implementation is assessed. Systematic comparisons with available experimental results, and especially with the results from other implementation schemes are carried out, which demonstrate the validity of our NAO-based DFT+$U$ formalism and implementation.

arXiv comment19 pages, 7 figures, 7 tables

References in corpus (24)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• A linear response approach to the calculation of the effective interaction parameters in the LDA+U method
• Optimized norm-conserving Vanderbilt pseudopotentials
• Electronic Structure Calculations with Dynamical Mean-Field Theory: A Spectral Density Functional Approach
• High-throughput electronic band structure calculations: challenges and tools
• Optimization Algorithm for the Generation of ONCV Pseudopotentials
• Hubbard-corrected DFT energy functionals: the LDA+U description of correlated systems
• Frequency-dependent local interactions and low-energy effective models from electronic structure calculations
• Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions
• Ab initio calculations of quasiparticle band structure in correlated systems: LDA++ approach
• Implementation of the Projector Augmented Wave LDA+U Method: Application to the Electronic Structure of NiO
• Implementation of the LDA+U method using the full potential linearized augmented plane wave basis
• Screened Coulomb interaction in the maximally localized Wannier basis
• Ab-initio procedure for effective models of correlated materials with entangled band structure
• Anisotropy and Magnetism in the LSDA+U Method
• Localized Orbital Scaling Correction for Systematic Elimination of Delocalization Error in Density Functional Approximations
• Systematically improvable optimized atomic basis sets for {\it ab inito} calculations
• Multipole decomposition of LDA+$U$ energy and its application to actinides compounds
• All-electron periodic $G_0{W}_0$ implementation with numerical atomic orbital basis functions: algorithm and benchmarks
• Effect of spin orbit coupling and Hubbard $U$ on the electronic structure of IrO${}_2$
• Calculation of Effective Coulomb Interaction for $P{r}^{3+}$, $U^{4+}$, and $UP{t}_3$
• Comparing particle-particle and particle-hole channels of random-phase approximation

Cited by in corpus (1)

• 1D Transition Metal Oxide Chains as a Challenging Model for Ab Initio Calculations