SHARC meets TEQUILA: Mixed Quantum-Classical Dynamics on a Quantum Computer using a Hybrid Quantum-Classical Algorithm

Eduarda Sangiogo Gil,Markus Oppel,Jakob S. Kottmann,Leticia González
DOI: https://doi.org/10.26434/chemrxiv-2024-dkv5r
2024-07-29
Abstract:Recent developments in quantum computing are highly promising, particularly in the realm of quantum chemistry. Due to the noisy nature of currently available quantum hardware, hybrid quantum-classical algorithms have emerged as a reliable option for near-term simulations. Mixed quantum-classical dynamics methods effectively capture nonadiabatic effects by integrating classical nuclear dynamics with quantum chemical computations of the electronic properties. However, these methods face challenges due to the high computational cost of the quantum chemistry part. To mitigate the computational demand, we propose a method where the required electronic properties are computed through a hybrid quantum-classical approach on combination of classical and quantum hardware. This framework employs the variational quantum eigensolver and variational quantum deflation algorithms to obtain ground and excited state energies, gradients, nonadiabatic coupling vectors, and transition dipole moments. These quantities are used to propagate the nonadiabatic molecular dynamics using the Tully’s fewest switches surface hopping method, although the implementation is also compatible with other molecular dynamics approaches. The approach, implemented by integrating the molecular dynamics program package SHARC with the TEQUILA quantum computing framework, is validated by studying the cis-trans photoisomerization of methanimine and the electronic relaxation of ethylene. The results show qualitatively accurate molecular dynamics that align with experimental findings and other computational studies. This work is expected to mark a significant step towards achieving a "quantum advantage" for realistic chemical simulations.
Chemistry
What problem does this paper attempt to address?
The main goal of this paper is to propose a hybrid quantum-classical algorithm to address computational challenges in chemical dynamics simulations, particularly for non-adiabatic process dynamics simulations. Specifically, the research aims to leverage the advantages of quantum computers to accelerate the computation of electronic properties in molecular dynamics simulations. To achieve this goal, the researchers combined the molecular dynamics package SHARC with the quantum computing framework TEQUILA, developing a method that can compute key physical quantities such as ground and excited state energies, gradients, non-adiabatic coupling vectors, and transition dipole moments on quantum hardware using the Variational Quantum Eigensolver (VQE) and the Variational Quantum Deflation (VQD) algorithms. These physical quantities are then used to drive Tully's fewest switches surface hopping method to simulate non-adiabatic molecular dynamics. Key points mentioned in the paper include: 1. **Background**: Current quantum computing technology is in the "Noisy Intermediate-Scale Quantum" (NISQ) era, where hardware is noisy and error-prone. Therefore, hybrid quantum-classical algorithms have become an important means for effective quantum chemical calculations in the near term. 2. **Methods**: - Using the VQE algorithm to compute electronic ground and excited state properties on a quantum computer. - Ensuring wavefunction orthogonality and computing multiple excited states through the VQD algorithm. - Calculating non-adiabatic coupling vectors and transition dipole moments for molecular dynamics simulations. 3. **Validation**: The proposed method's effectiveness was validated by studying the cis-trans isomerization of methanimine and the ultrafast electronic relaxation process of ethylene. The results showed that the method could produce qualitatively accurate molecular dynamics results consistent with experimental data and other computational studies. In summary, this paper presents a significant advancement in using hybrid quantum-classical algorithms on quantum computers for non-adiabatic dynamics simulations, opening new possibilities for future research in chemistry and materials science.