This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Unconventional Spin-Orbit Torques from Sputtered MoTe2 Films

Shuchen Li [email protected] Department of Materials Science and Engineering and Materials Research Laboratory, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, USA    Jonathan Gibbons Department of Materials Science and Engineering and Materials Research Laboratory, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, USA Department of Physics, University of California – San Diego, La Jolla, California 92093, USA    Stasiu Chyczewski Department of Electrical and Computer Engineering, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, USA    Zetai Liu Department of Electrical and Computer Engineering, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, USA    Hsu-Chih Ni Department of Materials Science and Engineering and Materials Research Laboratory, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, USA    Jiangchao Qian Department of Materials Science and Engineering and Materials Research Laboratory, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, USA    Jian-Min Zuo Department of Materials Science and Engineering and Materials Research Laboratory, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, USA    Jun-Fei Zheng Entegris Inc. Danbury, Connecticut 06810    Wenjuan Zhu Department of Electrical and Computer Engineering, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, USA    Axel Hoffmann [email protected] Department of Materials Science and Engineering and Materials Research Laboratory, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, USA
Abstract

Materials with strong spin-orbit coupling and low crystalline symmetry are promising for generating large unconventional spin-orbit torques (SOTs), such as in-plane field-like (FL) torques and out-of-plane damping-like (DL) torques, which can effectively manipulate and deterministcally switch an out-of-plane magnetization without the need for additional external in-plane magnetic fields. Here, we report SOTs generated by magnetron-sputtered 1T1T^{\prime} MoTe2/Permalloy (Py; Ni80Fe20)/MgO heterostructures using both spin-torque ferromagnetic resonance (ST-FMR) and second harmonic Hall measurements. We observed unconventional FL and DL torques in our samples due to spins polarized normal to the interface of MoTe2 and Py layers, and studied the influence of crystallographic order and MoTe2 layer thickness on the SOTs. By comparing the Raman spectra of 1T1T^{\prime} MoTe2 samples prepared in different ways, we found a tensile strain in sputtered MoTe2 films, which might further enhance the generation of unconventional torques by reducing the symmetry of 1T1T^{\prime} MoTe2.

I Introduction

Spin-orbit torques (SOT) are promising for novel information technologies, such as non-volatile magnetoresistive random-access memory [1, 2, 3] as they can efficiently manipulate magnetization dynamics by transferring angular momentum to the magnetic thin films [4] with demonstrated low power (130 pJ) and high speed (210 ps) [5]. When applying a charge current through the spin source material, which has large spin-orbit couplings, SOTs can be generated from spin currents and concomitant spin accumulations at material boundaries, e.g., through spin Hall effects (SHE) in bulk materials [6, 7, 8, 9, 10, 11, 12] and directly from electric current generated spin accumulations, e.g., from interfacial Rashba-Edelstein effects [13, 14, 15]. Such spin source materials can be heavy metals like Pt [8, 16] and Ta [9], topological insulators [15], and Weyl semimetals like WTe2 [17, 18]. Lately, there has been a growing interest in exotic in-plane FL torques, denoted as τFLz\vec{\tau}{{}_{FL}^{z}}, which are proportional to m^×z^\hat{m}\times\hat{z}, as well as out-of-plane DL torques τDLz\vec{\tau}{{}_{DL}^{z}} proportional to m^×(m^×z^)\hat{m}\times(\hat{m}\times\hat{z}), due to zz-polarized spins (where zz is the interface normal direction), which can deterministacally switch magnetizations with directions m^\hat{m} pointing out-of-plane without an external symmetry-breaking in-plane field. However, τFLz\vec{\tau}{{}_{FL}^{z}} and τDLz\vec{\tau}{{}_{DL}^{z}} are generally forbidden due to symmetry restrictions, but can be generated when the mirror symmetry in the plane perpendicular to the electric current direction is broken. Recent experiments have focused on exploring exotic SOTs generated from materials with intrinsically low symmetries, like WTe2/Permalloy (Py; Ni80Fe20) [17] and MnPd3/CoFeB [19], and from materials systems with artificial symmetry breakings, like strain-induced antisymmetry in NbSe2/Py devices [20]. Nevertheless, many questions still remain open with respect to the exact mechanisms behind the generation of those exotic torques.

In this work, we studied the SOTs generated from magnetron-sputtered MoTe2 films by using both spin-torque ferromagnetic resonance (ST-FMR) [8, 9] and second harmonic Hall measurements [21, 22]. We have observed sizable unconventional FL torques due to zz-polarized spins (τFLz\vec{\tau}{{}_{FL}^{z}}) and xx-polarized spins (τFLxm^×x^\vec{\tau}{{}_{FL}^{x}}\propto\hat{m}\times\hat{x}), and DL torques due to zz spins (τDLz\vec{\tau}{{}_{DL}^{z}}) in our MoTe2/Py/MgO devices. To investigate the origins of these torques, we studied their thickness and current direction dependencies. Through careful crystal structure characterizations, we found our sputtered MoTe2 has both preferential in- and out-of-plane alignments. We also found a strain in the MoTe2 films, which could further reduce the symmetry and contribute to the presence of the exotic torques.

II Sample fabrication and structural characterization

We synthesized MoTe2 films by magnetron co-sputtering elemental Te and Mo targets onto AA-plane sapphire (0.5-mm thick) substrates at temperatures ranging from 100–300C{}^{\circ}{C}, followed by 1 h of annealing ranging from 350–500C{}^{\circ}{C} in vacuum [23]. After the samples were cooled down, we deposited 10-nm Py and 2-nm MgO capping layer in-situ before bringing the samples to ambient conditions. We performed X-ray diffraction (XRD) and scanning transmission electron microscopy (STEM) to characterize the sputtered MoTe2 films. Fig. 1(a) shows the Raman measurements on the sputtered MoTe2 films and exfoliated 1T1T^{\prime} MoTe2 flakes (on AA-plane sapphire), indicating the sputtered MoTe2 to be 1T1T^{\prime} phase, a monoclinic structure with a single mirror plane along (perpendicular to) its a(b) axis and the c axis to be out-of-plane [shown in Fig. 1(b)].

However, interestingly, for all MoTe2 films (7 nm, 15 nm, and 40 nm) grown by magnetron sputtering, the characteristic Raman peak around 168 cm-1 shifted to lower energies compared to the same peak from the exfoliated 1T1T^{\prime} MoTe2 flake. Such a red shift of the Raman peak is generally caused by a strain from a tensile stress, which will lead to the elongation of the lattice and a decrease in the bond strength between neighboring atoms [24, 25, 26]. We think an in-plane strain with a perpendicular component to the mirror plane of 1T1T^{\prime} MoTe2 was induced during high-temperature processing, such that a red shift of 168 cm-1 Raman mode vibrating mostly in-plane perpendicular to the mirror plane of 1T1T^{\prime} MoTe2 [27] was observed across the samples. Fig. 1(d) shows the STEM image of our sputtered 15-nm MoTe2 sample. We calculated the distance between two adjacent bright fringes to be 1.39 nm-1 by performing the fast Fourier transform of the blue circled area, which corresponds to 0.719 nm in real space and is close to the distance between the two MoTe2 layers held by the van der Waals forces (0.693 nm), indicating the cc-axis of MoTe2 is mostly aligned vertically. We further explored the in-plane orientation of MoTe2 films to investigate if there is any preferential alignment with the AA-plane sapphire substrate by using polarized Raman spectroscopy. We illuminated sputtered MoTe2/sapphire samples with linearly polarized light at different ϕR\phi_{R}, where ϕR\phi_{R} is the angle between the polarization direction and the [11¯001\bar{1}00] direction of the AA-plane sapphire substrate. As shown in Fig. 1 (c), the polarization-dependent Raman peaks and intensities (green arrows) demonstrate there is a crystallographic texture within the sputtered MoTe2, with the a and b axis of the MoTe2 being preferentially aligned along the [00010001] and [11¯001\bar{1}00] directions of the AA-plane sapphire, respectively [28, 5].

Refer to caption
Figure 1: (a) Raman spectra of mangetron-sputtered 40-nm (dark blue), 15-nm (blue), 7-nm (light blue), and exfoliated (grey) MoTe2 samples. The spectra are shifted by an offset of 10 with respect to each other. The green dashed line and the red arrows indicate the theoretical and the measured MoTe2 Raman shift position. (b) Crystal structure of 1T1T^{\prime} MoTe2 with the only mirror plane (red dashed line) along the a-axis. On the left shows the AA-plane sapphire substrate orientations. (c) Polarized Raman spectra for 15-nm MoTe2/sapphire. A linearly polarized 633-nm light illuminates the sample with the polarization angle ϕR\phi_{R} with respect to the [11¯001\bar{1}00] direction of the AA-plane sapphire substrate [shown in (b)]. ϕR\phi_{R} = 0 means the polarization direction is parallel to the [11¯001\bar{1}00] direction. The green arrows indicate different Raman modes and different peak intensities. (d) The STEM image of the sputtered MoTe2 film and the fast Fourier transform of the blue circled area.

III Spintransport measurements

III.1 Spin-Torque Ferromagnetic Resonance

To study the characteristics of spin-orbit torques generated from sputtered MoTe2 layer, we performed spin-torque ferromagnetic resonance (ST-FMR) measurements [8] on our MoTe2(15)/Py(10)/MgO(2) devices (the individual layer thicknesses in parenthesis are in nm). As indicated in Figs. 2(a) and (b), the samples were patterned and integrated into ground-signal-ground coplanar waveguides by ion-milling and lithography, and a radio-frequency current IrfI_{rf} at 5–9 GHz along the [11¯001\bar{1}00] direction was applied through our devices with an external magnetic field HextH_{ext} sweeping in the xyxy plane from 0.1 T to 0 T and from -0.1 T to 0 T. We then measured the resultant homodyne dc mixing voltages, VmixV_{mix}, due to the coupling of the rf current and the anisotropic magnetoresistance (AMR) of Py modulated by SOTs from the MoTe2 layer through a lockin amplifer. Fig. 2(c) shows the measured mixing voltages VmixV_{mix} of device 1 for positive and negative field scans at ϕH\phi_{H} = 45, where ϕH\phi_{H} is the angle between the external magnetic field and the applied IrfI_{rf} current. VmixV_{mix} can be fitted by the sum of symmetric and antisymmetric Lorentzian functions

Vmix\displaystyle V_{mix} =\displaystyle= SΔ2(HextH0)2+Δ2\displaystyle S\dfrac{{\Delta}^{2}}{(H_{ext}-H_{0})^{2}+{\Delta}^{2}} (1)
+AΔ(HextH0)(HextH0)2+Δ2,\displaystyle+A\dfrac{{\Delta}(H_{ext}-H_{0})}{(H_{ext}-H_{0})^{2}+{\Delta}^{2}}~{},

where HextH_{ext} is the applied field, H0H_{0} is the ferromagnetic resonance field of permalloy, Δ{\Delta} is the half width at half maximum resonance linewidth, and SS and AA are the fitting parameters representing the sizes of the symmetric and antisymmetric Lorentzians VSV_{S} and VAV_{A}, which correspond to the sizes of the in-plane τ\tau_{\parallel} and out-of-plane τ\tau_{\bot} torques respectively.

As shown in Fig. 2(c), the different line shapes of symmetric (VSV_{S}, blue) and antisymmetric (VAV_{A}, red) Lorentzians for positive and negative field scans indicate the presence of unconventional SOTs due to zz-polarized spins. We added the mixing voltages measured from positive field and negative field scans, which cancels voltages contributed by spins polarized in-plane, while voltages due to spins polarized out-of-plane (along zz) will add up constructively, due to their distinct symmetry characteristics in relation to the field direction (Eqs. 2 and 3). We then divided the added voltages by two and obtained the voltage (Vmix_zV_{mix\_z}), due to torques from zz-polarized spins [Fig. 2(d)]. Through fitting Vmix_zV_{mix\_z} using Eq. 1, we confirmed the existence of both τDLz\vec{\tau}{{}_{DL}^{z}} and τFLz\vec{\tau}{{}_{FL}^{z}} within device 1 of our MoTe2(15)/Py(10)/MgO(2).

Refer to caption
Figure 2: (a) Diagram of our measurement setup for ST-FMR. A signal generator injects a GHz rf current whose amplitude is modulated by the reference signal of a lock-in amplifier into the device through the rf port of a bias tee. The mixing dc voltage is measured by the lock-in amplifier through the dc port of the bias tee. The dimension of the device is 80–130 μ\mum in length and 20–40 μ\mum in width. (b) A schematic of the spin-torque ferromagnetic resonance measurements on MoTe2/Py/MgO devices. (c) The measured dc mixing voltages of sample 1 device 1 of MoTe2(15)/Py/MgO at ϕH\phi_{H} = 45 for positive and negative field scans. The power and frequency of the current is 4 dBm and 6 GHz, and the current directin is along [11¯001\bar{1}00]. The fit for the mixing voltage is the green curve, which is the sum of VSV_{S} (blue) and VAV_{A} (red). (d) The mixing voltages Vmix_zV\textsubscript{mix\_z} with contributions solely from zz-polarized spins, and we found Sz=0.313S_{z}=-0.313 and Az=0.140A_{z}=0.140, which are proportional to the sizes of τFLz\vec{\tau}{{}_{FL}^{z}} and τDLz\vec{\tau}{{}_{DL}^{z}}.

To extract different components of SOTs from τ\tau_{\parallel} and τ\tau_{\bot} that are related to VSV_{S} and VAV_{A}, we varied the angle ϕH\phi_{H} between the current and field, and plot the extracted SS and AA as a function of angle ϕH\phi_{H}. SOTs generated by xx, yy and zz-polarized spins have distinct angular dependencies for both SS and AA described by the following equations,

S\displaystyle{S} =\displaystyle= SsDLyin(2ϕH)cos(ϕH)+SsFLzin(2ϕH)\displaystyle S{{}_{DL}^{y}}sin(2\phi_{H})cos(\phi_{H})+S{{}_{FL}^{z}}sin(2\phi_{H}) (2)
+SsDLxin(2ϕH)sin(ϕH)\displaystyle+S{{}_{DL}^{x}}sin(2\phi_{H})sin(\phi_{H})
A\displaystyle{A} =\displaystyle= AsFLyin(2ϕH)cos(ϕH)+AsDLzin(2ϕH)\displaystyle A{{}_{FL}^{y}}sin(2\phi_{H})cos(\phi_{H})+A{{}_{DL}^{z}}sin(2\phi_{H}) (3)
+AsFLxin(2ϕH)sin(ϕH),\displaystyle+A{{}_{FL}^{x}}sin(2\phi_{H})sin(\phi_{H}),

and we can obtain the sizes of τDLy\vec{\tau}{{}_{DL}^{y}}, τFLz\vec{\tau}{{}_{FL}^{z}}, τDLxm^×(m^×x^)\vec{\tau}{{}_{DL}^{x}}\propto\hat{m}\times(\hat{m}\times\hat{x}) ,τFLy\vec{\tau}{{}_{FL}^{y}}, τDLz\vec{\tau}{{}_{DL}^{z}}, τFLx\vec{\tau}{{}_{FL}^{x}} by fitting their angular dependencies using Eqs. 2 and 3, which are proportional to the fit values of SDLyS{{}_{DL}^{y}}, SFLzS{{}_{FL}^{z}}, SDLxS{{}_{DL}^{x}}, AFLyA{{}_{FL}^{y}}, ADLzA{{}_{DL}^{z}}, and AFLxA{{}_{FL}^{x}}, respectively.

As shown in Figs. 3(a) and (b), in device 1 of our MoTe2(15)/Py/MgO sample, we found AFLyA{{}_{FL}^{y}} = 0.624, which we assume is mainly contributed by the Oersted field generated by the rf currents, and SDLyS{{}_{DL}^{y}} = 0.155, related to conventional τDLy\vec{\tau}{{}_{DL}^{y}} in our sample. Also, we found SFLzS{{}_{FL}^{z}} = -0.240 and ADLzA{{}_{DL}^{z}} = 0.104, indicating sizeable unconventional τFLz\vec{\tau}{{}_{FL}^{z}} and τDLz\vec{\tau}{{}_{DL}^{z}} due to zz-polarized spins. In addition, we noticed that the polarity of τFLz\vec{\tau}{{}_{FL}^{z}} is always opposite to that of τDLz\vec{\tau}{{}_{DL}^{z}} for all the devices measured from MoTe2(15)/Py/MgO. This indicates that the mechanisms behind the generation of the two torques are the same or are strongly correlated. The exact mechanism still remains unclear but this phenomenon has also been reported in the exfoliated single crystalline 1T1T^{\prime} MoTe2/Py samples [5].

Refer to caption
Figure 3: (a), (c) and (e) antisymmetric components VAV_{A} as a function of angle ϕH\phi_{H} for ϕI\phi_{I} = 0, 30 and 90. (b), (d) and (f) Symmetric components VSV_{S} as a function of angle ϕH\phi_{H} for ϕI\phi_{I} = 0, 30 and 90. The red and blue dots are extracted from the measured VmixV_{mix}, and the black lines are the fitted curves using Eqs. 2 and 3.

III.2 Current Direction Dependence

We studied the effects of crystallographic order on different torque components by applying the current at different angles ϕI\phi_{I} through our MoTe2(15)/Py devices, where ϕI\phi_{I} is the angle between the current direction II and [11¯001\bar{1}00] of the AA-plane substrate, with ϕI\phi_{I} = 0 being II parallel to [11¯001\bar{1}00]. As shown in Figs. 3(c) and (d), we found sizeable τFLx\vec{\tau}{{}_{FL}^{x}} due to xx-polarized spins in the antisymmetric component when ϕI\phi_{I} = 30, but no τDLz\vec{\tau}{{}_{DL}^{z}} from the antisymmetric component. Moreover, we found the sizes of all unconventional torques diminished at ϕI\phi_{I} = 90, and only the conventional torques τDLy\vec{\tau}{{}_{DL}^{y}} remained [Figs.3(e) and (f)].

The spin-orbit torque (SOT) efficiency quantifies the ability of the spin-orbit coupling to convert an applied electric current into a torque that influences the magnetization dynamics, and we calculated the SOT efficiencies ξDLy\xi{{}_{DL}^{y}}, ξFLz\xi{{}_{FL}^{z}} and ξDLz\xi{{}_{DL}^{z}} using the equation

ξ=(S,DLyS,FLzorA)DLzAFLyeμ0MStd[1+Meff/H0]1/2,\xi=\dfrac{(S{{}_{DL}^{y}},{S{{}_{FL}^{z}}},or{A{{}_{DL}^{z}}})}{A{{}_{FL}^{y}}}\dfrac{e\mu_{0}M_{S}td}{\hbar}[1+M_{eff}/H_{0}]^{1/2}, (4)

where tt and dd are the thicknesses of Py and MoTe2 layers, MSM_{S} is the saturation magnetization of Py, which is approximately μ0MS0.8\mu_{0}M_{S}\approx 0.8 T. Figs. 4(a)–(c) show the calculated SOT efficiencies measured across different devices from MoTe2(15)/Py/MgO heterostructures at different ϕI\phi_{I}.

Refer to caption
Figure 4: (a) ξDLy\xi{{}_{DL}^{y}} at different ϕI\phi_{I}, (b) and (c) are absolute values of ξFLz\xi{{}_{FL}^{z}} and ξDLz\xi{{}_{DL}^{z}} for better study the trend with respect to ϕI\phi_{I}. (d), (e) and (f) are ξDLy\xi{{}_{DL}^{y}}, |ξ|FLz|\xi{{}_{FL}^{z}}| and |ξ|DLz|\xi{{}_{DL}^{z}}| for devices with different MoTe2 thicknesses (7, 15 and 40 nm).

We found ξDLy\xi{{}_{DL}^{y}} did not display a clear dependence on ϕI\phi_{I}, whereas ξFLz\xi{{}_{FL}^{z}} and ξDLz\xi{{}_{DL}^{z}} both showed strong relations to current directions. When the current was applied at an angle to the [11¯001\bar{1}00] direction, ξDLz\xi{{}_{DL}^{z}} vanished, while ξFLz\xi{{}_{FL}^{z}} decreased at ϕI\phi_{I} = 30 and diminished further at ϕI\phi_{I} = 60 and 90. For ξDLz\xi{{}_{DL}^{z}} and ξFLz\xi{{}_{FL}^{z}}, the strong ϕI\phi_{I} dependence indicates symmetry-related origin for the generation of τDLz\vec{\tau}{{}_{DL}^{z}} and τFLz\vec{\tau}{{}_{FL}^{z}}, which would emerge when a current is applied perpendicular to the mirror plane of 1T1T^{\prime} MoTe2, and become symmetry-forbidden when the current is applied along the mirror plane.

III.3 Thickness Dependence

To better understand origins of different spin-orbit torque components, we investigated the thickness dependencies of different torques. Similarly, we performed ST-FMR measurements on MoTe2(40 or 7)/Py/MgO heterostructures and calculated ξDLy\xi{{}_{DL}^{y}}, ξDLz\xi{{}_{DL}^{z}} and ξFLz\xi{{}_{FL}^{z}} using Eq. 4 for all the devices measured. As shown in Figs. 4(d)–(f), ξDLy\xi{{}_{DL}^{y}} and ξFLz\xi{{}_{FL}^{z}} increased when the MoTe2 thickness dd increased from 7 nm to 15 nm, and became saturated when dd = 40 nm, indicating possible origns to be bulk effects. The ϕI\phi_{I} independence of ξDLy\xi{{}_{DL}^{y}} suggests that τDLy\vec{\tau}{{}_{DL}^{y}} originates from the spin Hall effect (SHE). For ξτDLz\xi{\tau}{{}_{DL}^{z}}, we did not observe any clear trend with respect to dd [Fig. 4(f)], which would imply the origin of τDLz\vec{\tau}{{}_{DL}^{z}} to be related to the interfacial properties between MoTe2 and Py. Finally, for MoTe2(40) and MoTe2(7), we found the polarities of ξFLz\xi{{}_{FL}^{z}} and ξDLz\xi{{}_{DL}^{z}} were always opposite for all the devices measured, just as what has been found in MoTe2(15).

III.4 Second Harmonic Hall Measurements

Second harmonic Hall measurements are useful for analyzing in-plane τ\tau_{\parallel} and out-of-plane τ\tau_{\bot} spin-orbit torques, and can be complimentary to ST-FMR. For our MoTe2(15)/Py/MgO sample, we applied a low frequency (ω\omega = 1131 kHz) ac current with 4 mA amplitude along the xx direction, with the magnetization of Py defined in-plane by an external magnetic field (HextH_{ext}), and measured the transverse second harmonic Hall voltage V2ωV^{2\omega} along yy. τ\tau_{\parallel}(τ\tau_{\bot}) generated by MoTe2 has an effective out-of-plane(in-plane) field that rotates the magnetization of Py out-of-plane(in-plane) such that it modulates the change of the anomalous Hall resistance RAHER_{AHE} (planar Hall resistance RPHER_{PHE}) at a frequency ω\omega of the ac current. The change of RAHER_{AHE} and RPHER_{PHE} couples to the applied current and gives rise to Hall voltages V2ωV^{2\omega} on the second harmonic 2ω2\omega that can be detected through a lock-in amplifier. Figure 5(a) shows V2ωV^{2\omega} as a function of ϕH\phi_{H}, where ϕH\phi_{H} is the angle between HextH_{ext} and the applied current, at various fields. For an in-plane magnetization system, different spin-orbit torque components can be extracted by their distinct angular dependencies described in the following equation [22, 29, 17]

V2ω\displaystyle V^{2\omega} =\displaystyle= VPHEcos(2ϕH)(H+FLyHOeHextcos(ϕH)+HFLxHextsin(ϕH)+HDLzHext)\displaystyle V_{PHE}cos(2\phi_{H})(\dfrac{H{{}_{FL}^{y}}+H_{Oe}}{H_{ext}}cos(\phi_{H})+\dfrac{H{{}_{FL}^{x}}}{H_{ext}}sin(\phi_{H})+\dfrac{H{{}_{DL}^{z}}}{H_{ext}}) (5)
+12VAHE(HDLyHext+Hkcos(ϕH)+HDLxHext+Hksin(ϕH)+HFLzHext+Hk)+C,\displaystyle+\dfrac{1}{2}V_{AHE}(\dfrac{H{{}_{DL}^{y}}}{H_{ext}+H_{k}}cos(\phi_{H})+\dfrac{H{{}_{DL}^{x}}}{H_{ext}+H_{k}}sin(\phi_{H})+\dfrac{H{{}_{FL}^{z}}}{H_{ext}+H_{k}})+C~{},

where HFLyH{{}_{FL}^{y}}, HFLxH{{}_{FL}^{x}}, HDLzH{{}_{DL}^{z}}, HDLyH{{}_{DL}^{y}}, HDLxH{{}_{DL}^{x}} and HFLzH{{}_{FL}^{z}} are the spin-orbit fields corresponding to the respective torques, and CC is an offset constant. Through fitting, we observed sizable contributions from τFLy\vec{\tau}{{}_{FL}^{y}} + τOe\vec{\tau}_{Oe}, τDLy\vec{\tau}{{}_{DL}^{y}} and τDLz\vec{\tau}{{}_{DL}^{z}} to V2ωV^{2\omega}, and plot the extracted components VDL,z2ωV^{2\omega}_{DL,z} and VFL,y+Oe2ωV^{2\omega}_{FL,y+Oe} as a function of 1/Hext1/H_{ext} [see Figs. 5(c) and (d)]. Figure 5(b) shows the angular dependencies of the transverse VωV^{\omega}, from which we derived the planar Hall effect (VPHEV_{PHE}) to be 0.45 mV. Thus, we calculated H+FLyHOeH{{}_{FL}^{y}}+H_{Oe} and HDLzH{{}_{DL}^{z}} to be 8.0 A/m and 1.67 A/m. Finally, we calculated the SOT efficiency ξDLz\xi^{z}_{DL} = 0.068 using

ξDLz=2eMstFMHDLzJMoTe2,\xi^{z}_{DL}=\dfrac{2eM_{s}t_{FM}}{\hbar}\dfrac{H^{z}_{DL}}{J_{MoTe_{2}}}, (6)

where ee, tt and \hbar are the electron charge, thickness of permalloy and reduced Planck’s constant. JMoTe2=IMoTe2/(dw)J_{MoTe_{2}}=I_{MoTe_{2}}/(dw) is the current density flowing through the MoTe2 layer, where IMoTe2I_{MoTe_{2}} is estimated to be 0.27 mA (see Supplemental Material S1), dd = 15 nm and ww = 20 μm\mu{m} is the width of the device. We also estimated the Oersted field generated by IMoTe2I_{MoTe_{2}} to be 6.75 A/m using Ampere’s law HOe=IMoTe2/2wH_{Oe}=I_{MoTe_{2}}/2w assuming the sample to be an infinitely wide plate. This confirms that the Oersted field dominates VFL,y+Oe2ωV^{2\omega}_{FL,y+Oe}.

Refer to caption
Figure 5: (a) V2ωV^{2\omega} (dots) of MoTe2(15)/Py/MgO as a function of ϕH\phi_{H} for various fields and the fit curves (lines) using Eq. 5. (b) VωV^{\omega} as a function of ϕH\phi_{H} under HextH_{ext} = 0.22 T. (c) and (d) Components of V2ωV^{2\omega}, VDL,z2ωV^{2\omega}_{DL,z} and VFL,y+Oe2ωV^{2\omega}_{FL,y+Oe}, contributed by HDLzH{{}_{DL}^{z}} and H+FLyHOeH{{}_{FL}^{y}}+H_{Oe}, with linear fit red lines.

IV Discussion

The observations of unconventional spin-orbit torques in magnetron-sputtered 1T1T^{\prime} MoTe2 films are possibly caused by two effects. First, as it can seen from Figs. 1(c) and (d), the sputtered MoTe2 films have a textured structure with the c-axis mostly aligned vertically and a(b) axis preferentially aligned along the [00010001]([11¯001\bar{1}00]) direction of the AA-plane sapphire substrate, such that it preserves some macroscopic low symmetries that lack inversion symmetry and has only one mirror plane, enabling the generation of symmetry-forbidden SOTs when a current is applied at an angle to the mirror plane breaking this symmetry.

Second, a strain identified from the red shift of the 168 cm-1 Raman mode for all sputtered MoTe2 films can further reduce the symmetry and enhance the effects generating unconventional spin-orbit torques. Such a strain can be caused by a tensile stress introduced by the mismatch between the thermal expansion coefficients of the substrates [30, 31] and MoTe2 during high temperature growth. It is possible that the large SOT efficiency range in 1T1T^{\prime} MoTe2 films is caused by the relative alignment of the strain axis with the local grain orientations, analogous to the strain-induced τFLz\vec{\tau}{{}_{FL}^{z}} found in single crystalline NbSe2 [20], and τDLz\vec{\tau}{{}_{DL}^{z}} in exfoliated single crystalline 1T1T^{\prime} MoTe2 [5].

Finally, we performed similar measurements on 1T1T^{\prime} MoTe2 grown by metal-organic chemical vapor deposition[32], and observed different results [see S4 in Supplemental Materials].

V Conclusion

In conclusion, we have studied and reported spin-torques found in magnetron-sputtered MoTe2/Py devices by performing ST-FMR and second harmonic Hall measurements on MoTe2/Py(10)/MgO(2) samples. We explored the origins of different SOTs by studying the influence of the crystallographic order and the MoTe2 thickness on the SOT efficiencies. Based on the ϕI\phi_{I} and thickness dd dependence of SOT efficiencies of different torques, the origins for τDLy\vec{\tau}{{}_{DL}^{y}} and τFLz\vec{\tau}{{}_{FL}^{z}} are probably due to bulk spin-orbit coupling, specifically the spin Hall effect for τDLy\vec{\tau}{{}_{DL}^{y}}. For τDLz\vec{\tau}{{}_{DL}^{z}}, the origins may be related to interfacial properties between the sputtered MoTe2 films and the magnetic layer. Furthermore, the origins of τFLz\vec{\tau}{{}_{FL}^{z}} and τDLz\vec{\tau}{{}_{DL}^{z}} are closely related as their polarities are always opposite to each other. Our findings show that strains can be promising for generating large exotic in-plane FL and out-of-plane DL torques by introducing additional symmetry reduction, and single crystallinity are not necessary in generating exotic torques.

VI Acknowledgement

This work was supported by the National Science Foundation under Grant No. ECCS-2031870. In addition data analysis and discussion by J.G. was supported as part of Quantum Materials for Energy Efficient Neuromorphic Computing (Q-MEEN-C), an Energy Frontier Research Center funded by the US Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES), under Award No. DE-SC0019273. S.C., Z.L., and W.Z. would like to acknowledge the support from Entegris Incorporated under grant No. Entegris 108252.

References

  • Brataas et al. [2012] A. Brataas, A. D. Kent, and H. Ohno, Current-induced torques in magnetic materials, Nat. Mater. 11, 372 (2012).
  • Natsui et al. [2021] M. Natsui, A. Tamakoshi, H. Honjo, T. Watanabe, T. Nasuno, C. Zhang, T. Tanigawa, H. Inoue, M. Niwa, T. Yoshiduka, Y. Noguchi, M. Yasuhira, Y. Ma, H. Shen, S. Fukami, H. Sato, S. Ikeda, H. Ohno, T. Endoh, and T. Hanyu, Dual-port SOT-MRAM achieving 90-MHz read and 60-MHz write operations under field-assistance-free condition, IEEE Journal of Solid-State Circuits 56, 1116 (2021).
  • Shao et al. [2021] Q. Shao, P. Li, L. Liu, H. Yang, S. Fukami, A. Razavi, H. Wu, K. Wang, F. Freimuth, Y. Mokrousov, M. D. Stiles, S. Emori, A. Hoffmann, J. Åkerman, K. Roy, J.-P. Wang, S.-H. Yang, K. Garello, and W. Zhang, Roadmap of spin–orbit torques, IEEE Trans. Magn, 57, 1 (2021).
  • Manchon et al. [2019] A. Manchon, J. Železný, I. M. Miron, T. Jungwirth, J. Sinova, A. Thiaville, K. Garello, and P. Gambardella, Current-induced spin-orbit torques in ferromagnetic and antiferromagnetic systems, Rev. Mod. Phys. 91, 035004 (2019).
  • Stiehl et al. [2019] G. M. Stiehl, R. Li, V. Gupta, I. E. Baggari, S. Jiang, H. Xie, L. F. Kourkoutis, K. F. Mak, J. Shan, R. A. Buhrman, and D. C. Ralph, Layer-dependent spin-orbit torques generated by the centrosymmetric transition metal dichalcogenide βMoTe2\beta-{\mathrm{MoTe}}_{2}Phys. Rev. B 100, 184402 (2019).
  • Zhang [2000] S. Zhang, Spin Hall effect in the presence of spin diffusion, Phys. Rev. Lett. 85, 393 (2000).
  • Sinova et al. [2015] J. Sinova, S. O. Valenzuela, J. Wunderlich, C. H. Back, and T. Jungwirth, Spin Hall effects, Rev. Mod. Phys. 87, 1213 (2015).
  • Liu et al. [2011] L. Liu, T. Moriyama, D. C. Ralph, and R. A. Buhrman, Spin-torque ferromagnetic resonance induced by the spin Hall effect, Phys. Rev. Lett. 106, 036601 (2011).
  • Liu et al. [2012a] L. Liu, C.-F. Pai, Y. Li, H. W. Tseng, D. C. Ralph, and R. A. Buhrman, Spin-torque switching with the giant spin Hall effect of tantalum, Science 336, 555 (2012a).
  • Liu et al. [2012b] L. Liu, O. J. Lee, T. J. Gudmundsen, D. C. Ralph, and R. A. Buhrman, Current-induced switching of perpendicularly magnetized magnetic layers using spin torque from the spin Hall effect, Phys. Rev. Lett. 109, 096602 (2012b).
  • Hoffmann [2013] A. Hoffmann, Spin Hall effects in metals, IEEE Trans. Magn, 49, 5172 (2013).
  • Zhang et al. [2016] W. Zhang, J. Sklenar, B. Hsu, W. Jiang, M. B. Jungfleisch, J. Xiao, F. Y. Fradin, Y. Liu, J. E. Pearson, J. B. Ketterson, Z. Yang, and A. Hoffmann, Research Update: Spin transfer torques in permalloy on monolayer MoS2, APL Mater. 4, 032302 (2016).
  • Shao et al. [2016] Q. Shao, G. Yu, Y.-W. Lan, Y. Shi, M.-Y. Li, C. Zheng, X. Zhu, L.-J. Li, P. K. Amiri, and K. L. Wang, Strong Rashba-Edelstein effect-induced spin–orbit torques in monolayer transition metal dichalcogenide/ferromagnet bilayers, Nano Lett. 16, 7514 (2016).
  • Emori et al. [2016] S. Emori, T. Nan, A. M. Belkessam, X. Wang, A. D. Matyushov, C. J. Babroski, Y. Gao, H. Lin, and N. X. Sun, Interfacial spin-orbit torque without bulk spin-orbit coupling, Phys. Rev. B 93, 180402 (2016).
  • Mellnik et al. [2014] A. R. Mellnik, J. S. Lee, A. Richardella, J. L. Grab, P. J. Mintun, M. H. Fischer, A. Vaezi, A. Manchon, E.-A. Kim, N. Samarth, and D. C. Ralph, Spin-transfer torque generated by a topological insulator, Nature 511, 449 (2014).
  • Mihai Miron et al. [2010] I. Mihai Miron, G. Gaudin, S. Auffret, B. Rodmacq, A. Schuhl, S. Pizzini, J. Vogel, and P. Gambardella, Current-driven spin torque induced by the Rashba effect in a ferromagnetic metal layer, Nat. Mater. 9, 230 (2010).
  • MacNeill et al. [2017] D. MacNeill, G. M. Stiehl, M. H. D. Guimaraes, R. A. Buhrman, J. Park, and D. C. Ralph, Control of spin–orbit torques through crystal symmetry in WTe2/ferromagnet bilayers, Nat. Phys. 13, 300 (2017).
  • Li et al. [2018] P. Li, W. Wu, Y. Wen, C. Zhang, J. Zhang, S. Zhang, Z. Yu, S. A. Yang, A. Manchon, and X.-x. Zhang, Spin-momentum locking and spin-orbit torques in magnetic nano-heterojunctions composed of Weyl semimetal WTe2Nat. Commun. 9, 3990 (2018).
  • Dc et al. [2023] M. Dc, D.-F. Shao, V. D.-H. Hou, A. Vailionis, P. Quarterman, A. Habiboglu, M. B. Venuti, F. Xue, Y.-L. Huang, C.-M. Lee, M. Miura, B. Kirby, C. Bi, X. Li, Y. Deng, S.-J. Lin, W. Tsai, S. Eley, W.-G. Wang, J. A. Borchers, E. Y. Tsymbal, and S. X. Wang, Observation of anti-damping spin–orbit torques generated by in-plane and out-of-plane spin polarizations in MnPd3Nat. Mater. , 1 (2023).
  • Guimarães et al. [2018] M. H. D. Guimarães, G. M. Stiehl, D. MacNeill, N. D. Reynolds, and D. C. Ralph, Spin–orbit torques in NbSe2/permalloy bilayers, Nano Lett. 18, 1311 (2018).
  • Avci et al. [2014] C. O. Avci, K. Garello, M. Gabureac, A. Ghosh, A. Fuhrer, S. F. Alvarado, and P. Gambardella, Interplay of spin-orbit torque and thermoelectric effects in ferromagnet/normal-metal bilayers, Phys. Rev. B 90, 224427 (2014).
  • Hayashi et al. [2014] M. Hayashi, J. Kim, M. Yamanouchi, and H. Ohno, Quantitative characterization of the spin-orbit torque using harmonic Hall voltage measurements, Phys. Rev. B 89, 144425 (2014).
  • Pace et al. [2021] S. Pace, L. Martini, D. Convertino, D. H. Keum, S. Forti, S. Pezzini, F. Fabbri, V. Mišeikis, and C. Coletti, Synthesis of large-scale monolayer 1T’-MoTe2 and its stabilization via scalable hBN encapsulation, ACS Nano 15, 4213 (2021).
  • Karki et al. [2020] B. Karki, B. Freelon, M. Rajapakse, R. Musa, S. M. S. Riyadh, B. Morris, U. Abu, M. Yu, G. Sumanasekera, and J. B. Jasinski, Strain-induced vibrational properties of few layer black phosphorus and MoTe2 via Raman spectroscopy, Nanotechn. 31, 425707 (2020).
  • Imajo et al. [2021] T. Imajo, T. Suemasu, and K. Toko, Strain effects on polycrystalline germanium thin films, Sci. Rep. 11, 8333 (2021).
  • Angel et al. [2019] R. J. Angel, M. Murri, B. Mihailova, and M. Alvaro, Stress, strain and Raman shifts, Z. Kristallogr. Cryst. Mater. 234, 129 (2019).
  • Ma et al. [2016] X. Ma, P. Guo, C. Yi, Q. Yu, A. Zhang, J. Ji, Y. Tian, F. Jin, Y. Wang, K. Liu, T. Xia, Y. Shi, and Q. Zhang, Raman scattering in the transition-metal dichalcogenides of 1TMoTe2,TdMoTe21{T}^{{}^{\prime}}\text{$-$}\mathrm{MoT}{\mathrm{e}}_{2},{T}_{d}\text{$-$}\mathrm{MoT}{\mathrm{e}}_{2}, and TdWTe2{T}_{d}\text{$-$}\mathrm{WT}{\mathrm{e}}_{2}Phys. Rev. B 94, 214105 (2016).
  • Beams et al. [2016] R. Beams, L. G. Cançado, S. Krylyuk, I. Kalish, B. Kalanyan, A. K. Singh, K. Choudhary, A. Bruma, P. M. Vora, F. Tavazza, A. V. Davydov, and S. J. Stranick, Characterization of few-layer 1T’ MoTe2 by polarization-resolved second harmonic generation and Raman scattering, ACS Nano 10, 9626 (2016).
  • Mosendz et al. [2010] O. Mosendz, J. E. Pearson, F. Y. Fradin, G. E. W. Bauer, S. D. Bader, and A. Hoffmann, Quantifying spin hall angles from spin pumping: Experiments and theory, Phys. Rev. Lett. 104, 046601 (2010).
  • Tian et al. [2016] S. Tian, Y. Yang, Z. Liu, C. Wang, R. Pan, C. Gu, and J. Li, Temperature-dependent Raman investigation on suspended graphene: Contribution from thermal expansion coefficient mismatch between graphene and substrate, Carbon 104, 27 (2016).
  • Yoon et al. [2011] D. Yoon, Y.-W. Son, and H. Cheong, Negative thermal expansion coefficient of graphene measured by Raman spectroscopy, Nano Lett. 11, 3227 (2011).
  • Chyczewski et al. [2024] S. T. Chyczewski, H. Lee, S. Li, M. Eladl, J.-F. Zheng, A. Hoffmann, and W. Zhu, Strong damping-like torques in wafer-scale mote2 grown by mocvd (2024), arXiv:2404.19135 .