Abstract
A cluster of a few magnetic atoms on the surface of a nonmagnetic substrate is one suitable realization of a bit for spin-based information technology. The prevalent approach to achieve magnetic stability is decoupling the cluster spin from substrate conduction electrons in order to suppress destabilizing spin-flips. However, this route entails less flexibility in tailoring the coupling between the bits needed for spin-processing. Here, we use a spin-resolved scanning tunneling microscope to write, read, and store spin information for hours in clusters of three atoms strongly coupled to a substrate featuring a cloud of non-collinearly polarized host atoms, a so-called non-collinear giant moment cluster. The giant moment cluster can be driven into a Kondo screened state by simply moving one of its atoms to a different site. Using the exceptional atomic tunability of the non-collinear substrate mediated Dzyaloshinskii–Moriya interaction, we propose a logical scheme for a four-state memory.
Similar content being viewed by others
Introduction
Information technology is currently changing the paradigm from the separate use of the charge and spin degrees of freedom of the electrons toward a combination of both properties in so-called spin-electronic elements1. Using the scanning tunneling microscope as a tool, arrays of a well-defined number of atoms and shape can be assembled on a nonmagnetic substrate2, and their spin state3 and excitations4 can be read out atom-by-atom. Recently, corresponding arrays consisting of magnetic atoms have been intensely studied as the basic constituents for future spin-based information storage5, 6 and processing7 schemes.
One of the main hurdles in the former respect is to achieve long spin-energy relaxation and -decoherence times8, 9 of the spin states, which is usually established by two ingredients. First, a large energy barrier separating the spin states, provided by the so-called magnetic anisotropy energy (MAE)10, is needed, in order to prevent the system from undergoing thermally induced fluctuations at sufficiently low temperatures. Second, the spin states need to be protected from quantum tunneling11, and from scattering by electrons (the so-called Kondo scattering12) or by phonons13 from the substrate, which can shortcut the MAE barrier. To this end, an appropriate tuning of the MAE14 and a decoupling of the spins of the array from the substrate electrons using rare earth atoms15, 16 in combination with thin decoupling layers13, 17, 18, superconducting19, or semiconducting substrates20, was usually pursued. However, with respect to spin processing, a tunable communication between the array-spins is essential. Following the decoupling layer approach, communication with stable spin states is so far restricted to dipolar and superexchange coupling which offers limited flexibility.
A much larger flexibility can be achieved by the contrary, i.e., strongly coupling the array-spin to substrate electrons leading to Ruderman–Kittel–Kasuya–Yosida (RKKY) interactions21, which can transfer spin information between two arrays7, and can be tuned in strength, sign and even non-collinearity22. In view of applications, it is, therefore, a formidable task to realize an array of few atoms with a long spin-energy relaxation time useable for information storage, which at the same time, is sufficiently coupled to substrate conduction electrons in order enable strong RKKY coupling to the array. Moreover, the use of heavy substrate materials would additionally feature strong spin-orbit coupling, leading to highly non-collinear Dzyaloshinskii–Moriya contributions to the RKKY interaction22, which can potentially add functionality to the system.
Long-range interactions are typically promoted in a substrate material that almost fulfills the Stoner criterion for ferromagnetism, such as Rh23, Pd24, or Pt10, 21, 22, 25,26,27,28. As a consequence, such materials display a large cloud of magnetically polarized host atoms, leading to a strongly enhanced total magnetic moment of the impurity-polarization cloud system, which was, therefore, referred to as a giant moment (GM) in the early days of research on magnetic alloys29. Theses heavy substrate materials also support non-collinear magnetization states of the GM systems and Dzyaloshinskii–Moriya interactions22, and are thus ideal candidates for the pursuit of a magnetically stable cluster, which can be accessed by strong non-collinear RKKY interactions. With respect to information storage and processing applications, the main open questions are, thus, whether the MAE of such a non-collinear GM system can be tuned such that Kondo screening30 and quantum tunneling processes13 are avoided and the system is driven to support stable spin states as in a classical magnet31.
In order to realize a system that encompasses all these requirements, we built giant moment clusters (GMCs), each consisting of only three iron (Fe) atoms that are assembled into a triangular array on the (111) surface of a platinum (Pt) single crystal. We show that we can indeed store information for hours in the non-collinear spin state of such a GMC, enabling a tuneable non-collinear interaction between atomic scale bits for flexible schemes of information processing.
Results
GMC fabrication
The assembly of the GMCs is shown in Fig. 1 (see Methods). Each Fe atom can sit on one of the two hollow adsorption sites, face centered cubic (fcc) or hexagonally closed packed (hcp), on the lattice of the Pt surface atoms32. The Fe atoms are free from adsorbed hydrogen33 as indicated by their small apparent height (see Supplementary Fig. 1). By lateral atom manipulation three Fe atoms are placed artificially on next nearest neighboring adsorption sites of the same type in the most compact triangular geometry (Fig. 1a). In this way, four different species of GMCs can be built (Fig. 1b, see Supplementary Note 1 and Supplementary Fig. 1 for the unambiguous determination of the stacking). For each of the two adsorption sites, fcc and hcp, the cluster can be either positioned such that a Pt atom is underneath the center of the trimer, referred to as top, or such that a hollow site is underneath the center, referred to as hollow. The two hollow site trimers have a larger apparent height (≈250 pm) in comparison to the top site trimers (≈200 pm). The electronic properties of all four trimer types are characterized by (inelastic) scanning tunneling spectroscopy ((I)STS, see Methods) as shown in Fig. 1c. Their spectra show distinct excitations (arrows), which are markedly different from the excitations of isolated Fe atoms outside of trimers32, indicating the emergence of a coherent electron state after the formation of the trimer.
Density functional theory calculations
Before we proceed with the detailed experimental investigation of these excitations, we describe spin-polarized density functional theory (DFT) calculations within the Korringa-Kohn-Rostoker Green function (KKR-GF) method (Fig. 2, see Methods, Supplementary Note 2 and Supplementary Fig. 2). Such calculations were performed in order to study the non-collinear GMC character of the trimers, which is expected from the large atomic mass and the Stoner enhancement of the substrate22, 28. We find, that, for all four cluster types, the magnetic moments of the three iron atoms are coupled approximately ferromagnetic with an exchange interaction J of several ten to hundred millielectron-volts (see Eq. 4 in Supplementary Note 2 for the used Hamiltonian and Supplementary Table 2 for the parameters). Moreover, there is a large cloud of more than a hundred spin-polarized Pt atoms in the vicinity of each cluster (Fig. 2a), adding up to a total magnetic moment of the whole cluster/substrate system of m tot ~ 12 μ B for fcc hollow/top and hcp top, and of m tot ~ 13 μ B for hcp hollow (see Supplementary Table 1 for the distribution of spin and orbital moments). Closer inspection shows that, while there are considerable non-collinearities of the spins in the polarization cloud, the Fe atom spins have tilting angles of only θ ≤ 4°, since the Dzyaloshinskii–Moriya contribution D || to the exchange is an order of magnitude smaller than J. The only exception is the hcp top cluster, where \(\frac{{{D_{||}}}}{J} \sim 0.5\) leads to tilting angles of up to θ ~ 17° (see Supplementary Table 2 and Supplementary Figs. 3 and 4). The clusters’ magnetizations have strong MAEs K on the order of one milli-electronvolt per Fe atom (Fig. 2b, c). Interestingly, the polarization cloud has an RKKY-like contribution to the MAE26, as revealed by its behavior as function of the number of Pt atoms considered in the collinear calculation (Fig. 2b). For a sufficient number of Pt atoms taken into account, the calculations predict a preferred orientation of the magnetization perpendicular to the surface (“out-of-plane”) for both fcc clusters, and a preferred orientation of the magnetization in the surface plane (“easy-plane”) for the hcp hollow cluster. Interestingly, the MAE favors the same orientation as in the system where the constituent atom is individually adsorbed on the surface28. The case of hcp top cluster is peculiar in that, although the MAE favors an easy-plane orientation by approximately 1.5 meV per adatom in the collinear calculation (see Fig. 2b), the strong non-collinearity emerging from the Pt substrate reverses the trend and favors instead a nearly out-of-plane orientation by ~ 0.1 meV per adatom in the non-collinear calculation (Fig. 2c and Supplementary Fig. 5). The condition for the occurrence of this remarkable effect can be cast into the simple expression \(\left| {\frac{{9D_{||}^2}}{{4J}}} \right| >3K\) (see Eq. 21 in Supplementary Note 2), which shows that the effective anisotropy of the full trimer results from the competition between the non-collinear Dzyaloshinskii–Moriya contribution and the MAE. All in all, our DFT calculations reveal considerable MAEs of all four GMCs and a strong Dzyaloshinskii–Moriya–type coupling to the substrate conduction electrons, which is reflected by the non-collinear substrate polarization cloud.
Effective spin model
In order to get a rough idea of the spin degrees of freedom of these GMCs, we map the DFT results onto an effective spin Hamiltonian, which will be later used to link the experiments to the ab-initio results. Considering the application of a magnetic field B in the out-of-plane direction, the lowest-order approximate effective spin Hamiltonian reads \(\hat H = g{\mu _{\rm{B}}}{\hat S_z}B + {\cal D}\hat S_z^2\) 11. Here, \({\hat S_z}\) is the operator of the out-of-plane component of the total cluster effective spin, g is the Landé g-factor, μ B is the Bohr magneton, and \({\cal D}\) is the longitudinal magnetic anisotropy constant. Since the Fe atoms and their nearest neighbor Pt atoms are strongly exchange-coupled among each other (see Supplementary Table 2), the GMC is treated in this Hamiltonian as a single object with a well-defined macro spin. Neglecting the slight non-collinearity of the three ferromagnetically coupled Fe atoms plus the neighboring Pt atoms in the cluster, the DFT calculated total magnetic moments result in a mapping to a total spin with quantum number of \(S = - \frac{1}{2} + \sqrt {{{\left( {\frac{{{m_{{\rm{tot}}}}}}{{g{\mu _{\rm{B}}}}}} \right)}^2} + \frac{1}{4}} \approx 5 \) or \(\frac{{11}}{2}\) for the fcc hollow/top and hcp top GMCs, and of \(S \approx 7\) or \(\frac{{15}}{2}\) for the hcp hollow GMC. Here, we considered the g-factors extracted from ISTS as given in Table 1. Please note, that the effective spin model is a crude approximation for the case of magnetic impurities strongly interacting with a metallic substrate, since charge fluctuations can inhibit well-defined spin quantum numbers33. In the following, we choose the closest half-integer spin values (\(\frac{{11}}{2}\) for fcc hollow/top and hcp top, and \(\frac{{15}}{2}\) for hcp hollow) as we will see below that the hcp hollow cluster shows the Kondo effect, which strongly indicates a half integer spin ground state. However, the following discussion remains valid for the other choices of S by adjusting \({\cal D}\) accordingly. Taking into account the positive and negative values of \({\cal D}\) as estimated from the DFT-derived K-values for the easy-plane and out-of-plane clusters (see Table 1), the resulting qualitative energy level diagrams of the spin Hamiltonian for the two cases are shown in Fig. 2d, e, respectively, in zero magnetic field. We expect a degenerate \({S_z} = \pm \frac{1}{2}\) ground state for the hcp hollow cluster, and a degenerate \({S_z} = \pm \frac{{11}}{2}\) ground state for the fcc hollow/top and hcp top clusters. The former case is amenable to a quantum-mechanical superposition of the two \(\frac{1}{2}\) spin states of the effective spin with the spin states of the substrate conduction electrons involving Kondo-correlations33, 34, while the latter typically favors magnetic bi-stability, similar to a classical magnet5. We, therefore, expect that the magnetic properties of the hcp hollow cluster will be drastically different from that of the other clusters.
GMC spin excitations
First, we experimentally investigated the excitations (Fig. 1c) of all four clusters by (I)STS as a function of magnetic field B, and for the hcp hollow GMC also as a function of temperature T. Indeed, for the hcp hollow cluster (Fig. 3a, c), we find a Fano-resonance right at the Fermi energy E F (i.e., at V = 0 meV) that splits and broadens as a function of B and T, respectively. These are clear hallmarks for an emergent Kondo screening of the spin of this GMC33,34,35,36,37. The asymmetric lineshape of the Fano resonance is caused by the interference of a direct tunneling into the Kondo resonance, which arises as a consequence of the zero-energy spin flip scattering38, and other tunneling channels39. We can adequately fit the spectra to a sum of two Fano functions using the same form factor q and full width at half maximum Γ, but allowing for different energetic positions E i of the Kondo resonances (lines in Fig. 3a, c, see the used Fano functions and parameters in Supplementary Note 3, Supplementary Tables 3 and 4, and Supplementary Fig. 6). The extracted values are used to quantify the field-induced splitting and temperature-driven broadening, as shown in Fig. 3b, d. Indeed, the resonance shifts linearly with B (Fig. 3b) with a g-factor close to 2 (g hcp hollow ~ 1.7 ± 0.2), which we expect from the Kondo screening of the spin \(\frac{1}{2}\) ground state of the easy-plane anisotropic hcp hollow cluster (Fig. 2d)33, 34. The dependence of Γ on temperature (Fig. 3d) nicely fits to a power law and to numerical renormalization group calculations for a spin-1/2 impurity40 (see Supplementary Note 3) using a Kondo temperature of T K ~ 4.5 K. This value is one to two orders of magnitude smaller as compared to usual Kondo temperatures of single impurities on noble metal surfaces36, but has a similar size as that of atoms which are decoupled from the conduction electrons by thin insulators34 or that of Fe-hydrogen complexes on Pt(111)33. Our experimental findings thus prove Kondo correlations of the hcp hollow GMC with the substrate conduction electrons.
The magnetic-field evolution of the excitations of the fcc hollow/top and hcp top clusters behave drastically different (Fig. 1c, Supplementary Note 4, Supplementary Figs. 7 and 8). The spectra on these GMCs reveal a pair of features at energies symmetric to E F, which linearly shift away from E F with increasing magnetic field. This behavior is very reminiscent of systems with an out-of plane MAE28, which is excited from the ground state into the first excited state over the zero field energy splitting of \({{\it{\Delta}} _{01}} = {\cal D}\left( {2S - 1} \right)\) (Fig. 2e) that is linearly increasing by Zeeman energy. Note, however, that the shapes of the features deviate from simple steps, which probably indicates the limitation of the effective spin model (see Supplementary Note 4). By measuring the energetic positions of the excitations as a function of B (Supplementary Figs. 7 and 8 and Supplementary Table 5), we extract the corresponding longitudinal anisotropy constants \({\cal D}\) and g-factors close to 2, which are given in Table 1. The anisotropy constants for the two fcc and the hcp top clusters fit reasonably well with the constants estimated from the DFT calculations (see Table 1), substantiating the out-of-plane character of these GMCs.
GMC spin dynamics
Due to the out-of-plane MAE of the fcc top/hollow and hcp top GMCs, we expect a bi-stable switching between the \({S_z} = \pm \frac{{11}}{2}\) ground states, which is observable if there are not too frequent tunneling electrons with an energy larger than Δ01, which will drive the system across the barrier between these two states by sequential excitations (Fig. 2e)5. On the other hand, the emergent Kondo behavior of the hcp hollow GMC should lead to a quantum mechanical superposition of the \({S_z} = \pm \frac{1}{2}\) states of the effective spin with the spin states of substrate conduction electrons such that there is no observable switching of the magnetization. Indeed, using Cr coated tungsten tips, which are sensitive to the out-of-plane component of the GMC magnetization (Methods and Supplementary Fig. 9), we observe striped patterns in constant-current images taken on fcc top/hollow and hcp top clusters (shown exemplarily for the fcc top cluster in Fig. 4a), but not on the hcp hollow GMC. Such stripe patterns result from a random two-state telegraph noise5, 31, which is visible if the tip is positioned on-top of those clusters and the height is measured as a function of time (Fig. 4b). The telegraph signal can be unambiguously assigned to a statistical switching of the GMC between two spin states (0) and (1) of opposing out-of-plane magnetizations (Fig. 4d) by verifying that the asymmetry of the time-averaged occupational lifetimes \({\bar \tau _0}\) and \({\bar \tau _1}\) of the two states, defined by \({\cal A} = \frac{{{{\bar \tau }_1}}}{{{{\bar \tau }_0} + {{\bar \tau }_1}}}\), reverses by either reversing the out-of-plane oriented magnetic field B, or by reversing the out-of-plane spin-polarized tunneling current (bias polarity), which both drive the spin into one of the two states5, 31. We prove, that this is indeed the case for the fcc top/hollow and hcp top GMCs as shown in Fig. 4b, c. Similarly, full B-dependent loops of the asymmetry for two bias polarities recorded for all three clusters (Fig. 4e) show all the characteristics of a system randomly switching between two out-of-plane spin states and driven by the spin-pumping effect of the spin-polarized tunneling electrons5. Note, that there is a considerable shift of these asymmetry curves toward negative magnetic fields, which indicates an effective stray or exchange field of the tip3, 41 on the order of +0.1T. Moreover, the asymmetry curve recorded on the fcc hollow GMC is apparently inverted with respect to the vertical axis, which might be related with an inverted vacuum spin-polarization of this cluster as compared to the other two clusters, as sensed by the tip apex.
We consequently identify the two spin states (0) and (1) observed in the telegraph noise with the two ground states \({S_z} = + \frac{{11}}{2}\) and \({S_z} = - \frac{{11}}{2}\) of the fcc top/hollow and hcp top GMCs (Fig. 2e). In order to study the dynamics in these spin states in detail, we measure the telegraph noise as a function of bias voltage, tunneling current, and temperature (Fig. 4f–h). The temperature dependence of the average lifetime \(\bar \tau = {\left( {\frac{1}{{{{\bar \tau }_0}}} + \frac{1}{{{{\bar \tau }_1}}}} \right)^{ \!\! - 1}}\) (Fig. 4f) is slightly different for all three clusters and on the order of a few tens of seconds at T = 300 mK. It reveals an exponential decay with increasing temperature, which can be assigned to a quasi-classical Arrhenius-like behavior that has been observed for larger scale Fe islands31 and non-C 3v -symmetric clusters at larger temperatures5. The absence of any shoulder in the temperature dependence down to our lowest measurement temperature indicates, that transversal anisotropies, which could induce quantum tunneling processes, are negligible here5. Furthermore, as shown in Fig. 4g, \(\bar \tau \) increases over more than two decades if the bias voltage approaches the excitation thresholds \(\frac{{{{\it{\Delta}} _{01}}}}{e}\) determined from ISTS (Supplementary Table 5), and finally gets too long for a statistically relevant measurement. This behavior is indeed expected, considering that, without quantum tunneling through the barrier, the spin state of the GMC can only overcome the barrier by sequential excitation processes, which need either tunneling or substrate electrons with energies larger than Δ01. Indeed, we can consistenly reproduce all measurements by using a master equation approach (see Methods) in conjuction with the effective spin Hamiltonian given above, with \(S = \frac{{11}}{2}\), the g-factors extracted from ISTS (see Table 1) and fitting the longitudinal magnetic anisotropy parameter \({\cal D}\) (lines in Fig. 4e–g). The resulting magnetic anisotropies given in Table 1 agree with the parameters determined from ISTS remarkably well, giving a final support for the quasi-classical spin behavior of the fcc top/hollow and hcp top GMCs, in stark contrast to the Kondo character of the hcp hollow GMC.
GMCs as two-bit registers and four-state memories
As suggested by Fig. 4g, a spin-polarized tunneling current of sufficiently small bias below the excitation threshold \(\frac{{{{\it{\Delta}} _{01}}}}{e}\) could be used to “read” the quasi-classical spin state of the fcc top/hollow and hcp top GMCs with negligible disturbance. On the other hand, a current biased above this threshold could be used to switch the spin state and thereby “write” information. This is verified by building a two-bit register from two of the magnetically most stable fcc top GMCs located next to each other in a distance of only 2.5 nm (Fig. 5a–d). Indeed, we were able to write the state (0) or (1) into the desired GMC by application of a spin-polarized tunneling current driven by a voltage-ramp crossing \(\frac{{{{\it{\Delta}} _{01}}}}{e}\) that was switched off when the GMC switched from (0) to (1), or vice versa, at wish, without changing the state of the neighboring GMC. The consecutive read out of the four possible states (01), (11), (10), and (00) of the two-bit register written like this is shown in Fig. 5a–d. Furthermore, we demonstrated, that such prepared states are stable in time by reading the spin state of one of the CMCs prepared in state (0), while the other was in state (1), over time. Figure 5e shows that the state is stable for at least 10 h, after which we had to stop the measurement due to restrictions imposed by the experimental facility.
Since the GMCs are strongly coupled to a heavy metallic substrate that was experimentally shown to mediate Dzyaloshinskii–Moriya interactions tunable by the distance22, the realization of a two-bit register of spins, demonstrated here, offers further perspectives for spin-based information processing schemes. One possible scheme made from the two-bit register in Fig. 5a–d and an additional hcp-hollow cluster (or alternatively an Fe atom on hcp adsorption site which has easy-plane anisotropy33) in close to perpendicular geometry is shown in Fig. 5f. The hcp-hollow cluster should be assembled at a distance where the substrate mediated interaction is strong enough to quench the Kondo screening, permitted by its small Kondo temperature (k B T K ~ 0.4 meV). Since the Dzyaloshinskii–Moriya vector D has a strong contribution in the plane and perpendicular to the dotted line (see Fig. 5f)22, the spin of the hcp hollow cluster (or hcp atom) can be forced into the four different orientations shown in Fig. 5f–i depending on the four possible states of the two-bit register.
We finally compare the magnetic stability of the three-atomic fcc top GMC investigated here with the system of five Fe atoms constructed on the (111) surface of copper5. As shown in Fig. 4h, under writing conditions, a number of about \({\cal N} = \bar \tau \cdot I = {10^8}\) tunneling electrons (at 0.75 nA current and 5 mV bias, slightly depending on the current) is needed for a single switching event of the fcc top GMC. We investigated other fcc top clusters using different tips and found that \({\cal N}\) can be as large as 109 (0.75 nA, 10 mV), where the variation between clusters results from small variations in the MAE due to electronic inhomogeneity of the substrate5. This number is about five to ten times smaller than that needed to switch Fe5 on Cu(111) at comparable tunneling parameters. The increased sensitivity under writing conditions might be explained by the smaller number of atoms in the GMC investigated here, leading to a smaller spin and consequently to a smaller number of states separating the two spin states (0) and (1). Fully understanding why Fe3 on Pt(111) still is a magnetically stable GMC requires a detailed analysis of diverse effects, such as symmetry protection14 or quantum spin-fluctuations. For instance, the higher symmetry of the Fe3 on Pt(111) GMC leeds to negligible transversal anisotropies efficiently reducing quantum mechanical tunneling of the spin through the MAE barrier, which is not the case for the less symmetric Fe5 on Cu(111) system. Moreover, quantum spin-fluctuations can be suppressed if the hybridization with the substrate is less effective42, which appears to be the case for Pt(111) in comparison to Cu(111), according to our preliminary results. These effects are likely to be the reason for the fact, that the Fe3 on Pt(111) GMC can still serve as a stable magnetic bit under reading conditions, while it needs less electrons to switch its state under writing conditions.
Discussion
In summary, we experimentally demonstrated that a GMC of only three Fe atoms strongly coupled to a heavy-element substrate constitutes a stable quasi-classical magnet suitable for storing information for hours. Due to the strong coupling of the magnetic bit to a substrate that features non-collinear RKKY-interactions22, the spin-information can be processed with large variability, as we demonstrate by the scheme of a four-state memory. Moreover, the chosen substrate platinum is a prototypical strong spin-orbit coupling material widely used in the research on spintronics technology concepts. Therefore, the few-atom magnet we have realized here is an ideal model system to study the down-scaling of such concepts, as, for example, the application of spin-orbit torques for writing information43, 44, to the limit of single atoms.
Methods
Experimental procedures
All measurements have been performed under ultra-high vacuum conditions in a home-built low-temperature scanning tunneling microscope facility were the magnetic field B is applied perpendicular to the sample surface45. The Pt(111) single crystal was cleaned in situ by argon ion sputtering and annealing cycles and a final flash as described in ref. 28. While the sample cooled down to room temperature after the flash, a fraction of a monolayer Co was deposited from an e-beam heated rod resulting in the decoration of the Pt step edges and terraces with Co stripes and islands, respectively, of one atomic layer height3, 27 (see Supplementary Fig. 9). After the sample was cooled to T ~ 4 K, ~ 1% of a monolayer Fe was co-deposited from an e-beam heated rod. During deposition, the sample temperature did not exceed T ~ 10 K resulting in a statistical distribution of single Fe atoms on the Pt terraces (Fig. 1a and Supplementary Fig. 9).
Constant-current images were recorded at a tunnelling current I with a bias voltage V applied to the sample. The Fe clusters have been formed by lateral manipulation of Fe atoms using a current of I = 25–35 nA with voltages of V = 0.9–1.3 mV2, 7. For spin-resolved scanning tunneling microscopy (SPSTM), we coated flashed tungsten tips with about 50 monolayers of Cr32. Spin contrast with a sensitivity to the out-of-plane component of the sample magnetization was then achieved by picking up Fe atoms until a strong magnetic contrast was observed on the remanently out-of-plane magnetized Co stripes or islands3, 27 (Supplementary Fig. 9). For the investigation of the spin excitations by spin-averaged ISTS we dipped the tip into the Pt substrate until the telegraph signal on the out-of-plane clusters vanished and the spectrum taken on the Pt was feature-less within the voltage range used for ISTS. Then, ISTS was performed by positioning the tip on top of the GMC, adding a modulation voltage V mod = 80 μV (r.m.s.) of frequency f mod = 4142 Hz to V, stabilizing the tip at I stab = 2 nA and V stab = 5 mV, switching the feedback off, ramping the bias voltage and recording the \(\frac{{{\rm{d}}I}}{{{\rm{d}}V}}\) signal using a lock-in amplifier. For the investigation of the spin dynamics of the GMCs by SPSTM, we moved an out-of-plane spin-sensitive tip on top of the cluster, waited until the scanner creep was negligible, and recorded the tip height z as a function of time in constant-current mode5.
Ab initio DFT calculations
For the ab initio DFT calculations, we have used the KKR-GF method including spin-orbit coupling within the local-spin-density approximation46, 47. The Pt(111) hollow and top clusters have been modeled by a slab containing 24 and 64 Pt atoms, respectively, in both cases augmented by two vacuum regions. The relaxed distances of the Fe trimers toward the surface have been calculated by means of the QUANTUM-ESPRESSO package48. The MAE and the non-collinear energy landscape of the four types of trimers have been evaluated by band energy differences following the magnetic force theorem49, while the magnetic interactions among the Fe atoms have been computed using the generalized Liechtenstein formula50,51,52 and a fine mesh of 200 × 200 × 1 k points. For further details, we refer the reader to the Supplementary Note 2.
Master equation model
In order to simulate the GMC dynamics, we used the master equation model in conjunction with the effective spin Hamiltonian as described in ref. 5. Additionally, we took into account a bias voltage dependent change in the ratio between the coupling of the cluster to the tip, v tip, and to the surface, v surface, which we observed for the clusters investigated here. To this end, we assumed an exponential dependence on the bias voltage V via \({\left( {\frac{{{v_{{\rm{tip}}}}}}{{{v_{{\rm{surface}}}}}}} \right)_V} = {\left( {\frac{{{v_{{\rm{tip}}}}}}{{{v_{{\rm{surface}}}}}}} \right)_{V_T}} \times {{\rm e}^{ - \kappa \cdot \left( {V - {V_T}} \right)}}\), where \({\left( {\frac{{{v_{{\rm{tip}}}}}}{{{v_{{\rm{surface}}}}}}} \right)_{{V_T}}}\) is the ratio determined from the fit of the temperature dependence of the lifetime measured at voltage V T (Fig. 4f). The tip retraction parameter κ then follows from the fitting of the voltage dependence of the lifetime (Fig. 4g), resulting in \({\left( {\frac{{{v_{{\rm{tip}}}}}}{{{v_{{\rm{surface}}}}}}} \right)_{{V_T}}} = 0.085\) and κ = 0.15 for hcp top, \({\left( {\frac{{{v_{{\rm{tip}}}}}}{{{v_{{\rm{surface}}}}}}} \right)_{{V_T}}} \!\! = 0.095\) and κ = 0.1 for fcc hollow, and \({\left( {\frac{{{v_{{\rm{tip}}}}}}{{{v_{{\rm{surface}}}}}}} \right)_{{V_T}}} \!\!= 0.07\) and κ = 0.08 for the fcc top cluster. Since \(\frac{{{v_{{\rm{tip}}}}}}{{{v_{{\rm{surface}}}}}}\) additionally depends on the current, the fitted voltage dependence slightly deviates from the experimental data for the hcp top cluster. \(\frac{{{v_{{\rm{tip}}}}}}{{{v_{{\rm{surface}}}}}}\) varies between 0.01 and 0.022 for the field dependent asymmetry fits (Fig. 4e).
Data availability
The authors declare that the main data supporting the findings of this study are available within the article and its Supplementary Information files. Extra data are available from the corresponding author upon reasonable request.
References
Chappert, C., Fert, A. & Van Dau, F. N. The emergence of spin electronics in data storage. Nat. Mater. 6, 813–823 (2007).
Eigler, D. M. & Schweizer, E. K. Positioning single atoms with a scanning tunnelling microscope. Nature 344, 524–526 (1990).
Meier, F., Zhou, L., Wiebe, J. & Wiesendanger, R. Revealing magnetic interactions from single-atom magnetization curves. Science 320, 82–86 (2008).
Heinrich, A. J., Gupta, J. A., Lutz, C. P. & Eigler, D. M. Single-atom spin-flip spectroscopy. Science 306, 466–469 (2004).
Khajetoorians, A. A. et al. Current-driven spin dynamics of artificially constructed quantum magnets. Science 339, 55–59 (2013).
Loth, S., Baumann, S., Lutz, C. P., Eigler, D. M. & Heinrich, A. J. Bistability in atomic-scale antiferromagnets. Science 335, 196–199 (2012).
Khajetoorians, A. A., Wiebe, J., Chilian, B. & Wiesendanger, R. Realizing all-spin–based logic operations atom by atom. Science 332, 1062–1064 (2011).
Baumann, S. et al. Electron paramagnetic resonance of individual atoms on a surface. Science 350, 417–420 (2015).
Delgado, F. & Fernández-Rossier, J. Spin decoherence of magnetic atoms on surfaces. Prog. Surf. Sci. 92, 40–82 (2017).
Gambardella, P. et al. Giant magnetic anisotropy of single cobalt atoms and nanoparticles. Science 300, 1130–1133 (2003).
Gatteschi, D. & Sessoli, R. Molecular Nanomagnets (Oxford Univ. Press, 2006)
Kondo, J. Resistance minimum in dilute magnetic alloys. Prog. Theor. Phy. 32, 37–49 (1964).
Donati, F. et al. Magnetic remanence in single atoms. Science 352, 318–321 (2016).
Hübner, C., Baxevanis, B., Khajetoorians, A. A. & Pfannkuche, D. Symmetry effects on the spin switching of adatoms. Phys. Rev. B 90, 155134 (2014).
Singha, A. et al. Magnetic hysteresis in Er trimers on Cu(111). Nano Lett. 16, 3475–3481 (2016).
Balashov, T. et al. Dynamic magnetic excitations in 3d and 4f atoms and clusters. Surf. Sci. 630, 331–336 (2014).
Rau, I. G. et al. Reaching the magnetic anisotropy limit of a 3d metal atom. Science 344, 988–992 (2014).
Paul, W. et al. Control of the millisecond spin lifetime of an electrically probed atom. Nat. Phys. 13, 403–407 (2016).
Heinrich, B. W., Braun, L., Pascual, J. I. & Franke, K. J. Protection of excited spin states by a superconducting energy gap. Nat. Phys. 9, 765–768 (2013).
Khajetoorians, A. A. et al. Detecting excitation and magnetization of individual dopants in a semiconductor. Nature 467, 1084–1087 (2010).
Zhou, L. et al. Strength and directionality of surface Ruderman-Kittel-Kasuya-Yosida interaction mapped on the atomic scale. Nat. Phys. 6, 187–191 (2010).
Khajetoorians, A. A. et al. Tailoring the chiral magnetic interaction between two individual atoms. Nat. Commun. 7, 10620 (2016).
Błoński, P. et al. Magnetocrystalline anisotropy energy of Co and Fe adatoms on the (111) surfaces of Pd and Rh. Phys. Rev. B 81, 104426 (2010).
Oswald, A., Zeller, R. & Dederichs, P. H. Giant moments in palladium. Phys. Rev. Lett. 56, 1419–1422 (1986).
Šipr, O., Bornemann, S., Minár, J. & Ebert, H. Magnetic anisotropy of Fe and Co adatoms and monolayers: need for a proper treatment of the substrate. Phys. Rev. B 82, 174414 (2010).
Bouhassoune, M., Dias, Md. S., Zimmermann, B., Dederichs, P. H. & Lounis, S. RKKY-like contributions to the magnetic anisotropy energy: 3d adatoms on Pt(111) surface. Phys. Rev. B 94, 125402 (2016).
Meier, F. et al. Spin polarization of platinum (111) induced by the proximity to cobalt nanostripes. Phys. Rev. B 83, 075407 (2011).
Khajetoorians, A. A. et al. Spin excitations of individual Fe atoms on Pt(111): Impact of the site-dependent giant substrate polarization. Phys. Rev. Lett. 111, 157204 (2013).
Mydosh, J. A. Spin Glasses: An Experimental Introduction (Taylor & Francis Ltd, 1993)
Jamneala, T., Madhavan, V. & Crommie, M. F. Kondo response of a single antiferromagnetic chromium trimer. Phys. Rev. Lett. 87, 256804 (2001).
Krause, S., Berbil-Bautista, L., Herzog, G., Bode, M. & Wiesendanger, R. Current-induced magnetization switching with a spin-polarized scanning tunneling microscope. Science 317, 1537–1540 (2007).
Khajetoorians, A. A. et al. Itinerant nature of atom-magnetization excitation by tunneling electrons. Phys. Rev. Lett. 106, 037205 (2011).
Khajetoorians, A. A. et al. Tuning emergent magnetism in a Hund’s impurity. Nat. Nanotechnol. 10, 958–964 (2015).
Otte, A. F. et al. The role of magnetic anisotropy in the Kondo effect. Nat. Phys. 4, 847–850 (2008).
Madhavan, V., Chen, W., Jamneala, T., Crommie, M. F. & Wingreen, N. S. Tunneling into a single magnetic atom: Spectroscopic evidence of the Kondo resonance. Science 280, 567–569 (1998).
Nagaoka, K., Jamneala, T., Grobis, M. & Crommie, M. F. Temperature dependence of a single Kondo impurity. Phys. Rev. Lett. 88, 077205 (2002).
Zhang, Y.-h. et al. Temperature and magnetic field dependence of a Kondo system in the weak coupling regime. Nat. Commun. 4, 2110 (2013).
Suhl, H. Dispersion theory of the Kondo effect. Phys. Rev. 138, A515–A523 (1965).
Újsághy, O., Kroha, J., Szunyogh, L. & Zawadowski, A. Theory of the Fano resonance in the STM tunneling density of states due to a single Kondo impurity. Phys. Rev. Lett. 85, 2557–2560 (2000).
Costi, T. A. Kondo effect in a magnetic field and the magnetoresistivity of Kondo alloys. Phys. Rev. Lett. 85, 1504–1507 (2000).
Yan, S., Choi, D.-J., Burgess, J. A. J., Rolf-Pissarczyk, S. & Loth, S. Control of quantum magnets by atomic exchange bias. Nat. Nanotechnol. 10, 40–45 (2015).
Ibañez-Azpiroz, J., dos Santos Dias, M., Blügel, S. & Lounis, S. Zero-point spin-fluctuations of single adatoms. Nano Lett. 16, 4305–4311 (2016).
Brataas, A. & Hals, K. M. D. Spin-orbit torques in action. Nat. Nanotechnol. 9, 86–88 (2014).
Miron, I. M. et al. Perpendicular switching of a single ferromagnetic layer induced by in-plane current injection. Nature 476, 189–193 (2011).
Wiebe, J. et al. A 300 mk ultra-high vacuum scanning tunneling microscope for spin-resolved spectroscopy at high energy resolution. Rev. Sci. Instrum. 75, 4871–4879 (2004).
Papanikolaou, N., Zeller, R. & Dederichs, P. H. Conceptual improvements of the KKR method. J. Phys. Condens. Matter 14, 2799 (2002).
dos Santos Dias, M., Schweflinghaus, B., Blügel, S. & Lounis, S. Relativistic dynamical spin excitations of magnetic adatoms. Phys. Rev. B 91, 075405 (2015).
Paolo, G. et al. Quantum espresso: a modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 21, 395502 (2009).
Weinert, M., Watson, R. E. & Davenport, J. W. Total-energy differences and eigenvalue sums. Phys. Rev. B 32, 2115–2119 (1985).
Liechtenstein, A. I., Katsnelson, M. I. & Gubanov, V. A. Exchange interactions and spin-wave stiffness in ferromagnetic metals. J. Phys. F: Metal Phys. 14, L125 (1984).
Polesya, S. et al. Finite-temperature magnetism of FexPd1-x and CoxPt1-x alloys. Phys. Rev. B 82, 214409 (2010).
Udvardi, L., Szunyogh, L., Palotás, K. & Weinberger, P. First-principles relativistic study of spin waves in thin magnetic films. Phys. Rev. B 68, 104436 (2003).
Acknowledgements
We thank Kirsten von Bergmann for stimulating dicussions about non-collinear spin-processing schemes, Markus Ternes for valuable suggestions concerning the fitting of the Kondo resonances, and Roman Kovacic for technical help and discussions. J.H., A.S., M.S., R.W., and J.W. acknowledge funding from the SFB668 and the GrK1286 of the DFG. A.A.K. acknowledges financial support from the Emmy Noether Program (KH324/1-1) via the Deutsche Forschungsgemeinschaft, the Foundation of Fundamental Research on Matter (FOM), which is part of The Netherlands Organization for Scientific Research (NWO), and the VIDI project: ‘‘Manipulating the interplay between superconductivity and chiral magnetism at the single atom level’’ with project number 680-47-534 which is financed by NWO. J.I.-A., M.d.S.D., and S.L. acknowledge funding from Helmholtz Gemeinschaft Deutscher-Young Investigators Group Program No. VH-NG-717 (Functional Nanoscale Structure and Probe Simulation Laboratory), the Impuls und Vernetzungsfonds der Helmholtz-Gemeinschaft Postdoc Programme, and from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (ERC-consolidator grant 681405—DYNASORE).
Author information
Authors and Affiliations
Contributions
J.H., A.S., and J.W. designed the experiments. A.A.K. and M.S. helped with the design of the initial experiments. J.H., A.S., and K.T.T. carried out the measurements. J.H. and A.S. did the analysis of the experimental data. J.I.-A., M.d.S.D., and S.L. performed and analyzed the DFT calculations. Chr.H. and B.B. did the master equation modeling. J.H. and J.W. wrote the manuscript, to which all authors contributed via discussions and corrections.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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/.
About this article
Cite this article
Hermenau, J., Ibañez-Azpiroz, J., Hübner, C. et al. A gateway towards non-collinear spin processing using three-atom magnets with strong substrate coupling. Nat Commun 8, 642 (2017). https://doi.org/10.1038/s41467-017-00506-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-017-00506-7
This article is cited by
-
Trends in the hyperfine interactions of magnetic adatoms on thin insulating layers
npj Computational Materials (2021)
-
Magnetic resonance imaging of single atoms on a surface
Nature Physics (2019)
-
Stabilizing spin systems via symmetrically tailored RKKY interactions
Nature Communications (2019)
-
Spin-driven electrical power generation at room temperature
Communications Physics (2019)
-
Creating designer quantum states of matter atom-by-atom
Nature Reviews Physics (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.