Observation of spin-orbit magnetoresistance in metallic thin films on magnetic insulators

See allHide authors and affiliations

Science Advances  05 Jan 2018:
Vol. 4, no. 1, eaao3318
DOI: 10.1126/sciadv.aao3318


A magnetoresistance (MR) effect induced by the Rashba spin-orbit interaction was predicted, but not yet observed, in bilayers consisting of normal metal and ferromagnetic insulator. We present an experimental observation of this new type of spin-orbit MR (SOMR) effect in the Cu[Pt]/Y3Fe5O12 (YIG) bilayer structure, where the Cu/YIG interface is decorated with nanosize Pt islands. This new MR is apparently not caused by the bulk spin-orbit interaction because of the negligible spin-orbit interaction in Cu and the discontinuity of the Pt islands. This SOMR disappears when the Pt islands are absent or located away from the Cu/YIG interface; therefore, we can unambiguously ascribe it to the Rashba spin-orbit interaction at the interface enhanced by the Pt decoration. The numerical Boltzmann simulations are consistent with the experimental SOMR results in the angular dependence of magnetic field and the Cu thickness dependence. Our finding demonstrates the realization of the spin manipulation by interface engineering.


Relativistic spin-orbit interaction (SOI) plays a critical role in a variety of interesting phenomena, including the spin Hall effect (SHE) (13), topological insulators (4), and the formation of skyrmions (5, 6). In SHE, a pure spin current transverse to an electric current can be generated in conductors with strong SOI, such as Pt and Ta (7, 8). The inverse SHE (ISHE) is generally used to detect the spin current electrically by converting a pure spin current into a charge current (9, 10). It was recently discovered that the interplay of the SHE and ISHE in a nonmagnetic heavy metal (NM) with strong SOI in contact with a ferromagnetic insulator (FI) leads to an unconventional magnetoresistance (MR)—the spin Hall magnetoresistance (SMR), in which the resistance of the NM layer depends on the direction of the FI magnetization M (1113). SMR has been observed in several NM/FI systems and even in metallic bilayers (1417). However, it has been argued that SMR may originate from the magnetic moment in the NM layer induced by the magnetic proximity effect (MPE) (18). These two mechanisms were proposed to be distinguished by the angular-dependent MR measurements (11, 13). Very recently, another type of MR, the Hanle MR (HMR), was demonstrated in a single metallic film with strong SOI owing to the combined actions of SHE and Hanle effect (19). HMR depends on the direction and the strength of the external magnetic field H, rather than that of M in SMR. Within the framework of SMR, because of the negligible SOI in Cu (20), one would not expect any MR effect in a Cu/FI bilayer.

Recently, Grigoryan et al. (21) predicted a new type of MR effect in the NM/FI systems when a Rashba-type SOI is present at the interface between NM and FI. This new spin-orbit MR (SOMR) works even with light metals such as Cu or Al, with negligible bulk SOI, provided that the Rashba SOI is present at the NM/FI interface. However, because of the identical angular dependence on M direction for SOMR and SMR, it is difficult to distinguish SOMR from SMR in systems like Pt/Y3Fe5O12 (YIG), where both SOMR and SMR are present in principle. Here, we report the first observation of SOMR in a Cu/YIG bilayer, where the Rashba SOI at the Cu/YIG interface is enhanced by an ultrathin Pt layer (<1 nm). We also confirmed that SOMR almost disappears when Pt is placed inside or on the other side of the Cu layer, indicating that SMR from the ultrathin Pt layer cannot be the origin of the observed MR and that the Pt decoration of the Cu/YIG interface is crucial for SOMR. The observed SOMR has the same angular dependence as the SMR in Pt/YIG, in agreement with the SOMR prediction (21). The monotonous Cu thickness dependence of SOMR is different from the nonmonotonous dependence of SMR (13, 18). Both the angular and Cu thickness dependences of the observed MR are in good agreement with our Boltzmann simulations based on the SOMR mechanism. In addition, the MR shows two maxima as the Pt layer thickness increases, in sharp contrast to that of SMR (13, 22).


Sample morphology and structure

The YIG films used in this study are 10 nm thick, unless otherwise stated, and grown by pulsed laser deposition (PLD) on Gd3Ga5O12 (GGG) (111) substrates. The surface morphology of the YIG films was characterized by atomic force microscopy (AFM), as shown in Fig. 1A. The film is fairly smooth with a root mean square (RMS) roughness of 0.127 nm and a peak-to-valley fluctuation of 0.776 nm. The 0.4-nm-thick Pt layer, thinner than the peak-to-valley value of the YIG film, deposited on YIG by magnetron sputtering forms the nanosize islands with an RMS roughness of ~0.733 nm, as shown in Fig. 1B. This discontinuous Pt layer is nonconductive with a resistance over the upper limit of a multimeter. The surface roughness is reduced after the deposition of Cu onto Pt, as shown in fig. S1. Figure 1C presents the cross-sectional high-resolution transmission electron microscopy (HRTEM) image of the Au(3)/Cu(4)[Pt(0.4)]/YIG films, where the numbers are the thicknesses in nanometers. The YIG film is single-crystalline and smooth. The lattice constant of the YIG film is determined to be 1.2234 nm, as compared to 1.2366 nm for the bulk YIG. A clear interface is observed between the metallic films and the YIG film. The metallic films are polycrystalline.

Fig. 1 Sample characterization.

(A) AFM image of YIG(10)/GGG (RMS roughness, 0.127 nm). (B) AFM image of Pt(0.4)/YIG(10)/GGG (RMS roughness, 0.733 nm). (C) HRTEM image of the Au(3)/Cu(4)[Pt(0.4)]/YIG heterostructure, where Au is used to prevent the oxidation.

Field-dependent magnetization and transport measurements

Here, all the measurements were performed at room temperature. The YIG film is almost isotropic in the film plane, with a coercivity of about 0.4 Oe, as shown in Fig. 2A. Because of the large paramagnetic background of the GGG substrate, the magnetization of a thin YIG/GGG film in the out-of-plane geometry is difficult to measure. We measured the 400-nm-thick YIG/GGG(111) film instead. As shown in Fig. 2B, the magnetization is saturated at ~1800 Oe. The saturation magnetization Ms of our YIG film is determined to be 164.5 emu/cm3, as measured by ferromagnetic resonance (FMR) (see the Supplementary Materials). In comparison, Ms of bulk YIG is 140 emu/cm3.

Fig. 2 Field-dependent magnetization and transport measurements.

Magnetic hysteresis loops of (A) YIG(10)/GGG with field in-plane and (B) YIG(400)/GGG with field out-of-plane. ρ measured on the Cu(2)[Pt(0.4)]/YIG(10)/GGG sample for H applied along the (C) x axis, y axis, and (D) z axis, respectively.

Figure 2 (C and D) presents the resistivity ρ as a function of H for the Cu(2)[Pt(0.4)]/YIG(10) sample. In experiments, H was applied (i) along the direction of the current I (x axis), (ii) in the sample plane and perpendicular to the current direction (y axis), and (iii) perpendicular to the sample plane (z axis). The MR effects are present in all measurements. For H along x and y directions, ρ shows two peaks around the coercive fields of YIG. For H along z direction, ρ shows a minimum at H = 0 and remains almost a constant value above the saturation field. These features indicate that the MR effects are intimately correlated with M, meaning that the observed MR effects are not HMR.

Angular-dependent MR measurements

To further study the anisotropy of the MR effects in Cu[Pt]/YIG, we performed the angular-dependent MR measurements. Figure 3A shows Δρ/ρ of the Cu(3)[Pt(0.4)]/YIG(10) sample with rotation of H in the xy (α scan), yz (β scan), and xz (γ scan) planes, where α, β, and γ are the angles between H and x, z, and z directions, respectively, as defined in the inset of Fig. 3A. The applied magnetic field strength (H = 1.5 T) is large enough to align M with H. The MR effect is anisotropic. The MR ratio, defined as Δρ/ρ = [ρ(angle) − ρ(angle = 90°)]/ρ(angle = 90°), in α and β scans is about 0.012%, which is comparable to the SMR ratio in Pt/YIG (see fig. S3) (11, 13, 23).

Fig. 3 Angular-dependent MR measurements.

(A) Angular-dependent MR measurements in the xy, yz, and xz planes for Cu(3)[Pt(0.4)]/YIG. The solid lines are the Boltzmann simulation results. (B) Angular-dependent MR measurements in the xy plane for several control samples.

Next, we investigated the origin of the observed MR effect. Considering that Pt on YIG may suffer from the MPE-induced ferromagnetic moment and the corresponding anisotropic MR (AMR) (24), we replaced Pt by a 0.4-nm-thick Au layer, which is well known to have a negligible MPE (25). The MR effect of 0.002% still appears as shown in Fig. 3B, comparable to the SMR ratio in Au/YIG (see fig. S4), ruling out MPE as the origin of the observed MR. Furthermore, the MR ratios of the Cu(3)[Pt(0.4)]/YIG(10) sample in α and β scans are comparable and almost one order of magnitude larger than the MR ratio in the γ scan. This is different from AMR of a ferromagnetic metal, where the MR ratio in α and γ scans is much larger than that in the β scan (11, 14, 24). Therefore, the MPE-induced AMR can be ruled out.

The behaviors of the MR angular dependence follow the SMR scenario well (11, 1315, 17). However, with several control experiments, we can unambiguously exclude SMR as the explanation for our observations.

First, the observed MR amplitude cannot be explained by SMR. In our samples, the 0.4-nm-thick ultrathin Pt layer is nonconductive, and the conductivity of bulk Pt is about one order of magnitude smaller than that of bulk Cu, meaning that the current mainly passes through the Cu layer. We prepared a 3-nm-thick single-layer Cu on YIG without interface decoration and performed the MR angular-dependent measurement in the α scan. MR is not observed, as shown in Fig. 3B, evidencing that the Pt-decorated interface is indispensable. A conductive 0.4-nm-thick Pt layer is not available experimentally. Considering that a small fraction of current may flow in the Pt islands, there is a possibility of the occurrence of SMR from the Pt islands. According to the reported SMR results in Pt/YIG bilayers, the SMR ratio in Pt/YIG decreases rapidly with decreasing Pt thickness when the Pt thickness is less than about 3 nm (13, 18). The SMR ratio of Pt(0.4)/YIG is extrapolated to be well below 0.01% from the previously reported Δρ/ρ versus Pt thickness data (13, 18). Considering the pronounced shunting current of the highly conductive Cu layer, the SMR ratio should be significantly reduced in Cu[Pt]/YIG, that is, much less than 0.01%. In comparison, the MR ratio is as large as ~0.012% in Cu(3)[Pt(0.4)]/YIG (see Fig. 3A). Therefore, the SMR mechanism cannot explain our observations.

Second, the potential enhancement of SMR caused by intermixing or alloying between a strong SOI material and a weak SOI material can be excluded (17, 26, 27). For this purpose, we prepared two types of control samples, with the 0.4-nm-thick Pt layer either on top of or inserted inside the Cu layer: [Pt(0.4)]Cu(3)/YIG and Cu(1)[Pt(0.4)]Cu(3)/YIG. Because both samples are fabricated under the same condition as the Cu[Pt]/YIG samples, the intermixing of Pt and Cu should be similar. The MR vanishes in the [Pt(0.4)]Cu(3)/YIG and Cu(1)[Pt(0.4)]Cu(3)/YIG samples, as shown in Fig. 3B. These results rule out the Pt-Cu alloying–induced SMR. Thus, we conclude that the observed MR effect is not SMR.

Cu thickness–dependent transport measurements

To identify the physical origin of the observed unusual MR, we carried out the Cu thickness–dependent measurements. Figure 4A presents the angular-dependent MR measurements of Cu(tCu)[Pt(0.4)]/YIG in α scans for various Cu thicknesses (tCu). Obviously, the MR ratio steadily decreases with increasing tCu, highlighting the importance of the Pt-decorated Cu/YIG interface. The monotonous NM thickness dependence of this MR is in sharp contrast to the nonmonotonous behavior of SMR, which peaks at ~3 nm for Pt/YIG (13, 18). The Cu thickness dependence of ρ and the MR ratio extracted from Fig. 4A are shown in Fig. 4B. For very thin Cu film (tCu ≤ 5 nm), ρ dramatically increases with decreasing tCu, indicating that ρ is dominated by the interface/surface scatterings.

Fig. 4 Cu thickness– and Pt thickness–dependent transport measurements.

(A) Angular-dependent MR measurements in the xy plane for Cu(tCu)[Pt(0.4)]/YIG and Cu(3)/Pt(tPt)/YIG samples. (B) Cu thickness dependence of the MR ratio and ρ for Cu(tCu)[Pt(0.4)]/YIG. The solid lines are the Boltzmann simulation results. (C) Pt layer thickness dependence of the MR ratio and ρ for Cu(3)/Pt(tPt)/YIG. The solid lines are guide to the eyes.

Besides SMR, there is another type of MR recently predicted, which has the same angular dependence as we found (see Fig. 3A) (21). It originates from the Rashba SOI at the interface of an NM/FI bilayer. By comparing the samples of Cu[Pt]/YIG, [Pt]Cu/YIG, and Cu[Pt]Cu/YIG, one can see that only the Cu[Pt]/YIG samples exhibit a significant MR (see Fig. 3B). It strongly suggests that the MR observed in our experiments is the SOMR predicted by Grigoryan et al. (21) and that the Pt decoration enhances the Rashba SOI at the Cu/YIG interface.

First-principles calculations and Boltzmann simulations

To prove that the Pt decoration can induce Rashba SOI at the Cu/YIG interface, we carried out first-principles band structure calculations based on (i) a Cu ultrathin film of 14 monolayers, (ii) the same Cu film as (i) but covered by Au on surfaces on both sides, (iii) the same Cu film as (i) but covered by Pt on both surfaces, and (iv) the Pt layer inside the Cu film. By comparing these four different scenarios, we can see that there is no clear Rashba effect in the bare Cu film and the one covered by Au. A strong Rashba effect appears only for Pt on the Cu film surface (the details of calculations are given in the Supplementary Materials).

For a quantitative analysis, we use a Boltzmann formalism to calculate the charge and spin transport in an NM/FI bilayer structure. We solve the following spin-dependent Boltzmann equation in the NM layerEmbedded Image(1)where fα = 0,x,y,z(r,k) is the four-component distribution function denoting the charge/spin occupation at position r and wave vector k. The interface at FI z = z+ contains a Rashba-type SOI described by the Hamiltonian: Embedded Image, where η is the strength of the Rashba SOI, Embedded Image is the normal direction of the interface, and Embedded Image is the momentum operator. HR gives rise to an anomalous velocity localized at the interface. The Boltzmann equation is solved by discretizing the spherical Fermi surface of Cu and the real space in the z direction of Cu film. Using the full distribution function, we calculate all charge/spin transport properties, including the longitudinal and transverse conductivities. This method extends the earlier Boltzmann method developed for current perpendicular-to-plane structure–like spin valves to current in-plane structure–like NM/FI bilayers (2831) by taking into account the surface roughness and Rashba SOI at the interface. The details of the simulations are given in the Supplementary Materials.

In the numerical Boltzmann calculation, there are only two fitting parameters: the surface roughness and the Rashba coupling constant. All other parameters either are given by the experiment (such as the film thickness) or can be determined otherwise (such as the bulk relaxation time). Using the quantum description of the rough surface (3234), we are able to fit the thickness dependence of ρ in the ultrathin Cu film to a reasonably good precision, as shown in Fig. 4B. It is quite surprising considering that there is only one fitting parameter—the surface roughness. Once the surface roughness is determined, we can calculate the magnetization angular dependence of ρ, in good agreement with the experimental results (see Fig. 3A), from which we obtain the SOMR ratio. The calculated SOMR ratio is shown in Fig. 4B, which demonstrates monotonic decreasing behavior as a function of Cu film thickness, consistent with our experimental results but very different from the nonmonotonic behavior observed in SMR (13, 18).

Pt thickness–dependent MR measurements

Finally, to further differentiate SOMR from SMR, we carried out the Pt thickness (tPt)–dependent measurements. To reduce the sample fluctuation, we fabricated the YIG films successively under the same condition. Figure 4A shows the angular-dependent MR measurements of Cu(3)/Pt(tPt)/YIG in α scans, with H = 2000 Oe. The MR ratio extracted from Fig. 4A exhibits nonmonotonous behavior with increasing tPt, as shown in Fig. 4C. Two separate regimes can be identified: (i) the SOMR regime for tPt < 1 nm and (ii) the conventional SMR regime for tPt > 2.2 nm (see Fig. 4C). For tPt < ~0.6 nm, ρ and Δρ/ρ increase with increasing tPt because the Pt islands not only introduce the interface scattering but also enhance the Rashba SOI. For ~0.6 nm < tPt < ~1 nm, the Pt islands start to form a complete layer, leading to the reduction of the interface roughness and the rapid decrease of ρ as seen in Fig. 4C. The MR ratio continues to increase in this region because of the enhanced Rashba SOI with increasing Pt coverage on YIG. For ~1 nm < tPt < ~2 nm, ρ is smaller than the resistivity of Cu/YIG, suggesting that the interface scattering has minor contribution to ρ. Because SOMR is caused by the interface scattering, the MR ratio rapidly drops in this region. A sizable SMR ratio only appears when tPt > 2 nm in Pt/YIG (13, 18, 22). Therefore, around tPt ~ 2 nm, both SOMR and SMR are small, resulting in a minimum in MR. In the SMR regime, the SMR ratio exhibits a maximum, as expected for SMR (13, 18, 22). This result demonstrates the differences between SOMR and SMR.

A theoretical calculation shows that a rough interface can enhance SHE (34). To understand the role of the roughness to the SOMR, we fabricated the control sample of Cu(3)[Ag(0.7)]/YIG. The RMS roughness of Ag(0.7)/YIG is 0.797 nm, as shown in fig. S9A, similar to that of Pt(0.4 nm)/YIG. Owing to the weak SOI in Ag, the Rashba SOI in Cu[Ag]/YIG is expected to be weak. We do not observe any MR effect down to 5 × 10−6 in Cu(3)[Ag(0.7)]/YIG, as shown in fig. S9B. This result means that the rough surface alone cannot cause the SOMR.


In conclusion, we report the first observation of the SOMR effect predicted recently (21) at room temperature in Cu/YIG films with the Pt decoration at the interface. We show that this MR effect is caused by the enhanced Rashba SOI at the Pt-decorated interface. The angular dependence of SOMR is similar to that of SMR, but all other features are different, such as the increasing MR with decreasing Cu thickness. The amplitude of the SOMR ratio is comparable to that of the SMR ratio in Pt/YIG, highlighting the importance of the NM/FI interfaces. Our finding demonstrates the possibility of realizing spin manipulation by interface decoration.


The single-crystalline YIG films were epitaxially grown on GGG (111) substrates by the PLD technique using a KrF excimer laser with a wavelength of 248 nm. The PLD system was operated at a laser repetition rate of 4 Hz and an energy density of 10 J/cm2. The distance between the substrate and the target is 50 mm. Before film deposition, the chamber was evacuated to a base pressure of 1 × 10−7 torr. The YIG films were deposited at ~730°C in an oxygen pressure of 0.05 torr. The growth of the YIG films was monitored by in situ reflection high-energy electron diffraction. The structure was further examined by x-ray diffraction and HRTEM. The magnetic properties of all YIG films were characterized using a vibration sample magnetometer. Then, we used magnetron sputtering to deposit polycrystalline metallic films onto the YIG films via dc sputtering at room temperature with a shadow mask to define 0.3-mm-wide and 3-mm-long Hall bars. The deposition rate was calibrated by x-ray reflectivity. After the metallic film deposition, the samples were immediately mounted and transferred into a vacuum chamber for the transport measurements to minimize the metal oxidation. The resistance was measured by a Keithley 2002 multimeter in a four-probe mode. For the angular-dependent MR measurements with a magnetic field of less than 5000 Oe, the resistance was monitored as the magnet was rotated. The angular-dependent MR measurements with a magnetic field larger than 5000 Oe were performed in a physical property measurement system equipped with a rotatory sample holder.


Supplementary material for this article is available at

section S1. AFM images of Cu(tCu)[Pt (0.4)]/YIG(10)/GGG(111)

section S2. Magnetic properties of the YIG films

section S3. SMR in Pt/YIG

section S4. SMR and AFM image of Au/YIG

section S5. First-principles calculations

section S6. Boltzmann simulations

fig. S1. AFM images of Cu(tCu)[Pt(0.4)]/YIG(10)/GGG (111).

fig. S2. FMR of the YIG films.

fig. S3. SMR in Pt/YIG.

fig. S4. SMR and AFM image of Au/YIG.

fig. S5. The band structures of Cu, Au/Cu/Au, and Pt/Cu/Pt.

fig. S6. The spin textures of outer and inner bands.

fig. S7. The Rashba splitting in Cu/Pt/Cu.

fig. S8. Specular and diffusive interface scattering in the NM/FI bilayer.

fig. S9. AFM image of Ag(0.7)/YIG and MR of Cu(3)[Ag(0.7)]/YIG.

References (3551)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: Funding: L.Z. and D.W. were supported by the National Key R&D Program of China (2017YFA0303202), the National Natural Science Foundation of China (NSFC) (11674159, 51471086, and 11727808), and the National Basic Research Program of China (2013CB922103). H.S. and J.X. acknowledge the support of NSFC (11474065 and 11722430) and the National Key R&D Program of China (2016YFA0300702). K.X. was supported by the National Key R&D Program of China (2017YFA0303300) and NSFC (11734004, 61774017, and 21421003). Author contributions: J.X. and D.W. designed and supervised the project. L.Z. and Z.L. prepared the samples. L.Z. performed the transport measurements with support from Z.L., P.W., and S.J. H.S. performed the Boltzmann simulations under the supervision of J.X. and K.X. K.L. performed the first-principles calculations under the supervision of H.X. L.S. and Y.C. were responsible for the HRTEM characterization. J.X., D.W., and L.Z. wrote the manuscript, and J.D., H.X., H.D., and K.X. commented on the manuscript. All authors discussed the results and reviewed the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.
View Abstract

Navigate This Article