2026-07-08 | | Total: 12
Roaming reactions are organized not by potential-energy saddles but by transition states that are unstable invariant objects in phase space, periodic orbits in the two degrees of freedom studied here. To ask what controls roaming, we modify the Chesnavich model of a barrierless ion--molecule dissociation: its orientation-dependent angular hindrance is replaced by a transverse-stiffness ridge whose angular frequency peaks at an interior radius, and the classical dynamics are studied at a fixed energy just above the dissociation threshold. Comparing two ensembles that differ only in this angular interaction (same radial channel, energy, and inward initial conditions) isolates its effect. The ridge gates entry into the inner well, cutting inner capture from $57\%$ to $15\%$ and returning most of the incoming flux directly to reactants; it does not eliminate roaming but relocates it outward, suppressing it inside the ridge and switching it on farther out. The model retains analogues of the original model's three transition states (tight, free-rotor, and outer orbiting orbits), which we locate as unstable periodic orbits. The tight orbit spans a dividing surface that coincides, within numerical accuracy, with the variational minimum-flux surface, and it carries no barrier along the reaction coordinate: a deep entropic bottleneck placed at an interior radius by the stiffness maximum. Its entropic character is shared with the original model. Strength-matched monotone controls show that the gating tracks the hindrance strength at the bottleneck radius; what the interior maximum supplies is placement, concentrating that strength where it gates most effectively. The trajectories it admits roam nonstatistically, with nonexponential gap-time distributions: the entropic bottleneck governs how much is captured, not the dynamics that follow.
The Hessian of the ring polymer spring potential in the standard Trotter path integral is a $P\times P$ symmetric circulant matrix with a centroid eigenvalue of zero. All such matrices commute and are diagonalised by the same bead to normal mode transformation matrix, and their eigenvalues contain $\lceil P/2\rceil-1$ degenerate pairs by symmetry. However, this still leaves some freedom to improve on the Trotter approximation: one can optimise the remaining $\lfloor P/2\rfloor$ independent non-zero normal mode frequencies to fit the exact quantum mechanical radii of gyration of harmonic ring polymers with frequencies in the range $0\leω\leω_{\rm max}$, where $ω_{\rm max}$ is the maximum physical frequency in the problem of interest. The optimisation involves solving a simple least squares problem for the optimum (economised or "Eco") internal mode frequencies. The remainder of the calculation then proceeds in the same way as a Trotter path integral calculation. An example application to hexagonal ice shows that the convergence of the Eco path integral is comparable to that of the 4th order Suzuki-Chin path integral, but with purely 2nd order Trotter effort. There is no need to calculate the projected Hessians that arise in the Suzuki-Chin method by finite differences, there is no need to develop any new estimators for observables, and once the Eco frequencies have been calculated the implementation of the Eco path integral involves changing just a few lines of a Trotter path integral code. To provide a more impressive example we have implemented the Eco method in GPUMD and used it to converge the (negative) thermal expansion coefficient and the constant pressure heat capacity of MOF-5 with a machine-learned neuroevolution potential.
Characterization of the decomposition products of energized Criegee intermediates is essential for assessing their impact on the chemical evolution of the atmosphere. Here, a generic and microscopically resolved approach is used to determine the molecular fragmentation pathways and products for CF$_3$CHOO. They include, among others, direct formation of CO$_2$ + CHF$_3$ (HFC-23), HF + CO$_2$ + CF$_2$, and fragmentation routes that are not evident from static reaction path calculations alone. The computed probability for formation of HFC-23 of 14 \% qualitatively agrees with a value of $(7.9^{+0.4}_{-0.2})$ \% from recent measurements, given the differences in the two approaches. Non-statistical dynamics is found for almost all decomposition pathways and the simulations show that excess energy can redirect reaction outcomes away from minimum-energy pathways. The results highlight the power of machine-learned PESs to elucidate multi-step reaction mechanisms of atmospherically relevant intermediates beyond traditional Master equation/electronic structure approaches to provide molecular-level understanding of the role of dynamics.
We present a multiconfigurational binary-encounter (MC-BE) framework for direct analytical evaluation of electron-impact electronic-excitation cross sections for dipole-allowed transitions from ab initio excited-state data. The method combines the threshold-modified Mott-Massey (TMMM) approximation with binary-encounter (BE/BE$f$) scaling. The effective binding energy in the BE/BE$f$ prefactor is obtained from amplitude-weighted occupied-orbital contributions computed by linear-response time-dependent density functional theory (LR-TDDFT), without system-specific fitting parameters. For benzene, MC-BE/TMMM cross sections for the dominant $1{}^{1}\!E_{\mathrm{1u}}$ ($π\!\to\!π^{\ast}$) band agree well with experiment at incident energies $T=10$-$20$ eV and improve on the Schwinger multichannel/truncated configuration-interaction singles results of Falkowski et al. [J. Chem. Phys. 159, 194301 (2023)] for this band and energy range. For naphthalene, the calculated total excitation cross section reproduces the onset and principal maximum of the gas-phase apparent fluorescence cross section, used as an emission-based proxy under dipole-dominated conditions, without empirical energy shifts or intensity scaling. Analytic peak-position and peak-height expressions, parameterized by $r=\langle B\rangle/ΔE$, show that typical valence excitations peak at incident energies $T\simeq 1.5$-$1.6ΔE$ with substantial BE attenuation, providing a diagnostic for relating measured cross-section profiles to excitation energies. Although demonstrated with LR-TDDFT, the framework is transferable to other excited-state theories that provide compatible amplitudes and well-defined orbital energies. These results support MC-BE/TMMM as a practical, inexpensive route for modeling electron-impact excitation of polyatomic molecules with finite oscillator strength.
We present a unified, memory-bounded GPU realization of equation-of-motion coupled-cluster (EOM-CC) excited-state gradients and interstate nonadiabatic couplings (NACMEs) on a single 8\,GB consumer GPU. Both are built from one contraction directed acyclic graph: the EOM-CC relaxation is the reverse-mode transpose of the forward density build rather than a per-state re-derivation, and an atomic-orbital-direct Laplace-transform $J/K$ kernel, made non-symmetric ($J^x(A,B)\neq J^x(B,A)$) by the transition densities, resolves every energy denominator with no four-index molecular-orbital tensor; a two-sided Davidson returns both eigenvectors from one device-resident, spin-pure solve. The pipeline is \emph{validated end to end at small scale}: gradients and NACMEs match finite differences across four spin multiplicities and full configuration interaction to $<\!10^{-12}$ for two electrons, and the excited-state gradient matches the independent \textsc{Psi4} code to $\le\!4.6\times10^{-7}~E_h/a_0$ from \ce{H2O} to aromatic benzene. The kernels and the ground-state solve reach chromophores ($\le\!730$ AO) in 8\,GB, and a frozen-natural-virtual compression lets the eigensolver \emph{execute} a complete excited-state gradient and $Q$--$B$ NACME of the chlorophyll-core chromophore \ce{Mg}-porphine (def2-SVP, $439$ AO) on the card. We present that run as a \emph{capability demonstration} -- executed and translationally invariant to machine zero, but anchored only piece-wise and bounded by a direct convergence study at ${\sim}10^{-2}~E_h/a_0$ -- not a converged spectroscopic result. The validated small-scale capability and the memory-bounded implementation are the contribution.
Spherical colloids with harsh repulsive forces have long been used as experimental analogs of the hard sphere model, with demonstrated good agreement with computer simulations for bulk and structural properties of the fluid, glass and crystal phases. However, an enigmatic discrepancy remains for the crystal-melt stiffness coefficient. Here we perform computer simulations of colloidal hard spheres under tunable buoyant mass and show that the long-standing discrepancy can be traced to a hitherto unrecognized gravity dependent contribution of the stiffness coefficient. This effect is one practical realization of a more general result for the external field dependence of stiffness coefficients of arbitrary interfaces.
We present an efficient implementation of the Bethe--Salpeter equation (BSE) based on numerical atomic orbitals (NAOs) and norm-conserving pseudopotentials within the ABACUS+LibRPA framework. By exploiting the localized resolution-of-identity (LRI) technique, the screened Coulomb interaction is cast into a real-space, unit-cell-indexed form $W_{μν}(\boldsymbol R)$ that is inherently short-ranged and well localized. This spatial locality enables an efficient Fourier interpolation of the BSE kernel from the coarse $\boldsymbol k$-mesh used in the preceding $GW$ calculation to an arbitrarily dense $\boldsymbol k$-mesh on which the BSE Hamiltonian is assembled and diagonalized, thereby giving rise naturally to a dual-$\boldsymbol k$-mesh workflow. Building on this scheme, we systematically examine the convergence of the absorption spectra with respect to the NAO basis set, the auxiliary basis set, and the $\boldsymbol k$-point sampling. Benchmark calculations for both molecular and periodic systems collectively validate the accuracy of the present implementation and establish the dual-$\boldsymbol k$-mesh strategy as a practical and reliable approach for $GW$+BSE calculations.
We introduce population-transfer dynamics as a practical validation observable for active-space-derived reduced Hamiltonians in multistate reaction-center chemistry. Using a cytochrome P450-inspired Fe-complex model, we construct a reaction-coordinate-dependent effective Hamiltonian from state-averaged complete active-space self-consistent field (SA-CASSCF) calculations, map it to a quantum-circuit representation suitable for current hardware, and propagate dynamics from the reactant-side ground state. The reduced Hamiltonian reproduces the SA-CASSCF reference with an RMS deviation of 0.030 eV and a maximum absolute deviation of 0.143 eV. As a dynamics-based diagnostic, the product-manifold population p_P(t) identifies a pronounced near-degeneracy region around x = 0.3, where state mixing is strongest. Classical exact time evolution yields a product population of 0.488 at x = 0.3 after 10 fs, compared with 7.26 x 10^-2 at x = 0.2 and 5.90 x 10^-3 at x = 0.0. To enable execution on current trapped-ion hardware, we examine the trade-off between dynamical fidelity and circuit resources through coupling pruning and first-order Trotterization. A coupling cutoff of 0.02 eV reduces the non-zero coupling set from 32 to 7 while preserving the dominant transfer pathways, and M = 30 provides the best practical operating point. Finally, we demonstrate the workflow on Quantinuum's trapped-ion quantum computer Reimei. The hardware reproduces the key reaction-coordinate trend identified by the classical model, including the maximum at x = 0.3, where the measured product population is 0.42 on hardware and 0.43 on the matched emulator. This work establishes a dynamics-based framework for assessing active-space-derived reduced Hamiltonians and demonstrates chemically interpretable multistate electronic dynamics on current trapped-ion hardware.
Rapidly scaling autonomous science is limited not only by algorithms, compute or data volume, but by which physical records a platform exposes before action. We formulate physically accessible decision-making (PADM) and a measurement-access risk frontier: the Bayes-optimal target risk minimized over records realizable under cost, bandwidth, latency, disturbance, memory and actuation constraints. The frontier gives a no-free-autonomy limit: automation cannot collapse decision uncertainty by computation alone; an optimal controller cannot remove target components absent from its record, and closing that gap requires expanded access, auditing, tolerated disturbance, slower operation or restricted deployment. In monitored feedback, displacement-only control remains exposed to a hidden switching force, whereas a finite-bandwidth cue recovers part of the missing projection before action. A chemistry-aware candidate-ranking audit with a 1000-target stress panel, Gaussian sensing, hidden-regime decisions and cost-aware/thermodynamic channel selection provide reproducible checks. PADM identifies target-specific audit value and residual oracle gaps before deployment.
More than 80 years ago Kramers published a paper calculating how fast a Brownian particle escapes from a potential well over an activation barrier. Since then Kramers' model has been widely adopted by nuclear physics, biophysics and chemical physics communities as a description of activated barrier crossing. From a chemical kinetics perspective, Kramers' theory provides a mapping from continuous dynamics to discrete-state chemical kinetics. Motivated by recent developments, this Perspective provides a rigorous way of performing such a mapping, explaining why and how Kramers' theory works from several points of view. Specifically, we consider transitions of a Brownian particle between two potential wells corresponding to the ``reactant'' and the ``product'' of a chemical reaction. A central unifying idea is to divide the equilibrium ensemble of possible states of the system into two sub-ensembles corresponding to the reactant and product states and then to consider fluxes between these sub-ensembles. Importantly, naive separation based on the location measured relative to the barrier top does not result in a mapping that is physically tenable, and instead the past of the trajectory should be considered. Thus constructed reactant and product ensembles provide an internally consistent description of the problem when also viewed from two different perspectives: one based on the definition of the rate as a conditional transition probability per unit time and the other based on the relaxation modes of the time-evolution operator governing the dynamics.
Separating the contributions of homogeneous dephasing from inhomogeneous broadening in spectral linewidths is essential for connecting optical spectra to microscopic dissipation and disorder mechanisms. Voigt fits to one-dimensional spectra, such as photoluminescence, yield strongly correlated Gaussian and Lorentzian widths, so that neither width can be determined independently with confidence. We quantitatively show that two-dimensional coherent spectroscopy (2DCS) reduces this degeneracy by providing orthogonal spectral slices, diagonal and cross-diagonal, with complementary sensitivity to homogeneous and inhomogeneous broadening processes. As a demonstration, we measure the exciton resonance in hBN-encapsulated MoSe$_2$ at 8 K. A joint uncertainty-weighted fit maps the full $χ^2(σ,γ)$ landscape to quantify parameter covariance. Compared with linear Voigt analysis, 2DCS yields more compact confidence regions and markedly reduced parameter correlation, enabling reliable separation of the contributions to the excitonic linewidth.
Foundation machine learning force fields (MLFFs) such as MACE-MP-0 and UMA cover broad chemical space at near density functional theory (DFT) accuracy. However, they assume equilibrium ground-state physics and do not natively handle externally induced changes to the electronic state, such as charging, applied fields, or electronic excitation, which limits their use for driven processes such as photoexcitation and charge injection. We propose EquiFiLM, a lightweight extension that adds continuous external conditioning to any equivariant foundation MLFF via a per-layer Feature-wise Linear Modulation (FiLM) block, learning externally driven changes to the potential energy surface from minimal training data. The block modulates only scalar channels and preserves E(3)-equivariance exactly. We demonstrate the recipe on charged liquid water with the foundation model MACE-MatPES as the backbone, yielding E-MACE. On the four training charges, E-MACE delivers a $3.1\times$ reduction in force RMSE ($21.3$ to $6.96$ meV/$\mathring{A}$) and a $61\times$ reduction in per-atom energy RMSE ($6.1$ to $0.1$ meV/atom) over a baseline without EquiFiLM trained on the same data, at indistinguishable inference cost. Across seven held-out interpolation and extrapolation charges, force RMSE stays within $18-61$ meV/$\mathring{A}$ and energy RMSE within $0.7-5.4$ meV/atom. The model runs stable molecular dynamics across the full range tested and predicts the charge-dependent first-shell response of the reduced pair distribution function probed by ultrafast electron diffraction. Adding this conditioning axis to the foundation requires only a few thousand DFT-labeled frames, against the $\approx 10^8$ structures of a charge-aware foundation trained from scratch. The recipe is backbone- and conditioning-agnostic: it applies without architectural change to any equivariant MLFF with scalar interaction-layer channels.