Optimizing the number of measurements for vibrational structure on quantum computers: coordinates and measurement schemes

Vibrational ADAPT-VQE: Critical points lead to problematic convergence

Quantum chemistry is one of the most promising applications for which quantum computing is expected to have a significant impact. Despite considerable research in the field of electronic structure, calculating the vibrational properties of molecules on quantum computers remains a relatively unexplored field. In this work, we develop a vibrational Adaptive Derivative-Assembled Pseudo-Trotter Variational Quantum Eigensolver (vADAPT-VQE) formalism based on an infinite product representation (IPR) of anti-Hermitian excitation operators of the Full Vibrational Configuration Interaction (FVCI) wavefunction, which allows for preparing eigenstates of vibrational Hamiltonians on quantum computers. In order to establish the vADAPT-VQE algorithm using the IPR, we study the exactness of disentangled Unitary Vibrational Coupled Cluster (dUVCC) theory and show that dUVCC can formally represent the FVCI wavefunction in an infinite expansion. To investigate the performance of the vADAPT-VQE algorithm, we numerically study whether the vADAPT-VQE algorithm generates a sequence of operators that may represent the FVCI wavefunction. Our numerical results indicate frequent appearance of critical points in the wavefunction preparation using vADAPT-VQE. These results imply that one may encounter diminishing usefulness when preparing vibrational wavefunctions on quantum computers using vADAPT-VQE and that additional studies are required to find methods that can circumvent this behavior.

Time-dependent coupled cluster with orthogonal adaptive basis functions: General formalism and application to the vibrational problem

We derive equations of motion for bivariational wave functions with orthogonal adaptive basis sets and specialize the formalism to the coupled cluster Ansatz. The equations are related to the biorthogonal case in a transparent way, and similarities and differences are analyzed. We show that the amplitude equations are identical in the orthogonal and biorthogonal formalisms, while the linear equations that determine the basis set time evolution differ by symmetrization. Applying the orthogonal framework to the nuclear dynamics problem, we introduce and implement the orthogonal time-dependent modal vibrational coupled cluster (oTDMVCC) method and benchmark it against exact reference results for four triatomic molecules as well as a reduced-dimensional (5D) trans-bithiophene model. We confirm numerically that the biorthogonal TDMVCC hierarchy converges to the exact solution, while oTDMVCC does not. The differences between TDMVCC and oTDMVCC are found to be small for three of the five cases, but we also identify one case where the formal deficiency of the oTDMVCC approach results in clear and visible errors relative to the exact result. For the remaining example, oTDMVCC exhibits rather modest but visible errors.

Pitfalls in the n-mode representation of vibrational potentials

Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)-effectively a many-body expansion in the space of molecular vibrations-is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by "painting in" the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes-and restore SCF convergence-in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.