Pre-Born-Oppenheimer Quantum Simulation

• Pre-Born-Oppenheimer quantum simulation is a framework treating electrons and nuclei as coupled quantum particles that bypasses the Born–Oppenheimer approximation to capture nonadiabatic interactions.
• It employs explicitly correlated Gaussian functions and variational techniques to achieve spectroscopic accuracy in modeling few-body and highly correlated molecular systems.
• Quantum algorithms and analog simulation platforms integrated within this framework enable efficient resource scaling and precision computation for complex molecular dynamics.

Pre-Born-Oppenheimer (pre-BO) quantum simulation refers to the ab initio treatment of molecular systems in which all electrons and nuclei are treated as quantum particles on the same footing. Unlike conventional methods that rely on the Born–Oppenheimer (BO) approximation to separate fast electronic motion from slower nuclear dynamics, pre-BO formulations avoid this decoupling. This enables direct access to genuine electron-nuclear and nonadiabatic effects—critical in resonant, ultrafast, or highly correlated systems. Recent advances have driven the development of pre-BO theory, specialized numerics, and quantum algorithms for practical simulation of molecules, paving the way for high-precision computation and experimental quantum simulation beyond the reach of BO-based methods.

1. Fundamental Formalism of Pre-Born-Oppenheimer Theory

The pre-BO molecular Hamiltonian incorporates both electronic and nuclear degrees of freedom as quantum variables, eschewing any a priori partition into electronic and nuclear subspaces. For an $(n+1)$-particle, nonrelativistic system (electrons and nuclei), the Hamiltonian in atomic units is

$$\hat H = -\sum_{i=1}^{n+1} \frac{1}{2m_i} \Delta_{r_i} + \sum_{i=1}^{n+1}\sum_{j>i}^{n+1} \frac{q_i q_j}{|r_i - r_j|}$$

where $r_i$ are Cartesian coordinates, $m_i$ and $q_i$ are masses and charges, and $\Delta_{r_i}$ denotes the Laplacian with respect to $r_i$ (Simmen et al., 2014).

Pre-BO theory requires fully correlated wavefunctions,

$$\Psi(r_1, \ldots, r_{n+1}) = \sum_I c_I\, \phi_{L_I M_{L,I}}(r; A_I, u_I, K_I) \otimes \chi_{\rm spin}$$

using explicitly correlated Gaussian (ECG) basis functions in the laboratory-fixed frame. Angular momentum, parity, and (anti-)symmetry are enforced via the global vector representation and spin eigenfunctions (Matyus, 2018, Simmen et al., 2014).

Physical observables, transition amplitudes, and response functions must be evaluated with respect to these fully entangled states, yielding results that reflect true coupled electron–nuclear physics, including non-adiabatic couplings, rovibronic resonances, and quantum entanglement between all particles.

2. Numerical Methods: Explicitly Correlated Gaussians, Projections, and Variational Approaches

Most high-accuracy pre-BO structure calculations utilize ECGs, which provide a flexible, analytical machinery for few-body systems. An ECG with a shifted center is defined as

$$\varphi(\mathbf{r}) = \exp\left[-(\mathbf{r} - \mathbf{s})^T A (\mathbf{r} - \mathbf{s})\right], \quad A \succ 0,~\mathbf{s} \in \mathbb{R}^{3N_p}$$

where the shift vector allows basis functions to follow nuclear positions, critical for flexible description of floppy or highly delocalized systems (Muolo et al., 2018).

Specialized projection schemes numerically enforce total angular momentum and parity by integrating over rotational group elements, allowing construction of basis sets with well-defined quantum numbers: $$\varphi^{[N]}_{M_N}(\mathbf{r}) = \hat{P}^{[N]}_{M_N M_N} \varphi(\mathbf{r}), \qquad \hat{P}^{[N]}_{M_1 M_2} = \int \frac{d\Omega}{4\pi^3} D_{M_1 M_2}^{[N]}(\Omega)^* \hat{R}(\Omega)$$ While these methods increase computational cost due to required angular integrations, the approach is systematically improvable and essential for accurate rovibronic energy levels and spectra (Muolo et al., 2018).

Variational methods using large non-orthogonal ECG sets, optimized via stochastic or gradient-free procedures, achieve spectroscopic accuracy for bound states and resonances. Complex coordinate rotation techniques enable the calculation of resonance positions and lifetimes for metastable rovibronic states (Mátyus, 2018).

3. Quantum Computing Protocols for Pre-Born-Oppenheimer Simulation

Quantum simulation of pre-BO dynamics leverages various algorithmic paradigms:

• First-Quantized Real-Space/DVR Encoding: The universal form of the pre-BO Hamiltonian allows direct grid encoding. Recent block-encoding schemes using swap networks and alternating-sign representations of the Coulomb operator achieve linear scaling in particle number for the block-encoding Toffoli count, enabling efficient simulation of large systems such as $\rm NH_3+BF_3$ with $8.7 \times 10^9$ Toffolis per femtosecond and 1362 logical qubits (Pocrnic et al., 11 Feb 2026). These approaches are best-in-class for cost to date.
• Plane-Wave and Interaction-Picture Simulation: Gate complexity as low as $\tilde{O}(N^{1/3} \eta^{8/3})$ is possible with plane-wave basis simulations in the interaction picture, advantageous in regimes where fine basis discretization is critical to resolve both electron and nuclear cusps ($N \gg \eta$) (Babbush et al., 2018).
• Quantum Variational Methods: Highly compact ansätze such as NI-DUCC-VQE with Minimal Complete Pool (MCP) of Lie algebraic excitations avoid barren plateaus and achieve rapid convergence to $\sim 10^{-11}$ atomic units energy errors in three-body pre-BO calculations, with resource requirements scaling only logarithmically in basis size (Haidar et al., 20 Oct 2025).
• Second-Quantized and NEO Approaches: Multicomponent (electron + nucleus) exact Hamiltonians in a finite orbital basis are mapped via fermion-to-qubit transformations (Jordan–Wigner, Bravyi–Kitaev, parity mapping), enabling multicomponent coupled-cluster and unitary coupled-cluster (mcUCC) simulations (Kovyrshin et al., 2023, Cabral et al., 14 Nov 2025). Resource reduction by two-qubit tapering and symmetry exploitation enables full pre-BO treatments on near-term devices, with errors $<10^{-6}$ Ha in ground-state energies for small molecules such as $\mathrm{H}_2$ and malonaldehyde.

4. Analog and Hybrid Quantum Simulation Platforms

Analog quantum simulation in devices with bosonic modes (trapped ions, cavity/circuit QED, optical lattices) enables direct mapping of the pre-BO Hamiltonian—including linear and quadratic vibronic couplings, conical intersections, and nonadiabatic dynamics—onto hardware. Key architectures include:

• Mixed Qudit–Boson (MQB) and Coupled Qubit–Bosonic (cMQB): Simulation of coupled electron-nuclear dynamics is realized by encoding electrons in qubits or qudits and nuclear vibrations in bosonic modes. The mapping preserves entanglement and supports Trotterized time evolution or continuous evolution under the full vibronic Hamiltonian. Device resource scaling is linear in orbitals and vibrational modes, achieving exponential savings over classical and even digital quantum methods (MacDonell et al., 2020, Ha et al., 2024).
• Optical Lattice Emulation: Utilizing two rotational states of ultracold fermionic molecules, electron and nuclear degrees of freedom are mapped onto fermionic modes. Tunable hopping and dipolar interactions realize electron–electron and electron–nuclear couplings. These platforms qualitatively capture non-adiabatic scattering, exchange, and ionization phenomena, and the mass ratio of “nuclei” to “electrons” can be tuned dynamically (Argüello-Luengo et al., 30 Mar 2025).

5. Observables, Transition Dipoles, and Benchmark Results

Pre-BO simulations provide direct access to physically meaningful observables:

• Transition Dipole Moments: Explicit non-adiabatic evaluation of transition dipole moments between fully correlated rovibronic states is achieved by deriving closed-form expressions for matrix elements in the ECG basis, both in length and velocity gauges. Benchmark calculations for $\mathrm{H}_2$ yield $|\mu_{X\to B}| = 0.078414$ a.u., agreeing with clamped-nuclei treatments to better than 1%—demonstrating the precision gain when full electron–nuclear coupling is included (Simmen et al., 2014).
• Resonance Energies and Widths: Resonance (quasi-bound) rovibronic states embedded in dissociation continua can be characterized by complex-coordinate-rotation techniques, giving resonance energies $E$ and lifetimes $\tau = \hbar / \Gamma$. Bound and resonance energies for systems such as $\mathrm{Ps}_2$, $\mathrm{H}_2$ reproduce or exceed the precision of adiabatic methods (Mátyus, 2018).
• Dynamical Observables: Protocols for extracting one- and two-particle densities, current correlations, and dynamical structure factors have been established in both quantum trajectory-based frameworks and grid-based quantum simulations, allowing real-time tracking of non-adiabatic charge, energy, and coherence transfer processes (Larder et al., 2018, Ha et al., 2024).

6. Error Mitigation, Resource Scaling, and Implementation Challenges

Pre-BO quantum simulation faces significant challenges due to the exponential growth of Hilbert space with particle and basis size, requirement for antisymmetrization over all quantum particles, and the need for high-precision arithmetic. To address these:

• Error Mitigation: Physics-Inspired Extrapolation (PIE), symmetry-based qubit tapering, and low-depth ansätze such as Local Unitary Cluster Jastrow (LUCJ) compress circuit and measurement overhead while maintaining chemical accuracy—demonstrated in hardware with $\mathcal{O}(10^{-6})$ Ha errors (Cabral et al., 14 Nov 2025).
• Resource Efficiency: Leading algorithms achieve linear or even sublinear scaling in key bottlenecks such as number of two-body interactions, gate count per Trotter step, ancilla overhead, and qubit requirements (e.g., Toffoli-depth $O(\eta^3 t)$ and logical qubit count of order $10^3$–$10^4$ for 50-particle systems (Pocrnic et al., 11 Feb 2026)).
• Measurement and Decoherence: Advanced strategies for symmetry exploitation, state preparation (parameter transfer, VQE-based initialization), and error modeling in open system dynamics (Lindblad techniques) permit scaling to larger and more complex systems (Kovyrshin et al., 2023, MacDonell et al., 2020).

7. Scope, Applications, and Outlook

Pre-Born–Oppenheimer quantum simulation is essential for:

• Describing photochemical reactions, ultrafast dynamics, and proton-coupled electron transfer where non-adiabatic effects dominate.
• Benchmarking and validating new quantum devices—three-body and few-electron/multinuclear systems serve as stringent tests for algorithmic and hardware fidelity (Haidar et al., 20 Oct 2025).
• Developing and validating new approximations (e.g., coupled electron-nuclear Kohn–Sham or multi-component CC) that systematically recover pre-BO physics.
• Laying a foundation for quantum-enhanced metrology, precision spectroscopy, and fundamental studies in chemistry and physics (e.g., positronium complexes, isotope effects, molecule–positron interactions) (Matyus, 2018).

The major outstanding challenges are the extension to larger molecular systems (necessitating sublinear-scaling methods and hardware-efficient encodings), efficient handling of nuclear bosonic modes in strongly anharmonic systems, and automated selection of active quantum subspaces for scalable embedded simulation. Nevertheless, the rapid progress in algorithmic innovation, resource optimization, and device development strongly indicates that pre-Born–Oppenheimer quantum simulation will become a central methodology for first-principles molecular science in the coming decade.