The inverse thermal spinorbit torque and the relation of the DzyaloshinskiiMoriya interaction to groundstate energy currents
Abstract
Using the Kubo linearresponse formalism we derive expressions to calculate the heat current generated by magnetization dynamics in magnets with broken inversion symmetry and spinorbit interaction (SOI). The effect of producing heat currents by magnetization dynamics constitutes the Onsager reciprocal of the thermal spinorbit torque (TSOT), i.e., the generation of torques on the magnetization due to temperature gradients. We find that the energy current driven by magnetization dynamics contains a contribution from the DzyaloshinskiiMoriya interaction (DMI), which needs to be subtracted from the Kubo linear response of the energy current in order to extract the heat current. We show that the expressions of the DMI coefficient can be derived elegantly from the DMI energy current. Guided by formal analogies between the Berry phase theory of DMI on the one hand and the modern theory of orbital magnetization on the other hand we are led to an interpretation of the latter in terms of energy currents as well. Based on abinitio calculations we investigate the heat current driven by magnetization dynamics in Mn/W(001) magnetic bilayers. We predict that fast domain walls drive strong ITSOT heat currents.
pacs:
72.25.Ba, 72.25.Mk, 71.70.Ej, 75.70.Tj1 Introduction
The interaction of heat current with electron spins is at the heart of spin caloritronics [1]. It leads to thermal spintransfer torques (STTs) on the magnetization in spin valves, magnetic tunnel junctions, and domain walls when a temperature gradient is applied [2, 3, 4, 5, 6, 7, 8]. While the thermal STT does not require spinorbit interaction (SOI), it only exists in noncollinear magnets. In spin valves and magnetic tunnel junctions this noncollinearity arises when the magnetizations of the free and fixed layers are not parallel, while in domain walls it arises from the continuous rotation of magnetization across the wall.
In the presence of SOI electric currents and heat currents can generate torques also in collinear magnets: In ferromagnets with broken inversion symmetry the socalled spinorbit torque (SOT) acts on the magnetization when an electric current is applied (Figure 1a) [9, 10, 11, 12, 13, 14, 15, 16, 17] . The inverse spinorbit torque (ISOT) consists in the production of an electric current due to magnetization dynamics (Figure 1b) [18, 19, 20]. The application of a temperature gradient results in the thermal spinorbit torque (TSOT) (Figure 1c) [21]. TSOT and SOT are related by a Mottlike expression [22].
In this work we discuss the inverse effect of TSOT, i.e., the generation of heat current due to magnetization dynamics in ferromagnets with broken inversion symmetry and SOI (Figure 1d). We refer to this effect as inverse thermal spinorbit torque (ITSOT). While the SOT is given directly by the linear response of the torque to an applied electric field [16], expressions for the ITSOT are more difficult to derive because the energy current obtained from the Kubo formalism contains also a groundstate contribution that does not contribute to the heat current. Analogous difficulties are known from the case of the inverse anomalous Nernst effect, i.e., the generation of a heat current transverse to an applied electric field [23]. In this case the energy current obtained from the Kubo formalism contains besides the heat current also the materialdependent part of the Poynting vector, where is the orbital magnetization. This energy magnetization does not contribute to the heat current and needs to be subtracted from the Kubo linear response [23, 24, 25].
When inversion symmetry is broken in magnets with SOI the expansion of the free energy in terms of the magnetization direction and its gradients contains a term linear in the gradients of magnetization, the socalled DzyaloshinskiiMoriya interation (DMI) [26, 27]:
(1) 
where is the position and the index runs over the three cartesian directions, i.e., . The DMI coefficients can be expressed in terms of mixed Berry phases [22, 28]. DMI does not only affect the magnetic structure by energetically favoring spirals of a certain handedness but also enters spin caloritronics effects [29, 30]. Here, we will show that DMI gives rise to the groundstate energy current when magnetization precesses. This DMI energy current needs to be subtracted from the linear response of the energy current in order to obtain the ITSOT heat current.
This work is structured as follows. In section 2 we show that magnetization dynamics drives a groundstate energy current associated with DMI and we highlight its formal similarities with the materialdependent part of the Poynting vector. In section 3 we develop the theory of ITSOT. We derive the energy current based on the Kubo linearresponse formalism and subtract in order to extract the heat current. In section 4 we show that the expressions of DMI and orbital magnetization can also be derived elegantly by equating the energy currents obtained from linear response theory to and , respectively. In section 5 we present abinitio calculations of TSOT and ITSOT in Mn/W(001) magnetic bilayers.
2 Groundstate energy current associated with the DzyaloshinskiiMoriya interaction
To be concrete, we consider a flat cycloidal spin spiral propagating along the direction. The magnetization direction is given by
(2) 
where is the spinspiral wavenumber, i.e., the inverse wavelength of the spin spiral multiplied by . The free energy contribution given in (1) simplifies for the spin spiral of (2) as follows:
(3)  
where is the unit vector pointing in direction and we defined . Whether is nonzero or not depends on crystal symmetry. The tensor is axial and of second rank like the SOT torkance tensor [22]. Additionally, it is even under magnetization reversal, i.e., . Therefore, has the same symmetry properties as the even SOT torkance [16]. According to (3) the cycloidal spiral of (2) is affected by DMI if is nonzero. This is the case e.g. for magnetic bilayers such as Mn/W(001) and Co/Pt(111) (the interface normal points in direction), where also the component of the even SOT torkance tensor is nonzero [16, 22, 31, 32].
We consider now a Neeltype domain wall that moves with velocity in direction. The magnetization direction at time , which we denote by , is illustrated in Figure LABEL:figure1a. can be interpreted as a modification of ((2)), where the vector depends on position:
(4) 
Since the domain wall moves with velocity , the magnetization direction at position and time is given by
(5) 
In Figure LABEL:figure1 we discuss the magnetization direction at position at the three times , and . At time the domain wall is far away from . Therefore, the magnetization is collinear at and . At time the domain wall starts to arrive at . Consequently, the magnetization gradient becomes nonzero and thus . Due to the motion of the domain wall the DMI contribution to the free energy is time dependent: How much DMI free energy is stored at a given position in the magnetic structure is determined by the local gradient of magnetization, which moves together with the magnetic structure. The partial derivative of with respect to time is given by
(6) 
where in the last line is the component of the DMI energy current density
(7)  
By considering additionally spirals propagating in and direction we find that the general form of (6) is the continuity equation
(8) 
of the DMI energy current . According to (7) and (8) the energy current is driven by magnetization dynamics and its sources and sinks signal the respective decrease and increase of DMI energy density. When we compute the energy current driven by magnetization dynamics in section 3 we therefore need to be aware that this energy current contains in addition to the ITSOT heat current that we wish to determine. Thus, we need to subtract from the energy current in order to extract the ITSOT heat current.
It is reassuring to verify that the materialdependent part of the Poynting vector, which needs to be subtracted from the energy current to obtain the heat current in the case of the inverse anomalous Nernst effect [23], can be identified by arguments analogous to the above. We sketch this in the following. The energy density due to the interaction between orbital magnetization and magnetic field is given by
(9) 
We assume that the magnetic field is of the form
(10) 
i.e., the magnetic field at time can be obtained from the magnetic field at time by shifting it by , as illustrated in Figure LABEL:figure2. Additionally, we assume that the orbital magnetization is of the same form, i.e., . Consequently, also . can be expressed as in terms of the vector potential
(11) 
Due to the motion of the profile of the energy density in (9) changes as a function of time. The partial derivative of with respect to time is
(12)  
where we used the Maxwell equation and valid in Weyl’s temporal gauge with scalar potential set to zero. Thus,
(13) 
with
(14) 
as expected.
In the following we discuss several additional formal analogies and similarities between DMI, classical electrodynamics and orbital magnetization. We introduce the tensors and with elements
(15) 
and
(16) 
to quantify the noncollinearity of . and are related through the matrix
(17) 
as . The free energy can be expressed in terms of and as follows:
(18)  
where we defined . Similarly, in (7) can be expressed in terms of as
(19) 
The energy density in (9) involves the curl of the vector potential , while the materialdependent part of the Poynting vector, i.e., , involves the timederivative of . Similarly, the spatial derivatives enter in (18) via the tensor while the temporal derivative enters in (19). Thus, in the theory of DMI the magnetization direction plays the role of an effective vector potential.
The curl of orbital magnetization constitutes a bound current that does not contribute to electronic transport. Hence it needs to be subtracted from the linear response electric current driven by gradients in temperature or chemical potential in order to obtain the measurable electric current [23]. Similarly, the spatial derivatives that result from the presence of gradients in temperature or chemical potential constitute torques that are not measurable and need to be subtracted from the total linear response to temperature or chemical potential gradients in order to obtain the measurable torque [22]. Table 1 summarizes the formal analogies and similarities between the orbital magnetization and DMI.

OM DMI ’vector potential’ ’magnetic’ field energy density ’electric’ field energy current bound property
3 Inverse thermal spinorbit torque (ITSOT)
In ferromagnets with broken inversion symmetry and SOI, a gradient in temperature leads to a torque on the magnetization, the socalled thermal spinorbit torque (TSOT) [22, 21]:
(20) 
The inverse thermal spinorbit torque (ITSOT) consists in the generation of heat current by magnetization dynamics in ferromagnets with broken inversion symmetry and SOI. The effect of magnetization dynamics can be described by the timedependent perturbation to the Hamiltonian [16]
(21) 
where is the torque operator. is the exchange field, i.e., the difference between the potentials of minority and majority electrons. is the spin magnetic moment operator, is the Bohr magneton and is the vector of Pauli spin matrices. The energy current driven by magnetization dynamics is thus given by
(22) 
where the tensor with elements
(23) 
describes the Kubo linear response of the energy current operator
(24) 
to magnetization dynamics. is the chemical potential, is the velocity operator and the retarded energycurrent torque correlationfunction is given by
(25) 
In (23) we take the limit frequency , which is justified when the frequency is small compared to the inverse lifetime of electronic states, which is satisfied for magnetic bilayers at room temperature and frequency in the GHz range.
Within the independent particle approximation (23) becomes , with
(26) 
where and are the retarded and advanced singleparticle Green functions, respectively. is the Fermi function. contains scatteringindependent intrinsic contributions and, in the presence of disorder, additional disorderdriven contributions. The intrinsic Berrycurvature contribution is given by
(27) 
where
(28) 
and
(29) 
and are the Bloch wavefunctions with corresponding band energies , , and is the number of points.
As discussed in section 2 we subtract ((7)) from in order to obtain the heat current :
(30) 
with
(31) 
Inserting the Berrycurvature expression of DMI [22, 28]
(32) 
we obtain for the intrinsic contribution
(33)  
where . Using
(34) 
where is the Heaviside unit step function, we can rewrite (33) as
(35) 
Here,
(36) 
is the intrinsic SOT torkance tensor [16, 22] at zero temperature as a function of Fermi energy and is the elementary positive charge.
The intrinsic TSOT and ITSOT are even in magnetization, i.e., . (26) contains an additional contribution which is odd in magnetization, i.e., , and which is given by
(37) 
where is the odd contribution to the SOT torkance tensor as a function of Fermi energy [16]. The total coefficient, i.e., the sum of all contributions, is related to the total torkance for magnetization in direction by
(38) 
It is instructive to verify that the ITSOT described by (38) is the Onsagerreciprocal of the TSOT ((20)), where [22]
(39) 
Comparison of (38) and (39) yields
(40) 
and thus
(41) 
where is the thermal conductivity tensor and describes Gilbert damping and gyromagnetic ratio [18]. As expected, the response matrix
(42) 
satisfies the Onsager symmetry .
4 Using the groundstate energy currents to derive expressions for DMI and orbital magnetization
The expression (32) for the DMIspiralization tensor was derived both from semiclassics [28] and static quantum mechanical perturbation theory [22]. Alternatively, the expression of can also be obtained elegantly and easily by invoking the third law of thermodynamics: For the ITSOT must vanish, , because otherwise we could pump heat at zero temperature and thereby violate Nernst’s theorem. Hence, according to (31). In other words, at the energy current density in (22) is identical to the DMI energy current density because the heat current is zero. Thus, at we obtain from (27)
(43) 
which agrees with (32) at .
Similarly, we can derive the expression of orbital magnetization from the energy current discussed in (14): For the inverse anomalous Nernst effect (i.e., the generation of a transverse heat current by an applied electric field) has to vanish according to the third law of thermodynamics. Hence, the energy current driven by an applied electric field at does not contain any heat current and is therefore identical to . We introduce the tensor to describe the linear response of the energy current to an applied electric field , i.e., . We describe the effect of the electric field by the vector potential and take the limit later. The Hamiltonian density describing the interaction between electric current density and vector potential is , from which we obtain the timedependent perturbation
(44) 
Introducing the retarded energycurrent velocity correlationfunction
(45) 
we can write the elements of the tensor as
(46) 
This allows us to determine as , where the intrinsic Berrycurvature contribution to the response tensor is given by
(47)  
with
(48) 
and
(49) 
From we obtain
(50) 
It is straightforward to verify that given by (50) agrees to the expressions for orbital magnetization derived from quantum mechanical perturbation theory [33], from semiclassics [23], and within the Wannier representation [34, 35].
Combining the third law of thermodynamics with the continuity equations (8) and (13) provides thus an elegant way to derive expressions for and at . We can extend these derivations to if we postulate that the linear response to thermal gradients is described by Mottlike expressions. In the case of the TSOT this Mottlike expression is (39), while it is [23, 36, 37]
(51) 
in the case of the anomalous Nernst effect, where is the zerotemperature anomalous Hall conductivity as a function of Fermi energy and the anomalous Nernst current due to a temperature gradient in direction is . While (39) and (51) were, respectively, derived in the previous section and in [23], we now instead consider it an axiom that within the range of validity of the independent particle approximation the linear response to thermal gradients is always described by Mottlike expressions. Thereby, the derivation in the present section becomes independent from the derivation in the preceding section. Applying the Onsager reciprocity principle to (39) and (51) we find that the ITSOT and the inverse anomalous Nernst effect are, respectively, described by (38) and by
(52) 
Employing the general identity (34) (but in contrast to section 3 we now use it backwards) we obtain
(53)  
from (35) and, similarly, (52) can be written as
(54)  
The finite expressions of and are now easily obtained, respectively, by subtracting the ITSOT heat current given by (53) from the energy current (27) and by subtracting the heat current (54) from . This leads to (32) for the DMI spiralization tensor and to
(55) 
for the orbital magnetization. (55) agrees to the finite expressions of derived elsewhere [33, 23].
5 Abinitio calculations
We investigate TSOT and ITSOT in a Mn/W(001) magnetic bilayer composed of one monolayer of Mn deposited on 9 layers of W(001). The ground state of this system is magnetically noncollinear and can be described by the cycloidal spin spiral (2) [31]. Based on phenomenological grounds [38, 19] we can expand torkance as well as TSOT and ITSOT coefficients locally at a given point in space in terms of and :
(56)  
The coefficients , , ,…in this expansion can be extracted from magnetically collinear calculations. Analogous expansions of the TSOT and ITSOT coefficients are of the same form. Here, we consider only and , which give rise to the following contribution to the torque :
(57)  
where we used that for magnetization direction along it follows from symmetry considerations that , , and . The SOT in this system has already been discussed by us [16]. In order to obtain TSOT and ITSOT, we calculate the torkance for the magnetically collinear ferromagnetic state with magnetization direction set along