Abstract
Controlling magnetization dynamics is imperative for developing ultrafast spintronics and tunable microwave devices. However, the previous research has demonstrated limited electric-field modulation of the effective magnetic damping, a parameter that governs the magnetization dynamics. Here, we propose an approach to manipulate the damping by using the large damping enhancement induced by the two-magnon scattering and a nonlocal spin relaxation process in which spin currents are resonantly transported from antiferromagnetic domains to ferromagnetic matrix in a mixed-phased metallic alloy FeRh. This damping enhancement in FeRh is sensitive to its fraction of antiferromagnetic and ferromagnetic phases, which can be dynamically tuned by electric fields through a strain-mediated magnetoelectric coupling. In a heterostructure of FeRh and piezoelectric PMN-PT, we demonstrated a more than 120% modulation of the effective damping by electric fields during the antiferromagnetic-to-ferromagnetic phase transition. Our results demonstrate an efficient approach to controlling the magnetization dynamics, thus enabling low-power tunable electronics.
INTRODUCTION
Damping, which is a fundamental parameter that defines the magnetization relaxation process, plays a crucial role in the performances of spintronics, magnetic sensors, and magnon devices (1–6). For example, the spin transfer torque magnetic random access memory (STT-MRAM) devices favor a low damping parameter for a small critical switching current, while a high damping parameter is also desired for reaching an ultrafast switching speed (3). To meet both criteria in STT-MRAM and to achieve tunable microwave and magnon devices, dynamic manipulation of magnetic damping is critical.
The current-induced spin-orbit torque generated by spin-Hall source materials can be used to modulate the effective magnetic damping parameter in ferromagnets (4, 7, 8). The more energy-efficient approach, namely, voltage control of magnetic damping, has been demonstrated with electric field effects (9, 10) and piezo strain effects (11–14) by modulating the intrinsic part of the damping (Gilbert damping) in various material systems; however, their tunability has been limited. The extrinsic part of the damping can be contributed from the spin pumping in ferromagnets/spin sink bilayer structures. In such damping process, the angular momentum is transferred from the ferromagnet into the spin sink at ferromagnetic resonance (FMR) (2, 15). It has been shown recently that spin pumping can occur laterally in FeRh (16), a metallic alloy that has ferromagnetic and antiferromagnetic phases coexisting during the metamagnetic phase transition, where the spin current is transported from its ferromagnetic domain into the surrounding antiferromagnetic spin sink (as shown schematically in Fig. 1A). During the antiferromagnetic-to-ferromagnetic phase transition of FeRh, the damping would change drastically because the lateral spin pumping is sensitive to the fraction of the ferromagnetic phases. This suggests an approach to modulate the damping dynamically by small external perturbations. Here, we demonstrate the modulation of the lateral spin pumping in epitaxial thin films of FeRh on piezoelectric substrates by controlling the ferromagnetic phase fraction using the electric field–induced strain, which results in a reversible modulation of the effective magnetic damping of more than 120% during the phase transition. This is evidenced by the strong correlation between the electric-field–induced magnetic damping and the change of ferromagnetic domain fraction, both probed by the FMR. Our analytical model suggests that the spin pumping is the dominant contribution to the damping enhancement, demonstrating the role of mixed-phase coexistence, more importantly the coupling between the two phases, in achieving high tunability of damping. This study not only provides an efficient approach to control the magnetic damping that is useful for tunable magnetic devices but also demonstrates the FeRh/piezoelectric heterostructure to be a good platform to explore the spin dynamics and the interfacial spin transport.
(A) Schematic illustration of the FeRh/PMN-PT heterostructure where the antiferromagnetic (AFM) domains embedded in the ferromagnetic (FM) matrix of FeRh were driven into FMR by microwave magnetic fields hrf. The electric field (E) was applied across the thickness direction of the PMN-PT substrate to induce piezo strains. Bottom schematics show the isothermal growth of antiferromagnetic domains in the FeRh under relatively large electric fields. The curves in the blue region represents the profile of spin current density across the lateral direction of the antiferromagnetic domain, where |
FeRh with the ordered CsCl-type structure exhibits a first-order metamagnetic phase transition from an antiferromagnetic order with a G-type spin structure to a ferromagnetic order above the room temperature (17). The transport (18, 19), spin and orbital moment (20–22), FMR (16, 23, 24), and the ultrafast dynamics (25, 26) have been investigated in FeRh thin films across the phase transition, in which the unique magnetic phase transition and the coexisting phases trigger tremendous interest in heat-assisted magnetic recording (27–29), voltage-controlled magnetism (30), and antiferromagnetic spintronics (31–36). The magnetic phase transition is accompanied by a volume expansion of 1%, suggesting strong coupling between the magnetic ordering and the unit cell structure (37–41). By using piezoelectric substrates, electric-field control of the magnetic phase transition in FeRh has been realized (42), probed by the electric field dependence of magnetization (43–47) and resistivity measurements (48, 49).
RESULTS
In this study, highly [001]-oriented FeRh thin films were grown on (001)-oriented PMN-PT (0.72PbMg1/3Nb2/3O3-0.28PbTiO3) substrates by DC magnetron sputtering (Materials and Methods), which is confirmed by the out-of-plane x-ray diffraction with the observation of only (001) and (002) reflections from the FeRh thin film (see the Supplementary Materials).The details of the film growth and structural characterizations can be found in Materials and Methods and the Supplementary Materials (48). The phase transition temperature of the FeRh on PMN-PT substrates was ~380 K upon heating (from antiferromagnetic to ferromagnetic) and ~360 K upon cooling (from ferromagnetic to antiferromagnetic), confirmed by both resistivity and magnetization measurements (see the Supplementary Materials). We performed the FMR measurements on the FeRh/PMN-PT bilayers at various temperatures across the phase transition and under different applied electric fields. The FMR spectrum was acquired by measuring the field derivative of the power absorption intensity at a fixed frequency of 9.5 GHz while sweeping the external magnetic field along the film in-plane [100] direction. The PMN-PT substrates were poled before the FMR measurements.
Figure 1B (red curve) shows a typical FMR spectrum of the 50-nm FeRh thin film on PMN-PT at 380 K with 0 V μm−1 applied electric field, which was fitted with the derivative of the Lorentzian function. From the fitting, we determined the resonance peak-to-peak linewidth μ0ΔH, which is related to the effective magnetic damping parameter αeff through
The link between the αeff and the antiferromagnetic/ferromagnetic domain formation was further confirmed by estimating the ferromagnetic phase fraction from the FMR absorption curve obtained by integrating the FMR spectrum over the field (see the Supplementary Materials). The area under the FMR absorption curve represents the total magnetic moment in the sample and, therefore, the fraction of ferromagnetic domains, given that the saturation magnetization of each domain only decreases slightly with increasing temperature. Figure 1D shows the estimated ferromagnetic phase fraction as a function of the temperature, where the drastic increase of the ferromagnetic phase fraction is coincident with the decrease of the αeff. Below 340 K upon cooling (365 K upon heating), although the ferromagnetic phase still exists, the weak FMR amplitude and the wide linewidth make the fitting of the FMR spectrum unreliable.
Having established the FeRh magnetization dynamics across the magnetic phase transition at zero electric field, we then investigate their dependence on electric fields. The blue curve in Fig. 1B shows the FMR spectrum at 380 K upon heating with an applied electric field of 0.67 V μm−1 on the PMN-PT substrate. Noticeably, the broadening of the FMR linewidth and the decrease of the FMR amplitude are prominent. This is attributed to the electric-field–induced in-plane biaxial compressive strain that is expected to promote the formation of the antiferromagnetic phase (48, 49) and shift the magnetic phase transition to higher temperatures given the fact that the phase transition is extremely sensitive to the strain. Figure 1 (C and D) summarizes the temperature dependence of αeff (Fig. 1C) and the ferromagnetic phase fraction (Fig. 1D) with an applied electric field of 0.67 V μm−1, both showing the electric-field–induced shift of the magnetic phase transition temperature of ~15 K, which is consistent with our resistivity measurements (see the Supplementary Materials). Through this tuning mechanism, the modulation of effective damping (∆αeff) and the ferromagnetic phase fraction becomes more prominent around the phase transition temperatures.
For better visualization, we extracted the modulation of αeff (Fig. 2A) and the ferromagnetic phase fraction (Fig. 2B) by the electric field as a function of temperature. Both ∆αeff and the ferromagnetic phase modulation peak at the magnetic phase transition temperatures, demonstrating again the correlation between the magnetic damping enhancement and the mixture of ferromagnetic and antiferromagnetic phases, which causes the lateral spin pumping and the two-magnon scattering. It should also be noted that the ∆αeff is nonzero (Fig. 2A) even further above the phase transition temperature (>400 K) when the FeRh ferromagnetic phase modulation by electric field is essentially negligible (Fig. 2B). This small change in ∆αeff at a full ferromagnetic phase in FeRh could be attributed to the strain modulation of the Gilbert damping parameter or the piezoelectric surface morphology–induced inhomogeneous linewidth broadening, which has been reported in various magnetoelectric heterostructures (11, 13).
(A and B) Change of the effective magnetic damping (A) and the ferromagnetic phase fraction induced by the applied electric field of 0.67 V μm−1 as a function of the temperature. The solid curve in (A) represents the analytical model fitting.
We further studied the isothermal electric-field modulation of αeff and the ferromagnetic phase fraction at 380 K during heating (where ∆αeff peaks), shown in Fig. 3 (A and B). The strain-mediated ferromagnetic phase fraction in FeRh gives rise to a record high isothermal electric-field control of a relative change in damping (∆αeff/αeff,min) of ~120%. This tunability is higher than that reported in magnetic semiconductors (10% tunability) (10) or in ferromagnetic metals (50% tunability) (11). The butterfly-like feature shown in both curves mimics the behavior of the electric field–induced strain in piezoelectric substrates, evidencing a strain-mediated effect. At large applied electric fields (±0.67 V μm−1), the compressive piezo strain favors the antiferromagnetic phase (Fig. 3B), enhancing αeff in FeRh. While at the coercive fields (0.1 V μm−1) of the PMN-PT substrate, its polarization and induced strain are both zero, which promotes more uniform ferromagnetic domains (Fig. 3B) and therefore leads to a decrease in αeff.
Effective magnetic damping (A) and the ferromagnetic phase fraction (B) as a function of applied electric fields at 380 K during heating. The solid curve in (A) represents the analytical model fitting.
DISCUSSION
As the ferromagnetic phase fraction decreases, the feature size of ferromagnetic domains δ shrinks (bottom schematic of Fig. 1A). We assume the change of intrinsic Gilbert damping to be independent of δ, as the Gilbert damping is shown to be proportional to the resistivity in FeRh (16), and the change of the resistivity (see the Supplementary Materials) is much smaller than that of the damping. Then, the enhancement of αeff, which is mainly contributed from the spin pumping and the two-magnon scattering, can be expressed as
The present results demonstrate an effective approach to dynamically and reversibly modulate the effective magnetic damping up to 120% via the electric field–driven magnetic phase transition in FeRh/PMN-PT heterostructures. This large electric field control of spin dynamics is desirable for low-power spintronics, magnonics, and magnetic microwave devices, which may open opportunities, i.e., in STT-MRAM with tunable switching speed and critical current and in voltage-tunable spin waves for logic applications (5). Furthermore, we anticipate that the FeRh/piezoelectric hybrid structures to be a versatile platform for the fundamental understanding of damping mechanism in ferromagnets.
MATERIALS AND METHODS
Film growth and structural characterization
FeRh thin films were grown on the (001)-oriented PMN-PT (0.72PbMg1/3Nb2/3O3-0.28PbTiO3) substrates by DC magnetron sputtering using a stoichiometric Fe0.5Rh0.5 target in a vacuum chamber with a base pressure of ~10−8 torr. The film deposition was performed at an Ar gas pressure of 3 mtorr and at a substrate temperature of 375°C, followed with an in situ annealing at 500°C for 1 hour. The epitaxial arrangement and crystalline orientation of the thin films were confirmed by x-ray diffraction (see fig. S1). The thickness of the film was determined by x-ray reflectivity. Rutherford backscattering spectrometry indicates that the chemical composition of FeRh is Fe0.51Rh0.49. Further details of FeRh thin film deposition and structural characterization can be found in our previous report (48).
Electrical transport measurements
Electrical resistivity measurements were carried out on as-grown full films in a van der Pauw geometry. DC electric fields of up to ±0.67 V μm−1 were applied between the top FeRh film and the bottom Au electrode on the underside of the PMN-PT substrate with the bottom electrode grounded.
Magnetic property measurements
FMR measurements were performed using a Bruker EMX electron spin resonance spectrometer with a TE102 cavity and a field sweeping mode at a fixed frequency of 9.5 GHz. Static magnetization measurements were carried out in a superconducting quantum interference device magnetometer.
Analytical model fitting and interpretation
Equation 1 depicts how the modulation of effective damping parameter ∆αeff varies with the feature size of ferromagnetic domain δ. However, Figs. 2 and 3 only provide data of ∆αeff versus the ferromagnetic phase modulation (denoted as x ∈ [0, −100%]). Therefore, to fit these data using Eq. 1, a knowledge of the relationship between δ and x is necessary. Although quantitative determination of such relationship may require detailed magnetic imaging studies [e.g., (20) and (22)] as a function of both temperature and electric field, qualitatively, it is clear that δ should decrease when the magnitude of x increases. Guided by this principle, we tested different types of mathematical functions including exponential, logarithmic, and asymptotic functions. It is found that the simple asymptotic function δ = − L/x yields the most accurate fit (with coefficient of determination R2 = 0.98). By further considering that antiferromagnetic domain size [denoted as δAFM, from tens of nanonmeters to hundreds of nanometers; (20) and (22)] is substantially larger than the spin diffusion length in antiferromagnetic FeRh [λ is typically small in heavy metals because of the strong spin-orbit coupling, normally below 10 nm (53)], the bulk spin conductance of the antiferromagnetic phase
On the basis of the arguments above, Eq. 1 can now be rewritten as
Equation 2 can be directly used to fit the experimental data of ∆αeff versus x in Figs. 2 and 3. The fitting yields ∆αeff = Ax + Bx2 (A = −9.489 × 10−2; B = 2.861 × 10−2), which allows us to evaluate the individual contributions of two-magnon scattering and lateral spin pumping as a function of x. Specifically, the contribution of lateral spin pumping is Ax/(Ax + Bx2), which remains above 81% even when x reduces to −77% (the largest modulation among our data). Moreover, A = −9.489 × 10−2 means that
For pure antiferromagnetic phase of FeRh, resistivity ρ = 150 × 10−8 ohms·m is extracted from fig. S2; saturation magnetization Ms increases moderately from 1.252 to 1.303 MA/m when temperature cools from 385 to 345 K (16). For simplicity, Ms is fixed at 1.25 MA/m, which is the value at 380 K under which we performed the isothermal electric field tuning of effective damping experiment. Rearranging Eq. 3 gives rise to the relation among
Equation 4 indicates that
Equation 4 also indicates that 8.6151 × 1018(m−2)Lλ < 1 to ensure
Following the abovementioned analysis of estimating the lower bound of L, one can meanwhile estimate the upper bound of the spin diffusion length λ via Eq. 5, λ < 1.15 nm. Overall, the analyses above reasonably arrive at two conclusions. First, the real component of the spin-mixing conductance at the ferromagnetic/antiferromagnetic interface in FeRh has a lower bound of
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/40/eabd2613/DC1
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).