Arxiv:1904.04922V1 [Cond-Mat.Mtrl-Sci] 9 Apr 2019

Total Page:16

File Type:pdf, Size:1020Kb

Arxiv:1904.04922V1 [Cond-Mat.Mtrl-Sci] 9 Apr 2019 First-principles investigation of spin-phonon coupling in vanadium-based molecular spin qubits 1Andrea Albino, 2Stefano Benci, 1Lorenzo Tesi, 3Matteo Atzori, 2Renato Torre, 4Stefano Sanvito, 1Roberta Sessoli,∗ and 4 Alessandro Lunghi† 1 Dipartimento di Chimica ”Ugo Schiff” and INSTM RU, Universit´adegli Studi di Firenze, I50019 Sesto Fiorentino, Italy 2 Dipartimento di Fisica ed Astronomia and European Laboratory for Nonlinear Spectroscopy, Universit´adegli Studi di Firenze, I50019 Sesto Fiorentino, Italy 3 Dipartimento di Chimica ”Ugo Schiff” and INSTM RU, Universit´adegli Studi di Firenze, I50019 Sesto Fiorentino, Italy. Present address: Laboratoire National des Champs Magnetiques Intenses, UPR 3228 - CNRS, F38042 Grenoble, France. and 4 School of Physics, CRANN and AMBER, Trinity College, Dublin 2, Ireland Paramagnetic molecules can show long spin-coherence times, which make them good candidates as quantum bits. Reducing the efficiency of the spin-phonon interaction is the primary challenge towards achieving long coherence times over a wide temperature range in soft molecular lattices. The lack of a microscopic understanding about the role of vibrations in spin relaxation strongly undermines the possibility to chemically design better performing molecular qubits. Here we report a first-principles charac- terization of the main mechanism contributing to the spin-phonon coupling for a class of vanadium(IV) molecular qubits. Post Hartree Fock and Density Functional Theory are used to determine the effect of both reticular and intra-molecular vibrations on the modulation of the Zeeman energy for four molecules showing different coordination geometries and ligands. This comparative study provides the first insight into the role played by coordination geometry and ligand field strength in determining the spin- lattice relaxation time of molecular qubits, opening the avenue to a rational design of new compounds. I. INTRODUCTION from the environment, while keeping the possibility to interact with each others. In practice three physical sys- Quantum information science deals with the represen- tems implement spin qubits: nitrogen-vacancy centres tation, storage and processing of information by means in diamond (NVC),[9] atomic impurities in semiconduc- of a quantum mechanical system. The basic element is tors, such as P implanted in Si,[10] and paramagnetic the quantum bit or qubit, namely the quantum analo- molecules.[11–13] In contrast to solid-state spin qubits gous of the classical bit (0, 1). Quantum computation based on dopant atoms, such as the NVCs, molecules exploits the quantum properties of the qubit, such as show a significant advantage, namely the chemical sys- superposition and entanglement, [1] providing an ideal tems hosting the spin can be tailored to tune the quan- platform for improving algorithms’ efficiency. In particu- tum properties and the coupling to other qubits. This can lar, one can design a range of quantum algorithms, which create quantum platforms,[8, 14, 15] providing a bottom- scale with the complexity of the problem in a much more up route to large-scale quantum register fabrication. favourable way than their classical counterparts. Several A spin lifetime at least 104 times longer than the time systems are currently investigated for the practical real- needed for an individual quantum operation[16, 17] rep- ization of quantum devices. For example: superconduct- resents the minimal requirement for the development of ing circuits, [2] trapped ions [3] and polarized photons.[4] a qubit. Accordingly, the figures of merit to consider in Among the various physical systems with potential for the design of electronic spin molecular qubits are funda- developing quantum technologies the spin, with its intrin- mentally two:[18] i) the longitudinal (or spin-lattice[19]) sic two-levels qubit structure, occupies a special place. relaxation time, T1, which corresponds to the lifetime arXiv:1904.04922v1 [cond-mat.mtrl-sci] 9 Apr 2019 Both nuclear spins[5] and electron spins,[6] as well as of a classical bit; ii) the coherence (or spin-spin relax- electron-nuclear hybrid systems,[7, 8] can be exploited for ation) time, T , which is the time characteristic for a this purpose. Electron spins interact more strongly with 2 spin to loose memory of a coherent quantum superpo- the environment, compared to the nuclear ones, and thus sition state. In the last few years, remarkable results they are easier to read out. In contrast, electrons have have been achieved against both T and T by investigat- shorter spin lifetimes and must be carefully protected 1 2 ing mononuclear transition-metal complexes.[13, 18, 20– 28] These results place molecules containing light met- als back in the quantum race. [29] Vanadium(IV)-based ∗ roberta.sessoli@unifi.it complexes represent a promising class of compounds to † [email protected] be used as fundamental components in quantum tech- 2 nologies. These spin 1/2 systems show long spin-spin lation found between the spin-phonon interaction ampli- relaxation time T2, up to 1 millisecond at liquid helium tudes and experimental spin-lattice relaxation times for temperature, when complexed with nuclear-spin free lig- the four compounds is then discussed. ands and diluted in nuclear-spin free solvents.[30] This property makes VIV-based compounds very attractive for further development. II. METHODS A common trend found in molecular spins[27, 31] is the rapid decrease of T2 on raising the temperature, a A. Vibrational properties calculations feature that limits their potential use at room temper- ature. The interaction with lattice vibrations, typically The simulations of the crystals’ vibrational properties connected to the T1-type relaxation, also contributes to are performed with a Gaussian and plane waves (GPW) the T2-type one and becomes the predominant relaxation formalism as implemented in the Quickstep module[38] mechanism when increasing temperature.[24] In solids, of the CP2K[39, 40] package. The GGA[41] functional in thermal motion is usually described by phonons, which the Perdew-Burke-Ernzerhof approximation[42] (PBE) are energy quanta of lattice vibrations. Spin-lattice re- is chosen for all the calculations. Van-der-Waals inter- laxation is caused by the absorption/emission of phonons actions are taken into account with the non local rVV10 by the spin system. This process is possible due to the correction scheme.[43] The calculation of the Hessian ma- presence of spin-orbit coupling,[32] an interaction that trix is performed by finite differences, the Hessian matrix couples atomic and spin degrees of freedom and enables is symmetrized averaging it with its transpose and the the energy exchange among the two systems. acoustic sum rule is applied to impose the translational Although the experimental investigations in this field invariance.[44] Each atomic coordinate is displaced by are numerous,[18, 26, 28, 31, 33] the theoretical descrip- 0.01A around the equilibrium configuration, and such tion is still at an early stage.[14, 34–37] In particular, displacements± are used to compute energies and forces. the possibility to include molecules in a solid ab initio When displacing the atoms we have taken into account computational framework has been made possible only the crystal’s symmetries and only inequivalent displace- in the last few years, enabled by an extensive work of in- ments have been considered (see Figure S1 in Supple- tegration of density functional theory (DFT), post Hatree mentary Information - SI). The exploitation of symme- Fock (postHF) methods and spin dynamics for the calcu- try allows us to reduce the computational overheads by lation of the dynamical magnetic properties of multi-spin a factor of four for all the systems considered. The diag- systems. onalization of the Hessian matrix provides the phonons The present work introduces a comparative theoret- frequencies, ωα (Table S1 in SI), and the normal modes ical investigation of the spin-lattice relaxation in four of vibrations, qα. The latter are stored in the columns of VIV molecular complexes (Figure 1). In particular, we the Hessian’s eigenvectors matrix L. have focused on penta-coordinated vanadyl (VO2+) and hexa-coordinated VIV molecules, where the coordination is obtained by cathecolate[31] and dithiolene[28] ligands, B. Spin-phonon coupling coefficients calculations namely [PPh4]2[VO(cat)2](1), [PPh4]2[V(cat)3](2), [PPh4]2[VO(dmit)2](3), [PPh4]2[V(dmit)3](4) (cat = The calculation of the spin and the spin-phonon catecholate, dmit = 1,3-dithiole-2-thione-4,5-dithiolate, Hamiltonian parameters has been carried out with the PPh4 = tetraphenylphosphonium). These complexes ORCA[45, 46] package. The level of theory used is Com- have already experimentally shown to possess long co- plete Active Space Self Consistent Field plus second order herence times and remarkable differences in the temper- perturbation theory (CASSCF+NEVPT2), with a def2- ature dependence of the spin-lattice relaxation time. Our TZVP basis set for V, O and S, whereas def2-SVP is used approach consists in modelling their magnetic properties for C and H. The active space includes one electron and by first-principles, when perturbing the molecular struc- the five d-orbitals of the molecule. The molecular geome- tures along the normal modes of vibration, following a try used for these simulations is obtained by the periodic strategy adopted in previous works.[34, 36] In order to DFT calculation for the optimized crystal cell. correctly account for the vibrational properties of solid The calculation of the spin-phonon coupling coeffi- state systems, where intermolecular interactions become cients is performed following a tensor differentiation relevant, DFT calculations are performed in the crys- procedure as described in a previous report on Single tal phase, in contrast to previous studies where the gas Molecule Magnets (SMMs).[36] This procedure is here phase was taken into consideration.[34] The quality of applied to the Land´e g tensor that describes the coupling our computed vibrational properties is ascertained by IR between the spins and an external magnetic field.
Recommended publications
  • Molecular Spintronics
    Molecular Spintronics Gabriel Aeppli 1, Andrew Fisher 1, Nicholas Harrison 2, Sandrine Heutz 3, Tim Jones 4, Chris Kay 5 and Des McMorrow 1 1 Department of Physics and Astronomy, London Centre for Nanotechnology, University College London, London WC1E 6BT, U.K. 2 Department of Chemistry, London Centre for Nanotechnology, Imperial College London, London SW7 2AZ, U.K. 3 Department of Materials, London Centre for Nanotechnology, Imperial College London, London SW7 2AZ, U.K. 4 Department of Chemistry, University of Warwick, Coventry, CV4 7AL, U.K. 5 Department of Biology, London Centre for Nanotechnology, University College London, London WC1E 6BT, UK. Project Organisation Project Summary and Context Combine cheap organic electronics and high performance spintronics • Funded through the Basic Technology programme to develop molecular spintronics with outcomes in IT and biosensing. • November 2008 start, duration 4 years. Use expertise in small molecule film growth, magnetism, theory, • Includes 3 institutions (Warwick, UCL and Imperial) and 7 investigators. optoelectronics, device engineering and spin resonance applied to • Crossing boundaries: PIs experts in different branches of Science (Chemistry, Physics, Biology) and Engineering (Materials, EE). biology. • Directly employs 4 PDRAs and 3 PhD students Molecular Electronics Spintronics • Project extends boundaries: Additional academics (a.o. Hirjibehedin, OPV, OLED, GMR, MRAM Curson, Nathan, Ryan), more than 7 PhD students and PDRAs closely Transistors Magnetic HJ Semicond. polymers Magnetic linked to the project and molecules Nelson, Durrant (IC), Forrest Baibich, PRL88 Semiconductors Organic Spintronics Molecular Magnetism Molecular films Magnetic switching, as tunnelling layers spin-crossover Visibility and Outcomes Molecules BT Molecular Molecular powder, on magn. surfaces Xiong, Nature 04 Spintronics Verdaguer, Science 96 e.g.
    [Show full text]
  • New Material for Quantum Computing Discovered out of the Blue 27 October 2013
    New material for quantum computing discovered out of the blue 27 October 2013 a significant advantage over similar materials that have been studied previously. Now, researchers from the London Centre for Nanotechnology at UCL and the University of British Columbia have shown that the electrons in CuPc can remain in 'superposition' – an intrinsically quantum effect where the electron exists in two states at once - for surprisingly long times, showing this simple dye molecule has potential as a medium for quantum technologies. The development of quantum computing requires precise control of tiny individual "qubits", the quantum analogs of the classical binary bits, '0' and '1', which underpin all of our computation and communications technologies today. What distinguishes the "qubits" from classical bits is their ability to exist in superposition states. The decay time of such superpositions tells us how useful a candidate qubit could be in quantum Phthalocyanine thin film on a flexible plastic substrate, technologies. If this time is long, quantum data showing the coexistence of long-lived "0" and "1" qubits storage, manipulation and transmission become on the copper spin. The molecules form a regular array possible. together with the metal-free analogues, and the background represents the lattice fringes of the Lead author Marc Warner from the London Centre molecular crystals obtained by transmission electron for Nanotechnology, said: "In theory, a quantum microscopy. Credit: Phil Bushell, Sandrine Heutz and computer can easily solve problems that a normal, Gabriel Aeppli classical, computer would not be able to answer in the lifetime of the universe. We just don't know how to build one yet.
    [Show full text]
  • Ultralong Copper Phthalocyanine Nanowires with New Crystal
    Ultralong Copper Phthalocyanine Nanowires with New Crystal Structure and Broad Optical Absorption Hai Wang 1, ◊, Soumaya Mauthoor 2, Salahud Din 2, Jules A. Gardener 1, †, Rio Chang 2, Marc Warner 1, Gabriel Aeppli 1, David W. McComb 2, Mary P. Ryan 2, Wei Wu 1, Andrew J. Fisher 1, A. Marshall Stoneham 1, Sandrine Heutz 2, * 1 Department of Physics and Astronomy and London Centre for Nanotechnology, University College London, London WC1E 6BT, UK. 2 Department of Materials and London Centre for Nanotechnology, Imperial College London, London SW7 2AZ, UK. ◊ Present Address: Department of Physics, Kunming University, Kunming, 650214, P. R. China. † Present Address: Department of Physics, Harvard University, Cambridge, MA 02138, USA. * e-mail: [email protected] The development of molecular nanostructures plays a major role in emerging organic electronic applications, as it leads to improved performance and is compatible with our increasing need for miniaturisation. In particular, nanowires have been obtained from solution or vapour phase and have displayed high conductivity1, 2 or large interfacial areas in solar cells. 3 In all cases however, the crystal structure remains as in films or bulk, and the exploitation of wires requires extensive post-growth manipulation as their orientations are random. Here we report copper phthalocyanine (CuPc) nanowires with diameters of 10-100 nm, high directionality and unprecedented aspect ratios. We demonstrate that they adopt a new crystal phase, designated η-CuPc, where the molecules stack along the long axis. The resulting high electronic overlap along the centimetre length stacks achieved in our wires mediates antiferromagnetic couplings and broadens the optical absorption spectrum.
    [Show full text]
  • C Copyright 2021 Torin F. Stetina
    c Copyright 2021 Torin F. Stetina Advancing Quantum Chemistry Simulations with Relativistic Electronic Structure and Quantum Computation Torin F. Stetina A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy University of Washington 2021 Reading Committee: Xiaosong Li, Chair Cody Schlenker Stefan Stoll Program Authorized to Offer Degree: Chemistry University of Washington Abstract Advancing Quantum Chemistry Simulations with Relativistic Electronic Structure and Quantum Computation Torin F. Stetina Chair of the Supervisory Committee: Harry and Catherine Jaynne Boand Endowed Professor Xiaosong Li Department of Chemistry This dissertation covers two main topics: relativistic Hamiltonians for quantum chemistry, and quantum computation. The main theme of the work presented within is the development of higher accuracy ab initio quantum chemistry simulations within the context of relativis- tic theory, and the theoretical analysis of quantum computational algorithms. In the first chapter, approximate density functionals are calibrated with a variational two-component relativistic Hamiltonian for the prediction of excited state X-ray absorption spectra. The second chapter presents a new method for computing four-component Dirac Hartree{Fock ground states of molecules, with an emphasis on its practical efficiency for large heavy ele- ment molecular systems in the future. The third and final chapter goes beyond the realm of Dirac Hartree{Fock theory, and analyzes how one could model the fully correlated molecu- lar Hamiltonian including second order effective quantum electrodynamic effects within the digital model of quantum computation. TABLE OF CONTENTS Page List of Figures . iii Glossary . vi Chapter 1: Theoretical Background and Conceptual Overview . .1 1.1 Electronic Structure Theory .
    [Show full text]
  • Department of Materials
    Department of Materials Annual Report and Research in Progress Annual Report and Research Department of Materials Imperial College London Exhibition Road London SW7 2AZ Department of Materials UK Annual Report and Tel: +44 (0)20 7594 6734 2009–10 Fax: +44 (0)20 7594 6757 Research in Progress 2009–10 www.imperial.ac.uk/materials 1 Department of Materials Research in Progress and Annual Report 2009–10 www.imperial.ac.uk/materials Contents The information in this publication was compiled by: 5 Introduction by the Head of 45 Undergraduate students 111 Institute for Biomedical Dagmar Durham and Sonia Tomasetig Engineering Department 45 Undergraduate courses Designed by: Helen Davison 112 Centre for Advanced Stuctural 46 Sources of support for Printed by: Shanleys 9 Annual Report 2009–10 Ceramics Copyright © June 2011 undergraduates 113 Institute for Security Science Department of Materials, Imperial College London 11 People 46 Industrial experience and work Technology 11 Our academic, research and placements All rights reserved. No part of this publication support staff 113 Materials Characterisation may be reproduced, stored in a retrieval system or 47 Final year undergraduate projects transmitted in any form or by any means, electronic, 15 Visiting researchers 49 Postgraduate school 117 International links mechanical, photocopying, recording or otherwise, 18 Our students 49 Postgraduate masters courses 117 KAUST without the permission of the publisher. 25 Summary of staff, visiting 50 Postgraduate research students 117 IDEA League Department
    [Show full text]
  • Exchange Interaction Between the Triplet Exciton and the Localized Spin
    Exchange interaction between the triplet exciton and the localized spin in copper-phthalocyanine Wei Wu∗ London Centre for Nanotechnology, University College London, Gower Street, London, WC1E 6BT, United Kingdom Triplet excitonic state in the organic molecule may arise from a singlet excitation and the following inter-system crossing. Especially for a spin-bearing molecule, an exchange interaction between the triplet exciton and the original spin on the molecule can be expected. In this paper, such 1 exchange interaction in copper-phthalocyanine (CuPc, spin- 2 ) was investigated from first-principles by using density-functional theory within a variety of approximations to the exchange correlation, ranging from local-density approximation to long-range corrected hybrid-exchange functional. The magnitude of the computed exchange interaction is in the order of meV with the minimum value (1.5 meV, ferromagnetic) given by the long-range corrected hybrid-exchange functional CAM-B3LYP. This exchange interaction can therefore give rise to a spin coherence with an oscillation period in the order of picoseconds, which is much shorter than the triplet lifetime in CuPc (typically tens of nanoseconds). This implies that it might be possible to manipulate the localised spin on Cu experimentally using optical excitation and inter-system crossing well before the triplet state disappears. PACS numbers: 71.10.Aw,71.15.Mb,71.35.Gg,71.70.Gm I. INTRODUCTION nipulate the exchange interaction by external stimulus such as optical excitation. The metal-phthalocyanine (MPc), in which a metallic The optically excited state in molecule could not only ion is centred at a conjugated Pc ring, has attracted much facilitate the control of interaction in QIP [13] but also attention recently due to its fascinating physical prop- play an important role in organic spintronics [14, 15].
    [Show full text]