One of the primary challenges prohibiting demonstrations of practical quantum advantages for near-term devices amounts to excessive measurement overheads for estimating relevant physical quantities such as ground state energies. However, with major differences between the electronic and vibrational structures of molecules, the question of how the resource requirements of computing anharmonic, vibrational states can be reduced remains relatively unexplored compared to its electronic counterpart. Importantly, bosonic commutation relations, distinguishable Hilbert spaces and vibrational coordinates allow manipulations of the vibrational system that can be exploited to minimize resource requirements. In this work, we investigate the impact of different coordinate systems and measurement schemes on the number of measurements needed to estimate anharmonic, vibrational states for a variety of three-mode (six-mode) molecules. We demonstrate an average of 3-fold (1.5-fold), with up to 7-fold (2.5-fold), reduction in the number of measurements required by employing appropriate coordinate transformations, based on an automized construction of qubit Hamiltonians from a conventional vibrational structure program.
One of the primary challenges prohibiting demonstrations of practical quantum advantages for near-term devices amounts to excessive measurement overheads for estimating relevant physical quantities such as ground state energies. However, with major differences between the electronic and vibrational structures of molecules, the question of how the resource requirements of computing anharmonic, vibrational states can be reduced remains relatively unexplored compared to its electronic counterpart. Importantly, bosonic commutation relations, distinguishable Hilbert spaces and vibrational coordinates allow manipulations of the vibrational system that can be exploited to minimize resource requirements. In this work, we investigate the impact of different coordinate systems and measurement schemes on the number of measurements needed to estimate anharmonic, vibrational states for a variety of three-mode (six-mode) molecules. We demonstrate an average of 3-fold (1.5-fold), with up to 7-fold (2.5-fold), reduction in the number of measurements required by employing appropriate coordinate transformations, based on an automized construction of qubit Hamiltonians from a conventional vibrational structure program.
Abstract Quantum computing has recently been emerging in theoretical chemistry as a realistic avenue meant to offer computational speedup to challenging eigenproblems in the context of strongly-correlated molecular systems or extended materials. Most studies so far have been devoted to the quantum treatment of electronic structure, for which the transformation of fermionic operators into the qubit space is quite transparent. In contrast, only a few were directed to the quantum treatment of vibrational structure, which at the moment remains not devoid of unknowns. When simulating an anharmonic vibrational mode under bosonic second quantization, the biggest problem, which we analyze in detail in the present work, is the disruption of the resolution of the identity when truncating the infinite harmonic-oscillator basis set of primitive modals and how this alters the canonical commutation relation. This may lead to serious incorrectness in the evaluation of the Hamiltonian matrix elements when assembling them from finite matrix products, which eventually entail a nonvariational behavior with respect to the finite size of the computational basis. As we show here, a simple cure occurs to be obtained upon using Wick's normal order, for reasons that are not standardly evoked. We also provide a detailed comparison between the boson-to-qubit unary and binary mappings under two different representations: one for the bosonic ladder operators within the harmonic-oscillator primitive basis, and one for the so-called n-mode representation within any type of computational basis. In addition, we discuss the impact of choosing an adequate primitive basis set in terms of quantum computing with respect to its variational convergence efficiency (number of basis functions, hence of qubits) and as regards the magnitude of the 1-norm of the encoded Hamiltonian (a measure of the computational complexity of the quantum algorithm). Such fundamental aspects are illustrated numerically on a one-dimensional anharmonic Hamiltonian model corresponding to a symmetric double-well potential, of interest both for vibrational spectroscopy and chemical reactivity, and which is a challenging situation for numerical convergence due to fine tunneling splitting.
Accurate simulations of vibrational molecular spectra are expensive on conventional computers. Compared to the electronic structure problem, the vibrational structure problem with quantum computers is less investigated. In this work we accurately estimate quantum resources, such as number of logical qubits and quantum gates, required for vibrational structure calculations on a programmable quantum computer. Our approach is based on quantum phase estimation and focuses on fault-tolerant quantum devices. In addition to asymptotic estimates for generic chemical compounds, we present a more detailed analysis of the quantum resources needed for the simulation of the Hamiltonian arising in the vibrational structure calculation of acetylene-like polyynes of interest. Leveraging nested commutators, we provide an in-depth quantitative analysis of trotter errors compared to the prior investigations. Ultimately, this work serves as a guide for analyzing the potential quantum advantage within vibrational structure simulations.
Quantum computation of vibrational properties of molecules is a promising platform to obtain computational advantages for computational chemistry. However, fault-tolerant quantum computations of vibrational properties remain a relatively unexplored field in quantum computing. In this work, we present different algorithms for efficient encodings of vibrational Hamiltonians using qubitization. Specifically, we investigate different encoding representations, high-order tensor decomposition to obtain low rank approximations for the vibrational Hamiltonian, rectilinear, and polyspherical coordinate systems, parallelization, and grouping algorithms. To investigate the performance of the different methods, we perform benchmark computations for both small and large molecules with more than one hundred vibrational modes with two- and three-mode coupled Hamiltonian. We present the first resource estimates for qubitization using multimode vibrational Hamiltonians with 1 cm –1 spectral resolution for molecules with up to one hundred vibrational modes. Additionally, we show that tensor decomposition and operator parallelization effectively reduce the T gate depth of the quantum algorithm.
Determining the vibrational structure of a molecule is central to fundamental applications in several areas, from atmospheric science to catalysis, fuel combustion modeling, biochemical imaging, and astrochemistry. However, when significant anharmonicity and mode coupling are present, the problem is classically intractable for a molecule of just a few atoms. Here, we outline a set of quantum algorithms for solving the molecular vibrational structure problem for both near- and long-term quantum computers. There are previously unaddressed characteristics of this problem which require approaches distinct from most instances of the commonly studied quantum simulation of electronic structure: many eigenstates are often desired, states of interest are often far from the ground state (requiring methods for ``zooming in'' to some energy window), particle conservation is not required, and transition amplitudes with respect to a nonunitary Hermitian operator must be calculated. We address each of these hurdles and consider problem instances of four molecular vibrational Hamiltonians. Finally and most importantly, we give analytical and numerical results which suggest that, to a given energy precision, a vibrational problem instance will be simulatable on a quantum computer before an electronic structure problem instance. These results imply that more focus in the quantum information community ought to shift toward scientifically and industrially important quantum vibrational problems.
We developed an efficient approach to study the coherent control of vibrational state-to-state transitions. The approximations employed in our model are valid in the regime of the low vibrational excitation specific to the vibrational quantum computer. Using this approach we explored how the vibrational properties of a two-qubit system affect the accuracy of subpicosecond quantum gates. The optimal control theory and numerical propagation of laser-driven vibrational wave packets were employed. The focus was on understanding the effect of the three anharmonicity parameters of the system. In the three-dimensional anharmonicity parameter space we identified several spots of high fidelity separated by low fidelity planar regions. The seemingly complicated picture is explained in terms of interferences between different state-to-state transitions. Very general analytic relationships between the anharmonicity parameters and the frequencies are derived to describe the observed features. Geometrically, these expressions represent planes in the three-dimensional anharmonicity parameter space. Results of this work should help to choose a suitable candidate molecule for the practical implementation of the vibrational two-qubit system.
Simulating vibrational dynamics is essential for understanding molecular structure, unlocking useful applications such as vibrational spectroscopy for high-fidelity chemical detection. Quantum algorithms for vibrational dynamics are emerging as a promising alternative to resource-demanding classical approaches, but this domain is largely underdeveloped compared with quantum simulations of electronic structure. In this work, we describe in detail three distinct forms of the vibrational Hamiltonian: canonical bosonic quantization, real space representation, and the Christiansen second-quantized form. Leveraging the Lie algebraic properties of each, we develop efficient fragmentation schemes to enable the use of Trotter product formulas for simulating time evolution. We introduce circuits required to implement time evolution in each form and highlight factors that contribute to the simulation cost, including the choice of vibrational coordinates. Using a perturbative approach for the Trotter error, we obtain tight estimates of the T gate cost for the simulation of time evolution in each form, enabling their quantitative comparison. Combining tight Trotter error estimates and efficient fragmentation schemes, we find that for the CH 4 molecule with 9 vibrational modes, time evolution for approximately 1.8 ps may be simulated using as few as 36 qubits and approximately 3 × 10 8 T gates─an order-of-magnitude speedup over state-of-the-art algorithms. Finally, we present calculations of vibrational spectra using each form to demonstrate the fidelity of our algorithms. This work presents a unified and highly optimized framework that makes simulating vibrational dynamics an attractive use case for quantum computers.