Dynamical amplification of magnetoresistances and Hall currents up to the THz regime

TL;DR

This paper demonstrates that dynamical effects in Co/Pt and Fe/W bilayers can significantly amplify magnetoresistances and Hall currents at THz frequencies, with phase differences playing a crucial role in these phenomena.

Contribution

It reveals how phase differences influence dynamical currents and shows giant enhancements of spin Hall angles near ferromagnetic resonance in bilayer systems.

Findings

Dynamical spin Hall angles can change sign and increase over 500%.

Charge and spin pumping can significantly modify magnetoresistive and Hall effects.

Giant dynamical effects occur near ferromagnetic resonance in bilayers.

Abstract

Spin-orbit-related effects offer a highly promising route for reading and writing information in magnetic units of future devices. These phenomena rely not only on the static magnetization orientation but also on its dynamics to achieve fast switchings that can reach the THz range. In this work, we consider Co/Pt and Fe/W bilayers to show that accounting for the phase difference between different processes is crucial to the correct description of the dynamical currents. By tuning each system towards its ferromagnetic resonance, we reveal that dynamical spin Hall angles can non-trivially change sign and be boosted by over 500%, reaching giant values. We demonstrate that charge and spin pumping mechanisms can greatly magnify or dwindle the currents flowing through the system, influencing all kinds of magnetoresistive and Hall effects, thus impacting also dc and second harmonic experimental measurements.

Tables

Table 1. Table 1: Magnetoresistances Δ ​ I α , β , γ = I ω , ∥ ′ ​ ( 90 ∘ ) − I ω , ∥ ′ ​ ( 0 ∘ ) \Delta I^{\alpha,\beta,\gamma}=I^{\prime}_{\omega,\|}(90^{\circ})-I^{\prime}_{\omega,\|}(0^{\circ}) for applied dc and resonant electric fields for the different rotation angles represented in the central panel of Fig.  3 .

	Frequency	$\Delta I^{\alpha }$ (xy plane)	$\Delta I^{\beta }$ (zy plane)	$\Delta I^{\gamma }$ (zx plane)
Co/Pt(001)	$\omega =0$ (dc)	-1.86	-0.99	0.89
	$\omega ={\omega }_{\text{res}}$ (ac)	-13.22	-7.37	5.85
Fe/W(110)	$\omega =0$ (dc)	0.52	0.91	0.39
	$\omega ={\omega }_{\text{res}}$ (ac)	-0.26	5.63	5.89

Table 2. Table 2: Maximum amplitude of Hall currents for applied dc and resonant electric fields. When the magnetization is rotated in the x ​ y 𝑥 𝑦 xy plane (i.e., PHE symmetry), we also indicate in parenthesis the angle for which the current is maximum. The minus sign represents the change from the reference case of static Co/Pt(001) currents.

	Frequency	$I_{\omega ,⟂}^{\prime }\left({𝐦}_0\in xy\right)$	$I_{\omega ,⟂}^{\prime }\left({𝐦}_0\parallel \widehat{𝐳}\right)$
Co/Pt(001)	$\omega =0$ (dc)	0.05 ($\varphi ={35}^{\circ }$)	0.44
	$\omega ={\omega }_{\text{res}}$ (ac)	0.22 ($\varphi ={42}^{\circ }$)	0.54
Fe/W(110)	$\omega =0$ (dc)	0.09 ($\varphi ={55}^{\circ }$)	0.22
	$\omega ={\omega }_{\text{res}}$ (ac)	0.36 ($\varphi ={69}^{\circ }$)	-0.29

Equations

$$\begin{array}{rl}{\mathbf{I}}_{\text{C}}(t)& =\mathrm{Re}\{(I_{\omega ,\Vert }\hat{\mathbf{x}}+I_{\omega ,\perp }\hat{\mathbf{y}})E_0{e}^{i\omega t}\}.\end{array}$$

$$\begin{split}\mathbf{I}^{\text{dyn}}_{\text{C}}(t)=&\left\{\Delta I_{\text{1}}(\omega)\left(m^{2}_{0}-m^{y2}_{0}\right)\mathbf{\hat{x}}+\Big{[}\Delta I_{\text{2}}(\omega)m^{z}_{0}+\Delta I_{\text{1}}(\omega)m^{x}_{0}m^{y}_{0}\Big{]}\mathbf{\hat{y}}\right\}E_{0}\cos(\omega t)\end{split}$$

Full text

Introduction

The interrelation between magnetic properties and electric currents has been studied for more than a century Thomson (1856); Jan (1957). The electric resistance of magnetically ordered systems, for example, depends on the direction of their magnetization, a phenomenon given the general name of magnetoresistance. Systems exhibiting relatively large magnetoresistance are used in sensing devices for reading information stored in magnetic media. Understanding how the spin of the electron can affect its propagation in materials with strong spin-orbit interaction is the main focus of various recent developments. A plethora of novel magnetoresistive effects have been discovered in the last years: spin Hall Nakayama et al. (2013); Chen et al. (2013), unidirectional spin Hall Avci et al. (2015) and, independently, linear spin Hall Olejník et al. (2015), as well as Hanle Vélez et al. (2016), and Rashba-Edelstein effects Nakayama et al. (2016). A common characteristic shared by all of them is the appearance of non-equilibrium spin polarizations accumulated at the interfaces/edges of those systems due to the flow of an electric current.

This spin accumulation can also be used to excite spin waves in adjoining magnetic units, by means of the so-called spin-orbit torques, when time-dependent (ac) currents are applied Fang et al. (2011); Garello et al. (2013); Freimuth et al. (2014, 2015); Guimarães et al. (2015). These torques can be effective even at room temperatures and may be used for magnetization switching purposes Miron et al. (2011); Pai et al. (2012); Liu et al. (2012); Garello et al. (2014); Ciccarelli et al. (2016); Fukami et al. (2016). In certain metallic systems, the induced spin-orbit torques and magnetoresistive effects can reach fairly high intensities Miron et al. (2010); Kim et al. (2016). Although the connection between the magnetoresistance and static torques is well-established, a dynamical connection between these two effects is missing. Unraveling how the induced dynamics influences the magnetoresistance and the Hall currents are essential questions, especially for prospective spintronics applications within the THz frequency range.

The effective magnetic fields that drive the magnetization into precession in those systems are obtained experimentally through second order signals, either dc (rectified) Fang et al. (2011); Liu et al. (2012); Kurebayashi et al. (2014) or doubled-frequency (second-harmonic) voltages Garello et al. (2014); Avci et al. (2014); Hayashi et al. (2014). These measurements depend both on the transverse magnetization components, and on the magnetoresistance (or Hall resistance). The latter is assumed to be independent of both the frequency of the oscillatory electric field and intensity of the bias magnetic field that may be swept through the ferromagnetic resonance. The measurements are usually described by simple functional forms Jan (1957); McGuire and Potter (1975), and the characteristics of the magnetization dynamics are studied within a semi-classical approach. In addition, even though the effect of the magnetization precession on second order signals is commonly accessed, measurements of its feedback in the first order signal are scarce Kupferschmidt et al. (2006); Kovalev et al. (2007), specially for the in-plane geometry Tserkovnyak and Bender (2014); Ciccarelli et al. (2014), where the pumped currents may give rise to important contributions yet to be explored.

In this work, we investigate the intrinsic feedback of the current-induced magnetization dynamics to the electric current that drives it. By applying a time-dependent electric field to ferromagnetic/heavy metal bilayer systems such as Co/Pt(001) and Fe/W(110), we demonstrate the emergence of an enhanced dynamical magnetoresistance originating from spin and charge pumping. This can be one order of magnitude larger than the dc values. The magnetization dynamics also induces currents along the transverse directions that substantially influence Hall effects. We show that charge to spin conversion in the dynamical regime hinges crucially on signal dephasing, implying that the spin Hall angle is a complex quantity. It is characterized by an amplitude and a phase difference, contrary to what is assumed in common phenomenological approaches. This complex nature allows striking variations — up to 500% — on the generated transverse spin currents. Our approach is based on realistic descriptions of the electronic structure of the system. The theory describes collective spin excitations in the presence of spin-orbit interaction, and naturally takes into account all the intrinsic spin-orbit-related mechanisms, with no need for any adjustable parameters.

In the systems we have considered, we further show that magnetic excitations are the main source — but not the only one — of the dynamic contributions. They can be manipulated by varying the frequency and intensity of the applied fields, which can also affect second order signals measured in experiments. Moreover, since the spin-pumping contribution is switched on or off by changing the magnetization direction, it leads to an angular dependence of the currents that deviates starkly from the static counterpart. This change in current signals provides an efficient way of manipulating magnetic states and establish a powerful ingredient for the development of spintronic devices.

Results

Conceptual framework

We consider metallic bilayer heterostructures consisting of a ferromagnetic metal deposited on a heavy metal substrate. The systems are subjected to a time-dependent uniform electric field given by $\mathbf{E}(t)=\operatorname{Re}[E_{0}e^{-i\omega t}\mathbf{\hat{x}}]$, where $E_{0}$ is the field amplitude, and $\omega$ its angular frequency. An electric current may be measured along the longitudinal (x) or along the transverse (y) direction. The induced first order electric currents flowing through the system can be summarized in vector form as

[TABLE]

where $I_{\omega}=I^{\prime}_{\omega}+iI^{\prime\prime}_{\omega}$, with $I^{\prime}_{\omega}$ and $I^{\prime\prime}_{\omega}$ expressing the in-phase and out-of-phase currents per electric field, respectively. When $\omega\rightarrow 0$, we have $I^{\prime\prime}_{\|}=I^{\prime\prime}_{\perp}=0$. The quantity $I_{\omega}$ is directly related to the conductance $G$ and conductivity $\sigma_{0}$ of the system. Assuming a uniform in-plane electric field, $I_{\omega}=G.d=\int_{S}\sigma_{0}dA$, where $d$ is the length of the sample along the direction of the applied electric field, and $S$ is the cross-sectional area.

A schematic setup and our choice of axes are illustrated in Fig. 1. The total electric current flowing across the system can be split into two parts: $\mathbf{I}_{\text{C}}(t)=\mathbf{I}_{\text{C}}^{0}(t)+\mathbf{I}^{\text{dyn}}_{\text{C}}(t)$, where $\mathbf{I}^{0}_{\text{C}}(t)$ is the current driven directly by the electric field, and $\mathbf{I}^{\text{dyn}}_{\text{C}}(t)$ is the dynamical part produced by the magnetization precession. In the presence of spin-orbit coupling, $\mathbf{I}^{0}_{\text{C}}(t)$ depends upon the magnetization direction, giving rise to the anisotropic magnetoresistance (AMR) and also to transverse components due to the planar Hall effect (PHE) and anomalous Hall effect (AHE) McGuire and Potter (1975); Nagaosa et al. (2010); Seemann et al. (2011); Bouaziz et al. (2016). Here, the spin Hall magnetoresistance (SMR) Nakayama et al. (2013); Chen et al. (2013) is also included in $\mathbf{I}^{0}_{\text{C}}(t)$. The frequency dependence of $\mathbf{I}^{0}_{\text{C}}(t)$ can be well described by the Drude model in metallic materials, and exhibits very little variation in the range of frequencies we are interested in.

The dynamical contribution, $\mathbf{I}^{\text{dyn}}_{\text{C}}(t)$, originates from a combination of several mechanisms as follows: The spin Hall and inverse spin galvanic effects induce an oscillatory spin accumulation at the heavy metal surfaces, as illustrated on Fig. 1a. Then, spin-orbit torques may set the magnetization into precessional motion, depicted in Figs. 1(b-d). Note that when $\mathbf{m}_{0}\|\hat{\mathbf{y}}$, no precession is induced. Finally, spin pumping followed by inverse spin Hall effect and charge pumping generated by spin galvanic effects produce an additional charge current, as shown in Figs. 1(c,d).

To reveal the spin-orbit origin of these effects, they are usually given in terms of the spin Hall angle $\theta_{\text{SH}}$. This quantity is defined as $\theta_{\text{SH}}\propto I_{\perp}^{\text{S}}/I_{\|}^{\text{C}}$, where $I_{\perp}^{\text{S}}$ is the transverse spin current with out-of-plane spin polarization and $I_{\|}^{\text{C}}$ is the longitudinal charge current. It is a real number for dc currents, and measures how good the material is as a charge-to-spin converter. Phenomenologically, the symmetry of the dynamical contribution is similar to the SMR, and can be written as

[TABLE]

This formula was obtained by considering the charge current converted from the pure spin current ($\propto\theta_{\text{SH}}$) that was pumped away from the magnetic layer due to the precessing magnetization Tserkovnyak et al. (2002). The motion of the magnetization is described by introducing the even and odd spin-orbit torques Freimuth et al. (2014) (also $\propto\theta_{\text{SH}}$) into the Landau-Lifshitz-Gilbert equation. Therefore, $\Delta I_{\text{1,2}}(\omega)$ is proportional to $\theta_{\text{SH}}^{2}$, representing the two charge$\leftrightarrow$spin conversions involved in the effect (see Ref. 28 and Supplementary Note S1). In this approach, $\theta_{\text{SH}}$ is assumed to be real.

The longitudinal and transverse contributions obtained in Eq. 2 are similar to the ones obtained for the SMR Chen et al. (2013). While the current flowing in the direction of the applied electric field ($\mathbf{\hat{x}}$) depends on the $y$ component of the magnetization (i.e., the direction of the spin accumulation induced by the spin-orbit coupling), the transverse current (flowing along the $\mathbf{\hat{y}}$ direction) has symmetries similar to the AHE ($\propto m_{z}$) and to the PHE ($\propto m_{x}m_{y}$). It is important to notice that, since this contribution involves the magnetization precession, $\Delta I_{1,2}$ are frequency-dependent. A detailed explanation of this contribution can be found in Supplementary Note S1.

To investigate the features of the dynamical effects described above, we use an ab-initio-based approach in which we map the electronic structure obtained in a density functional theory calculation into a multi-orbital tight-binding hamiltonian. The magnetic excitations and currents are obtained using linear response theory (see Methods and Supplementary Note S2).

Dynamical currents and spin Hall angles

We chose two bilayer metallic systems, Co/Pt(001) and Fe/W(110), manifesting different aspects of the phenomena due to the character of their charge to spin conversion. There are two main differences between these systems. First, Co/Pt(001) presents an uniaxial perpendicular anisotropy, while Fe/W(110) has a biaxial in-plane anisotropy with the long axis being the lowest-energy direction. The second difference is the band filling of the substrate, which, according to Hund’s rule, leads to parallel (anti-parallel) orbital and spin moments in Pt (W) due to the spin-orbit coupling Hoffmann (2013). In our approach, the position of the resonance frequency is intrinsically determined by the electronic structure, without any adjustable parameter. The magnetocrystalline anisotropy energies can be estimated through the modified Kittel’s formula Smit and Beljers (1955); Farle (1998) or the static transverse magnetic susceptibility dos Santos Dias et al. (2015). The values obtained using the latter and further details about the ground state properties of these systems are listed in Supplementary Note S3. The large anisotropies of the systems we investigate are mainly due to the fact that we are considering a single atomic ferromagnetic layer (See, for example, Refs. 38; 39). The strength and sign of the anisotropy energy can be tailored by changing the thickness and/or the composition of the two constituents forming the bilayer. For example, increasing the thickness of the FM layer lowers the effective anisotropy, which moves the resonance frequency to the GHz range.

Pt and W are known to have strong spin-orbit coupling and relatively high spin Hall angles. We have calculated $\theta_{\text{SH}}$ for slabs containing 10 layers of pure Pt(001) and pure W(110) and obtained the values 0.07 and -0.03, respectively, in accordance with experimental observations Rojas Sánchez et al. (2014); Hoffmann (2013). For a monolayer of Co deposited on Pt(001), we have found values of $\theta_{\text{SH}}\sim 0.06$ for Pt thicknesses ranging from 4 to 9 atomic planes, in fair agreement with recent experimental estimatesGarello et al. (2014).

We start by exploring the behavior of the longitudinal and transverse currents, calculated as functions of the applied electric field frequency, for different thicknesses of the heavy metal substrate. The time-dependent electric field is applied along the [110] and [1$\bar{1}$0] directions for the Co/Pt(001) and Fe/W(110) systems, respectively, which defines the corresponding $x$ axis in each case. In Fig. 2, we plot the total first order in-phase currents along the directions parallel and perpendicular to the applied field, $I^{\prime}_{\omega,\|}$ and $I^{\prime}_{\omega,\perp}$, respectively. As expected from symmetry considerations, for Fe/W(110) with the electric field parallel to the magnetization direction, $I^{\prime}_{\omega,\perp}=0$ (as illustrated in Fig. 1c). A drastic variation, however, is observed for all the other currents as the frequency is varied. The structures manifested in the currents are related to the ferromagnetic resonances of the magnetic layers, located at frequencies determined by the magnetic anisotropy (calculations for Fig. 2 were performed for zero external magnetic field). They induce substantial changes owing to the pumping mechanisms described earlier, which clearly affect both the longitudinal and transverse currents. We observe that the dynamical contributions for the longitudinal currents displayed in Figs. 2(a,c) exhibit symmetric line shapes in the vicinity of the resonance frequencies, whilst the line shapes of the transverse currents presented in Fig. 2b are clearly asymmetric. This can be explained by the properties of the charge-spin response function (see Supplementary Note S4). The thickness-dependent variations on the features shown in Fig. 2 reflect quantum well effects present in our atomically thin films. These could be explored experimentally, e.g., by having ultrathin films of Pt or W deposited on an insulator.

Surprisingly, from the results of Figs. 2(a,c), we note that the longitudinal currents calculated for Co/Pt(001) exhibit a maximum, whereas the corresponding ones for Fe/W(110) display a minimum. This means that the dynamical contribution to $I^{\prime}_{\omega,\|}$ has opposite signs for Co/Pt(001) and Fe/W(110). Such a sign change has been experimentally observed in Ref. 41 for permalloy deposited on Pt and Ta, and seems to contradict the phenomenological expectations from Eq. 2. It turns out that the issue is not the postulated quadratic dependence on the spin Hall angle, but neglecting its dynamical nature: we must take into account the dephasing between the charge and spin currents, i.e., the spin Hall angle must be defined as a complex number. Figures 2(d,e) show the amplitude and phase of the spin Hall angle, respectively, obtained for Co/Pt(001), and Figs. 2(f,g) represents the same quantities for Fe/W(110). They reveal that the magnitude of the spin currents flowing perpendicularly to the electric field can be controlled by changing either its frequency or the strength of an applied external magnetic field (i.e., by tuning the ferromagnetic resonance frequency). On the one hand, for Co/Pt(001) the spin Hall angle amplitude decays to rather low values, close to turning off the flow of spin currents (see Fig. 2d). On the other hand, for Fe/W(110) it can reach giant values of $|\theta_{\text{SH}}|\sim 0.5$ (see Fig. 2f) — an enhancement of 500% with respect to the dc value. Complementary results displaying in-phase and out-of-phase spin Hall angles can be found in the Supplementary Figure S1. We remark that the previous considerations on the importance of dephasing also applies to the other processes involved in the generation of the ac currents (such as the spin accumulation and the spin-orbit torques).

An advantage of our approach is that all dephasing effects are naturally taken into account, and based on the underlying electronic structure. Since these factors are also layer-dependent and include the oscillatory induced magnetic moments, they may give different contributions to the total currents. For instance, our results show that increasing the thickness of the Pt(001) substrate leads to an increase of the current flowing through the Co layer, while the current on the Fe layer is quite insensitive on the thickness of W(110). The importance of how the current is distributed between the ferromagnetic and heavy metal layers has been raised in Refs. 42; 43. Due to the non-homogeneous distribution of the currents along the layers, a time-dependent Oersted field can also be generated inside the sample. Nevertheless, in ultrathin films like the ones we investigate, it was shown and observed that this contribution is negligible compared to the spin-orbit torques Miron et al. (2011); Garello et al. (2013); Ciccarelli et al. (2014).

Angular dependence of magnetoresistances and Hall effects

The excitation of the spin waves depends on the relative angle between the magnetization and the induced spin accumulation at the interface of the bilayer. Hence, the angular dependence of the currents with respect to the magnetization direction — which defines the magnetoresistance ratios and the Hall currents — can be drastically modified. To investigate this effect, we focus on bilayers consisting of a ferromagnetic monolayer deposited on four layers of heavy metal substrate. As in experimental setups, we consider a magnetic field larger than the anisotropy, to orient the magnetization along any chosen direction. Since increasing the static effective magnetic field diminishes the ability of the ac electric field to drive the magnetization precession, the dynamical charge currents are even stronger for smaller and more realistic external fields (see Supplementary Figure S2).

To demonstrate the dynamical effects on the magnetoresistances, we plot in Figs. 3(a,b) the angular dependencies of the in-phase longitudinal currents for an applied dc electric field ($\omega=0$) and for a resonant ac field ($\omega=\omega_{\text{res}}$). In Table 1, we summarize the values of the variations $\Delta I^{\alpha,\beta,\gamma}=I^{\prime}_{\omega,\|}(90^{\circ})-I^{\prime}_{\omega,\|}(0^{\circ})$, where the angles $\alpha,\beta$ and $\gamma$ are illustrated on the central panel of Fig. 3. The large magnetoresistance enhancements (with respect to the static case) in the vicinities of the resonant frequencies are evident. As expected from the results shown in the previous section, the average current increases for Co/Pt(001) and decreases for Fe/W(110). Current variations that were almost vanishing in the dc limit (as the SMR and AMR changes depicted by the dashed blue curve in Fig. 3a, and red curve in Fig. 3b, respectively) become much larger than any static alterations. It is worth noticing that the ac electric field can even change the sign of some of the magnetoresistance components.

The naive expectation is that the change in the currents is due to a modulation of the static magnetoresistance by the amplitude and shape of the magnetization precession Kim et al. (2012). This is true for frequencies away from the resonance, when the amplitude of precession and the pumped currents are small. Nevertheless, this is clearly not the case for the ac currents presented in Fig. 3, otherwise the upper limit for the effect amplitude would be the dc magnetoresistance.

The curves displayed in Fig. 3 highlight the importance of higher order terms in the pumped currents, especially in the dynamical regime. Even though some of the obtained results for the $C_{4v}$ symmetry of Co/Pt(001) follow the behavior predicted by Eq. 2, those for the $C_{2v}$ symmetry of Fe/W(110) definitely do not. Large angular dependencies on the first order voltage were measured on strained (Ga,Mn)As, also $C_{2v}$Ciccarelli et al. (2014). However, they showed marked deviations with respect to predictions using phenomenological models. These were attributed to the presence of higher order contributions on the pumped currents, which are intrinsically included in our results shown in Fig. 3b.

As expected from the phenomenology, pumped currents also influence Hall effects. In Figs. 3(c,d), we show the in-phase currents that flow along the transverse direction, and we collect their maximum amplitude in Table 2. In the dc-limit, we find that the PHE (dashed black curves) is smaller than the AHE (blue and red dashed curves) in both cases, as obtained in experimental measurements Garello et al. (2013). The dynamical contribution changes this picture and brings the maximum currents of all angular sweeps to the same order of magnitude. The angular dependence presented in Fig. 3d can be understood by analyzing the precessional motion of the magnetization (see Supplementary Note S4).

Other contributions not coming from the magnetization precession may also be relevant. As illustrated in Fig. 1b, spin waves are not excited when the magnetization points along the $y$ direction and, therefore, no dynamical contribution is expected. However, a substantial change from dc to ac excitations can be seen on the black and blue curves of Fig. 3b at $90^{\circ}$ or $270^{\circ}$. We relate this effect to a possible precession of the spin current polarization around the magnetic field, an aspect not included in Fig. 1b and in Eq. 2.

The angular variations discussed above depend on the chosen frequency of the applied field. Figure 4 presents a complete map of the angular and frequency dependencies of all the dynamical effects described above, calculated for a monolayer of Co on 4 layers of Pt(001). The black dashed lines indicate the frequency in which the angular dependencies of the dynamical currents plotted in Figs. 3(a,c) were obtained. They provide a clear picture of how the precessional motion of the magnetization influences longitudinal and charge currents. Different magnetoresistances and Hall currents can be obtained depending on the anisotropy of the system and the applied magnetic fields. Therefore, by tuning the frequency of the oscillatory field to a suitable value, one can obtain large current variations as a result of slight changes in the magnetization direction.

Discussion

We have investigated the dynamical effects on the linear response charge currents due to spin-orbit-related phenomena. They involve the simultaneous action of both spin-Hall and inverse spin galvanic effects together with their reciprocal counterparts, as in the SMR, with an additional contribution arising from the excitation of the ferromagnetic resonance by the ac spin-orbit torques. We have presented results for Co/Pt(001) and Fe/W(110) that unveil the complex nature of spin-charge excitations in these systems. They show that relevant physical mechanisms are being neglected by current phenomenological approaches, which consequently fail to properly describe the amplitude and angular dependence of the ac currents induced in such setups. Our material-dependent first-principles approach avoid these pitfalls by taking into account the lattice structure and all the dephasing factors between different physical quantities, which are intrinsically determined by the electronic structure. Although our main focus is metallic multilayers, the underlying processes discussed here also apply to other systems such as magnetic insulators in contact with heavy metals (e.g. YIG/Pt), and magnetic semiconductors with relatively large SOC, where the spin pumping contribution to the magnetoresistance can be also relevant Ciccarelli et al. (2014).

We propose that the dynamical contributions, that are present not only in the charge but also in the spin currents, can be used to manipulate the spin Hall angle of the system. This means that the flow of spin currents can be controlled by changing either the frequency of the electric field or the intensity of the magnetic field. Recently, giant spin Hall angles were measured for a few systems in the dc limit Liu et al. (2012); Pai et al. (2012); Zhang et al. (2016). Our work indicates a huge potential for magnifying even further these angles in the ac regime.

The dynamical processes we describe manifest prominently in the ac currents flowing with the same frequency as the external perturbation (first harmonic signal). Nevertheless, the frequency and field dependence of the current response may also affect dc and second harmonic signals, especially for frequencies close to the ferromagnetic resonance of the system (see Supplementary Note S5). This may be related to why the second harmonic signal, measured in the quasi-static limit, provides more reliable results than the dc measurements (see Supplemental Material of Ref. 10).

Since the main contribution to the charge current is generated by the magnetization precession, it can be deactivated by rotating the magnetization direction, substantially modifying dynamical magnetoresistances. Hall currents can also be sensitively controlled by tuning the frequency of the applied electric field in the vicinity of the ferromagnetic resonance. These mechanisms can be used to engineer drastic variations of the ac currents upon rotation of the magnetization by only a few degrees. Conversely, ac currents can be employed for manipulating magnetic states due to the enhancement of the spin-orbit torque, in particular at the ferromagnetic resonance frequency. This change in current signals provides an efficient way of manipulating and detecting magnetic states and may be an alternative to the microwave-assisted magnetic recording. A hypothetical protocol exploiting these effects could be: first apply a short weak current pulse that sets the magnetization into precessional motion, immediately followed by a larger pulse that would induce the switching. Current pulses may be used to incorporate high frequency components into the actuating current

The control of the current signal may help to unravel mechanisms to improve and measure spin-orbit torques, providing an efficient way to manipulate magnetic states, e.g., through microwave assisted recording technologies Okamoto et al. (2015). The dynamical physical concepts described above can be adjusted by choosing lattice symmetries, material elements and layer thicknesses of multilayer structures. This engineering procedure already proved effective in building efficient THz devices by taking advantage of the different spin-orbit properties of Pt and W Seifert et al. (2016). Our results demonstrate that the exploration of dynamical contribution gives an extra degree of freedom to manipulate future spintronic devices and may support ultrafast reading and writing mechanisms.

Methods

The electronic structures of the adsorbate/substrate combinations are described by a multi-orbital tight-binding hamiltonian with nine orbitals (the five $d$ and the $sp$ complex) and the two spin components included per atom. The hopping matrices are obtained from density functional calculations based on the real-space linear-muffin tin orbitals method, within the atomic sphere approximation (RS-LMTO-ASA) Andersen and Jepsen (1984); Peduto et al. (1991); Frota-Pessôa (1992). Mean field theory is used to describe the ground state of the system, considering an effective intra-atomic Coulomb interaction $U=$1\text{,}\mathrm{eV}$$ within the d shells of both the ferromagnetic and heavy metal layers. The ground-state magnetization of the ferromagnetic layer is calculated self-consistently, assuming a fixed value of the Fermi energy and adjusting the atomic $d$-orbital energy levels to reproduce the electronic occupancy obtained by the RS-LMTO-ASA calculations for each layer Barbosa et al. (2001). At the end of this process the atomic planes of the non-magnetic substrate close to ferromagnetic layer become also slightly spin polarized due to the proximity effect. Spin-orbit interaction is incorporated within the local (on-site) approximation $\hat{H}_{\text{SO}}=\sum_{i}\lambda_{i}\mathbf{L}_{i}\cdot\mathbf{S}_{i}$, explicitly involving the $d$-band complex only. All coupling constants $\lambda_{i}$ are obtained from RS-LMTO-ASA first principles calculations. Spin rotational symmetry is broken in the presence of SOC, and in the cases we consider, it leads to a perpendicular magnetic anisotropy for the Co/Pt(001) system and the [1$\bar{1}$0] direction being the easy axis for Fe/W(110).

Our description of the excitations is based on linear response theory using the Kubo formalism Kubo (1957), and follows through application of the RPA. Due to spin-orbit coupling, the responses of the transverse components of the magnetization are coupled to the longitudinal ones, and we end up with a $4\times 4$ susceptibility matrix in spin space whose elements need to be calculated simultaneously Costa et al. (2006, 2010).

A more detailed description of the theory is given in the Supplementary Note S2.

Acknowledgments

We are very grateful to A. B. Klautau for providing the tight-binding parameters, to J. R. Suckert and S. Blügel for discussions, and for computing time on JURECA and JUQUEEN supercomputers at Jülich Supercomputing Centre. We thank the support CNPq (Brazil) and Alexander von Humboldt foundation (Germany). This work is supported by the HGF-YIG Programme VH-NG-717 (Functional Nanoscale Structure and Probe Simulation Laboratory-Funsilab) and the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (ERC-consolidator grant 681405 – DYNASORE).

Author contributions

F.S.M.G. performed calculations and wrote the manuscript. F.S.M.G., M.d.S.D., J.B., A.T.C., R.B.M. and S.L. analyzed and discussed the results. All authors reviewed the manuscript.

Competing interests

The authors declare no competing financial interests.