arxiv: v1 [cond-mat.mtrl-sci] 29 Sep 2010

Transcription

1 Spin-spiral inhomogeneity as origin of ferroelectricity in orthorhombic manganites I. V. Solovyev National Institute for Materials Science, -- Sengen, Tsukuba, Ibaraki , Japan (Dated: December, 08) arxiv: v [cond-mat.mtrl-sci] 9 Sep 00 We argue that homogeneous distribution of spins in state of orthorhombic manganites is deformed by relativistic spin-orbit interaction. The thus induced inhomogeneity gives rise to ferroelectric response of purely electronic origin in even-periodic magnetic structures. The mechanism is generic and explains appearance of ferroelectricity in twofold periodic structure of HoMnO as well as fourfold periodic structure of TbMnO. Neverless, odd-periodic magnetic structures preserve inversion symmetry and thus are not ferroelectric. Our analysis is based on low-energy model derived from first-principles electronic structure calculations and Berry-phase formalism for electronic polarization. The possibility of switching electric polarization P by means of magnetic field and spin magnetization by means of electric field has attracted enormous attention to multiferroic materials because of ir potential applicability in new electronic devises as well as fundamental interest to problem of coupling between magnetic and electronic degrees of freedom []. The orthorhombic rare-earth manganites RMnO ( space group Pbnm) are regarded as key materials for understanding details of such coupling. After discovering switching phenomena in TbMnO [], an unprecedented number of investigations has been carried out in order to clarify magnetic, structural, and ferroelectric properties of various RMnO compounds (e.g., Refs. [ ]). All of m are improper ferroelectrics, where inversion symmetry is broken due to some complex magnetic ordering []. Since Mn-atoms in Pbnm structure occupy inversion centers, latter can be destroyed only if magnetic unit is larger than crystallographic one. There are two types of compounds with different periodicity, which are typically regarded as representatives of two main ories of ferroelectricity in RMnO. The first one is (nearly) fourfold periodic TbMnO,whereferroelectricactivityisbelievedtobe due to spin-orbit interaction (SOI) related spiral spin alignment [6]. Anor one is twofold periodic HoMnO, where P is typically ascribed to (independent on SOI) magnetostriction effect in noncentrosymmetric E-type antiferromagnetic (AFM) structure [7]. Neverless, many questions remain. () There is no unique ory of multiferroicity in RMnO. Do HoMnO and TbMnO really behave as completely different systems, where ferroelectricity is caused by different microscopic mechanisms? Is it possible to unify se two cases? () When temperature decreases, many manganites, including TbMnO, exhibit a lock-in transition into a fourfold periodic commensurate phase []. What is origin of se commensurability? Is it different from twofold periodicity in HoMnO? () The experimental magnetic structure, derived for TbMnO and some related compounds by assuming model, rises many questions: even at low temperature, in order to fit experimental data, one had to rely on elliptical deformation of spin spiral [4]. However, in elliptical distribution, magnetic moments at certain Mn-sites are substantially reduced (up to about µ B [4]), which clearly contradicts to Hund s rule physics. Thus, model is probably incomplete. If so, what is true magnetic ground state of RMnO and how is it related to ferroelectric activity of se systems? In this work we will rationalize some of se questions. We will argue that ground state of RMnO is not spin spiral. Instead, we will introduce concept of inhomogeneous, where inhomogeneity is driven by relativistic SOI and is actually responsible for ferroelectric activity in RMnO. We will argue that this concept is applicable to all even-periodic systems, including TbMnO and HoMnO. The magnetic properties of manganites can be linked to electronic structure of Mnd-bands located near Fermi level. Thus, se bands, after transformation to real space, can serve as Wannierbasis for an effective low-energy model. For RMnO this basis includes three t g and two e g orbitals per spin for each of four Mn-sites in unit cell. The model Hamiltonian is taken in Hubbard form, where all parameters, such as crystal field, SOI, transfer integrals, and effective Coulomb interactions are calculated rigorously by starting from local-density approximation (LDA) and using experimental crystal structure [8]. The rare-earth 4f states are treated as core, which does not contribute to low-energy properties. The details of computational procedure can be found in review article[9]. The results of such calculationsfor wholeseriesofrmno compoundswithout SOI have been reported in Ref. [0]. Here, we only emphasize two points, which are important for magnetic inversion symmetry breaking: () Jahn-Teller distortion(jtd) gives rise to large(. ev) crystalfield splitting between e g levels, which manifests itself in orbital ordering (Fig. ); () The on-site Coulomb repulsion U is not particularly strong (. ev) due to veryefficient screeningby Opband[0], which isim-

2 J a J J J J y -r J b x -r E(q) (mev) 0 0 (,0) q x FIG.. (Color online) Left panel: schematic view on orbital ordering and magnetic interactions in ab-plane of RMnO. Typical values of magnetic interactions for TbMnO are J =.7, J a = 0., J b =., and J =. mev [0, ]. Rigt panel: dependence of total energy (per one Mn atom) on uniform vector q = (q x,q y,) obtained in HF calculations for TbMnO without SOI. portant for right balance between nearest-neighbor (NN) and some longer-range (LR) magnetic interactions, whose form is controlled by JTD [0, ]. Besides JTD, NN interactions in ab-plane depend on bucklingof Mn-O-Mnbonds. In least distorted LaMnO, J is ferromagnetic (FM) [0, ]. However, as hybridization between t g and e g states increases in more distorted compounds, J can become AFM [0]. Since on-site Coulomb repulsion is not particularly strong, besides conventional superexchange, re are or interactions, which appear in higher orders of /U, connect more remote sites, and compete with J. Among m, third-neighbor interaction J, operating via unoccupied e g states of intermediate Mn-sites, is strongest one. The second-neighbor interactions J a and Jb are weaker. Neverless, what important is anisotropy J b Ja, which predetermines direction of propagation (b) of and AFM E-state [0]. The NN interactions between adjacent ab-planes are strongly AFM (about 8 mev). The numerical values of J s for TbMnO are listed in Fig.. Similar behavior was found for HoMnO [0]. Without SOI, competition of isotropic interactions in ab-plane gives rise to incommensurate spinspiral state (Fig. ), which can be obtained by applying generalized Bloch orem [] and solving model in Hartree-Fock(HF) approximation[4]. Due to AFM interactions between planes, vector q, which specifies phase ϕ τ R = q (τ+r) + ατ of direction of spin e τ R = (cosϕτ R,sinϕτ R,0) at Mn-site τ of unit cell R, can be searched in form q=(q x,q y,), in units of reciprocal lattice translations. The total-energy minimum corresponds to homogeneous(α τ =0)propagationalongorthorhombicb-axis with q y =0.68 for TbMnO (and q y =0.7 for HoMnO ), which exceeds experimental values q y =0.8 and 0., reportedfor bc[]andab[]helix, respectively. This discrepancy will be resolved later by considering rel- (0,0) q y (0,) P a ( C/cm ) a b (deg) FIG.. (Color online) Nonrelativistic electronic polarization forgeneral periodicq y= structureoftbmno depending on angle φ between spin magnetic moments in two Mnsublattices. Relative directions of spins in ab plane are explained in inset. The numerical value of P in E- state (φ= 0 and 80 ) of TbMnO is. µc/cm (and.7 µc/cm in E-state of HoMnO ). ativistic SOI, but first we turn to analysis of electronic polarization, which is believed to be closely related to alignment [6]. For se purposes we compute P from HF eigenvectors C nk in Wannier-basis by using Berry-phase formalism on discrete grid of k-points []. Since Wannierbasis (and model itself) was constructed by starting from LDA bandstructure, which respects inversion symmetry, it does not contribute to P. The details can be found in Ref. [6]. First, we confirmed that without SOI, homogeneous state does not produce electronic polarization. A finite P can be indeed obtained by switching on relativistic SOI and performing one iteration by starting from self-consistent nonrelativistic HF potential. As a test example, let us consider q y = spiral in TbMnO, which was intensively discussed in literature [7]. Then, for bc and ab helix structure, we obtain P c µc/m and P a 470 µc/m, respectively. Thus, we confirm that our minimal model successfully reproduces results of first-principles LDA+U calculations for state [7]: () inequality P a P c, which holds forabandbchelixstructures; ()absolutevalues of P, which are not particularly large. Being encouraged by this good agreement, we turn to central part of our work, where we will show that () homogeneous state in RMnO is magnetically unstable as it tends to deform to some new inhomogeneous state under SOI; () this inhomogeneity gives rise to new electronic contribution to P. The simplest example, which illustrates how ferroelectricity can be induced by inhomogeneity without SOI is general q y = periodic structure (Fig. ). In this caseone can realize collinear(and in-

3 &$'(&)*+,-. & & &' & & % &'( ' & ( &' &() &() &()% &()' /0'& &)*+,-. &'& ()( ' L C/cm ) ( 0.4 L Energy gain (mev/mn) P a 0.0 Total Energy (mev/mn) 0.0 cterized by two $ $ $ $ E-state and & on angle %& %& lane: =0 and th opposite po& $ $ e homogeneous %&' ization. Thus, & nt to perturb direction of ar E-state this & t al. [7]. Then, & & polarization for nd on basis re is a large (Color online)spin Spin patterns patterns in of of..(color online) in ababplane plane s semiquantita-fig.fig. TbMnO obtained in HF calculations with SOI for = ons in Fig.. TbMnO obtained in HF calculations with SOI for qy = (a), (b), (c), (d), (e), and (f). The numbers (b), (c), 4 between (d), Mn-spins (e), and in6 different (f). The numbers this inhomoge-(a), indicate angles Mn-O-Mn angles between Mn-spins bonds along orthorhombic axis. in different Mn-O-Mn ity is of courseindicate ich leads to bonds along orthorhombic b axis. s and stabilizes ack to analmechanisms at homogeneous in sense that it is characterized by two o show that Mn-O-Mn bonds) E-state and ch gives rise tomagnetically different $ % $ % e naturally sta-homogeneous state, depending on angle φ hanism is verybetween two Mn-sublattices in ab-plane: φ=0 and riodic HoMnO 80 give rise to two E-phase domains with opposite pofig. 4. Left panel: electronic polarization in TbMnO de nce general-larization, corresponds. homogeneous pending while on φ=90 size of magnetic unitto cell Right panel: relativistic dependence of total energy for homogestructure with zeroonnet=/q polarization. Thus, ies correspondneous to TbMnO without SOI cirin order obtain states finiteofp, it is sufficient to(filled perturb e homogeneous cles), stabilization energy ( ) caused by SOI-induced state energies direction of magnetic relaxation (open circles), andintotal for magher iterate homogeneous E-state. (squares). collinear E-state this netic superstructures shownfor in Fig. ically, this pro-inhomogeneous s of iterations,mechanism was proposed by Sergienko et al. [7]. Then, of total en-picozzi et al. considered behavior of polarization for or spins compromise between NN, LR, and and relativistic r TbMnO are more general noncollinear alignment on basis anisotropic and Dzyaloshinsky-Moriya interactions. For magnetic struc-of first-principles calculations argued that re is a large structure one can clearly distinguish different n of angles electronic contribution to P [8], which is semiquantita domains formed by NN interactions and singlends, which indiby our model calculations inall Fig.. ionreproduced anisotropy, respectively. Most importantly, oddpin patterns re-tively periodic restore inversion centers magnetic questionstructures is how to stabilize this inhomoges in system: Now (associated withstructure? central Mn-sites in Fig. )isand he uniform spinneous One possibility of thus course exclude any ferroelectric activity. ples spin mag- nonrelativistic exchange striction, which leads to Thus, we predictatomic oscillatory behavior ofanddependestore com-off-centrosymmetric displacements stabilizes ing on =/L (Fig. 4). The inhomogeneity tic moments lie AFM E-state [7, 8]. We will come back to analin even-periodic structures gives rise to electronic l canting in ysis of relative roles played by different mechanisms polarization, similar to E-state, which decreases with at mainly formed of paper, but canting first weofwant show that increase of. The spins to away from anisotropy [], end magnetic structure,reduces which gives rise collinear arrangement substantially values of to e due to pe-inhomogeneous ture, spins finite. For example, electronic polarization =.staferroelectric polarization, can be naturally bonds and min-bilized C/cm, obtained for SOI. This structure of TbMnO by relativistic mechanism is very Moreover, generic (Fig.and 4), takes is reduced by factor two in comparison with place in twofold periodic HoMnO nimizes en collinear E-state (Fig. ). Taking into account as well as fourfold periodic TbMnO. Since general structures, numerical values of (Fig. ), this reduction is readily Bloch orem is no longer valid spins minimizeized explained by angle-dependence of in in relativistic nonrelwe consider ).supercell geometries he directions ofcase,ativistic case (Fig. Thus, SOI only correspondforms. ing to qy =/L for which we start from homogeneous state, turn on SOI and furr iterate HF equations until self-consistency. Typically, this procedure requires several tens of thousands of iterations, which are accompanied by decrease of total energy. The results of such calculations for TbMnO are summarized in Fig.. The equilibrium magnetic struc- FIG. 4. Left panel: electronic polarization in TbMnO depending on size of magnetic unit cell L. Right panel: dependence of total energy on L=/qy for homogeneous states of TbMnO without SOI (filled circles), stabilization energy ( E) caused by SOI-induced magnetic relaxation (open circles), and total energies for magnetic superstructures shown in Fig. (squares). tures are characterized by an alternation of angles formed by spins in different Mn-O-Mn bonds, which indicates inhomogeneity. The obtained spin patterns reflect competition of many interactions in system: isotropic interactions tend to form uniform spin spiral, while relativistic SOI, which couples spin magnetic moments to lattice, tends to restore commensurability in system. The magnetic moments lie mainly in ab-plane (apart from small canting in c-direction) [9]. The qy = structure is mainly formed by NN interactions and magnetic anisotropy [], while LR interactions are not affective due to periodicity constraint. For qy = structure, spins are aligned parallel to longest Mn-O bonds and minimize single-ion anisotropy energy. Moreover, magnetic coupling along b-axis minimizes energy of LR interactions. For large-l structures, situation is even more complex: certain spins minimize single-ion anisotropy energy, while directions of or spins compromise between NN, LR, and relativistic anisotropic and Dzyaloshinsky-Moriya interactions. For qy = 6 structure one can clearly distinguish different domains formed by NN interactions and singleion anisotropy, respectively. Most importantly, all oddperiodic magnetic structures restore inversion centers (associated with central Mn-sites in Fig. ) and thus exclude any ferroelectric activity. Thus, we predict oscillatory behavior of P depending on qy =/L (Fig. 4). The inhomogeneity in even-periodic structures gives rise to electronic polarization, similar to E-state, which decreases with increase of L. The canting of spins away from collinear arrangement substantially reduces values of P. For example, electronic polarization P =. µc/cm, obtained for qy = structure of TbMnO (Fig. 4), is reduced by factor two in comparison with collinear E-state (Fig. ). Taking into account numerical values of φ (Fig. ), this reduction is readily explained by angle-dependence of P in nonrelativistic case (Fig. ). Thus, SOI only forms

4 4 inhomogeneous state, while microscopic mechanism yielding finite P is essentially nonrelativistic. The existing ories overestimate values of P (sometimes by an order of magnitude) [7, 8], which is one of unresolved problems for multiferroic manganies. In our calculations we obtain P a=.4 µc/cm for q y = structure of HoMnO (φ=0 ) and 0. µc/cm for q y = 4 structure of TbMnO shown in Fig. 4. Both values are about three times larger than experimental ones []. Neverless, our oretical approach also suggests that it is difficult to obtain a good quantitative agreement because P is sensitive to many factors: () precise value of φ, which itself depends on fragile balance of many magnetic interactions; () structural relaxation[7, 8]; () due to oscillatory behavior (Fig. 4), P can be reduced by possible deviations from even-periodic commensurate alignment. Finally, we comment on relative roles played by magnetic and structural relaxation effects in formation of inhomogeneous magnetic structures. The behavior of magnetic stabilization energies E, defined as energy difference between fully optimized magnetic structure with SOI and homogeneous state without SOI, is explained in Fig. 4. Since each increase of supercell provides additional degrees of freedom for magnetic relaxation, E decreases with increase of L. The absolute value E 6 mev/mn obtained for q y = structure of HoMnO is comparable with energy gain caused by noncentrosymmetric atomic displacements in E-state ( 8 mev/mn [8]). Thus, even for q y =, magnetic relaxation cannot be neglected and should be considered on an equal footing with structural relaxation. Moreover, transition to inhomogeneous state is driven by SOI, prior to structural relaxation, and in this sense re is no conceptual difference between HoMnO and TbMnO. Since E increases with L, relative role of magnetic relaxation is also expected to increase. The magnetic relaxation alone readily explains experimentally observed q y 4 periodicity in TbMnO : E of relativistic origin shifts total energy minimum of homogeneous state (q y =0.68, Fig. ) towards smaller q y. Thus, new oretical minimum corresponds to q y = 4 (Fig. 4). Similar behavior was found for HoMnO, which is expected to form similar fourfold periodic structure. The disagreement with experimentally observed twofold periodicity in HoMnO is probably caused by neglect of structural relaxation, which is relatively more important for q y =. Quantitative aspects of interplay between magnetic and structural relaxation effects should be addressed in future oretical studies. At present, structural relaxation cannot be easily implemented in present model analysis. In summary, we propose that homogeneous spinspiral state in orthorhombic manganites is deformed by relativistic SOI. The thus induced inhomogeneity is responsible for ferroelectric activity of RMnO. It would be interesting to check our finding experimentally. This work is partly supported by Grant-in-Aid for Scientific Research (C) No from Ministry of Education, Culture, Sport, Science and Technology of Japan and by Russian Federal Agency for Science and Innovations, grant No SOLOVYEV.Igor@nims.go.jp [] Y. Tokura, Science, 48 (006); S.-W. Cheong and M. Mostovoy, Nature Materials 6, (007); D. Khomskii, Physics, 0 (009). [] T. Kimura et al., Nature 46, (00). [] T. Arima at al., Phys. Rev. B 7, 000 (00). [4] M. Kenzelmann at al., Phys. Rev. Lett. 9, (00); T. Arima at al., ibid. 96, 0970 (006); Y. Yamasaki at al., ibid. 98, 4704 (007); Y. Yamasaki at al., ibid. 0, (008). [] J. Hemberger et al., Phys. Rev. B 7, 08 (007); S. Ishiwata et al., ibid. 8, 004 (00). [6] H. Katsura, N. Nagaosa, and A. V. Balatsky, Phys. Rev. Lett. 9, 070 (00); M. Mostovoy, ibid. 96, (006); I. A. Sergienko and E. Dagotto, Phys. Rev. B 7, (006). [7] I. A. Sergienko, C. Şen, and E. Dagotto, Phys. Rev. Lett. 97, 704 (006). [8] J. Blasco et al., Phys. Rev. B 6, 609 (000); A. Muñoz et al., Inorg. Chem. 40, 00 (00). [9] I. V. Solovyev, J. Phys.: Cond. Matter 0, 90(008). [0] I. Solovyev, J. Phys. Soc. Jpn. 78, 0470 (009). [] The interatomic interactions are defined as one-half of HF energy difference between AFM and FM configurations in each bond. [] I. Solovyev, N. Hamada, and K. Terakura, Phys. Rev. Lett. 76, 48 (996). [] L. M. Sandratskii, Adv. Phys. 47, 8 (998). [4] Once orbital degeneracy is lifted by JTD, HF ory provides a good approximation for groundstate properties. The effect of correlation interactions, which can be treated as a perturbation to HF ground state [9], is partially compensated by magnetic polarization of Op states: if former tend to stabilize AFM structures, latter favors FM alignment. Details can be found in Ref. [0]. [] D. Vanderbilt and R. D. King-Smith, Phys. Rev. B 48, 444 (99); R. Resta, J. Phys.: Cond. Matter, 0 (00). For practical purposes we use grid of k-points (and proportionally decrease number of divisions in supercell calculations). [6] I. V. Solovyev and Z. V. Pchelkina, Phys. Rev. B 8, 0944 (00). [7] A. Malashevich and D. Vanderbilt, Phys. Rev. Lett. 0, 070 (008); Phys. Rev. B 80, 4407 (009); H. J. Xiang et al., Phys. Rev. Lett. 0, 0709 (008). [8] S. Picozzi et al., Phys. Rev. Lett. 99, 70 (007). [9] When we start with bc-helix structure, magnetic moments eir flop into ab plane (and converge to same solution as for ab-helix) or form anor inhomogeneous structure with higher energy. The ab-

5 planar structure minimizes anisotropy energy in Mn-sublattice, while bc-planar structures are probably stabilized by magnetism of R-ions [], which is not considered in present work.