The interfacial Dzyaloshinskii-Moriya Interaction (iDMI) is an antisymmetric exchange interaction that is induced by the broken inversion symmetry at the interface of, e.g., a ferromagnet/heavy metal. Thus, the presence of iDMI is not expected in symmetrical multilayer stacks of such structures. Here, we use thermal annealing to induce the iDMI in a [Py/Pt]×10 symmetrical multilayer stack. Brillouin light scattering spectroscopy is used to directly evidence the iDMI induction in the annealed sample. Structural characterizations highlight the modified crystallinity as well as a higher surface roughness of the sample after annealing. First principles electronic structure calculations demonstrate a monotonic increase of the iDMI with the interfacial disorder due to the interdiffusion of atoms, depicting the possible origin of the induced iDMI. The presented method can be used to tune the iDMI strength in symmetric multilayers, which are the integral part of racetrack memories, magnonic devices as well as spin-orbitronic elements.
Article
Inducing Dzyaloshinskii-Moriya interaction in symmetrical multilayers using post annealing
https://doi.org/10.21203/rs.3.rs-1528012/v1
This work is licensed under a CC BY 4.0 License
You are reading this latest preprint version
Dzyaloshinskii-Moriya Interaction
post annealing
interface
symmetrical structure
The Dzyaloshinskii-Moriya Interaction (DMI) is an antisymmetric exchange interaction that appears in inversion asymmetric structures and favors perpendicular alignment of neighboring spins in a magnetic material1. The DMI can lead to the formation of stable chiral spin textures such as Néel-type domain walls or skyrmions 2,3, which are promising candidates for memory and logic applications 4,5. In ferromagnet/heavy metal (FM/HM) heterostructures, the large spin–orbit coupling (SOC) of the HM and the broken inversion symmetry at the interface result in an interfacial Dzyaloshinskii–Moriya interaction (iDMI) 6.
For real applications e.g. racetrack memory devices, the DMI is demanded in multilayers or superlattices with both symmetric and asymmetric structures7,8 .However, the DMI vanishes in the symmetric structures due to the symmetry of the stacks, therefore, finding ways to induce the DMI in such structures is highly desired.
So far, several methods have been used for this purpose. For example, it has been reported that the iDMI is realized in such magnetic superlattices due to the gradience in the interface quality which is induced during the growth of the layers 2,3. Moreover, it has been shown that a finite iDMI can be observed if the upper and lower interfaces are of different crystallographic quality in a symmetric Pt/Co/Pt stack 9,10,11. This highlights the importance of the structure and the relative interface morphology of the upper and lower FM interfaces in determining the iDMI in such systems 12.
In this context, several studies have focused on the induction or tuning of the DMI in thin film structures via surface engineering. For example, Zimmermann et al. 13 used density-functional theory calculations to show the possibility to modify the magnetic properties at Co/Pt bilayer interfaces with chemical disorder based on the intermixing of atoms at the interface. Furthermore, the influence of the interface quality, considering the surface roughness and atomic intermixing on the iDMI has been studied in symmetric Pt/Co/Pt structure by Wells et al. 12. They found that the net iDMI increases as the difference between the interface quality of the top and bottom Co layer increases. In addition, Samardak et al. 14 studied the effect of atomic-scale surface modulation on the iDMI in ultrathin films composed of 5d heavy metal/ferromagnet/4d (5d) heavy metal or oxide interface. They predicted that a significant enhancement of the iDMI by up to 2.5 times can be achieved through interface roughness engineering on the atomic scale. Moreover, Torrejon et al. 15 have found that the strength of the iDMI varies by changing the 5d orbitals filling, in other words, the electronegativity of the heavy-metal layer that is used in the structure.
In this article, we use thermal annealing to induce the iDMI in a [HM/FM/HM]10 symmetrical multilayer stack. Structural characterizations show that the thermal annealing results in the interfacial microstructure evolution and roughness which causes the induction of the iDMI. Direct observation of the iDMI is performed by BLS spectroscopy measurement. To support the measurements, we carried out first-principles calculations which show a monotonic variation of the iDMI with the interfacial disorder due to interdiffusion of atoms. In order to study effects of the interface, we suggested substitutional disorder at the interface by switching the positions of Ni and Pt atoms and also the concentration of Fe atoms at the interface
A multilayer stack with Ta(5nm)/Pt(6nm)/[Ni80Fe20(3nm)/Pt(1nm)]×10/Au(3nm) structure was deposited on a thermally oxidized Si substrate at room temperature by DC magnetron sputtering. We first investigate the structural characteristics of the as-deposited and the annealed samples to understand the impact of the annealing on the microstructure of the stacks. Figure 1(a-b) show the atomic force microscopy (AFM) images of the as-deposited and annealed samples, respectively. The surface roughness is about 0.5 nm for the as-deposited while it is increased to approximately 5 nm for the annealed sample. The magnetic domain patterns which are measured by the magnetic force microscopy (MFM) for the as-deposited and annealed samples without an external field, are presented in Fig. 1(c-d), respectively. In the as-deposited sample, as demonstrated by the uniform color contrast, the magnetization is in the plane, while in the annealed sample, the magnetic domains formed in a labyrinthine pattern which can indicate that domains do not orientate purely in the plane 17. Field emission scanning electron microscopy (FESEM) was employed to characterize the structural and morphological properties. In Fig. 1(e-f) which display the FESEM images, one can observe that the as-deposited sample is uniform, however after annealing, the film surface is rough and planar nano-sized grains are formed. In general, from surface morphology characterizations, we conclude that a higher surface roughness can be observed for the annealed sample which can affect the iDMI consequently.
The X-ray pattern of the samples, presented in Fig. 1(g) indicate that Ni80Fe20 (Py) forms a face-centered-cubic structure of the type Ni3Fe (Reference code: 96-90-3942, Space Group: Pm-3m, 221) in both as-deposited and annealed samples. The intensity of Py and Pt peaks are more pronounced in the annealed samples. This indicates that the degree of crystallinity for the Py and Pt layers are appreciably improved after annealing. The improved crystallinity of the annealed sample is important for the presence of iDMI at the interfaces, which will be discussed later.
The hysteresis loops of the samples which are measured for the in-plane (IP) and out-of-plane (OP) configurations using a Vibrating-sample magnetometer (VSM) system are presented in Fig. 1(h). It can be seen that neither the as-deposited nor the annealed samples have a purely in-plane or an out-of-plane magnetic anisotropy, that highlights the reorientation of the magnetization. Moreover, the coercivity of the sample is increased after annealing which can be due to the higher crystallinity of this sample as well as the interface sharpening 18,19. In these multilayers, the sharp interfaces between the magnetic and non-magnetic lattices leads to the hybridization of the magnetic metal’s 3d and non-magnetic metal’s 5d orbitals, which can tune the direction of the magnetic anisotropy depending on the orbital bindings 20,21,22,23. In our sample, this anisotropy axis is more tilted towards the film plane after annealing, which is evidenced by the OP/IP relative saturation fields (Hs) observed from the VSM data. This relative increase is depicted at the inset of Fig. 1(h).
We now turn our focus to magnetic properties of the samples. First, the SW spectroscopy using BLS is used to measure the iDMI strength in both samples. An in-plane bias magnetic field sufficient to saturate the sample is applied to allow the SWs to propagate in the film plane. The iDMI leads to a frequency non-reciprocity of counter-propagating magnetostatic surface SWs (aka Damon-Eshbach mode). Such a frequency non-reciprocity can be observed by the difference between the Stokes (S) and anti-Stokes (AS) frequencies in the BLS spectroscopy 24,25. The same argument applies once the direction of the applied field is inverted. The measured frequency non-reciprocity with respect to the wave vector can be easily characterized, which depends on the DMI constant (D) based on the \(\varDelta\)f =\([2\gamma /{\pi M}_{s}\)] D ksw, where \(\gamma\) is the gyromagnetic ratio and MS is the saturation magnetization 25. It is required to probe surface waves, for which iDMI-induced nonreciprocity becomes maximal and, in particular, linear in k with DMI constant.
The BLS spectra of the SWs having wavenumber of ksw = 20.44 rad/µm in the presence of an applied field of \({\mu }_{0}\)H = ±50 mT (with opposite direction) are shown for the as-deposited and annealed samples in Fig. 2 (a-b) respectively. The BLS spectrum of both field directions shows a SW mode. The annealed sample shows a different SW characteristic. The results from the Lorentzian fits show that the frequency non-reciprocity for this wavenumber is as large as Δf = 490 MHz.
In order to understand the nature of the observed SW mode, we present the measured SW dispersion relation. Figure 2(c-d) shows frequency of the S and AS peaks as a function of the wave vector ksw, at \({\mu }_{0}\)H = ±50 mT for the AS and S modes for the as-deposited and annealed samples respectively. Indeed, the SW dispersion relation is reciprocal, i.e. the AS and S peaks for a given wavenumber have similar frequency for the as-deposited sample. From Fig. 2d, there is frequency difference between counter-propagating SWs for both orientations of the applied magnetic field for the annealed sample.
Figure 2(e) shows the frequency non-reciprocity (∆f) with respect to the wave vector for the as-deposited sample. A vanishing of the frequency difference is observed, indicating that the iDMI is absent in this sample. The sign and the amplitude of the iDMI constant D, have been obtained from the measured frequency non-reciprocity which is presented in Fig. 2(f) for the annealed sample 24,26. The results show a linear variation of the frequency non-reciprocity (∆f) with respect to the wave vector. The slopes obtained from a linear fit of this variation is 0.07881\(\pm\)0.011 GHz/radµm−1. From the slope of the curve, we estimate the value of the D by assuming \(\gamma\) =185 rad/T−1ns−1 27, 28 and MS = 700 kA/m for Py in our samples 29. Therefore, the estimated value of the iDMI constant from BLS measurements is D = 0.47\(\pm\) 0.02 mJ/m2. The obtained iDMI constant is positive (D \(>\) 0), thus, it favors a right-handed chirality 30.
As the structural characterizations show the crystal related evolution by annealing and also the BLS results show the evolution in iDMI, therefore, for a better evaluation on the possible origin of the induced iDMI, computational study were performed for different interface conditions as the following.
The iDMI is particularly sensitive to the structural changes of the layers, especially at the interfaces. This includes processes such as crystallization or interdiffusion at the interfaces. Therefore, the improved ordering of the atoms at the Pt/Py interface, due to the annealing treatment, can be the reason for the modification of the iDMI 9, 31. In order to support this hypothesis, we have carried out first-principle calculations, discussed in detail below.
The first principles calculations of DMI coefficients are carried out using relativistic exchange couplings 32 calculated from the Hamiltonian determined from density functional theory calculations employing the OpenMX 33, 34 ab initio package. We adopted Troullier-Martins type norm-conserving pseudopotentials 35 with partial core correction. We used a 24\(\times\)24\(\times\)1 k-point mesh for the first Brillouin zone (BZ) integration, and an energy cutoff of 500 Ry for numerical integrations in the real space grid. The localized orbitals were generated with radial cutoffs of 6.0 a.u, 6.0 a.u and 7.0 a.u. for Ni, Fe and Pt, respectively 33, 28.
Structural relaxations were carried using the Vienna ab initio simulation package (VASP) 36, 37 within the generalized gradient approximation as parameterized by Perdew et al 38. The pseudopotential and wave functions are treated within the projector-augmented wave method 39. A 5 
thick vacuum region is introduced to separate the periodic slabs along the stacking direction. The plane wave cutoff energy was set to 500 eV and a 6\(\times\)6\(\times\)1 k-points mesh was used in the 2D Brillouin zone sampling.
The crystal structure of the Ni3Fe (1 nm) /Pt (0.85 nm) bilayer device and the calculated DMI values are presented in Fig. 3. In Figs. 3(a) and (b) we show the L12 crystal structure of Ni3Fe and the FCC crystal structure of Pt grown along the \([\text{1,1},\stackrel{-}{1}]\) and \([3,\stackrel{-}{1},\stackrel{-}{1}]\) directions, respectively. We build up supercells by aligning the in-plane [101] and [011] directions of Ni3Fe with the \([0,\stackrel{-}{1},2]\) and [011] directions of Pt as shown in Fig. 3(a) and (b), respectively. This choice of growth direction results in a bilayer film with hexagonal crystal structure and a minimized lattice mismatch of 2% between the Ni3Fe and Pt films.
We have investigated the effect of two types of interfacial disorder on the iDMI. The first type corresponds to interdiffusion of Ni atoms in the Pt thin film, while the second corresponds to variation of the interfacial Fe concentration across the multilayer. In the first type of the disorder, shown in Figs. 3(c)-(e), we introduced substitutional interfacial disorder by switching the positions of Ni and Pt atoms at the interface. Figure 3 (c) shows the ideal interface, while Figs. 3(d) and (e) illustrate lower and higher levels of disordered interfaces with one and two interfacial Ni atoms switched with the interfacial Pt, respectively.
We find that the iDMI for the ideal, one and two substitutional atoms is 0.05 pJ/m, 0.7 pJ/m and 2.5 pJ/m, respectively, indicating a monotonic increase of the DMI with interfacial disorder, which is opposite to the case of Pt/Co bilayer devices13. It should be noted that, since the interfacial disorder breaks the in-plane isotropy, we consider the average values of the DMI, i.e., DMI= (\({D}_{xy}-{D}_{yx}\))/2.
The second type of interfacial disorder [Fig. 3(f)] in the [Ni3Fe/Pt] multilayer corresponds to substitution of interfacial Ni atoms with Fe atoms. Figure 3(f) displays the calculated DMI versus concentration of Fe atoms at the interface, which shows an almost linear decrease of DMI with the interfacial Fe concentration, resulting in a sign reversal of the DMI around 35% Fe concentration.
In conclusion, we have studied the magnetic properties of sputter-deposited [Py(3nm)/Pt(1nm)]×10 symmetrical multilayer stacks. Thermal annealing at 300°C and in inert atmosphere is used to induce interfacial Dzyaloshinskii–Moriya interaction in such symmetrical structures. Experiments based on Brillouin light scattering spectroscopy have been carried out to evidence the induction of the iDMI in the annealed samples. The estimated value of the induced iDMI in the annealed samples is approximately D = 0.47 mJ/m2. Structural characterizations showed that the annealing leads to a modified crystallinity and surface roughness, which can explain the underlying mechanism of the iDMI induction using this method. To further support our experimental observations, first principles calculations have been performed to show that increase of atomic interdiffusion at the interfaces of the magnetic/heavy metal layers enhances the iDMI from 0.05 pJ/m to 2.5 pJ/m. The presented method can be used to tune the asymmetric exchange interactions in multilayers, which is an essential ingredient for spintronic and magnonic applications such as memory devices as well as nonreciprocal spin-wave elements.
A multilayer stack with Ta(5nm)/Pt(6nm)/[Ni80Fe20(3nm)/Pt(1nm)]×10/Au(3nm) structure was deposited on a thermally oxidized Si substrate at room temperature by DC magnetron sputtering. The base chamber pressure was 10− 8 Torr, and samples were grown under an Ar pressure of 4 mTorr. The annealing of the sample has been carried out under Ar atmosphere (100 mTorr) at 300°C for 1 hour.
BLS measurements are performed using a λ = 532 nm laser beam by rotating the sample around a planar axis covering angles of incidence θ = 2.5°- 75°, corresponding to the spin-wave (SW) wavenumber of ksw = 1.03–22.8 rad/µm given by ksw = 4πsin (θ)/λ 16. The SWs frequency is determined by analyzing the backscattered light with a (3 + 3)-pass tandem Fabry-Pérot interferometer.
ACKNOWLEDGEMENT
Majid Mohseni acknowledges foundation from Iran Science Elites Federation (ISEF) and Iran National Science Foundation (INSF) under Grant No. 97001930. This project is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - TRR 173 - 268565370 (“Spin + X”, Project B01), the DFG Priority Programme “SPP2137 Skyrmionics”. The work at CSUN is supported by NSF ERC-Translational Applications of Nanoscale Multiferroic Systems (TANMS) Grant No. 1160504
AUTHOR CONTRIBUTIONS
K.A., L.J. and S.M.M. conceived and designed the study. S.S.M. fabricated the samples. K.A. carried out the experiments. K.A., M.M. and L.J. analyzed the data and interpreted the results. F.M. and N.K. performed the first-principles calculations and the analysis of relevant data. M.M., T.B. and P.P. carried out BLS analysis. J.A. commented on the manuscript. K.A. wrote the paper with L.J., M.M. and F.M. All authors discussed and contributed to the interpretation of the results, as well as to the writing of the manuscript.
ADDITIONAL INFORMATION
Competing Interests
The authors declare no competing interests.
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding authors upon reasonable request.
- Cao, A. et al. Enhanced interfacial Dzyaloshinskii—Moriya interactions in annealed Pt/Co/MgO structures. Nanotechnology 31, 155705 (2020).
- Bode, M. et al. Chiral magnetic order at surfaces driven by inversion asymmetry. Nature 447, 190–193 (2007).
- Rößler, U. K., Bogdanov, A. N. & Pfleiderer, C. Spontaneous skyrmion ground states in magnetic metals. Nature 442, 797–801 (2006).
- Parkin, S. & Yang, S. H. Memory on the racetrack. Nat. Nanotechnol. 10, 195–198 (2015).
- Tomasello, R. et al. A strategy for the design of skyrmion racetrack memories. Sci. Rep. 4, 1–7 (2014).
- Kim, D. H. et al. Bulk Dzyaloshinskii–Moriya interaction in amorphous ferrimagnetic alloys. Nat. Mater. 18, 685–690 (2019).
- Moreau-Luchaire, C. et al. Additive interfacial chiral interaction in multilayers for stabilization of small individual skyrmions at room temperature. Nat. Nanotechnol. 11, 444–448 (2016).
- Soumyanarayanan, A. et al. Tunable room-temperature magnetic skyrmions in Ir/Fe/Co/Pt multilayers. Nat. Mater. 16, 898–904 (2017).
- Hrabec, A. et al. Measuring and tailoring the Dzyaloshinskii-Moriya interaction in perpendicularly magnetized thin films. Phys. Rev. B 90, 020402 (2014).
- Je, S. G. et al. Asymmetric magnetic domain-wall motion by the Dzyaloshinskii-Moriya interaction. Phys. Rev. B - Condens. Matter Mater. Phys. 88, 1–5 (2013).
- Petit, D., Seem, P. R., Tillette, M., Mansell, R. & Cowburn, R. P. Two-dimensional control of field-driven magnetic bubble movement using Dzyaloshinskii-Moriya interactions. Appl. Phys. Lett. 106, 0–4 (2015).
- Wells, A. W. J., Shepley, P. M., Marrows, C. H. & Moore, T. A. Effect of interfacial intermixing on the Dzyaloshinskii-Moriya interaction in Pt/Co/Pt. Phys. Rev. B 95, 054428 (2017).
- Zimmermann, B. et al. Dzyaloshinskii-Moriya interaction at disordered interfaces from ab initio theory: Robustness against intermixing and tunability through dusting. Appl. Phys. Lett. 113, 5–10 (2018).
- Samardak, A. S. et al. Enhancement of perpendicular magnetic anisotropy and Dzyaloshinskii–Moriya interaction in thin ferromagnetic films by atomic-scale modulation of interfaces. NPG Asia Mater. 12, (2020).
- Torrejon, J. et al. Interface control of the magnetic chirality in CoFeB/MgO heterostructures with heavy-metal underlayers. Nat. Commun. 5, 4–11 (2014).
- Sebastian, T., Schultheiss, K., Obry, B., Hillebrands, B. & Schultheiss, H. Micro-focused Brillouin light scattering: imaging spin waves at the nanoscale. Front. Phys. 3, 1–23 (2015).
- Hellwig, O., Berger, A., Kortright, J. B. & Fullerton, E. E. Domain structure and magnetization reversal of antiferromagnetically coupled perpendicular anisotropy films. 319, 13–55 (2007).
- Cui, B. et al. Perpendicular magnetic anisotropy in CoFeB / X (X = MgO, Ta, W, Ti, and Pt) multilayers. J. Alloys Compd. 559, 112–115 (2013).
- Den Broeder, F. J. A., Kuiper, D., van de Mosselaer, A., P.and Hoving, W. Perpendicular Magnetic Anisotropy of Co-Au Multilayers Induced by Interface Sharpening. Phys. Rev. Lett. 60, 2769–2772 (1988).
- Daalderop, G.H.O., Kelly, P.J. and Schuurmans, M. F. H. Magnetic anisotropy of a free-standing Co monolayer and of multilayers which contain Co monolayers. Phys. Rev. B 50, (1994).
- Nakajima, N. et al. Perpendicular Magnetic Anisotropy Caused by Interfacial Hybridization via Enhanced Orbital Moment in Co Í Pt Multilayers: Magnetic Circular X-Ray Dichroism Study. Phys. Rev. B 81, 5229–5232 (1998).
- P. F. Carcia, A. D. Meinhaldt, and A. S. Perpendicular magnetic anisotropy in Pd / Co thin film layered structures. Appl. Phys. Lett. 47, 178–180 (1985).
- Jamilpanah, L., Hajiali, M. & Mohseni, S. M. Interfacial magnetic anisotropy in Py/MoS2 bilayer. J. Magn. Magn. Mater. 514, (2020).
- Di, K. et al. Direct observation of the Dzyaloshinskii-Moriya interaction in a Pt/Co/Ni film. Phys. Rev. Lett. 114, 1–15 (2015).
- Soucaille, R. et al. Probing the Dzyaloshinskii-Moriya interaction in CoFeB ultrathin films using domain wall creep and Brillouin light spectroscopy. Phys. Rev. B 94, 104431 (2016).
- Belmeguenai, M. et al. Interfacial Dzyaloshinskii-Moriya interaction in perpendicularly magnetized Pt/Co/AlOx ultrathin films measured by Brillouin light spectroscopy. Phys. Rev. B 91, 180405 (2015).
- Stashkevich, A. A. et al. Experimental study of spin-wave dispersion in Py/Pt film structures in the presence of an interface Dzyaloshinskii-Moriya interaction. Phys. Rev. B 91, 214409 (2015).
- Shaw, J. M., Nembach, H. T., Silva, T. J. & Boone, C. T. Precise determination of the spectroscopic g -factor by use of broadband ferromagnetic resonance spectroscopy. J. Appl. Phys. 114, 243906 (2013).
- Nembach, H. T., Shaw, J. M., Weiler, M., Jué, E. & Silva, T. J. Linear relation between Heisenberg exchange and interfacial Dzyaloshinskii-Moriya interaction in metal films. Nat. Phys. 11, 825–829 (2015).
- Rowan-Robinson, R. M. et al. The interfacial nature of proximity-induced magnetism and the Dzyaloshinskii-Moriya interaction at the Pt/Co interface. Sci. Rep. 7, 16835 (2017).
- Lavrijsen, R. et al. Asymmetric magnetic bubble expansion under in-plane field in Pt / Co / Pt: Effect of interface engineering. Phys. Rev. B 104414, 1–8 (2015).
- Mahfouzi, F. & Kioussis, N. First Principles Calculation of Dzyaloshinskii-Moriya Interaction: A Green’s function Approach. (2021).
- Ozaki, T. Variationally optimized atomic orbitals for large-scale electronic structures. Phys. Rev. B - Condens. Matter Mater. Phys. 67, 155108 (2003).
- Ozaki, T. & Kino, H. Efficient projector expansion for the ab initio LCAO method. Phys. Rev. B - Condens. Matter Mater. Phys. 72, 045121 (2005).
- Troullier, N. & Martins, J. L. Efficient pseudopotentials for plane-wave calculations. Phys. Rev. B 43, 1993–2006 (1991).
- Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B - Condens. Matter Mater. Phys. 54, 11169–11186 (1996).
- Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).
- Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).
- Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B - Condens. Matter Mater. Phys. 59, 1758–1775 (1999).
No competing interests reported.