Transport Property and Spin-Orbit Torque in 2D Rashba Ferromagnetic Electron Gas.

In this paper, we investigate the spin-orbit torque and transport property in a 2D Rashba ferromagnetic electron gas. The longitudinal conductivity can be divided into two parts: the first term is determined by the charge density and is independent of the spin degrees of freedom. The second term depends on the two bands that spin in the opposite directions, and it is directly proportional to spin-orbit torque regardless of the band structure and temperature. This is a general and underlying relation between the transport property and spin-orbit torque. Moreover, we show the impacts of the spin-orbit coupling constant and Fermi energy on transverse conductivity and spin-orbit torque, which is helpful for relevant experiments.

1. Introduction

In recent years, spin–orbit torque has become one of the most attractive topics in the field of spintronics because of its great application prospects in magnetic information storage technology [1,2,3]. spin–orbit torque (SOT) is based on spin–orbit interaction (SOI), which uses the non-equilibrium spin accumulation induced by charge flow to generate a torque on a local magnetic moment [4,5], thus achieving the purpose of regulating magnetic storage units. SOT-based magnetic random access memory (SOT-MRAM) overcomes the disadvantages of STT-MRAM, especially the fact that it separates the read and write paths, so it has a higher read and write speed and lower power consumption than STT-MRAM [6,7]. Relevant studies show that SOT-MRAM can achieve ultra-fast information writing, which is reduced from the tens of nanoseconds required by the original STT-MRAM to less than 10 ns, while the power consumption of the device is further reduced [8,9].

The essence of spin–orbit torque is that the directional movement of electrons produces a non-equilibrium spin accumulation in the spin–orbit coupling system, by means of the s–d interaction, the spin accumulation of conducting electrons applies a torque on the local magnetic moment [10,11]. Meanwhile, this directional movement of electrons also generates a charge current in general. Thus, the current and spin–orbit torque are both manifestations of the non-equilibrium transport properties of the SOC system, there must be some relations between them. Literature have showed that the magnitude of spin–orbit torque is often proportional to the current density [12,13,14,15]. Because the spin accumulation and charge current are proportional to the external electric field under linear transport conditions, the ratio depends on the band structure of the system [4,5]. There is no clear conclusion about the underlying relationship between spin–orbit torque and conductivity.

Two-dimensional Rashba ferromagnetic electron gas is an ideal platform for investigating spin–orbit torque [16,17], and the latter is often realized at heavy metal/ferromagnetic interfaces [18,19]. These systems have both Rashba spin–orbit interaction and ferromagnetism, and interesting transport phenomena such as spin–orbit torque and anomalous Hall effect have been found [20,21]. In this paper, we use 2D Rashba ferromagnetic electron gas as an example to explore the relationship between spin–orbit torque and transport properties, as well as the regulation of spin–orbit torque. In Section 2, we show the energy splitting caused by spin–orbital interactions; in Section 3, we study the transport properties of the system, including longitudinal conductance and intrinsic anomalous Hall conductivity; in Section 4, we investigate the spin–orbit torque, and show its regulation in terms of the spin–orbit coupling constant and Fermi energy; Section 5 is the conclusion.

2. Model

The Hamiltonian of an electron in a 2D Rashba ferromagnetic electron gas is , where is the effective mass of the electron, is the direction of magnetization and we take an out-of-plane magnetization as: , J and are the constant of the s–d interaction and Rashba spin–orbit coupling. To be sure, the s–d interaction is generated between the electron gas (s electrons) and the local magnetic moments (d electrons), where the latter form a three-dimensional ferromagnetic layer [4,5]. Therefore, the Hamiltonian can be written as [16,17,21]:Solving the Schrodinger’s equation, we obtain where is the energy splitting of these two bands. Correspondingly, the wavefunctions are , where Then, the average spin of the state is , and the average spin of the state is .

Figure 1 shows the bands splitting caused by spin–orbit interaction. The spins of electrons in these two bands are in opposite directions, and the energy difference is . The band structure depends on the relative magnitudes of and J. If , where , and are both parabolic. If , has a maximum at and two bottoms. The specific band structures will affect transport property.

Figure 1

(a) Bands splitting with opposite spins; and bands diagram when (b) , (c) .

3. Transport Property

In a semiclassical transport equation, the velocity of electrons in band n can be expressed as , where is the traditional velocity and is the anomalous velocity arising from Berry curvature [22]. In our system, the band splits into , therefore the velocities of electrons are: where is the external electric field and the Berry curvature [22] can be calculated from wavefunctions: .

In a weak external electric field, the distribution functions are , where are the equilibrium distribution functions, and is the first-order perturbation caused by an external electric field. By the relaxation time approximation, can be expressed as: and is the momentum relaxation time [4]. By the definition of current density , longitudinal and transverse current densities are: and To simplify matters, the equilibrium distribution functions can be expressed via the step function at 0K: , where is the Fermi energy. Thus, the longitudinal conductivity is where the first term of Equation (7) can be calculated as: where is calculated with the formulas: . Allowing for the charge density , can be expressed as which has the same form with the conductivity of free electron gas. This conductivity is determined by the charge density of electrons, and independent of the spin degree of freedom.

The second term of Equation (7) is derived from the different velocities in two bands:This conductivity results from the different velocities of spin-up and spin-down electrons, and reveals the spin-polarized transport property. Considering the change in Fermi surface, we discuss the results in categories:

When , the results of equations are and , respectively. Thus, When , the result of equation is . Thus, When , only if , the equation has the solutions and (). Thus,

Figure 2 shows the longitudinal conductivity as a function of Fermi energy . With the increase in , the charge density n increases gradually. Thus, increases monotonously. At , the number of crossover points of the Fermi surface changes. When , the Fermi surface intersects with the top of at . Below this point, there are two Fermi wave vectors and ; above that point, increases with the only Fermi wave vector . Furthermore, when , the Fermi-surface intersects with the bottom of . Above this point, electrons in band participate in the conduction process. From Equation (7), we can see a negative contribution made by . Thus, decreases gradually.

Figure 2

Longitudinal conductivities and vs. the Fermi energy , where , the solid line is , and the dashed line is .

The transverse conductivity is just the intrinsic anomalous Hall conductance, which originates from the Berry curvature [17]. Accordingly, the transverse conductance is

Figure 3 shows the transverse conductivity as a function of Fermi energy for different band structures. According to Equation (6), depends on the equilibrium distribution functions , which means that all states below Fermi energy contribute to transverse conductivity. For , and are both parabolic. When , begins to make a negative contribution, and thus decreases. For , has a top at and two bottoms. The maximum of Berry curvature is just at , which means a maximum of transverse velocity. Therefore, the transverse conductivity mainly depends on the states around . When , increases with the increase in , because the states in around gradually contribute to . When , is almost unchanged. When , decreases because of the negative contribution by states in .

Figure 3

Transverse conductivity vs. the Fermi energy , where (a) and (b) , the red line and blue line are bandstructures of and respectively.

4. Spin–Orbit Torque

The spin–orbit torque is the result of s–d interaction between the local magnetic moment and spin accumulation of conduction electrons, and it is always expressed as [23,24]: where is the spin accumulation and can be calculated by definition [4,15]:According to the distribution function above, the spin accumulation along the y axis is Compared to Equation (9), is directly proportional to . , and this specific value is a constant which has no connection with the band structure and temperature. Thus, the spin accumulation and have the same variation trend in terms of .

Figure 4 shows the spin accumulation as a function of the Rashba SOC constant . When is small, the band structure is still parabolic. Therefore, researchers usually treat the spin–orbit interaction as a perturbation [4,5,25]. The interaction between conduction electrons and the local magnetic moment can be regarded as a spin–orbit effective field acting on local magnetic moment [17,26]: , where is the drift of in an external electric field. is calculated based on the ferromagnetic electron gas model regardless of the band splitting caused by spin–orbit interaction. Consequently, the spin–orbit torque is a direct proportional function of . When is large, is smaller than the result based on the ferromagnetic electron gas model. Our numerical results show the deviation from the proportional relation.

Figure 4

Spin accumulation along the y axis vs. Rashba SOC constant , where the solid line is the result based on the split bands and the dashed line is based on the band of ferromagnetic electron gas.

Similarly, the spin accumulation along the z axis is This spin accumulation depends on the equilibrium distribution function and is parallel to the local magnetic moment; thus, it makes no contribution to the spin–orbit torque.

5. Conclusions

In this paper, we investigated the spin–orbit torque and transport property in a 2D Rashba ferromagnet. The main conclusions are as follows:

The longitudinal conductivity can be divided into two parts: the first term is determined by the charge density and is independent of the spin degrees of freedom. The second term depends on the two bands that spin in opposite directions, and this reflects the spin-polarized transport property of the system.

The spin–orbit torque is directly proportional to the second term of longitudinal conductivity, because spin–orbit torque is just caused by the spin-polarized transport. Although this conclusion is obtained at 0K and using parabolic bands, it is suitable for general linear transport cases. This proportionate relation reveals the underlying connection between spin–orbit torque and conductivity.

Moreover, we demonstrate the impacts of the spin–orbit coupling constant and Fermi energy on spin–orbit torque. When these constants are experimentally changed, the spin–orbit torque can be adjusted accordingly. The results are helpful for relevant experiments.