Discontinuous Galerkin discretization for quantum simulation of chemistry

Methods for electronic structure based on Gaussian and molecular orbital discretizations offer a well established, compact representation that forms much of the foundation of correlated quantum chemistry calculations on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from $O(N^{4.5})$ in molecular orbital representations to $O(N^{2.6})$ in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1-2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.

[1]  D. Hamann Optimized norm-conserving Vanderbilt pseudopotentials , 2013, 1306.4707.

[2]  George Em Karniadakis,et al.  The Development of Discontinuous Galerkin Methods , 2000 .

[3]  Chao Yang,et al.  Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework II: Force, vibration, and molecular dynamics calculations , 2012, J. Comput. Phys..

[4]  Jeremy L O'Brien,et al.  Quantum simulation of Hamiltonian spectra on a silicon chip , 2016 .

[5]  Scott B. Baden,et al.  Parallel implementation of γ‐point pseudopotential plane‐wave DFT with exact exchange , 2011, J. Comput. Chem..

[6]  J. Carter,et al.  Hybrid Quantum-Classical Hierarchy for Mitigation of Decoherence and Determination of Excited States , 2016, 1603.05681.

[7]  Alexei Y. Kitaev,et al.  Quantum measurements and the Abelian Stabilizer Problem , 1995, Electron. Colloquium Comput. Complex..

[8]  Alán Aspuru-Guzik,et al.  Exploiting Locality in Quantum Computation for Quantum Chemistry. , 2014, The journal of physical chemistry letters.

[9]  M. Hastings,et al.  Gate count estimates for performing quantum chemistry on small quantum computers , 2013, 1312.1695.

[10]  S. Goedecker Linear scaling electronic structure methods , 1999 .

[11]  Ryan Babbush,et al.  Exponentially more precise quantum simulation of fermions in the configuration interaction representation , 2015, 1506.01029.

[12]  Chao Yang,et al.  Chebyshev polynomial filtered subspace iteration in the discontinuous Galerkin method for large-scale electronic structure calculations. , 2016, The Journal of chemical physics.

[13]  Luke Schaeffer,et al.  Trading T-gates for dirty qubits in state preparation and unitary synthesis , 2018, 1812.00954.

[14]  D. Abrams,et al.  Simulation of Many-Body Fermi Systems on a Universal Quantum Computer , 1997, quant-ph/9703054.

[15]  Garnet Kin-Lic Chan,et al.  Efficient tree tensor network states (TTNS) for quantum chemistry: generalizations of the density matrix renormalization group algorithm. , 2013, The Journal of chemical physics.

[16]  J. Pask,et al.  Finite element methods in ab initio electronic structure calculations , 2005 .

[17]  M. Hastings,et al.  Progress towards practical quantum variational algorithms , 2015, 1507.08969.

[18]  Chao Yang,et al.  DGDFT: A massively parallel method for large scale density functional theory calculations. , 2015, The Journal of chemical physics.

[19]  I. Chuang,et al.  Hamiltonian Simulation by Qubitization , 2016, Quantum.

[20]  A. Kitaev,et al.  Fermionic Quantum Computation , 2000, quant-ph/0003137.

[21]  Chao Yang,et al.  Edge reconstruction in armchair phosphorene nanoribbons revealed by discontinuous Galerkin density functional theory. , 2015, Physical chemistry chemical physics : PCCP.

[22]  T. H. Dunning Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[23]  F. Verstraete,et al.  Tensor product methods and entanglement optimization for ab initio quantum chemistry , 2014, 1412.5829.

[24]  H. Lischka,et al.  PNO–CI (pair natural orbital configuration interaction) and CEPA–PNO (coupled electron pair approximation with pair natural orbitals) calculations of molecular systems. I. Outline of the method for closed‐shell states , 1975 .

[25]  Vikram Gavini,et al.  Large-scale all-electron density functional theory calculations using an enriched finite-element basis , 2016, 1610.02121.

[26]  White,et al.  Density matrix formulation for quantum renormalization groups. , 1992, Physical review letters.

[27]  J. Gambetta,et al.  Hardware-efficient variational quantum eigensolver for small molecules and quantum magnets , 2017, Nature.

[28]  Bryan O'Gorman,et al.  Generalized swap networks for near-term quantum computing , 2019, ArXiv.

[29]  E Weinan,et al.  Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation , 2011, J. Comput. Phys..

[30]  Frank Neese,et al.  An efficient and near linear scaling pair natural orbital based local coupled cluster method. , 2013, The Journal of chemical physics.

[31]  Aihui Zhou,et al.  Adaptive Finite Element Approximations for Kohn-Sham Models , 2013, Multiscale Model. Simul..

[32]  M. Head‐Gordon,et al.  Simulated Quantum Computation of Molecular Energies , 2005, Science.

[33]  R. Pooser,et al.  Cloud Quantum Computing of an Atomic Nucleus. , 2018, Physical Review Letters.

[34]  Jiangfeng Du,et al.  NMR implementation of a molecular hydrogen quantum simulation with adiabatic state preparation. , 2010, Physical review letters.

[35]  Baldereschi,et al.  Self-consistent Hartree-Fock and screened-exchange calculations in solids: Application to silicon. , 1986, Physical review. B, Condensed matter.

[36]  S. Lloyd,et al.  Quantum Algorithm Providing Exponential Speed Increase for Finding Eigenvalues and Eigenvectors , 1998, quant-ph/9807070.

[37]  J. McClean,et al.  Application of fermionic marginal constraints to hybrid quantum algorithms , 2018, 1801.03524.

[38]  Alán Aspuru-Guzik,et al.  The theory of variational hybrid quantum-classical algorithms , 2015, 1509.04279.

[39]  Nathan Wiebe,et al.  Hamiltonian Simulation in the Interaction Picture , 2018, 1805.00675.

[40]  David M. Ceperley,et al.  Towards the solution of the many-electron problem in real materials: equation of state of the hydrogen chain with state-of-the-art many-body methods , 2017, 1705.01608.

[41]  Matthew B. Hastings,et al.  Improving quantum algorithms for quantum chemistry , 2014, Quantum Inf. Comput..

[42]  P. Coveney,et al.  Scalable Quantum Simulation of Molecular Energies , 2015, 1512.06860.

[43]  E. Miles Stoudenmire,et al.  Multisliced gausslet basis sets for electronic structure , 2019, Physical Review B.

[44]  D. Arnold An Interior Penalty Finite Element Method with Discontinuous Elements , 1982 .

[45]  Alexandru Paler,et al.  Encoding Electronic Spectra in Quantum Circuits with Linear T Complexity , 2018, Physical Review X.

[46]  Alán Aspuru-Guzik,et al.  On the Chemical Basis of Trotter-Suzuki Errors in Quantum Chemistry Simulation , 2014, 1410.8159.

[47]  Annie Y. Wei,et al.  Exponentially more precise quantum simulation of fermions in second quantization , 2015, 1506.01020.

[48]  Alán Aspuru-Guzik,et al.  Quantum Simulation of Electronic Structure with Linear Depth and Connectivity. , 2017, Physical review letters.

[49]  T. Monz,et al.  Quantum Chemistry Calculations on a Trapped-Ion Quantum Simulator , 2018, Physical Review X.

[50]  Hartmut Neven,et al.  Quantum simulation of chemistry with sublinear scaling in basis size , 2018, npj Quantum Information.

[51]  Hartmut Neven,et al.  Low-Depth Quantum Simulation of Materials , 2017, 1706.00023.

[52]  David,et al.  Gaussian basis sets for use in correlated molecular calculations . Ill . The atoms aluminum through argon , 1999 .

[53]  Sandeep Sharma,et al.  The density matrix renormalization group in quantum chemistry. , 2011, Annual review of physical chemistry.

[54]  M. Ratner Molecular electronic-structure theory , 2000 .

[55]  M. Yung,et al.  Quantum implementation of the unitary coupled cluster for simulating molecular electronic structure , 2015, 1506.00443.

[56]  Ryan Babbush,et al.  The theory of variational hybrid quantum-classical algorithms , 2015, 1509.04279.

[57]  Frank Neese,et al.  Natural triple excitations in local coupled cluster calculations with pair natural orbitals. , 2013, The Journal of chemical physics.

[58]  Alán Aspuru-Guzik,et al.  A variational eigenvalue solver on a photonic quantum processor , 2013, Nature Communications.

[59]  Ryan Babbush,et al.  Qubitization of Arbitrary Basis Quantum Chemistry Leveraging Sparsity and Low Rank Factorization , 2019, Quantum.

[60]  Yingzhou Li,et al.  Globally Constructed Adaptive Local Basis Set for Spectral Projectors of Second Order Differential Operators , 2019, Multiscale Model. Simul..

[61]  Jarrod R. McClean,et al.  Computation of Molecular Spectra on a Quantum Processor with an Error-Resilient Algorithm , 2017, 1707.06408.

[62]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[63]  Ryan Babbush,et al.  Low rank representations for quantum simulation of electronic structure , 2018, npj Quantum Information.

[64]  Juan C. Meza,et al.  KSSOLV—a MATLAB toolbox for solving the Kohn-Sham equations , 2009, TOMS.

[65]  E. Knill,et al.  Quantum algorithms for fermionic simulations , 2000, cond-mat/0012334.

[66]  Angela K. Wilson,et al.  Gaussian basis sets for use in correlated molecular calculations. IX. The atoms gallium through krypton , 1993 .

[67]  Steven R White,et al.  Sliced Basis Density Matrix Renormalization Group for Electronic Structure. , 2017, Physical review letters.

[68]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

[69]  B. Lanyon,et al.  Towards quantum chemistry on a quantum computer. , 2009, Nature chemistry.

[70]  Phanish Suryanarayana,et al.  Two-Level Chebyshev Filter Based Complementary Subspace Method: Pushing the Envelope of Large-Scale Electronic Structure Calculations. , 2017, Journal of chemical theory and computation.

[71]  Steven R White,et al.  Hybrid grid/basis set discretizations of the Schrödinger equation. , 2017, The Journal of chemical physics.

[72]  Andrew M. Childs,et al.  Simulating Hamiltonian dynamics with a truncated Taylor series. , 2014, Physical review letters.

[73]  Phanish Suryanarayana,et al.  Discrete discontinuous basis projection method for large-scale electronic structure calculations. , 2018, The Journal of chemical physics.

[74]  S. Lehtola A review on non‐relativistic, fully numerical electronic structure calculations on atoms and diatomic molecules , 2019, International Journal of Quantum Chemistry.

[75]  Nathan Wiebe,et al.  Hamiltonian simulation using linear combinations of unitary operations , 2012, Quantum Inf. Comput..

[76]  H. Neven,et al.  Low-Depth Quantum Simulation of Materials , 2018 .

[77]  Jonathan Carter,et al.  Computation of Molecular Spectra on a Quantum Processor with an Error-Resilient Algorithm , 2018 .