Molecular quantum dynamics simulations that treat both electrons and nuclei quantum mechanically are crucial for predicting chemical reactions. With classical computation, a full wave-function representation of both types of particles requires resources that grow exponentially with their number. Oscillator-qubit processors, including trapped-ion and circuit QED devices, represent electronic states with qubits and nuclear vibrations with oscillators. Their native operations produce linear vibronic interactions, from which programmable anharmonic interactions can be built. Product-formula approaches build the corresponding phase gates by repeating small noncommuting phase steps, and quantum signal processing (QSP) approaches approximate each gate directly with a Fourier series. Simulations apply many such gates in sequence, so accumulated cost and error limit how long the dynamics can be followed. We propose near-optimal quantum primitives for molecular quantum dynamics simulation, covering the representations of both types of particles and their interactions. We synthesize nonlinear and multimode bosonic phase gates using generalized Jacobi-Anger (GJA) expansions within QSP, with a cost that grows near-optimally in both interaction strength and precision. Our interaction model and gate synthesis each have their own accuracy setting. Chosen together, these settings make our multimode benchmark circuits about one-third shorter than those of the previous direct QSP approach. We then simulate energy exchange between the stretching and bending vibrations of a simplified carbon dioxide model. With circuit-QED noise, a short simulation narrowly misses our accuracy criteria at qubit and oscillator lifetimes near 2 ms and 5 ms. These lifetimes mark a concrete checkpoint on the way to simulating anharmonic multimode molecular quantum dynamics on near-term oscillator-qubit processors.
Molecular quantum dynamics simulations that treat both electrons and nuclei quantum mechanically are crucial for predicting chemical reactions. With classical computation, a full wave-function representation of both types of particles requires resources that grow exponentially with their number. Oscillator-qubit processors, including trapped-ion and circuit QED devices, represent electronic states with qubits and nuclear vibrations with oscillators. Their native operations produce linear vibronic interactions, from which programmable anharmonic interactions can be built. Product-formula approaches build the corresponding phase gates by repeating small noncommuting phase steps, and quantum signal processing (QSP) approaches approximate each gate directly with a Fourier series. Simulations apply many such gates in sequence, so accumulated cost and error limit how long the dynamics can be followed. We propose near-optimal quantum primitives for molecular quantum dynamics simulation, covering the representations of both types of particles and their interactions. We synthesize nonlinear and multimode bosonic phase gates using generalized Jacobi-Anger (GJA) expansions within QSP, with a cost that grows near-optimally in both interaction strength and precision. Our interaction model and gate synthesis each have their own accuracy setting. Chosen together, these settings make our multimode benchmark circuits about one-third shorter than those of the previous direct QSP approach. We then simulate energy exchange between the stretching and bending vibrations of a simplified carbon dioxide model. With circuit-QED noise, a short simulation narrowly misses our accuracy criteria at qubit and oscillator lifetimes near 2 ms and 5 ms. These lifetimes mark a concrete checkpoint on the way to simulating anharmonic multimode molecular quantum dynamics on near-term oscillator-qubit processors.
Accurately simulating molecular dynamics during chemical reactions remains a central challenge in computational chemistry. These processes are inherently quantum mechanical, involving coupled electronic and nuclear motion and non-adiabatic effects. A fully quantum description is therefore required, yet the associated computational cost grows exponentially with system size. Even with modern supercomputers, storing and evolving a molecular wavefunction rapidly becomes intractable, limiting exact simulations of chemical dynamics. Quantum simulation offers a promising alternative by using controllable quantum systems to directly reproduce molecular dynamics. However, fully digital, qubit-based approaches remain impractical for chemistry due to limited qubit numbers, coherence times, and gate fidelities. This motivates the exploration of more resource-efficient quantum architectures. In this thesis, we demonstrate that mixed-qudit-boson (MQB) quantum simulators provide a practical platform for simulating chemical dynamics in the quantum regime. In MQB devices, multilevel qudits encode electronic states while bosonic modes natively represent vibrational degrees of freedom, avoiding costly qubit encodings of nuclear motion and substantially reducing resource requirements. The thesis makes three main contributions. First, it reports the direct observation of geometric-phase interference in non-adiabatic dynamics using an MQB simulation of the Jahn–Teller model. Second, it extends the MQB framework to open-system chemical dynamics through programmable Lindblad dissipation. Third, it presents a resource analysis showing that MQB simulators can reduce resource costs by several orders of magnitude relative to qubit-only approaches. Together, these results establish MQB quantum simulators as a promising route for studying chemical dynamics beyond the reach of classical computation.
Simulating chemical dynamics is computationally challenging, especially for nonadiabatic dynamics, where numerically exact classical simulations scale exponentially with system size, becoming intractable for even small molecules. On quantum computers, chemical dynamics can be simulated efficiently using either universal, qubit-only devices or specialized mixed-qudit-boson (MQB) simulators, which natively host electronic and vibrational degrees of freedom. Here, we compare the quantum resources required for a qubit-only approach to achieve the same accuracy as an MQB device at simulating nonadiabatic molecular dynamics. We find that MQB simulations require orders-of-magnitude fewer quantum operations than qubit-only simulations, with a one-gate MQB circuit requiring a qubit-equivalent circuit volume of over 400,000 when simulating an isolated molecule, which increases to over ten million when environmental effects are included. These estimates assume perfect qubits and gates, and would increase by additional orders of magnitude if error correction were used for fault tolerance. When errors are small, the advantage of MQB simulators becomes even larger as system size increases. Our results highlight the enormous resource advantages of representing non-qubit chemical degrees of freedom natively, rather than encoding them into qubits.
Hybrid oscillator-qubit processors have recently demonstrated high-fidelity control of both continuous- and discrete-variable information processing. However, most quantum algorithms remain limited to homogeneous quantum architectures. Here, we present a compiler for hybrid oscillator-qubit processors that implements state preparation and time evolution. In this setting, the compiler invokes generalized quantum signal processing (GQSP) to synthesize arbitrary analytic bosonic phase gates in a constructive manner, with circuit depth scaling as . The approximation cost scales with the Fourier bandwidth of the target bosonic phase, rather than with the degree of nonlinearity. Armed with OQ-GQSP, multi-state vibronic coupling Hamiltonian dynamics can be decomposed into state-dependent arbitrary-phase potential propagators, which also enable extension to multi-state systems through the parity-measurement technique. Compared to fully discrete encodings, our approach avoids the overhead of truncating continuous variables, resulting in linear dependence on the number of vibrational modes. We validate our method on the uracil cation, a canonical system whose accurate modeling requires anharmonic vibronic models, and estimate the cost of state preparation and time evolution.
Key to solving many societally relevant problems in materials science, quantum chemistry, and nuclear physics is understanding their fundamental microscopic quantum behavior. All of these problems can be phrased as fermionic (e.g., electrons, quarks) and bosonic particles (e.g., lattice vibrations, photons) interacting with each other. However, solving these systems of strongly interacting fermions and bosons is out of reach for classical computers and also non-error-corrected qubit-only quantum computers. In this paper, we introduce a hybrid qubit and bosonic-oscillator quantum-simulation framework, which is so efficient that it could simulate some of these strongly correlated problems even without error correction. Specifically, we leverage the available universal bosonic gate set with O ( 1 ) operations per gate as a function of the bosonic cutoff. This enables our framework to provide exact decompositions of higher-order fermion-boson interactions, including notoriously challenging terms such as ( 2 + 1 ) D gauge theories. We show how to compile dynamics, ancilla-free partial error detection, nonlocal observable measurements, and ground-state energy estimation using both qubit-boson variational quantum eigensolvers and quantum signal processing. We introduce a “post-qubit quantum advantage” experiment which is out of reach for both near-term qubit-only quantum computers and classical computers. We justify this proposal by estimating bosonic gate fidelities through high–boson-occupation Lindblad simulation and by comparing to qubit-only algorithms. For the latter, we develop a new qubit-based algorithm for simulating a beam splitter with a boson number cutoff N max , which scales as O ˜ ( N max 1 + o ( 1 ) ) , asymptotically outperforming previously known qubit methods. However, to reach the same circuit fidelity as qubit-boson hardware with 99.9% per-gate fidelity, a qubit quantum computer using the best-known method would still need fidelities of at least 99.9999%, likely requiring substantial error correction. Our framework can be used in superconducting hardware, which we focus on, as well as trapped-ion and neutral-atom hardware. This work establishes digital hybrid qubit-boson quantum simulation as a promising and advantageous method for solving the many challenging, strongly correlated, fermion-boson models of relevance to science and technology.
Modern computational techniques used to simulate quantum chemistry are on the boundary of tractability due to the exponential growth of the molecular wavefunction, requiring a careful balance between molecule size and simulation accuracy. In recent years, quantum computing has risen in popularity as a potential alternative to conventional (classical) techniques; however, most methods rely on access to "digital" quantum computers composed of qubits and quantum gates, which at present are severely limited by noise. We have developed a real-time, analog approach to simulate vibronic chemical dynamics with existing quantum technology. Our approach uses an intuitive mapping of molecular electronic and vibrational degrees of freedom onto quantum resonators and qudit (d-level system) states, with controllable couplings between degrees of freedom. The measurement output can be mapped onto different time-dependent observables, including the time-domain simulation of vibronic spectra. Our approach can also incorporate controlled sources of noise to simulate system-bath interactions and dissipative dynamics at a minimal cost. We present experimental results using a trapped-ion device, thus showing the potential for near-term simulation of chemical dynamics in complex environments beyond the abilities of classical computers.
A key requirement to perform simulations of large quantum systems on near-term quantum hardware is the design of quantum algorithms with short circuit depth that finish within the available coherence time. A way to stay within the limits of coherence is to reduce the number of gates by implementing a gate set that matches the requirements of the specific algorithm of interest directly in hardware. Here, we show that exchange-type gates are a promising choice for simulating molecular eigenstates on near-term quantum devices since these gates preserve the number of excitations in the system. Complementing the theoretical work by Barkoutsos et al. [PRA 98, 022322 (2018)], we report on the experimental implementation of a variational algorithm on a superconducting qubit platform to compute the eigenstate energies of molecular hydrogen. We utilize a parametrically driven tunable coupler to realize exchange-type gates that are configurable in amplitude and phase on two fixed-frequency superconducting qubits. With gate fidelities around 95% we are able to compute the eigenstates within an accuracy of 50 mHartree on average, a limit set by the coherence time of the tunable coupler.