1 Introduction
Thermal spin transfer, which is the spin transfer created by thermoelectric effects, has attracted considerable attention in the field of “spin caloritronics” [
1]. The large thermally driven spin transfer torque (STT) in all-metallic spin valves was first predicted by Hatami
et al. [
2], and it was experimentally found in a Co|Cu|Co nanowire by Yu
et al. [
3]. Subsequently, thermal torque has been predicted and found in spin valves [
2-
5], domain walls [
6-
9], magnetic tunnel junctions (MTJs) [
10-
13], and even normal-metal|ferromagnet (NM|FM) interfaces [
14-
18]. Among them, the NM|FM interface, especially NM|ferromagnetic insulator (NM|FI) interface, has the potential to serve as the basic unit in spin caloritronics.
The thermal torque at NM|FI interfaces is considered to be interesting for its pure spin-driven process. Two mechanisms are proposed for the interface [
18]: the first one is the local effect of torque on the magnetization via magnon-magnon scattering driven by a spin Seebeck effect (SSE) [
5], the other one is the nonlocal effect of torque on the magnetization via absorption of spin current generated by misaligned magnetization. These mechanisms are based on SSE, which is driven by the temperature difference between the magnon of FM(FI) and the lattice of NM. This difference is large in spin valves, in which the thermal torque is carried by the spin polarizing electrons and driven by the temperature difference between the electrons in FMs. To classify the two types of thermal torques, the torque from the SSE is named as spin Seebeck torque (
) and that from the nonlocal electrons as thermal spin pumping torque (
). The latter can be calculated using the spin pumping theory [
19], which is based on the scattering theory. Recent experiments revealed that a small temperature bias of
can build a spin torque of the order 1.0 μJ/(K⋅m
2) [
20,
21], which is several times larger than that [
10,
22] in ultrathin MgO-based magnetic tunnel junctions (MTJs).
When a ferromagnet is attached to the NM|FI interface, forming a FM|NM|FI multilayer, the spin current pumped out of the FI/FM generates a torque on the FM/FI. Both SEE and anomalous Nernst effect (ANE) are present under the temperature bias. A recent experiment in Co|Cu|YIG demonstrates that SSE dominates over ANE [
23]. As the spin torque overcomes the damping torque, a stable magnetic oscillation or magnetic switching would be generated. For an NM|FI interface, the spin mixing conductance is a good parameter to estimate the power required to generate a spin torque [
20]. Through experiments, it is easier to measure the spin torque than the spin mixing conductance; therefore, it is convenient to derive the spin mixing conductance from the measured spin torque.
Here, we study the thermal torque in Ag|Fe|Ag|YIG multilayers using first principles. Owing to the ultralow intrinsic magnetic damping
, Yttrium iron garnet (Y
3Fe
5O
12) is a promising material for spin-wave applications [
5,
24-
32]. A low magnetic damping [
33-
37]
was reported in YIG thin films, which is close to that in bulk YIG. In our calculations, an enhanced magnetic damping
from interface scattering and
TSSE∼1.0 μJ/(K·m
2) at 300 K were estimated. Based on the Landau–Lifshitz–Gilbert (LLG) equation [
38], we predict that a temperature bias of
across the sandwiched NM is sufficient to reverse the magnetic configurations circularly in a multilayer at room temperature.
This article is organized as follows: In Section 2, we introduce the spin Seebeck torque induced by magnetic precession. In Section 3, we present our results on Ag|Fe|Ag|YIG multilayers, where two methods are used to estimate the spin Seebeck torque from the spin mixing conductance and the magnetic damping parameter. The study is summarized in Section 4.
2 Enhanced magnetic damping and spin Seebeck torque
In this section, we explain the use of scattering theory to calculate the magnetic damping and spin Seebeck torque in NM|FM|NM|FI|NM multilayers, as shown in Fig. 1. The same NM is used for both the right and left leads, and when the FI is sufficiently thick, the right lead shows less effect on the spin transport. Magnetization dynamics in the FI is usually described by the phenomenological LLG equation [
40]:
where
,
is the unit vector of the magnetization texture at position
and time
t,
and
are the effective magnetic field and random magnetic field [
5] due to the magnetization, respectively, and
is the gyromagnetic ratio. The magnetic damping
is a 3×3 tensor, which depends on the orientation of magnetization.
The time derivative of the magnetic energy
E, obtained from the above equation is [
40]
where
Ms is the saturation magnetization. Therefore, the magnetic damping is directly related to the energy pumping process. The rate of energy loss in the scattering region in terms of the scattering matrix [
41,
42] at finite temperatures is given by
where is the scattering matrix with transmission/reflection matrices for the Bloch waves incident from the left (without prime) and right (with prime) leads. Comparing Eqs. (2) and (3), we obtain the magnetic damping tensor expressed by the scattering matrix:
where
V is the volume of magnetization with section area
A and thickness
df. For the FI|NM interface, neglecting the spin-orbit coupling, the magnetic damping tensor can be reduced to a magnetic damping parameter
, where
gr is the real part of the spin-mixing conductance [
41,
42].
At finite temperatures,
α includes the value of temperature via two parameters: the Fermi–Dirac distribution
, and the temperature dependent magnetic coherence volume [
20]
, where
D is the the spin-wave stiffness constant and
ξ is the Riemann Zeta function. The thermal magnetic damping parameter can be expressed as
At zero temperature, Va is reduced to V, and is reduced to at the Fermi energy .
When there is a temperature difference between two magnetic layers, as in the FM|NM|FI multilayer shown in Fig. 1, a thermal spin current arising from the thermally induced spin fluctuation is generated and it flows from the high-temperature layer
across the NM [
5]
The spin current is polarized along the equilibrium magnetization . The thermally induced spin current enters the magnetic layer and generates a spin-transfer torque. When the working frequencies of FM and FI are different, there is a weak interaction between the spin currents being pumped from FI and FM. Therefore, the total spin current in the sandwiched NM would be the sum of the spin currents from both FI and FM. Considering that the spin current pumping out of the NM|FI interface is fully absorbed by the FM, the spin torque exerted on is
which is referred to as the spin Seebeck torque. When the quantum effect in the sandwiched NM is negligible, the total spin Seebeck torque is reduced to [
5]
.
3 Spin Seebeck torque in Fe|Ag|YIG multilayers
Here, we study the enhanced magnetic damping and thermal torque in Ag|Fe|Ag|YIG multilayers, as shown in Fig. 1. We use a
lateral supercell for the fcc Ag, so that it matches the cubic YIG unit cell. The electronic structure of the YIG is calculated using a tight-binding linear-muffin-tin-orbital (TB-LMTO) code in the augmented spherical wave approximation, as implemented in the Stuttgart code [
43-
45] using the generalized gradient approximation (GGA) to the local density approximation (LDA); and the bulky potentials are added into a wave-function matching (WFM) transport package [
46]. During transport calculations, we consider the particle current along the (010) material growth direction, and a
k-mesh in the two-dimensional Brillouin zone (2d-BZ) is used to ensure good transport convergence. Additional numerical details about the electronic structure and transport calculations can be found in the reference [
47].
First, let us consider the
in the Ag|Fe|Ag|YIG multilayers.
can be estimated from the scattering matrix [
48], where
summarize the modes of the left/right (L/R) lead and
is the voltage bias estimated from the local chemical potentials
. In the case of a thicker YIG, the right-going electrons would be reflected to left completely, and no pure electrons will penetrate through the YIG. Here, a
YIG was used for the calculations, and the transmission was on the order of
. Consequently, the in-plane STT from Fe/YIG to YIG/Fe would be zero for the nonlocal nature. However, the out-of-plane STT is non-zero, as shown in the inset of Fig. 2(a), which is due to the local exchange interaction between Fe and YIG. From the energy dependent STT, we obtain the
[
10]. A
of 71 nJ/(K·m
2) was found in Ag|Fe(16)|Ag(16)|YIG at 300 K, where the numbers in brackets indicate the thickness of the atomic layers. However, when Pt is considered as a separate metal, a large
of
nJ/(K·m
2) is observed, as listed in Table 1, which is even larger than that in Fe|MgO(3)|Fe MTJs [
10]. The large difference in
results from the large symmetric difference in energy-dependent STT as shown in the inset of Fig. 2(a), which is related to the large difference in the symmetry of the density of states of the sandwiched NM at Fermi energy in the two multilayers, as shown in the Fig. 2(b).
Though the out-of-plane spin torque shows less effect on the magnetic switching [
22], it can change the resistance in the structure via inverse spin Hall effect (ISHE). The out-of-plane
originates from a spin current along the
y axis with a spin polarization in the
x direction. When the spin current is absorbed by the sandwiched NM via spin-orbital coupling, it would produce an electron flow along the
z axis via ISHE. Therefore, the lateral resistance can be varied in the presence of thermal pumping torque.
The spin Seebeck torque can be considered as a local effect, which is directly related to the thermal magnetic damping
. Figure 3 shows
and
as functions of temperature in the Ag|Fe(16)|Ag(16)|YIG multilayer. At high temperatures, both
and
are proportional to
, and a deviation is observed at low temperatures. The temperature dependence is consistent with that estimated from Eq. (5). Therein, for
, where
is the magnetic coherence length, we obtain
; consequently,
and
. That is,
can decrease to one-thirty, when the temperature decreases one order. We found that
and
nJ/(K·m
2) in Ag|Fe(16)|Ag(16)|YIG at 300 K, considering
and taking [
5]
J·m
2. When Pt is the separate metal,
and
nJ/(K·m
2), as listed in Table 1. Comparatively,
would dominate over
in both the spin valves. By including the thermal torque in the LLG equation [
38], we obtain a temperature gradient Δ
TSW=12/14 K that is sufficient to switch the magnetization in Fe between antiparallel (AP) and parallel (P) configurations in the NM|Fe(16)|NM(16)|YIG multilayer with Ag/Pt spacer, as listed in Table 1, which is of the same order as that in the MgO-based MTJs [
10,
19,
22].
The origin of large and can be found from the inset of Fig. 3, which shows the energy-dependent α. We find that: i) is dominated by the diagonal part, and the off-diagonal part is several times smaller than the diagonal part. ii) is considerable larger than both and , indicating that the magnetic damping is anisotropic. The break in rotational symmetry is related to the structural symmetry and the existence of quantum states in the spin valves. , which corresponds to M2 rotates slightly towards x, and this rotation pumps spin current into M1 with a spin polarization along the y-axis. The pumped spin current can be absorbed by M1, which is along the x-axis, leading to an enhancement of the magnetic damping. When M2 rotates towards the y-axis, the pumped spin current is polarized along the x direction, and is hardly absorbed by M1, leading to small magnetic damping. As a result, the damping component becomes asymmetric.
When the sandwiched NM is sufficiently thick, the contribution from quantum effects would be very weak, and the magnetic damping α would be dominated by the spin mixing conductance at the NM|FI interface, as spin flip in the sandwiched NM is negligible. In Fig. 3, the and estimated from the spin mixing conductance are also given. Both and estimated from the spin mixing conductance are about smaller than that estimated from the calculations. This indicates that the quantum effect in multilayers cannot be ignored.
Interlayer exchange coupling is usually observed in multilayers. The effects of interlayer exchange coupling can be observed from Fig. 4, which shows the thickness dependence of
α on the sandwiched NM in Ag|Fe(16)|Ag(
x)|YIG multilayers. A slowly damped oscillation is observed in
, as the function of Ag thickness varied from 1 to
. As the thickness of the sandwiched Ag is below
, the oscillation period is approximately
, which is close to those in FM|Ag|FM and FM|Cu|FM [
49]. This result indicates that spin pumping is affected by the interlayer exchange coupling in the Ag spacer. At the same time, the oscillatory behavior in
α, also exists in
, thus providing another method for measuring the interlayer exchange coupling through the calculation of spin torque.
4 Summary
We investigated the enhanced magnetic damping and thermal spin transfer effects in NM|Fe|NM|YIG multilayer from first principles studies. The spin Seebeck torque of nJ/(K·m2) is found in the multilayers Ag/Pt spacer at 300 K, which is considerably larger than the thermal spin pumping torque of nJ/(K·m2). Based on the LLG equation, we estimated that a temperature bias of across the sandwiched NM would be sufficient to reverse the magnetic configurations circularly in a multilayer at room temperature. Moreover, estimated on the basis of first-principles, is considerably larger than that estimated from spin mixing conductance, indicating that the quantum effect in spin valves cannot be ignored. It is easy to measure the spin torque experimentally, therefore, quantum effects such as interlayer exchange, can be estimated from the spin torque in multilayers.
Higher Education Press and Springer-Verlag Berlin Heidelberg 2017