Interfacial-hybridization-modified Ir Ferromagnetism and Electronic Structure in LaMnO3/SrIrO3 Superlattices
Abstract
Artificially fabricated 3/5 superlattices (SLs) involve both strong electron correlation and spin-orbit coupling in one material by means of interfacial 3-5 coupling, whose mechanism remains mostly unexplored. In this work we investigated the mechanism of interfacial coupling in LaMnO3/SrIrO3 SLs by several spectroscopic approaches. Hard x-ray absorption, magnetic circular dichroism and photoemission spectra evidence the systematic change of the Ir ferromagnetism and the electronic structure with the change of the SL repetition period. First-principles calculations further reveal the mechanism of the SL-period dependence of the interfacial electronic structure and the local properties of the Ir moments, confirming that the formation of Ir-Mn molecular orbital is responsible for the interfacial coupling effects. The SL-period dependence of the ratio between spin and orbital components of the Ir magnetic moments can be attributed to the realignment of electron spin during the formation of the interfacial molecular orbital. Our results clarify the nature of interfacial coupling in this prototypical 3/5 SL system and the conclusion will shed light on the study of other strongly correlated and spin-orbit coupled oxide hetero-interfaces.
Introduction
Entanglement of charge, spin, lattice and orbital degrees of freedom in transition metal oxides (TMOs) has attracted a great amount of research attentions recently [1_Imada, 2_Dagotto, 3_Quintanilla, 4_Keimer, 5_Ramirez]. Strong electron correlation in TMOs is a necessity to support the existence of local magnetic moments and various magnetic/charge/orbital orderings. Meanwhile, novel topics related to spin-orbit coupling (SOC) have also become research hot spots in condensed matter physics [6_Soumyanarayanan, 7_Qi]. SOC can work in real space and give rise to a variety of non-collinear magnetic structures such as skyrmions [6_Soumyanarayanan, 8_Muhlbauer], or in reciprocal -space to produce topologically non-trivial band structures [6_Soumyanarayanan, 7_Qi, 9_Bansil]. Consequently, fabrication and investigation of systems involving both of these two interactions become of not only scientific but also technical interest [10_Pesin, 11_Kim]. Local magnetic moments and orderings within a strong SOC regime still remain mostly unexplored. TMOs with heavy transition metals, such as iridates, are ideal hosts for coexistence of both significant electron correlation and SOC. An unprecedented eff=1/2 Mott insulating state with canted antiferromagnetic (AFM) ordering was realized in layered perovskite iridates such as Sr2IrO4 [11_Kim, 12_Moon, 13_Kim] and Sr3Ir2O7 [12_Moon, 14_Fujiyama], where the collaboration of strong electron correlation and SOC plays a crucial role. Nevertheless, generally 5 TMOs are not capable to sustain magnetic orderings due to the large spatial extension of the 5 electrons [15_Cao]. While their well-investigated 3 TMO counterparts usually possess weak SOC, even though the strength of electron correlation is always sufficient to support magnetism.
Interfaces between dissimilar materials can provide an intriguing playground for manipulation of various physical properties [16_Hwang, 17_Zubko, 18_Matsuno, 19_Ohmoto, 20_Hellman]. Great improvement of thin-film fabrication techniques enables accurately controlled design of epitaxial TMO heterostructures and SLs with atomically abrupt interfaces. It appears to be a natural strategy that artificial 3/5 TMO heterostructures or SLs are promising candidates to involve both significant electron correlation and SOC simultaneously. Pioneering research about 3/5 SLs was triggered by investigation on SrIrO3/SrTiO3 (SIO/STO) perovskite SLs [21_Matsuno], as a comparison with the Ruddlesden-Popper series iridates Srn+1IrnO3n+1. Meanwhile, strong interfacial 3-5 coupling was reported in La1-xSrxMnO3/SrIrO3 (0x1, LSMO/SIO) SLs [22_Nichols, 23_Yi, 24_Yi, 25_Kim, 26_Huang]. Emergent Ir ferromagnetic (FM) moments can be induced by the interfacial coupling with Mn FM moments, and in turn the magnetic properties of Mn layers can be significantly modified as well. Perpendicular magnetic anisotropy and concomitant anomalous Hall effect were observed in SLs [22_Nichols], and modulation of magnetic anisotropy in LSMO layers was also studied [23_Yi, 24_Yi]. Recent reports claim that interfacial hybridization between Ir and Mn is responsible for the charge transfer in SLs [27_Okamoto] and spectroscopic properties of SLs [25_Kim].
Perovskite SLs are ideal systems for investigation of interfacial coupling mechanisms thanks to their high interface quality. The modification of electronic structure at the interfaces can lead to consequent change of magnetic properties. Research on the SL-period dependent evolution of the interfacial electronic structure will be informative for understanding the role of interfacial 3-5 coupling to affect the Ir magnetism, which has not been systematically investigated so far. For this purpose, we fabricated LaMnO3/SrIrO3 (LMO/SIO) SLs with different repetition periods. X-ray magnetic circular dichroism (XMCD) were employed to study how the SL period and the interfacial coupling can affect the properties of FM Ir moments. X-ray absorption spectra (XAS) and hard x-ray photoemission spectroscopy (HAXPES) were carried out to characterize the SL-period dependence of the electronic structure. A systematic SL-period dependent trend of the ratio between orbital and spin magnetic moments of Ir as well as the electronic structure were observed. First-principles calculations demonstrate a satisfactory consistency with the experimental results and reveal that the formation of the interfacial Ir-Mn molecular orbital associated with concomitant electronic-structure change is the pivotal mechanism behind this interfacial coupling.
Methods
[(LMO)a/(SIO)a]b ((,)=(1,24), (2,12) and (8,3), where is counted in unit cells, SL in abbreviation) SL samples as well as LMO and SIO reference samples (24 unit cells) were fabricated by laser molecular beam epitaxy. A KrF excimer pulsed laser (=248 nm) with a repetition rate of 2 Hz and an energy density of 1.5 J/cm2 was employed. The sample temperature and ambient oxygen pressure were controlled at 720 oC and 16 Pa, respectively. The distance between the stoichiometric LMO or SIO targets and STO(001) single crystal substrates was set at 6 cm. The layer-by-layer growth of the SLs was guaranteed by monitoring the oscillation of the reflection high-energy electron diffraction signal (as reported in Ref. [26_Huang]). The sample structure is schematically displayed in Fig. 1(a). The crystal structure of the SLs was characterized by an x-ray diffractometer with Cu radiation (XRD, Rigaku RINT-2200). Basic magnetic properties of the SLs were characterized by a superconducting quantum interference device (SQUID, Quantum Design).
The Ir edge XAS/XMCD measurements were conducted at BL39XU of SPring-8. A He-flowing cryostat was used to cool the samples to a lowest temperature of 30 K. In-plane magnetic field up to 2 T along the x-ray propagation was applied by an electromagnet. The Ir edge XAS/XMCD spectra were collected by standard helicity reversal technique [28_Suzuki] with a grazing incidence geometry (5.5o incidence angle) and partial fluorescence yield (PFY) mode. For PFY detection of XAS/XMCD at the Ir and edges, Ir and emissions were collected and energy-analyzed respectively by a four-element silicon drift detector (Sirius 4, SGX Sensortech Inc.). The positive magnetic field direction is defined opposite to the x-ray propagation. HAXPES measurements were carried out at BL47XU of SPring-8. The incidence angle of 7.94 keV hard x-ray was set at 1o and the emitted photoelectrons were detected by a Scienta R-4000 electron energy analyzer, whose energy resolution was 280 meV.
First-principles calculations were carried out within the framework of density functional theory (DFT) [29_Hohenberg, 30_Kohn] using the generalized gradient approximation (GGA) in the parameterization of Perdew-Burke-Ernzerhof (PBE) format exchange-correlation functional [31_Perdew], as implemented in the Vienna Simulation Package (VASP) [32_Kresse, 33_Kresse, 34_Kresse, 35_Kresse]. SOC is implemented in the projector augmented wave (PAW) method [36_Blochl, 37_Kresse] which is based on a transformation that maps all electron wave functions to smooth pseudowave functions to describe the interaction between electrons and ions. The corresponding electronic configurations for each element are Sr: 445; Ir: 56; O: 22; La: 5565; Mn: 334. The cutoff energy is set to 500 eV. To account for strong correlation effects [38_Liechtenstein], we included the Hubbard correction for Ir and Mn states with Ir5d=2 eV and Mn3d=3 eV [39_Caviglia, 27_Okamoto]. We used K-points following the Monkhorst-Pack scheme in our systems. The convergence criterion for the electronic relaxation is eV. In this calculation, we relaxed the SL cell parameters and atomic positions with the in-plane lattice constant constrained to that of STO. The doubled unit cell has been used with the experimental in-plane lattice constant of STO, ==3.905 Å (see Fig. 7 in the Appendix A). Optimized SL structures were achieved when forces on all the atoms were 0.01 eV/Å. We calculated SL11, SL22, SL33 and SL44 rather than SL11, SL22 and SL88 investigated in our experiments since the supercell of SL88 is too large for DFT-based calculations.
Results and Discussions

The XRD scan spectra in (0 0 ) direction of the SLs are presented in Fig. 1(b). SL satellite peaks can be clearly observed beside the (0 0 1) diffraction of the STO substrate, confirming the high quality of the SLs. In-plane magnetization-field (-) and magnetization-temperature (-) curves shown in Fig. 1(c,d) indicate FM behaviors of all the three SLs, as well as the pure LMO reference sample. The FM Curie temperatures () are around 150 K, 170 K and 200 K for SL11, SL22 and SL88, respectively. The saturated magnetization and of SLs change systematically with the SL period, which should be attributed to the enhancement of Mn-Mn electron hopping within LMO layers as the SL period increases. The reference SIO sample remains paramagnetic down to 10 K as previously reported [15_Cao].

To detect and comprehensively analyze the SL-period dependence of the Ir magnetism, hard x-ray XAS/XMCD measurements at the Ir edges were conducted and the main results are displayed in Fig. 2. The XAS (defined as and averaged for positive and negative magnetic fields, where and are PFY-XAS measured by x-rays with left and right helicities) in Fig. 2(a,b) of all the SLs have similar line shape with white-line peaks at around 11.223 keV() and 12.834 keV(), indicating a localized feature of Ir 5 states. No obvious peak shift induced by interfacial charge transfer between Ir and Mn was observed, which is consistent with previous reports [40_Zhong]. Ir-Mn charge transfer will lead to deviation of the Ir valence from the nominal valence state of 4+ and add complexity to the investigation. So we chose pure LMO to construct the SLs since it is reported that with the increase of the Sr%, significant charge transfer between Ir and Mn will appear in LSMO/SIO SL system [22_Nichols, 27_Okamoto]. The XAS peak intensity slightly changes in different SLs, which suggests possible modification of the spin-orbital states of Ir and will be detailedly discussed in the latter parts.
The XMCD signal in Fig. 2(c,d) is defined as (averaged for positive and negative magnetic fields). All three SLs exhibit clear XMCD signal at 30 K. The XMCD peak positions are located at 2 eV lower than the XAS white-line peaks (about 11.221 keV for and 12.832 keV for ), which is a sign of intrinsic XMCD signal mainly originating from states at lower energy, rather than artifacts from the XAS measurement. The positive XMCD signals at both the and edges and the larger intensity of the XMCD peak indicate that the net magnetization of Ir is antiparallel to the external field and antiferromagnetically coupled to the Mn magnetization as previously reported [22_Nichols, 23_Yi, 24_Yi, 25_Kim]. The XAS/XMCD spectra of SL11 measured at temperature above the FM transition (200 K) are also displayed in Fig. 2. XMCD signal is negligible at both the and edges, evidencing the disappearance of FM ordering of Ir moments.
Temperature dependent XMCD measurements of SL11 at the Ir edge (Fig. 2(e)) shows that the integrated XMCD intensity decreases with the increase of temperature and vanishes at 150 K, in consistency with the of the SL11 sample. The - curve of SL11 is also plotted into Fig. 2(e) for comparison. Since the total magnetization is dominated by the Mn moments, this result confirms that the emergent FM Ir moments originate from the interfacial coupling and rely on the existence of Mn FM ordering. Fig. 2(f) shows the magnetic field dependence of Ir XMCD signal. Clear FM behaviors were observed for all the SLs.

Based on the successful observation of FM XMCD signal at the Ir edges, sum-rules analysis can provide important information about the local properties of the Ir magnetic moments. By applying sum-rules analysis with the following formulas [41_Thole, 42_Carra], it can be obtained that at 30 K the orbital magnetic moments of Ir () are -0.0235, -0.0170 and -0.0055 B/atom, and the effective spin magnetic moments of Ir () are -0.0190, -0.0125 and -0.0024 B/atom, for SL11, SL22 and SL88 respectively. With the assumption of the negligible charge transfer between Ir and Mn according to the previous theoretical report [40_Zhong], Ir has a nominal valence state of 4+ in our SL system, so the number of 5 holes is estimated as .
Fig. 3(a) displays the sum-rules analysis results of our samples as well as previous reports of other LSMO/SIO SL and magnetic iridate systems [22_Nichols, 23_Yi, 25_Kim, 43_Haskel, 44_LagunaMarco]. The size of measured Ir FM moments is quite consistent with the previous reports [22_Nichols, 23_Yi, 24_Yi, 25_Kim] that Mn moments dominate the total magnetization and Ir magnetization is 1 or 2 orders of magnitude smaller than Mn. Since FM Ir moments are antiparallel to the external field, the signs of the Ir magnetic moments are negative, which is opposite to the perovskite iridates. The size of the Ir moments decreases with the SL period, which suggests that the FM Ir moments mainly distribute near the interfaces (see more detailed comments on this point in the Appendix B).
Since Ir moments often exhibit canted AFM ordering in perovskite iridates [45_Kim, 21_Matsuno, 25_Kim], the net Ir moment size depends on both the absolute size of the local Ir moment and the canting angle between moments in different AFM sublattices. Hence the net moment size evaluated by XMCD sum-rules analysis varies in different systems (Fig. 3(a)). While the / ratio can be compared among different systems (Fig. 3(b)) and reflects the local properties of the Ir moments. First it can be noticed that / ratio of SL11 is quite consistent with a previous report of the same SL [24_Yi], indicating a satisfactory reproducibility of the XMCD measurements. Moreover, the / ratios of LSMO/SIO SL systems are generally much larger than magnetic iridates such as Sr2IrO4 and BaIrO3. This should be attributed to the different origin of Ir moments. The FM Ir moments mainly originate from the interfacial coupling with Mn moments in LSMO/SIO SLs while in magnetic iridates they mainly originate from the electron correlation of the eff=1/2 band [11_Kim].
The / ratio of ideal eff=1/2 model is 2 [11_Kim], as indicated in Fig. 3(b). The contributions of both and magnetic dipole term are included in the (), but they can not be easily separated by experiment. To estimate the / ratio of our samples and compare with the ideal eff=1/2 model, we used the estimation of =0.18 (similar to the values obtained by cluster-model calculations for Sr2IrO4 [43_Haskel] and BaIrO3 [44_LagunaMarco], is the spin angular momentum) to evaluate the term. The estimated / ratio of our samples are roughly consistent with the eff=1/2 scheme.
Interfacial Ir-Mn coupling should significantly affect the local properties of the Ir moments, which can be evidenced by the SL-period dependence of the / ratio of Ir. The / ratio systematically decreases with decreasing the SL period. The decrease of the SL period leads to the enhancement of the interfacial coupling effects since the volume ratio of interfacial layers increases with the decrease of the SL period. In other words, the / ratio can be effectively decreased by the interfacial coupling. Remarkably, our DFT calculation results (Fig. 3(c)) also show a similar trend that / ratio systematically decrease with the decreasing of the SL period, which is in good agreement with the experiments. Although the calculated absolute values of / ratio deviates from the experimental values, the consistent trend indicates that the local properties of Ir magnetic moments indeed vary systematically with the SL period. It is worth mentioning that the calculated is more strongly dependent on the SL period than . The strong SL-period dependence of should be attributed to the orbital reconstruction at the interface, which can make the originally antiparallel electron spins of Ir align parallel. This point will be discussed in detail in the following parts.

To understand this SL-period dependent behaviors of the Ir magnetic moments, the electronic structure of the SLs should be investigated in detail. Fig. 4(a) shows schematic orbital energy level of the LMO/SIO SL system [27_Okamoto, 46_Okamoto]. The octahedral crystal field splits both the Ir 5 and Mn 3 levels into and states. According to the relaxed crystal structures in our DFT calculations, the O-Ir-O bond is compressed while O-Mn-O bond is elongated along the -axis in the SLs. Therefore in SIO layer the orbital lies above the orbital and vice versa in LMO layer due to the Jahn-Teller effect. The orbitals of Ir and Mn can hybridize with each other along the -axis and form molecular orbitals [25_Kim, 27_Okamoto]. The formation of the molecular orbitals can be visualized by the calculated charge density difference in Fig. 4(b). Electrons are spatially redistributed from Ir/Mn atoms to the interfacial region near the O atoms. The physical picture of the interfacial molecular orbital can also be evidenced by the partial density of states (PDOS), as shown in Fig. 5. With the decrease of SL period, the orbital (mainly located above 1 eV) of Ir obviously shifts to higher energy above the Fermi level (), while the orbital of Mn exhibits a trend of PDOS redistribution from higher to lower energy. The SL-period dependence of PDOS is consistent with the formation of molecular orbital that anti-bonding molecular orbital is mainly contributed by Ir and lies at higher energy than the original orbital of Ir, while bonding molecular orbital is mainly contributed by Mn and lies at lower energy than the original orbital of Mn, as schematically shown in Fig. 4(a).

Experimental results also show some clues to understand the SL-period dependent electronic structure of the SLs. Details in XAS at the Ir edge can provide information about the unoccupied Ir 5 states. By taking a closer look at the white-line regions of the XAS (Fig. 6(a,b)), the white-line intensity at the edge () obviously increases with the decrease of the SL period, while the white-line intensity of the edge () keeps nearly constant and slightly decreases for SL11. This variation of the white-line intensity induces the systematic change of the branching ratio () and the expectation value of SOC operator (Fig. 6(c)). Here , [47_Laan] and . It can be observed that the of SL11 is 4.61, which is similar to that of previous reported (LSMO)1/(SIO)1 SL () [25_Kim]. The and value systematically increase with the SL period. Since the edge corresponds to the electric dipole transition from 2 to 5 states while the edge corresponds to the electric dipole transition from 2 to 5 states, the decrease of indicates less occupation of 5 states and/or more occupation of 5 states. In perovskite iridates, the octahedral crystal field splits the 5 () states into states and eff=1/2 states, and eff=3/2 states originate from the atomic 5 () states [11_Kim, 44_LagunaMarco]. The change of indicates that the occupation of eff=3/2 states is decreased by the interfacial coupling and the occupation of states or eff=1/2 increases when the interfacial coupling is present.
We conducted bulk-sensitive HAXPES measurements to further characterize the valence band structure of these SLs (Fig. 6(d)). It can be clearly observed in the inset of Fig. 6(d) that a feature at a binding energy of 1.2 eV is enhanced with the decrease of the SL period, which should be the interfacial-coupling-enhanced occupation of the bonding molecular orbital. Simultaneously, the intensity near the decreases with the decrease of the SL period, showing clear evidence that some density of states (DOS) near the was transferred to deeper levels due to the Ir-Mn interfacial coupling.
Enhanced occupation of the states rather than eff=3/2 states induced by the interfacial coupling is more likely due to the following reasons. As displayed in the schematic in Fig. 4(a), the bonding molecular orbital appears below the and changes the relative occupation of different orbitals. In maganite/iridate SL systems, rather than an ideal eff=1/2 scheme, mixed occupation of both eff=1/2 and eff=3/2 states can often occur [23_Yi]. When the molecular orbital is formed by interfacial coupling, some of the eff=1/2 and eff=3/2 electrons of Ir near the will be transferred into the bonding molecular orbital. In particular, the electron transfer from eff=3/2 states to the molecular orbital will lead to the change of and consequently the local property change of the Ir moments. One may argue that when the SL period increases, the change of Ir 5 bandwidth induced by the dimensionality of SIO layers may also account for the electron redistribution among eff=3/2 and eff=1/2 states and give rise to the change of . However, as the SL period increases, the enhanced Ir-Ir hopping will result in simultaneous increase of the bandwidth of both eff=3/2 and eff=1/2 states. Since eff=3/2 states are nearly fully occupied, the center of eff=3/2 states lies deep below the . While the center of eff=1/2 states lies close to the , as schematically shown in Fig. 4(a). With the center of mass of the states fixed, widening of the eff=3/2 states should lead to more DOS above and less occupation of itself, which does not agree with the change of XAS white-line intensity. Consequently, the experimentally observed change should be mainly attributed to the electron transfer between eff=3/2 states and the states, but not simply the redistribution of electrons within the (eff=3/2 and eff=1/2) states.

Up to now, we illustrated the effect of interfacial Ir-Mn coupling on the electronic structure in LMO/SIO SL system. Based on this, we can discuss the consequent effect on the / ratio. Since eff=3/2 states originate from the atomic 5 () states [11_Kim], whose spin and orbital moments are antiparallel. While eff=1/2 and states originate from the atomic 5 () states [11_Kim], whose spin and orbital moments are parallel. When Ir-Mn interfacial coupling transfers some electrons from eff=3/2 states to states, the SOC sign of these electrons are effectively changed. Spin and orbital components of the Ir moments become more parallel. As a result, we observed smaller / ratio in SLs with shorter SL period, in which the interfacial coupling is more efficient. From another angle of view, the consistent trend of / ratio in DFT calculations further reveals that is more sensitive to the interfacial coupling than . This is due to the fact that spin is directly carried by the redistributed electrons from states to states. The shorter SL period, the more electrons transferred from Ir to molecular orbitals, and the more Ir spins which are originally antiparallel to the total Ir spin tend to reverse its direction and align parallel to the total Ir spin due to the Hund’s coupling [27_Okamoto]. On the other hand, the is relatively robust to the interfacial coupling, which drives the / ratio to be smaller when the SL period decreases.
Conclusions
In conclusion, by comprehensive experimental and theoretical investigations of the LMO/SIO SL system, we conclude that the local properties of Ir FM moments and interfacial electronic structure can be modified by the SL period. Hybridization between the orbitals of Ir and Mn along the -axis of the SLs can form a bonding molecular orbital which lies below the , so that electrons from Ir states near the are pulled down into this bonding molecular orbital. The / ratio of Ir can be modified by the interfacial coupling due to the electron transfer between Ir and orbitals, which have different spin alignments. Our results demonstrate a clear physical picture of the Ir-Mn interfacial coupling in manganite-iridate SL system. The conclusions could also be generalized to other TMO-based perovskite heterostructures and SL systems, such as SrRuO3/SrIrO3 heterostructures [18_Matsuno], etc.
Acknowledgments
This work was supported by Grant-in-Aid for JSPS fellows (No. 17F17327) and JSPS KAKENHI (Grant No. 19H05824), Japan, as well as the Natural Science Foundation of China (Grant No. 51729201) . The synchrotron radiation experiments at SPring-8 were performed under the approval of the Japan Synchrotron Radiation Research Institute (Proposal No. 2018A1232, 2018B1449 and 2019A1239). X. R. Wang acknowledges supports from the Nanyang Assistant Professorship grant from Nanyang Technological University and Academic Research Fund Tier 1 (RG108/17 and RG177/18) from Singapore Ministry of Education. L. Shen acknowledges the support from Singapore MOE Tier 1 (Grant R-265-000-615-114). We also acknowledge the enlightening discussion with Prof. M. van Veenendaal, as well as the experimental supports and discussions provided by Dr. K. Ishii and Dr. Y. Takeda.
Appendix A
Fig. 7 shows the detailed lattice structure used for the DFT calculations. SL33 is taken as an example. Our theoretical approaches have been validated by comparing with previously reported experimental benchmarks [27_Okamoto, 48_Elemans, 49_Lee].
Appendix B
To get a feeling about the spatial range that the interfacial coupling can influence, we averaged the magnetic moments obtained by sum-rules analysis to each interface. As shown in Fig. 8(a), the size of Ir magnetic moment per interface increases with the SL period, which indicates that the effect of the interfacial coupling may not be restricted only in the unit cells adjacent to the interface, especially by comparing SL22 and SL88. On the other hand, the difference between SL11 and SL22/SL88 can be attributed to other factors. As depicted in Fig. 8(b), in SL11 every SIO unit cell is sandwiched by two LMO unit cells, and in SL22 and SL88, every SIO unit cell is only adjacent to one Ir-Mn interface. So in SL11, every Ir orbital is directly hybridizing with two Mn orbitals simultaneously. While in SL22 and SL88, every Ir orbital has only one counterpart for interfacial hybridization. This factor may cause different electronic structure and local properties of Ir moments at the interface.

