Spin dependence of ferroelectric polarization in the double exchange model for manganites

The double exchange (DE) model is systematically applied for studying the coupling between ferroelectric (FE) and magnetic orders in several prototypical types of multiferroic manganites. The model was constructed for the magnetically active Mn $3d$ bands in the basis of Wannier functions and include the effect of screened on-site Coulomb interactions. The essence of our approach for the FE polarization is to use the Berry phase theory, formulated in terms of occupied Wannier functions, and to evaluate the asymmetric spin-dependent change of these functions in the framework of the DE model. This enables us to quantify the effect of the magnetic symmetry breaking and derive several useful expressions for the electronic polarization ${\bf P}$, depending on the relative directions of spins. The proposed theory is applied to the solution of three major problems: (i) The magnetic-state dependence of ${\bf P}$ in hexagonal manganites; (ii) The microscopic relationship between canted ferromagnetism and ${\bf P}$ in monoclinic BiMnO$_3$; (iii) The origin of FE activity in orthorhombic manganites. We show that for an arbitrary noncollinear magnetic structure, propagating along the orthorhombic $\boldsymbol{b}$ axis and antiferromagnetically coupled $\boldsymbol{c}$, ${\bf P}$ can be obtained by scaling the one of the E-phase with the prefactor depending only on the relative directions of spins and being the measure of the spin inhomogeneity. This picture works equally well for the twofold (HoMnO$_3$) and fourfold (TbMnO$_3$) periodic manganites. The basic difference is that the twofold periodic magnetic structure is strongly inhomogeneous, that leads to large ${\bf P}$. On the contrary, the fourfold periodic magnetic structure can be viewed as a moderately distorted homogeneous spin spiral, which corresponds to weaker ${\bf P}$.Comment: 32 pages, 7 figure

Double-exchange theory of ferroelectric polarization in orthorhombic manganites with twofold periodic magnetic texture

We argue that many aspects of improper ferroelectric activity in orthorhombic manganites can be rationalized by considering the limit of infinite intra-atomic splitting between the majority- and minority-spin states (or the double exchange limit), which reduces the problem to the analysis of a spinless double exchange (DE) Hamiltonian. We apply this strategy to the low-energy model, derived from the first-principles calculations, and combine it with the Berry-phase theory of electric polarization. We start with the simplest two-orbital model, describing the behavior of the eg bands, and apply it to the E-type antiferromagnetic (AFM) phase, which in the DE limit effectively breaks up into one-dimensional zigzag chains. We derive an analytical expression for the electronic polarization (Pel) and explain how it depends on the orbital ordering and the energy splitting Delta between eg states. Then, we evaluate parameters of this model, starting from a more general five-orbital model for all Mn 3d bands and constructing a new downfolded model for the eg bands. From the analysis of these parameters, we conclude that the behavior of Pel in realistic manganites corresponds to the limit of large Delta. We further utilize this property in order to derive an analytical expression for Pel in a general two-fold periodic magnetic texture, based on the five-orbital model and the perturbation-theory expansion for the Wannier functions in the first order of 1/Delta. This expression explains the functional dependence of Pel on the relative directions of spins. Furthermore, it suggests that Pel is related to the asymmetry of the transfer integrals, which should simultaneously have symmetric and antisymmetric components. Finally, we explain how the polarization can be switched between orthorhombic directions a and c by inverting the zigzag AFM texture in every second ab plane.Comment: 41 page, 10 figure

Lattice Distortion and Magnetism of 3d-t2gt_{2g} Perovskite Oxides

Several puzzling aspects of interplay of the experimental lattice distortion and the the magnetic properties of four narrow $t_{2g}$-band perovskite oxides (YTiO$_3$, LaTiO$_3$, YVO$_3$, and LaVO$_3$) are clarified using results of first-principles electronic structure calculations. First, we derive parameters of the effective Hubbard-type Hamiltonian for the isolated $t_{2g}$ bands using newly developed downfolding method for the kinetic-energy part and a hybrid approach, based on the combination of the random-phase approximation and the constraint local-density approximation, for the screened Coulomb interaction part. Then, we solve the obtained Hamiltonian using a number of techniques, including the mean-field Hartree-Fock (HF) approximation, the second-order perturbation theory for the correlation energy, and a variational superexchange theory. Even though the crystal-field splitting is not particularly large to quench the orbital degrees of freedom, the crystal distortion imposes a severe constraint on the form of the possible orbital states, which favor the formation of the experimentally observed magnetic structures in YTiO$_3$, YVO_, and LaVO$_3$ even at the HF level. Beyond the HF approximation, the correlations effects systematically improve the agreement with the experimental data. Using the same type of approximations we could not reproduce the correct magnetic ground state of LaTiO$_3$. However, we expect that the situation may change by systematically improving the level of approximations for dealing with the correlation effects.Comment: 30 pages, 17 figures, 8 tables, high-quality figures are available via e-mai

Microscopic analysis of the magnetic form factor in low-dimensional cuprates

We analyze the magnetic form factor of Cu$^{2+}$ in low-dimensional quantum magnets by taking the metal-ligand hybridization into account explicitly. In this analysis we use the form of magnetic Wannier orbitals, derived from the first-principles calculations, and identify the contributions of different atomic sites. Having performed local density approximation calculations for cuprates with different types of ligand atoms, we discuss the influence of the on-site Coulomb correlations on the structure of the magnetic orbital. The typical composition of Wannier functions for copper oxides, chlorides and bromides is defined and related to features of the magnetic form factor. We propose easy-to-use approximations of the partial orbital contributions to the magnetic form factor in order to give a microscopic explanation for the results obtained in previous first-principles studies.Comment: 5 pages, 4 figure

Validity and limitations of the superexchange model for the magnetic properties of Sr2IrO4 and Ba2IrO4 mediated by the strong spin-orbit coupling

Layered perovskites Sr2IrO4 and Ba2IrO4 are regarded as the key materials for understanding the properties of magnetic relativistic insulators, mediated by the strong spin-orbit (SO) coupling. One of the most fundamental issues is to which extent these properties can be described by the superexchange (SE) model, formulated in the limit of the large Coulomb repulsion. In the present work we address this issue by deriving the relevant models and extracting parameters of these models from the first-principles calculations. First, we construct the effective Hubbard-type model for the t2g bands, by recasting the problem in the language of Wannier orbitals. Then, we map the obtained electron model onto the pseudospin model by applying the theory of SE interactions. We discuss the microscopic origin of anisotropic SE interactions, inherent to the compass Heisenberg model, and the appearance of the antisymmetric Dzyaloshinskii-Moriya term, associated with the additional rotation of the IrO6 octahedra in Sr2IrO4. In order to evaluate the Neel temperature (TN), we employ the non-linear sigma model. While for Sr2IrO4 our value of TN agrees with the experimental one, for Ba2IrO4 it is overestimated by a factor two. We argue that this discrepancy is related to limitations of the SE model: while for more localized t2g states in Sr2IrO4 it works reasonably well, the higher-order terms, beyond the SE model, play a more important role in the more "itinerant" Ba2IrO4, giving rise to the new type of isotropic and anisotropic exchange interactions. This conclusion is supported by unrestricted Hartree-Fock calculations for the same electron model, where in the case of Ba2IrO4, already on the mean-field level, we were able to reproduce the experimentally observed magnetic ground state, while for Sr2IrO4 the main results are essentially the same as in the SE model.Comment: 37 pages, 9 figure

Hybridization and spin-orbit coupling effects in quasi-one-dimensional spin-1/2 magnet Ba3Cu3Sc4O12

We study electronic and magnetic properties of the quasi-one-dimensional spin-1/2 magnet Ba3Cu3Sc4O12 with a distinct orthogonal connectivity of CuO4 plaquettes. An effective low-energy model taking into account spin-orbit coupling was constructed by means of first-principles calculations. On this basis a complete microscopic magnetic model of Ba3Cu3Sc4O12, including symmetric and antisymmetric anisotropic exchange interactions, is derived. The anisotropic exchanges are obtained from a distinct first-principles numerical scheme combining, on one hand, the local density approximation taking into account spin-orbit coupling, and, on the other hand, projection procedure along with the microscopic theory by Toru Moriya. The resulting tensors of the symmetric anisotropy favor collinear magnetic order along the structural chains with the leading ferromagnetic coupling J1 = -9.88 meV. The interchain interactions J8 = 0.21 meV and J5 = 0.093 meV are antiferromagnetic. Quantum Monte Carlo simulations demonstrated that the proposed model reproduces the experimental Neel temperature, magnetization and magnetic susceptibility data. The modeling of neutron diffraction data reveals an important role of the covalent Cu-O bonding in Ba3Cu3Sc4O12.Comment: 11 pages, 12 figure

Construction of Wannier functions from localized atomic-like orbitals

The problem of construction of the Wannier functions (WFs) in a restricted Hilbert space of eigenstates of the one-electron Hamiltonian $\hat{H}$ (forming the so-called low-energy part of the spectrum) can be formulated in several different ways. One possibility is to use the projector-operator techniques, which pick up a set of trial atomic orbitals and project them onto the given Hilbert space. Another possibility is to employ the downfolding method, which eliminates the high-energy part of the spectrum and incorporates all related to it properties into the energy-dependence of an effective Hamiltonian. We show that by modifying the high-energy part of the spectrum of the original Hamiltonian $\hat{H}$, which is rather irrelevant to the construction of WFs in the low-energy part of the spectrum, these two methods can be formulated in an absolutely exact and identical form, so that the main difference between them is reduced to the choice of the trial orbitals. Concerning the latter part of the problem, we argue that an optimal choice for trial orbitals can be based on the maximization of the site-diagonal part of the density matrix. The main idea is illustrated for a simple toy model, consisting of only two bands, as well as for a more realistic example of $t_{2g}$ bands in V$_2$O$_3$. An analogy with the search of the ground state of a many-electron system is also discussed.Comment: 13 pages, 6 figure

Realization of anisotropic compass model on the diamond lattice of Cu2+^{2+} in CuAl2_2O4_4

Spin-orbit (SO) Mott insulators are regarded as a new paradigm of magnetic materials, whose properties are largely influenced by SO coupling and featured by highly anisotropic bond-dependent exchange interactions between the spin-orbital entangled Kramers doublets, as typically manifested in $5d$ iridates. Here, we propose that a very similar situation can be realized in cuprates when the Cu$^{2+}$ ions reside in a tetrahedral environment, like in spinel compounds. Using first-principles electronic structure calculations, we construct a realistic model for the diamond lattice of the Cu$^{2+}$ ions in CuAl$_2$O$_4$ and show that the magnetic properties of this compound are largely controlled by anisotropic compass-type exchange interactions that dramatically modify the magnetic ground state by lifting the spiral spin-liquid degeneracy and stabilizing a commensurate single-$\boldsymbol{q}$ spiral