In this work we demonstrate a practical prospect of using quantum annealers for simulation of molecular dynamics. A methodology developed for this goal, dubbed Quantum Differential Equations (QDE), is applied to propagate classical trajectories for the vibration of the hydrogen molecule in several regimes: nearly harmonic, highly anharmonic, and dissociative motion. The results obtained using the D-Wave 2000Q quantum annealer are all consistent and quickly converge to the analytical reference solution. Several alternative strategies for such calculations are explored and it was found that the most accurate results and the best efficiency are obtained by combining the quantum annealer with classical post-processing (greedy algorithm). Importantly, the QDE framework developed here is entirely general and can be applied to solve any system of first-order ordinary nonlinear differential equations using a quantum annealer.
more » « less- Award ID(s):
- 2102465
- PAR ID:
- 10373223
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
We report the first-ever formulation of the molecular dynamics problem as a Quadratic Unconstrained Binary Optimization (QUBO) problem and demonstrate its successful solution on a D-Wave quantum annealer. This methodology, named a Quantum Differential Equation (QDE) solver, is applied to propagate trajectories for the vibrational motion of hydrogen molecule, H2, in three different energy regimes: nearly harmonic, highly anharmonic, and above dissociation threshold. The results obtained using the D-Wave 2000Q quantum annealer are mutually consistent and quickly converge to the exact solution. Several alternative strategies for such calculations are explored and it is shown that the most accurate result and the best efficiency are obtained by combining quantum annealer with classical post-processing (known as “greedy” algorithm). Importantly, the QDE framework developed here is entirely general and can be applied to solve any system of first-order ordinary nonlinear differential equations using a quantum annealer.more » « less
-
null (Ed.)Quantum computing is a new and rapidly evolving paradigm for solving chemistry problems. In previous work, we developed the Quantum Annealer Eigensolver (QAE) and applied it to the calculation of the vibrational spectrum of a molecule on the D-Wave quantum annealer. However, the original QAE methodology was applicable to real symmetric matrices only. For many physics and chemistry problems, the diagonalization of complex matrices is required. For example, the calculation of quantum scattering resonances can be formulated as a complex eigenvalue problem where the real part of the eigenvalue is the resonance energy and the imaginary part is proportional to the resonance width. In the present work, we generalize the QAE to treat complex matrices: first complex Hermitian matrices and then complex symmetric matrices. These generalizations are then used to compute a quantum scattering resonance state in a 1D model potential for O + O collisions. These calculations are performed using both a software (classical) annealer and hardware annealer (the D-Wave 2000Q). The results of the complex QAE are also benchmarked against a standard linear algebra library (LAPACK). This work presents the first numerical solution of a complex eigenvalue problem of any kind on a quantum annealer, and it is the first treatment of a quantum scattering resonance on any quantum device.more » « less
-
One of the most fundamental problems that has no efficient solutions on classical computers is simulation of quantum systems. It has been long hypothesized that quantum computing devices are naturally more suitable for this task, but many aspects of practical implementations of such simulations remain unknown. One particularly important kind of these simulations is the simulation of molecular dynamics, i.e. prediction of time evolution for a system of interacting particles. In this work we show how a quantum annealer can be used to carry out such simulations by solving differential equations of motion, on the example of the hydrogen molecule. Although the considered system is simple, our method is well scalable and can be readily applied to more complicated systems as annealers with larger number of qubits become available. Importantly, the method is general and can be used to solve arbitrary systems of ordinary non-linear differential equations, which can be helpful not only in the field of computational chemistry, but in many other fields as well.more » « less
-
Efficiency optimization in quantum computing: balancing thermodynamics and computational performance
Abstract We investigate the computational efficiency and thermodynamic cost of the D-Wave quantum annealer under reverse-annealing with and without pausing. Our demonstration on the D-Wave 2000
Q annealer shows that the combination of reverse-annealing and pausing leads to improved computational efficiency while minimizing the thermodynamic cost compared to reverse-annealing alone. Moreover, we find that the magnetic field has a positive impact on the performance of the quantum annealer during reverse-annealing but becomes detrimental when pausing is involved. Our results, which are reproducible, provide strategies for optimizing the performance and energy consumption of quantum annealing systems employing reverse-annealing protocols. -
Quantum computers can produce a quantum encoding of the solution of a system of differential equations exponentially faster than a classical algorithm can produce an explicit description. However, while high-precision quantum algorithms for linear ordinary differential equations are well established, the best previous quantum algorithms for linear partial differential equations (PDEs) have complexity p o l y ( 1 / ϵ ) , where ϵ is the error tolerance. By developing quantum algorithms based on adaptive-order finite difference methods and spectral methods, we improve the complexity of quantum algorithms for linear PDEs to be p o l y ( d , log ( 1 / ϵ ) ) , where d is the spatial dimension. Our algorithms apply high-precision quantum linear system algorithms to systems whose condition numbers and approximation errors we bound. We develop a finite difference algorithm for the Poisson equation and a spectral algorithm for more general second-order elliptic equations.more » « less