Long-range magnetic coupling across a polar insulating layer.

Magnetic interactions in solids are normally mediated by short-range exchange or weak dipole fields. Here we report a magnetic interaction that can propagate over long distances (∼10 nm) across a polar insulating oxide spacer. Evidence includes oscillations of magnetization, coercivity and field-cooled loop shift with the thickness of LaAlO3 in La0.67Sr0.33MnO3/LaAlO3/SrTiO3 heterostructures. Similar modifications of the hysteresis loop appear when two coupled films of La0.67Sr0.33MnO3 are separated by LaAlO3, or another polar insulator, but they are absent when the oxide spacer layer is nonpolar. The loop shift is attributed to strong spin-orbit coupling and Dzyaloshinskii-Moriya interaction at the interfaces. There is evidence from inelastic light scattering that the polar spacer mediates long-range transmission of orbital magnetization. This coupling mechanism is expected to apply for any conducting ferromagnetic oxide with mixed valence; in view of electron hopping frequency involved, it raises the prospect of terahertz tunability of magnetic coupling.

Exchange interactions in insulators usually involve virtual transfer of spin-polarized electrons in a process known as superexchange. The magnetic interaction between an atom with spin Si at the origin and another Sj at distance r is described by the Heisenberg Hamiltonian =−2Jij Si·Sj, where the exchange constant Jij are positive for ferromagnetic interactions and negative for antiferromagnetic interactions. Exchange interactions in oxides are often negative1. Since Jij depends on the overlap integrals between the wave functions of neighbouring atoms, superexchange may be relatively strong (JS2≤10 meV) for nearest-neighbour cations, but the interaction falls off rapidly when r>0.5 nm (ref. 2).

Longer-ranged exchange is possible when the spin-polarized electrons are delocalized. Examples are metallic oxides, where ferromagnetic direct or double exchange is important3, and thin-film structures such as spin valves, where two ferromagnetic metallic layers are separated by a metallic spacer and the exchange interaction between them oscillates in sign with spacer thickness4. For oxide-based systems, oscillatory exchange coupling as a function of the thickness of the spacer has been found, for example, in superlattices involving Fe3O4 and TiN (ref. 5). The TiN spacer layers are conducting, and the coupling is attributed to oscillatory spin polarization of the magnetite conduction electrons as a function of TiN layer thickness. The theory of exchange coupling in thin-film sandwich structures with conducting or insulating spacers has been developed by Bruno6.

Here we are concerned with the coupling of films of ferromagnetic conducting oxides across relatively thick layers of insulating oxide, where exchange is not anticipated. The spacers are (001) polar films of LaAlO3 (LAO) or NdGaO3 (NGO) where the formal layer charges are ±e, with e the elementary charge, and (001) SrTiO3 (STO), which is formally nonpolar since the structure is a stack of uncharged layers of SrO and TiO2. The ferromagnetic mixed-valence manganite La0.67Sr0.33MnO3 (LSMO) is also polar with a formal layer charge of ±0.67e, although the charge may be screened by the Mn conduction electrons as the Thomas–Fermi screening length, lTF=0.17 nm for carrier density ne=1 × 1027 electrons per m3, is comparable to half a unit-cell layer thickness, 0.19 nm.

Results

The configuration and growth of the heterostructures

Our heterostructures usually involve a thickness parameter tf=10 nm film of LSMO on top of a polar or nonpolar oxide spacer. On the lower side, there may be a similar LSMO layer, or else a quasi two-dimensional electron gas (2DEG) induced in an STO substrate when the spacer is LAO7. The upper LSMO can be replaced by a different ferromagnetic oxide, SrRuO3 (SRO) or La0.67Sr0.33CoO3 (LSCO). It is convenient to express the layer thickness t in these heterostructures as a multiple n of the cubic perovskite unit cell (u.c.) parameter a0=0.39 nm. The heterostructures grow epitaxially, with little interdiffusion of the constituents (Supplementary Figs 1 and 2)

The 2DEG that forms at the LAO/STO interface has been shown to exhibit magnetic order891011. The coexistence of ferromagnetic, and paramagnetic or diamagnetic phases below 100 K is associated with nanoscale electronic phase separation8. Ferromagnetism persists up to room temperature in samples grown at high oxygen pressure. Sufficient conditions for electronic surface reconstruction and emergence of the 2DEG are that the STO should be TiO2 terminated and that the LAO should be >4-u.c. thick12.

Magnetic coupling between LSMO and 2DEG at LAO/STO

Figure 1 shows the 10 K M–H hysteresis loops of 10 nm LSMO on LAO/STO, including the effect of field cooling in an in-plane field of 7 T. All the magnetic parameters (spontaneous magnetization Ms, coercivity Hc and the horizontal field-cooled shift of the hysteresis loop Hshift) depend on the thickness of the LAO spacer layer, provided it exceeds the threshold n=4. There are no changes when n≤4. This demonstrates that the charge transfer to the underlying STO that creates the 2DEG also modifies the magnetism of the LSMO. The magnetization of a 10-nm layer of LSMO prepared without LAO, n=0, is Ms=350 kA m−1 or 2.2 μB/Mn, where μB is Bohr magneton, compared with the bulk value of 3.6 μB/Mn. The reduction of the magnetization in thin (100) LSMO films13 is due to a noncollinear ferromagnetic structure. The loop shift when n≥4 is reminiscent of exchange bias at a ferromagnetic/antiferromagnetic interface14, although it can be of either sign, depending on the spacer thickness. When the 7 T cooling field is reversed, so is the M–H loop shift (Fig. 1a). The maximum value μ0Hshift=28 mT corresponds to an asymmetric coupling energy σa=μ0HshiftMstf of 0.11 mJ m−2. Moreover, the sign of the shift changes from negative to positive when n=11, and back again when n=16. The oscillations in magnetization, with a maximum at n=7 and minimum at n=12 are reminiscent of the oscillations in exchange found in spin valve structures as a function of the spacer layer thickness4, except that the period of the oscillatory behaviour in our insulating system is three to eight times greater than that found for Co/Cr/Co spin valves or magnetic tunnel junctions with an MgO barrier thinner than 3 nm (refs 15, 16). None of the explanations proposed for metallic spacers such as Ruderman–Kittel–Kasuya–Yosida or other exchange coupling via interlayer conduction electrons17, Fermi velocity mismatch18, hybridization of the interlayer conduction electrons with the interface d-states of the magnetic films19 applies to an insulating spacer. Nor can tunnelling or dipole interactions associated with correlated surface roughness (Supplementary Fig. 3) account for the magnitude of the effect20, or its change of sign. We propose a new magnetic interaction mechanism.

Figure 1

Oscillatory behaviour in La0.67Sr0.33MnO3/LaAlO3/SrTiO3 heterostructures.

(a) At 10 K, hysteresis loops with different LaAlO3 thickness after 7 T field cooling; a reverse case is shown in second panel. (b–d) The horizontal shift of the M–H loop (Hshift), coercivity (Hc) and saturation magnetization (Ms) as a function of insulator thickness, measured as a multiple n of the unit cell (uc) parameter (0.39 nm). The solid lines in b–d are all fits to a exponentially damped sinusoidal oscillation with period n−9 u.c. and attenuation length n=9 u.c.

Transport properties of 2DEG under the magnetic coupling

The reciprocal interaction of LSMO with the 2DEG across the LAO spacer is demonstrated using a 3-μm strip of LSMO laid across the 5 × 5 mm2 (001) STO substrate. The conduction of the 2DEG was measured either parallel or transverse to the LSMO strip with four in-line contacts (Fig. 2). The insulating character of the LAO spacer was confirmed by direct electrical measurement. In the parallel case, the usual temperature dependence of conductance was found but, in the transverse geometry, a resistance minimum appears due to a marked increase of scattering at low temperature. Such minima are associated with magnetic (Kondo) scattering21, and are accompanied by negative magnetoresistance (Supplementary Note 1; Supplementary Fig. 4). There is no perceptible diffusion of Mn across the LAO barriers (Supplementary Fig. 2), so we are seeing an effect of the distant LSMO layer on the 2DEG, mediated by the LAO. It is the counterpart of the effect of the 2DEG on the LSMO hysteresis loop, shown in Fig. 1. The mutual interaction of the 2DEG and the LSMO is confirmed by gating experiments2223 which show striking changes in loop shift, coercivity and resistance minimum with bias voltage (Supplementary Note 2; Supplementary Fig. 5). The changes of both magnetization and transport properties on varying the LAO thickness or gate potential indicate correlated charge transfer on opposite sides of the LAO. Electron transfer to the 2DEG, which improves its conductivity and increases the magnetic scattering, is coupled with hole depletion in the LSMO2425, which increases its ferromagnetic moment. The magnitude of the moments on either side are coupled, but the relative orientations cannot be inferred because the 2DEG moment is only ∼1% of that of the LSMO, and too small to measure directly.

Figure 2

Transport properties of the 2DEG at the LaAlO3/SrTiO3 with La0.67Sr0.33MnO3 strip.

(a–c) The resistance versus temperature for 10, 13 and 17 uc LaAlO3 spacers, measured in the two configurations illustrated in d. (e) The temperature of resistance minimum measured in configuration (1) as a function of LaAlO3 thickness.

An estimate of the symmetric coupling energy between the ferromagnet and the 2DEG can be derived from the maximum change of LSMO coercivity μ0ΔHc when the 2DEG forms. It is ∼40 mT. The change in anisotropy field μ0ΔHa will be roughly an order of magnitude greater. The corresponding change of anisotropy energy ΔK1=μ0MsΔHa/2 expressed per unit area is σs≈0.8 mJ m−2 or 9 Ka0−2. Since 0.5 electrons per unit cell area at the STO interface are needed to avert the polarization catastrophe, and <0.5 are ferromagnetic26, the energy per electron is >18 K. The additional scattering in the transverse strip geometry of Fig. 2 sets in above this temperature. A plot of the temperature of the resistivity minimum as a function of LAO thickness reveals an oscillatory variation quite similar to that of Hshift or Ms (Fig. 2e). The curves do not overlap at zero bias, but their similarity reinforces the idea that the LAO is transmitting magnetic information from the LSMO to the 2DEG.

The effect of polar or nonpolar insulator

To investigate the magnetization further, sandwich structures consisting of two 10 nm films of LSMO separated by a layer of LAO, NGO or STO of varying thicknesses were prepared to compare the effects of polar and nonpolar barriers. Results shown in Fig. 3 for a polar insulator, LAO or NGO, resemble those already described for the LSMO/LAO/STO heterostructures, except there is now no threshold thickness. The two 10 nm films of LSMO switch together, although thinner films show evidence of independent switching (Supplementary Note 3; Supplementary Fig. 6), which demonstrates that they are not identical. The M–H loop again varies with barrier thickness in an oscillatory manner. The behaviour with a nonpolar STO spacer is quite different. Then, there are no oscillations and no loop shift, only a monotonic increase of coercivity and decrease of magnetization with n. To mediate oscillatory coupling, it seems that the spacer must be polar. The essential differences between polar and nonpolar oxide layers are the possibilities of charge transfer with electronic interface reconstruction, and oscillations or standing waves involving the sheets of positive and negative charge, that distinguish the (001) direction normal to the sheets in a thin film from the in-plane directions.

Figure 3

Heterostructures with polar and nonpolar insulators.

(a) LaAlO3. (b) NdGaO3 and (c) SrTiO3. (d–f) The variation of horizontal shift of the M–H loop (Hshift), coercivity (Hc) and saturation magnetization (Ms) as a function of insulator thickness. The data are measured at 10 K.

Discussion

We associate the loop shift with the interface charge that is responsible for the variation of the moment with spacer thickness2425. The symmetry-breaking interaction is likely to be the Dzyaloshinskii–Moriya (D–M) interaction Dij·(Si × Si), where Dij is the D–M constant, associated with enhanced spin–orbit interaction at the charged interface27. Orbital currents associated with geometric frustration of triangles of spins at the manganite interface28 may be involved. The shift oscillates as a function of LAO thickness with the oscillating interfacial charge, being greatest when electron depletion at the LAO/LSMO interfaces is the greatest and the most screening charge appears in the LSMO2425. The direction of the loop shift depends on both the sign of the interfacial charge and the direction of magnetization of the LSMO. If the sign of the charge at the two LAO interfaces is opposite, Dij is directed in the same sense perpendicular to the interfaces. The interface charging effects are maximum at about 5 u.c. and minimum near 10 u.c., as seen in Fig. 3, but the effects are weaker than those for a single LSMO layer and the strongly polar LAO/STO interface that we discussed first. The opposite sign of the charge at the two LSMO interfaces is attributed to the different growth conditions of the upper and lower LSMO films.

The fact that we never observe a loop shift for a single 10 nm LSMO layer in the absence of interface charge, no matter how we prepare it (Supplementary Fig. 6), points to the existence of some mechanism that couples not only the magnitudes, but also the directions of the LSMO moments across the polar spacer. Two mechanisms for transmitting magnetic interactions can be envisaged. We ruled out dipolar interactions due to correlated interface roughness because they are much too weak (Supplementary Table 1; Supplementary Fig. 3). A long-range exchange interaction is implausible. Any exchange between partly filled d-orbitals in the two ferromagnetic layers due to virtual charge transfer involving the filled 2p band should be very weak and ferromagnetic, regardless of the polar nature of the oxide; coupling by spin-polarized electron tunnelling is expected to decay exponentially over ∼2–3 u.c. (refs 6, 15, 16).

An alternative idea is a coupling mechanism that is not based on spin, but depends on the transmission of orbital correlations across the polar barrier. Electron hopping in mixed-valence3 manganites is a pair-localized process involving the fourth, majority spin Mn eg electron, with an associated polaronic lattice distortion29. The hopping frequency f is of order 1013 Hz (refs 3, 30), and an orbital moment morb≈πd2ef/8 is created in LSMO as the electron hops back and forth in the exchange field, where d is the Mn–Mn separation, a0/√2=0.28 nm. Planar charge excitations in the LSMO or LAO can be decomposed into two counter-rotating orbital motions, one of which could resonate with the electron hopping and couple to the Mn spins in LSMO via spin–orbit interaction. The asymmetric coupling energy of 0.11 mJ m−2 is equivalent to ∼1 K per interface Mn. Taking S=1.85 for Mn in LSMO and the spin–orbit coupling constant for Mn as Λ=150 K in the interaction ΛLS with a, we deduce L=3.7 × 10−3, where L and S are the orbital and spin angular momentum, respectively. This value of L corresponds to the orbital moment morb/μB induced by electron hopping at a frequency of 7 THz.

Some further experiments were conducted to validate these ideas. First, the upper LSMO was replaced by SRO, which is a nonpolar ferromagnetic metal that does not exhibit pairwise electron hopping. In that case, there was no loop shift on field cooling for 5, 10 or 16 u.c. of LAO (Fig. 4), and the coercivity of the bottom LSMO layer remained essentially constant at 10–13 mT. Second, the top manganite electrode was replaced by a similar cobaltite (LSCO) electrode. The sign of the spin–orbit coupling for cobalt Λ=−272 K is opposite to that of manganese, so we would expect the sign of the loop shift to change. This is exactly what we observe (Fig. 4).

Figure 4

The loop shift for a 10 nm La0.67Sr0.33MnO3 bottom layer with different top layers.

There is no shift for SrRuO3, which is a nonpolar metal. The shift for La0.67Sr0.33CoO3, which is polar but with the opposite signs of spin–orbit coupling, is opposite to that of La0.67Sr0.33MnO3.

To shed some light on the role of the oxide in transmitting magnetic coupling, we have studied the effect of a magnetic field on inelastic light scattering. Results shown in Fig. 5b–d are for a LAO crystal, where a number of both Raman and luminescent modes are found to exhibit Zeeman splitting. Figure 5e compares the temperature dependences of the loop shift and the intensity of the magnetic field-sensitive lines in LAO. They are similar, especially for the defect-related luminescent transitions at 272 and 290 cm−1, and the Raman excitations at 1,100 cm−1, and its harmonics which the field splits into three, corresponding to g=1.2±0.2. As shown in Supplementary Figs 7 and 8, the observation of similar field-sensitive excitations in other oxides correlates with the observation of a loop shift when the oxide is used as a polar spacer. The requirements for a suitable spacer oxide are that it should include a polar plane containing a heavy (rare earth) cation and oxygen, but no light cation (Sr, Al…) (Supplementary Note 4). These observations support the idea that magnetic-field-dependent excitations in a polar oxide are associated with the transmission of magnetic information across the insulating polar layer; this could involve electrons associated with Schottky defects or metal in gap states31 in LAO, which couple to electrons hopping at the Fermi level of LSMO. The phase of these states, possibly influenced by standing waves in the LSMO6 controls the oscillatory behaviour.

Figure 5

Inelastic light scattering from an LaAlO3 crystal at 5 K under magnetic field.

(a) The zero magnetic field spectrum of LaAlO3. (b–d) The field-dependent splitting of the spectral lines up to 9 T. R denotes a Raman mode and L a luminescent transition. The transitions in d are third order of Raman lines in c. (e) The coincident temperature dependence of the magnetization oscillations including horizontal shift of the M–H loop (Hshift) and the intensity of the magnetically split components in the luminescence/Raman spectrum.

If the interface charge is critical, why is the oscillatory effect manifested only below 80 K, although the 2DEG is stable at room temperature? The onset temperature reflects the interlayer coupling that is mediated by orbital excitation via the LAO layer, and establishment of the phase of the gap states in the LAO. The appearance of the loop shift and magnetic oscillation only below 80 K is evidence for a role of the magnetic-field-dependent Raman and luminescent modes in the polar layer in somehow providing a coupling mechanism between the two ferromagnetic layers.

In conclusion, we have set out a range of evidence for magnetic coupling that acts across relatively thick polar layers of a nonmagnetic insulating oxide. The magnitudes of the moments are coupled by the correlated interface charge on either side of the barrier, and there is evidence that the directions of the moments are also aligned. The proposed coupling mechanism is not via interlayer exchange or dipolar interactions, but it involves orbital magnetism, which is excited at surface or defect-related optical modes in the oxide spacer. We have shown that they couple both to LSMO and to the 2DEG at the STO interface, but the modes do not couple to the nonpolar metallic ferromagnet SRO nor do they propagate in a nonpolar insulator.

Methods

Sample fabrication

For the LSMO/LAO/STO heterostructures, first, LAO was deposited on a TiO2-terminated STO substrate under 5 × 10−2 torr at 780 °C using pulsed laser deposition. The thickness of LAO was monitored by in situ reflection high-energy electron diffraction, which revealed layer-by-layer growth. Then, 10 nm of LSMO was grown under the same oxygen pressure at 720 °C. For the LSMO/oxide barrier/LSMO sandwich heterostructures, (La0.3Sr0.7)(Al0.65Ta0.35)O3 (LSAT) substrates were first annealed at 1,050 °C in air for 2.5 h to obtain atomically flat surfaces. Then, the LSMO (10 nm)/oxide barrier/LSMO (10 nm) heterostructures with different barrier thickness were prepared under 5 × 10−2 torr oxygen pressure at 720 °C. The same growth condition was also applied for SrRuO3 and La0.67Sr0.33CoO3 layer deposition. LSAT was used as the substrate because it is lattice matched with LSMO, but similar results were obtained with LAO, STO and NGO substrates (Supplementary Note 5; Supplementary Fig. 9). Exact same fabrication condition was applied in the reproducibility test (Supplementary Note 6; Supplementary Fig. 10).

Interface structure and elementary diffusion

The interface structure was studied by high-resolution transmission electron microscopy, and a lattice image for a 10 u.c. spacer is shown in Supplementary Fig. 1. Energy-dispersive X-ray spectroscopy was scanned for elemental identification of Mn, La and Sr.

Magnetic and transport measurement

Magnetization measurements were made using a Quantum Design superconducting quantum interference device magnetometer. Transport measurements were carried out by the Quantum Design Physical Properties Measurement System. The MR anisotropy was measured in a linear geometry as well with two different directions of the applied magnetic field, out of plane and in plane, but perpendicular to the current in each case.

Additional information

How to cite this article: Lü, W. M. et al. Long-range magnetic coupling across a polar insulating layer. Nat. Commun. 7:11015 doi: 10.1038/ncomms11015 (2016).