Thank you for visiting nature.com. You are using a browser version with
limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off
compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site
without styles and JavaScript.

Subjects

Abstract

Chiral spin textures of a ferromagnetic layer in contact to a heavy non-magnetic metal, such as Néel-type domain walls and skyrmions, have been studied intensively because of their potential for future nanomagnetic devices. The Dyzaloshinskii–Moriya interaction (DMI) is an essential phenomenon for the formation of such chiral spin textures. In spite of recent theoretical progress aiming at understanding the microscopic origin of the DMI, an experimental investigation unravelling the physics at stake is still required. Here we experimentally demonstrate the close correlation of the DMI with the anisotropy of the orbital magnetic moment and with the magnetic dipole moment of the ferromagnetic metal in addition to Heisenberg exchange. The density functional theory and the tight-binding model calculations reveal that inversion symmetry breaking with spin–orbit coupling gives rise to the orbital-related correlation. Our study provides the experimental connection between the orbital physics and the spin–orbit-related phenomena, such as DMI.

Introduction

Chiral interaction between two atomic spins owing to a strong spin–orbit coupling (SOC), which is known as the Dzyaloshinskii–Moriya interaction (DMI), has attracted intense interest1,2. In particular, it has been demonstrated that the DMI at the interface between ferromagnetic (FM) and non-magnetic heavy metals (HMs) plays a major role for the formation of chiral spin textures, such as skyrmions3,4 and homochiral Néel-type domain walls (DWs)5,6,7, which are attractive for the development of future information storage technology8. Understanding the microscopic origin of the DMI is indispensable for the realization of such chiral spin textures9,10. It has been reported that the proximity-induced magnetic moment in HM layers is critical to promote the DMI11. However, this proximity effect is still controversial because it has been also reported that the induced magnetic moment has no direct correlation with the DMI in the case of the Co/Pt system12,13. The scattering of spin-polarized electrons on spin–orbit coupled impurities is known to give rise to the DMI in spin-glass systems as a microscopic viewpoint14,15. However, the orbital hybridization between the spin–orbit coupled ions and the magnetic matrix was explicitly neglected in the theory, which becomes problematic when considering transition metal interfaces. On the other hand, theories have predicted that SOC combined with inversion symmetry breaking (ISB) naturally introduces a chirality to conduction electron spins in equilibrium and the interfacial DMI at an FM/HM interface is related to this spin chirality16,17,18,19,20. It has also been reported that the spin chirality is a manifestation of the chirality of the orbital magnetism in strongly spin–orbit coupled systems with ISB21,22. These previous studies suggest a possible microscopic origin of the interfacial DMI, which has remained experimentally unaddressed so far.

Here we discuss the microscopic origin of the interfacial DMI with experimental and theoretical studies as follows: First, we show the temperature dependence of the DMI for a Pt/Co/MgO trilayer, which is one of the standard structures used for the studies of the DMI7,12,23, using the extended droplet model24. We find that the DMI increases with decreasing temperature in a range from 300 to 100 K. In general, the electron–phonon interaction promotes thermally induced hopping between nearest neighbours when increasing the temperature25. As a result, it is expected that the difference between in-plane and out-of-plane hopping energies is reduced upon temperature increase. Therefore, changing the temperature of the system allows for charge redistribution between in-plane and out-of-plane orbitals while preserving the integrity of electronic states of the trilayer unlike other interface control methods such as ion irradiation26 or thermal annealing technique27, which may cause a permanent atomic rearrangement and thus induce undesired extrinsic effects. To discuss this temperature dependence of the DMI, that of the spin (ms) and orbital (mo) magnetic moments of Co and Pt is studied by X-ray magnetic circular dichroism (XMCD) spectroscopy. We find that ms values of Co and Pt show temperature dependences due to change in Heisenberg exchange. Furthermore, the intra-atomic magnetic dipole moment (mD), which is due to the asymmetric spin density distribution (fSD)28,29, displays strong temperature dependence, suggesting a sizable modification of the charge distribution between the in-plane and the out-of-plane d-orbitals under temperature variation. We also find that the out-of-plane orbital moment (\({\boldsymbol{m}}_{\mathrm{o}}^ \bot\)) shows large temperature dependence while in-plane orbital moment (\({\boldsymbol{m}}_{\mathrm{o}}^\parallel\)) does not, revealing a close connection between the anisotropy of m0 (orbital anisotropy) and the DMI. The ab initio and the tight-binding model calculations suggest that the ISB-dependent electron hopping, which gives rise to the asymmetric charge distribution at the interface of the FM/HM, is a possible microscopic origin of the correlation between the orbital anisotropy and the DMI.

Results

Temperature dependence of the DMI in the Pt/Co/MgO trilayer

The temperature dependence of the DMI-induced effective field (HDMI) of the Pt (2 nm)/Co (0.5 nm)/MgO (2 nm) trilayer is determined by measuring the nucleation field (Hn) applied along the in-plane (Hx) and out-of-plane (Hz) direction as schematically shown in Fig. 1a. The measured Hn values are analysed with the extended droplet model24. Since Hn is proportional to square of DW energy (\({\mathrm{\sigma }}_{{\mathrm{DW}}}^2\)), the σDW is a crucial parameter for the nucleation24,30. In case of the droplet with sizable DMI, DW magnetizations are aligned in the radial direction. When we consider the two DW magnetizations with respect to Hx, DW1 and DW2 (see the inset of Fig. 1b), \({\mathrm{\sigma }}_{{\mathrm{DW}}}^2\) of the DW1 and the DW2, respectively, follows the blue and the red curves in terms of Hx due to DMI as illustrated in Fig. 1b31. Assuming that total DW energy mainly depends on the two magnetizations at the DW of a droplet with respect to Hx, the total \({\mathrm{\sigma }}_{{\mathrm{DW}}}^2\) follows the trend of the purple-dotted line. Note that there is a threshold where the \({\mathrm{\sigma }}_{{\mathrm{DW}}}^2\) starts to become Hx dependent. Our previous report showed that the extended droplet model allows us to estimate the DMI-induced field from the threshold field24. In spite of such one-dimensional approximation, we have demonstrated that the fitting gives reliable value of the DMI energy density (D).

Fig. 1

Droplet nucleation with DMI. a Schematic image of a magnetic droplet in a ferromagnetic medium under Hz and Hx. Inset shows a magneto-optical Kerr effect (MOKE) image to prove the droplet nucleation inside the Pt/Co/MgO microstrip. The white bar is a scale bar of 5 μm. b Schematic diagrams of \({\mathrm{\sigma }}_{{\mathrm{DW}}}^2\) in terms of Hx. Inset shows the magnetization alignments of the DW1 and DW2 in the DW of the droplet. The thresholds in the curve of averaged \({\mathrm{\sigma }}_{{\mathrm{DW}}}^2\) are highlighted with blue shades

The measured Hn as a function of Hx shows a threshold arising from the DMI as predicted by the extended droplet model24 (see also Supplementary Note 1, 2 and Method section for details about the analysis and the measurement). Figures 2a–d show Hn/HSW(Hx=0) as a function of Hx/HK at various temperatures (a, T=300 K; b, 200 K; c, 150 K and d, 100 K). Here Hn and Hx are normalized by the switching field at Hx=0 [HSW(Hx=0)] and by the anisotropy field (HK), respectively. Those normalized values allow us to clearly confirm the DMI-dependent threshold with ruling out the temperature-dependent characteristics of HK and the HSW(Hx=0). The temperature-dependent HDMI can be determined from the best fitting using the extended droplet model; HDMI at 300, 200, 150 and 100 K are 166±50, 245±45, 324±15 and 372±30 mT, respectively (Fig. 2e). The temperature-dependent D is readily calculated from HDMI and Δvia \({\boldsymbol{H}}_{{\mathrm{DMI}}} = {\boldsymbol{D}}/\mu _0M_{\mathrm{s}}\varDelta\), where μ0 is the permeability, Ms is the saturation magnetization and Δ is the DW width6. Δ and Ms values in terms of the temperature are listed in Table 1. We find that D has a strong temperature dependence as shown in Fig. 2e; D increases by a factor of 2.2 as the temperature decreases from 300 to 100 K.

Fig. 2

HDMI measurement from the Hn of the droplet. a–d The Hn/Hsw(Hx=0) vs Hx/HK plots measured at 300, 200, 150 and 100 K, respectively. The grey solid lines are the best fitting results using the droplet model. Here the vertical axis is normalized by the nucleation field Hsw(Hx=0) at Hx=0 and the horizontal axis is normalized by the effective perpendicular anisotropy field HK. e Plots of HDMI and D in terms of T. The error bars are based on the standard deviation of the HDMI distribution from the best fitting

Heisenberg exchange is one of key parameters to understand the microscopic origin of DMI2,32. From the temperature dependence of Ms, the exchange stiffness constant (A) values were quantitatively obtained as detailed in Supplementary Note 3. As listed in Table 1, we find a clear correlation of the Heisenberg exchange with DMI; the A value increases by 57% when the temperature decreases from 300K to 100 K while DMI shows 100% increase. Based on Moriya’s theory, DMI requires three ingredients: (i) SOC, (ii) magnetic exchange, and (iii) ISB. Therefore, we can consider two parameters in addition to Heisenberg exchange: one is the proximity-induced magnetism of the non-magnetic element, which is related to the source ii, and the other is the hybridization between 3d and 5d orbitals, which is related to the sources i and iii. In order to find how DMI correlated with the induced moment in the Pt layer and the orbital structure of Co, XMCD studies were performed using the soft and hard X-ray as discussed in the following sections.

Temperature dependence of the proximity-induced magnetic moment of Pt

In this section, we first investigate the role of the induced magnetic moment in the Pt layer as mentioned in previous section. The temperature dependence of the Pt-induced magnetic moment was measured using the XMCD method, which enables element-specific analyses of spin and orbital magnetism33,34. Figure 3a presents the XMCD and integration of XMCD spectra measured at the Pt L2,3 edge. The intensities of XMCD are ~3% of the X-ray absorption spectra (XAS) edge heights. At both the L3 and L2 edges (around 11.57 keV and 13.28 keV, respectively), the integrated XMCD spectra show temperature dependences. In contrast, there is no temperature dependence in XAS spectra (see inset of Fig. 3a). The total magnetic moment (mtotal), which is the sum of an effective spin magnetic moment \(\left( {{\boldsymbol{m}}_{\mathrm{s}}^{{\mathrm{eff}}} = {\boldsymbol{m}}_{\mathrm{s}} + {\boldsymbol{m}}_{\mathrm{D}}} \right)\) and mo, was estimated by a sum rule calculation35 (see the details about the sum rule calculation in Supplementary Note 4). In this study, we use moment values per single hole rather than those per atom because it is difficult to precisely determine the hole number (nh) of Pt in the Pt/Co/MgO. The moment values normalized by nh are directly obtained from the sum rule formula without considering nh35. The changes in the induced magnetic moments with temperature are small (~15%) and comparable to the error range of the analysis as shown in Fig. 3b. This suggests that there is a weak correlation between the temperature dependences of the proximity-induced magnetic moments of the Pt layer and the DMI in the Pt/Co/MgO system11.

Fig. 3

Temperature dependence of the proximity-induced moment in the Pt layer. a The XMCD and integrated XMCD spectra at the Pt L3 and L2 edges in terms of temperature. The XMCD spectra were fitted by the Lorentzian function, and the fitted curves were integrated to get the integrated XMCD spectra. The insets are the XAS spectra at the Pt L3 and L2 edges. b Temperature dependence of the \({\boldsymbol{m}}_{\mathrm{s}}^{{\mathrm{eff}}}\) and \({\boldsymbol{m}}_{{\mathrm{total}}}\)=\({\boldsymbol{m}}_{\mathrm{s}}^{{\mathrm{eff}}}\)+mo. The error bars are based on the standard deviation of the Lorentzian-fitting residual

Correlation between asymmetric orbital structure of the FM layer and the DMI

In this section, we study the correlation of DMI with various physical quantities associated with the magnetism of Co such as ms, \({\boldsymbol{m}}_{\mathrm{o}}^\parallel\) and \({\boldsymbol{m}}_{\mathrm{o}}^ \bot\) and the intra-atomic dipole moment mD. As explained in the previous section, ms,eff is the sum of ms and mD. The ms and mD values in terms of temperature were obtained from the relation \({\boldsymbol{m}}_{\mathrm{D}}\left( \theta \right) = {\boldsymbol{m}}_{\mathrm{D}}(0^\circ ) \cdot (1 - 3{\mathrm{cos}}^2\theta )\), thereby \({\boldsymbol{m}}_{{\mathrm{s}},{\mathrm{eff}}}\left( \theta \right) = {\boldsymbol{m}}_{\mathrm{s}} + {\boldsymbol{m}}_{\mathrm{D}}(0^\circ ) \cdot (1 - 3{\mathrm{cos}}^2\theta )\)28,34.

Details about the sum rule calculation are explained in Supplementary Note 4. The XAS and XMCD of the film are measured at 100, 200 and 300 K. Two incident angles (θ=0° and 70° with respect to the film normal) were used to separately estimate \({\boldsymbol{m}}_{\mathrm{o}}^\parallel\) and \({\boldsymbol{m}}_{\mathrm{o}}^ \bot\) values using the relation \({\boldsymbol{m}}_{\mathrm{o}}\left( \theta \right) = {\boldsymbol{m}}_{\mathrm{o}}^ \bot {\mathrm{cos}}^2\theta + {\boldsymbol{m}}_{\mathrm{o}}^\parallel {\mathrm{sin}}^2\theta\)36. Figure 4a, b are typical XAS and XMCD spectra at the Co L2,3 edges obtained at 0° and 300 K. The XMCD spectra show a clear temperature dependence of their intensity at both 0° and 70° as shown in Fig. 4c–f. In case of the ms, it increases by 20% (from 0.87 μB/μh to 1.11 μB/nh) as temperature decreases from 300 to 100 K (see Fig. 4g), which is consistent with the observed temperature dependence of Ms.

Fig. 4

Temperature dependence of the Co magnetic moments. a XAS spectra for positive (σ+) and negative (σ−) X-ray helicities. b XMCD and integrated XMCD spectra at 0° with 300 K. The temperature dependence of XMCD spectra at the Co L3 and L2 edges measured at (c, d) 0° and (e, f) 70°. g Plots of ms and mD vs temperature as a function of the X-ray incident angle. h\({\boldsymbol{m}}_{\mathrm{o}}^ \bot\), \({\boldsymbol{m}}_{\mathrm{o}}^\parallel\) and \({\boldsymbol{m}}_{\mathrm{o}}^ \bot /{\boldsymbol{m}}_{\mathrm{o}}^\parallel\) as functions of temperature. The error bars are based on the standard deviation of the integral distribution after subtracting the backgrounds of XAS spectra

In addition to ms, mD shows a strong temperature dependence. The mD reflects the anisotropy of fSD distorted by the SOC or the crystal-field effect.28,29. When the electron occupation becomes asymmetric, for instance, more electrons in the in-plane orbitals than the out-of-plane orbitals, the situation gives rise to both non-zero mD term and asymmetry in the charge distribution. As a result, asymmetry of the charge distribution naturally implies non-vanishing fSD at transition metals' interfaces. As a result, the value of mD determined from the sum rule analysis increases from 0.014 μB/nh to 0.094 μB/nh (Fig. 4g). In addition, the orbital anisotropy also increases as temperature decreases; \({\boldsymbol{m}}_o^ \bot\) increases from 0.058 μB/nh to 0.080μB/nh, whereas \({\boldsymbol{m}}_{\mathrm{o}}^\parallel\) slightly decreases from 0.044 μB/nh to 0.039 μB/nh (Fig. 4h). These results imply a correlation of DMI with mD and orbital anisotropy. Since the variation of the perpendicular component of the mo is also related to magnetocrystalline anisotropy (KU) based on the Bruno theory37, the correlation between DMI and KU is also reasonable (see the Supplementary Note 5). Figure 5a shows \({\boldsymbol{m}}_{\mathrm{o}}^ \bot /{\boldsymbol{m}}_{\mathrm{o}}^\parallel\)(=[\({\boldsymbol{m}}_{\mathrm{o}}^ \bot (T)/{\boldsymbol{m}}_{\mathrm{o}}^\parallel (T)\)]/[\({\boldsymbol{m}}_{\mathrm{o}}^ \bot (300{\mathrm{K}})/{\boldsymbol{m}}_{\mathrm{o}}^\parallel (300{\mathrm{K}})\)]), and \({\boldsymbol{m}}_{\mathrm{D}}(T)/{\boldsymbol{m}}_{\mathrm{D}}(300{\mathrm{K}})\) plots as a function of the normalized D=D=(T)/D(300K) where T inside of parentheses is the measurement temperature. The ratio \({\boldsymbol{m}}_{\mathrm{o}}^ \bot /{\boldsymbol{m}}_{\mathrm{o}}^\parallel\) increases by 53% as D increases, and mD also shows a clear correlation with the DMI; \({\boldsymbol{m}}_{\mathrm{D}}(100{\mathrm{K}})/{\boldsymbol{m}}_{\mathrm{D}}(300{\mathrm{K}})\)≈ 6.4. Because both orbital anisotropy and mD are closely related to the orbital occupation with ISB28,37,38, these results suggest that the temperature dependence of DMI is governed by the change in asymmetric electron occupation in orbitals in addition to spin magnetic moments as we discuss with the following theoretical studies. We also note that peak intensity of the XAS spectra reflects the hole number in 3d orbitals. We found small change in the absorption intensity when varying the temperature. Quantitatively, difference in the XAS intensity integral (IXAS) between 300 and 100 K is about 4% \(\{ [I_{{\mathrm{XAS}}}\left( {300{\mathrm{K}}} \right) - I_{{\mathrm{XAS}}}\left( {100{\mathrm{K}}} \right)]/ I_{{\mathrm{XAS}}}\left( {300{\mathrm{K}}} \right) = 0.04\}\). Hence, this indicates that there is a small change in the hole number within the temperature range, reflecting temperature dependence of the charge distribution in orbitals. Nonetheless, the value is so small that it does not affect the trend in the respective quantities normalized by nh from the sum rule calculation.

Fig. 5

Correlation of DMI with \({\boldsymbol{m}}_{\mathbf{o}}^ \bot /{\boldsymbol{m}}_{\mathbf{o}}^\parallel\) and mD, and theoretical calculations based on the tight binding model and DFT. a Normalized \({\boldsymbol{m}}_{\mathrm{o}}^ \bot /{\boldsymbol{m}}_{\mathbf{o}}^\parallel\) and mD vs D. Values of \({\boldsymbol{m}}_{\mathrm{o}}^ \bot /{\boldsymbol{m}}_{\mathbf{o}}^\parallel\) and mD for all temperatures are normalized by the values measured at 300 K. b Calculated \({\boldsymbol{m}}_{\mathrm{o}}^ \bot\), \({\boldsymbol{m}}_{\mathbf{o}}^\parallel\) and EDMI values based on the tight-binding model as a function of temperature. Inset shows \(f_{{\mathrm{SD}}}^ \bot\)/\(f_{{\mathrm{SD}}}^\parallel\) as a function of temperature. c Physical parameters such as total D (Dtot), \({\boldsymbol{m}}_{\mathrm{o}}^ \bot /{\boldsymbol{m}}_{\mathbf{o}}^\parallel\), mD and ms obtained by DFT calculation. Strength and sign of Dtot are calculated around their magnetic ground state using the combination of the relativistic effect spin–orbit coupling with the spin spirals. A positive sign of Dtot indicates a left-rotational sense or left chirality

Theoretical consideration about the microscopic origin of the DMI

In order to support the insights obtained from our experimental study, we carry out two types of theoretical calculations: a tight-binding model calculation and ab initio calculation based on density function theory (DFT), and the main result is shown in Fig. 5b, c. The details of these two complementary theoretical studies can be found in Supplementary Notes 6 and 7.

For the tight-binding model, we extend the trimer model suggested by Kashid et al.10, which contains two magnetic atoms coupled to a spin–orbit coupled non-magnetic ion. Compared to Kashid’s trimer model, we add one more orbital (dxz) on the non-magnetic site, enabling the computation of the orbital anisotropy (see details in Supplementary Note 6). In order to describe the temperature dependence of parameters in our trimer model, we assume the level broadening increases with the temperature. This broadening can be phenomenologically explained by magnetization fluctuations and the electron–phonon interaction given by atomic vibrations. As this model calculation is too simple, we do not aim to give a quantitative explanation of the experimental data but provide a qualitative understanding of the experimental result.

Figure 5b shows that the tight-binding model calculation reproduces qualitatively our experimental observations; both the DMI and \({\boldsymbol{m}}_{\mathrm{o}}^ \bot\) decrease with temperature, while the \({\boldsymbol{m}}_{\mathrm{o}}^\parallel\) is almost temperature independent. Note that the \({\boldsymbol{m}}_{\mathrm{o}}^ \bot\) and \({\boldsymbol{m}}_{\mathrm{o}}^\parallel\) cases are related to electron hopping with and without ISB, respectively (Supplementary Note 6). In addition, the ratio of the spin density distribution (\(f_{{\mathrm{SD}}}^ \bot\)/\(f_{{\mathrm{SD}}}^\parallel\)) between the in-plane (xy) and out-of-plane (yz, xz) orbital states decreases as temperature increases as shown in the inset in Fig. 5b (In our tight-binding model, the spin distribution is defined as SDn=fn,↑−fn,↓, where fn,↑(fn,↓) is the occupation of spin up (down) state and n represents the orbital quantum state \(d_{xy}^A,d_{xy}^B,d_{xy}^C,d_{yz}^C,d_{xz}^C\).). This result implies that the spin distribution variation of the orbital states with ISB is key to change of the orbital moment and the DMI, which is in line with the following DFT result.

To examine the correlation between the DMI and the orbital anisotropy in realistic structures, we also perform ab initio calculations for Pt(111)/X ultra-thin films, where X is a 3d transition metal (X=V, Cr, Mn, Fe, Co, Ni). Based on the experimental observation, the microscopic origin of the correlation between the DMI and the orbital anisotropy involves the impact of the temperature on 3d-orbital magnetization and their electron filling. In this respect, it is instructive to vary the 3d transition metals on the Pt substrate and examine the general chemical trend. Details about this study are also discussed in Supplementary Note 7. Figure 5c shows the summary of the DFT calculation. Here changing the overlayer results in modification of the relative alignment between 3d and 5d orbitals and thereby in a modification of the charge distribution: more hybridization results in less asphericity of the charge distribution between the in-plane and the out-of-plane, as reflected by the change in mD. Our calculation results provide a physical insight on how the asymmetric charge distribution induces DMI. As we mentioned in the second section, DMI requires three ingredients: (i) SOC, (ii) magnetic exchange, and (iii) symmetry breaking. While the first two aspects are easily quantifiable, the latter is more difficult to apprehend. It has been reported that the larger the asymmetric spin distribution with symmetry breaking and SOC, the larger intra-dipole spin moment29. Therefore, our experimental and theoretical studies propose a method to address the symmetry breaking term. Figure 5c explains this aspect explicitly. The Mn 3d orbitals and the Pt 5d orbitals provide large magnetic exchange and SOC, respectively. In addition, the dipole moment term becomes maximized with the hybridization between Mn 3d and Pt 5d orbitals. As a result, DMI is maximal. Furthermore, both the DMI and the orbital anisotropy follow the same trend in their signs: Dtot and \({\boldsymbol{m}}_{\mathrm{o}}^ \bot /{\boldsymbol{m}}_{\mathrm{o}}^\parallel\) at the Pt/V interface have negative sign, while those of other interfaces have positive sign. We also find the correlation between ms and DMI as confirmed experimentally in this study and in ref. 32. In contrast, our calculations do not reproduce the experimental correlation reported above between DMI and the induced moment. This may be because of the difference between the actual film stack (~3 monolayers of Co) and the simulated system (1 monolayer of Co). In fact, the total DMI is dominated by the contribution of the first Co monolayer at the Co/Pt interface, while the Heisenberg exchange, which is a source of the induced moment, is affected by the total number of Co overlayers. This is why our calculation, limited to one Co overlayer, well reproduces the correlation of DMI with spin and orbital moments, but not with the induced moment.

Now, the remaining question is how the orbital asphericity can give rise to DMI. At the interface, in-plane orbitals of Co and out-of-plane orbitals of Pt are mixed through inter-atomic hybridization (in our tight-binding model: dxy of Co and dyz of Pt). In the presence of ISB, this hybridization results in a non-vanishing orbital angular momentum: the hopping trajectory of an electron produces an effective orbital angular momentum as illustrated in Fig. 6. Moreover, the large SOC on Pt mixes in-plane and out-of-plane orbitals, which is accompanied by spin canting (in our tight-binding model: dxy and dyz orbitals on Pt site). In other words, the SOC on Pt converts the interfacial orbital angular momentum of hopping electrons into spin canting, resulting in a non-collinear magnetic texture. In spite of its simplicity, this picture intuitively explains how the orbital asphericity generates the spin canting. A similar idea has been reported in ref. 21: the interfacial electric dipole induced by the asymmetric charge distribution results in chiral orbital angular momentum, and thereby spin canting via SOC. Although the latter approach is based on p-orbitals, it agrees with the intuitive picture based on d-orbitals.

Fig. 6

Schematic illustration of the spin canting. Explicit illustration of the onset of the orbital angular momentum and how SOC induces spin canting from it. The red and grey spheres represent the top ferromagnetic and bottom normal metal atoms in each layer, while the yellow and blue spheres represent the dxy (top) and dyz (bottom) orbitals, respectively. Dotted red lines and black curved arrow indicate the electron hopping path and effective electron orbiting motion due to the hopping, respectively

Discussion

Our experimental study on the temperature dependence of the DMI suggests that the interfacial DMI in FM/HM bilayers originates from Heisenberg exchange and the asymmetric charge distribution caused by the ISB, as evidenced by the orbital anisotropy and magnetic dipole moment (\({\boldsymbol{m}}_{\mathrm{o}}^ \bot /{\boldsymbol{m}}_{\mathrm{o}}^\parallel ,\)mD). Our DFT simulation and tight-binding calculation provide a clear evidence of the close link between the DMI and orbital physics. Based on the theoretical discussions, the temperature-dependent mD and mo indicate that increase in the temperature promotes the phonon-induced electron hopping between the out-of-plane and in-plane orbitals, thereby resulting in a reduced asphericity (or reduced asymmetry) of fSD over the d-orbitals and thus quenching mD.

In our experiment, however, the correlation between \({\boldsymbol{m}}_{\mathrm{o}}^ \bot /{\boldsymbol{m}}_{\mathrm{o}}^\parallel\) and the DMI is only semi-quantitative, i.e. the temperature-dependent change in the DMI does not perfectly scale with the temperature-dependent change in \({\boldsymbol{m}}_{\mathrm{o}}^ \bot /{\boldsymbol{m}}_{\mathrm{o}}^\parallel\). This semi-quantitative correlation between the DMI and mo demands a more detailed discussion. In systems with ISB, \({\boldsymbol{m}}_{\mathrm{o}}^ \bot\) can be decomposed into ISB-independent part and ISB-dependent part. Given that ISB is an essential ingredient for the DMI, there should be a direct correlation between the DMI and ISB-dependent mo, as evidenced by our tight-binding model calculation. However, \({\boldsymbol{m}}_{\mathrm{o}}^ \bot\) also has an ISB-independent part, which precludes a direct and quantitative correlation between the DMI and \({\boldsymbol{m}}_{\mathrm{o}}^ \bot\). This statement can be rephrased technically as the DMI involves only off-diagonal elements of the SOC operator10, while mo involves all of them. We note, however, that even with this uncertainty, the experimentally observed correlation between the DMI and mo anisotropy for the Pt/Co/MgO structure is rather clear, implying that the ISB-dependent mo would dominate over the ISB-independent one in this structure. Therefore, our findings based on both experimental and theoretical studies provide a link between orbital physics and spin–orbit-related phenomena such as the DMI, which are essential for spin–orbitronic devices.

Method

Film preparation and device fabrication

Si/Ta (1.5)/Pt (2)/Co (0.5)/MgO (2)/HfO (5) (in nm) film with perpendicular magnetic anisotropy was deposited on an undoped Si substrate by direct current magnetron sputtering and the atomic layer deposition technique. A 5-μm-wide Hall cross-structure were fabricated using the photo lithography and the Ar ion milling. For the XMCD measurement, the same stack film was prepared.

Nucleation field measurement

Angular-dependent coercivity of the Co/Pt Hall device was measured to estimate the Hn of the magnetic droplet at 300, 200, 150 and 100 K. The angle between magnetic field and the sample normal was varied from 0° to 89° rotating the electromagnet. At each angle, magnetic field was swept with in ±0.5 T to observe the coercivity. Details about this measurement are also discussed in ref. 24.

XMCD measurement

Soft X-ray: Soft X-ray absorption spectra were measured using the total electron yield method with 96% circularly polarized incident X-rays at the BL25SU at SPring-8. XMCDs at the Co L3 and L2 edges (in a range between 770 and 840 keV) were recorded in the helicity-switching mode with an applied magnetic field of 1.9 T. Homogeneity of the magnetic field was better than 99% for ɸ10 mm at the sample position. The incident light direction was inclined by 10° with respect to the magnetic field direction. Temperature was varied from 300 to 100 K using a continuous liquid He flow-type cryostat.

Hard X-ray: XMCD experiments using hard X-rays were carried out at BL39XU of SPring-8. A circularly polarized X-ray beam with a high degree of circular polarization (>95%) was produced with a transmission-type diamond X-ray phase retarder of 1.4-mm thickness. XAS of the film were observed at a 0.6 T magnetic field applied parallel to the X-ray propagation direction of which an incident angle was 0° with respect to the surface normal. The X-ray fluorescence yield mode was used to record the spectra. The X-ray energy was scanned around the Pt L3 and L2 edges in a range between 11.5 and 3.5 keV, reversing the X-ray photon helicity at 0.5 Hz. In this manner, two helicity dependent spectra I+ and I− were recorded simultaneously. Here I+ and I− denote the intensities when the incident photon momentum and the magnetization vectors are parallel and antiparallel, respectively. The XMCD spectrum, ΔI = I+−I−, is given by the difference of the two spectra. Detailed experimental set-ups for the soft and hard X-ray MCD measurements are described elsewhere39,40.

DFT calculation

To understand the behaviour of DMI in 3d/Pt(111) ultrathin films and its correlation with the orbital moment anisotropy (OMA) and magnetic dipole moment (Tz), we have performed DFT calculations in the local density approximation41 to the exchange correlation functional, using the full potential linearized augmented plane wave method in film geometry42 as implemented in the Jülich DFT code FLEUR43. Both collinear and non-collinear magnetic states have been studied employing an asymmetric film consisting of six substrate layers of Pt covered by a pseudomorphic 3d monolayer on one side of the film at the distance optimized for the lowest collinear magnetic states. For the non-collinear calculations, we used p(1x1) unit cell applying the generalized Bloch theorem (http://www.flapw.de). The OMA is obtained at relaxed geometry as total energy difference for two different magnetization directions employing the force theorem (the principle axes point along hard and easy axis). We considered 512 and 1024 k-points in the two-dimensional Brillouin zone (2D-BZ) for the scalar relativistic computation and the calculation with SOC treated within first-order perturbation theory, respectively. Note that high computational accuracy is required since energy differences between different magnetic configurations are tiny (~meV) in the present case.

In order to investigate the DMI, first we self-consistently calculate the total energy of homogeneous magnetic spin spirals employing the generalized Bloch theorem within the scalar-relativistic approach44. We have considered the energy dispersion E(q) of planar spin spirals, which are the general solution of the Heisenberg Hamiltonian, i.e. states in which the magnetic moment of an atom site Ri is given by Mi=M[cos(q·Ri), sin(q·Ri), 0] where q is the wave vector propagation of the spin spiral. By imposing the Néel spin spirals along the high symmetry lines of 2D-BZ, we can scan all possible magnetic configurations that can be described by a single q-vector. So, varying the q-vector with small steps along the paths connecting the high symmetry points, we find the well-defined magnetic phases of the hexagonal lattice: FM state at \({\bar{\mathrm \Gamma }}\)-point (q=0), anti-FM state at the \({\bar{\mathrm M}}\)-point, and periodic 120° Néel state the \({\bar{\mathrm K}}\)-point. When the energy E(q) along the high symmetry lines of 2D-BZ is lower than any of the collinear magnetic phases studied previously, the system most likely adopts an incommensurate spin-spiral magnetic ground state structure.

In a second step, we evaluate the DMI contribution from the energy dispersion of spin spirals by applying the SOC treated within first-order perturbation theory combined with the spin spirals5,42. Phenomenologically, the antisymmetric exchange interaction DMI has the typical form \(E_{{\mathrm{DM}}} = \mathop {\sum }\limits_{i,j} {\boldsymbol{D}}_{i,j} \cdot ({\boldsymbol{S}}_i \times {\boldsymbol{S}}_j)\), where Di,j is the DM vector that determine the strength and sign of DMI and Si and Sj are magnetic spin moments located on neighbouring atomic sites i and j, respectively. Considering the Néel-type out-of-plane configuration, the Di,j vector should be oriented in plane and normal to the q-vector. Note that DMI term must vanish for both configurations Néel-type in-plane and Bloch-type spin spirals due to symmetry arguments44. According to our definition, the vector chirality reads \(\boldsymbol{C}= C\hat{\boldsymbol c} = {\boldsymbol{S}}_i \times {\boldsymbol{S}}_{i + 1}\), where the direction of the vector spin chirality \(\hat{\boldsymbol c}\) is considered as spin rotation axis. Thus the left-handed (right-handed) spin spiral correspond to C = +1 (C = −1)45.

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgements

We also thank H.-W. Lee for fruitful discussion about the relation between the orbital magnetism and the DMI. This work was partly supported by JSPS KAKENHI Grant Numbers 15H05702, 26870300, 26870304, 26103002 and 25220604; JSPS Postdoctoral Fellowship program (Grant Number 2604316, P16314); Collaborative Research Program of the Institute for Chemical Research, Kyoto University; R & D project for ICT Key Technology of MEXT from the Japan Society for the Promotion of Science (JSPS) and the Cooperative Research Project Program of the Research Institute of Electrical Communication, Tohoku University. This work has also been performed with the approval of the SPring-8 Program Advisory Committee (Proposal Nos. 2015A0117, 2015A0125). A.M. and A.B. acknowledge support from King Abdullah University of Science and Technology (KAUST) and fruitful discussiojns with S. Blügel and G. Bihlmayer. G.G., P.-H.J. and K.-J.L. also acknowledge support from the National Research Foundation of Korea (NRF-2015M3D1A1070465, 2017R1A2B2006119). K.-J.K. was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea Government (MSIP) (Nos. 2017R1C1B2009686, NRF-2016R1A5A1008184) and by the DGIST R&D Program of the Ministry of Science, ICT and Future Planning (17-BT-02). S.K. was supported by Priority Research Centers Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2009-0093818).

Search for Takahiro Moriyama in:

Search for Kab-Jin Kim in:

Search for Teruo Ono in:

Contributions

S.K., K.U., K.-J.K., T.M. and T.O. conceived and designed the study. S.K., K.U., T.K. and D.C. fabricated the device and performed the nucleation field measurement. S.K., D.-H.K., P.-H.J., T.M., K.-J.K. and T.O. contributed to determine the DMI-induced effective field using the droplet model. S.K., M.S., Y.K., T.N. and K.-J.K. designed and performed XMCD experiment and analysed the data. K.U. and K.Y characterizes the magnetic properties of films. A.B., A.M. and K.N. supported theoretical analyses with ab initio calculation. G.G. and K.-J.L. supported the tight-binding model for the temperature dependence of the DMI and orbital magnetism. All authors discussed the results and wrote the manuscript.

Competing interests

The authors declare no competing interests.

Corresponding authors

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.