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.

Spin dynamics of current-driven single magnetic adatoms and molecules

A scanning tunneling microscope can probe the inelastic spin excitations of a single magnetic atom in a surface via spin-flip assisted tunneling in which transport electrons exchange spin and energy with the atomic spin. If the inelastic transport time, defined as the average time elapsed between two inelastic spin flip events, is shorter than the atom spin-relaxation time, the scanning tunnel microscope (STM) current can drive the spin out of equilibrium. Here we model this process using rate equations and a model Hamiltonian that describes successfully spin-flip-assisted tunneling experiments, including a single Mn atom, a Mn dimer, and Fe Phthalocyanine molecules. When the STM current is not spin polarized, the nonequilibrium spin dynamics of the magnetic atom results in nonmonotonic dI/dV curves. In the case of spin-polarized STM current, the spin orientation of the magnetic atom can be controlled parallel or antiparallel to the magnetic moment of the tip. Thus, spin-polarized STM tips can be used both to probe and to control the magnetic moment of a single atom.

Cotunneling theory of atomic spin inelastic electron tunneling spectroscopy

We propose cotunneling as the microscopic mechanism that makes possible inelastic electron tunneling spectroscopy of magnetic atoms in surfaces for a wide range of systems, including single magnetic adatoms, molecules, and molecular stacks. We describe electronic transport between the scanning tip and the conducting surface through the magnetic system (MS) with a generalized Anderson model, without making use of effective spin models. Transport and spin dynamics are described with an effective cotunneling Hamiltonian in which the correlations in the magnetic system are calculated exactly and the coupling to the electrodes is included up to second order in the tip MS and MS substrate. In the adequate limit our approach is equivalent to the phenomenological Kondo exchange model that successfully describes the experiments. We apply our method to study in detail inelastic transport in two systems, stacks of cobalt phthalocyanines and a single Mn atom on Cu2N. Our method accounts for both the large contribution of the inelastic spin exchange events to the conductance and the observed conductance asymmetry.

Theory of single-spin inelastic tunneling spectroscopy

I show that recent experiments of inelastic scanning tunneling spectroscopy of single and a few magnetic atoms are modeled with a phenomenological spin-assisted tunneling Hamiltonian so that the inelastic dI/dV line shape is related to the spin spectral weight of the magnetic atom. This accounts for the spin selection rules and dI/dV spectra observed experimentally for single Fe and Mn atoms deposited on Cu2N. In the case of chains of Mn atoms it is found necessary to include both first and second-neighbor exchange interactions as well as single-ion anisotropy.

Efficient intrinsic spin relaxation in graphene due to flexural distortions

We propose an intrinsic spin scattering mechanism in graphene originated by the interplay of atomic spin-orbit interaction and the local curvature induced by flexural distortions of the atomic lattice. Starting from a multiorbital tight-binding Hamiltonian with spin-orbit coupling considered non-perturbatively, we derive an effective Hamiltonian for the spin scattering of the Dirac electrons due to flexural distortions. We compute the spin lifetime due to both flexural phonons and ripples and we find values in the 1-10 ns range at room temperature. The proposed mechanism dominates the spin relaxation in high mobility graphene samples and should also apply to other planar aromatic compounds.

Spin alignment of extra electrons in K-phenanthrene clusters taken from the crystalline tripotassium-intercalated phenanthrene structure

The appearance of ferromagnetic correlations among π electrons of phenanthrene (C14H10) molecules in the herringbone structure is proven for K doped clusters both by ab initio quantum-chemistry calculations and by the direct solution of the many-body Pariser-Parr-Pople Hamiltonian. Magnetic ground states are predicted for one or three additional electrons per phenanthrene molecule. These results are a consequence of the small overlap between the lowest unoccupied molecular orbitals (and lowest unoccupied molecular orbitals + 1) of neutral neighboring phenanthrene molecules, which makes the gain in energy by delocalization similar to the corresponding increase due to the Coulomb interaction.

Anisotropic intrinsic spin relaxation in graphene due to flexural distortions

We propose an intrinsic spin scattering mechanism in graphene originated by the interplay of atomic spin-orbit interaction and the local curvature induced by flexural distortions of the atomic lattice. Starting from a multiorbital tight-binding Hamiltonian with spin-orbit coupling considered nonperturbatively, we derive an effective Hamiltonian for the spin scattering of the Dirac electrons due to flexural distortions. We compute the spin lifetime due to both flexural phonons and ripples and we find values in the microsecond range at room temperature. Interestingly, this mechanism is anisotropic on two counts. First, the relaxation rate is different for off-plane and in-plane spin quantization axis. Second, the spin relaxation rate depends on the angle formed by the crystal momentum with the carbon-carbon bond. In addition, the spin lifetime is also valley dependent. The proposed mechanism sets an upper limit for spin lifetimes in graphene and will be relevant when samples of high quality can be fabricated free of extrinsic sources of spin relaxation.

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.