# Full text of "Ferromagnetic insulating state in tensile-strained LaCoO$_3$ thin films"

## See other formats

(N oo c o o 00 en o (N X Ferromagnetic insulating state in tensile-strained LaCoOs thin films Han Hsu, 1 Peter Blaha, 2 and Renata M. Wentzcovitch 1 1 Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota, USA O - CS| ' 6 Institute of Materials Chemistry, Vienna University of Technology, & . 7 A-1060 Vienna, Getreidemarkt 9/165-TC, Austria (Dated: February 25, 2013) Abstract io With local density approximation + Hubbard U (LDA+C7) calculations, we show that the fer- ii romagnetic (FM) insulating state observed in tensile-strained LaCo03 epitaxial thin films is most 12 likely a mixture of low-spin (LS) and high-spin (HS) Co, namely, a HS/LS mixture state. Com- 13 pared with other FM states, including the intermediate-spin (IS) state (metallic within LDA+U), 14 which consists of IS Co only, and the insulating IS/LS mixture state, the HS/LS state is the most 15 favorable one. The FM order in HS/LS state is stabilized via the superexchange interactions be- i6 tween adjacent LS and HS Co. We also show that Co spin state can be identified by measuring the 17 electric field gradient (EFG) at Co nucleus via nuclear magnetic resonance (NMR) spectroscopy. is PACS numbers: 75.50.Cc, 75.70.Ak, 76.60.Gv 19 Perovskite-structure oxides have been proven a fertile area in condensed matter physics. 20 They exhibit amazing properties, including ferroelectricity, ferromagnetism, colossal magne- 21 toresistance (CMR), and multiferroics (simultaneous ferroelectricity and ferromagnetism), 22 as a consequence of their spin, lattice, charge, and orbital degree of freedom. Advances in 23 thin-film growth techniques have even brought more promising potentials for their future 24 application, as their properties can be engineered via epitaxial strains. A few examples 25 include strontium titanate (SrTiOs), ferroelectric in tensile-strained thin film while para- 26 electric in bulk, 1 lanthanum titanate (LaTiOs), conducting in compressive-strained thin 27 film while insulating in bulk, 2 and Europium titanate (EuTiOs), in which multiferroics in- 28 duced by tensile strains has been observed. 3 As to lanthanum cobaltite (LaCoOs), a diamag- 29 netic insulator in bulk at low temperatures (T < 35 K), a ferromagnetic (FM) insulating 30 state has been observed in tensile-strained thin films, e.g. LaCoOs grown on SrTiOs or 31 (LaA103)o.3(Sr 2 AlTa0 6 )o.7, at T < 85 K,-~— while the ferromagnetism induced by com- 32 pressive strains, e.g. LaCoOs grown on LaAlOs, is not conclusive .— — « Two questions arise 33 immediately: (1) Given that there are six 3d electrons in Co 3+ , which can thus have a to- 34 tal electron spin S = 0, 1, or 2, refereed to as low-spin (LS), intermediate-spin (IS), and 35 high-spin (HS) state, respectively, what is the spin state of Co in FM LaCoO"3 thin films, in 36 contrast to the LS Co in diamagnetic bulk? (2) What leads to the formation of FM order in 37 LaCoOs thin films? After all, FM insulators are rarely seen. So far, all first-principles cal- 38 culations have only found FM metallic LaCoOs thin films with all Co ions in IS state,— ^^ 39 a prediction clearly inconsistent with transport measurement s.- 40 While finite Co spin induced by tensile strains in LaCoOs thin films has just started 4i attracting attention, finite Co spin induced by thermal excitation in bulk LaCoOs has been 42 a highly controversial issue for decades.— 1 ^ With LS Co at T < 35 K, bulk LaCoOs becomes 43 a paramagnetic insulator with finite Co spin at about 90 K. Such a spin-state crossover in 44 the temperature range of 35-90 K was first suggested to be a LS-HS crossover— ~— but was 45 later suggested to be LS-IS based on a local density approximation + Hubbard U (LiDA+U) 46 calculation.— Since then, both scenarios have received supports from various experimental 47 and theoretical works, but a consensus is not yet achieved (see Ref. l23| for a brief review). 48 A study regarding LaCoOs thin films may also help understanding LaCoOs bulk from a 49 different perspective. In this paper, we investigate the Co spin state in tensile-strained thin 50 films and the formation of FM order via a series of LDA+U calculations. While LDA+U has 51 been frequently used to study cobaltites and Co spin state, the choice of Hubbard U can be 52 an issue. It has been shown that under the same lattice parameter, the Hubbard U affects the 53 total energy and the determination of ground state.— A well justified Hubbard U determined 54 by first principles would thus be necessary for finding out the the actual ground state. In 55 this paper, we compute the Hubbard U parameters of Co in all spin states self-consistently 56 with a linear response approach.— ~— This method has successfully found the ground state 57 of iron-bearing magnesium silicate (MgSiOs) perovskite at a wide range of pressure.— Both 58 the plane-wave pseudopotential (PWPP) method 2 - implemented in Quantum ESPRESSO 59 codes 28 and the augmented plane wave + local orbitals (APW+lo) method 29 implemented eo in WIEN2k codes 3 -^ are used. As shall be pointed out later, the orbital occupancies of Co in ei thin films are different from those in bulk, due to their different symmetries. We therefore 62 compute the electric field gradient (EFG) tensor at Co nucleus, V zz , with WIEN2k, to see 63 whether the Co spin state in thin films can be identified via EFG, as demonstrated in bulk. 23 64 The pseudocubic lattice parameter of bulk LaCo0 3 (R3c symmetry) is about 3.81 A 65 at T ~ 5 K.— »2S To model tensile-strained LaCoOs thin films via bulk calculations, we ee constrain the in-plane pseudocubic lattice parameters a pc and b pc of the hypothetical bulk 67 to 3.899 A (the lattice constant of cubic SrTiOs at low temperatures),- 3 set a = (3 = 7 = 6B 90°, and optimize the out-of-plane pseudocubic lattice parameter (c pc ). Due to the lack of 69 accurate information regarding C0O6 octahedral rotation in thin-film LaCoOs in the low- 70 temperature FM phase, we consider two extreme cases shown in Fig. [TJ (a) cube-on-cube, 71 namely, no CoOg octahedral rotation degree of freedom, and (b) full CoOg rotation degree 72 of freedom subject to the above-mentioned constraints. We also consider several magnetic 73 configurations shown in Fig. [2j (a) all Co ions in LS state, (b) all Co ions have the same 74 magnetic moment aligned in FM order, and (c) a mixture state with LS Co surrounded by 75 magnetic Co (and vice versa) aligned in FM order. The configuration shown in Fig. 2(c) is 76 a legitimate postulate, as the observed magnetization in LaCoOs thin films rarely exceeds 77 0.85 hb/Co.—~— For the configuration in Fig. [2(b), a convergent wave function for HS state 7B cannot be obtained; only IS state can be found. For the configuration in Fig. [2(c), both 79 HS/LS and IS/LS mixture state can be obtained. The self-consistent Hubbard U parameter so (U sc ) of LS and IS Co are 7.0 eV, while the HS/LS state has f/ s ( c HS) = 5.4 and u£ S) = 7.2 si eV.— The dependence of U sc on c pc is negligible. To demonstrate how the choice of Hubbard 82 U can affect the determination of ground state, we also present the result obtained using 83 a constant U = 7 eV for all Co in our PWPP calculations. In tensile-strained LaCo03 84 thin films, CoOg octahedra possess tetragonal symmetry, namely, longer Co-0 distance on 85 the xy plane, regardless of CoOg rotation, as will be discussed in the next paragraph. In 86 tetragonal symmetry (Ai/i), the spin-down electron of HS Co occupies d xy orbital, and the 87 spin-down electrons of IS Co occupy d xz and d yz orbitals, as shown in Fig. [2](d) . Such orbital as occupancies are very different from those in bulk LaCoOa, in which Co06 octahedra have 89 trigonal symmetry (D^d), with the [111] direction being the high-symmetry axis. In bulk 90 LaCoOa, the spin-down electron of HS Co occupies the 0^2-like orbital oriented along the 91 [111] direction, and the IS Co spin-down electrons occupy the doublet with 3-fold rotation 92 symmetry about the [111] direction.— 93 The optimized out-of-plane pseudocubic lattice parameter (c pc ) of each FM state and 94 associated relative energy (AE) and band gap (E gap ) are listed in Tables I-II. Regardless 95 of CoOq rotation, the HS/LS mixture state (with U sc ) is the most stable FM state given 96 by the PWPP method (Table l)M While the choice of U = 7.0 eV makes HS/LS state 97 less favorable in PWPP calculations, APW+lo calculations still find HS/LS the most stable 98 FM state (Table II). Both PWPP and APW+lo methods open an energy gap for HS/LS 99 state, consistent with transport measurements. 6 Also, the presence of HS Co is consistent ioo with recent x-ray magnetic circular dichroism (XMCD) and x-ray absorption spectroscopy ioi (XAS) spectra.- 1 ^ In contrast, the IS state is never the most favorable FM state, regardless 102 of the computation method and CoOg rotation. Its partially filled bands formed by partially 103 occupied d z 2 and d x 2_ y 2 orbitals in IS Co lead to a nonzero density of state at the Fermi 104 level. When the IS Co concentration is reduced to 50%, an energy gap is opened in APW+lo 105 calculations (Table II), while PWPP calculations still give a conducting IS/LS state (Table we I). Such a difference is likely to result from the way that the Hubbard U is applied. In 107 PWPP, the Hubbard U is applied to the projection of the total wave function onto Co 3d io8 orbitals; 24 in APW+lo, the Hubbard U is directly applied to the 3d orbitals within the 109 muffin-tin radius of Co (1.9 bohr). Insulating or not, the IS/LS state is highly unlikely; no its energy is even higher than that of IS state. One more thing worthy of mention is that in the Co-0 distance (rfco-o) an d Co-O-Co angle obtained in our calculation are different 112 from those estimated in Ref. |5|, where a constant dc -o = 1-93 A is assumed (regardless n3 of compressive or tensile strains), and an Co-O-Co angle of 176° is estimated for LaCoOs H4 grown on SrTi03. In our PWPP calculation with U sc and CoOg rotation, the HS/LS state us has dco(HS)-o — 2.015 and 1.872 A, g?co(ls)-o = 1.922 and 1.886 A, and Co-O-Co angles of n6 163.8° and 156.6°, on the xy plane and along the z axis, respectively. Even for the IS state H7 suggested by Ref. |5|, we have (Zco(iS)-O — 1-973 and 1.939 A, and Co-O-Co angles of 162.3° us and 154.9°, on the xy plane and along the z axis, respectively. n9 Other than the total energy, structural properties can be a useful criterion to determine 120 which FM state favors tensile strains (c pc /a pc < 1). Starting with the structures listed in 121 Table I with CoOg rotation, we perform full structural optimization (at constrained volume) 122 via variable cell-shape damped molecular dynamics. 3 - All lattice parameters, including a, 123 /3, and 7, are optimized, so the final structures only experience hydrostatic pressures. With 124 a, (3, and 7 slightly deviated from 90°, the IS state has c pc /a pc > 1, while all other states 125 remain c pc /a pc < 1 (but IS/LS state still has a c pc /a pc larger than that of HS/LS state), as 126 shown in Table III. The larger c pc /a pc ratio associated with IS Co is a direct consequence of 127 its occupied d xz and d yz orbitals (by spin-down electrons), which elongate the Co-0 distance i2B along the z-direction. In contrast, the fully optimized HS/LS state has c pc /a pc = 0.969, in 129 great agreement with c pc /a pc = 0.967 observed in experiments- 130 A significant part of cobalt-spin controversy arises from the difficulty in directly measuring 131 the total electron spin of Co. Such difficulty, also appearing in other spin systems, can be 132 resolved by comparing the calculated and measured EFGs.— -2SM. So far, insulating FM state 133 has been observed in LaCoOa thin films with a pc ranging from 3.84 to 3.90 A.-<£ In these 134 thin films, the magnetic Co concentration and Co-0 distance may be different, which can 135 lead to slightly different EFG for Co in the same spin state. To find out possible upper 136 and lower limits of HS and IS Co EFG, we compute them in two extreme cases: (1) thin 137 films with a pc = 3.899 A and 50% of magnetic Co, namely, the HS/LS and IS/LS states i3B listed in Tables I and II, and (2) single isolated HS or IS Co in an array of LS Co in a fully 139 relaxed structure with a pc ~ 3.81 A, where the orbital occupancies of isolated HS and IS Co 140 are maintained in tetragonal symmetry [Fig. |2{d)]. In these APW+lo calculations, IS Co 141 does not lead to a metallic state, in contrast to bulk LaCoOs {D^d symmetry).— Different 142 choices of Hubbard U have been adopted as well (5 eV, 7 eV, and U sc ). The results of all 143 these calculations show that the EFG mainly depends on the spin state: 14.7 < V zz /(10 21 144 V/m 2 ) < 19.9 and -14.6 < \4 IS) /(10 21 V/m 2 ) < -8.0. The quadrupole frequency, v Q = 145 3eQ\V xx \/2I (21 - 1) h, can thus be easily predicted, with Q = 0.42 x 10~ 28 m 2 and I = 7/2 we for 59 Co nucleus. Based on the range of V zz and V zz , we conclude that in insulating 147 LaCoOs thin films, a measured vq via nuclear magnetic resonance (NMR) spectroscopy us within ~ 8.2 ± 2.4 MHz can be a strong evidence for IS Co, and a measured vq within 149 ~ 12.6 ± 1.9 MHz should indicate HS Co. 150 Analysis of electronic structures can help developing a physical understanding for the FM 151 order in HS/LS state, whose projected density of states (PDOS) are shown in Fig. |3j The 152 case with U = 7 eV and no CoOg rotation is presented, as the main features in PDOS are 153 not sensitive to the choice of U and CoOg rotation. Extracted from Fig. [3]^a), both HS and 154 LS Co have nonzero magnetic moment: 2.97 and 0.56 p#, respectively. The FM order is 155 established via the superexchange interaction between HS and LS Co, as described by the 156 Goodenough-Kanamori rule,— >^~— which states that the superexchange interaction between 157 two cations (with or without a shared anion) is ferromagnetic if the electron transfer is from 158 a filled to a half-filled orbital or from a half-filled to an empty orbital. Indeed, for the HS/LS 159 state, electrons transfer from the filled d xz and d yz orbitals of LS Co to the half-filled d xz wo and d yz orbitals of HS Co via the oxygen in between, and also from the half-filled e g (d z 2 lei and d x 2_ y 2) orbitals of HS Co to the empty e g orbitals of LS Co, as depicted in the inset of 162 Fig. [3]^b). The PDOS shown in Fig. EJ^b) confirms this model: the finite spin- up e g electrons 163 localized at the LS Co site (transferred from the HS Co site) and the finite spin-down d xz 164 and d yz electrons localized at the HS Co site (transferred from the LS Co site). Such electron 165 transfers have been also described via a configuration fluctuation model,— ^ which further we suggests that the interchange of spin states (without net transfer of charge) led by electron 167 transfers stabilizes the FM order in this HS/LS states. lea The above-mentioned superexchange interaction can be visualized via electron spin den- 169 sity s(r) = p-|-(r) — p±(r), where Pt( r ) an d Pi( r ) are spin-up and spin-down electron density, 170 respectively. Figure Hl^a) shows s(r) corresponding to the configuration with all HS Co 171 magnetic moments aligned (same as the configuration in Fig. [3]). The nonzero magnetic 172 moments localized at LS Co sites (with e g character), aligned with the HS Co magnetic 173 moments, are consistent with the PDOS shown in Fig. [3j When the magnetic moment of 174 one HS Co in a 40-atom supercell is nipped [Fig. 11(b)], the alignment of magnetic moments 175 is altered, and so is the condition that allows configuration fluctuation [inset of Fig. 13(b)]. 176 The spin density at surrounding LS Co sites is thus significantly affected. One flipped HS 177 Co spin (in a 40-atom cell) increases the total energy by 195 meV/supercell. Flipping one i7B more HS Co spin, so the total magnetization per supercell becomes zero, further increases 179 the total energy by 78 meV/supercell. With CoOe rotation, the energy increases associated i8o with one and two flipped HS Co spins are 96 and 34 meV/supercell, respectively. These i8i results indicates that the magnetic moment of HS Co in the HS/LS state should align at 182 low temperatures. 183 While our calculations have shown that HS/LS state is the most favorable state among 184 the ferromagnetic states being considered, magnetic state in actual LaCo03 thin films 185 can be more complicated. The magnetization observed in experiments rarely exceeds 0.85 186 fiB /Co,— — smaller than that of the HS/LS mixture (2 //b/Co). Such a magnetization sug- 187 gests that the HS Co population should be smaller than 50%. In fact, XAS spectra combined 188 with atomic multiplet calculations have suggested that LaCoOs thin film on SrTiOs consists 189 of about 64% of LS Co and 36% of HS Co.— Given that the FM order is achieved via the 190 superexchange interaction within the HS-LS-HS Co configuration shown in Fig. E(c), one 191 can thus expect ferromagnetic HS/LS domain and nonmagnetic LS domain coexist in tensile 192 strained thin films, as observed using magnetic force microscopy (MFM). 7 Also, since HS/LS 193 state favors larger in-plane lattice parameter, thin films with larger in-plane lattice parame- 194 ters can be expected to have larger HS/LS domain, and thus larger magnetization, consistent 195 with the increase of magnetization with lattice parameter observed in experiment .— we In summary, we use LDA+U calculations to show that the ferromagnetic insulating state 197 in tensile-strained LaCoOa thin films is most likely a mixture of HS and LS Co. Among 198 all the ferromagnetic states studied in this paper (HS/LS, IS/LS, and IS), the insulating 199 HS/LS mixture state is the most favorable one, energetically and structurally. Its FM order 200 is established via the superexchange interaction between LS and HS Co. We also show that 201 cobalt spin states in LaCoOa thin films could be identified via NMR spectroscopy. 202 This work was primarily supported by the MRSEC Program of NSF grants DMR-02 12302 203 and DMR-0819885, and partially supported by EAR-081272 and EAR-1047629. P.B. was 204 supported by the Austrian Science Fund (SFB F41, " ViCoM"). Calculations were performed 205 at the Minnesota Supercomputing Institute (MSI). We thank C. Leighton for valuable dis- 206 cussions. 1 J. H. Haeni et al, Nature 430, 758 (2004). 2 F. J. Wong et al, Phys. Rev. B 81, 161101(R) (2010). 3 J. H. Lee et al, Nature 466, 954 (2010). 4 D. Fuchs et al, Phys. Rev. B 75, 144402 (2007). 5 D. Fuchs et al, Phys. Rev. B 77, 014434 (2008). 6 J. W. Freeland et al, Appl. Phys. Lett. 93, 212501 (2008). 7 S. Park et al, Appl. Phys. Lett. 95, 072508 (2009). 8 A. Herklotz et al, Phys. Rev. B 79, 092409 (2009). 9 V. V. Mehta et al, J. Appl. Phys. 105, 07E503 (2009). 10 M. Merz et al, Phys. Rev. B 82, 174416 (2010). 11 A. Posadas et al, Appl. Phys. Lett. 98, 053104 (2011). 12 C. Pinta et al, Phys. Rev. B 78, 174402 (2008). 13 D. Fuchs et al, Phys. Rev. B 79, 024424 (2009). 14 V. Mehta and Y. Suzuki, J. Appl. Phys. 109, 07D717 (2011). 15 K. Gupta and P. Mahadevan, Phys. Rev. B 79, 020406(R) (2009). 16 J. M. Rondinelli and N. A. Spaldin, Phys. Rev. B 79, 054409 (2009). 17 J. B. Goodenough, Localized to Itinerant Electronic Transition in Perovskite Oxides (Springer, 2001). 18 C. N. R. Rao et al, Top. Curr. Chem. 234, 1 (2004). 19 P. M. Raccah and J. B. Goodenough, Phys. Rev. 155, 932 (1967). 20 M. A. Senaris-Rodriguez and J. B. Goodenough, J. Solid State Chem. 116, 224 (1995). 21 M. A. Senaris-Rodriguez and J. B. Goodenough, J. Solid State Chem. 118, 323 (1995). 22 M. A. Korotin et al, Phys. Rev. B 54, 5309 (1996). 23 Han Hsu et al, Phys. Rev. B 82, 100406(R) (2010). 24 M. Cococcioni and S. de Gironcoli, Phys. Rev. B 71, 035105 (2005). 25 V. L. Campo Jr and M. Cococcioni, J. Phys.: Condens. Matter 22, 055602 (2010). 26 Han Hsu et al, Phys. Rev. Lett. 106, 118501 (2011). 27 Pseudopotentials used in this paper have been reported in H. Hsu et al., Phys. Rev. B 79, 125124 (2009). 28 P. Giannozzi et al, J. Phys.: Condens. Matter 21, 395502 (2009). 29 G. Madsen et al, Phys. Rev. B 64, 195134 (2001). 30 P. Blaha et al, WIEN2k, An Augmented Plane Wave Plus Local Orbitals Program for Calcu- lating Crystal Properties, edited by K. Schwarz, Techn. Universitat Wien, Vienna (2001). 31 P. G. Radaelli and S.-W. Cheong, Phys. Rev. B 66, 094408 (2002). 32 Whi Ref. iDA+C/ gives a pseudocubic lattice parameter of about 3.78 A for bulk LaCo03 (see 271 ). inclusion of zero point motion energy can improve this 1% underestimate, and better agreement with experiments can be achieved, as reviewed by R. M. Wentzcovitch et al, Rev. Min. Geochem. 71, 59 (2010). 33 F. W. Lytle, J. Appl. Phys. 35, 2212 (1964). 34 In this paper, f7 sc is extracted from a series of LDA+U ground states associ- ated with different trial U, as detatiled in the supplemental material of Ref. [26 24B (http://link.aps.org/supplemental/10.1103/PhysRevLett.106.118501 ). In our earlier work on LS 249 LaCo03 (Ref. 271 ). the Hubbard U for LS Co (~ 8.3 eV) was extracted from the LDA ground 24 250 state, as detailed in Ref. 251 35 In our calculations with constrained in-plane lattice parameters (a pc = b pc = 3.899 A), a HS 252 state with G-type antiferromagnetic (AFM) order can be obtained, and its energy is lower than 253 that of the HS/LS state. Given the neglegct of finite thickness, the possible uncertainty of energy 254 given by LDA+U method, and the fact that AFM thin film is not observed in experiments, we 255 limit our discussions on the available FM states (IS, IS/LS, and HS/LS). 256 36 R. M. Wentzcovitch et al, Phys. Rev. Lett. 70, 3947 (1993). 257 37 Han Hsu et al, Earth Planet. Sci. Lett. 294, 19 (2010). 25B 38 J. B. Goodenough, Phys. Rev. 100, 564 (1955). 259 39 J. B. Goodenough, J. Phys. Chem. Solids 6, 287 (1958). 260 40 J. Kanamori, J. Phys. Chem. Solids 10, 87 (1959). TABLE I. Optimized out-of-plane pseudocubic lattice parameter (c pc ) and associated relative en- ergy (AE) and energy gap (E gSu p) of each FM state in tensile-strained LaCo03 thin film (PWPP method). No CoOg rotation Cpc (A) AE (eV/f.u.) E ga ,p (eV) LS 3.865 0.35 0.54 IS 3.785 0.20 metal IS/LS 3.720 0.35 metal HS/LS (U sc ) 3.685 0.00 0.92 Full CoOg rotation Cpc (A) AE (eV/f.u.) £ga P (eV) 3.660 0.32 1.24 3.785 0.19 metal 3.680 0.29 metal 3.680 0.00 0.90 HS/LS (U = 7 eV) 3.700 0.29 1.12 3.695 0.29 0.90 TABLE II. Optimized c pc and associated AE and E gap of each FM state (APW+lo method, with CoOg rotation). Cpc (A) AE ' (eV/f.u.) £gap (eV) LS 3.660 0.37 1.72 IS 3.741 0.18 metal IS/LS 3.672 0.29 0.59 HS/LS (U = 7 eV) 3.686 1.52 TABLE III. Fully optimized pseudocubic lattice parameters of each FM state at the volume as in Table I (with CoOg rotation). %>C7 "pC (A) Cpc (A) Cpc/Opc IS 3.848 3.893 1.012 IS/LS 3.847 3.778 0.982 HS/LS (Use ) 3.863 3.745 0.969 HS/LS (U = 7 eV) 3.865 3.757 0.972 10 FIG. 1. (color online). Possible atomic structures of LaCoO"3 thin film (La is not shown) subject to constrained in-plane lattice parameters, (a) Cube on cube, no C0O6 octahedral rotation; (b) full octahedral rotation degree of freedom. FIG. 2. (color online), (a)-(c) Possible magnetic configurations in LaCo03 epitaxial thin film (La is not shown). The arrows denote for nonzero magnetic moments, either IS or HS. (a) LS state; (b) HS or IS state in FM order; (c) HS/LS or IS/LS mixture state in FM order, (d) The 3d orbitals occupied by the spin-down electrons of HS and IS Co. 11 Energy (eV) FIG. 3. (color online). Projected density of states of ferromagnetic HS/LS state (no CoOg rotation, £7 = 7 eV) onto (a) each atomic site, and (b) some of the Co 3d orbitals. The inset in (b) shows the electron transfer between HS and LS Co. FIG. 4. (color online). Spin density, s(r), of HS/LS state (no Co06 rotation, U = 7 eV) with (a) all HS Co magnetic moments aligned, and (b) one HS Co magnetic moment flipped downward (indicated by arrow). The isosurface values are 0.02 (red) and —0.02 (blue). 12