2026-04-15 | | Total: 110
We introduce a new equation-of-motion coupled-cluster method based on a pair coupled-cluster doubles (pCCD) reference, termed frozen-pair EOM-CCSD (EOM-fpCCSD). This approach combines the computational efficiency of the pCCD ansatz with a dynamical correlation correction, enabling a reliable description of electronically excited states within the EOM framework. The method has been implemented in the open-source PyBEST software package. Its performance is systematically benchmarked against standard EOM-CCSD and its pair-tailored variant (EOM-ptCCSD), using both canonical Hartree-Fock and pCCD natural orbitals. For charge-transfer (CT) excitations taken from the QUEST database, EOM-fpCCSD yields excitation energies very close to those of EOM-CCSD, outperforming EOM-ptCCSD, as well as to the theoretical best estimates (TBEs). Working within the localized pCCD natural orbital basis allows us to determine the directed CT character, which quantifies the directed charge flow from one molecular domain to another. Numerical results show that EOM-fpCCSD, EOM-CCSD, and EOM-ptCCSD provide nearly identical descriptions of the directed CT character, despite changes in excitation energies. The true advantage of EOM-fpCCSD becomes evident for the challenging QUEST subset of doubly excited states. While EOM-ptCCSD performs similarly to standard EOM-CCSD, EOM-fpCCSD significantly outperforms both methods for these problematic states compared to TBEs. In addition to improving the accuracy of excitation energies, EOM-fpCCSD also converges for several states that standard EOM-CCSD and EOM-ptCCSD fail to converge. These results demonstrate that EOM-fpCCSD offers a promising and computationally efficient route toward a more accurate description of complex electronic excitations.
A screening of oxygen profiles in mid-T treated SRF cavities is crucial, in order to infer physical correlations between the microscopic cavity lattice and cavity performance - a problem concerning acceleration physicists for years. This thesis provides an analysis of oxygen diffusion profiles for three differently treated samples: Two mid-T baked and the third with the standard EuXFEL recipe. The measurement method utilizes EXAFS spectroscopy and was carried out at the DELTA facility in Dortmund. The result suffers heavily under noise, making the quantity of the result barely useable. Qualitatively, no deviations of current models regarding the profiles, could be proven, and no results of previous studies were contradicted. The experimental analysis is described in precision, interpretations of the possible are undertaken and theoretical considerations regarding error estimation and possible EXAFS simulations for future attempts are provided. A repetition at the PETRA III facility is indicated.
This paper examines the polarization-dependent output of various types of 3-mode photonic lanterns fabricated using double-clad fibers. We explore the sensitivity of the modal response across several types of photonic lanterns, from the fully symmetric and strongly coupled structure of regular photonic lanterns to the fully asymmetric structure of mode-selective photonic lanterns. We demonstrate the high sensitivity of the output of photonic lanterns with strong coupling between their ports to the polarization of the input state. In contrast, ports with high isolation or low coupling, such as in mode-selective photonic lanterns, exhibit responses that are almost polarization independent.
Recent advances in data-driven modeling have shown that diffusion models can successfully generate synthetic Lagrangian trajectories in turbulent flows. Building on this progress, we extend the method to the joint generation of pairs of Lagrangian velocity trajectories, enabling a fully data-driven representation of turbulent pair dispersion, a long-standing fundamental problem with broad relevance in fluid dynamics. We demonstrate that diffusion models accurately reproduce the evolution of particle-pair separation, including deviations from Richardson's classical scaling law, while simultaneously preserving all key single-particle statistical properties reported in previous studies. These findings underscore the potential of diffusion-based generative models to emulate high-dimensional, multi-scale turbulent dynamics, further establishing them as a powerful tool for scientific modeling and for future geophysical and astrophysical applications.
We present an efficient implementation of the low-cost linear-response coupled-cluster singles and doubles (LR-CCSD) method for computing static and frequency-dependent polarizabilities in systems with significant relativistic and electron-correlation effects. The approach employs X2C-based Hamiltonians (X2CAMF and X2CMP) and incorporates Cholesky decomposition to reduce memory requirements. In the current implementation, costly three- and four-external index integrals are generated on the fly, eliminating the need for their storage. Benchmark results indicate that the X2CMP Hamiltonian provides more consistent performance than X2CAMF, particularly for large and highly augmented basis sets. The proposed FNS++CD-X2CMP-LR-CCSD method shows excellent agreement with four-component reference values across a wide range of systems. Additionally, different strategies for constructing the FNS++ basis were assessed, and an averaged density approach was found to offer a favorable balance between accuracy and computational cost. On average, about 73% of the virtual spinor space is removed, demonstrating the efficiency and consistency of the FNS++ density-based truncation approach. The present implementation enables accurate and scalable relativistic response calculations for large molecular systems, as demonstrated by the calculation of the static polarizability of the Uranium Hexafluoride complex with a triple-zeta basis set more than 1400 basis functions.
Stimulated Brillouin scattering (SBS), an optical nonlinearity arising from photon-phonon interactions, has formed the basis for a large class of optical signal processing devices, including Brillouin amplifiers. A limiting factor of such amplifiers is the noise due to thermal-mechanical fluctuations that the phonons imprint on the optical signal. Prior work has either inferred or experimentally observed a noise factor ($F$) that depends only on the thermal occupation of the phonons ($F\approx 1+n_{th}$). We show that this noise factor results naturally from a Hamiltonian-based spatio-temporal coupled mode treatment in the limit of large Brillouin amplification and when phonon propagation is neglected. Moreover, this theoretical framework allows us to extend our treatment to a much larger and more representative parameter space for emerging SBS systems; specifically, this analysis accounts for the forward or backward nature of the scattering process and the effects of phonon propagation, optical loss, and small Brillouin gains. Our results demonstrate that the noise factor can deviate radically from $F\approx 1+n_{th}$ for a host of modern SBS devices, especially those in which phonon propagation significantly changes the coupled mode dynamics.
The term 'neutrinoless' is a cornerstone of modern particle physics, yet it defines a fundamental process by what is missing rather than what is created. We trace the origins of this privative neologism to a 1953 experimental claim and show how a 'sociology of suspicion' transformed Ettore Majorana's affirmative ontology into an agnostic shorthand. By examining this linguistic shift, we argue that our current terminology may obscure the profound physical meaning of the search. Reclaiming the language of 'matter creation' is not merely a semantic choice, but a timely conceptual shift to bridge the gap between experimental caution and the radical character of the laws of nature we aim to uncover.
Neuromorphic photonics is emerging as a powerful platform for fast and efficient optical information processing and sensing. However, future brain-inspired photonic systems require compact and scalable light sources, capable of generating the neuro-mimetic optical signals needed for their operation. This work demonstrates a single-stack laser that delivers optical and electrical neural-like spiking emission under solitary operation. Termed the Neuron Surface-Emitting Laser (NeuronSEL), this compact, multi-junction Vertical-Cavity Surface Emitting Laser (VCSEL) exhibits non-linear Negative Differential Resistance (NDR), similar to that observed in memristive devices. Leveraging this NDR behaviour enables the novel demonstration of multiple neuronal features in the NeuronSEL including refractoriness and threshold-/integrate-and-fire dynamics. We demonstrate the NeuronSEL's behaviour as an optical spiking neuron and its ability to perform processing functions, such as coincidence detection and exclusive OR operations. Its scalability is illustrated by proposing a network based on an array of NeuronSELs, able to perform classification tasks. The NeuronSEL emerges as a strong candidate for practical and scalable neuromorphic photonic hardware, with potential impact across a range of applications in optical sensing, communications and computing technologies, whilst benefitting from the inherent advantages of VCSEL technology -low manufacturing cost, compactness, efficiency, vertical emission, and straightforward integration into large arrayed-structures and networks.
Three-dimensional (3D) photonic integrated circuits (PIC) are emerging as an indispensable scheme for high density and multifunctional photonic systems. However, the wafer-scale scaling of PICs towards a 3D configuration is constrained by two key factors: (i) the trade-off between inter-layer taper efficiency and footprint, and (ii) wafer-scale uniformity of inter-layer transition loss. In this work, we introduce etch-back assisted chemical mechanical polishing (E-CMP) to achieve high wafer-scale uniformity of the spacer layer. Moreover, we break the efficiency-footprint trade-off by demonstrating a novel $κ$-engineered taper, achieving a reliability metric that is 75\% higher than the traditional linearly tapered structure. Building on these design and fabrication developments, we enable reliable 3D PICs with typical loss of 0.077 and 0.068 dB/cm on two silicon nitride (SiN) waveguide layers and typical 3D transition loss as low as 6 mdB. Furthermore, the low 3D transition loss enables the first class of 3D high-Q optical cavities occupying two distinct device layers, providing new design space for high-Q optical cavities. The scalable fabrication process and design methodology provide routes for wafer-scale reliable 3D PICs that are promising in a series of applications ranging from photonic interconnects and computing networks to high-density photonic sensors and nonlinear photonics.
There has been a veritable explosion of methods and software to perform machine-learned regression on datasets of electronic energies and forces to develop high-dimensional machine learned potential energy surfaces (ML-PESs). A major, but not deeply-studied aspect is how well different ML-PESs represent the same dataset on which they are trained, beyond the standard fitting precision metrics. Here, this is examined in detail using several ''stress tests'', for two widely applied machine-learned potential approaches. One is based on permutationally invariant polynomial (PIP) linear least square regression and the other is the message-passing neural network PhysNet approach. These potentials and dipole moment surfaces are used in VSCF/VCI calculations of vibrational energies and wavefunctions. The energies from the two PESs are directly compared as are the IR spectra. In addition, tunneling splittings for the hydrogen transfer between two equivalent structures are reported from using three methods: ring polymer instanton theory, diffusion Monte Carlo simulations, and the $Q_{im}$ path method. These calculations require the evaluation of on the order of one billion energies that are widely dispersed in the 15-dimensional configurational space. The two PESs yield results for these quantities in excellent agreement with each other.
Gas hydrates are potential candidates in future energy sources while simultaneously providing structures with extensive applications in carbon capture and storage, gas transport, and important separation processes. Prior research in the field considers the dynamics of the water molecule backbone in particular. We investigated the pressure-enthalpy landscape and mechanical stability envelope of sI methane and carbon dioxide hydrates simulated using DFT. We investigated the effect of the revPBE + DFT-D2 and the SCAN + rVV10 and their treatment of the exchange correlation interactions. We examined the zero pressure material properties, finding that revPBE comparatively underbinds the interactions, causing more flexible structures with large equilibrium volumes. Under pressure, the carbon dioxide molecule was found to align itself parallel to the hexagonal faces of the large cage despite the functional used. Additionally, the property differences are caused by the ability of the carbon dioxide molecule to rotate and disperse the changes in the energy landscape in ways that methane molecules cannot. This computational methodology describes the elastic stability of gas hydrate, marginal stability, and critical differences across important molecular interactions, confirming experimentally observed restrictions in guest molecule rotations and novel pressure behaviors under hydrostatic loads
The integration of active light-emitting elements into planar photonic circuits on a silicon nitride platform remains challenging due to material incompatibilities and high-temperature processing. Proposed hybrid method embeds monodisperse luminescent particles into lithographically defined wells above a 200 nm-thick silicon nitride taper coupler. A fabrication process involving wells etching, particle deposition, and planarization enables precise integration while maintaining waveguide integrity. When pumped at 950 nm the particles emit across 1500--1600 nm, peaking at 1532 nm (FWHM 60 nm), covering the optical telecommunication C-band. Numerical simulations yield an average coupling efficiency of 0.25% into the fundamental waveguide mode, suggesting significant potential for further device optimization. The approach provides a scalable route for integrating broadband telecommunications emitters on a silicon nitride platform.
Neural operators have been increasingly used as data-driven surrogates for time-marching predictions of turbulent flows. However, long-horizon autoregressive prediction is sensitive to error accumulation and the choice of prediction interval. Excessively small time increments may increase temporal redundancy and lengthen rollouts, which can degrade the stability of neural operators in turbulence forecasting. This work pursues a unified objective: stable long-horizon autoregressive prediction at fine temporal resolution for three-dimensional turbulence. We propose a multi-stepsize mixture-of-experts (Ms-MoE) neural operator built on an implicit factorized Transformer (IFactFormer) backbone. The model conditions on a requested relative stride and uses a time-step router to activate scale-specific routed experts together with a shared expert, yielding a single architecture that represents a family of stride-parameterized time-advancement operators. We evaluate the approach on forced homogeneous isotropic turbulence (HIT) and turbulent channel flow using filtered direct numerical simulation datasets. Relative to sampling intervals used in previous studies, we construct training datasets with up to 20 times finer temporal resolution and report long-horizon autoregressive rollouts using qualitative time-slice comparisons and long-time-averaged statistics. Ms-MoE-IFactFormer yields more stable long-horizon rollouts and improved agreement with long-time-averaged statistics on both HIT and turbulent channel flow, suggesting potential for stable time-marching at fine temporal resolution in more complex turbulent flows.
We report on the shape taken by the interface of a liquid bath when hit by a smooth oblique steady jet. When the angle between the jet and the bath decreases below $50^\circ$, a cavity is formed in front of the jet. In the inertial regime we explore, the jet boundary layer detaches in the impact region, thereby delimiting a core jet region outside of which the liquid is mainly in hydrostatic equilibrium. The shape of the outer meniscus is shown to be related to the one outside a tilted fiber piercing the fluid interface. In order to unravel the flow features and separation, we perform direct numerical simulations and show that the flow detachment displays an asymmetry, which results in the acceleration of the liquid below the surface, thereby creating a depression. With this observation, we propose a model balancing the suction force of this depression with the weight of the displaced water and the surface tension force to obtain a prediction for the typical width of the cavity.
Sub-wavelength thick photonic crystal (PhC) slabs coupled to 2D excitonic materials, such as transition metal dichalcogenides (TMDs), are a promising platform for highly tunable, room-temperature, on-chip optoelectronic devices. Unlike conventional Fabry-Perot microcavities, these compact open cavities exhibit non-trivial electric field profiles, leading to spatially distinct regions of weak and strong coupling with excitons within the PhC unit cell. Using coupled mode theory and rigorous solutions to Maxwell's equations, we investigate how the PhC geometry can be used to control these coexisting exciton/polariton contributions and tailor the resulting optical spectra. For large filling factors, i.e., small air gaps, we show that PhC polaritons can be modeled as dark waveguide modes brightened via the periodicity of the PhC slab. Furthermore, by spatially patterning the TMD monolayer based on the local field intensity, we reveal the simultaneous presence of excitons in both the weak and strong coupling regimes. Overall, this work provides fundamental insights into the strong light-matter coupling regime in structured photonic environments, offering a pathway to design and optimize metal-free, ultra-compact polaritonic devices.
Purpose: Accurate dose calculation is essential in radiotherapy for precise tumor irradiation while sparing healthy tissue. With the growing adoption of MRI-guided and real-time adaptive radiotherapy, fast and accurate dose calculation on CT and MRI is increasingly needed. The DoseRAD2026 dataset and challenge provide a public benchmark of paired CT and MRI data with beam-level photon and proton Monte Carlo dose distributions for developing and evaluating advanced dose calculation methods. Acquisition and validation methods: The dataset comprises paired CT and MRI from 115 patients (75 training, 40 testing) treated on an MRI-linac for thoracic or abdominal lesions, derived from the SynthRAD2025 dataset. Pre-processing included deformable image registration, air-cavity correction, and resampling. Ground-truth photon (6 MV) and proton dose distributions were computed using open-source Monte Carlo algorithms, yielding 40,500 photon beams and 81,000 proton beamlets. Data format and usage notes: Data are organized into photon and proton subsets with paired CT-MRI images, beam-level dose distributions, and JSON beam configuration files. Files are provided in compressed MetaImage (.mha) format. The dataset is released under CC BY-NC 4.0, with training data available from April 2026 and the test set withheld until March 2030. Potential applications: The dataset supports benchmarking of fast dose calculation methods, including beam-level dose estimation for photon and proton therapy, MRI-based dose calculation in MRI-guided workflows, and real-time adaptive radiotherapy.
This Perspective is part of a Special Topic that explored the maturity of nonadiabatic molecular dynamics for predicting photochemical processes. In 2023, a prediction challenge was issued to the community of computational photochemists to simulate the photochemistry of cyclobutanone, photoexcited at 200 nm, and the resulting time-resolved MeV-UED signal. The challenge attracted 15 theoretical predictions from more than 70 researchers, employing a wide range of strategies for electronic structure and nonadiabatic molecular dynamics to predict the time-resolved MeV-UED signal before the experiment had been conducted at SLAC (Stanford, USA). The MeV-UED instrument at Shanghai Jiao Tong University was also used to provide a second independent time-resolved MeV-UED signal for the photochemistry of cyclobutanone. This Perspective discusses the various approaches and strategies used by the participants to predict the photochemistry of cyclobutanone. This work also summarizes the strengths and weaknesses of various methods used for photoexcitation, electronic structure, nonadiabatic dynamics, and calculation of observables, as agreed by the participants during a CECAM workshop dedicated to the results of the challenge and organized in Lausanne in April 2025. This Perspective also collects all the predicted time-resolved MeV-UED signals into a single figure, together with the experimental signal. This challenge (i) demonstrated the qualitative predictive power of nonadiabatic molecular dynamics and (ii) underscore the impact of electronic-structure theory on the outcome of the excited-state dynamics and the need for its careful benchmarking. This effort allowed the community to share practical strategies to perform nonadiabatic dynamics (discussed in the present Perspective) and constitutes a 'calibration' exercise for computational photochemistry.
Thin sheets can be assembled into tubular origami structures that combine deployability with pronounced anisotropic stiffness, enabling applications ranging from robotics to deployable systems. However, most existing tubular origami designs remain limited to degree-four vertex topologies and are characterized primarily in axial and radial loading modes, without a full assessment of anisotropic stiffness. Here, we present an automated design framework for tubular origami that jointly explores local vertex topology through generalized degree-$n$ vertices and global tube topology through the polygonal cross-section, for the systematic design and optimization of anisotropic stiffness. Using a calibrated bar-and-hinge model together with experimental validation, we quantify large-deformation stiffness responses in axial translation, in-plane translation, torsion about the tube axis, and rotation about in-plane axes, thereby characterizing the anisotropic stiffness of the tube across its compliant and constrained deformation modes. The resulting design-space exploration showed that the polygonal cross-sectional topology is the primary factor governing the anisotropic stiffness. We further show that increasing the local vertex degree can improve global structural performance, particularly for tubes with a small number of cross-sectional vertices, demonstrating that higher local kinematic freedom does not necessarily compromise stiffness at the structural scale. Compared with a benchmark design, the optimized architectures achieve more than 50 times higher constrained rotational stiffness. Together, these results highlight higher-degree vertices and polygonal cross-sectional topology as powerful design variables for tailoring anisotropic stiffness in tubular origami.
Transferable excited-state dynamics offer a route to efficient screening of photophysical behavior across molecular systems, but conventional nonadiabatic simulations remain prohibitively expensive. Here we introduce X-MACE, a transferable machine-learning potential for excited-state dynamics that predicts multiple potential energy surfaces, forces and oscillator strengths, and combine it with curvature-driven surface hopping to enable data-efficient screening of photochemical pathways. We apply this framework to fluorescent chromophores as an example application, using green fluorescent protein chromophore variants to demonstrate how subtle structural modifications reshape excited-state relaxation, lifetimes and photoisomerization yields. Fine-tuning a single pretrained model with fewer than 100 reference geometries per derivative yields accurate dynamics across a chemically diverse set of analogues. The screening reveals two governing design principles: steric crowding on the phenolate ring lowers the torsional barrier and accelerates access to twisted conical intersections, whereas conjugation extension stabilizes planar excited-state configurations, suppresses non-radiative decay and prolongs fluorescence. More broadly, this workflow provides a general framework for scalable excited-state screening and interpretable design of photophysical properties.
This paper investigates DeGroot-Friedkin (DF) dynamics over signed influence networks with interdependent topics. We propose a multi-topic signed framework that combines repelling interpersonal interactions with cross-issue self-appraisal, examining how antagonism and topic interdependence shape the evolution of agent-level social power. When the logic matrices (for topic interdependence) of all agents share a common dominant left eigenvector, we identify structural conditions under which the original dynamics admit an exact reduction to an explicit scalar DF map. This yields a complete classification of limiting social power configurations into pluralistic, mixed, and vertex-dominant types. In all three cases, the dynamics are globally convergent, and in the first two the ordering induced by the interaction centrality is preserved. We further show local robustness under small heterogeneous perturbations of the logic matrices. We also clarify what changes when this common-eigenvector structure is lost. These results extend signed social power dynamics beyond the standard nonnegative scalar setting and shed light on the robustness and scope of centrality-based social power formation in multi-topic signed influence systems.
We present a quantitative framework for on-axis optical transduction of vertical InP nanowire resonators, correlating laser position to signal amplitude, calibration, and frequency stability. Photothermal resonance detuning is used to reconstruct the local beam intensity profile and to calibrate the photodetector signal using the thermomechanical noise. A noise model incorporating shot noise and spatial variation in substrate reflectance predicts the position-dependent Allan deviation. We find that the optimal detection position lies near the steepest intensity gradient, and that increasing laser power does not significantly improve frequency stability, because the accompanying temperature rise enhances thermomechanical noise and offsets the signal gain. These results establish design guidelines for optimizing nanowire-based sensors in on-axis optical detection schemes.
This work proposes a statistically enhanced framework to address the instability and limited predictive capability of conventional Galerkin-Proper Orthogonal Decomposition (Galerkin-POD) models. The method reformulates the correction of the Galerkin-projected ODE system as a statistical inverse problem, in which the coefficients are inferred through Bayesian inference. By accounting for model uncertainty arising from POD mode truncation and data uncertainty introduced by data noise and numerical postprocessing, the framework systematically updates the ODE system coefficients using an analytical, sampling-free solution based on Gaussian likelihood and inverse-Gamma priors. The approach is first validated using a self-sustained oscillating flow over a dimpled surface at a moderate Reynolds number (Re=3000), demonstrating stable and accurate reproduction of the temporal dynamics and phase trajectories of coherent structures when compared with direct numerical simulation (DNS). It is then applied to a centrifugal compressor featuring strong tip-leakage vortex breakdown and impeller-diffuser interactions at Re=100000, where the model successfully captures dominant unsteady structures and frequency characteristics despite limited mode retention. Overall, the results show that Bayesian inference substantially enhances the robustness, stability, and predictive fidelity of Galerkin-POD models for compressible flow systems. The proposed methodology combines the physical interpretability of Galerkin projection with the statistical rigour of Bayesian inference, offering a general, computationally efficient, and uncertainty-aware reduced-order modelling framework for complex fluid dynamic applications.
Controlling the statistical properties of light, namely the fluctuations in photon arrival, entropy and number, is essential for both classical and quantum photonics. While integrated systems provide tunable control over amplitude, phase, and wavelength, real-time modulation of photon statistics has remained a long-standing challenge. Herein, we introduce the concept and experimental realization of a photon statistics transducer: a high-extinction, broadband electro-optic device capable of deterministically shaping photon-number distributions at nanosecond timescales. Our approach employs a cascaded thin-film lithium niobate (TFLN) Mach-Zehnder amplitude modulator delivering more than 50 dB extinction, enabling precise suppression and release of coherent seed light from an integrated InP laser. By exploiting the interplay between seed suppression and erbium-doped fiber amplifier dynamics, we demonstrate smooth, voltage-controlled switching between Poissonian and super-Poissonian photon statistics, with second-order coherence g2(0) tunable from 1.0 to 1.7. Complementary measurements with superconducting nanowire single-photon detectors further show photon-flux control down to sub-photon levels, highlighting the potential for future operation with non-classical sources. The photon statistics transducer thus establishes statistical modulation as a new functional primitive in integrated photonics. Applications range from entropy generation and secure communication to neuromorphic and hybrid quantum-classical processing, where controlled randomness and entropy are essential resources. By enabling programmable transitions between statistical regimes using only electronic drive signals, our work lays the foundation for adaptive, entropy-aware photonic systems that bridge classical and quantum domains.
A path-integral hybrid Monte Carlo approach with enveloping bridging potentials (PIHMC-EBP) is proposed for calculating numerically exact tunneling splittings in molecular systems. The central idea is to construct an approximately barrierless bridging potential that smoothly connects symmetry-related regions of ring-polymer phase space, enabling direct sampling of the free-energy profile from which the relevant splittings are obtained. Two tailored nonlocal updates are designed to enhance the sampling of slow collective motions. Compared with path-integral molecular dynamics using thermodynamic integration, PIHMC-EBP requires neither quadrature nor time-step convergence checks, thereby substantially reducing the manual effort required to analyze the results. Applications to malonaldehyde (and its deuterated isotopologue) and the HCl dimer using state-of-the-art potential energy surfaces provide the most precise tunneling splittings reported to date for both systems, while simultaneously reducing the overall computational cost by several times and three orders of magnitude, respectively. Finally, application to the water dimer yields the first numerically exact path-integral calculations of the ground-state tunneling splittings on three different potential energy surfaces, all obtained simultaneously by reweighting a single set of trajectories.
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$.