Engineering of spin mixing conductance in Ru/FeCo/Ru interfaces: Effect of Re Doping
Rahul Gupta Nilamani Behera Vijay A. Venugopal Swaraj Basu Mark A. Gubbins Lars Bergqvist Rimantas Brucas Peter Svedlindh Ankit Kumar






frame=tb,
  	language=Matlab,
  	aboveskip=3mm,
  	belowskip=3mm,
  	showstringspaces=false,
  	columns=flexible,
  	basicstyle=,
  	numbers=none,
  	numberstyle=,
 	keywordstyle=,
	commentstyle=,
  	stringstyle=,
  	breaklines=true,
  	breakatwhitespace=true
  	tabsize=3








^1Department of Engineering Sciences, Uppsala University, Box 534, SE-751 21 Uppsala, Sweden
 
^2Seagate Technology, BT47 6SE, Londonderry, United Kingdom,

^3Department of Applied Physics, School of Engineering Sciences, KTH Royal Institute of Technology, Electrum 229, SE-16440 Kista,
 
^4Sweden Swedish e-Science Research Center, KTH Royal Institute of Technology, SE-10044 Stockholm, Sweden
 








§ ABSTRACT

We have deposited polycrystalline Re doped (Fe_65Co_35)_100-xRe_x (0 ≤ x ≤ 12.6 at%) thin films grown under identical conditions and sandwiched between thin layers of Ru in order to study the phenomenon of spin pumping as a function of Re concentration. In-plane and out-of-plane ferromagnetic resonance spectroscopy results show an enhancement of the Gilbert damping with an increase in Re doping. We found evidence of an increase in the real part of effective spin mixing conductance [Re(g^↑↓_eff)] with the increase in Re doping of 6.6 at%, while a decrease is evident at higher Re doping. The increase in Re(g^↑↓_eff) can be linked to the Re doping induced change of the interface electronic structure in the non-magnetic Ru layer and the effect interfacial spin-orbit coupling has on the effective spin-mixing conductance. The lowest and highest values of Re(g^↑↓_eff)  are found to be 9.883(02) nm^-2 and 19.697(02) nm^-2 for 0 at% and 6.6 at% Re doping, respectively. The saturation magnetization decreases with increasing Re doping, from 2.362(13) T for the undoped film to 1.740(03) T for 12.6 at% Re doping. This study opens a new direction of tuning the spin-mixing conductance in magnetic heterostructures by doping of the ferromagnetic layer, which is essential for the realization of energy efficient operation of spintronic devices.

^*corresponding author: ankit.kumar@angstrom.uu.se, peter.svedlindh@angstrom.uu.se




Engineering of spin mixing conductance in Ru/FeCo/Ru interfaces: Effect of Re Doping
    Rahul Gupta[itm:Rahul]^1, Nilamani Behera[itm:Rahul]^1, Vijay A. Venugopal[itm:VJ]^2, Swaraj Basu[itm:VJ]^2, Mark A. Gubbins[itm:VJ]^2, Lars Bergqvist[itm:Lars]^3,4, Rimantas Brucas[itm:Rahul]^1, Peter Svedlindh[itm:Rahul]^1[itm:corrs]^*, Ankit Kumar[itm:Rahul]^1[itm:corrs]^*

    December 30, 2023
=========================================================================================================================================================================================================================================================================================



§ INTRODUCTION

Generation, transportation and detection of pure spin currents play a fundamental role in spintronic devices. Pure spin currents can be generated in ferromagnet and normal metal (FM/NM) bilayer structures, which are governed by the magnetization dynamics via spin pumping [PhysRevLett.88.117601]. With the uniform ferromagnetic resonance (FMR) mode excited in the FM/NM structure, transfer of spin angular momentum occurs from the FM layer to the NM layer via the interface, due to magnetization precession. This phenomenon is called spin pumping, which can be measured by the FMR technique [PhysRevB.19.4382,PhysRevB.66.104413,PhysRevLett.97.216603,saitoh2006conversion]. Quantitatively, the spin pumping is determined by the spin mixing conductance (SMC) and damping (α) of the spin dynamics [PhysRevB.66.224403,brataas2004spin].

The generalized form of the pure spin current (a 2 × 2 matrix in spin space) is governed by the spin-dependent conductance, which depends on the reflection (r) and transmission (t) matrices. The SMC is defined as g^↑↓ =  e^2/h [M-Σ r^↑ (r^↓)^*], where M is the number of propagating quantum channels at the Fermi level and r^↑(↓) is the reflection matrix for spin-up (spin-down) electrons. The first part of this equation (e^2/h M) is called the Sharvin conductance, where M is proportional to the Fermi surface-averaged density of states at the interface [li2018effect]. The SMC is a complex parameter g^↑↓ = Re(g^↑↓) + ι 
Im(g^↑↓) [in units of Ω^-1 m^-2 (or per unit quantum conductance per unit area)], where Re(g^↑↓) and Im(g^↑↓) are the real and imaginary parts of the SMC, respectively [PhysRevB.66.224403]. However, for realistic interfaces, as per first principle calculations, the imaginary part of g^↑↓ is negligible for metallic systems in the absence of interfacial Rashba states [PhysRevB.65.220401]. This approximation has been considered throughout the paper. The SMC also varies between different combinations of FM/NM layers [gerrits2006enhanced,carva2007ab,PhysRevLett.107.046601,deorani2013role,weiler2014detection]. For example, g^↑↓ is found to be  about six times smaller for low conductivity FMs (e.g. Fe_3O_4) than for highly conductive FMs (3d transition metal alloys; e.g. permalloy and Heusler compounds) with Pt as NM layer [PhysRevLett.107.046601]. Moreover, M. Schoen  et. al. [schoen2016ultra] (in supplementary information) reported the compositional dependence of the SMC in iron-cobalt (FeCo) alloy thin films. Studies of FeCo have led to many advances in the fundamental understanding of magnetic devices because of its metallic nature and ultra-low damping [schoen2016ultra], providing a low threshold current density in spin-transfer-torque magnetic-random-access-memories. However, interfacial properties of FeCo based heterostructures are still to be investigated in detail. In addition, engineering the properties of FM layer and its interface with NM layer are required to achieve high efficient operation of spin-logic devices, e.g., spin torque MRAM [PhysRevApplied.5.014002]. The next generation spin torque MRAM for high-density catch memory requires high magnetic anisotropy for thermal stability and high damping for ultrafast operation [6374706,bhatti2017spintronics]. The benchmark for efficient operation of spin-logic devices is the product of the switching power and the delay, where the delay represents the switching time. The power-delay product describes the average energy consumed per switching event, and hence balances the performance and energy consumption. However, there is trade-off between ultrafast switching and energy efficient operations as an increase in damping enhances the writing energy. A reduction of the power-delay can be achieved by enhancing the SMC of the FM/NM interface.
The energy-delay reduction scaling can be achieved by enhancing the SMC of the FM/NM layers interface. The SMC governs the efficiency of spin angular momentum transfer across the interface, providing spin torques to switch the magnetic state of the FM layer in MRAMs. Manipatruni et al. have reported that by increasing the SMC at the FM/NM interface, the power-delay product of in-plane all spin-logic devices can be reduced significantly [PhysRevApplied.5.014002].

Doping the FM layer with 5d transition elements can be used to tune the magnetic properties [PhysRevB.87.014430], providing a way of optimizing  magnetic heterostructures for efficient spin logic devices. For example, it is expected that damping relaxation due to two-magnon-scattering (TMS) increases with non-magnetic doping. Moreover, other contributions to the damping, such as radiative and eddy current contributions, owing to the inductive coupling between the sample and the co-planar waveguide (CPW), also have to be considered when estimating the SMC of the heterostructures; the extrinsic contributions have to be subtracted from the total measured damping in order to reliably extract the interfacial properties.

To tune the interfacial and dynamic properties, as per theory and experimental suggestions [schoen2016ultra], we used low damping Fe_65Co_35 thin films with Ru as seed and capping layers, as Ru can help to decrease the coercivity and the effective damping parameter of Fe_65Co_35 thin films [akansel2018effect]. In the present work, by 6.6 at% Re-doping of technologically important Fe_65Co_35 interfaced with Ru layers, we have evidenced a 100% increase of the SMC along with 135% increase of the damping, while the saturation magnetization decreases marginally. These properties make Re-doped Fe_65Co_35 thin films interesting for ultrafast and energy efficient spin torque logic devices.





§ EXPERIMENTAL DETAILS

Using DC magnetron sputtering, rhenium (Re) doped polycrystalline Fe_65Co_35 alloy thin films with different thickness were deposited at room temperature on Si(100)/SiO_2 substrates with Ru as seed and capping layers. The Re concentration (x) in the (Fe_65Co_35)_100-xRe_x (0 ≤ x ≤ 12.6 at%) thin films was varied by changing the deposition rate of Re. The target power, base pressure and Ar working pressure in the vacuum chamber were the same during all the depositions; for further details see Ref. [akansel2019enhanced]. The nominal thicknesses of the FeCo films was 5, 10, 15, 20 and 30 nm and the nominal thickness of the Ru layer was 3 nm. The sample configuration is shown in the inset of Fig. [fig:fig1]1.

Our polycrystalline films have the body-centered cubic crystal structure. The lattice constant increases with increasing Re doping, which follows the trend from theoretical calculations [akansel2019enhanced]. Rutherford backscattering spectrometry was used to discover the absolute elemental composition of Co, Fe, Ru and Re in our films. The detailed analysis can be found in our previous work [akansel2019enhanced].

Film thickness, surface/interface roughness, and density were measured by X-ray reflectivity (XRR) using the X’Pert-Pro system. The scan angle (2θ) was set to 0-8^∘ during the measurements. Fittings were done in specular mode using the HighScore software. The XRR spectra of the 20 nm thick of FeCo films for all Re concentrations are shown in Fig. [fig:fig1]1. Resistivity measurements were performed in a Quantum Design Physical Property Measurement System (PPMS) using the four-probe method.

The spin dynamic properties were studied using IP- and OP-FMR spectroscopy. The direction of the external magnetic field was parallel and perpendicular to the film plane in IP- and OP-FMR geometry, respectively, as shown in the inset of Fig. [fig:fig1]1. The IP- and OP-FMR measurements were based on lock-in amplifier and vector-network-analyzer techniques, respectively. The samples (4 × 4 mm^2 in size) were mounted face-down on the CPW and a fixed frequency (f) and amplitude microwave field was passed through the CPW using a microwave source [PhysRevB.96.224425]. The measurements were performed in field-sweep mode; the uniform FMR resonance mode was observed in the presence of the magnetic field, fulfilling the resonant condition for the material. The resonance field will hereafter be referred to as H_r and the linewidth of the resonance will be denoted Δ H.










§  RESULT AND DISCUSSION

Results from XRR measurements provide accurate estimation of thickness, roughness and density of the different layers in the SiO_2/Ru/FeCo(t_FeCo)/Ru/Oxide samples. Fig. [fig:fig1]1 shows observed (open symbols) along with simulated spectra (red lines) for FeCo (20 nm) with different Re concentrations; the results from the analysis are presented in Table [tab:table1]I. The model used for precise fitting, Ru/(Fe_65Co_35)_100-xRe_x/Ru/Oxide, includes 4 layers, where the FeCo layer can have different Re concentration. We use the XRR determined thicknesses of the FeCo layers for further calculations, represented by t_FeCo^E.




 §.§ IP-FMR




To understand the fundamental behavior of the magnetization dynamics, IP-FMR spectra were recorded using a signal generator–detector diode based lock-in technique. The microwave frequency range used in these measurements was from 8 to 20 GHz. Fig. [fig:fig2]2 shows typical IP-FMR spectra of a FeCo(20 nm) sample at different frequencies. To find the resonance field and linewidth (half width at half maximum, HWHM), the observed IP-FMR spectra were fitted with the field derivative of the susceptibility function [χ(H,t)] as given below [woltersdorf2004spin],



    dχ/dH∼ -S[Δ H^2-(H-H_r)^2]/[Δ H^2+(H-H_r)^2]^2 -A2(H-H_r)Δ H/[Δ H^2+(H-H_r)^2]^2+Dt,



where the first and second terms on the right-hand side represent the symmetric and the asymmetric contributions to the spectra, respectively. The asymmetry arises because of the electromagnetic coupling between the experimental setup and the low frequency, low amplitude magnetic modulation field. A and S are the asymmetric and symmetric coefficients, H is the applied magnetic field and Dt is linear drift term.



According to the Landau-Lifshitz-Gilbert model, the linewidth is directly proportional to microwave frequency, and the proportionality constant gives the intrinsic damping parameter (α) in the absence of inhomogeneity. However, in reality, the linewidth exhibits different contributions, e.g. from sample inhomogeneity, spin pumping, radiative, eddy-current, and TMS. As per Arias and Mills model [PhysRevB.60.7395], TMS can be caused by lattice geometrical defects and surface-induced inversion symmetry breaking. The spin pumping contribution to the linewidth is due to the adjacent NM layer, and radiative and eddy-current contributions are due to the inductive coupling between the sample and the CPW. These other contributions will be discussed separately in the next sections. Therefore, in IP-FMR measurements, the total linewidth is a function of all contributions (i.e.  μ_0 Δ H = μ_0 Δ H_0 + μ_0 Δ H^(G+SP) + μ_0 Δ H^rad + μ_0 Δ H^eddy + μ_0 Δ H^TMS, and is defined as,


    ΔH = hα_eff^∥/g^⊥μ_Bμ_0f + ΔH_0,


where h is the Planck's constant, g^⊥ is the Landé spectroscopic splitting factor, μ_B is the Bohr magneton, ΔH_0 is the frequency independent inhomogeneity contribution to the linewidth, and α_eff^∥ = [α^(G+SP) + α^rad + α^eddy + α^TMS] is known as the effective damping for the IP-FMR geometry, where α^G, α^SP, α^rad, α^eddy and α^TMS correspond to Gilbert damping, spin-pumping damping, radiative damping, eddy-current damping and two-magnon-scattering damping, respectively. The total linewidth vs. f is shown in Fig. [fig:fig3]3 for different Re concentrations and different FeCo thickness, and the thickness dependent effective damping is shown by red symbols in Fig. [fig:fig6]6. The extracted μ_0 Δ H_0  values from IP-FMR vary in the range 2.38-4.29 mT for t_FeCo^N = 5nm for all Re concentration, being smaller for t_FeCo^N > 5nm; 0.51-1.20 mT.

The easy axis IP-FMR Kittel equation, given below, is used to fit the experimental data [PhysRev.73.155]. The fits, shown in the insets of Fig. [fig:fig3]3, give the effective magnetization (M_eff) and the IP uniaxial anisotropic field (H_u),


    f=g^⊥μ_Bμ_o/h√((H_r+H_u)(H_r+H_u+M_eff)),



In Eqs. [eq:eq2]2 & [eq:eq3]3, the average value of g^⊥ = 2.07 extracted from the OP-FMR results using asymptotic analysis of the FMR data [shaw2013precise], was used.



 §.§ OP-FMR





Magnon scattering is the process where transitions between the uniform spin wave mode (k = 0) and degenerate spin wave modes (k ≠ 0) occur in the sample. According to Hurben and Patton [hurben1998theory], this process is strongly dependent on the wavelength of spin waves, and the coupling between k = 0 with k ≠ 0 modes are called two-magnon scattering, three-magnon scattering and so on. The level of three and higher order magnon scattering is expected to be much less than for two-magnon scattering. Therefore, it is assumed that only TMS is present in the IP-FMR results and that this contribution is minimum for the OP-FMR geometry [hurben1998theory]. Broadband OP-FMR measurements were performed where the complex transmission parameter (S_21) was recorded in field-sweep mode using VNA-technique. The used frequency range was 14-30 GHz. A typical result of the real and imaginary parts of the complex S_21 parameter is shown in Fig. [fig:fig4]4. To analyze the data, we used the fitting model by Nembach et.al. [PhysRevB.84.054424]. The complex S_21  function is used as follows,


    S_21(H,t) = S_21^0 + Dt - χ(H)/χ̃_̃õ,



    χ(H) = M_eff(H-M_eff)/(H-M_eff)^2 -H_eff^2 -ιΔH/2 (H-M_eff),


where S_21^0  is the nonmagnetic contribution to S_21 [ding2004coplanar], Dt is the first order drift correction, χ̃_̃õ is an imaginary function of frequency, film thickness and effective length of the sample on the CPW, μ_0 H_eff = hf/(g^⊥μ_B), and μ_0 Δ H is the HWHM linewidth, free from the TMS contribution. To find the precise μ_0 M_eff and g^⊥ values, considering both as fitting parameters, the OP Kittel equation can be used,


    f = g^⊥μ_B μ_0/h (H_r-M_eff),


where M_eff = M_s- H_k^⊥, μ_0 H_k^⊥ = (2k^⊥)/(M_s t_FeCo^E) is the OP anisotropy field and μ_0 M_s is the saturation magnetization, which can be found from extrapolation of the μ_0 M_eff vs. 1/t_FeCo^E curve (not shown here) to 1/t_FeCo^E = 0. The fitted μ_0 M_s values are included in Table [tab:table3]III for all Re concentrations, and the extracted μ_0 M_eff  and g^⊥ values are given in Table [tab:table1]I.





The effective damping (α_eff^⊥) can be extracted using,


    ΔH = hα_eff^⊥/g^⊥μ_Bμ_0f + ΔH_0,

where α_eff^⊥ = α^(G+SP) + α^rad + α^eddy is known as the OP effective damping parameter, free from the TMS contribution. The μ_0 Δ H_0 contribution to the linewidth is found to be less than 0.58 (02) mT for all samples except for 5 nm FeCo with the highest Re concentration, for which the value increased to 2.58 (07) mT. The μ_0 Δ H vs. f results are shown in Fig. [fig:fig5]5 for all Re concentrations and all FeCo thicknesses. A comparison of the effective damping, with (α_eff^∥) and without TMS contribution (α_eff^⊥), is presented in Fig. [fig:fig6]6 and in Table [tab:table2]II.

In our studied films, the thickness of the Ru layer, t_Ru^N (=3 nm)  is less than the spin diffusion length λ_Ru^sd (=15 nm) [behera2015effect]. Therefore, significant back flow of the spin angular momentum from the Ru layer to the FeCo layer is expected. The ratio of the IP to OP effective damping (ζ) can be calculated using the effective damping expression derived in Ref. [PhysRevLett.114.126602],  


    ζ = α_eff^∥/α_eff^⊥∝1+6ηξ/1+ξ+η/2(1+ξ)^2/1+4ηξ/1+ξ,


where η = (α_R k_F/E_F)^2 is an interface parameter, which depends on the Rashba coefficient (α_R), Fermi wave vector (k_F) and Fermi energy (E_F), and ξ is the back flow factor. In the case of complete back flow (ξ → ∞) ζ → 1.5, while for no back flow (ξ = 0) ζ → 1.0. According to this model, ζ values in the range 1-1.5 are expected, while our experimentally determined ζ values for the Re-doped FeCo films are found to be greater than 1.5, which infers that α_eff^∥ contains other contributions. The most common contribution to the effective damping enhancement in IP-FMR measurements is the TMS contribution. Therefore, in order to calculate the true value of the effective SMC, OP-FMR data devoid of TMS contributions can be used. Furthermore, extrinsic contributions, i.e. radiative and eddy contributions, to the OP effective damping must be subtracted before determining Re-doping dependent changes in the effective SMC.



 §.§ Radiative and Eddy contributions

According to Faraday's law, an AC voltage is generated in a metallic magnetic material when it is exposed to a time varying magnetic flux created by the microwave magnetic field and the precessing magnetization, which induces eddy-currents in the material and in the CPW. Therefore, damping due to eddy-currents in the sample is called eddy-current damping (α^eddy), whereas the damping in the CPW is referred to as radiative damping (α^rad). Both contributions to the effective damping depend on the wave-guide dimensions and the magnetic properties of the sample.
The eddy-current damping for the lowest order mode in FMR can be defined as [schoen2016ultra,PhysRevB.92.184417],


    α_eddy = C/16g^⊥μ_B μ_0^2 M_s (t^E)^2/ħρ,


where C describes the eddy-current profile. In our case, we used C = 0.5 for all samples and t^E is the total thickness of the full stack determined from XRR analysis. The α_eddy values are found to be (0.16-3.32) × 10^-5, (0.13-2.81) × 10^-5, (0.09-0.67) × 10^-5 and (0.08-0.45) × 10^-5 for the 0 at%, 3 at%, 6.6 at% and 12.6 at% Re doped samples, respectively. The first-last values of α_eddy within parentheses correspond to the smallest-largest thicknesses. The Landé g^⊥ factor and μ_0 M_s from the OP-FMR fitting were used in the calculations. ρ is the resistivity of the full stack sample, which takes values in the range (1.243-0.611) Ω-μm, (2.094-2.432) Ω-μm, (1.734-2.290) Ω-μm and (1.743-2.915) Ω-μm for the 0 at%, 3 at%, 6.6 at% and 12.6 at% Re doped samples, respectively. The first-last values of ρ within parentheses correspond to the smallest-largest thicknesses. α_eddy decreases with increasing Re doping, while α_eddy increases with an increase of the FeCo thickness.

Apart from eddy current damping, the radiative contribution to the damping also depends on the magnetic properties of the sample and the CPW parameters. The radiative damping is defined as [schoen2016ultra,PhysRevB.92.184417], 


    α_rad = ηg^⊥μ_B μ_0^2 M_s t^El/2ħ Z_0 W,


where η is a dimensionless parameter, which describes the actual FMR mode profile in the sample. In our case, we assume η = 0.25, considering the uniform mode of the FMR. μ_0 M_s and g^⊥ are used as obtained from the OP-FMR fitting. Z_0 = 50 Ω, W = 300 μm, t^E and l are the CPW impedance, central width of the CPW, total thickness and length of the sample on the CPW, respectively. The α_rad values (smallest - largest) are found to be (2.99-9.24) × 10^-4, (2.61-8.29) × 10^-4, (2.37-7.48) × 10^-4 and (2.1-6.42) × 10^-4 for the 0 at%, 3 at%, 6.6 at% and 12.6 at% Re doped samples, respectively. α_rad decreases with an increase of Re doping, while α_rad increases with an increase of FeCo thickness.







 §.§ Spin pumping (SP) and Gilbert contribution

The precessing magnetization transfers spin angular momentum to the NM layer in a FM/NM bilayer, which causes enhancement of the damping. According to the spin pumping model [PhysRevB.66.224403], the effective damping is governed by the real part of effective SMC and is expressed as,


    α_eff(α_G, α_SP)=α_G+μ_B/4πRe(g^↑↓_eff)/ M_sg^⊥/t_FeCo^E,


where α_eff is the effective damping without TMS, radiative  and eddy current contributions to the total damping, and hence only a function of the Gilbert and spin pumping contributions. The real part of the effective SMC for both the interface is obtained from thickness-dependent experimental effective damping, by linear fitting of α_eff vs. g^⊥/t_FeCo^E. The slope gives Re(g^↑↓_eff) and the interception with the y-axis gives the Gilbert damping parameter. α_G is found to be 4.5 times higher with an increase of Re doping from 0 at% to 12.6 at%. The α_G, Re(g^↑↓_eff) and μ_0 M_s values are presented in Table [tab:table3]III, and α_eff vs. g^⊥/t_FeCo^E plots are presented in Fig. [fig:fig6]6 together with fits according to Eq. [eq:eq11]11. 
As shown in Fig. [fig:fig7]7, Re(g^↑↓_eff) (per unit quantum conductance per unit area) exhibits non-monotonic behavior as a function of Re doping. The observed lowest and highest values of Re(g^↑↓_eff) are found to be 9.883(02) nm^-2 and 19.697(02) nm^-2 for Re 0 at% and Re 6.6 at% doping, respectively. A comparison with previously reported Re(g^↑↓_eff) values is presented in Table [tab:table4]IV. 

According to existing theory, Re(g^↑↓_eff) can be enhanced either by increasing the intrinsic SMC or by decreasing interfacial spin-memory-loss (SML) [PhysRevLett.114.126602]. A decrease of the SML with an increase of Re doping is unlikely; rather the opposite effect can be expected. The observed increase in Re(g^↑↓_eff) with the increase of Re doping from 0 at% to 6 at% should therefore according to existing theory be linked to the increase of the intrinsic SMC. However, this would require an increase of the number of quantum mechanical conductance channels in the NM layer, which again seems unlikely considering that the Re doping takes place in the FM layer. As we will argue below, the explanation for the increase of Re(g^↑↓_eff) can instead be found in the influence of the spin-orbit coupling (SOC) on the effective SMC. 

The intrinsic SMC calculated using the existing theory [PhysRevB.66.224403,brataas2004spin] is only well defined without SOC. It is nevertheless possible from theory to identify an effective SMC [Re(g^↑↓_eff)] by calculating the interface-enhanced damping parameter. Theoretical calculations for Pt/FeCo varying the concentration of Co show that the magnitude of Re(g^↑↓_eff) exhibits a non-monotonic behaviour with Co-concentration and that Re(g^↑↓_eff) can be increased by a factor of two in the presence of interfacial SOC on the NM side [PhysRevB.98.174412]. This change in Re(g^↑↓_eff) results from the partially filled d-bands and half-filled s-band of Pt (5d^96s^1). The interaction between the d-bands of FeCo and Pt near the Fermi level results in the non-monotonic Co-concentration dependent behaviour of Re(g^↑↓_eff). Similar calculations for Au(5d^106s^1)/FeCo show that  Re(g^↑↓_eff) is less influenced by SOC due to the mainly s-electron character of the conduction channel in Au. In addition, calculations of the Re(g^↑↓_eff) for Pt/Py indicated an increase by 25% in the presence of interfacial SOC, while no significant change was observed for Cu(3d^104s^1)/Py [PhysRevLett.113.207202]. For comparison, the experimentally observed enhancement in Re(g^↑↓_eff) for Ru(4d^75s^1)/FeCo [cf. Fig. [fig:fig7]7] can also be attributed to interfacial SOC. As the Re-concentration increases, the density of states at the Fermi level increases, which will also modify the interaction between the d-bands of FeCo and Ru and the interfacial SOC in the Ru layer. The interpretation for low Re-concentrations is that the interfacial SOC in the Ru layer increases with increasing Re-concentration. This interpretation might also explain the non-monotonic behaviour of Re(g^↑↓_eff) with increasing Re-concentration. However, the decrease of Re(g^↑↓_eff) at higher Re doping may be due to SML, since an increase of the SML for large enough Re-concentration cannot be ruled out [PhysRevLett.114.126602].




§ CONCLUSION

The dynamics properties of Re-doped (0-12.6 at%) polycrystalline Fe_65Co_35 thin films with Ru as capping and seed layers have been investigated using room temperature IP- and OP-FMR measurements. Comparison of the effective damping in IP- and OP-FMR indicates that the IP damping is affected by TMS, which overestimates the real part of the effective spin mixing conductance. The thickness dependent results for the effective damping, after subtracting the radiative and the eddy-current damping contributions, indicate that the enhancement of the damping is due to the spin pumping contribution. By further analyses of the OP-FMR results, free from TMS, radiative and eddy-current contributions, a non-monotonic dependence of the real part of the effective spin-mixing conductance on Re concentration is found. The increase of Re(g^↑↓_eff) with increasing doping from 0 at% to 6 at% Re is tentatively explained by a corresponding increase of the interfacial SOC in the Ru layer. Apart from this, an enhancement of the bulk Gilbert damping is found with an increase in Re doping, while μ_0 M_s decreases with increasing Re concentration. This study opens an entirely new direction of tuning the spin-mixing conductance in magnetic heterostructures by doping the ferromagnetic layer, thus providing a method for optimizing the design of spintronic devices.



§ ACKNOWLEDGEMENT

This work is supported by the Swedish Research Council (VR; grant no 2017-03799), Olle Engkvist Byggmästare (project number 182-0365), and by the Marie Curie Action ‘Industry-Academia Partnership and Pathways’ (ref. 612170, FP7-PEOPLE-2013-IAPP).

