## Quantifying exchange forces of a spin spiral on the atomic scale

Nadine  Hauptmann 1* ,  Soumyajyoti  Haldar 2 ,  Tzu-Chao  Hung 1 ,  Wouter  Jolie ,  Mara  Gutzeit , 1 2 Daniel Wegner 1 , Stefan Heinze 2 , Alexander A. Khajetoorians 1

- 1 Institute for Molecules and Materials, Radboud University, 6525 AJ Nijmegen, Netherlands
- 2 Institut fŸr Theoretische  Physik  und  Astrophysik,  Christian-Albrechts-UniversitŠt, Kiel, Germany

* Correspondence to: n.hauptmann@science.ru.nl

The large interest in chiral magnetic structures for realization of nanoscale magnetic storage or logic devices has necessitated methods which can quantify magnetic interactions at the atomic  scale.  To  overcome  the  limitations  of  the  typically  used  current-based  sensing  of atomic-scale exchange interactions, a force-based detection scheme is highly advantageous. Here, we quantify the atomic-scale exchange force field between a ferromagnetic tip and a cycloidal  spin  spiral  using  our  developed  combination  of  current  and  exchange  force detection. Compared to the surprisingly weak spin polarization, the exchange force field is more sensitive to atomic-scale variations in the magnetization. First-principles calculations reveal  that  the  measured  atomic-scale  variations  in  the  exchange  force  originate  from different contributions of direct and indirect (Zener) type exchange mechanisms, depending on the chemical tip termination. Our work opens the perspective of quantifying different exchange mechanisms of chiral magnetic structures with atomic-scale precision using 3D magnetic exchange force field measurements.

Chiral nanoscale spin structures such as magnetic skyrmions 1, 2 , as well as individual magnetic adatoms 3 , have raised a lot of interest in recent years due to their potential for magnetic storage or logic devices at the atomic scale 4, 5 . This has necessitated methods that can quantify and manipulate magnetic interactions between coupled magnetic atoms with atomic-scale precision. Spin-sensitive approaches based on scanning tunneling microscopy (STM) are sensitive to local variations in the exchange interaction at the atomic scale by detecting inelastic tunneling 6, 7, 8, 9 , spin resonance 10, 11 , or magnetization between a spin-sensitive tip and a surface 12, 13 . Most notable of these methods is spin-polarized scanning tunneling microscopy (SP-STM), in which the total detected current is sensitive to the relative magnetization between the surface and a magnetized tip 14, 15, 16 . SP-STM has  been  widely  used  to  characterize  non-collinear  magnetization  textures 17,  18 driven  by  the Dzyaloshinskii-Moriya  interaction  (DMI) 19,  20 that  occurs  in  structures  with  broken  inversion symmetry such as surfaces due to spin-orbit coupling. In metals the DMI is mediated by conduction electrons 21 .  The  Hamiltonian  of  the  DMI  between  two  spins 𝐒" and 𝐒 # is 𝐻 = -𝐃"# ∙ (𝐒 " ×𝐒#) where 𝑫"# is the DM vector the direction of which is given by symmetry 19  and the cross-product implies that canted spin structures are preferred. However, the tunneling current is also sensitive to the tunneling anisotropic magnetoresistance 22 , the non-collinear magnetoresistance 23 , as well as structural and electronic variations, which are all convoluted in the overall measured signal.

Magnetic exchange force microscopy (MExFM) is a method, based on non-contact atomic force microscopy  (nc-AFM),  complementary  to  spin-sensitive  tunneling  spectroscopy.  It  directly measures  the  local  exchange  force  between  a  magnetic  tip  and  the  surface 24,  25 .  Force-based detection  of  magnetic  interactions  presents  numerous  advantages  over  current-based  detection. Most prominently, MExFM can address magnetic insulators 24, 26 , and it can be combined with simultaneous spin-polarized current detection (SPEX) 27  to disentangle the geometrical structure

from  electronic  properties  of  magnetic  structures 28, 29 . Moreover,  magnetic  exchange  force spectroscopy (MExFS) can directly quantify exchange forces 27, 30 , which, combined with firstprinciples calculations, allows for determining the interplay of various exchange mechanisms 31  as well as the role of chemical functionalization of the tip 32, 33 . Nevertheless, this method has been scarcely applied, in contrast to the non-magnetic 3D chemical force field method that has been widely  used  to  quantify  chemical  forces  between  the  tip  and  various  surfaces 34,  35 as  well  as molecules 36 .  However, 3D-MExFS contributes to answer open questions in surface magnetism, for example, what the role of spin-orbit coupling on the total force is, in analogy to the modification of the spin-averaged electronic structure known as the tunneling anisotropic magnetoresistance. A further question is if it is possible to detect and quantify various exchange mechanisms, including Zener-type exchange and anisotropic exchange, of chiral magnetic structures with atomic-scale resolution.

Here, we utilize SPEX to investigate a single atomic layer of Mn on W(110), which exhibits a cycloidal spin-spiral ground state driven by the DMI 17 . We demonstrate that the spatially resolved exchange  force  is  much  more  sensitive  to  the  atomic-scale  magnetic  exchange  force  field, including the local atomic sites within the magnetic unit cell, when directly compared to the spin polarization.  Utilizing  simultaneous  MExFS  and I z -spectroscopy  (MExFIS),  we  map  out  the evolution  of  the  exchange  force  field  and  spin  polarization  between  various  non-collinear orientations  of  the  tip  magnetization  relative  to  the  sample.  Using  first-principles  calculations based on density functional theory, taking into account both the tip and the surface, we understand the experimental findings of the exchange force field, and we quantify the influence of chemically different tip terminations on the overall distance dependence and the various possible exchange pathways. The surprisingly weak spin polarization of the tunneling current is explained by the

dominant contribution of low spin-polarized p z -states near the Fermi energy, E F, due to the large exchange splitting of d -bands in Mn.

## Results

## Atomic-resolution imaging of the Mn spin spiral.

We sense the atomic-scale spatial variation of the magnetic exchange interaction of the cycloidal spin  spiral  utilizing  constant-height  (ch)  imaging.  Fig.  1(d)  sketches  the  cycloidal  spin  spiral magnetic ground state along the [ 11 . 0 ] direction of the Mn monolayer, which has a periodicity of 𝜆		~ 12 nm, as shown previously by SP-STM imaging 17 . The angle between the magnetic moments of adjacent Mn rows is about 173¡ resulting in a locally almost antiferromagnetic rectangular c(2 x 2) unit cell (Fig. 1(a)). Fig. 1 shows simultaneously measured maps of (b) the current I T ch , and (c) the frequency shift Δ f ch and at a constant height z 1 that is determined by opening the feedback loop at  height z 0,  at  which  constant-current  (cc) SP-STM  images  (Supplementary Fig. 1) were acquired ( V s = -10 mV; I T cc = -2 nA), and moving the tip by a tip displacement of Δ z = -0.29 nm. We note that negative Δ z reflects smaller tip-surface distances and Δ z = 0 nm corresponds to our reference height z 0). I T ch shows a stripe pattern along the [001] direction (Fig. 1(b)), illustrating spin contrast consistent with the spin spiral order reported in Ref. 17 . A long-range modulation of the stripe-pattern contrast is observed with a periodicity of Δ x m =  (5.70 ± 0.03) nm reflecting half the periodicity ( 𝜆/2 ) of the spin spiral 17 . We refer to the position with maximum contrast of the stripes as 𝜑 = 0° and 𝜑 = 180° (Fig. 1(d)), which correspond to positions with reverse alignments of the out-of-plane magnetic moments. We have carefully acquired and analyzed all data to make sure that the tip has a dominantly magnetic out-of-plane sensitivity (see Supplementary Note 1 and

Supplementary Fig. 1 for more details). We therefore can assign the darker and brighter stripes of I T ch to Mn rows exhibiting an antiparallel ( ap ) alignment at 𝜑 = 0° and a parallel ( p ) alignment at 𝜑 = 180° with respect to the out-of-plane projection of the tip magnetization (see Fig. 1(d)). The periodicity of Δ x  = 1 (0.43 ± 0.03) nm along the [ 11 . 0 ] corresponds to the lattice constant of the magnetic c(2 x 2) unit cell as sketched in Fig. 1(a) 14 . In agreement with previous works 14, 17 , the lattice constant along the [001] direction ( Δ x 2) is not resolved in I T ch . Imaging at constant height ( z 1) gives spin contrast in the I T ch image comparable to SP-STM images at z 0 (Supplementary Fig. 2 and Supplementary Note 2), but does reveal more information on the magnetism within the spin spiral.

In contrast, the simultaneously probed Δ f ch resolves the atomic lattice for both the [ 11 . 0 ] and [001] direction. The Δ f ch image in Fig. 1(c) reflects the force gradient averaged within z 1 ± z mod due to the oscillation of the tuning fork ( z mod = 50 pm). The total force that gives rise to the contrast in the image, is composed of chemical and magnetic exchange forces, and results in the observation of well-separated bright protrusions in the Δ f ch image, which we assign to the positions of the Mn atoms  (see  Supplementary  Note  3  and  Supplementary  Fig.  3  for  more  information  on  the assignment of atomic sites). We note that electrostatic forces due to local charges are not expected to  be  significant  because  of  the  short  screening  length  of  the metallic  substrate.  Our  extracted lattice constants ( Δ x 1 = (0.43 ± 0.03) nm and Δ x 2 = (0.30 ± 0.03) nm agree well with the predicted ones 14 . We refer to the top sites,  , within each c(2 x 2) unit cell, t i (Fig. 1(a)), within the spin spiral as t i ,  where i = 1 and   = 14 correspond to i 𝜑 = 0° and 𝜑 = 180° of the spin-spiral period (Fig. 1(d)), respectively. An analogous definition is used for the hollow sites hi (Fig. 1(a)). As the Mn layer  is  planar  due  to  pseudomorphic  growth  on  W(110) 14 ,  atomic  contrast  emerging  from chemical interaction is expected to result in the same variation of Δ f ch at all atomic sites ( t i and tÕi )

along the spin spiral.  However, we observe a modulation Δ f ex, i = Δ f i ch (z 1)  Δ fÕi ch (z 1)  of  the contrast between t i and tÕi along the [ 11 . 0 ] direction (cf. solid vs. dashed circles at position   = 13 i in Fig. 1(c)).  We obtain Δ f ex,1 ( z 1) = -(0.8 ± 0.2) Hz, Δ f ex,7 ( z 1) = (0.0 ± 0.2) Hz and Δ f ex,14 ( z 1) = (0.8 ± 0.2) Hz, reflecting the reversal of the out-of-plane projection of the surface magnetization within the spin spiral as seen with SP-STM (see also sketch in Fig. 1(d)). We therefore assign Δ f ex, i to variations of the magnetic exchange force along the spin spiral. We note that for different tips we have also observed an opposite correlation between I T ch and Δ f ch than that seen in Fig. 1(b) and Fig.  1(c).  Our  measurements  show  that  the  exchange  force  component  of  the  measured  force directly correlates with the out-of-plane projection of the surface magnetization. Therefore, this indicates that Δ f ch directly probes not only the local antiferromagnetic order, but is also sensitive to the variations of the magnetic exchange interaction of the non-collinear magnetic moments on the surface resulting from the Dzyaloshinskii-Moriya interaction driven spin spiral.

## Quantifying the atomic-scale magnetic exchange force field.

Next,  we  quantify  the  exchange  force  field  at  each  atomic  site  as  a  function  of  tip-sample separation. We address both the top and hollow sites to study the influence of the local lattice configuration on the detected magnetic exchange interaction between each site and the tip, which we were not able to resolve in SP-STM mode. First, we focus on out-of-plane magnetic orientation of the atoms in the spin spiral with p and ap alignment (Fig. 1(d)). We record Δ f i (from which we calculate Δ f ex, i ) together with I T, i with respect to Δ z using the same atomic-scale tip (see Fig. 2(a) and  the  insets  in  (b)  and  (c),  more  details  are  available  in  the  Supplementary  Note  4  and Supplementary Fig. 4). The error margin reflects an uncertainty of ±0.2 Hz for Δ f ex, i (Fig. 2(b,c), see Supplementary Note 2 for discussion of the influence of a surface tilt). Using the procedure

described  in  Ref. 37 ,  we  derive  the  magnetic  exchange  forces F ex, i from Δ f ex, i (Fig.  2(d,e)).  As expected, we observe a sign reversal of the distance dependence when comparing Δ f ex,1 ( Δ z ) and Δ f ex,14 ( Δ z ), as well as for F ex,1 ( Δ z ) and F ex,14 ( Δ z . ) Moreover, we observe a significant difference regarding the distance dependence of the exchange force field between top and hollow sites: for the hollow sites, the curves show a monotonous increase of the magnetic interaction (Fig. 2(c,e)). This is different from the curves acquired at the top sites (Fig. 2(b,d)) which show a change of the slope. The maximum absolute Δ f ex values are ~ 3 Hz resulting in magnetic exchange forces up to 40 pN, i.e. , smaller by two orders of magnitude compared to the total forces (Supplementary Note 4 and Supplementary Fig. 4). However, the observed exchange forces here are larger by about a factor of three than probed on the non-collinear magnetic nano-skyrmion lattice of a monolayer Fe on Ir(111) using also a bulk Fe tip 27 . The distance-dependent variations of the exchange force field suggest  that  the  measured  exchange  force  results  from  different  exchange  mechanisms  when comparing  data  for  the  top  and  hollow  sites,  as  a  result  of  the  different  nearest-neighbor configurations. Moreover, this evidence illustrates that F ex,i is highly sensitive to variations of the magnetic exchange force contributions on different atomic sites.

In order to have a measure for the distance-dependent spin-polarization, we consider the current asymmetry 𝐴" (𝛥𝑧) = ; &lt;,&gt; (?@)	-	;' &lt;,&gt; (?@) ; &lt;,&gt; (?@)C	;' &lt;,&gt; (?@) with I T, i ( Δ z ) acquired at t i or hi and IÕ T, i ( Δ z ) acquired at tÕi or hÕi .  The  error  margins reflect  an  uncertainty  of  ±3%  for I T, i .  The  different  signs  of A 1( Δ z )  and A 14( Δ z ) (Fig. 2(f, g)) reflect the reverse relative alignments of the magnetic moments of the tip and Mn atoms (see insets in Fig. 2(b,c)). We observe that A 1( Δ z ) and A 14( Δ z ) do not show a significant variation in amplitude regardless of the atomic site (hollow or top) or distance. This is different from previous distance-dependent measurements of the Fe/Ir(111) skyrmion lattice 27 , where A ( Δ z ) showed a monotonous increase. In addition, the absolute values of A observed here are one order

of magnitude smaller than for one monolayer Fe on Ir(111). Our findings here further show that detection  of F ex( Δ z ),  is  more  sensitive  to  atomic-scale  variations, i.e. top  and  hollow  sites, compared to A ( Δ z ), and provides more information about different exchange mechanisms.

So far, we have quantified the magnetic exchange force field at positions of the spin spiral with dominantly p and ap alignment (i.e., for   = 1 and 14). Next, we probe the exchange force field at i different  sites i along  the  spin  spiral  using  MExFIS,  in  order  to  quantify  the  variations  of  the exchange forces for the non-collinear orientations relative to the tip magnetization (Fig. 3a). We determine F ex, i ( Δ z ) and Ai ( Δ z ) (Figs. 3(b,c)) at the hollow sites ( hi and hÕi ) at which we probe larger magnitudes of the exchange forces than on the top sites for the p and ap alignment (Fig. 2(e)). The gradual color change from blue to olive reflects the position within the spin spiral as sketched in Fig. 3(a). We observe a direct correlation between the magnetization of the spin spiral and the measured magnitude of the magnetic exchange force as shown in Fig. 3(b) for ( 𝜆/4 ) of the spin spiral. Given our error margins, we observe a gradual decrease of the absolute values of F ex, i ( Δ z ) (as well as of Ai ( Δ z )) from   = 1 to   = 7. We assign the change of i i F ex, i ( Δ z ) to the variation of the magnetic exchange interaction due to the different directions of the magnetic moments between 𝜑 = 0° and 𝜑 = 90° (Fig.  3(a)).  We  note  that  we  also  acquired  data  for i =  7  to i =  14  (see Supplementary Note 5, Supplementary Fig. 5 and Supplementary Fig. 6). We carefully acquired and analyzed all data to make sure that the tip has a dominantly magnetic out-of-plane sensitivity and that we can rule out measurement artifacts, e.g. piezo creep (for more information on the various types of force curves acquired, see Supplementary Note 6 and Supplementary Fig. 7). Nonetheless, we cannot exclude a misalignment of the tip magnetization by a few degrees from the out-of-plane direction and small changes of absolute distance due to a slight tilt of the sample (see Supplementary Note 5 for more details). Furthermore, we determined a spontaneous change

of the tip magnetization during acquisition of the data for   = 7 to   = 14 (see Supplementary Fig. i i 6 for more information). Therefore, we have excluded this data from our discussion. Nevertheless, our results  show  that  we  can  probe  the  variation  of  the  magnetic  exchange  forces  and  current asymmetry along the spin spiral.

## First-principles calculations.

In order  to  understand  the  experimental  observations,  we  have  performed  first-principles calculations  for  Mn/W(110)  and  for  the  interaction  with  a  magnetic  tip  (see  Methods  for computational  details  and  Supplementary  Note  7  and  Supplementary  Fig.  8  for  the  studied geometries). Although an Fe bulk tip was used in the experiments, the tip was prepared by gentle dipping into the Mn layer with co-adsorbed Co adatoms. Therefore, we consider the influence of three different tip terminations (Fe, Co, Mn) of Fe-based tips. In our calculations we approached the tip on both the top and the hollow site of the Mn monolayer (cf. Fig. 1(a)) and considered a parallel ( p ) and an antiparallel ( ap ) alignment of the magnetic moments between the tip apex atom and the Mn surface atom (top site) or nearest-neighbor Mn surface atoms (hollow site). The effect of structural relaxations due to tip-sample interaction does not change our conclusions, as shown for  the  Fe  tip  (see  Supplementary  Note  8  and  Supplementary  Fig.  9). The  magnetic  exchange energy, 𝐸GH (𝑑) , can be calculated from the total energies of the ap - and p - configurations:

$$E _ { \text{ex} } ( d ) = \, E _ { \text{ap} } ( d ) - \, E _ { \text{p} } ( d ), \quad ( 1 )$$

where 𝐸JK (𝑑) and 𝐸K (𝑑) denotes the total energies of the ap - and p - configurations, respectively. Therefore, 𝐸GH (𝑑) &gt; 0 and 𝐸GH (𝑑) &lt; 0 indicate ferromagnetic  and  antiferromagnetic  coupling, respectively. Similarly, the magnetic exchange force 𝐹 GH (𝑑) is defined as:

$$F _ { \text{ex} } ( d ) = \ F _ { \text{ap} } ( d ) - \ F _ { \text{p} } ( d ). \quad \ ( 2 )$$

Here, 𝐹 JK (𝑑) and 𝐹 K (𝑑) represent the total forces on the tip along the [110] (i.e., the z -) direction of Mn/W(110) for ap - and p - configurations, respectively.

Fig. 4 shows the calculated magnetic exchange energies E ex for top (a,c) and hollow sites (e,g) and the  corresponding  exchange  forces F ex  (b,d,f,h)  as  a  function  of Δ d that  is  related  to  the experimentally determined Δ z by an offset as defined further below. Δ d is defined as d Ð d 0, where d 0 = 0.5 nm. Already at first glance, we observe that the exchange energy and forces qualitatively differ for the two sites on the Mn monolayer Ð in agreement with the experimental observation Ð as well as for different tip terminations. For Mn-terminated tips, we find a regime of ferromagnetic coupling on the top site at large distances, i.e. E ex( Δ d ) &gt; 0 for Δ d &gt; -0.16 nm (Fig. 4(a)). For smaller Δ d , the coupling  is  antiferromagnetic.  In  contrast,  for  a  Fe-  or  Co-terminated  tip,  only antiferromagnetic  coupling  is  observed  (Fig.  4(c)).  We  attribute  the  change  of  sign  for  Mnterminated tips to a transition  from  an  indirect  long-range  ferromagnetic  exchange  mechanism between the Mn d -states mediated by s electrons (Zener-type exchange), to a short-range direct antiferromagnetic exchange between d -orbitals of Mn apex and surface atom (as suggested in Ref. 31 ). The Zener-type exchange mechanism is much weaker for Fe and Co due to their smaller s-d coupling 38, 39, 40 , and no transition occurs, in agreement with a previous work 30 . The different types of  exchange  mechanisms  can  be  seen  in  spin-resolved  charge-density  difference  plots  (see Supplementary  Note  9  and  Supplementary  Fig.  10).  For  hollow  sites,  the  coupling  is  always ferromagnetic for the Mn-terminated tip (Fig. 4(e)). For Feand Co-terminated tips, antiferromagnetic coupling at large tip-sample separations changes to a ferromagnetic coupling at small tip-sample separations (Fig. 4(g)). The relatively small values of E ex and F ex for hollow sites and the change of sign for Fe- and Co-terminated tips can be attributed to competing exchange interactions between the tip apex atom and nearest neighbor and next-nearest neighbor surface Mn

atoms (cf. Fig. 1(a)). For Mn-terminated tips the Zener-type exchange dominates. Our calculations reveal  that  the  magnetic  exchange  force  is  highly  sensitive  to  competing  magnetic  exchange interactions between the tip and the Mn layer depending on the chemical tip termination as well as on the atomic site.

The calculated exchange force curves for a Mn-terminated tip (Fig. 4(f)) agree qualitatively with the experimental data shown in Fig. 2(d,e) while Fe- and Co-terminated tips (Fig. 4(h)) cannot explain the observations. From a comparison between the onset of magnetic exchange forces in the calculations and in the experimental curves we estimate that Δ d = 0.0 nm corresponds to a tip displacement of about Δ z = -0.16 nm. This estimate is also in reasonable agreement regarding the total measured and calculated forces (see Supplementary Fig. 3 and Supplementary Fig. 4). For the  discussion  of  the  calculated  and  experimental F ex,  we  focus  on  the  tip-sample  separations between Δ d =  -0.15  nm and Δ d = 0.0 nm corresponding to Δ z =  -0.31  nm and Δ z =  -0.16  nm, respectively. The experimental exchange force F ex ( Δ z )  on  top  sites  first  rises  up  to  an  absolute maximum value and changes its sign when the tip is brought closer to the surface (Fig. 2(d)). At the hollow sites, F ex( Δ z ) exhibits a monotonous increase of the absolute values (Fig. 2(e)). This agrees  with  the  change  in  the  calculations  from  positive  exchange  force  for Δ d &gt;  -0.1  nm  to negative exchange force for Δ d &lt; 0.1 nm on the top site for the Mn-terminated tip (Fig. 4(b)), while there is a monotonous rise of the exchange force on the hollow sites (Fig. 4(f)). We point out that we have observed many different types of distance-dependent curves experimentally (see Supplementary  Note  10  and  Supplementary  Fig. 11),  which  can  result  from  a  combination  of different tip terminations and tip geometries, as well as the presence of irreversible dissipative mechanisms such as surface-induced magnetization flipping of the tip or hydrogen diffusion (see

Supplementary  Note  6  and  Supplementary  Fig.  7).  For  our  discussion  here,  we  excluded  all experimental data showing any indications of irreversible dissipative processes.

To  understand  the  experimentally  observed  small  and  site-independent  spin  polarization,  we analyze the vacuum local density of states (LDOS) at the top and hollow sites. Within the TersoffHamann model of STM 41, 42 , the differential conductance is proportional to the LDOS a few GLYPH&lt;129&gt; above  the  surface,  and  the A ( Δ z )  curves  in  spin-polarized  measurements  (Fig.  2(f,g))  can  be compared to the spin polarization of the vacuum LDOS 43, 44 . As shown in the Supplementary Fig. 12, there is very little difference between the spin polarization on the hollow and the top site, in agreement with the experimental observation. In addition, we find quite a small value on the order of 20% around E F. As shown in the orbitally decomposed LDOS (see Supplementary Fig. 13 and Supplementary Note 11) the small spin polarization is due to p z states which dominate the vacuum LDOS. Typically, the spin-polarized LDOS on 3 d transition-metal films and surfaces is dominated by d states with large values of the spin polarization on the order of 50 to 80% 43, 45 . However, since Mn exhibits a very large exchange splitting and a large magnetic moment of 3.6 𝜇Q the d bands are far from the Fermi energy and p states dominate which are much less spin polarized. This is a rather different finding than for the recently studied Fe monolayer on Ir(111), where the magnetic moment is significantly smaller and the d -states strongly contribute to the spin polarization close to the Fermi energy 18 .

## Discussion.

In conclusion, we quantify the atomic-scale variation of the magnetic exchange force field between a  ferromagnetic  tip  and  a  chiral  magnetic  surface.  Using  constant-height  SPEX  imaging,  we spatially resolve the magnetic exchange interaction of the tip and the cycloidal spin spiral in one

monolayer Mn on the W(110) surface. While the constant-height current image shows the expected magnetic  contrast,  but  only  resolves  the  atomic  lattice  in  one  direction,  the  constant-height frequency shift image resolves the magnetism within the entire antiferromagnetic c(2 x 2) unit cell. Using MExFIS, we quantify the exchange force field along the spin spiral with atomic resolution and  observe  a  strong  influence  of  the  exchange  forces  which  depend  on  the  local  spin configurations on hollow and top sites, indicating different exchange regimes. In contrast, the sitedependent variations are absent in the current, reflecting that magnetic exchange forces are more sensitive  to  changes  of  the  atomic  environment.  Utilizing  first-principles  calculations,  we understand  those  site-dependent  differences  of  the  exchange  forces  and  quantify  the  chemical contribution of the tip apex to the overall distance dependence of the force field and the various possible  exchange  pathways.  The  experimental  findings  qualitatively  agree  with  the  exchange interactions probed by a Mn-terminated tip that reveals an interplay between antiferromagnetic and Zener-type ferromagnetic exchange mechanisms. The experimentally observed, surprisingly weak spin polarization is explained by the dominant contribution of p -states near E F due to the large exchange splitting in Mn.

Our work provides the first combined experimental and theoretical study to quantify the different exchange force regimes of a chiral magnetic structure on the atomic scale. As a future perspective, SPEX can be used to simultaneously quantify the influence of spin-orbit coupling on both the current and the magnetic interactions in chiral magnetic structures, as well as characterize them with antiferromagnetic probes or functionalized tips. Furthermore, 3D-MExFIS opens a promising route to quantify the various direct and indirect exchange forces, with atomic-scale resolution, in addition to extracting dissipation mechanisms, providing direct access to interaction energies in comparison to state-of-the-art methods based solely on tunneling current 7, 46 .

## Methods

## Experimental procedures.

We performed SPEX imaging using a bulk Fe tip with a magnetization out of plane to the surface (see  Supplementary  Note  1  for  more  details)  in  a  modified  commercial  ultra-high  vacuum STM/AFM 47   operating at a base temperature of T =  5.8  K.  Force detection was done in noncontact  frequency-modulation mode, utilizing a tuning fork-based qPlus sensor 48  with the free prong oscillating at its resonance frequency ( f 0 = 30.8 kHz) and with a quality factor larger than 25000. The bias voltage V s was applied to the sample, and we performed both constant-current (cc) and constant-height (ch) imaging. The W(110) surface was cleaned by annealing ( T ~ 1800 K) in an oxygen atmosphere ( p ~ 3 x 10 -7 mbar) followed by a final flash to T ~ 2800 K. The 0.8 monolayer of Mn were deposited from Mn granules (99.99% purity, MaTeck) by using an e-beam evaporator. The W(110) was kept at room temperature during deposition, and was annealed ( T ~ 800 K) after the Mn deposition. To characterize our magnetic tips, we utilize the spin-direction dependent appearance of co-adsorbed Co adatoms as a reference 49, 50 , which were deposited onto the cold sample ( T ~ 6  K) from a Co rod source (99.99% purity, MaTeck) by using an e-beam evaporator, with 1% of coverage (see Supplementary Note 1 for more details). Tips with spin contrast were obtained by gently dipping a bulk Fe tip into the Mn/W(110) monolayer. The tip was brought closer to the Mn layer by values between -600 pm to -900 pm with respect to the tunneling position z0, which is defined by the constant-current feedback setpoint V s = -10 mV and I T cc = -2 nA. We used the Δ f values at this feedback setpoint to characterize the bluntness of our tips, which were typically between -10 Hz and -17 Hz.

## Density functional theory.

We applied first-principles calculations using a plane wave-based density functional theory (DFT) code  VASP 51,  52 within  the  projector  augmented  wave  method  (PAW) 53,  54 .  We  used  the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange correlation 55, 56 . The combined system of tip and sample was calculated in a supercell geometry (see Supplementary Fig. 8). The Mn/W(110) surface was modelled using by a symmetric slab with 5 W layers and 1 layer of Mn on each side. The local magnetic order of the Mn/W(110) surface can be approximated as collinear antiferromagnetic (AFM) due to the long periodicity of the spin spiral ground state. We have used a c(6 × 6) AFM surface unit cell with the GGA lattice constant of W (3.17 GLYPH&lt;129&gt;), which is in good agreement with the experimental value of 3.165 GLYPH&lt;129&gt;. The tip was modelled by a 14-atom pyramid in a bcc(001) orientation (see Supplementary Fig. 8). We have used ferromagnetic Fe to model the tip. For the alloyed tips, we have used Co and Mn atoms as the apex atom. For a Mn-terminated tip apex, antiferromagnetic coupling occurs between the Mn and Fe atoms. The tip and the surface were initially relaxed independently before considering the coupled system, i.e. , the tip-sample interaction. For the isolated tips, i.e., without any tip-sample interaction, the apex atom and the four atoms of the adjacent layer have been relaxed. The in-plane interatomic distances between the base atoms are kept constant at the GGA lattice constant of Fe (2.85 GLYPH&lt;129&gt;). Note that for Fe tips we have also considered the effect of structural relaxations due to tip-sample interaction (see Supplementary Note 8). We have used a 450 eV energy as plane wave basis set cutoff. A (5 × 5 × 1) Γ . centered k -grid mesh was used for the Brillouin zone integration. For every tip-sample separation d (defined as the distance between the centers of the tip apex atom and  the  surface  Mn  atom  underneath  without  structural  relaxations  due  to  the  tip-sample interaction), we carried out two sets of calculation for p and ap alignment of the magnetic moments between the tip apex atom and the Mn surface atom (top site) or nearest-neighbor Mn surface

atoms (hollow site). The magnetic exchange energy was calculated from the total energies in the p and ap alignment. The forces were calculated by using the Hellmann-Feynman theorem 57 .

## Acknowledgements

We  acknowledge  financial  support  from  the  NWO  VIDI  project:  ÔManipulating  the  interplay between superconductivity and chiral magnetism at the single atom levelÕ with project number 680-47-534 and the NWO Physics/f grant (project: 680.91.003 SPEX). This project has received funding from the European Research Council (ERC) under the European UnionÕs Horizon 2020 research and innovation programme (SPINAPSE, grant agreement No 818399). S.Ha. and S.He. acknowledge financial  support from  the  Deutsche  Forschungsgemeinschaft  via  SFB677.  S.Ha. M.G. and S.He. thank the Norddeutscher Verbund fŸr Hoch-und Hšchleistungsrechnen (HLRN) for  providing  computational  resources.  W.J.  acknowledges  support  from  the  Alexander  von Humboldt Foundation via the Feodor Lynen Research Fellowship.

## Author contributions

N.H.,  S.He.  and  A.A.K.  designed  the  experiment.  T.-C.H.,  W.J.  and  N.H.  carried  out  the measurements and did the analysis of the experimental data. M. G. and S. Ha. performed the DFT calculations  for  the  clean  Mn/W(110)  surface.  S.  Ha.  performed  the  DFT  calculations  for  the interaction of the tip with the Mn/W(110) surface. S. Ha. and S. He. analyzed the results of the DFT calculations. N.H., T.-C.H., S.Ha., S.He., D.W. and A.A.K. wrote the manuscript. All authors contributed  to  the  discussion  and  interpretation  of  the  results  as  well  as  the  discussion  of  the manuscript.

## Additional information

Supplementary Information accompanies this paper at É

## Data availability

The authors declare that all relevant data are included in the paper and its Supplementary

Information file. Additional data are available from the corresponding author upon reasonable request.

## Competing interests

The authors declare no competing interests.

## References

- 1. Goncalves M.A.P., Escorihuela-Sayalero C., Garca-Fernandez P., Junquera J., Iniguez J. Theoretical guidelines to create and tune electric skyrmion bubbles. Sci. Adv. 5 , (2019), doi: 10.1126/sciadv.aau7023.
- 2. Bogdanov A.N., Yablonsky D.A. Thermodynamically Stable Vortexes in Magnetically Ordered Crystals - Mixed State of Magnetics. Zh. Eksp. Teor. Fiz 95 , 178-182, (1989), doi:
- 3. Donati F. , et al. Magnetic remanence in single atoms. Science 352 , 318-321, (2016), doi: 10.1126/science.aad9898.
- 4. Fert A., Cros V., Sampaio J. Skyrmions on the track. Nat. Nanotechnol. 8 , 152-156, (2013), doi: 10.1038/nnano.2013.29.
- 5. Khajetoorians A.A., Wiebe J., Chilian B., Wiesendanger R. Realizing All-Spin-Based Logic Operations Atom by Atom. Science 332 , 1062-1064, (2011), doi: 10.1126/science.1201725.
- 6. Heinrich A.J., Gupta J.A., Lutz C.P., Eigler D.M. Single-atom spin-flip spectroscopy. Science 306 , 466-469, (2004), doi: 10.1126/science.1101077.
- 7. Czap G. , et al. Probing and imaging spin interactions with a magnetic single-molecule sensor. Science 364 , 670-673, (2019), doi: 10.1126/science.aaw7505.
- 8. Hirjibehedin C.F., Lutz C.P., Heinrich A.J. Spin coupling in engineered atomic structures. Science 312 , 1021-1024, (2006), doi: 10.1126/science.1125398.
- 9. Verlhac B. , et al. Atomic-scale spin sensing with a single molecule at the apex of a scanning tunneling microscope. Science 366 , 623-627, (2019), doi: 10.1126/science.aax8222.

|   10. | Yang K. , et al. Tuning the Exchange Bias on a Single Atom from 1 mT to 10 T. Phys. Rev. Lett. 122 , 227203 (2019), doi: 10.1103/PhysRevLett.122.227203.                                                                        |
|-------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|    11 | Baumann S., Paul W., Choi T., Lutz C.P., Ardavan A., Heinrich A.J. Electron paramagnetic resonance of individual atoms on a surface. Science 350 , 417-420, (2015), doi: 10.1126/science.aac8703.                               |
|    12 | Meier F., Zhou L.H., Wiebe J., Wiesendanger R. Revealing magnetic interactions from single-atom magnetization curves. Science 320 , 82-86, (2008), doi: 10.1126/science.1154415.                                                |
|    13 | Zhou L.H., Meier F., Wiebe J., Wiesendanger R. Inversion of spin polarization above individual magnetic adatoms. Phys. Rev. B 82 , 012409, (2010), doi: 10.1103/PhysRevB.82.012409.                                             |
|    14 | Heinze S. , et al. Real-Space Imaging of Two-Dimensional Antiferromagnetism on the Atomic Scale. Science 288 , 1805-1808, (2000), doi: 10.1126/science.288.5472.1805.                                                           |
|    15 | Wortmann D., Heinze S., Kurz P., Bihlmayer G., BlŸgel S. Resolving complex atomic-scale spin structures by spin- polarized scanning tunneling microscopy. Phys. Rev. Lett. 86 , 4132, (2001), doi: 10.1103/PhysRevLett.86.4132. |
|    16 | Bode M. Spin-polarized scanning tunnelling microscopy. Rep. Prog. Phys. 44 , 523Ð582, (2003), doi: 10.1088/0034- 4885/66/4/203.                                                                                                 |
|    17 | Bode M. , et al. Chiral Magnetic Order at Surfaces driven by Inversion Asymmetry. Nature 447 , 190-193, (2007), doi: 10.1038/nature05802.                                                                                       |
|    18 | Heinze S. , et al. Spontaneous atomic-scale magnetic skyrmion lattice in two dimensions. Nat. Phys. 7 , 713-718, (2011), doi: 10.1038/nphys2045.                                                                                |
|    19 | Dzyaloshinskii I.E. Thermodynamic Theory of Weak Ferromagnetism in Antiferromagnetic Substances. Sov. Phys. Jetp-Ussr 5 , 1259-1272, (1957), doi:                                                                               |
|    20 | Moriya T. Anisotropic Superexchange Interaction and Weak Ferromagnetism. Phys. Rev. 120 , 91-98, (1960), doi: 10.1103/PhysRev.120.91.                                                                                           |
|    21 | Fert A., Reyren N., Cros V. Magnetic skyrmions: advances in physics and potential applications. Nat Rev Mater 2 , 17031, (2017), doi: 10.1038/natrevmats.2017.31.                                                               |
|    22 | Bode M. , et al. Magnetization-direction-dependent local electronic structure probed by scanning tunneling spectroscopy. Phys. Rev. Lett. 89 , 237205, (2002), doi: 10.1103/PhysRevLett.89.237205.                              |
|    23 | Hanneken C. , et al. Electrical detection of magnetic skyrmions by tunnelling non-collinear magnetoresistance. Nat. Nanotechnol. 10 , 1039-1042, (2015), doi: 10.1038/Nnano.2015.218.                                           |
|    24 | Kaiser U., Schwarz A., Wiesendanger R. Magnetic exchange force microscopy with atomic resolution. Nature 446 , 522-525, (2007), doi: 10.1038/nature05617.                                                                       |
|    25 | Schmidt R. , et al. Probing the Magnetic Exchange Forces of Iron on the Atomic Scale. Nano Lett. 9 , 200-204, (2009), doi: 10.1021/nl802770x.                                                                                   |

|   26. | Pielmeier F., Giessibl F.J. Spin Resolution and Evidence for Superexchange on NiO(001) Observed by Force Microscopy. Phys. Rev. Lett. 110 , 266101, (2013), doi: 10.1103/PhysRevLett.110.266101.                                                                                            |
|-------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|    27 | Hauptmann N., Gerritsen J., Wegner D., Khajetoorians A.A. Sensing Noncollinear Magnetism at the Atomic Scale Combining Magnetic Exchange and Spin-Polarized Imaging. Nano Lett. 17 , 5660-5665, (2017), doi: 10.1021/acs.nanolett.7b02538.                                                  |
|    28 | Hauptmann N. , et al. Revealing the correlation between real-space structure and chiral magnetic order at the atomic scale. Phys. Rev. B 97 , 100401, (2018), doi: 10.1103/PhysRevB.97.100401.                                                                                              |
|    29 | Grenz J., Kšhler A., Schwarz A., Wiesendanger R. Probing the Nanoskyrmion Lattice on Fe/Ir(111) with Magnetic Exchange Force Microscopy. Phys. Rev. Lett. 119 , 047205, (2017), doi: 10.1103/PhysRevLett.119.047205.                                                                        |
|    30 | Schmidt R., Lazo C., Kaiser U., Schwarz A., Heinze S., Wiesendanger R. Quantitative Measurement of the Magnetic Exchange Interaction across a Vacuum Gap. Phys. Rev. Lett. 106 , 257202, (2011), doi: 10.1103/PhysRevLett.106.257202.                                                       |
|    31 | Tao K., Stepanyuk V.S., Hergert W., Rungger I., Sanvito S., Bruno P. Switching a Single Spin on Metal Surfaces by a STM Tip: Ab Initio Studies. Phys. Rev. Lett. 103 , 057202 (2009), doi: 10.1103/PhysRevLett.103.057202.                                                                  |
|    32 | Lazo C., Caciuc V., Hšlscher H., Heinze S. Role of tip size, orientation, and structural relaxations in first-principles studies of magnetic exchange force microscopy and spin-polarized scanning tunneling microscopy. Phys. Rev. B 78 , 214416, (2008), doi: 10.1103/PhysRevB.78.214416. |
|    33 | Lazo C., Heinze S. First-principles study of magnetic exchange force microscopy with ferromagnetic and antiferromagnetic tips. Phys. Rev. B 84 , 144428, (2011), doi: 10.1103/PhysRevB.84.144428.                                                                                           |
|    34 | Hšlscher H., Langkat S.M., Schwarz A., Wiesendanger R. Measurement of three-dimensional force fields with atomic resolution using dynamic force spectroscopy. Appl. Phys. Lett. 81 , 4428-4430, (2002), doi: 10.1063/1.1525056.                                                             |
|    35 | Albers B.J. , et al. Three-dimensional imaging of short-range chemical forces with picometre resolution. Nat. Nanotechnol. 4 , 307-310, (2009), doi: 10.1038/Nnano.2009.57.                                                                                                                 |
|    36 | Braun D.A., Weiner D., Such B., Fuchs H., Schirmeisen A. Submolecular features of epitaxially grown PTCDA on Cu(111) analyzed by force field spectroscopy. Nanotechnology 20 , 264004, (2009), doi: 10.1088/0957- 4484/20/26/264004.                                                        |
|    37 | Sader J.E., Jarvis S.P. Accurate Formulas for Interaction Force and Energy in Frequency Modulation Force Spectroscopy. Appl. Phys. Lett. 84 , 1801-1803, (2004), doi: 10.1063/1.1667267.                                                                                                    |
|    38 | Zener C. Interaction between the d-Shells in the Transition Metals. II. Ferromagnetic Compounds of Manganese with Perovskite Structure. Phys. Rev. 82 , 403-405, (1951), doi: 10.1103/PhysRev.82.403.                                                                                       |
|    39 | Zener C. Interaction Between the d-Shells in the Transition Metals. Phys. Rev. 81 , 440-444, (1951), doi: 10.1103/PhysRev.81.440.                                                                                                                                                           |
|    40 | Zener C. Interaction between the d-Shells in the Transition Metals. III. Calculation of the Weiss Factors in Fe, Co, and Ni. Phys. Rev. 83 , 299-301, (1951), doi: 10.1103/PhysRev.83.299.                                                                                                  |

- 41. Tersoff J., Hamann D.R. Theory of the scanning tunneling microscope. Phys. Rev. B 31 , 805-813, (1985), doi: 10.1103/PhysRevB.31.805.
- 42. Tersoff J., Hamann D.R. Theory and Application for the Scanning Tunneling Microscope. Phys. Rev. Lett. 50 , 19982001, (1983), doi: 10.1103/PhysRevLett.50.1998.
- 43. NŽel N., Ferriani P., Ziegler M., Heinze S., Kršger J., Berndt R. Energy-resolved spin-polarized tunneling and exchange coupling of Co and Cr atoms on Fe islands on W(110). Phys. Rev. B 85 , 155406, (2012), doi: 10.1103/PhysRevB.85.155406.
- 44. Schšneberg J., Caffrey N.M., Ferriani P., Heinze S., Berndt R. Distance- and Spin-resolved spectroscopy of iridium atoms on an iron bilayer. Phys. Rev. B 94 , 115418, (2016), doi: ARTN 115418 10.1103/PhysRevB.94.115418.
- 45. Ferriani P., Lazo C., Heinze S. Origin of the spin polarization of magnetic scanning tunneling microscopy tips. Phys. Rev. B 82 , 054411, (2010), doi: 10.1103/PhysRevB.82.054411.
- 46. Willke P., Yang K., Bae Y., Heinrich A.J., Lutz C.P. Magnetic resonance imaging of single atoms on a surface. Nat. Phys. 5 , 1005-1010, (2019), doi: 10.1038/s41567-019-0573-x.
- 47. CreaTec Fischer &amp; Co GmbH, http://www.createc.de, 2019.
- 48. Giessibl F.J. High-speed force sensor for force microscopy and profilometry utilizing a quartz tuning fork. Appl. Phys. Lett. 73 , 3956-3958, (1998), doi: 10.1063/1.122948.
- 49. Serrate D. , et al. Imaging and manipulating the spin direction of individual atoms. Nat. Nanotechnol. 5 , 350, (2010), doi: 10.1038/nnano.2010.64.
- 50. Haldar S., Heinze S. Noncollinear spin density of an adatom on a magnetic surface. Phys. Rev. B 98 , 220401(R) (2018), doi: 10.1103/PhysRevB.98.220401.
- 51. Kresse G., FurthmŸller J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54 , 11169-11186, (1996), doi: 10.1103/PhysRevB.54.11169.
- 52. The Vienna Ab initio Simulation Package: atomic scale materials modelling from first principles, https://www.vasp.at, 2019.
- 53. Blšchl P.E. Projector Augmented-Wave Method. Phys. Rev. B 50 , 17953-17979, (1994), doi: 10.1103/PhysRevB.50.17953.
- 54. Kresse G., Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59 , 17581775, (1999), doi: 10.1103/PhysRevB.59.1758.
- 55. Perdew J.P., Burke K., Ernzerhof M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77 , 3865-3868, (1996), doi: 10.1103/PhysRevLett.77.3865.
- 56. Perdew J.P., Burke K., Ernzerhof M. Generalized gradient approximation made simple (vol 77, pg 3865, 1996). Phys. Rev. Lett. 78 , 1396-1396, (1997), doi: 10.1103/PhysRevLett.78.1396.

Figure 1: Atomic-resolution SPEX imaging of the Mn spin spiral. (a) Sketches of the antiferromagnetic c(2 x 2) unit cell which locally approximates the spin spiral. The extracted lattice constants from the experiment ( Δ x 1 and Δ x 2), the tip-Mn distance d , the different top ( t i and tÕi ) and hollow ( hi and hÕi ) sites are indicated for the tip positioned at tÕi (left) and hi (right). (b) Current I T ch and (c) frequency shift Δ f ch images of one half of the spin-spiral period (see (d), both images line-flattened, frequency shift image Gaussian-smoothed by two points) measured at constant height ( z 1) that is by 0.29 nm closer to the surface than the height z 0 at which the current feedback loop was opened ( V s = -10 mV and I T cc = -2 nA). Parameters: oscillation amplitude z mod = 50 pm, V s = -0.1 mV, tip magnetization normal to the surface. The arrow in (c) depicts the contrast variation due to the reversal of a single spin along the [ 11 . 0 ] direction. The scale bars in (c) and (d) correspond to 0.5 nm. (d) Side-view sketch of one half of the cycloidal spin spiral along the [ 11 . 0 ] direction in one monolayer Mn on W(110) together with the experimental measurement scheme. The magnetic moments of the spin spiral have a parallel ( p ) alignment with respect to the tip magnetization for tÕ 1 and  14, and an anti-parallel ( t ap ) alignment for  1 and t tÕ 14. Positions t7 and tÕ7 indicate the position with in-plane magnetization.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00959v2-figure-d256bfac-4b42-48fe-93b5-8ac7254f6dfa.jpg)

Figure  2: Quantifying  the atomic-scale magnetic exchange  force  field  at  top  and  hollow  sites. (a)  Sketch  of  the distance-dependent measurement scheme for the top and hollow sites along the spin spiral. The upper and lower sketch of the spiral reflects the color code for top and hollow sites, respectively. Data is acquired at top ( t 1, tÕ 1 and  14, t tÕ 14) and hollow ( h 1, hÕ 1  and h 14, hÕ 14)  sites  within  the  c(2  x  2)  unit  cell  (Fig.  1(a)).  (b,c)  Frequency  shift  difference Δ f ex  versus  tip displacement Δ z ,  reflecting  the  magnetic  exchange  force  contribution  ( Δ f ex( Δ z )  = Δ f ( Δ z )  Ð Δ fÕ ( Δ z )).  (d,  e)  Magnetic exchange force F ex( Δ z) derived from Δ f ex( Δ z ). (f,g) Current asymmetry A ( Δ z ). The grey error margins reflect an uncertainty of ±0.2 Hz and ±3% for Δ f ex and I T, respectively. All curves are smoothed by three points using a Savitzky-Golay filter. Prior to calculation of F ex, Δ f ex has been smoothed by five points. Parameters: oscillation amplitude z mod = 50 pm , V s = -0.05 mV, out-of-plane tip magnetization. The position Δ z = 0 nm corresponds to our reference height z 0, that is defined by a setpoint of the current feedback loop of V s = -10 mV; I T cc = -2 nA. Negative Δ z reflects smaller tip-surface distances.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00959v2-figure-f25b80fe-5f53-4a5f-971f-3bbbbdcac235.jpg)

Figure 3: (a) Quantifying the atomic-scale magnetic exchange force field along the spin spiral. Sketch of the cycloidal  spin  spiral  for  the  hollow  sites  and  illustration  of  the  distance-dependent  measurement  scheme.  (b) Magnetic exchange force F ex,i ( Δ z) derived from Δ f ex,i ( Δ z ), versus the tip displacement Δ z for different hollow sites hi / hÕi (with   = 1 to   = 7). (c) Current asymmetry i i 𝐴S (𝛥𝑧) . The grey error margins reflect an uncertainty of ±0.2 Hz and ±3% for Δ f i / Δ fÕi and I T,i / IÕ T,i, respectively. Parameters: oscillation amplitude z mod = 43 pm , V s = -0.1 mV. During the data acquisition, the tip magnetization may differ from the out-of-plane orientation which was characterized prior to and after the measurements. The position Δ z = 0 nm corresponds to our reference height z 0, that is defined by a setpoint of the current feedback loop of V s = -10 mV; I T cc = -2 nA. Negative Δ z reflects smaller tip-surface distances.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00959v2-figure-e6a53531-d551-406d-93de-9fada61020fd.jpg)

hollow site

Figure 4: First-principles calculations of the exchange interaction . Calculated magnetic exchange energies E ex on top (a,c) and hollow (e, g) sites for collinear alignments of the magnetic moments of tip and Mn atoms, for Fe-, Coterminated (c, g), and Mn-terminated (a, c) tips versus Δ d that is defined as d d Ð 0, where d 0 = 0.5 nm and d is the tipsample separation. Corresponding exchange forces F ex on top (b, d) and hollow sites (f, h) versus Δ d .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00959v2-figure-5007d4e7-ad78-4641-8e39-980524551207.jpg)
Nadine Hauptmann Soumyajyoti Haldar Tzu-Chao Hung Wouter Jolie Mara Gutzeit Daniel Wegner Stefan Heinze Alexander A. Khajetoorians
Quantifying Exchange Forces of a Non-Collinear Magnetic Structure on the Atomic Scale
The large interest in chiral magnetic structures for realization of nanoscale
magnetic storage or logic devices has necessitated methods which can quantify
magnetic interactions at the atomic scale. To overcome the limitations of the
typically used current-based sensing of atomic-scale exchange interactions, a
force-based detection scheme is highly advantageous. Here, we quantify the
atomic-scale exchange force field between a ferromagnetic tip and a cycloidal
spin spiral using our developed combination of current and exchange force
detection. Compared to the surprisingly weak spin polarization, the exchange
force field is more sensitive to atomic-scale variations in the magnetization.
First-principles calculations reveal that the measured atomic-scale variations
in the exchange force originate from different contributions of direct and
indirect (Zener) type exchange mechanisms, depending on the chemical tip
termination. Our work opens the perspective of quantifying different exchange
mechanisms of chiral magnetic structures with atomic-scale precision using 3D
magnetic exchange force field measurements.