We present a first-principles investigation of the structural, electronic, and magnetic properties of pyrolusite ($beta$-MnO$_2$) using conventional and extended Hubbard-corrected density-functional theory (DFT+$U$ and DFT+$U$+$V$). The onsite $U$ an
d intersite $V$ Hubbard parameters are computed using linear-response theory in the framework of density-functional perturbation theory. We show that while the inclusion of the onsite $U$ is crucial to describe the localized nature of the Mn($3d$) states, the intersite $V$ is key to capture accurately the strong hybridization between neighboring Mn($3d$) and O($2p$) states. In this framework, we stabilize the simplified collinear antiferromagnetic (AFM) ordering (suggested by the Goodenough-Kanamori rule) that is commonly used as an approximation to the experimentally-observed noncollinear screw-type spiral magnetic ordering. A detailed investigation of the ferromagnetic and of other three collinear AFM spin configurations is also presented. The findings from Hubbard-corrected DFT are discussed using two kinds of Hubbard manifolds - nonorthogonalized and orthogonalized atomic orbitals - showing that special attention must be given to the choice of the Hubbard projectors, with orthogonalized manifolds providing more accurate results than nonorthogonalized ones within DFT+$U$+$V$. This work paves the way for future studies of complex transition-metal compounds containing strongly localized electrons in the presence of pronounced covalent interactions.
Recent neutron-diffraction experiments in honeycomb CrI$_3$ 2D ferromagnets have evinced the existence of a gap at the Dirac point in their spin-wave spectra. The existence of this gap has been attributed to strong Dzyaloshinskii-Moriya or Kitaev (DM
/K) interactions and suggested to set the stage for topologically protected edge states to sustain non-dissipative spin transport. We perform state-of-the-art simulations of the spin-wave spectra in this system, based on time-dependent density-functional perturbation theory and fully accounting for spin-orbit couplings (SOC) from which DM/K interactions ultimately stem. Our results are in fair overall agreement with experiments, but fail to account for the observed gap. The magnitude of the DM/K interaction parameters, determined from constrained density-functional theory (CDFT), is not consistent with observations either. Lattice-dynamical calculations, performed within density-functional perturbation theory (DFpT), indicate that a substantial degeneracy and a strong coupling between vibrational and magnetic excitations exist in this system, thus providing an alternative path to understanding the observed gap. In order to pursue this path, we introduce an interacting magnon-phonon Hamiltonian featuring a linear coupling between lattice and spin fluctuations, enabled by the magnetic anisotropy induced by SOC. Upon determination of the relevant interaction constants by DFpT and CDFT, this model allows us to propose magnon-phonon interactions as the most likely mechanism to explain the observed gap.
In Peierls-distorted materials, photoexcitation leads to a strongly coupled transient response between structural and electronic degrees of freedom, always measured independently of each other. Here we use transient reflectivity in the extreme ultrav
iolet to quantify both responses in photoexcited bismuth in a single measurement. With the help of first-principles calculations based on density-functional theory (DFT) and time-dependent DFT, the real-space atomic motion and the temperature of both electrons and holes as a function of time are captured simultaneously, retrieving an anticorrelation between the $A_{1g}$ phonon dynamics and carrier temperature. The results reveal a coherent, bi-directional energy exchange between carriers and phonons, which is a dynamical counterpart of the static Peierls-Jones distortion, providing first-time validation of previous theoretical predictions.
Electron-phonon ($e$-ph) interactions are pervasive in condensed matter, governing phenomena such as transport, superconductivity, charge-density waves, polarons and metal-insulator transitions. First-principles approaches enable accurate calculation
s of $e$-ph interactions in a wide range of solids. However, they remain an open challenge in correlated electron systems (CES), where density functional theory often fails to describe the ground state. Therefore reliable $e$-ph calculations remain out of reach for many transition metal oxides, high-temperature superconductors, Mott insulators, planetary materials and multiferroics. Here we show first-principles calculations of $e$-ph interactions in CES, using the framework of Hubbard-corrected density functional theory (DFT+$U$ ) and its linear response extension (DFPT+$U$), which can describe the electronic structure and lattice dynamics of many CES. We showcase the accuracy of this approach for a prototypical Mott system, CoO, carrying out a detailed investigation of its $e$-ph interactions and electron spectral functions. While standard DFPT gives unphysically divergent and short-ranged $e$-ph interactions, DFPT+$U$ is shown to remove the divergences and properly account for the long-range Frohlich interaction, allowing us to model polaron effects in a Mott insulator. Our work establishes a broadly applicable and affordable approach for quantitative studies of e-ph interactions in CES, a novel theoretical tool to interpret experiments in this broad class of materials.
The self-consistent evaluation of Hubbard parameters using linear-response theory is crucial for quantitatively predictive calculations based on Hubbard-corrected density-functional theory. Here, we extend a recently-introduced approach based on dens
ity-functional perturbation theory (DFPT) for the calculation of the on-site Hubbard $U$ to also compute the inter-site Hubbard $V$. DFPT allows to reduce significantly computational costs, improve numerical accuracy, and fully automate the calculation of the Hubbard parameters by recasting the linear response of a localized perturbation into an array of monochromatic perturbations that can be calculated in the primitive cell. In addition, here we generalize the entire formalism from norm-conserving to ultrasoft and projector-augmented wave formulations, and to metallic ground states. After benchmarking DFPT against the conventional real-space Hubbard linear response in a supercell, we demonstrate the effectiveness of the present extended Hubbard formulation in determining the equilibrium crystal structure of Li$_x$MnPO$_4$ (x=0,1) and the subtle energetics of Li intercalation.
We present a derivation of the exact expression for Pulay forces in density-functional theory calculations augmented with extended Hubbard functionals, and arising from the use of orthogonalized atomic orbitals as projectors for the Hubbard manifold.
The derivative of the inverse square root of the orbital overlap matrix is obtained as a closed-form solution of the associated Lyapunov (Sylvester) equation. The expression for the resulting contribution to the forces is presented in the framework of ultrasoft pseudopotentials and the projector-augmented-wave method, and using a plane wave basis set. We have benchmarked the present implementation with respect to finite differences of total energies for the case of NiO, finding excellent agreement. Owing to the accuracy of Hubbard-corrected density-functional theory calculations - provided the Hubbard parameters are computed for the manifold under consideration - the present work paves the way for systematic studies of solid-state and molecular transition-metal and rare-earth compounds.
We present a joint theoretical and experimental study of the oxygen $K$-edge spectra for LaFeO$_3$ and homovalent Ni-substituted LaFeO$_3$ (LaFe$_{0.75}$Ni$_{0.25}$O$_3$), using first-principles simulations based on density-functional theory with ext
ended Hubbard functionals and x-ray absorption near edge structure (XANES) measurements. Ground-state and excited-state XANES calculations employ Hubbard on-site $U$ and inter-site $V$ parameters determined from first principles and the Lanczos recursive method to obtain absorption cross sections, which allows for a reliable description of XANES spectra in transition-metal compounds in a very broad energy range, with an accuracy comparable to that of hybrid functionals but at a substantially lower cost. We show that standard gradient-corrected exchange-correlation functionals fail in capturing accurately the electronic properties of both materials. In particular, for LaFe$_{0.75}$Ni$_{0.25}$O$_3$ they do not reproduce its semiconducting behaviour and provide a poor description of the pre-edge features at the O $K$ edge. The inclusion of Hubbard interactions leads to a drastic improvement, accounting for the semiconducting ground state of LaFe$_{0.75}$Ni$_{0.25}$O$_3$ and for a good agreement between calculated and measured XANES spectra. We show that the partial substitution of Fe for Ni affects the conduction-band bottom by generating a strongly hybridized O($2p$)-Ni($3d$) minority-spin empty electronic state. The present work, based on a consistent correction of self-interaction errors, outlines the crucial role of extended Hubbard functionals to describe the electronic structure of complex transition-metal oxides such as LaFeO$_3$ and LaFe$_{0.75}$Ni$_{0.25}$O$_3$ and paves the way to future studies on similar systems.
Contradictory theoretical results for oxygen vacancies in SrTiO$_3$ (STO) were often related to the peculiar properties of STO, which is a $d^0$ transition metal oxide with mixed ionic-covalent bonding. Here, we apply, for the first time, density fun
ctional theory (DFT) within the extended Hubbard DFT+$U$+$V$ approach, including on-site as well as inter-site electronic interactions, to study oxygen-deficient STO with Hubbard $U$ and $V$ parameters computed self-consistently via density-functional perturbation theory. Our results demonstrate that the extended Hubbard functional is a promising approach to study defects in materials with electronic properties similar to STO. Indeed, DFT+$U$+$V$ provides a better description of stoichiometric STO compared to standard DFT or DFT+$U$, the band gap and crystal field splitting being in good agreement with experiments. In turn, also the description of the electronic properties of oxygen vacancies in STO is improved, with formation energies in excellent agreement with experiments as well as results obtained with the most frequently used hybrid functionals, however at a fraction of the computational cost. While our results do not fully resolve the contradictory findings reported in literature, our systematic approach leads to a deeper understanding of their origin, which stems from different cell sizes, STO phases, the exchange-correlation functional, and the treatment of structural relaxations and spin-polarization.
We present in full detail a newly developed formalism enabling density functional perturbation theory (DFPT) calculations from a DFT+$U$ ground state. The implementation includes ultrasoft pseudopotentials and is valid for both insulating and metalli
c systems. It aims at fully exploiting the versatility of DFPT combined with the low-cost DFT+$U$ functional. This allows to avoid computationally intensive frozen-phonon calculations when DFT+$U$ is used to eliminate the residual electronic self-interaction from approximate functionals and to capture the localization of valence electrons e.g. on $d$ or $f$ states. In this way, the effects of electronic localization (possibly due to correlations) are consistently taken into account in the calculation of specific phonon modes, Born effective charges, dielectric tensors and in quantities requiring well converged sums over many phonon frequencies, as phonon density of states and free energies. The new computational tool is applied to two representative systems, namely CoO, a prototypical transition metal monoxide and LiCoO$_2$, a material employed for the cathode of Li-ion batteries. The results show the effectiveness of our formalism to capture in a quantitatively reliable way the vibrational properties of systems with localized valence electrons.
We propose a self-consistent site-dependent Hubbard $U$ approach for DFT+$U$ calculations of defects in complex transition-metal oxides, using Hubbard parameters computed via linear-response theory. The formation of a defect locally perturbs the chem
ical environment of Hubbard sites in its vicinity, resulting in different Hubbard $U$ parameters for different sites. Using oxygen vacancies in SrMnO$_3$ as a model system, we investigate the dependence of $U$ on the chemical environment and study its influence on the structural, electronic, and magnetic properties of defective bulk and strained thin-film structures. Our results show that a self-consistent $U$ improves the description of stoichiometric bulk SrMnO$_3$ with respect to GGA or GGA+$U$ calculations using an empirical $U$. For defective systems, $U$ changes as a function of the distance of the Hubbard site from the defect, its oxidation state and the magnetic phase of the bulk structure. Taking into account this dependence, in turn, affects the computed defect formation energies and the predicted strain- and/or defect-induced magnetic phase transitions, especially when occupied localized states appear in the band gap of the material upon defect creation.
