Arxiv:2009.11417V2 [Quant-Ph] 14 Jan 2021

Total Page:16

File Type:pdf, Size:1020Kb

Arxiv:2009.11417V2 [Quant-Ph] 14 Jan 2021 A state-averaged orbital-optimized hybrid quantum-classical algorithm for a democratic description of ground and excited states Saad Yalouz,1, 2, ∗ Bruno Senjean,2, 1, y Jakob G¨unther,1 Francesco Buda,3 Thomas E. O'Brien,4, 2 and Lucas Visscher1 1Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV, Amsterdam, The Netherlands 2Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 3Leiden Institute of Chemistry, Leiden University, Einsteinweg 55, P.O. Box 9502, 2300 RA Leiden, The Netherlands. 4Google Research, 80636 Munich, Germany In the Noisy Intermediate-Scale Quantum (NISQ) era, solving the electronic structure problem from chemistry is considered as the \killer application" for near-term quantum devices. In spite of the success of variational hybrid quantum/classical algorithms in providing accurate energy pro- files for small molecules, careful considerations are still required for the description of complicated features of potential energy surfaces. Because the current quantum resources are very limited, it is common to focus on a restricted part of the Hilbert space (determined by the set of active orbitals). While physically motivated, this approximation can severely impact the description of these complicated features. A perfect example is that of conical intersections (i.e. a singular point of degeneracy between electronic states), which are of primary importance to understand many prominent reactions. Designing active spaces so that the improved accuracy from a quantum com- puter is not rendered useless is key to finding useful applications of these promising devices within the field of chemistry. To answer this issue, we introduce a NISQ-friendly method called \State- Averaged Orbital-Optimized Variational Quantum Eigensolver" (SA-OO-VQE) which combines two algorithms: (1) a state-averaged orbital-optimizer, and (2) a state-averaged VQE. To demonstrate the success of the method, we classically simulate it on a minimal Schiff base model (namely the formaldimine molecule CH2NH) relevant also for the photoisomerization in rhodopsin | a crucial step in the process of vision mediated by the presence of a conical intersection. We show that merg- ing both algorithms fulfil the necessary condition to describe the molecule's conical intersection, i.e. the ability to treat degenerate (or quasi-degenerate) states on the same footing. I. INTRODUCTION a chemical reaction), several chemical reactions also de- pend on excited states (see for instance Ref. 21 and refer- Quantum computing promises to solve problems that ences therein). To date, most variational quantum algo- are beyond the capacities of classical devices. Technolog- rithms have been specifically designed to prepare ground ical advances in quantum computers occur rapidly [1{8], states [1, 22{37] while excited state calculations have re- and it is now of primary importance to develop quantum ceived less attention until only recently [23, 38{49]. It algorithms dedicated to specific tasks of high industrial should be noted that some excited states can also be ob- and societal impact. The electronic structure problem in tained as the lowest energy solutions of a given set of chemistry is considered as one of the killer applications symmetries [50{55]. of quantum computers [9{14], and it remains an open The landscape of ground- and excited-state potential question whether these applications will be achieved in energy surfaces (PES) can be quite complicated, making the noisy intermediate-scale quantum (NISQ) era [15]. even qualitative numerical description challenging. As While quantum phase estimation algorithms [16{18] can atomic nuclei move, ground- and excited-state energies in principle solve the electronic structure problem in its may tend towards each other, leading to either avoided entirety, they have received limited focus in the NISQ era crossings or true inversions in the ordering of the elec- due to their apparent longer circuit depth requirements. tronic states. These inversions are referred to as conical arXiv:2009.11417v2 [quant-ph] 14 Jan 2021 Instead, variational hybrid quantum/classical algorithms intersections and are ubiquitous in photochemistry and require much shallower circuits at the cost of additional photobiology [56, 57]. Their presence mediates impor- measurements, with as a spearhead the variational quan- tant reactions, for instance the photoisomerisation of the tum eigensolver (VQE) [7, 19, 20]. Those NISQ-friendly chromophore of the visual photoreceptor rhodopsin, en- algorithms are hoped to solve quantum chemical prob- suing absorption of light. Photoisomerisation is a crucial lems of moderate-large size in the near future. step in the process of vision [58{62]. Conical intersections While computing the ground state of a molecule gives also find their importance in photosynthesis [63], photo- information about the equilibrium geometries and tran- stabilization of DNA [64{67] and excitation-energy trans- sition states (essential to determine the energy barrier of fer [68] such as in supramolecular light-harvesting anten- nae [69]. Despite their indisputable importance, their description remains a challenge in the chemistry commu- nity [56, 57, 70]. Diagonalization of the electronic Hamil- ∗ [email protected] tonian matrix being intractable on classical devices, it is y [email protected] 2 of common use to reduce the Hilbert space size by consid- II. THEORY ering only a subspace of (so-called active) orbitals of the full problem. However, such a truncated Hilbert space is A. Electronic structure problem not guaranteed to be accurate enough to represent a con- ical intersection anymore. A possible solution is provided In quantum chemistry, the Born-Oppenheimer approx- by state-averaged orbital-optimization, which consists in imation assumes that molecular structure and dynamics modifying the truncated Hilbert space (without modify- can be treated in two parts: considering first the atomic ing its size) to hopefully recover the conical intersection. nuclei, and then the motion of the electrons around them. On classical computers, this orbital-optimization is part Atomic nuclei, being many times heavier than their elec- of the so-called state-averaged complete active space self- tronic counterparts, are often treated as classical point consistent field (SA-CASSCF) method [71, 72], known for particles with definite position RA and momentum PA its democratic treatment of multiple (possibly degener- (with A the index of the nucleus). As classical parti- ate) eigenstates. cles have no zero-point energy, the ground state of the For quantum computers in the NISQ era, the concept nuclei is motionless, i.e. PA = 0. By contrast, the of active space also represents an important issue. In- light electrons behave as fundamentally quantum parti- deed, the limitations of quantum resources (e.g. max- cles, and form a collective many-body wavefunction Ψ(r), imum number of qubits) constrain algorithms to focus where r = (r ; r ;:::; r ) is the 3N -dimensional global on small-sized active spaces, thus preventing any calcu- 1 2 Ne e position vector of the Ne electrons. This wavefunction lation on large systems. Designing Hilbert spaces so that evolves by the Schr¨odingerequation with the fixed nuclei the improved accuracy from a quantum computer is not providing a background potential field. The quantum op- rendered useless by the active space approximation itself erator describing this dynamics is the so-called electronic is key to finding useful applications of these promising structure Hamiltonian (in atomic units) devices within the field of chemistry. In this context, getting access to complicated chemical features such as Ne Ne Ne Na ^ 1 X 1 X 1 X X ZA avoided crossing or conical intersections can appear as a H = − rri + − ; (1) 2 2 rij riA real challenge for near-term quantum computers. i=1 i6=j i=1 A=1 Motivated by this issue, we introduce in this paper a NISQ-friendly variational hybrid quantum/classical al- where rij = jri − rjj, riA = jri − RAj and ZA are re- gorithm inspired by the SA-CASSCF method, so-called spectively the distance between electrons i and j, the state-averaged orbital-optimized VQE (SA-OO-VQE). distance between electron i and nuclei A, and the atomic The ability of SA-OO-VQE to capture a conical inter- number of atom A. In this context, solving the electronic section is demonstrated on a minimal Schiff base model structure problem means solving the time-independent of the rhodopsin chromophore, i.e. the formaldimine Schr¨odingerequation molecule. ^ The paper is organized as follows. After a brief descrip- H jΨ`i = E` jΨ`i (2) tion of the electronic structure problem in Sec. IIA, the concept of conical intersection is introduced in Sec.IIB. which yields a set of electronic eigenstates jΨ`i ≡ Ψ`(r) A particular attention is paid to the efficiency of classical with corresponding energies E`. These eigenstates and algorithms to represent this characteristic in chemistry energies depend on the atomic positions themselves | (Sec. IIC), as shown on the formaldimine molecule in E` = E`(R) and jΨ`i = jΨ`(R)i, where R is the 3Na- Sec. IID. Then, the SA-OO-VQE hybrid quantum/clas- dimensional position vector of the Na nuclei. (Note the sical algorithm is detailed in Sec. III, where each building change of notation: Ψ`(r) is a single number from evalu- blocks of the algorithm is described separately. Following ating the function Ψ` at the fixed electronic co-ordinates a brief introduction on quantum variational algorithms in r, while jΨ`(R)i is an entire function defined by a fixed Sec.IIIA, we detail how SA-OO-VQE mimics its classical set of nuclear co-ordinates R). In practice, it is conve- analog SA-CASSCF by combining a state-averaged VQE nient to work in a finite basis of No orthonormal molec- (Sec. IIIB1) with a state-averaged orbital-optimization ular orbitals (MO) fφp(r)g. The latter are usually so- procedure (Sec. IIIB2). The importance of choosing a lutions to the mean-field single-particle problem, deter- flexible wavefunction ansatz such as the generalized uni- mined by the Hartree{Fock (HF) method [72].
Recommended publications
  • BCG) Is a Global Management Consulting Firm and the World’S Leading Advisor on Business Strategy
    The Next Decade in Quantum Computing— and How to Play Boston Consulting Group (BCG) is a global management consulting firm and the world’s leading advisor on business strategy. We partner with clients from the private, public, and not-for-profit sectors in all regions to identify their highest-value opportunities, address their most critical challenges, and transform their enterprises. Our customized approach combines deep insight into the dynamics of companies and markets with close collaboration at all levels of the client organization. This ensures that our clients achieve sustainable competitive advantage, build more capable organizations, and secure lasting results. Founded in 1963, BCG is a private company with offices in more than 90 cities in 50 countries. For more information, please visit bcg.com. THE NEXT DECADE IN QUANTUM COMPUTING— AND HOW TO PLAY PHILIPP GERBERT FRANK RUESS November 2018 | Boston Consulting Group CONTENTS 3 INTRODUCTION 4 HOW QUANTUM COMPUTERS ARE DIFFERENT, AND WHY IT MATTERS 6 THE EMERGING QUANTUM COMPUTING ECOSYSTEM Tech Companies Applications and Users 10 INVESTMENTS, PUBLICATIONS, AND INTELLECTUAL PROPERTY 13 A BRIEF TOUR OF QUANTUM COMPUTING TECHNOLOGIES Criteria for Assessment Current Technologies Other Promising Technologies Odd Man Out 18 SIMPLIFYING THE QUANTUM ALGORITHM ZOO 21 HOW TO PLAY THE NEXT FIVE YEARS AND BEYOND Determining Timing and Engagement The Current State of Play 24 A POTENTIAL QUANTUM WINTER, AND THE OPPORTUNITY THEREIN 25 FOR FURTHER READING 26 NOTE TO THE READER 2 | The Next Decade in Quantum Computing—and How to Play INTRODUCTION he experts are convinced that in time they can build a Thigh-performance quantum computer.
    [Show full text]
  • Chapter 3 Quantum Circuit Model of Computation
    Chapter 3 Quantum Circuit Model of Computation 3.1 Introduction In the previous chapter we saw that any Boolean function can be computed from a reversible circuit. In particular, in principle at least, this can be done from a small set of universal logical gates which do not dissipate heat dissi- pation (so we have logical reversibility as well as physical reversibility. We also saw that quantum evolution of the states of a system (e.g. a collection of quantum bits) is unitary (ignoring the reading or measurement process of the bits). A unitary matrix is of course invertible, but also conserves entropy (at the present level you can think of this as a conservation of scalar prod- uct which implies conservation of probabilities and entropies). So quantum evolution is also both \logically" and physically reversible. The next natural question is then to ask if we can use quantum operations to perform computational tasks, such as computing a Boolean function? Do the principles of quantum mechanics bring in new limitations with respect to classical computations? Or do they on the contrary bring in new ressources and advantages? These issues were raised and discussed by Feynman, Benioff and Manin in the early 1980's. In principle QM does not bring any new limitations, but on the contrary the superposition principle applied to many particle systems (many qubits) enables us to perform parallel computations, thereby speed- ing up classical computations. This was recognized very early by Feynman who argued that classical computers cannot simulate efficiently quantum me- chanical processes. The basic reason is that general quantum states involve 1 2 CHAPTER 3.
    [Show full text]
  • Luigi Frunzio
    THESE DE DOCTORAT DE L’UNIVERSITE PARIS-SUD XI spécialité: Physique des solides présentée par: Luigi Frunzio pour obtenir le grade de DOCTEUR de l’UNIVERSITE PARIS-SUD XI Sujet de la thèse: Conception et fabrication de circuits supraconducteurs pour l'amplification et le traitement de signaux quantiques soutenue le 18 mai 2006, devant le jury composé de MM.: Michel Devoret Daniel Esteve, President Marc Gabay Robert Schoelkopf Rapporteurs scientifiques MM.: Antonio Barone Carlo Cosmelli ii Table of content List of Figures vii List of Symbols ix Acknowledgements xvii 1. Outline of this work 1 2. Motivation: two breakthroughs of quantum physics 5 2.1. Quantum computation is possible 5 2.1.1. Classical information 6 2.1.2. Quantum information unit: the qubit 7 2.1.3. Two new properties of qubits 8 2.1.4. Unique property of quantum information 10 2.1.5. Quantum algorithms 10 2.1.6. Quantum gates 11 2.1.7. Basic requirements for a quantum computer 12 2.1.8. Qubit decoherence 15 2.1.9. Quantum error correction codes 16 2.2. Macroscopic quantum mechanics: a quantum theory of electrical circuits 17 2.2.1. A natural test bed: superconducting electronics 18 2.2.2. Operational criteria for quantum circuits 19 iii 2.2.3. Quantum harmonic LC oscillator 19 2.2.4. Limits of circuits with lumped elements: need for transmission line resonators 22 2.2.5. Hamiltonian of a classical electrical circuit 24 2.2.6. Quantum description of an electric circuit 27 2.2.7. Caldeira-Leggett model for dissipative elements 27 2.2.8.
    [Show full text]
  • Evaluation of Gaussian Molecular Integrals I
    The Mathematica® Journal Evaluation of Gaussian Molecular Integrals I. Overlap Integrals Minhhuy Hô Julio Manuel Hernández-Pérez This article discusses the evaluation of molecular overlap integrals for Gaussian-type functions with arbitrary angular dependence. As an example, we calculate the overlap matrix for the water molecule in the STO-3G basis set. ‡ Introduction Computational quantum chemistry makes extensive use of various integrals (and their derivatives) of the general form [1, 2, 3] ¶ ‡ caHrL O cbHrL dr, (1) -¶ where caHrL is an unnormalized Cartesian Gaussian function centered at A = 9Ax, Ay, Az=: 2 ax ay az -a r-A caHr; a, A, aL = Hx - AxL Iy - AyM Hz - AzL e , (2) where A is normally taken at the nucleus, a is the orbital exponent, and the polynomial rep- resents the angular part, in that the sum of the Cartesian angular momenta ax + ay + az = 0, 1, 2, … corresponds to functions of type s, p, d, f, …. When the operator O is 1, one simply has the overlap/density integral; otherwise it can be the energy operator for kinetic energy - 1 2, electron-nuclear attraction r - R -1, or electron-electron 2 ! -1 repulsion ri - rj (which would involve double integrals). Other molecular properties involving external fields (response functions) or transition moments can also be computed from integrals of this form. The Mathematica Journal 14 © 2012 Wolfram Media, Inc. 2 Minhhuy Hô and Julio Manuel Hernández-Pérez · Gaussian-Type Functions Gaussian-type functions are not the most natural choice for expanding the wavefunction. Slater-type functions, where the exponent is -a r - A instead, can describe atomic systems more realistically.
    [Show full text]
  • Multidimensional Dynamics Involving a Conical Intersection: Wavepacket Calculations Using the MCTDH Method
    Multidimensional dynamics involving a conical intersection: Wavepacket calculations using the MCTDH method G. A. Worth†, H.-D. Meyer∗, and L. S. Cederbaum∗ ∗ Theoretical Chemistry, University of Heidelberg, INF 229, 69120 Heidelberg, Germany † Dept. of Chemistry, King’s College London, The Strand, London, WC2R 2LS, U.K. This article will appear in Conical Intersections, W. Domcke, D. R. Yarkony, and H. K¨oppel, Eds., World Scientific Co., (Singapore, 2004), pp. 573–607. June 1, 2004 1 Introduction The time-evolution of a molecular system through a conical intersection is multi- dimensional in nature, i.e. its description cannot be reduced to a single reaction coordinate. In addition to this, such processes are inherently quantum mechan- ical due to the mixing of electronic states as the intersection is traversed, and a powerful quantum dynamics method is required to treat these, often large, systems accurately. The most intuitive way of obtaining a numerically exact solution of the time- dependent Schr¨odinger equation is to use standard wavepacket propagation meth- ods. These represent the evolving system wavefunction using basis set expansion techniques or grid representations [1–3]. Unfortunately, the standard methods are restricted by an exponential increase in the computational resources required with the number of degrees of freedom treated, usually no more than 4-6 nuclear modes can be included. Systems in which vibronic coupling plays a role are often much larger than this. To treat these systems using the standard methods it is then necessary to turn to reduced dimensionality models in which only a few degrees of freedom are explicitly included, and the remainder either ignored or included in a phenomenological manner.
    [Show full text]
  • Optimization of Reversible Circuits Using Toffoli Decompositions with Negative Controls
    S S symmetry Article Optimization of Reversible Circuits Using Toffoli Decompositions with Negative Controls Mariam Gado 1,2,* and Ahmed Younes 1,2,3 1 Department of Mathematics and Computer Science, Faculty of Science, Alexandria University, Alexandria 21568, Egypt; [email protected] 2 Academy of Scientific Research and Technology(ASRT), Cairo 11516, Egypt 3 School of Computer Science, University of Birmingham, Birmingham B15 2TT, UK * Correspondence: [email protected]; Tel.: +203-39-21595; Fax: +203-39-11794 Abstract: The synthesis and optimization of quantum circuits are essential for the construction of quantum computers. This paper proposes two methods to reduce the quantum cost of 3-bit reversible circuits. The first method utilizes basic building blocks of gate pairs using different Toffoli decompositions. These gate pairs are used to reconstruct the quantum circuits where further optimization rules will be applied to synthesize the optimized circuit. The second method suggests using a new universal library, which provides better quantum cost when compared with previous work in both cost015 and cost115 metrics; this proposed new universal library “Negative NCT” uses gates that operate on the target qubit only when the control qubit’s state is zero. A combination of the proposed basic building blocks of pairs of gates and the proposed Negative NCT library is used in this work for synthesis and optimization, where the Negative NCT library showed better quantum cost after optimization compared with the NCT library despite having the same circuit size. The reversible circuits over three bits form a permutation group of size 40,320 (23!), which Citation: Gado, M.; Younes, A.
    [Show full text]
  • Multi-Mode Ultra-Strong Coupling in Circuit Quantum Electrodynamics
    www.nature.com/npjqi ARTICLE OPEN Multi-mode ultra-strong coupling in circuit quantum electrodynamics Sal J. Bosman1, Mario F. Gely1, Vibhor Singh2, Alessandro Bruno3, Daniel Bothner1 and Gary A. Steele1 With the introduction of superconducting circuits into the field of quantum optics, many experimental demonstrations of the quantum physics of an artificial atom coupled to a single-mode light field have been realized. Engineering such quantum systems offers the opportunity to explore extreme regimes of light-matter interaction that are inaccessible with natural systems. For instance the coupling strength g can be increased until it is comparable with the atomic or mode frequency ωa,m and the atom can be coupled to multiple modes which has always challenged our understanding of light-matter interaction. Here, we experimentally realize a transmon qubit in the ultra-strong coupling regime, reaching coupling ratios of g/ωm = 0.19 and we measure multi-mode interactions through a hybridization of the qubit up to the fifth mode of the resonator. This is enabled by a qubit with 88% of its capacitance formed by a vacuum-gap capacitance with the center conductor of a coplanar waveguide resonator. In addition to potential applications in quantum information technologies due to its small size, this architecture offers the potential to further explore the regime of multi-mode ultra-strong coupling. npj Quantum Information (2017) 3:46 ; doi:10.1038/s41534-017-0046-y INTRODUCTION decreasing gate times20 as well as the performance of quantum 21 Superconducting circuits such as microwave cavities and Joseph- memories. With very strong coupling rates, the additional modes son junction based artificial atoms1 have opened up a wealth of of an electromagnetic resonator become increasingly relevant, new experimental possibilities by enabling much stronger light- and U/DSC can only be understood in these systems if the multi- matter coupling than in analog experiments with natural atoms2 mode effects are correctly modeled.
    [Show full text]
  • Concentric Transmon Qubit Featuring Fast Tunability and an Anisotropic Magnetic Dipole Moment
    Concentric transmon qubit featuring fast tunability and an anisotropic magnetic dipole moment Cite as: Appl. Phys. Lett. 108, 032601 (2016); https://doi.org/10.1063/1.4940230 Submitted: 13 October 2015 . Accepted: 07 January 2016 . Published Online: 21 January 2016 Jochen Braumüller, Martin Sandberg, Michael R. Vissers, Andre Schneider, Steffen Schlör, Lukas Grünhaupt, Hannes Rotzinger, Michael Marthaler, Alexander Lukashenko, Amadeus Dieter, Alexey V. Ustinov, Martin Weides, and David P. Pappas ARTICLES YOU MAY BE INTERESTED IN A quantum engineer's guide to superconducting qubits Applied Physics Reviews 6, 021318 (2019); https://doi.org/10.1063/1.5089550 Planar superconducting resonators with internal quality factors above one million Applied Physics Letters 100, 113510 (2012); https://doi.org/10.1063/1.3693409 An argon ion beam milling process for native AlOx layers enabling coherent superconducting contacts Applied Physics Letters 111, 072601 (2017); https://doi.org/10.1063/1.4990491 Appl. Phys. Lett. 108, 032601 (2016); https://doi.org/10.1063/1.4940230 108, 032601 © 2016 AIP Publishing LLC. APPLIED PHYSICS LETTERS 108, 032601 (2016) Concentric transmon qubit featuring fast tunability and an anisotropic magnetic dipole moment Jochen Braumuller,€ 1,a) Martin Sandberg,2 Michael R. Vissers,2 Andre Schneider,1 Steffen Schlor,€ 1 Lukas Grunhaupt,€ 1 Hannes Rotzinger,1 Michael Marthaler,3 Alexander Lukashenko,1 Amadeus Dieter,1 Alexey V. Ustinov,1,4 Martin Weides,1,5 and David P. Pappas2 1Physikalisches Institut, Karlsruhe Institute of Technology,
    [Show full text]
  • HETEROLEPTIC PADDLEWHEEL COMPLEXES and MOLECULAR ASSEMBLIES of DIMOLYBDENUM and DITUNGSTEN: a Study of Electronic and Structural Control
    HETEROLEPTIC PADDLEWHEEL COMPLEXES AND MOLECULAR ASSEMBLIES OF DIMOLYBDENUM AND DITUNGSTEN: A Study of Electronic and Structural Control DISSERTATION Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the Graduate School of The Ohio State University By DOUGLAS J. BROWN, B.S., M.S. ***** The Ohio State University 2006 Dissertation Committee: Approved by Professor MALCOLM H. CHISHOLM, Advisor Professor CLAUDIA TURRO _________________________________ Professor CHRISTOPHER M. HADAD Advisor Graduate Program in Chemistry ABSTRACT Molybdenum and tungsten have a rich chemistry in quadruply-bonded paddlewheel complexes that are very similar in many ways, yet significantly different in other ways. Details and findings from Density Functional Theory studies of homoleptic and heteroleptic ligand complexes as well as simple pairs of dimers linked by a bifunctional oxalate bridge are discussed. Special interest is given to metal-to-ligand interactions in paddlewheel complexes and metal-to-bridge interactions in simple pairs of dimers. Computational studies provide great insight, reliable explanations, and qualitative predictions of physicochemical properties. A new series of dimolybdenum and ditungsten paddlewheel complexes with i i mixed amidinate-carboxylate ligands of the form trans-M2(O2CCH3)2( PrNC(R)N Pr)2 where M = Mo, W and R = Me, C≡CtBu, C≡CPh, or C≡CFc have been synthesized and characterized with the exception of the tungsten complex where R = Me. The synthesis of a series of complexes with each metal allowed for the analysis and comparison of molecular structural features, spectroscopic properties, and electrochemical behavior. The dinuclear centers exhibit redox chemistry which is dependent upon the nature of the ancillary ligand, and the electrochemical properties can be tuned through ligand ii modifications.
    [Show full text]
  • Orbital Exchange Calculations: an Alternative View of the Chemical Bond Paul B
    Orbital Exchange Calculations: An Alternative View of the Chemical Bond Paul B. Merrithewa PO Box 120, Amherst, New Hampshire 03031-0120 USA Abstract The purpose of this paper is to explore a model of the chemical bond which does not assume that the electrons of the chemical bonding electron pair can be unambiguously identified with either the left hand or right hand of the bonding atoms when their orbitals overlap to bond. In order to provide maximum flexibility in the selection of the electron’s orbitals, the orbitals have been represented as spatial arrays and the calculations performed numerically. This model of the chemical bond assumes that the identifiability of the bonding electrons is a function of 1-(overlap/(1+overlap)) where the overlap of the two bonding electron’s orbitals is calculated in the usual manner. The kinetic energy of the bonding electron pair and the energy required to meet the orthogonality requirements, mandated by the Pauli principle, are a function of overlap/(1+overlap). The model assumes that the bonding orbitals are straight-forward atomic orbitals or hybrids of these atomic orbitals. The results obtained by applying this simple approach to eleven di-atomics and seven common poly-atomics are quite good. The calculated bond lengths are generally within 0.005Å of the measured values and bond energies to within a few percent. Bond lengths for bonds to H are about 0.02 Å high. Except for H2, bond lengths are determined, independent of bond energy, at that point where overlap/(1+overlap) equals 0.5. I. INTRODUCTION The purpose of the work described here is to test a model for the chemical bond that assumes that the two electrons of the bonding pair are not completely identifiable, with respect to their source atom, when their orbitals overlap.
    [Show full text]
  • Chapter 1: Electronic Structure and Bonding a Brief Review of General Chemistry
    Chapter 1: Electronic Structure and Bonding Or A Brief Review of General Chemistry Part 2: Quantum Mechanics and Molecular Orbital Theory Review ideas from general chemistry: atoms, bonds, molecular geometry Hybrid Orbitals of Ethene, C2H4 Bonding in Ethene (Ethylene): A Double Bond The sigma bond is unaffected by rotation of one of the CH2 groups. The overlap of the p orbitals is disrupted by rotation of one of the CH2 groups. This would cause the double bond to break. Hybrid Orbitals of Ethene, C2H4 An sp2-Hybridized Carbon - The bond angle in the sp2 carbon is 120°. The sp2 carbon is the trigonal planar carbon. Molecular Orbital Theory So let’s look at Ethyne (Acetylene) . Draw a Lewis Dot Structure for Ethyne Hybrid Orbitals of Ethyne, C2H2 Bonding in Ethyne (Acetylene) : A Triple Bond A triple bond consists of one s bond and two p bonds. Bond angle of the sp carbon: 180°. Molecular Orbital Theory + - + Molecular Orbital Theory Summary The sharing of one pair of electrons is a single bond. The sharing of two pairs gives a double bond. The sharing of three pairs gives a triple bond. The greater the electron density in the region of orbital overlap, the stronger the bond. The more s character, the shorter and stronger is the bond. The more s character, the larger is the bond angle. A bond is weaker than a bond. Hybridization of Nitrogen and Oxygen Elements other than C can have hybridized orbitals H–N–H bond angle in ammonia (NH3) 107.3° C-N-H bond angle is 110.3 ° N’s orbitals (sppp) hybridize to form four sp3 orbitals One sp3 orbital is occupied by two nonbonding electrons, and three sp3 orbitals have one electron each, forming bonds to H and CH3.
    [Show full text]
  • European Conference on the Dynamics of Molecular Systems Book of Abstracts
    European conference on the dynamics of molecular systems Book of abstracts Department of Chemistry Jesus College University of Oxford 1 2 Plenary talks 3 PLENARY TALK P1 Imaging studies of inelastic collisions David H. Parker Dept of Molecular and Laser Physics, IMM, Radboud University, Nijmegen, the Netherlands. E‐mail address: [email protected] State-to-state imaging studies of the interaction of small molecules including O2,H2O, OH, CO and ammonia with H2 and with He will be described. These molecules are prominent components of molecular interstellar matter. In order to extract quantitative information on these molecules from telescope spectra, models are employed that depend critically on the rates of rotational energy exchange due to collisions with molecular hydrogen and helium. Collision rates are currently determined by theory from the multidimensional Potential Energy Surface (PES) describing the interaction of H2O and H2 [1] or He [2] and OH with H2 [3], or He. Our velocity map imaging [4] measurements of state-to-state differential and relative integral cross sections of rotational inelastic collisions, also as a function of collision energy, are used to test these PESs. For the H2O/H2, He system experiment is compared with state-of-the-art theoretical calculations by the group of L. Wiesenfeld (Grenoble) [5]. Our studies of astrochemistry relevant small molecules with the collision partners He and H2 at collision energy relevant to that of the interstellar media should place the theoretically determined PESs and the collision rates extracted from the PES on a firmer basis. References 1. P. Valiron M. Wernli, A.
    [Show full text]