Pre‑Born‑Oppenheimer Molecular Hamiltonian

• Pre‑Born‑Oppenheimer Hamiltonian is the full non‑relativistic many‑particle Coulomb operator that treats electrons and nuclei equally without any a priori separation.
• It incorporates all kinetic terms and Coulomb interactions, enabling simultaneous treatment of rovibrational, electronic, and non‑adiabatic effects while preserving symmetry constraints.
• Variational methods using explicitly correlated Gaussian functions and quantum algorithms demonstrate its potential for precise electron–nuclear simulations and benchmarking approximate models.

The pre‑Born–Oppenheimer molecular Hamiltonian is the full non‑relativistic many‑particle Coulomb Hamiltonian for a molecular system, acting on a single wavefunction of all electrons and nuclei without any a priori electronic–nuclear separation. In this formulation, electrons and nuclei are treated on equal footing as quantum particles with their physical masses and charges, and the Hamiltonian contains all kinetic terms and all Coulomb interactions in one operator. In laboratory‑fixed Cartesian coordinates, a standard form is

$$\hat H = -\sum_i \frac{1}{2m_i}\nabla_{\mathbf r_i}^2 + \sum_{i<j}\frac{q_i q_j}{|\mathbf r_i-\mathbf r_j|},$$

or, with explicit electronic and nuclear partitions,

$$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$

This is the Hamiltonian that underlies pre‑BO molecular structure theory, explicitly correlated Gaussian implementations, multicomponent orbital methods such as NOMO, and recent first‑quantized quantum algorithms for direct electron–nuclear simulation (Matyus, 2018, Muolo et al., 2018, Pocrnic et al., 11 Feb 2026).

1. Definition and relation to the Born–Oppenheimer framework

The defining feature of the pre‑BO Hamiltonian is the absence of any separation into “electronic” and “nuclear” problems. The exact stationary equation is

$\hat H \Psi = E \Psi,$

with $\Psi$ depending on all electronic and nuclear coordinates simultaneously. In the formulation used for rigorous large‑mass analysis, the same structure is written as

$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$

where $\mathcal V(\cdot,X)$ is the electronic operator at fixed nuclear configuration and the nuclear kinetic energy appears explicitly with the nuclear mass parameter $M$ (Bayer et al., 2011).

This differs fundamentally from the standard Born–Oppenheimer construction. In BO theory one first solves, for fixed nuclear positions,

$$\hat H_{\rm el}(\{\mathbf R_A\}) = \hat T_e+\hat V_{ee}+\hat V_{en}(\{\mathbf R_A\})+V_{nn}(\{\mathbf R_A\}),$$

and only afterwards solves a nuclear problem on the resulting potential energy surfaces. By contrast, in the pre‑BO Hamiltonian the nuclear degrees of freedom are dynamic at the Hamiltonian level, electron–nucleus coupling is explicit, and there is no potential‑energy‑surface reduction built into the operator itself (Veis et al., 2015, Matyus, 2018).

A useful consequence is conceptual completeness: rovibrational, electronic, and non‑adiabatic effects are not added in layers but are contained in a single spectrum. A corresponding difficulty is that the full Hamiltonian acts on a very high‑dimensional Hilbert space and must respect translation, rotation, inversion, and permutation symmetries at the many‑particle level (Muolo et al., 2018, Matyus, 2018).

2. Symmetry, coordinates, and internal motion

The pre‑BO Hamiltonian is invariant under overall translations and rotations, and, in the absence of external fields, under inversion. This allows eigenstates to be classified by total spatial angular momentum $N$, its projection $M_N$, parity $$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$0, and spin and permutation quantum numbers appropriate to the particle content (Muolo et al., 2018).

Because the center‑of‑mass motion is physically irrelevant for internal molecular structure, one introduces translationally invariant Cartesian coordinates together with the center‑of‑mass coordinate, so that

$$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$1

In coordinate‑transformed approaches this yields an exact internal Hamiltonian. In laboratory‑fixed Cartesian coordinate implementations, one may instead retain the simple laboratory form of the operator and remove translational contamination analytically from matrix elements. In the explicitly correlated Gaussian treatment of the singlet hydrogen molecule, the translational contamination term was identified as

$$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$2

which is subtracted to recover translationally invariant energies while preserving the simple LFCC form of the Hamiltonian and basis functions (Simmen et al., 2012).

Rotation is more delicate. In explicitly correlated Gaussian approaches, total angular momentum is enforced at the basis level, so the Hamiltonian need not be rewritten in body‑fixed form. In the NOMO translation‑ and rotation‑free construction, by contrast, the working Hamiltonian is

$$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$3

where the center‑of‑mass kinetic energy is “simply subtracted,” while rotational subtraction is only approximate because molecular rotations and vibrations are coupled. In the numerical study of H$$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$4 and HT, $$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$5 was treated by a zeroth‑order Taylor expansion, which was adequate for the lowest rotationless states considered (Veis et al., 2015).

These coordinate and symmetry issues are a common source of misunderstanding. Pre‑BO theory does not require internal coordinates or body‑fixed coordinates as a matter of principle; laboratory‑fixed Cartesian formulations are viable provided translational contamination is removed and rotational symmetry is enforced correctly (Simmen et al., 2012, Simmen et al., 2014).

3. Variational realizations in explicitly correlated bases

A central realization of the pre‑BO Hamiltonian is variational expansion in explicitly correlated Gaussian basis functions. In one widely used form, the spatial basis function is

$$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$6

and the angular‑momentum adapted global‑vector representation augments this with a polynomial and spherical harmonic,

$$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$7

with $$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$8 (Matyus, 2018).

The total wavefunction is then expanded as a linear combination of symmetrized products of spatial and spin functions,

$$\hat H = \hat T_e+\hat T_n+\hat V_{ee}+\hat V_{nn}+\hat V_{en}.$$9

or, in review notation,

$\hat H \Psi = E \Psi,$0

This permits direct variational solution of the full electron–nuclear Schrödinger equation while enforcing the Pauli principle and the permutation symmetry of identical nuclei (Muolo et al., 2018, Matyus, 2018).

Shifted‑center or floating ECGs increase flexibility but are not, in general, eigenfunctions of total angular momentum and parity. One remedy is numerical projection onto irreducible representations of the rotation–inversion group,

$\hat H \Psi = E \Psi,$1

combined with parity projection. In the five‑particle pre‑BO treatment of H$\hat H \Psi = E \Psi,$2, shifted FECGs plus projection yielded a variational upper bound that substantially improved earlier ECG results (Muolo et al., 2018).

The same Hamiltonian supports property calculations and resonance theory. Electric transition dipole moments have been evaluated directly from the full pre‑BO Hamiltonian in laboratory‑fixed Cartesian coordinates, in both length and velocity forms, without relying on clamped‑nuclei transition dipole surfaces (Simmen et al., 2014). Resonances have been accessed by complex coordinate rotation, with the transformed Hamiltonian

$\hat H \Psi = E \Psi,$3

allowing rovibronic resonances of H$\hat H \Psi = E \Psi,$4, Ps$\hat H \Psi = E \Psi,$5, and Ps$\hat H \Psi = E \Psi,$6 to appear as discrete complex eigenvalues (Mátyus, 2018).

4. Orbital, second‑quantized, and quantum‑algorithmic representations

A distinct practical realization is the nuclear‑orbital plus molecular‑orbital framework. In NOMO, electrons and nuclei are represented by one‑particle orbitals, and the translation‑ and rotation‑free Hamiltonian is written in second quantization as

$\hat H \Psi = E \Psi,$7

The indices $\hat H \Psi = E \Psi,$8 run over electronic spin orbitals and $\hat H \Psi = E \Psi,$9 over nuclear spin orbitals. In the complete‑basis and complete‑configuration limit, the authors state that “the NOMO/FCI theory for a complete configuration space is an exact theory” within the TRF framework (Veis et al., 2015).

This multicomponent second‑quantized form can be mapped to qubits. Fermionic particles use Jordan–Wigner or Bravyi–Kitaev mappings; bosonic nuclei admit either a direct boson mapping or the compact boson mapping proposed in that work; and distinguishable nuclei can be treated with one‑qubit occupation variables. In the proof‑of‑principle H$\Psi$0 and HT calculations, the nuclei were treated as distinguishable, which the authors state is justified because exchange interaction between nuclei is negligibly small for the states studied (Veis et al., 2015).

Recent first‑quantized quantum algorithms instead represent each particle on a real‑space grid and simulate the full Hamiltonian

$\Psi$1

In this setting, the Hamiltonian is block‑encoded using swap networks and an alternating‑sign implementation of the Coulomb interaction. For the $\Psi$2 reaction, the reported resource estimate is a Toffoli cost of $\Psi$3 per femtosecond with $\Psi$4 logical qubits, illustrating that the full pre‑BO Hamiltonian has become an explicit target for fault‑tolerant quantum simulation (Pocrnic et al., 11 Feb 2026).

5. Controlled reductions and alternative beyond‑BO Hamiltonians

Although the pre‑BO Hamiltonian is the fundamental operator, much work analyzes controlled reductions derived from it. One rigorous route starts from the full molecular Schrödinger operator and shows that, under a uniform electronic spectral gap and in the large nuclear mass limit, Born–Oppenheimer molecular dynamics approximates stationary Schrödinger observables with

$\Psi$5

for any bounded nuclear observable $\Psi$6 and any $\Psi$7 (Bayer et al., 2011). This does not redefine the pre‑BO Hamiltonian; it quantifies when a reduced BO description is accurate for specific observables.

A second route derives a second‑order effective nuclear Hamiltonian associated with one isolated electronic state. The resulting operator contains the BO potential energy surface, the diagonal Born–Oppenheimer correction, and a mass‑correction tensor. In compact form,

$\Psi$8

with

$\Psi$9

This Hamiltonian is explicitly presented as a reduction of the full electron–nuclear Hamiltonian, not as a replacement for pre‑BO theory itself (Mátyus, 2018).

More recently, phase‑space electronic Hamiltonians have been proposed that depend on both nuclear positions and momenta,

$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$0

These constructions are designed to recover beyond‑BO electronic momentum, current density, and vibrational circular dichroism while conserving total linear and angular momentum along quantum‑classical trajectories (Tao et al., 2024, Tao et al., 2024). They are not, however, fully quantum pre‑BO Hamiltonians: in the authors’ own characterization, nuclei enter as classical phase‑space parameters, so the method is best viewed as a mixed quantum–classical approximation motivated by pre‑BO structure rather than a full pre‑BO treatment.

6. Spectral consequences, observables, and conceptual significance

The spectrum of the pre‑BO Hamiltonian contains bound states, rovibronic excitations, and resonances without invoking potential energy surfaces. This has been demonstrated for systems as different as H$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$1, H$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$2, H$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$3, Ps$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$4, Ps$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$5, and hydrogen isotopomers such as HT (Matyus, 2018, Mátyus, 2018, Veis et al., 2015). In H$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$6 and HT, for example, vibrational excitation appears in the NOMO picture as excitation among nuclear orbitals, with dominant configurations

$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$7

for the ground and first excited rotationless vibrational states, respectively (Veis et al., 2015).

Pre‑BO calculations also serve as benchmarks for approximate non‑adiabatic models. For H$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$8, shifted FECGs with projection produced an extrapolated ground‑state energy close to a non‑adiabatic estimate based on BO, DBOC, and effective‑mass corrections, showing that well‑constructed BO plus non‑adiabatic schemes can approach full pre‑BO results while also clarifying their residual error (Muolo et al., 2018).

At the level of wavefunction structure, the pre‑BO Hamiltonian defines a genuine bipartite electron–nuclear problem. Schmidt decompositions and reduced density matrices therefore provide a natural measure of electronic–nuclear entanglement. In the one‑dimensional H$$\mathcal H(x,X)=\mathcal V(x,X)-\frac{1}{2M}\sum_{n=1}^N \Delta_{X^n},$$9 and Shin–Metiu analyses, the ground BO vibronic state could be “almost separable and non‑entangled,” whereas vibrational excitation increased the entanglement monotonically, avoided crossings strongly enhanced entanglement in the adiabatic BO picture, and Born–Huang superpositions could be more entangled than any individual BO component (Mosquera et al., 29 Aug 2025). A plausible implication is that the pre‑BO Hamiltonian is not only a route to non‑adiabatic spectra but also a direct object for quantifying electron–nuclear correlation.

A persistent conceptual issue is the status of molecular structure. Because exact eigenstates of the full Hamiltonian respect overall symmetry, familiar localized structures do not appear as fixed classical geometries in the operator itself. Instead, one extracts structural information from marginal probability densities and angle distributions; for H$\mathcal V(\cdot,X)$0 and H$\mathcal V(\cdot,X)$1D$\mathcal V(\cdot,X)$2 these reveal shell structures and characteristic angular correlations (Matyus, 2018). This suggests that classical molecular structure is not encoded by adding geometry parameters to the Hamiltonian, but emerges from symmetry‑resolved quantum states, measurement, and, in the broader view of the review literature, decoherence (Matyus, 2018).

The pre‑Born–Oppenheimer molecular Hamiltonian is therefore both a specific operator and a framework. As an operator, it is the full many‑particle Coulomb Hamiltonian with all electrons and nuclei quantum mechanical. As a framework, it encompasses exact internal‑coordinate and laboratory‑coordinate formulations, explicitly correlated variational methods, multicomponent orbital theories, resonance formalisms, and modern quantum algorithms, while also providing the reference point from which controlled BO and beyond‑BO reductions are derived (Matyus, 2018, Simmen et al., 2012, Pocrnic et al., 11 Feb 2026).