Spectroscopy and Scattering Studies Using Interpolated Ab Initio Potentials

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.

Orthogonal State Reduction Variational Eigensolver for the Excited-State Calculations on Quantum Computers

Variational quantum eigensolver (VQE) is a promising method for ground-state calculations on current noisy intermediate-scale quantum computers. However, the research progress of excited-state calculations on quantum computers is relatively slow. In order to extend the framework of VQE for excited-state calculations, we propose a new algorithm, orthogonal state reduction variational eigensolver (OSRVE), to determine the energies of excited states. Theoretical derivations prove that the optimized state in the OSRVE method can ensure the energy minimum and orthogonality constraint simultaneously, and OSRVE is also applicable for the degenerate state. The performance of OSRVE is demonstrated by numerical calculations of the H₄ and H₂O molecules. Compared with other excited-state calculation algorithms, OSRVE has obvious advantages in calculating lower-order excited states. This work extends the VQE algorithm to excited-state calculations, and OSRVE can be implemented on near-term quantum computers.