We study the transport of thermally excited non-equilibrium magnons through the ferrimagnetic insulator YIG using two electrically isolated Pt strips as injector and detector. The diffusing magnons induce a non-local inverse spin Hall voltage in the detector corresponding to the so-called non-local spin Seebeck effect (SSE). We measure the non-local SSE as a function of temperature and strip separation. In experiments at room temperature we observe a sign change of the non-local SSE voltage at a characteristic strip separation d0, in agreement with previous investigations. At lower temperatures however, we find a strong temperature dependence of d0. This suggests that both the angular momentum transfer across the YIG/Pt interface as well as the transport mechanism of the magnons in YIG as a function of temperature must be taken into account to describe the non-local SSE.
Magnons, the collective excitations in magnetically ordered systems, represent an attractive option for information transfer and processing. Using the ferrimagnetic insulator Yttrium Iron Garnet (YIG) as a model system, magnon-based information processing schemes have been put forward on the basis of coherently excited spin waves.1–4 Recent experiments in YIG/Pt heterostructures furthermore show that information can also be carried by incoherent non-equilibrium magnons5,6 diffusing in YIG. This approach even allows for the implementation of logic operations within the magnetic system.7 The non-equilibrium magnons can be excited and detected electrically via spin scattering mechanisms at the YIG/Pt interface. In a scheme referred to as magnon-mediated magnetoresistance (MMR), magnons are generated by driving a dc charge current through an injector Pt strip and detected as a non-local voltage in a second strip. The MMR effect has been studied as a function of the distance d between injector and detector,5 temperature6 and magnetic field magnitude and orientation,8 allowing for the extraction of the length scales involved in the magnon diffusion process. In addition to the electrical injection, non-equilibrium magnons can also be generated thermally, via local Joule heating in the injector strip. The ensuing thermal non-local voltage is called non-local spin Seebeck effect in analogy to the well established local (longitudinal) spin Seebeck effect (SSE).9 While the microscopic mechanisms and in particular the relevant length scales for the SSE have been investigated in quite some detail,10,11 the physics behind the non-local thermal signal is not well established. In this context, the non-local SSE has recently been studied at room temperature as a function of strip separation d and YIG thickness.12 While at short distances d the local and the non-local SSE signals have the same sign, for larger distances the non-local SSE amplitude is inverted, which was attributed to the profile of the non-equilibrium magnon distribution in the YIG film.
In this article, we systematically study the non-local spin Seebeck effect in YIG/Pt nanostructures as a function of temperature and strip separation and find that the non-local SSE voltage changes sign at a characteristic strip separation d0, which is strongly temperature dependent. We interpret our findings as evidence of a complex interplay between the temperature dependencies of the interfacial transparency, i.e. angular momentum transfer across the YIG/Pt interface, and the diffusive properties of the thermally excited non-equilibrium magnons.
We investigate the non-local spin Seebeck effect in YIG/Pt bilayers fabricated and nano-patterned at the Walther-Meißner-Institut (sample series A) and at Johannes Gutenberg-University Mainz (sample series B). Series A was fabricated starting from a commercially available 2 μm thick YIG film grown onto (111) oriented Gd3Ga5O12 (GGG) via liquid phase epitaxy (LPE). After Piranha cleaning and annealing (see Ref. 13 for details) to improve the interface quality, 10 nm of Pt were deposited onto the YIG film using electron beam evaporation. A series of nanostructures consisting of 2 parallel Pt strips with length , width and an edge-to-edge separation of were patterned using electron beam lithography followed by Ar ion etching. Series B was fabricated using a 3.3 μm thick LPE-YIG film grown onto a GGG substrate as well. A series of strips with width and edge-to-edge distance were patterned using electron beam lithography followed by a lift-off process with a Pt thickness of 7.5 nm deposited by dc sputtering. An optical micrograph of one of these nanostructures is depicted in Fig. 1 (a). Figure 1 (b) is a schematic representation of the YIG/Pt bilayer (not to scale). The color coding in the YIG film qualitatively represents the magnon accumulation profile as proposed in Ref. 12, where red/blue corresponds to a magnon depletion/accumulation with respect to the local thermal equilibrium value. Two detectors are exemplarily shown for different injector/detector separations d, i.e. above an area of magnon depletion or accumulation for detector 1 and 2, respectively. The implications of the magnon accumulation profile for the detected non-local thermal signal are discussed in the following.
(a) Optical micrograph of a non-local nanostructure with two Pt strips (bright) on a YIG film (dark), including the electrical wiring. (b) Sketch of the YIG/Pt heterostructure: a dc charge current Jc (not shown here) is applied to the injector strip (left) and the spin Seebeck signal is detected locally via the ISHE. The corresponding non-local thermal signal is measured along the detector strips for different strip separations. The color coding gives a qualitative profile of the magnon accumulation in the YIG film, as proposed in Ref. 12, where red corresponds to (magnon depletion) and blue to (magnon accumulation). In the short distance regime at the injector and detector, such that the same sign is expected for local and non-local SSE. With increasing distance from the injector, and consequently the spin current across the interface as well as the detected ISHE voltage change sign.
(a) Optical micrograph of a non-local nanostructure with two Pt strips (bright) on a YIG film (dark), including the electrical wiring. (b) Sketch of the YIG/Pt heterostructure: a dc charge current Jc (not shown here) is applied to the injector strip (left) and the spin Seebeck signal is detected locally via the ISHE. The corresponding non-local thermal signal is measured along the detector strips for different strip separations. The color coding gives a qualitative profile of the magnon accumulation in the YIG film, as proposed in Ref. 12, where red corresponds to (magnon depletion) and blue to (magnon accumulation). In the short distance regime at the injector and detector, such that the same sign is expected for local and non-local SSE. With increasing distance from the injector, and consequently the spin current across the interface as well as the detected ISHE voltage change sign.
In order to study the devices in series A and B, the samples were mounted in the variable temperature insert of a superconducting magnet cryostat (). For series A an external magnetic field was rotated in the thin film plane, while for series B the external magnetic field was applied along the y-direction and swept from to (, in Fig. 1 (b)). For local longitudinal spin Seebeck effect measurements in one single (injector) strip, we used the current heating method described in Ref. 14: a charge current is applied to the Pt strip along the x direction using a Keithley 2400 source meter, inducing Joule heating in the normal metal. The ensuing temperature gradient across the Pt/YIG interface gives rise to the spin Seebeck effect and generates a spin current Js flowing across the interface, with the spin current spin polarization s determined by the orientation of the magnetization M in YIG. This spin current is accompanied by a charge current JISHE flowing along the x-direction in the Pt, as shown in Fig. 1 (b). The voltage drop Vloc, which includes the local SSE and the resistive response of the injector strip is recorded in a 4-point measurement configuration using a Keithley 2182 nanovoltmeter. Since the spin Seebeck effect is of thermal origin, the SSE voltage is proportional to the Joule heating power in the Pt and therefore independent of the heating current direction. Using the switching scheme of Ref. 14, we extract Vtherm,loc = (Vloc(+Jc) + Vloc(−Jc))/2, and thereby eliminate additional resistive effects such as the spin Hall magnetoresistance. Note that we do not consider contributions from the unidirectional SMR,15 since to our knowledge there is no experimental evidence of this effect in ferrimagnetic insulator/normal metal bilayers to date.
Figure 2 (a) shows Vtherm,loc measured as a function of the magnetic field orientation with respect to the x-axis at 50 K for a device from series A. We observe the characteristic SSE dependence yielding a positive amplitude of the local SSE, as expected in YIG/Pt heterostructures for this field configuration.16
(a) Local spin Seebeck voltage detected at the injector strip at 50 K as a function of the in-plane magnetic field orientation with respect to the x axis. (b), (c) Non-local thermal voltage detected at the second strip at for a strip separation of and , respectively. (d) Non-local SSE amplitude ASSE,nl extracted from in-plane field rotations at temperatures between and as a function of the injector-detector separation d.
(a) Local spin Seebeck voltage detected at the injector strip at 50 K as a function of the in-plane magnetic field orientation with respect to the x axis. (b), (c) Non-local thermal voltage detected at the second strip at for a strip separation of and , respectively. (d) Non-local SSE amplitude ASSE,nl extracted from in-plane field rotations at temperatures between and as a function of the injector-detector separation d.
Using an additional nanovoltmeter, we simultaneously measure the voltage drop Vnl arising along the unbiased and electrically isolated second Pt strip. In analogy to the local thermal signal, the non-local thermal voltage is extracted as Vtherm,nl = (Vnl(+Jc) + Vnl(−Jc))/2 in order to distinguish it from resistive non-local effects such as the magnon mediated magnetoresistance.6 Vtherm,nl as a function of the external magnetic field orientation at 50 K is depicted in Fig. 2 (b) and (c) for strip separations of and (series A). In both devices we observe a sin dependence, with an amplitude ASSE,nl about one order of magnitude smaller than for the local SSE. In the following we do not discuss the absolute amplitude of the non-local thermal signal in detail, but focus on the sign of ASSE,nl as a function of separation between the Pt strips. While the signal is positive for the device, a negative ASSE,nl is observed in the device with a larger injector-detector separation of . In order to confirm this sign change, we extract the amplitude of the non-local SSE measured at 50 K in different devices with strip separations ranging from to . The resulting data is shown in Fig. 2 (d) as green symbols. Indeed, a sign change is observed at a strip separation . Repeating these measurements as a function of temperature in the range between 10 K and 300 K yields the data compiled in Fig. 2 (d). For all temperatures, a sign change in ASSE,nl is observed as a function of the strip separation. Invariably, for small gaps the local and non-local SSE are both positive, but for large gaps the non-local SSE becomes negative. The experimental data in Fig. 2 (d) show that the critical strip separation d0, which is defined by ASSE = 0, shifts to larger values as the temperature increases. The values d0 extracted from Fig. 2 for different temperatures are shown in Fig. 3 as red symbols for sample series A. With increasing temperature d0 increases monotonically and seems to saturate around .
Temperature dependence of the critical strip separation d0 at which the non-local SSE changes sign for sample series A (red) and B (blue).
Temperature dependence of the critical strip separation d0 at which the non-local SSE changes sign for sample series A (red) and B (blue).
Similar experiments as a function of temperature and strip separation were performed on devices from series B and the critical strip separation extracted from these measurements is included in Fig. 3 as blue squares. While the temperature dependence is much steeper, a qualitatively similar increase of d0 with temperature is observed in both series.
This characteristic sign change in the non-local SSE in YIG/Pt heterostructures above a particular separation d0 has been previously observed by Shan et al.12 at room temperature and was attributed to the spatial profile of the non-equilibrium magnon accumulation in the YIG film, as shown schematically in Fig. 1 (b). Magnons are thermally excited in the ferrimagnet due to Joule heating in the injector strip and diffuse vertically towards the GGG/YIG interface as well as laterally to the sample edges. According to the model proposed by Shan et al.,12 this leads to a depletion of magnons (, red in Fig. 1 (b)) compared to the thermal equilibrium population beneath the injector. On the other hand, diffusing magnons accumulate further away from the injector, giving rise to (blue in Fig. 1 (b)). This model is applied to describe the increase of d0 with increasing YIG thickness observed by Shan et al.:12 in contrast to phonons, the magnons cannot cross the YIG/GGG interface and accumulate there. As a consequence, d0 (which marks the sign change of ) shifts to smaller values for thinner YIG films. Since the overall profile of is governed by diffusive magnon transport, the corresponding length scales can reach several μm.12 Note that while the magnon accumulation profile is affected by changes in the magnon diffusion length , d0 is not directly proportional to but exhibits a more complex behavior as discussed in the following. Furthermore, we expect d0 to be independent of the Joule heating power to a first approximation, since in linear response the latter changes only the absolute number of non-equilibrium magnons which are thermally generated at the injector. The boundary conditions at the Pt/YIG and YIG/GGG interfaces, however, determine the qualitative profile of the magnon accumulation and remain independent of the heating power.
As shown in Fig. 1 (b), the sign of determines the direction of the interfacial spin current Js at the detector, i.e. towards (away from) the YIG for negative (positive) at detector 1 (detector 2), and consequently governs the sign of the measured non-local ISHE voltage. Non-local SSE measurements as a function of the strip separation therefore allow us to map out the non-equilibrium magnon distribution in the YIG film. In particular, the characteristic length d0 for the sign change of can be determined.
In order to rationalize the measured temperature dependence of d0, the parameters governing the angular momentum transfer across the YIG/Pt interface as well as the magnon diffusion process need to be analyzed. It has been shown that the transparency of the YIG/Pt interface, described by the effective spin-mixing conductance gs, influences the magnon accumulation and hence the sign-reversal distance d0.12 For a fully opaque interface (obtained using an Al2O3 interlayer) which suppresses angular momentum backflow into the injector and therefore preserves a strong magnon depletion, an increase of the sign-reversal distance d0 was observed.12 Previous measurements of the MMR effect in a YIG/Pt heterostructure - prepared using exactly the same fabrication method as for sample series A - as a function of temperature have shown that the MMR signal decreases with decreasing temperature,6 consistent with as predicted by theory.17,18 A similar behavior was observed in a Pt/YIG/Pt heterostructure which is based on dc sputtered Pt patterned by a lift-off process.19 We therefore assume a qualitatively similar temperature dependence of the interfacial transparency for evaporated (series A) and sputtered Pt (series B). However, with this temperature dependence we expect a decreasing transparency of the YIG/Pt interface with decreasing temperature, leading to an increase of d0 at low temperatures according to the model presented in Ref. 12. Since this is not consistent with our experimental observations depicted in Fig. 3, the interface properties alone are not sufficient to describe the temperature dependence of the non-local SSE.
In addition to the interfacial transparency, the magnon diffusion length and the magnon spin conductivity determine the spatial distribution of the non-equilibrium magnons in YIG. We extracted the magnon diffusion length from temperature dependent MMR measurements conducted in the sample series A and found an increase of with decreasing temperature by about a factor of 3 between 300 K and 50 K, following a 1/T dependence. This is different from the temperature independent diffusion length reported by Cornelissen et al.,20 who extracted together with a magnon spin conductivity vanishing at low temperatures. While the detailed evolution of and with T thus must be studied more systematically in future work, it is clear that these quantities depend on temperature. This implies that they can qualitatively impact the magnon diffusion process and consequently the non-local SSE. Indeed, the strong dependence of d0 on the relative amplitudes of gs, and at a fixed temperature has been demonstrated by Shan et al. using a one-dimensional analytical model for the spin Seebeck effect.12
Based on the available experimental data, we can conclude that for a quantitative modeling of the non-local SSE the temperature dependence of both the angular momentum transfer across the YIG/Pt interface and the magnon diffusion in YIG must be taken into account. We furthermore note that additional effects due to phononic heat transport, known to be of importance for the local SSE, cannot fully be excluded at this point as a source of the non-local SSE-like signal. The local SSE originates from a finite difference of the effective temperature of the magnon and phonon subsystems in YIG close to the YIG/Pt interface, giving rise to a magnon spin current.21 While this difference may also influence the thermal signal measured at the non-local detector strip, it was shown that the thermalization of the magnon and phonon subsystems takes place on a length scale of the order of a few nm at room temperature,22–24 due to a very efficient (magnon conserving) magnon-phonon scattering. The effective temperature model described in Ref. 21 is therefore applicable mainly in the local limit, i.e. close to the injector. The long-distance non-local limit, however, will be dominated by diffusing magnons and the magnon-phonon scattering which does not conserve the number of magnons and can reach a larger length scale of the order of several m,18 as discussed above. Further contributions from phonons, i.e. an enhancement of the SSE signal due to phonon drag25,26 might also play a role for the magnon accumulation profile (see supplementary material).
In summary, we have measured the non-local SSE in a YIG/Pt heterostructure as a function of injector-detector distance and temperature. The non-local SSE changes sign at a characteristic injector-detector separation d0, consistent with previous observations put forward by Shan et al.12 We furthermore observe a decrease of the characteristic separation d0 with decreasing temperature, which was confirmed by the independent study of two sample series in different setups. Our results suggest a complex dependence of the non-local SSE on interfacial transparency, magnon diffusion properties as well as phonon heat transport.
Note added. During the review process of this paper, we have become aware that Zhou et al.27 reported similar results on the temperature dependence of the non-local SSE sign change.
SUPPLEMENTARY MATERIAL
See supplementary material for details on the temperature dependence of the local and non-local SSE amplitude.
ACKNOWLEDGMENTS
KG acknowledges N. Vlietstra for discussions and JC thanks S. Kauschke for sample preparation. This work is financially supported by the Deutsche Forschungsgemeinschaft through the Priority Program Spin Caloric Transport (GO 944/4, GR 1132/18, KL 1811/7) and the SFB TRR 173 Spin+X, the Graduate School of Excellence Materials Science in Mainz (MAINZ) and EU projects (IFOX, NMP3-LA-2012 246102, INSPIN FP7-ICT-2013-X 612759).