Exact tunneling splittings of rotationally excited states from symmetrized path-integral molecular dynamics

Jeremy O. Richardson; Lea Zupan; Yu-Chen Wang

arxiv: 2604.12638 · v1 · submitted 2026-04-14 · ⚛️ physics.chem-ph

Exact tunneling splittings of rotationally excited states from symmetrized path-integral molecular dynamics

Lea Zupan , Yu-Chen Wang , Jeremy O. Richardson This is my paper

Pith reviewed 2026-05-10 14:22 UTC · model grok-4.3

classification ⚛️ physics.chem-ph

keywords tunneling splittingpath-integral molecular dynamicsrotational excitationammonia inversionEckart framesymmetrized ring polymerquantum dynamics

0 comments

The pith

A path-integral method projects ring polymers onto specific rotational states to extract exact tunneling splittings for multiple J values from single simulations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends symmetrized path-integral molecular dynamics by adding an Eckart spring that links the ring-polymer ends through a permutation-inversion-rotation operation. This projects the sampled configurations onto chosen rotational manifolds and symmetries, yielding tunneling splittings that are numerically exact once all parameters converge. The same simulation run simultaneously produces results for several total angular-momentum quantum numbers J without extra cost. Validation on water recovers rotational levels beyond the rigid-rotor limit, and application to ammonia reproduces exact benchmarks while reproducing the experimental trend of decreasing splitting with rising J.

Core claim

The symmetrized path-integral molecular dynamics formalism, augmented with an Eckart spring, projects the ring polymer onto selected rotational manifolds and states of chosen symmetry, delivering numerically exact tunneling splittings for rotationally excited states within statistical uncertainty; the same trajectories furnish splittings for multiple J values at no added cost and recover the observed decrease of splitting with J for ammonia.

What carries the argument

The Eckart spring, a harmonic restraint that connects the two end beads of the ring polymer through a chosen permutation-inversion-rotation operation to enforce projection onto a specific rotational manifold and symmetry.

If this is right

Splittings for all relevant J can be obtained from one simulation set rather than separate runs for each rotational level.
The method applies directly to other molecules exhibiting inversion or proton-transfer tunneling without requiring rigid-rotor approximations.
Rotationally resolved spectra become accessible for systems where full quantum variational calculations remain prohibitive.
Trends such as the systematic reduction of splitting with increasing J can be mapped systematically across isotopologues or vibrational states.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The approach may scale to larger polyatomic systems or higher excitations where matrix-diagonalization methods become intractable.
If the same Eckart-spring projection is combined with more accurate or ab-initio surfaces, it could predict splittings for states not yet measured experimentally.
The observed J-dependence may prove general for inversion tunneling in symmetric tops, offering a simple diagnostic for the quality of computed surfaces.

Load-bearing premise

The underlying potential energy surface must be accurate enough and all simulation parameters (beads, time step, spring constants) must be sufficiently converged for the target systems.

What would settle it

A set of converged simulations for ammonia that produce a tunneling splitting for any J differing from the exact variational benchmark by more than the reported statistical uncertainty would falsify the claim of numerical exactness.

Figures

Figures reproduced from arXiv: 2604.12638 by Jeremy O. Richardson, Lea Zupan, Yu-Chen Wang.

**Figure 2.** Figure 2: FIG. 2. Rotation–inversion energy levels of ammonia, labeled [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Difference between the rotationally excited and [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗

read the original abstract

We extend our previous symmetrized path-integral molecular dynamics approach to calculate tunneling splittings of molecules in rotationally excited states. In this new formalism, the system is rigorously projected onto selected rotational manifolds and states of a chosen symmetry through an Eckart spring, which connects the two end beads of the ring polymer via a permutation--inversion--rotation operation. This method is numerically exact within statistical uncertainty once convergence with respect to all simulation parameters has been achieved. Importantly, it enables the simultaneous extraction of tunneling splittings for multiple total angular-momentum quantum numbers $J$ from a single set of simulations, without additional computational cost relative to the original approach. After validating the formalism by computing the rotational levels of water (beyond the rigid-rotor approximation), we apply it to ammonia and obtain rotationally resolved tunneling splittings in excellent agreement with exact variational benchmarks. Except for small errors due to the underlying potential energy surface, the results capture the experimentally observed trend that the tunneling splitting decreases with $J$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The Eckart-spring projection extends symmetrized PIMD to rotationally excited tunneling splittings and extracts multiple J values from one run, with validation that holds up against benchmarks.

read the letter

The key advance here is the Eckart spring that projects the ring polymer onto specific rotational manifolds via a permutation-inversion-rotation operation. This lets them pull tunneling splittings for several J values simultaneously from the same simulation set, at no added cost over the prior symmetrized PIMD method. They first validate the rotational projection on water beyond the rigid-rotor limit, then apply it to ammonia and recover splittings that match exact variational results within PES errors while also capturing the experimental trend of decreasing splitting with J. The derivation stays grounded in symmetry and Eckart conditions rather than fitted parameters, which keeps the extension clean and independent of their earlier work. The claim that the method is numerically exact within statistical uncertainty once all parameters converge is supported by the reported agreement. What stands out is the practical payoff: one set of runs gives the full rotational manifold data. The soft spots are the standard ones for path-integral work. Results still depend on having an accurate enough PES and on reaching convergence in beads, timestep, and spring strength for the target system. The authors flag the PES contribution to small discrepancies, and nothing in the presented checks suggests the projection itself introduces bias or hidden fitting. For larger molecules the ergodic sampling could become the practical limit, but that is not a flaw in the formalism shown here. This paper is for computational chemists already using or building on path-integral methods for spectroscopy and tunneling. A reader in that area gets a usable extension with clear implementation details and benchmark numbers to check against. The thinking is direct, the checks are honest, and the central argument stands without internal contradictions, so it deserves a serious referee. I would send it out for peer review.

Referee Report

0 major / 3 minor

Summary. The paper extends symmetrized path-integral molecular dynamics to rotationally excited states by projecting the ring polymer onto selected rotational manifolds via an Eckart spring that connects the end beads through a permutation-inversion-rotation operation. It claims the resulting method is numerically exact within statistical uncertainty once all parameters (beads, timestep, spring strength) converge, enables simultaneous extraction of tunneling splittings for multiple J from a single simulation at no extra cost, and validates this on water rotational levels (beyond rigid-rotor) and ammonia, where the computed J-dependent splittings agree with exact variational benchmarks (within PES errors) and reproduce the experimental trend of decreasing splitting with increasing J.

Significance. If the convergence and projection claims hold, the work provides a parameter-free, symmetry-based route to rovibrational tunneling splittings that avoids separate simulations per J. The rigorous derivation from Eckart conditions and the direct validation against variational results for both water and ammonia are strengths; the method could be useful for larger systems where full variational calculations are prohibitive, provided the underlying PES is reliable.

minor comments (3)

[Abstract and §4.3] The abstract and results sections state agreement 'within PES errors' and 'excellent agreement' with variational benchmarks, but do not tabulate the separate contributions of statistical uncertainty versus PES discrepancy for each J; adding this breakdown (e.g., in Table 2 or 3) would strengthen the 'numerically exact upon convergence' claim.
[§4.2 and Figure 3] Convergence with respect to the number of beads and the Eckart-spring strength is asserted for the ammonia calculations, but the supplementary figures or main text do not show explicit plateauing for J>1; a dedicated convergence panel for the highest J reported would address the weakest assumption noted in the review.
[§2.2] The definition of the permutation-inversion-rotation operator in the projection (Eq. (X)) uses standard symmetry labels, but a brief explicit mapping to the molecular symmetry group for ammonia would prevent any ambiguity when readers implement the Eckart spring.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive and accurate summary of our manuscript, as well as for the recommendation of minor revision. The referee correctly identifies the key advances in extending symmetrized PIMD to rotationally excited states via Eckart projection and the ability to extract multiple J values from a single simulation. Since no specific major comments were raised in the report, we have no points requiring rebuttal or clarification at this time.

Circularity Check

0 steps flagged

Minor self-citation to prior method; new extension derived independently

full rationale

The paper extends its own prior symmetrized PIMD work but derives the rotational projection explicitly from symmetry operations, permutation-inversion-rotation, and Eckart conditions rather than by fitting or redefinition. The claim of numerical exactness within statistical error follows directly from standard PIMD convergence (beads, timestep, spring constants) once parameters are converged, without reducing to a fitted input or self-referential definition. Validation uses independent external benchmarks (variational results for ammonia tunneling splittings, rotational levels for water) and reproduces the experimental J-dependence trend. Self-citation exists but is not load-bearing for the central claims, as the new formalism stands on its own derivations and external checks. No step equates a prediction to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on the standard path-integral formalism plus a domain-specific projection; no new free parameters or invented entities are introduced in the abstract.

axioms (1)

domain assumption The ring-polymer system can be rigorously projected onto selected rotational manifolds and symmetry states via an Eckart spring connecting end beads with a permutation-inversion-rotation operation.
This is the central new device enabling the extension to excited rotational states.

pith-pipeline@v0.9.0 · 5480 in / 1100 out tokens · 31890 ms · 2026-05-10T14:22:20.930934+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages

[1]

Exact tunnelling splittings from symmetrized path integrals,

1G. Trenins, L. Meuser, H. Bertschi, O. Vavourakis, R. Fl¨ utsch, and J. O. Richardson, “Exact tunnelling splittings from symmetrized path integrals,” J. Chem. Phys.159, 034108 (2023)

work page 2023

[1] [1]

Exact tunnelling splittings from symmetrized path integrals,

1G. Trenins, L. Meuser, H. Bertschi, O. Vavourakis, R. Fl¨ utsch, and J. O. Richardson, “Exact tunnelling splittings from symmetrized path integrals,” J. Chem. Phys.159, 034108 (2023)

work page 2023