Evidence for non-conservative current-induced forces in the breaking of Au and Pt atomic chains

This experimental work aims at probing current-induced forces at the atomic scale. Specifically it addresses predictions in recent work regarding the appearance of run-away modes as a result of a combined effect of the non-conservative wind force and a ‘Berry force’. The systems we consider here are atomic chains of Au and Pt atoms, for which we investigate the distribution of break down voltage values. We observe two distinct modes of breaking for Au atomic chains. The breaking at high voltage appears to behave as expected for regular break down by thermal excitation due to Joule heating. However, there is a low-voltage breaking mode that has characteristics expected for the mechanism of current-induced forces. Although a full comparison would require more detailed information on the individual atomic configurations, the systems we consider are very similar to those considered in recent model calculations and the comparison between experiment and theory is very encouraging for the interpretation we propose.

Can model Hamiltonians describe the electron–electron interaction in π-conjugated systems?: PAH and graphene

Model Hamiltonians have been, and still are, a valuable tool for investigating the electronic structure of systems for which mean field theories work poorly. This review will concentrate on the application of Pariser–Parr–Pople (PPP) and Hubbard Hamiltonians to investigate some relevant properties of polycyclic aromatic hydrocarbons (PAH) and graphene. When presenting these two Hamiltonians we will resort to second quantisation which, although not the way chosen in its original proposal of the former, is much clearer. We will not attempt to be comprehensive, but rather our objective will be to try to provide the reader with information on what kinds of problems they will encounter and what tools they will need to solve them. One of the key issues concerning model Hamiltonians that will be treated in detail is the choice of model parameters. Although model Hamiltonians reduce the complexity of the original Hamiltonian, they cannot be solved in most cases exactly. So, we shall first consider the Hartree–Fock approximation, still the only tool for handling large systems, besides density functional theory (DFT) approaches. We proceed by discussing to what extent one may exactly solve model Hamiltonians and the Lanczos approach. We shall describe the configuration interaction (CI) method, a common technology in quantum chemistry but one rarely used to solve model Hamiltonians. In particular, we propose a variant of the Lanczos method, inspired by CI, that has the novelty of using as the seed of the Lanczos process a mean field (Hartree–Fock) determinant (the method will be named LCI). Two questions of interest related to model Hamiltonians will be discussed: (i) when including long-range interactions, how crucial is including in the Hamiltonian the electronic charge that compensates ion charges? (ii) Is it possible to reduce a Hamiltonian incorporating Coulomb interactions (PPP) to an 'effective' Hamiltonian including only on-site interactions (Hubbard)? The performance of CI will be checked on small molecules. The electronic structure of azulene and fused azulene will be used to illustrate several aspects of the method. As regards graphene, several questions will be considered: (i) paramagnetic versus antiferromagnetic solutions, (ii) forbidden gap versus dot size, (iii) graphene nano-ribbons, and (iv) optical properties.

Dynamic bonding of metallic nanocontacts: Insights from experiments and atomistic simulations

The conductance across an atomically narrow metallic contact can be measured by using scanning tunneling microscopy. In certain situations, a jump in the conductance is observed right at the point of contact between the tip and the surface, which is known as “jump to contact” (JC). Such behavior provides a way to explore, at a fundamental level, how bonding between metallic atoms occurs dynamically. This phenomenon depends not only on the type of metal but also on the geometry of the two electrodes. For example, while some authors always find JC when approaching two atomically sharp tips of Cu, others find that a smooth transition occurs when approaching a Cu tip to an adatom on a flat surface of Cu. In an attempt to show that all these results are consistent, we make use of atomistic simulations; in particular, classical molecular dynamics together with density functional theory transport calculations to explore a number of possible scenarios. Simulations are performed for two different materials: Cu and Au in a [100] crystal orientation and at a temperature of 4.2 K. These simulations allow us to study the contribution of short- and long-range interactions to the process of bonding between metallic atoms, as well as to compare directly with experimental measurements of conductance, giving a plausible explanation for the different experimental observations. Moreover, we show a correlation between the cohesive energy of the metal, its Young's modulus, and the frequency of occurrence of a jump to contact.

Large spin splitting in the conduction band of transition metal dichalcogenide monolayers

We study the conduction band spin splitting that arises in transition metal dichalcogenide (TMD) semiconductor monolayers such as MoS2, MoSe2, WS2, and WSe2 due to the combination of spin-orbit coupling and lack of inversion symmetry. Two types of calculation are done. First, density functional theory (DFT) calculations based on plane waves that yield large splittings, between 3 and 30 meV. Second, we derive a tight-binding model that permits to address the atomic origin of the splitting. The basis set of the model is provided by the maximally localized Wannier orbitals, obtained from the DFT calculation, and formed by 11 atomiclike orbitals corresponding to d and p orbitals of the transition metal (W, Mo) and chalcogenide (S, Se) atoms respectively. In the resulting Hamiltonian, we can independently change the atomic spin-orbit coupling constant of the two atomic species at the unit cell, which permits to analyze their contribution to the spin splitting at the high symmetry points. We find that—in contrast to the valence band—both atoms give comparable contributions to the conduction band splittings. Given that these materials are most often n-doped, our findings are important for developments in TMD spintronics.

Quantum fluctuations stabilize skyrmion textures

We study the quantum spin waves associated to skyrmion textures. We show that the zero-point energy associated to the quantum spin fluctuations of a noncollinear spin texture produce Casimir-like magnetic fields. We study the effect of these Casimir fields on the topologically protected noncollinear spin textures known as skyrmions. In a Heisenberg model with Dzyalonshinkii-Moriya interactions, chosen so the classical ground state displays skyrmion textures, we calculate the spin-wave spectrum, using the Holstein-Primakoff approximation, and the associated zero-point energy, to the lowest order in the spin-wave expansion. Our calculations are done both for the single-skyrmion case, for which we obtain a discrete set of skyrmion bound states, as well as for the skyrmion crystal, for which the resulting spectrum gives the spin-wave bands. In both cases, our calculations show that the Casimir magnetic field contributes up to 10% of the total Zeeman energy necessary to delete the skyrmion texture with an applied field.

Electronic properties of transition metal atoms on Cu2N/Cu(100)

We study the nature of spin excitations of individual transition metal atoms (Ti, V, Cr, Mn, Fe, Co, and Ni) deposited on a Cu2N/Cu(100) surface using both spin-polarized density functional theory (DFT) and exact diagonalization of an Anderson model derived from DFT. We use DFT to compare the structural, electronic, and magnetic properties of different transition metal adatoms on the surface. We find that the average occupation of the transition metal d shell, main contributor to the magnetic moment, is not quantized, in contrast with the quantized spin in the model Hamiltonians that successfully describe spin excitations in this system. In order to reconcile these two pictures, we build a zero bandwidth multi-orbital Anderson Hamiltonian for the d shell of the transition metal hybridized with the p orbitals of the adjacent nitrogen atoms, by means of maximally localized Wannier function representation of the DFT Hamiltonian. The exact solutions of this model have quantized total spin, without quantized charge at the d shell. We propose that the quantized spin of the models actually belongs to many-body states with two different charge configurations in the d shell, hybridized with the p orbital of the adjacent nitrogen atoms. This scenario implies that the measured spin excitations are not fully localized at the transition metal.

Noncollinear magnetic phases and edge states in graphene quantum Hall bars

Application of a perpendicular magnetic field to charge neutral graphene is expected to result in a variety of broken symmetry phases, including antiferromagnetic, canted, and ferromagnetic. All these phases open a gap in bulk but have very different edge states and noncollinear spin order, recently confirmed experimentally. Here we provide an integrated description of both edge and bulk for the various magnetic phases of graphene Hall bars making use of a noncollinear mean field Hubbard model. Our calculations show that, at the edges, the three types of magnetic order are either enhanced (zigzag) or suppressed (armchair). Interestingly, we find that preformed local moments in zigzag edges interact with the quantum spin Hall like edge states of the ferromagnetic phase and can induce backscattering.

Quantum theory of spin waves in finite chiral spin chains

We calculate the effect of spin waves on the properties of finite-size spin chains with a chiral spin ground state observed on biatomic Fe chains deposited on iridium(001). The system is described with a Heisenberg model supplemented with a Dzyaloshinskii-Moriya coupling and a uniaxial single ion anisotropy that presents a chiral spin ground state. Spin waves are studied using the Holstein-Primakoff boson representation of spin operators. Both the renormalized ground state and the elementary excitations are found by means of Bogoliubov transformation, as a function of the two variables that can be controlled experimentally, the applied magnetic field and the chain length. Three main results are found. First, because of the noncollinear nature of the classical ground state, there is a significant zero-point reduction of the ground-state magnetization of the spin spiral. Second, there is a critical external field from which the ground state changes from chiral spin ground state to collinear ferromagnetic order. The character of the two lowest-energy spin waves changes from edge modes to confined bulk modes over this critical field. Third, in the spin-spiral state, the spin-wave spectrum exhibits oscillatory behavior as function of the chain length with the same period of the spin helix.

Majorana Zero Modes in Graphene

A clear demonstration of topological superconductivity (TS) and Majorana zero modes remains one of the major pending goals in the field of topological materials. One common strategy to generate TS is through the coupling of an s-wave superconductor to a helical half-metallic system. Numerous proposals for the latter have been put forward in the literature, most of them based on semiconductors or topological insulators with strong spin-orbit coupling. Here, we demonstrate an alternative approach for the creation of TS in graphene-superconductor junctions without the need for spin-orbit coupling. Our prediction stems from the helicity of graphene’s zero-Landau-level edge states in the presence of interactions and from the possibility, experimentally demonstrated, of tuning their magnetic properties with in-plane magnetic fields. We show how canted antiferromagnetic ordering in the graphene bulk close to neutrality induces TS along the junction and gives rise to isolated, topologically protected Majorana bound states at either end. We also discuss possible strategies to detect their presence in graphene Josephson junctions through Fraunhofer pattern anomalies and Andreev spectroscopy. The latter, in particular, exhibits strong unambiguous signatures of the presence of the Majorana states in the form of universal zero-bias anomalies. Remarkable progress has recently been reported in the fabrication of the proposed type of junctions, which offers a promising outlook for Majorana physics in graphene systems.

Magnetic Edge Anisotropy in Graphenelike Honeycomb Crystals

The independent predictions of edge ferromagnetism and the quantum spin Hall phase in graphene have inspired the quest of other two-dimensional honeycomb systems, such as silicene, germanene, stanene, iridates, and organometallic lattices, as well as artificial superlattices, all of them with electronic properties analogous to those of graphene, but a larger spin-orbit coupling. Here, we study the interplay of ferromagnetic order and spin-orbit interactions at the zigzag edges of these graphenelike systems. We find an in-plane magnetic anisotropy that opens a gap in the otherwise conducting edge channels that should result in large changes of electronic properties upon rotation of the magnetization.

Imaging of spin waves in atomically designed nanomagnets

The spin dynamics of all ferromagnetic materials are governed by two types of collective phenomenon: spin waves and domain walls. The fundamental processes underlying these collective modes, such as exchange interactions and magnetic anisotropy, all originate at the atomic scale. However, conventional probing techniques based on neutron1 and photon scattering2 provide high resolution in reciprocal space, and thereby poor spatial resolution. Here we present direct imaging of standing spin waves in individual chains of ferromagnetically coupled S = 2 Fe atoms, assembled one by one on a Cu2N surface using a scanning tunnelling microscope. We are able to map the spin dynamics of these designer nanomagnets with atomic resolution in two complementary ways. First, atom-to-atom variations of the amplitude of the quantized spin-wave excitations are probed using inelastic electron tunnelling spectroscopy. Second, we observe slow stochastic switching between two opposite magnetization states3, 4, whose rate varies strongly depending on the location of the tip along the chain. Our observations, combined with model calculations, reveal that switches of the chain are initiated by a spin-wave excited state that has its antinodes at the edges of the chain, followed by a domain wall shifting through the chain from one end to the other. This approach opens the way towards atomic-scale imaging of other types of spin excitation, such as spinon pairs and fractional end-states5, 6, in engineered spin chains.

Derivation of the spin Hamiltonians for Fe in MgO

A method to calculate the effective spin Hamiltonian for a transition metal impurity in a non-magnetic insulating host is presented and applied to the paradigmatic case of Fe in MgO. In the first step we calculate the electronic structure employing standard density functional theory (DFT), based on generalized gradient approximation (GGA), using plane waves as a basis set. The corresponding basis of atomic-like maximally localized Wannier functions is derived and used to represent the DFT Hamiltonian, resulting in a tight-binding model for the atomic orbitals of the magnetic impurity. The third step is to solve, by exact numerical diagonalization, the N electron problem in the open shell of the magnetic atom, including both effects of spin–orbit and Coulomb repulsion. Finally, the low energy sector of this multi-electron Hamiltonian is mapped into effective spin models that, in addition to the spin matrices S, can also include the orbital angular momentum L when appropriate. We successfully apply the method to Fe in MgO, considering both the undistorted and Jahn–Teller (JT) distorted cases. Implications for the influence of Fe impurities on the performance of magnetic tunnel junctions based on MgO are discussed.

Edge states in graphene-like systems

The edges of graphene and graphene like systems can host localized states with evanescent wave function with properties radically different from those of the Dirac electrons in bulk. This happens in a variety of situations, that are reviewed here. First, zigzag edges host a set of localized non-dispersive state at the Dirac energy. At half filling, it is expected that these states are prone to ferromagnetic instability, causing a very interesting type of edge ferromagnetism. Second, graphene under the influence of external perturbations can host a variety of topological insulating phases, including the conventional quantum Hall effect, the quantum anomalous Hall (QAH) and the quantum spin Hall phase, in all of which phases conduction can only take place through topologically protected edge states. Here we provide an unified vision of the properties of all these edge states, examined under the light of the same one orbital tight-binding model. We consider the combined action of interactions, spin–orbit coupling and magnetic field, which produces a wealth of different physical phenomena. We briefly address what has been actually observed experimentally.

The emergence of classical behaviour in magnetic adatoms

A wide class of nanomagnets shows striking quantum behaviour, known as quantum spin tunnelling (QST): instead of two degenerate ground states with opposite magnetizations, a bonding-antibonding pair forms, resulting in a splitting of the ground-state doublet with wave functions linear combination of two classically opposite magnetic states, leading to the quenching of their magnetic moment. Here we study how QST is destroyed and classical behaviour emerges in the case of magnetic adatoms, where, contrary to larger nanomagnets, the QST splitting is in some instances bigger than temperature and broadening. We analyze two different mechanisms for the renormalization of the QST splitting: Heisenberg exchange between different atoms, and Kondo exchange interaction with the substrate electrons. Sufficiently strong spin-substrate and spin-spin coupling renormalize the QST splitting to zero allowing the environmental decoherence to eliminate superpositions between classical states, leading to the emergence of spontaneous magnetization. Importantly, we extract the strength of the Kondo exchange for various experiments on individual adatoms and construct a phase diagram for the classical to quantum transition.

Controlled Complete Suppression of Single-Atom Inelastic Spin and Orbital Cotunneling

The inelastic portion of the tunnel current through an individual magnetic atom grants unique access to read out and change the atom’s spin state, but it also provides a path for spontaneous relaxation and decoherence. Controlled closure of the inelastic channel would allow for the latter to be switched off at will, paving the way to coherent spin manipulation in single atoms. Here, we demonstrate complete closure of the inelastic channels for both spin and orbital transitions due to a controlled geometric modification of the atom’s environment, using scanning tunneling microscopy (STM). The observed suppression of the excitation signal, which occurs for Co atoms assembled into chains on a Cu2N substrate, indicates a structural transition affecting the dz2 orbital, effectively cutting off the STM tip from the spin-flip cotunneling path.