Abstract:The magnetic, noncollinear parametrization of Dudarev's DFT+$U$ method is generalized to fully-relativistic ultrasoft pseudopotentials. We present the definition of the DFT+$U$ total energy functional, and the calculation of forces and stresses in the case of orthogonalized atomic orbitals defining the localised Hubbard manifold, where additional contributions arising from the derivative of the inverse square root of the overlap matrix appear. We further extend the perturbative calculation of the Hubbard $U$ parameters within density-functional perturbation theory to the noncollinear relativistic case, by exploiting an existing and recently developed theoretical approach that takes advantage of the time-reversal operator to solve a second Sternheimer equation. We validate and apply the new scheme by studying the electronic structure and the thermodynamics of the binary compounds EuX (where X = O, S, Se, Te is a chalcogen atom), as representative simple crystals, and of the pyrochlore Cd$_2$Os$_2$O$_7$, representative of a more structurally complex oxide.

What problem does this paper attempt to address?

The problem this paper attempts to address is: how to accurately describe the electronic structure and thermodynamic properties of materials containing non-collinear magnetism and spin-orbit coupling (SOC) using first-principles methods. Specifically, the authors propose an extended DFT+U method that combines fully relativistic ultrasoft pseudopotentials (FR US-PP) to handle the effects of non-collinear magnetism and SOC on the electronic structure. ### Main Issues: 1. **Accuracy of Electronic Structure**: Traditional density functional theory (DFT) has limitations when dealing with materials containing transition metals or rare earth atoms, especially for strongly correlated electron systems such as 3d, 4d, and 5d electrons. These systems have significant self-interaction errors, leading to inaccurate predictions of the electronic structure. 2. **Non-collinear Magnetism and Spin-Orbit Coupling**: Many materials exhibit non-collinear magnetism and significant spin-orbit coupling effects, which further complicate the electronic structure calculations. Existing DFT+U methods are primarily designed for collinear magnetic systems and cannot be directly applied to non-collinear magnetic systems. 3. **Determination of Hubbard U Parameter**: The choice of the Hubbard U parameter significantly affects the results of the DFT+U method. Traditional DFT+U methods often rely on empirical or semi-empirical methods to determine the U parameter, lacking first-principles calculation methods. ### Solutions: 1. **Extended DFT+U Method**: The authors extend the DFT+U method to the framework of non-collinear magnetism and fully relativistic ultrasoft pseudopotentials by introducing new total energy functionals and methods for calculating forces and stresses, ensuring accuracy when dealing with non-collinear magnetism and SOC. 2. **First-Principles Calculation of Hubbard U Parameter**: Using density functional perturbation theory (DFPT), the authors propose a method to calculate the Hubbard U parameter under conditions of non-collinear magnetism and SOC. This method solves two independent Sternheimer equations, considering the breaking of time-reversal symmetry, thus avoiding the introduction of adjustable parameters. 3. **Application Verification**: The authors apply the new method to the study of EuX series compounds (EuO, EuS, EuSe, EuTe) and Cd2Os2O7, verifying the method's effectiveness. By calculating properties such as the electronic structure, lattice parameters, bulk modulus, and critical pressure of these materials, they demonstrate the method's performance in practical applications. ### Application Examples: - **EuX Series Compounds**: The electronic structure and thermodynamic properties of EuX compounds under different pressures were studied, including lattice parameters, bulk modulus, and critical pressure. The results show that the new method can accurately describe the insulating properties and magnetic behavior of these materials. - **Cd2Os2O7**: The electronic structure and energy of Cd2Os2O7 under different non-collinear magnetic arrangements were studied, analyzing the impact of magnetic modes on the crystal structure. Through these applications, the authors demonstrate the effectiveness and reliability of the new method in handling complex material systems.

Noncollinear DFT+$U$ and Hubbard parameters with fully-relativistic ultrasoft pseudopotentials

Orbital-resolved DFT+U for molecules and solids

Hubbard $U$ through polaronic defect states

Orbital-Resolved DFT +U for Molecules and Solids

Explicit demonstration of the equivalence between DFT+U and the Hartree-Fock limit of DFT+DMFT

High-throughput determination of Hubbard and Hund values for transition metal oxides via the linear response formalism

Communication: Comparing Ab Initio Methods of Obtaining Effective U Parameters for Closed-Shell Materials

Magnons from time-dependent density-functional perturbation theory and the noncollinear Hubbard formulation

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

DFT+U Type Strong Correlation Functional Derived from Multiconfigurational Wavefunction Theory

A Self-consistent DFT+DMFT scheme in the Projector Augmented Wave : Applications to Cerium, Ce2O3 and Pu2O3 with the Hubbard I solver and comparison to DFT+U

Hubbard-$U$ corrected Hamiltonians for non-self-consistent random-phase approximation total-energy calculations: A study of ZnS, TiO$_2$, and NiO

Computing total energies in complex materials using charge self-consistent DFT+DMFT

Minimum tracking linear response Hubbard and Hund corrected Density Functional Theory in CP2K

Self-consistent on-site and inter-site Hubbard parameters within DFT+U+V for UO$_2$ using density-functional perturbation theory

DFT+U+J with linear response parameters predicts non-magnetic oxide band gaps with hybrid-functional accuracy

Nonadiabatic Exchange-Correlation Potential for Strongly Correlated Materials in the Weak and Strong Interaction Limits

Fermi surface evolution of the 2D Hubbard model within a novel four-pole approximation

Parametrization of LSDA+$U$ for noncollinear magnetic configurations: Multipolar magnetism in UO$_2$

Dynamical Linear Response Of Tddft With Lda Plus U Functional: Strongly Hybridized Frenkel Excitons In Nio

Auxiliary Field Quantum Monte Carlo for Multiband Hubbard Models: Controlling the Sign and Phase Problems to Capture Hund's Physics