DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Calculation of Molecular Vibrational Spectra on a Quantum Annealer

Abstract

Until recently molecular energy calculations using quantum computing hardware have been limited to gate-based quantum computers. In this paper, a new methodology is presented to calculate the vibrational spectrum of a molecule on a quantum annealer. The key idea of the method is a mapping of the ground state variational problem onto an Ising or quadratic unconstrained binary optimization (QUBO) problem by expressing the expansion coefficients using spins or qubits. The algorithm is general and represents a new revolutionary approach for solving the real symmetric eigenvalue problem on a quantum annealer. The method is applied to two chemically important molecules: O2 (oxygen) and O3 (ozone). The lowest two vibrational states of these molecules are computed using both a hardware quantum annealer and a software based classical annealer. Extension of the algorithm to higher dimensions is explicitly demonstrated for an N-dimensional harmonic oscillator (N ≤ 5). The algorithm scales exponentially with dimensionality if a direct product basis is used but will exhibit polynomial scaling for a nondirect product basis.

Authors:
ORCiD logo [1];  [1]; ORCiD logo [2]
  1. Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
  2. Department of Chemistry, Marquette University, Milwaukee, Wisconsin 53021, United States
Publication Date:
Research Org.:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1545969
Alternate Identifier(s):
OSTI ID: 1558208
Report Number(s):
LA-UR-18-29586
Journal ID: ISSN 1549-9618
Grant/Contract Number:  
20170221ER; 89233218CNA000001
Resource Type:
Published Article
Journal Name:
Journal of Chemical Theory and Computation
Additional Journal Information:
Journal Name: Journal of Chemical Theory and Computation Journal Volume: 15 Journal Issue: 8; Journal ID: ISSN 1549-9618
Publisher:
American Chemical Society
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY; Inorganic and Physical Chemistry; quantum computing; quantum annealing; eigensolver; molecular spectrum; QAE

Citation Formats

Teplukhin, Alexander, Kendrick, Brian K., and Babikov, Dmitri. Calculation of Molecular Vibrational Spectra on a Quantum Annealer. United States: N. p., 2019. Web. doi:10.1021/acs.jctc.9b00402.
Teplukhin, Alexander, Kendrick, Brian K., & Babikov, Dmitri. Calculation of Molecular Vibrational Spectra on a Quantum Annealer. United States. https://doi.org/10.1021/acs.jctc.9b00402
Teplukhin, Alexander, Kendrick, Brian K., and Babikov, Dmitri. Wed . "Calculation of Molecular Vibrational Spectra on a Quantum Annealer". United States. https://doi.org/10.1021/acs.jctc.9b00402.
@article{osti_1545969,
title = {Calculation of Molecular Vibrational Spectra on a Quantum Annealer},
author = {Teplukhin, Alexander and Kendrick, Brian K. and Babikov, Dmitri},
abstractNote = {Until recently molecular energy calculations using quantum computing hardware have been limited to gate-based quantum computers. In this paper, a new methodology is presented to calculate the vibrational spectrum of a molecule on a quantum annealer. The key idea of the method is a mapping of the ground state variational problem onto an Ising or quadratic unconstrained binary optimization (QUBO) problem by expressing the expansion coefficients using spins or qubits. The algorithm is general and represents a new revolutionary approach for solving the real symmetric eigenvalue problem on a quantum annealer. The method is applied to two chemically important molecules: O2 (oxygen) and O3 (ozone). The lowest two vibrational states of these molecules are computed using both a hardware quantum annealer and a software based classical annealer. Extension of the algorithm to higher dimensions is explicitly demonstrated for an N-dimensional harmonic oscillator (N ≤ 5). The algorithm scales exponentially with dimensionality if a direct product basis is used but will exhibit polynomial scaling for a nondirect product basis.},
doi = {10.1021/acs.jctc.9b00402},
journal = {Journal of Chemical Theory and Computation},
number = 8,
volume = 15,
place = {United States},
year = {Wed Jul 17 00:00:00 EDT 2019},
month = {Wed Jul 17 00:00:00 EDT 2019}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record
https://doi.org/10.1021/acs.jctc.9b00402

Citation Metrics:
Cited by: 20 works
Citation information provided by
Web of Science

Figures / Tables:

Figure 1 Figure 1: Computed wave functions of the ground and first excited states. The molecules are (A) O2 and (B) O3, and the wave function is squared for ozone. Red and green bold curves are the results of a classical numerical eigensolver (LAPACK). Solid thin black curves were obtained using QAEmore » with a software classical QUBO solver (they lie on top of the LAPACK curves). Dashed black curves are the results obtained using QAE with a hardware quantum annealer (DWave machine). The blue curve is (A) oxygen potential V(r) and (B) the minimum energy path V(ρ) for ozone.« less

Save / Share:

Works referenced in this record:

Simulated Quantum Computation of Molecular Energies
journal, September 2005


Quantum computers
journal, March 2010


LAPACK Users' Guide
software, January 1999


Structure of small molecular hydrogen clusters
journal, March 1991

  • Rama Krishna, M. V.; Whaley, K. B.
  • Zeitschrift f�r Physik D Atoms, Molecules and Clusters, Vol. 20, Issue 1-4
  • DOI: 10.1007/BF01543978

Communication: An accurate global potential energy surface for the ground electronic state of ozone
journal, November 2013

  • Dawes, Richard; Lolur, Phalgun; Li, Anyang
  • The Journal of Chemical Physics, Vol. 139, Issue 20
  • DOI: 10.1063/1.4837175

Accurate localized and delocalized vibrational states of HCN/HNC
journal, March 1987

  • Bačić, Z.; Light, J. C.
  • The Journal of Chemical Physics, Vol. 86, Issue 6
  • DOI: 10.1063/1.452017

Adiabatic quantum computation
journal, January 2018


Quantum annealing with manufactured spins
journal, May 2011

  • Johnson, M. W.; Amin, M. H. S.; Gildert, S.
  • Nature, Vol. 473, Issue 7346
  • DOI: 10.1038/nature10012

The ratio of the number of states in asymmetric and symmetric ozone molecules deviates from the statistical value of 2
journal, March 2019

  • Gayday, Igor; Teplukhin, Alexander; Babikov, Dmitri
  • The Journal of Chemical Physics, Vol. 150, Issue 10
  • DOI: 10.1063/1.5082850

Quantum Monte Carlo simulations of selected ammonia clusters ( n  = 2–5): Isotope effects on the ground state of typical hydrogen bonded systems
journal, December 2010

  • Curotto, E.; Mella, Massimo
  • The Journal of Chemical Physics, Vol. 133, Issue 21
  • DOI: 10.1063/1.3506027

Efficient method for calculations of ro-vibrational states in triatomic molecules near dissociation threshold: Application to ozone
journal, September 2016

  • Teplukhin, Alexander; Babikov, Dmitri
  • The Journal of Chemical Physics, Vol. 145, Issue 11
  • DOI: 10.1063/1.4962914

Signatures of Large-Amplitude Vibrations in the Spectra of H 5 + and D 5 +
journal, December 2012

  • Lin, Zhou; McCoy, Anne B.
  • The Journal of Physical Chemistry Letters, Vol. 3, Issue 24
  • DOI: 10.1021/jz3017683

Isotope effect on the melting of para-hydrogen and ortho-deuterium clusters
journal, September 1992


The theory of variational hybrid quantum-classical algorithms
journal, February 2016


Diffusion Monte Carlo Calculations of Zero-Point Structures of Partially Deuterated Isotopologues of H 7 +
journal, April 2014

  • Qu, Chen; Bowman, Joel M.
  • The Journal of Physical Chemistry B, Vol. 118, Issue 28
  • DOI: 10.1021/jp501371z

A variational eigenvalue solver on a photonic quantum processor
journal, July 2014

  • Peruzzo, Alberto; McClean, Jarrod; Shadbolt, Peter
  • Nature Communications, Vol. 5, Issue 1
  • DOI: 10.1038/ncomms5213

MULTIMODE calculations of the infrared spectra of H 7 + and D 7 + using ab initio potential energy and dipole moment surfaces
journal, October 2013

  • Qu, Chen; Prosmiti, Rita; Bowman, Joel M.
  • Theoretical Chemistry Accounts, Vol. 132, Issue 12
  • DOI: 10.1007/s00214-013-1413-1

Unraveling Anharmonic Effects in the Vibrational Predissociation Spectra of H 5 O 2 + and Its Deuterated Analogues
journal, June 2011

  • Guasco, Timothy L.; Johnson, Mark A.; McCoy, Anne B.
  • The Journal of Physical Chemistry A, Vol. 115, Issue 23
  • DOI: 10.1021/jp109999b

Quantum Simulators
journal, October 2009


Electronic Structure Calculations and the Ising Hamiltonian
journal, October 2017

  • Xia, Rongxin; Bian, Teng; Kais, Sabre
  • The Journal of Physical Chemistry B, Vol. 122, Issue 13
  • DOI: 10.1021/acs.jpcb.7b10371

Diffusion Monte Carlo Approaches for Evaluating Rotationally Excited States of Symmetric Top Molecules: Application to H 3 O + and D 3 O +
journal, August 2009

  • Petit, Andrew S.; McCoy, Anne B.
  • The Journal of Physical Chemistry A, Vol. 113, Issue 45
  • DOI: 10.1021/jp905098k

Thermal effects on small para-hydrogen clusters
journal, November 2010

  • Navarro, Jesús; Guardiola, Rafael
  • International Journal of Quantum Chemistry, Vol. 111, Issue 2
  • DOI: 10.1002/qua.22575

Ab Initio Diffusion Monte Carlo Calculations of the Quantum Behavior of CH 5 + in Full Dimensionality
journal, May 2004

  • McCoy, Anne B.; Braams, Bastiaan J.; Brown, Alex
  • The Journal of Physical Chemistry A, Vol. 108, Issue 23
  • DOI: 10.1021/jp0487096

Deuteration Effects on the Structure and Infrared Spectrum of CH 5 +
journal, March 2006

  • Huang, Xinchuan; Johnson, Lindsay M.; Bowman, Joel M.
  • Journal of the American Chemical Society, Vol. 128, Issue 11
  • DOI: 10.1021/ja057514o

Investigation of the Structure and Spectroscopy of H 5 + Using Diffusion Monte Carlo
journal, May 2013

  • Lin, Zhou; McCoy, Anne B.
  • The Journal of Physical Chemistry A, Vol. 117, Issue 46
  • DOI: 10.1021/jp4014652

The structure of para-hydrogen clusters
journal, December 2009


Structural transitions in the 309-atom magic number Lennard-Jones cluster
journal, March 2006

  • Noya, Eva G.; Doye, Jonathan P. K.
  • The Journal of Chemical Physics, Vol. 124, Issue 10
  • DOI: 10.1063/1.2173260

On the solid- and liquidlike nature of quantum clusters in their ground state
journal, June 2008

  • Cuervo, Javier Eduardo; Roy, Pierre-Nicholas
  • The Journal of Chemical Physics, Vol. 128, Issue 22
  • DOI: 10.1063/1.2938369

Universal quantum simulation with prethreshold superconducting qubits: Single-excitation subspace method
journal, June 2015

  • Geller, Michael R.; Martinis, John M.; Sornborger, Andrew T.
  • Physical Review A, Vol. 91, Issue 6
  • DOI: 10.1103/PhysRevA.91.062309

Visualization of Potential Energy Function Using an Isoenergy Approach and 3D Prototyping
journal, November 2014

  • Teplukhin, Alexander; Babikov, Dmitri
  • Journal of Chemical Education, Vol. 92, Issue 2
  • DOI: 10.1021/ed500683g

Scalable Quantum Simulation of Molecular Energies
journal, July 2016

  • O’Malley, P. J. J.; Babbush, R.; Kivlichan, I. D.
  • Physical Review X, Vol. 6, Issue 3
  • DOI: 10.1103/PhysRevX.6.031007

Simple Proof of Equivalence between Adiabatic Quantum Computation and the Circuit Model
journal, August 2007


Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.