Writing using spin–orbit torque (SOT) has been widely investigated in the field of magnetic random-access memory (MRAM). Heavy metal (HM)/CoFeB/MgO is the core of this SOT-MRAM structure. The heterostructure consisting of Ta as the spin current source and CoFeB/MgO as the perpendicular magnetic anisotropy (PMA) material is the most researched structure, owing to its high tunneling magnetoresistance ratio. However, Ta is difficult to integrate into the CMOS process due to its poor thermal stability against annealing at temperatures greater than 350 °C. Currently, β-tungsten (W) is the only HM with the CoFeB/MgO system, which can provide both thermal stability and SOT switching. Nevertheless, to achieve the high resistive β phase of W is a challenging task. Here, we report another material rhenium (Re) capable of providing thermally stable PMA up to temperature 425 °C with a perpendicular anisotropic field greater than 5000 Oe; Re possesses a spin hall angle (ϴSH) of 0.065 ± 0.003, and SOT switching can be achieved with a current density around 1.36 × 1011 A/m2. Our findings pave a new avenue for the material design of perpendicular SOT-based MRAM.

Perpendicular magnetic tunnel junction (p-MTJ)-based magnetic random-access memory (MRAM) is the front-runner of the next generation non-volatile memory due to its robust magnetization states, which are undisturbed by external perturbation except for a strong magnetic field.1–3 Manipulation of magnetization in p-MTJ using electric current [either spin transfer torque (STT) or spin–orbit torque (SOT)] has already pushed the rise of p-MTJ MRAMs into commercial productions.2–5 Producing commercial MRAM products with high density requires the free layer of the MTJ to keep perpendicular magnetic anisotropy (PMA) after annealing up to a temperature of at least 400 °C, given that the back-end semiconductor processing is performed at a temperature of around 400 °C.6,7 Thus, finding material systems with thermally stable PMA is of supreme importance. So far, heavy metal (HM)/CoFeB/MgO-based p-MTJs have been commonly used due to their sufficiently high tunneling magnetoresistance (TMR) value.8–11 Although a range of materials, including Ta, Mo, W, Hf, and Pd, have been shown to have PMA, only Mo and β-W can still reveal PMA after annealing at 400 °C or higher temperature.6,7 The origin of PMA in HM/CoFeB/MgO is believed to be mainly originating from the Fe–O orbital hybridization at the CoFeB/MgO interface and boron (B) out-diffusion.6,7,12–15 Furthermore, the property of HM is another crucial factor affecting anisotropy and TMR, which needs to promote the CoFeB crystallization along MgO [002].16–19 Thermal stability of CoFeB against the annealing temperature of 400 °C is shown to depend strongly on the crystal structure and oxygen affinity of HM metals.6,7 It has been reported that in the case of Ta as the HM layer, the formation enthalpy of Ta oxide is comparable to that of MgO, which may lead to diffusion of Ta toward the CoFeB/MgO interface during high temperature annealing and thus causes severe dead layer formation.6,7 Moreover, thin Ta is an amorphous layer, which is a high energy state as compared to the crystalline ones, leading to further enhancement of Ta diffusion into CoFeB. Severe Ta diffusion destroys the PMA.6,7

Mo and β-W crystallize into a low energy BCC structure at the as-deposited state and also have lower formation enthalpy for their oxides as compared to that of Ta oxides. So, they facilitate the growth of MgO along the [002] direction during the deposition and allow CoFeB to crystallize along MgO [002] during the annealing process.6,7 Because of their low-energy stable BCC phase and lower oxygen affinity than that of Mg, they can suppress interdiffusion during high temperature annealing and still keep strong PMA.6,7 Moreover, W and Mo are well-known refractory metals with high melting points, which might be another factor to slow down the diffusion.

One of the critical factors for MRAM cells is thermal stability ∆ = KeffV/KBT, where Keff is the effective magnetic anisotropy energy density, V is the volume of the magnetic free layer, KB is Boltzmann’s constant, and T is the absolute temperature.20–22 In addition to enhancing Keff by the underlayer, another straightforward way to increase ∆ is to increase the volume of the free layer, so this is where the concept of composite free layers comes in. Two CoFeB layers separated by a very thin HM layer have been reported to substantially increase the volume of the free layer and, thus, the thermal stability.20–24 In addition to promoting perpendicular anisotropy, the HM may also need to play a role as the spin current source for writing in the SOT-MRAM.25–27 Writing using the three terminal SOT, which utilizes the spin hall effect (SHE) and/or the Rashba effect, provides protection for the insulator against the thin oxide barrier (MgO) breakdown due to the separate write/read paths.25–27 In addition, SOT writing may provide much faster switching than the STT writing.25–27 So, having material systems, which can simultaneously provide SHE and thermally stable PMA against high temperature annealing, is important for the development of commercial SOT-MRAMs. Currently, β-W is the only HM serving well for these purposes.28,29 The switching current density can be reduced to 1.1 × 1011 A/m228 However, achieving a highly resistive β phase may pose a different challenge for mass production.28,29 Moreover, very high resistivity of β-W (∼250 µΩ cm) may further lead to strong Joule heating during writing.

Here, we report another HM, Re, which provides both a thermally stable perpendicular CoFeB free layer with perpendicular anisotropy field HK larger than 5000 Oe after 400 °C annealing and SOT switching of the perpendicular CoFeB with a critical switching current density of 1.36 × 1011 A/m2 under the in-plane field HX of 50 Oe. Unlike W,28,29 Re does not require a lot of tuning of deposition conditions to achieve a specific phase to have PMA and SHE. Moreover, a perpendicular composite free layer with thin Re as the spacer separating two CoFeB layers is demonstrated to have double volume of the ferromagnet (FM) and even higher HK of 6000 Oe than that of a single FM. To further quantize the SHE in the Re/CoFeB/MgO structure, we employed a loop shift method30–32 to determine the spin hall angle of Re to be 0.065 ± 0.003. Our results provide an additional material option for SOT-MRAM application.

Si/SiO2/Re(5)/Co20Fe60B20(0.9–1.5)/MgO(2)/Ta(3.5) (in nm) samples were deposited to study the dependence of PMA on CoFeB thickness. Si/SiO2/Re(t)/CoFeB(1.1)/MgO(2)/Ta(3.5) samples with t = 1.1–10 nm were deposited to investigate how the Re thickness affects the PMA properties. Si/SiO2/Ta(0.4)/Re(t)/CoFeB(1.07)/MgO(2)/Ta(3.5) samples with t = 1.1–10 nm were studied to check the effects of the Ta seed layer on the magnetic properties of CoFeB. Finally, the composite free layer structures Si/SiO2/Ta(3.5)/MgO(1.5)/CoFeB(1.1)/Re(t)/CoFeB(1.1)/MgO(1.5)/Ta(3.5) were deposited. The base pressure of the sputtering system was 1.0 × 10−7 Torr, and the metallic layers were deposited using DC magnetron sputtering, and MgO was deposited using RF sputtering. All the samples were annealed at 300, 350, 400, and 425 °C under pressure lower than 1 × 10−5 Torr for 30 min. Magnetic properties were studied using the vibrating sample magnetometer (VSM). X-ray diffraction (XRD) analysis was performed with the x-ray wavelength of 0.254 nm. To study the SOT properties, we used photolithography and Ar ion milling to fabricate the typical Hall cross structure. The width and length of the current channel are 10 and 40 µm, respectively, while those of the voltage channel are 3 and 40 µm, respectively.

For annealing temperatures at 350 °C or higher, the samples Si/SiO2/Re(5)/Co20Fe60B20(t)/MgO(2)/Ta(3.5) (in nm) in Fig. 1(a) with the CoFeB thickness between 0.9 and 1.4 nm all exhibited perpendicular magnetic anisotropy. The variations in HK with the CoFeB thickness after annealing at 400 °C are shown in Fig. 1(b). The initial increase in HK with the thickness may be due to the existence of the dead layer. When the thickness further increases, HK is reduced, reflecting the perpendicular anisotropy that originates from the interface. The VSM measurements for the samples Re(5 nm)/CoFeB(1.1)/MgO annealed at 400 °C are shown in Fig. 1(c). HK can be as high as 5500 Oe after annealing at 350 °C and then gradually decreases with a further increase in the annealing temperature. Good PMA with HK of 4500 Oe can still be obtained after annealing at 425 °C [Fig. 1(d)]. Other properties including the dead layer thickness (tdead), saturation magnetization (MS), and interfacial anisotropy (KS) of samples annealed at 400 °C are shown in Figs. 2(a) and 2(b). tdead and MS were calculated from the plot MStCoFeB vs tCoFeB, as indicated in Fig. 2(a) for 400 °C. MS and tdead for the as-deposited samples were 939 emu/cc and 0.51 nm, respectively. The higher MS after annealing is due to the crystallization of CoFeB and B out-diffusion.

FIG. 1.

Tuning magnetic anisotropy of the samples Si/SiO2/Re(5)/Co20Fe60B20(t)/MgO(2)/Ta(3.5) (in nm). (a) Sample structures, (b) HK at different CoFeB thicknesses annealed at 400 °C, (c) in-plane and out-of-plane MH curves for the CoFeB thickness of 1.1 nm annealed at 400 °C, and (d) highest HK obtained at different annealing temperatures.

FIG. 1.

Tuning magnetic anisotropy of the samples Si/SiO2/Re(5)/Co20Fe60B20(t)/MgO(2)/Ta(3.5) (in nm). (a) Sample structures, (b) HK at different CoFeB thicknesses annealed at 400 °C, (c) in-plane and out-of-plane MH curves for the CoFeB thickness of 1.1 nm annealed at 400 °C, and (d) highest HK obtained at different annealing temperatures.

Close modal
FIG. 2.

Quantitative analysis of magnetic properties for samples annealed at different temperatures. (a) Areal saturation magnetization vs CoFeB thickness for 400 °C annealing to calculate MS and the dead layer at 400 °C, (b) calculation of KS after 400 °C annealing, the solid line is the linearly fitted plot, (c) dead layer thickness after different annealing temperatures, (d) saturation magnetization after different annealing temperatures, and (e) interfacial anisotropies after different annealing temperatures.

FIG. 2.

Quantitative analysis of magnetic properties for samples annealed at different temperatures. (a) Areal saturation magnetization vs CoFeB thickness for 400 °C annealing to calculate MS and the dead layer at 400 °C, (b) calculation of KS after 400 °C annealing, the solid line is the linearly fitted plot, (c) dead layer thickness after different annealing temperatures, (d) saturation magnetization after different annealing temperatures, and (e) interfacial anisotropies after different annealing temperatures.

Close modal

KS was calculated using the equation Keff = (Kb − 2πMS2) + KS/tCoFeB, where volume anisotropy (Kb) is assumed to be negligible, and KS at different temperatures was calculated from the plot of Keff × tCoFeB vs tCoFeB, as shown in Fig. 2(b) for 400 °C annealing as an example. Finally, the temperature dependences of tdead, MS, and KS are shown in Figs. 2(c)2(e), respectively. KS is about 1.98, 1.73, and 1.62 erg/cm2 for samples annealed at 350, 400, and 425 °C, respectively [Fig. 2(e)]. KS is reported to be 1.7 and 2.05 erg/cm2 for W and Mo buffers, respectively, after 400 °C annealing.6,7 Our results of HK, MS, and tdead are comparable to those reported in the case of W/CoFeB/MgO after 400 °C annealing.6 

Next, the Re thickness effect on the magnetic properties was studied. The hysteresis loops measured by the VSM of the samples Si/SiO2/Re(t)/CoFeB(1.1)/MgO(2)/Ta(3.5) in Fig. 3(a) revealed that the PMA only existed with the Re thickness ranging between 4 and 8 nm after annealing. For very thin Re (3 nm or less) samples, both in-plane and out-of-plane hysteresis curves look like S-shaped hard axes; the MH loop for 400 °C annealed Re(1.1)/CoFeB(1.1)/MgO(2)/Ta(5) sample is shown in Fig. 3(b) as an example. HK reduced significantly for samples with the Re thickness of 9 nm or higher. The HK dependence on the Re thickness is shown in Fig. 3(c) 

FIG. 3.

Magnetic properties for different Re thicknesses. (a) Structure, (b) M–H loop after 400 °C annealing for the samples Re(1.125)/CoFeB(1.1)/MgO(2)/Ta(3.5), (c) HK vs Re thickness for samples annealed at 350 and 400 °C, and (d) M–H loop for 400 °C annealed samples Ta(0.4 nm)/Re(1.1)/CoFeB(1.1)/MgO(2)/Ta(3.5). The magnetization is normalized with respect to Re(5)/CoFeB(1.1)/MgO(2)/Ta(3.5).

FIG. 3.

Magnetic properties for different Re thicknesses. (a) Structure, (b) M–H loop after 400 °C annealing for the samples Re(1.125)/CoFeB(1.1)/MgO(2)/Ta(3.5), (c) HK vs Re thickness for samples annealed at 350 and 400 °C, and (d) M–H loop for 400 °C annealed samples Ta(0.4 nm)/Re(1.1)/CoFeB(1.1)/MgO(2)/Ta(3.5). The magnetization is normalized with respect to Re(5)/CoFeB(1.1)/MgO(2)/Ta(3.5).

Close modal

To further understand why the thin Re layer cannot provide good PMA, we deposited a very thin Ta with the thickness of 0.4 nm as a seed layer for Re and revealed that even very thin Re (1.1–3 nm) can have thermally stable PMA with HK close to 5000 Oe after 400 °C annealing. The MH curve for 1.1 nm Re with the 0.4 nm Ta seed layer is shown in Fig. 3(d). However, for the samples with thicker Re (9 nm or higher), the perpendicular anisotropy is still poor even with the Ta seed layer.

To test if Re can be used for the interlayer of the composite free layer, we prepared a composite free layer with the structures Si/SiO2/Ta(3.5)/MgO(1.5)/CoFeB(1.1)/Re(t)/CoFeB(1.1)/MgO(1.5)/Ta(3.5) with t = 0.3–1 nm, as shown in Fig. 4(a). All samples show PMA and ferromagnetic coupling between the two CoFeB layers in this Re interlayer thickness range after 400 °C annealing. We obtained maximum net magnetization for the Re thickness of 0.47 nm [Fig. 4(b)] with nearly double net magnetization and slightly higher HK of ∼6000 Oe than a single CoFeB. The higher HK is possibly due to the double Fe–O interface.

FIG. 4.

Composite free layer with the interlayer Re thickness t = 0.3–1 nm annealed at 400 °C. The maximum MSt was obtained for Re = 0.47 nm. (a) Structure for the composite free layer with ferromagnetic coupling between the two CoFeBs and (b) M–H curves for the Re thickness of 0.47 nm; the magnetization is normalized with respect to a single ferromagnet thickness (1.1 nm CoFeB).

FIG. 4.

Composite free layer with the interlayer Re thickness t = 0.3–1 nm annealed at 400 °C. The maximum MSt was obtained for Re = 0.47 nm. (a) Structure for the composite free layer with ferromagnetic coupling between the two CoFeBs and (b) M–H curves for the Re thickness of 0.47 nm; the magnetization is normalized with respect to a single ferromagnet thickness (1.1 nm CoFeB).

Close modal

Next, SOT studies were done by fabricating the typical hall cross geometry from the samples Re(t)/CoFeB(1.1)/MgO(2)/Ta(5), t = 1.1–5 nm, and a Ta seed layer of 0.4 nm was used to promote PMA for t = 1.1–3 nm. The SOT measurement setup and spin current generation due to SHE are shown in Fig. 5(a) with 5 nm thick Re. We applied 0.1 ms current pulses with an increase in the magnitude to achieve SOT switching. As shown in Fig. 5(b), we achieved full magnetization switching under HX = 50 Oe with a critical current density of 3 × 1011 A/m2 for Re(5)/CoFeB(1.1)/MgO(2)/Ta(5). The critical current density was calculated by considering the proportion of the current flowing into Re. We also performed the HX field dependence of SOT switching by applying HX up to 800 Oe, as shown in Fig. 5(c). Under this field range, the critical current density is slightly reduced with an increase in HX, which is a typical SOT switching behavior.

FIG. 5.

SOT switching in Re(5)/CoFeB(1.1)/MgO(2)/Ta(5). (a) Current induced SOT switching using SHE and the measurement setup in our structure, (b) SOT switching under HX = ±50 Oe, (c) SOT under different HX, and (d) JC for SOT switching for different Re thicknesses.

FIG. 5.

SOT switching in Re(5)/CoFeB(1.1)/MgO(2)/Ta(5). (a) Current induced SOT switching using SHE and the measurement setup in our structure, (b) SOT switching under HX = ±50 Oe, (c) SOT under different HX, and (d) JC for SOT switching for different Re thicknesses.

Close modal

Based on the switching polarity, we may know that the spin hall angle of Re has the same sign as W, that is, it is negative. Finally, we investigated how the Re thickness affected SOT switching. We found that even 1.1 nm Re could switch CoFeB. The dependence of the critical current density JC on the Re thickness is shown in Fig. 5(d). JC is slightly increased with an increase in the Re thickness and ranges from 1.36 × 1011 to 3 × 1011 A/m2.

To further quantify the SOT behavior of Re, we performed the loop shift measurement to calculate the damping-like (DL) efficiency (ξDL) of Re.30–32 The measurement setup is shown in Fig. 6(a). The damping-like (DL) effective field generated due to SHE in heavy metals can be obtained by applying a series of different magnitudes of DC current pulses while sweeping HZ under different HX fields [Fig. 6(b)]. From the loop shift, we can calculate the DL efficiency using the formula

where s=tFMρHMtHMρFM. Here, e is the elementary charge, is reduced Planck’s constant, u0 is the permeability of vacuum, MS is the saturation magnetization of the ferromagnet, tFM is the thickness of the ferromagnet, w is the width of the current channel, tHM is the thickness of the HM, and HzEffI is the saturated effective field per unit current calculated from the loop shift measurement. Figure 6(b) shows the loop shift corresponding to IDC = ±9 mA for HX = 680 Oe, and Fig. 6(c) shows the effective field (net shift) generated by different IDC applied for HX = 680 Oe. From the plot in Fig. 6(c), we can calculate the effective field per unit current under different HXHzEffI=4Oe/mAatHX=680Oe. Finally, HzEffI calculated in Fig. 6(c) is plotted as a function of HX, as shown in Fig. 6(d). The slope will saturate at a given HX = HDMI and the saturated value of HzEffI can be used to calculate ξDL using the above equation. The calculated ξDL was −0.065 ± 0.003. The DL efficiency is comparable to Pt with ξDL + 0.07 − +0.1233,34 but less than that of Ta (∼−0.15)34,35 and β-W (∼−0.21).28 HDMI in the Re/CoFeB structure is about 400 Oe based on the measurement shown in Fig. 6(d). ϴSH can be calculated from ξDL using the equation ξDL = TintϴSH, where Tint is the interfacial spin transparency. When Tint is assumed to be unity, the ϴSH is equal to −0.065 ± 0.003. In reality, the Tint is typically less than 1, so this ϴSH value can be underestimated.

FIG. 6.

Loop shift measurement. (a) Optical micrograph of the Hall cross device and the setup for the loop shift measurement, (b) illustration of the loop shift measurement under HX = 680 Oe for the current ±9 mA, (c) HZeff/I calculation under HX = 680 Oe, and (d) HZeff/I vs HX. The saturated value of HZeff/I gives the DL efficiency, while HX at which HZeff/I is saturated gives the HDMI field.

FIG. 6.

Loop shift measurement. (a) Optical micrograph of the Hall cross device and the setup for the loop shift measurement, (b) illustration of the loop shift measurement under HX = 680 Oe for the current ±9 mA, (c) HZeff/I calculation under HX = 680 Oe, and (d) HZeff/I vs HX. The saturated value of HZeff/I gives the DL efficiency, while HX at which HZeff/I is saturated gives the HDMI field.

Close modal

We believe that the origin of PMA in Re/CoFeB/MgO still mainly comes from the Fe–O orbital hybridization at the CoFeB/MgO interface. The magnetic properties of CoFeB (Ms, dead layer, and HK) are affected by the interdiffusion between CoFeB and Re. The wide angle XRD and grazing incidence x-ray diffraction (GIXD) analysis of Re, as shown in Figs. 7(a) and 7(b), respectively, reveal that the as-deposited 5 nm Re is amorphous; therefore, during annealing, CoFeB can easily crystallize along the MgO [002], and boron can be absorbed by the Re underlayer. The GIXD measurements shown in Fig. 7(b) reveal that after annealing, 5 nm Re is partially crystallized into the hexagonal close-packed (HCP) crystal structure, which is a stable and lower energy state than the as-deposited amorphous state. Thus, we believe that this crystallization suppresses the interdiffusion between Re and CoFeB during high temperature annealing.

FIG. 7.

XRD analysis of the rhenium thin film. (a) Wide angle XRD analysis and (b) GIXD measurements at 0.5° x-ray incidence.

FIG. 7.

XRD analysis of the rhenium thin film. (a) Wide angle XRD analysis and (b) GIXD measurements at 0.5° x-ray incidence.

Close modal

Moreover, after W, Re is the second highest melting-point metal, which could further suppress interdiffusion of Re into CoFeB during high temperature annealing. In fact, a recent report also revealed good thermal stability for Co/Re multilayers.36 Consequently, we can get high MS and low dead layers after high temperature annealing for Re/CoFeB/MgO. In addition, the affinity between HM and oxygen is also critical. If the affinity of HM-O is high, the HM may have the tendency to diffuse to the MgO side and react with oxygen; therefore, the bonding of Fe–O is reduced, and the PMA is destroyed during high temperature annealing, as observed in the case of Ta.7 Regarding the oxygen affinity of Re to form oxides, the enthalpy of the formation of Re oxide is −144 kJ/mol, the least among Mo oxide (−196 kJ/mol), Ta oxide (−292 kJ/mol), and W oxide (−197 kJ/mol).37 Thus, among all these heavy metals, Re is least expected to “steal” oxygen from the CoFeB/MgO interface and least expected to diffuse toward the Fe–O interface, which makes CoFeB sustain its PMA after high temperature annealing.

A thin Re (≦3 nm) layer without the Ta seed layer shows the isotropic behavior of CoFeB for both in-plane and out-of-plane measurements, as shown in Fig. 3(b), but a thin Re layer with the Ta seed layer shows good PMA after high temperature annealing, as shown in Fig. 3(d). It is possible that a very thin Re film is not continuous, leading to loss of PMA similar to the results observed for thin W.6 On the other hand, when we use the Ta seed layer, the film becomes rather continuous38 probably because of much better adhesion of Ta on SiO2, so the PMA can be achieved. Thicker Re does not have good perpendicular anisotropy no matter if the Ta seed layer is used or not. The XRD analysis of the thick Re (10 nm), as shown in Fig. 7, shows that both as-deposited and annealed 10 nm Re films are well crystallized into the HCP structure. The well crystalized Re layer may hinder the B out-diffusion and suppress the crystallization of CoFeB along MgO (002). Another possible reason could be the increase in the roughness of Re with the increase in the thickness, which could further suppress the PMA.

The SOT in Re/CoFeB might be coming from the resistive amorphous/amorphous + HCP state of Re, as it has been shown that the highly resistive phase of HMs such as Ta and β-W possesses higher DL SOT efficiency.28,29 Moreover, when the Pt resistivity is tuned to a higher value by alloying or tuning deposition conditions, it can show larger ξDL32,39 Our critical current density of 1.34–3 × 1011 A/m2 and ξDL of −0.065 ± 0.003 are similar to Pt, but ξDL is smaller than Ta or W. Slightly lower JC was observed for a very thin Re that maybe ascribed to larger resistivity of Re at a very thin Re thickness. On the other hand, note that the requirement of HX in Re/CoFeB/MgO for switching is quite small, which can be ascribed to small HDMI of ∼400 Oe in the Re/CoFeB structure as compared to Pt/Co or Pt/CoFeB, which have HDMI of 5000 and 2500 Oe, respectively.30 The HDMI in Re/CoFeB/MgO is similar to that of typical β-W/CoFeB/MgO (∼300 Oe)40 and Ta/CoFeB/MgO (∼250–300 Oe),30 which can also be switched under rather small HX. Finally, the lower resistivity of Re (less than 100 µemu cm) as compared to that of Ta (∼200–270 µemu cm) and β-W (∼250 µemu cm) may relieve the concern of the Joule heating during writing.

In this paper, we present a new HM rhenium that could generate PMA in CoFeB/MgO systems after high temperature annealing up to 425 °C. With a thin Ta seed layer, even 1.1 nm Re can support the perpendicular CoFeB with HK close to 5000 Oe. The as-deposited thin Re layer is amorphous, which becomes partially crystallized during annealing. The initial amorphous Re may promote the B out-diffusion, and CoFeB crystallizes along MgO [002]. Later on, the partially crystallized Re and low oxygen affinity of Re may suppress the interdiffusion between CoFeB and Re, which sustains the PMA after high temperature annealing. Furthermore, we demonstrate that Re can be a spin current source for SOT switching with similar damping-like torque efficiency to that of Pt. Unlike β-W, Re does not need to form a specific phase, and it has lower resistivity to avoid Joule heating. So far, only β-W and Re have been demonstrated to simultaneously possess high thermal stability of perpendicular CoFeB up to 400 °C and generate the spin current for SOT switching. It may pave an alternative avenue for the material design of SOT-MRAM.

This work was supported by Solar Applied Materials Technology Corporation on sputtering targets and partially supported by the Ministry of Science and Technology of Taiwan (Grant Nos. 109-2218-E-007-018 and 110-2218-E-007-034). We also acknowledge the help of Professor Chi-Feng Pai, National Taiwan University, Taipei, Taiwan, and his Ph.D. student Yu-Hu Chen for helping out with the loop shift measurements.

The authors declare no conflicts of interest to disclose.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
S. A.
Wolf
,
D. D.
Awschalom
,
R. A.
Buhrman
,
J. M.
Daughton
,
S.
von Molnár
,
M. L.
Roukes
,
A. Y.
Chtchelkanova
, and
D. M.
Treger
,
Science
294
,
1488
(
2001
).
2.
K. C.
Chun
,
H.
Zhao
,
J. D.
Harms
,
T.-H.
Kim
,
J.-P.
Wang
, and
C. H.
Kim
,
IEEE J. Solid-State Circuits
48
,
598
(
2013
).
3.
S.
Mangin
,
D.
Ravelosona
,
J. A.
Katine
,
M. J.
Carey
,
B. D.
Terris
, and
E. E.
Fullerton
,
Nat. Mater.
5
,
210
(
2006
).
4.
J. C.
Slonczewski
,
J. Magn. Magn. Mater.
159
,
L1
(
1996
).
5.
T.
Kawahara
,
K.
Ito
,
R.
Takemura
, and
H.
Ohno
,
Microelectron. Reliab.
52
,
613
(
2012
).
6.
G.-G.
An
,
J.-B.
Lee
,
S.-M.
Yang
,
J.-H.
Kim
,
W.-S.
Chung
, and
J.-P.
Hong
,
Acta Mater.
87
,
259
(
2015
).
7.
T.
Liu
,
Y.
Zhang
,
J. W.
Cai
, and
H. Y.
Pan
,
Sci. Rep.
4
,
5895
(
2014
).
8.
S.-E.
Lee
,
T.-H.
Shim
, and
J.-G.
Park
,
NPG Asia Mater.
8
,
e324
(
2016
).
9.
W.
Kim
,
J. H.
Jeong
,
Y.
Kim
,
W. C.
Lim
,
J. H.
Kim
,
J. H.
Park
,
H. J.
Shin
,
Y. S.
Park
,
K. S.
Kim
,
S. H.
Park
,
Y. J.
Lee
,
K. W.
Kim
,
H. J.
Kwon
,
H. L.
Park
,
H. S.
Ahn
,
S. C.
Oh
,
J. E.
Lee
,
S. O.
Park
,
S.
Choi
,
H. K.
Kang
, and
C.
Chung
, in
2011 International Electron Devices Meeting
(
IEEE
,
2011
), p.
24.1.1
10.
H.
Sato
,
M.
Yamanouchi
,
K.
Miura
,
S.
Ikeda
,
R.
Koizumi
,
F.
Matsukura
, and
H.
Ohno
,
IEEE Magn. Lett.
3
,
3000204
(
2012
).
11.
B.
Dieny
and
M.
Chshiev
,
Rev. Mod. Phys.
89
,
025008
(
2017
).
12.
T.
Liu
,
J. W.
Cai
, and
L.
Sun
,
AIP Adv.
2
,
032151
(
2012
).
13.
Y.-W.
Oh
,
K.-D.
Lee
,
J.-R.
Jeong
, and
B.-G.
Park
,
J. Appl. Phys.
115
,
17C724
(
2014
).
14.
N.
Miyakawa
,
D. C.
Worledge
, and
K.
Kita
,
IEEE Magn. Lett.
4
,
1000104
(
2013
).
15.
J.
Sinha
,
M.
Gruber
,
M.
Kodzuka
,
T.
Ohkubo
,
S.
Mitani
,
K.
Hono
, and
M.
Hayashi
,
J. Appl. Phys.
117
,
043913
(
2015
).
16.
M.
Akyol
,
J. Supercond. Novel Magn.
32
,
457
(
2019
).
17.
H.
Almasi
,
D. R.
Hickey
,
T.
Newhouse-Illige
,
M.
Xu
,
M. R.
Rosales
,
S.
Nahar
,
J. T.
Held
,
K. A.
Mkhoyan
, and
W. G.
Wang
,
Appl. Phys. Lett.
106
,
182406
(
2015
).
18.
J.
Dwivedi
,
M.
Gupta
,
V. R.
Reddy
,
A.
Mishra
,
V.
Srihari
,
K. K.
Pandey
, and
A.
Gupta
,
J. Magn. Magn. Mater.
466
,
311
(
2018
).
19.
M.
Yamanouchi
,
R.
Koizumi
,
S.
Ikeda
,
H.
Sato
,
K.
Mizunuma
,
K.
Miura
,
H. D.
Gan
,
F.
Matsukura
, and
H.
Ohno
,
J. Appl. Phys.
109
,
07C712
(
2011
).
20.
J.-H.
Kim
,
J.-B.
Lee
,
G.-G.
An
,
S.-M.
Yang
,
W.-S.
Chung
,
H.-S.
Park
, and
J.-P.
Hong
,
Sci. Rep.
5
,
16903
(
2015
).
21.
H.
Cheng
,
J.
Chen
,
S.
Peng
,
B.
Zhang
,
Z.
Wang
,
D.
Zhu
,
K.
Shi
,
S.
Eimer
,
X.
Wang
,
Z.
Guo
,
Y.
Xu
,
D.
Xiong
,
K.
Cao
, and
W.
Zhao
,
Adv. Electron. Mater.
6
,
2000271
(
2020
).
22.
H.
Sato
,
M.
Yamanouchi
,
S.
Ikeda
,
S.
Fukami
,
F.
Matsukura
, and
H.
Ohno
,
Appl. Phys. Lett.
101
,
022414
(
2012
).
23.
H.
Sato
,
M.
Yamanouchi
,
S.
Ikeda
,
S.
Fukami
,
F.
Matsukura
, and
H.
Ohno
,
IEEE Trans. Magn.
49
,
4437
(
2013
).
24.
H.
Sato
,
E. C. I.
Enobio
,
M.
Yamanouchi
,
S.
Ikeda
,
S.
Fukami
,
S.
Kanai
,
F.
Matsukura
, and
H.
Ohno
,
Appl. Phys. Lett.
105
,
062403
(
2014
).
25.
C.-F.
Pai
,
L.
Liu
,
Y.
Li
,
H. W.
Tseng
,
D. C.
Ralph
, and
R. A.
Buhrman
,
Appl. Phys. Lett.
101
,
122404
(
2012
).
26.
L.
Liu
,
C.-F.
Pai
,
Y.
Li
,
H. W.
Tseng
,
D. C.
Ralph
, and
R. A.
Buhrman
,
Science
336
,
555
(
2012
).
27.
R.
Ramaswamy
,
J. M.
Lee
,
K.
Cai
, and
H.
Yang
,
Appl. Phys. Rev.
5
,
031107
(
2018
).
28.
S.
Cho
,
S.-h. C.
Baek
,
K.-D.
Lee
,
Y.
Jo
, and
B.-G.
Park
,
Sci. Rep.
5
,
14668
(
2015
).
29.
C.
Zhang
,
S.
Fukami
,
K.
Watanabe
,
A.
Ohkawara
,
S.
DuttaGupta
,
H.
Sato
,
F.
Matsukura
, and
H.
Ohno
,
Appl. Phys. Lett.
109
,
192405
(
2016
).
30.
C.-F.
Pai
,
M.
Mann
,
A. J.
Tan
, and
G. S. D.
Beach
,
Phys. Rev. B
93
,
144409
(
2016
).
31.
M.-H.
Nguyen
and
C.-F.
Pai
,
APL Mater.
9
,
030902
(
2021
).
32.
C. Y.
Hu
and
C. F.
Pai
,
Adv. Quantum Technol.
3
,
2000024
(
2020
).
33.
L.
Liu
,
T.
Moriyama
,
D. C.
Ralph
, and
R. A.
Buhrman
,
Phys. Rev. Lett.
106
,
036601
(
2011
).
34.
C. F.
Pai
,
Y.
Ou
,
L. H. V.
Leao
,
D. C.
Ralph
, and
R. A.
Buhrman
,
Phys. Rev. B
92
,
064426
(
2015
).
35.
K.
Garello
,
I. M.
Miron
,
C. O.
Avci
,
F.
Freimuth
,
Y.
Mokrousov
,
S.
Blügel
,
S.
Auffret
,
O.
Boulle
,
G.
Gaudin
, and
P.
Gambardella
,
Nat. Nanotechnol.
8
,
587
(
2013
).
36.
C.-L.
Yang
and
C.-H.
Lai
,
Sci. Rep.
11
,
15214
(
2021
).
37.
F. R.
de Boer
,
W. C. M.
Mattens
,
R.
Boom
,
A. R.
Miedema
, and
A. K.
Niessen
,
Cohesion in Metals: Transition Metal Alloys
(
North-Holland
,
The Netherlands
,
1988
).
38.
J.
Liu
,
T.
Ohkubo
,
S.
Mitani
,
K.
Hono
, and
M.
Hayashi
,
Appl. Phys. Lett.
107
,
232408
(
2015
).
39.
J. W.
Lee
,
Y. W.
Oh
,
S. Y.
Park
,
A. I.
Figueroa
,
G. V. D.
Laan
,
G.
Go
,
K. J.
Lee
, and
B. G.
Park
,
Phys. Rev. B
96
,
064405
(
2017
).
40.
R.
Soucaille
,
M.
Belmeguenai
,
J.
Torrejon
,
J.-V.
Kim
,
T.
Devolder
,
Y.
Roussigné
,
S.-M.
Chérif
,
A. A.
Stashkevich
,
M.
Hayashi
, and
J.-P.
Adam
,
Phys. Rev. B
94
,
104431
(
2016
).