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.

On the forbidden gap of finite graphene nanoribbons

The electronic structure of isolated finite graphene nanoribbons is investigated by solving, at the Hartree-Fock (HF) level, the Pariser, Parr and Pople (PPP) many-body Hamiltonian. The study is mainly focused on 7-AGNR and 13-AGNR (Armchair Graphene Nano-Ribbons), whose electronic structures have been recently experimentally investigated. Only paramagnetic solutions are considered. The characteristics of the forbidden gap are studied as a function of the ribbon length. For a 7-AGNR, the gap monotonically decreases from a maximum value of ~6.5 eV for short nanoribbons to a very small value of ~0.12 eV for the longer calculated systems. Gap edges are defined by molecular orbitals that are spatially localized near the nanoribbon extremes, that is, near both zig-zag edges. On the other hand, two delocalized orbitals define a much larger gap of about 5 eV. Conductance measurements report a somewhat smaller gap of ~3 eV. The small real gap lies in the middle of the one given by extended states and has been observed by STM and reproduced by DFT calculations. On the other hand, the length dependence of the gap is not monotonous for a 13-AGNR. It decreases initially but sharply increases for lengths beyond 30 Å remaining almost constant thereafter at a value of ~2.1 eV. Two additional states localized at the nanoribbon extremes show up at energies 0.31 eV below the HOMO (Highest Occupied Molecular Orbital) and above the LUMO (Lowest Unoccupied Molecular Orbital). These numbers compare favorably with those recently obtained by means of STS for a 13-AGNR sustained by a gold surface, namely 1.4 eV for the energy gap and 0.4 eV for the position of localized band edges. We show that the important differences between 7- and 13-AGNR should be ascribed to the charge rearrangement near the zig-zag edges obtained in our calculations for ribbons longer than 30 Å, a feature that does not show up for a 7-AGNR no matter its length.

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.

Earth’s Rotation: A Challenging Problem in Mathematics and Physics

A suitable knowledge of the orientation and motion of the Earth in space is a common need in various fields. That knowledge has been ever necessary to carry out astronomical observations, but with the advent of the space age, it became essential for making observations of satellites and predicting and determining their orbits, and for observing the Earth from space as well. Given the relevant role it plays in Space Geodesy, Earth rotation is considered as one of the three pillars of Geodesy, the other two being geometry and gravity. Besides, research on Earth rotation has fostered advances in many fields, such as Mathematics, Astronomy and Geophysics, for centuries. One remarkable feature of the problem is in the extreme requirements of accuracy that must be fulfilled in the near future, about a millimetre on the tangent plane to the planet surface, roughly speaking. That challenges all of the theories that have been devised and used to-date; the paper makes a short review of some of the most relevant methods, which can be envisaged as milestones in Earth rotation research, emphasizing the Hamiltonian approach developed by the authors. Some contemporary problems are presented, as well as the main lines of future research prospected by the International Astronomical Union/International Association of Geodesy Joint Working Group on Theory of Earth Rotation, created in 2013.

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.

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.

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.

Towards Accurate and Efficient Description of Excited States

Thesis (Ph.D.)--University of Washington, 2016-06

Quantum transport and integrability of the Anderson model for a quantum dot with multiple leads

We show that an Anderson Hamiltonian describing a quantum dot connected to multiple leads is integrable. A general expression for the nonlinear conductance is obtained by combining the Bethe ansatz exact solution with Landauer-Buttiker theory. In the Kondo regime, a closed form expression is given for the matrix conductance at zero temperature and when all the leads are close to the symmetric point. A bias-induced splitting of the Kondo resonance is possible for three or more leads. Specifically, for N leads, with each at a different chemical potential, there can be N-1 Kondo peaks in the conductance.

Lower bounds on the complexity of simulating quantum gates

We give a simple proof of a formula for the minimal time required to simulate a two-qubit unitary operation using a fixed two-qubit Hamiltonian together with fast local unitaries. We also note that a related lower bound holds for arbitrary n-qubit gates.

Nonlinear dynamics of a satellite with deployable solar panel arrays

The multibody dynamics of a satellite in circular orbit, modeled as a central body with two hinge-connected deployable solar panel arrays, is investigated. Typically, the solar panel arrays are deployed in orbit using preloaded torsional springs at the hinges in a near symmetrical accordion manner, to minimize the shock loads at the hinges. There are five degrees of freedom of the interconnected rigid bodies, composed of coupled attitude motions (pitch, yaw and roll) of the central body plus relative rotations of the solar panel arrays. The dynamical equations of motion of the satellite system are derived using Kane's equations. These are then used to investigate the dynamic behavior of the system during solar panel deployment via the 7-8th-order Runge-Kutta integration algorithms and results are compared with approximate analytical solutions. Chaotic attitude motions of the completely deployed satellite in circular orbit under the influence of the gravity-gradient torques are subsequently investigated analytically using Melnikov's method and confirmed via numerical integration. The Hamiltonian equations in terms of Deprit's variables are used to facilitate the analysis. (C) 2003 Published by Elsevier Ltd.

The effects of J-gate potential and interfaces on donor exchange coupling in the Kane quantum computer architecture

We calculate the electron exchange coupling for a phosphorus donor pair in silicon perturbed by a J-gate potential and the boundary effects of the silicon host geometry. In addition to the electron-electron exchange interaction we also calculate the contact hyperfine interaction between the donor nucleus and electron as a function of the varying experimental conditions. Donor separation, depth of the P nuclei below the silicon oxide layer and J-gate voltage become decisive factors in determining the strength of both the exchange coupling and hyperfine interaction-both crucial components for qubit operations in the Kane quantum computer. These calculations were performed using an anisotropic effective-mass Hamiltonian approach. The behaviour of the donor exchange coupling as a function of the parameters varied in this work provides relevant information for the experimental design of these devices.