2026-11-09 | | Total: 7
Quantum-state conversion between the ortho and para nuclear-spin states of H$_2$ is forbidden in the isolated molecule but these states can become connected through interactions with a surrounding environment. How that same environment governs the dynamics of the newly connected states remains poorly understood. Here, we investigate this problem using H$_2$ confined in crystalline molecular solids, where the host environment selectively enables conversion between a pair of otherwise disconnected ortho-H$_2$ and para-H$_2$ states. Rotationally resolved infrared spectroscopy directly determines the state-specific conversion dynamics in crystalline CO$_2$ and N$_2$O at 10 K, revealing that increasing H$_2$ concentration suppresses nuclear-spin conversion by nearly an order of magnitude, reaching complete suppression by 20% H$_2$ loading. Molecular dynamics and Monte Carlo simulations identify two complementary origins of this concentration dependence: changes in the thermally populated phonon modes available to mediate relaxation, and strain-driven H$_2$ clustering that progressively removes molecules from the host crystal field responsible for enabling conversion. The resulting picture connects a classical collective response of the host to microscopic quantum-state dynamics: H$_2$ loading generates strain that reorganizes the local molecular environment, thereby modifying both the phonon bath and the molecular coupling to that bath. These results establish a multiscale mechanism by which condensed-phase environments can govern the connectivity and dynamics of an effective two-level system.
Simulating molecules in quantum electrodynamic (QED) cavities requires a consistent description of electronic response, nuclear motion, and electromagnetic-field fluctuations. Here, I present Cavity-QED Molecular Dynamics, an interface connecting eT and i-PI for atomistic simulations of finite molecular systems coupled to a single photon mode. The code evaluates a self-consistent cavity Born-Oppenheimer Hartree-Fock (CBO-HF) energy surface and supports classical-nuclei dynamics under canonical (NVT) and microcanonical (NVE) ensembles, as well as canonical joint nuclear-photon path-integral sampling including nuclear quantum effects (NQEs). A fixed-displacement gradient correction provides nuclear and photon forces consistent with the same energy surface as the photon coordinate evolves independently. OpenMP threading, concurrent bead evaluations, and analysis utilities support practical simulations. Water-monomer benchmarks assess both dynamical propagation and equilibrium sampling. Energy-derivative checks and NVE timestep tests demonstrate consistent force evaluation and integration, while path integral molecular dynamics (PIMD) simulations with PILE-L thermostatting reproduce analytic finite-bead photon second moments, including their centroid and internal contributions, within sampling uncertainties. At the tested parameters, cavity coupling produces a small but statistically resolved increase in photon fluctuations, while changes in mean molecular geometry remain within sampling uncertainty. The resulting Cavity-QED Molecular Dynamics framework enables controlled investigations of how coupled quantum fluctuations of matter and the electromagnetic field influence molecular equilibrium properties within the CBO-HF approximation, providing a foundation for future studies of cavity-modified hydrogen bonding, vibrational spectra, and molecular polaritonic phenomena.
The hierarchical equations of motion (HEOM) are a powerful, numerically exact approach for simulating the time evolution of open quantum systems. Recently, several tensor-train- and tensor-network-based approaches have been suggested and realized to extend the range of applicability of the HEOM approach to larger and more realistic systems. In this paper, we build upon the existing twin-space formulation of the HEOM and present a novel approach to propagate the fermionic HEOM efficiently using the multi-configurational time-dependent Hartree method (MCTDH). We benchmark the resulting HEOM+MCTDH method with electron transport calculations for a Holstein-type nanojunction model. Furthermore, we demonstrate the capabilities of the HEOM+MCTDH method by applying it to a more complex vibronic coupling nanojunction model with strong molecule-electrode coupling for which fully quantum simulations have not been feasible previously.
Within classical nucleation theory, crystal nucleation requires sufficient supersaturation, but is suppressed at low temperatures due to reduced molecular mobility, leading to a maximum nucleation rate typically above the glass transition temperature $T_g$. Here, we experimentally show that this conventional picture breaks down in a biocompatible polymer poly(2-methoxyethyl acrylate)(PMEA)/water mixture. We find that the number of nuclei increases dramatically below $T_g$, with the effective nucleation rate exhibiting an unusual maximum at approximately $T_g-10$ K. Calorimetric experiments reveal a strong temperature-history dependence of cold crystallization of water in hydrated PMEA, which is explained by the observed nucleation behavior. We propose that the unusual enhancement of nucleation below $T_g$ provides a key mechanism linking high glass-forming ability and cold crystallization upon heating, thereby explaining the striking asymmetry in crystallization behavior between cooling and heating. Our findings suggest that a hidden nucleation-rate maximum below $T_g$ may also underlie cold crystallization in other glass-forming liquids.
Few molecular processes are more fundamental to chemistry than the autoionization of water, whose equilibrium constant, $K_w \approx 10^{-14}$, underlies acid-base equilibria and the chemistry of life. Yet the molecular origins of the free energy cost of forming and separating ions remain incompletely understood. Efforts to explain it have been limited by the difficulty of modeling water reactivity in sufficient detail. Here we combine machine learning interatomic potentials (MLIPs) with spectral adaptive biasing force enhanced sampling to connect proton transfer and water wire configurations to the underlying free energy landscape. We identify competing influences of long-range electrostatics, system size, temperature, and nuclear quantum motion on the predicted ionization thermodynamics. The free energy surface extends from neutral water to ions separated beyond the solvent-separated ion pair (SSIP). It reveals a transition state in which a proton is shared between the second and third molecules of a three-membered water wire. The contact ion pair is unstable and the SSIP only weakly metastable, so determining pK$_w$ requires sampling longer wires. Long-range electrostatics raise the free energy at ion separations beyond $3\,\mathrm{\mathring{A}}$. Models with comparable force errors give substantially different free energy profiles in this region. Increasing system size lowers the calculated pK$_w$, while nuclear quantum effects bring it closer to experiment and reproduce the higher pK$_w$ of D$_2$O relative to H$_2$O. These results establish a baseline for interpreting thermodynamic results from MLIP simulations in which charged and polar groups react and diffuse in finite cells. Our analysis also explains how cancellation among competing effects can yield ionization constants close to experiment even when simulations omit several physical contributions.
Infrared spectroscopy is a key tool to study complex structural dynamics in molecular systems. However, resolving weak or broad vibrational modes is challenging due to their rapid decay. Here, we demonstrate a technique to significantly extend the detectable lifetime - and thereby spectral resolution - of molecular signatures by employing synchronized time-gated optical parametric amplification of mid-infrared free-induction decay signals. Utilizing a Cr:ZnS laser system and a nonlinear ZnGeP2 crystal, we achieve a 10-fold increase in the electric field strength of the infrared molecular signal, yielding a two-orders of magnitude enhancement in dynamic range and a 2-fold increase in the signal-to-noise ratio. Using isopropyl alcohol as a model system, we demonstrate that OPA-enhanced detection extends the observable window of vibrational modes up to 72 %. This approach allows for the reconstruction of high-fidelity molecular fingerprints that resolve complex vibrational relaxation pathways and weak structural dynamics, previously hidden within the background and noise, providing a pathway for new high-sensitivity and time-resolved studies of complex chemical structures.
Significant challenges exist in the nonlinear extensional rheology of entangled polymers. With simulations, we show that the key to understanding this problem is to recognize the existence and importance of a strain-induced crossover from the entangled state to a packed state. This crossover, following the saturation of primitive chain stretch, takes place with massive release of entanglements via convective flow and progressive chain alignment. After the crossover, the disentangled, fully-aligned chain segments pack similarly to the random packing of rods. Meanwhile, the system enters the steady state. In this state, the tube model, built on the concept of entanglement, fails, while the stress can be quantitatively calculated by combining the intra-chain conformational contribution and a frictional contribution from the inter-chain separation along flow.