Ion-mediated hydrogen-bond rearrangement through tunnelling in the iodide–dihydrate complex

MBX V1.2: Accelerating Data-Driven Many-Body Molecular Dynamics Simulations

The MBX software provides an advanced platform for molecular dynamics simulations, leveraging state-of-the-art MB-pol and MB-nrg data-driven many-body potential energy functions. Developed over the past decade, these potential energy functions integrate physics-based and machine-learned many-body terms trained on electronic structure data calculated at the "gold standard" coupled-cluster level of theory. Recent advancements in MBX have focused on optimizing its performance, resulting in the release of MBX v1.2. While the inherently many-body nature of MB-pol and MB-nrg ensures high accuracy, it poses computational challenges. MBX v1.2 addresses these challenges with significant performance improvements, including enhanced parallelism that fully harnesses the power of modern multicore CPUs. These advancements enable simulations on nanosecond time scales for condensed-phase systems, significantly expanding the scope of high-accuracy, predictive simulations of complex molecular systems powered by data-driven many-body potential energy functions.

Current Status of the MB-pol Data-Driven Many-Body Potential for Predictive Simulations of Water Across Different Phases

Developing a molecular-level understanding of the properties of water is central to numerous scientific and technological applications. However, accurately modeling water through computer simulations has been a significant challenge due to the complex nature of the hydrogen-bonding network that water molecules form under different thermodynamic conditions. This complexity has led to over five decades of research and many modeling attempts. The introduction of the MB-pol data-driven many-body potential energy function marked a significant advancement toward a universal molecular model capable of predicting the structural, thermodynamic, dynamical, and spectroscopic properties of water across all phases. By integrating physics-based and data-driven (i.e., machine-learned) components, which correctly capture the delicate balance among different many-body interactions, MB-pol achieves chemical and spectroscopic accuracy, enabling realistic molecular simulations of water, from gas-phase clusters to liquid water and ice. In this review, we present a comprehensive overview of the data-driven many-body formalism adopted by MB-pol, highlight the main results and predictions made from computer simulations with MB-pol to date, and discuss the prospects for future extensions to data-driven many-body potentials of generic and reactive molecular systems.

Entropy of Liquid Water as Predicted by the Two-Phase Thermodynamic Model and Data-Driven Many-Body Potentials

We investigate the entropy of liquid water at ambient conditions using the two-phase thermodynamic (2PT) model, applied to both common pairwise-additive water models and the MB-pol and MB-pol(2023) data-driven many-body potentials. Our simulations demonstrate that the 2PT model yields entropy values in semiquantitative agreement with experimental data when using MB-pol and MB-pol(2023). Additionally, our analyses indicate that the entropy values predicted by pairwise-additive water models may benefit from error compensation between the inherent approximations of the 2PT model and the known limitations of these water models in describing many-body interactions. Despite its approximate nature, the simplicity of the 2PT model makes it a valuable tool for estimating relative entropy changes of liquid water across various environments, especially when combined with water models that provide a consistently robust representation of the structural, thermodynamic, and dynamical properties of liquid water.

Nuclear Quantum Effects in Hydroxide Hydrate Along the H-Bond Bifurcation Pathway

Path integral (PI) simulations are used to explore nuclear quantum effects (NQEs) in hydroxide hydrate and its perdeuterated isotopomer along the H-bond bifurcation pathway. Toward this, a new potential energy surface using the symmetric gradient domain machine learning method with ab initio data at the CCSD(T)/aug-cc-pVTZ level is built. From PI umbrella sampling (US) simulations, free energy profiles along the bifurcation coordinate are explored as a function of temperature. At ambient temperature, the bifurcation barrier is increased upon inclusion of NQEs. At low temperatures in the deep tunneling regime, the barrier is strongly decreased and flattened. These trends are examined, and the role of the O-O distance is also investigated through two-dimensional US simulations.

Nuclear quantum effects in gas-phase 2-fluoroethanol

Torsional motions along the FCCO and HOCC dihedrals lead to the five unique conformations of 2-fluoroethanol, of which the conformer that is gauche along both dihedrals has the lowest energy. In this work, we explore how nuclear quantum effects (NQEs) manifest in the structural parameters of the lowest energy conformer, in the intramolecular free energy landscape along the FCCO and HOCC dihedrals, and also in the infrared spectrum of the title molecule, through the use of path integral simulations. We have first developed a full dimensional potential energy surface using the reaction surface Hamiltonian framework. On this potential, we have carried out path integral molecular dynamics simulations at several temperatures starting from the minimum energy well to explore structural influences of NQEs including geometrical markers of the interaction between the OH and F groups. From the computed free energy landscapes, significant reduction of the torsional barrier is found at low temperature near the cis region of the dihedrals, which can be understood through the trends in the radii of gyration of the atomic ring polymers. We find that the inclusion of NQEs in the computation of infrared spectrum is important to obtain good agreement with the experimental band positions.

All-Round Ionic Liquids for Shuttle-Free Zinc-Iodine Battery

The practical implementation of aqueous zinc-iodine batteries (ZIBs) is hindered by the rampant Zn dendrites growth, parasite corrosion, and polyiodide shuttling. In this work, ionic liquid EMIM[OAc] is employed as an all-round solution to mitigate challenges on both the Zn anode and the iodine cathode side. First, the EMIM+ embedded lean-water inner Helmholtz plane (IHP) and inert solvation sheath modulated by OAc- effectively repels H2 O molecules away from the Zn anode surface. The preferential adsorption of EMIM+ on Zn metal facilitates uniform Zn nucleation via a steric hindrance effect. Second, EMIM+ can reduce the polyiodide shuttling by hindering the iodine dissolution and forming an EMIM+ -I3 - dominated phase. These effects holistically enhance the cycle life, which is manifested by both Zn || Zn symmetric cells and Zn-I2 full cells. ZIBs with EAc deliver a capacity decay rate of merely 0.01 ‰ per cycle after over 18,000 cycles at 4 A g-1 , and lower self-discharge and better calendar life than the ZIBs without ionic liquid EAc additive.

Consistent density functional theory-based description of ion hydration through density-corrected many-body representations

Delocalization error constrains the accuracy of density functional theory in describing molecular interactions in ion-water systems. Using Na+ and Cl- in water as model systems, we calculate the effects of delocalization error in the SCAN functional for describing ion-water and water-water interactions in hydrated ions, and demonstrate that density-corrected SCAN (DC-SCAN) predicts n-body and interaction energies with an accuracy approaching coupled cluster theory. The performance of DC-SCAN is size-consistent, maintaining an accurate description of molecular interactions well beyond the first solvation shell. Molecular dynamics simulations at ambient conditions with many-body MB-SCAN(DC) potentials, derived from the many-body expansion, predict the solvation structure of Na+ and Cl- in quantitative agreement with reference data, while simultaneously reproducing the structure of liquid water. Beyond rationalizing the accuracy of density-corrected models of ion hydration, our findings suggest that our unified density-corrected MB formalism holds great promise for efficient DFT-based simulations of condensed-phase systems with chemical accuracy.

Sequential flipping: the donor-acceptor exchange mechanism in water trimers

The donor-acceptor exchange (DAE) is a significant hydrogen bond network rearrangement (HBNR) mechanism because it can lead to the change of the hydrogen bond direction. In this work, we report a new DAE mechanism found in water trimers that is realized by sequential flipping (SF) of all molecules rather than the well-known proton transfer (PT) process. Meanwhile, the SF process has a much smaller potential barrier (0.262 eV) than the previously predicted collective rotation process (about 1.7 eV), implying that the SF process is the main flipping process that can lead to DAE. Importantly, high-precision ab initio calculations show that SF-DAE can make the water ring to show a clear chiral difference from PT-DAE, which brings the prospect of distinguishing the two confusing processes based on circular dichroism spectra. The reaction rate analysis including quantum tunneling indicates an obvious temperature-dependent competitive relationship between the SF and PT processes; specifically, the SF process dominates above 65 K, while the PT process dominates below 65 K. Therefore, in most cases, the contribution for DAE mainly comes from the flipping process, rather than the PT process as previously thought. Our work enriches the understanding of the DAE mechanism in water trimers and provides a piece of the jigsaw that has been sought for the HBNR mechanism.

MBX: A many-body energy and force calculator for data-driven many-body simulations

Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the "many-body energy" (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.

Beyond steric selectivity of ions using ångström-scale capillaries

Ion-selective channels play a key role in physiological processes and are used in many technologies. Although biological channels can efficiently separate same-charge ions with similar hydration shells, it remains a challenge to mimic such exquisite selectivity using artificial solid-state channels. Although there are several nanoporous membranes that show high selectivity with respect to certain ions, the underlying mechanisms are based on the hydrated ion size and/or charge. There is a need to rationalize the design of artificial channels to make them capable of selecting between similar-sized same-charge ions, which, in turn, requires an understanding of why and how such selectivity can occur. Here we study ångström-scale artificial channels made by van der Waals assembly, which are comparable in size with typical ions and carry little residual charge on the channel walls. This allows us to exclude the first-order effects of steric- and Coulomb-based exclusion. We show that the studied two-dimensional ångström-scale capillaries can distinguish between same-charge ions with similar hydrated diameters. The selectivity is attributed to different positions occupied by ions within the layered structure of nanoconfined water, which depend on the ion-core size and differ for anions and cations. The revealed mechanism points at the possibilities of ion separation beyond simple steric sieving.

Stabilization of the Alkylammonium Cations in Halide Perovskite Thin Films by Water-Mediated Proton Transfer

The development of alkylammonium lead trihalide perovskite (ALHP) photovoltaics has grown rapidly over the past decade. However, there are remaining critical challenges, such as proton defects, which can lead to the material instability of ALHPs. Although specific strategies, including the use of halide additives, have significantly reduced the defects, a fundamental understanding of the defect passivation mechanism remains elusive. Herein, an approach and mechanism for minimizing proton defects in ALHP crystals by adding ionized halides to the perovskite precursor solution are reported. This work clarifies that the ionized halides induced proton transfer from H₂ O to the alkylammonium cation in the precursor solution, stabilizing the ALHP crystals. The fundamental characteristics of ALHP and its precursors are examined by X-ray diffraction, transmittance electron microscopy, in situ extended X-ray absorption fine structure, Fourier transform NMR spectroscopy, and Fourier transform infrared spectroscopy. The findings from this work will guide the development of highly stable ALHP crystals, enabling efficient and stable optoelectronic ALHP devices.

Decreased Water Mobility Contributes To Increased α-Synuclein Aggregation

The solvation shell is essential for the folding and function of proteins, but how it contributes to protein misfolding and aggregation has still to be elucidated. We show that the mobility of solvation shell H2 O molecules influences the aggregation rate of the amyloid protein α-synuclein (αSyn), a protein associated with Parkinson's disease. When the mobility of H2 O within the solvation shell is reduced by the presence of NaCl, αSyn aggregation rate increases. Conversely, in the presence CsI the mobility of the solvation shell is increased and αSyn aggregation is reduced. Changing the solvent from H2 O to D2 O leads to increased aggregation rates, indicating a solvent driven effect. We show the increased aggregation rate is not directly due to a change in the structural conformations of αSyn, it is also influenced by a reduction in both the H2 O mobility and αSyn mobility. We propose that reduced mobility of αSyn contributes to increased aggregation by promoting intermolecular interactions.

Unraveling the effect of solvents on the excited state dynamics of C540A by experimental and theoretical study

In this work, the excited-state dynamics including intramolecular charge transfer (ICT) and the redshift of C540A have been investigated in a series of solvents on the basis of the Kamlet-Taft solvatochromic parameters (π*, α, β) using femtosecond transient absorption spectra and systematic theoretical calculation. We demonstrate that the redshift of the emission peak has a linear relationship with the α and π* scales and the effect of the π* scale is slightly stronger than that of the α scale. Meanwhile, the ICT rates can be suggested as relevant to not only the α scale but also the π* scale. Additionally, C540A-AN has proved that the excited state molecules have a unique inactivation mechanism because of the dark feature of the S1 (CT) state. The valuable mechanistic information gleaned from the excited-state dynamics by the experimental and theoretical study would facilitate the design of organic materials for prospective applications in photochemistry and photobiology.

Decreased Water Mobility Contributes To Increased α-Synuclein Aggregation

The solvation shell is essential for the folding and function of proteins, but how it contributes to protein misfolding and aggregation has still to be elucidated. We show that the mobility of solvation shell H2O molecules influences the aggregation rate of the amyloid protein α-synuclein (αSyn), a protein associated with Parkinson's disease. When the mobility of H2O within the solvation shell is reduced by the presence of NaCl, αSyn aggregation rate increases. Conversely, in the presence CsI the mobility of the solvation shell is increased and αSyn aggregation is reduced. Changing the solvent from H2O to D2O leads to increased aggregation rates, indicating a solvent driven effect. We show the increased aggregation rate is not directly due to a change in the structural conformations of αSyn, it is also influenced by a reduction in both the H2O mobility and αSyn mobility. We propose that reduced mobility of αSyn contributes to increased aggregation by promoting intermolecular interactions.

Water Complex of Imidogen

Imidogen (NH) is the simplest nitrogen hydride that plays an important role in combustion and interstellar chemistry, and its combination with H₂O is the prototypical amidation reaction of O-H bonds involving a nitrene intermediate. Herein, we report the observation of the elusive water complex of NH, a prereaction complex associated with the amidation reaction in a solid N₂ matrix at 10 K. The hydrogen-bonded structure of NH···OH₂ (versus HN···HOH) is confirmed via IR spectroscopy with comprehensive isotope labeling (D, ¹⁸O, and ¹⁵N) and quantum chemical calculations at the UCCSD(T)/aug-cc-pVQZ level of theory. In line with the observed absorption at 350 nm, irradiation of the complex at 365 nm leads to O-H bond insertion, yielding hydroxylamine NH₂OH.

Data-Driven Many-Body Potential Energy Functions for Generic Molecules: Linear Alkanes as a Proof-of-Concept Application

We present a generalization of the many-body energy (MB-nrg) theoretical/computational framework that enables the development of data-driven potential energy functions (PEFs) for generic covalently bonded molecules, with arbitrary quantum mechanical accuracy. The "nearsightedness of electronic matter" is exploited to define monomers as "natural building blocks" on the basis of their distinct chemical identity. The energy of generic molecules is then expressed as a sum of individual many-body energies of incrementally larger subsystems. The MB-nrg PEFs represent the low-order n-body energies, with n = 1-4, using permutationally invariant polynomials derived from electronic structure data carried out at an arbitrary quantum mechanical level of theory, while all higher-order n-body terms (n > 4) are represented by a classical many-body polarization term. As a proof-of-concept application of the general MB-nrg framework, we present MB-nrg PEFs for linear alkanes. The MB-nrg PEFs are shown to accurately reproduce reference energies, harmonic frequencies, and potential energy scans of alkanes, independently of their length. Since, by construction, the MB-nrg framework introduced here can be applied to generic covalently bonded molecules, we envision future computer simulations of complex molecular systems using data-driven MB-nrg PEFs, with arbitrary quantum mechanical accuracy.

Heavy-Atom Quantum Tunnelling in Spin Crossovers of Nitrenes

We simulate two recent matrix-isolation experiments at cryogenic temperatures, in which a nitrene undergoes spin crossover from its triplet state to a singlet state via quantum tunnelling. We detail the failure of the commonly applied weak-coupling method (based on a linear approximation of the potentials) in describing these deep-tunnelling reactions. The more rigorous approach of semiclassical golden-rule instanton theory in conjunction with double-hybrid density-functional theory and multireference perturbation theory does, however, provide rate constants and kinetic isotope effects in good agreement with experiment. In addition, these calculations locate the optimal tunnelling pathways, which provide a molecular picture of the reaction mechanism. The reactions involve substantial heavy-atom quantum tunnelling of carbon, nitrogen and oxygen atoms, which unexpectedly even continues to play a role at room temperature.

Toward High-level Machine Learning Potential for Water Based on Quantum Fragmentation and Neural Networks

Accurate and efficient simulation of liquids, such as water and salt solutions, using high-level wave function theories is still a formidable task for computational chemists owing to the high computational costs. In this study, we develop a deep machine learning potential based on fragment-based second-order Møller-Plesset perturbation theory (DP-MP2) for water through neural networks. We show that the DP-MP2 potential predicts the structural, dynamical, and thermodynamic properties of liquid water in better agreement with the experimental data than previous studies based on density functional theory (DFT). The nuclear quantum effects (NQEs) on the properties of liquid water are also examined, which are noticeable in affecting the structural and dynamical properties of liquid water under ambient conditions. This work provides a general framework for quantitative predictions of the properties of condensed-phase systems with the accuracy of high-level wave function theory while achieving significant computational savings compared to ab initio simulations.

The behavior of methane-water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane-water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane-water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane-water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. We find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

Instanton theory for Fermi's golden rule and beyond

Instanton theory provides a semiclassical approximation for computing quantum tunnelling effects in complex molecular systems. It is typically applied to proton-transfer reactions for which the Born-Oppenheimer approximation is valid. However, many processes in physics, chemistry and biology, such as electron transfers, are non-adiabatic and are correctly described instead using Fermi's golden rule. In this work, we discuss how instanton theory can be generalized to treat these reactions in the golden-rule limit. We then extend the theory to treat fourth-order processes such as bridge-mediated electron transfer and apply the method to simulate an electron moving through a model system of three coupled quantum dots. By comparison with benchmark quantum calculations, we demonstrate that the instanton results are much more reliable than alternative approximations based on superexchange-mediated effective coupling or a classical sequential mechanism. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

General Many-Body Framework for Data-Driven Potentials with Arbitrary Quantum Mechanical Accuracy: Water as a Case Study

We present a general framework for the development of data-driven many-body (MB) potential energy functions (MB-QM PEFs) that represent the interactions between small molecules at an arbitrary quantum-mechanical (QM) level of theory. As a demonstration, a family of MB-QM PEFs for water is rigorously derived from density functionals belonging to different rungs across Jacob's ladder of approximations within density functional theory (MB-DFT) and from Møller-Plesset perturbation theory (MB-MP2). Through a systematic analysis of individual MB contributions to the interaction energies of water clusters, we demonstrate that all MB-QM PEFs preserve the same accuracy as the corresponding ab initio calculations, with the exception of those derived from density functionals within the generalized gradient approximation (GGA). The differences between the DFT and MB-DFT results are traced back to density-driven errors that prevent GGA functionals from accurately representing the underlying molecular interactions for different cluster sizes and hydrogen-bonding arrangements. We show that this shortcoming may be overcome, within the MB formalism, by using density-corrected functionals (DC-DFT) that provide a more consistent representation of each individual MB contribution. This is demonstrated through the development of a MB-DFT PEF derived from DC-PBE-D3 data, which more accurately reproduce the corresponding ab initio results.

Data-driven many-body models enable a quantitative description of chloride hydration from clusters to bulk

We present a new data-driven potential energy function (PEF) describing chloride-water interactions, which is developed within the many-body-energy (MB-nrg) theoretical framework. Besides quantitatively reproducing low-order many-body energy contributions, the new MB-nrg PEF is able to correctly predict the interaction energies of small chloride-water clusters calculated at the coupled cluster level of theory. Importantly, classical and quantum molecular dynamics simulations of a single chloride ion in water demonstrate that the new MB-nrg PEF predicts x-ray spectra in close agreement with the experimental results. Comparisons with an popular empirical model and a polarizable PEF emphasize the importance of an accurate representation of short-range many-body effect while demonstrating that pairwise additive representations of chloride-water and water-water interactions are inadequate for correctly representing the hydration structure of chloride in both gas-phase clusters and solution. We believe that the analyses presented in this study provide additional evidence for the accuracy and predictive ability of the MB-nrg PEFs, which can then enable more realistic simulations of ionic aqueous systems in different environments.

Highly Accurate Many-Body Potentials for Simulations of N2O5 in Water: Benchmarks, Development, and Validation

Dinitrogen pentoxide (N₂O₅) is an important intermediate in the atmospheric chemistry of nitrogen oxides. Although there has been much research, the processes that govern the physical interactions between N₂O₅ and water are still not fully understood at a molecular level. Gaining a quantitative insight from computer simulations requires going beyond the accuracy of classical force fields while accessing length scales and time scales that are out of reach for high-level quantum-chemical approaches. To this end, we present the development of MB-nrg many-body potential energy functions for nonreactive simulations of N₂O₅ in water. This MB-nrg model is based on electronic structure calculations at the coupled cluster level of theory and is compatible with the successful MB-pol model for water. It provides a physically correct description of long-range many-body interactions in combination with an explicit representation of up to three-body short-range interactions in terms of multidimensional permutationally invariant polynomials. In order to further investigate the importance of the underlying interactions in the model, a TTM-nrg model was also devised. TTM-nrg is a more simplistic representation that contains only two-body short-range interactions represented through Born-Mayer functions. In this work, an active learning approach was employed to efficiently build representative training sets of monomer, dimer, and trimer structures, and benchmarks are presented to determine the accuracy of our new models in comparison to a range of density functional theory methods. By assessing the binding curves, distortion energies of N₂O₅, and interaction energies in clusters of N₂O₅ and water, we evaluate the importance of two-body and three-body short-range potentials. The results demonstrate that our MB-nrg model has high accuracy with respect to the coupled cluster reference, outperforms current density functional theory models, and thus enables highly accurate simulations of N₂O₅ in aqueous environments.

Instanton calculations of tunneling splittings in chiral molecules

We report the ground state tunneling splittings (ΔE_± ) of a number of axially chiral molecules using the ring-polymer instanton (RPI) method (J. Chem. Phys., 2011, 134, 054109). The list includes isotopomers of hydrogen dichalcogenides H₂ X₂ (X = O, S, Se, Te, and Po), hydrogen thioperoxide HSOH and dichlorodisulfane S₂ Cl₂ . Ab initio electronic-structure calculations have been performed on the level of second-order Møller-Plesset perturbation (MP2) theory either with split-valance basis sets or augmented correlation-consistent basis sets on H, O, S, and Cl atoms. Energy-consistent pseudopotential and corresponding triple zeta basis sets of the Stuttgart group are used on Se, Te, and Po atoms. The results are further improved using single point calculations performed at the coupled cluster level with iterative singles and doubles and perturbative triples amplitudes. When available for comparison, our computed values of ΔE_± are found to lie within the same order of magnitude as values reported in the literature, although RPI also provides predictions for H₂ Po₂ and S₂ Cl₂ , which have not previously been directly calculated. Since RPI is a single-shot method which does not require detailed prior knowledge of the optimal tunneling path, it offers an effective way for estimating the tunneling dynamics of more complex chiral molecules, and especially those with small tunneling splittings.

Microcanonical Tunneling Rates from Density-of-States Instanton Theory

Semiclassical instanton theory is a form of quantum transition-state theory which can be applied to the computation of thermal reaction rates in complex molecular systems including quantum tunneling effects. There have been a number of attempts to extend the theory to treat microcanonical rates. However, the previous formulations are either computationally unfeasible for large systems due to an explicit sum over states or they involve extra approximations, which make them less reliable. We propose a robust and practical microcanonical formulation called density-of-states instanton theory, which avoids the sum over states altogether. In line with the semiclassical approximations inherent to the instanton approach, we employ the stationary-phase approximation to the inverse Laplace transform to obtain the densities of states. This can be evaluated using only post-processing of the data available from a small set of instanton calculations, such that our approach remains computationally efficient. We show that the new formulation predicts results that agree well with quantum scattering theory for an atom-diatom reaction and with experiments for a photoexcited unimolecular hydrogen transfer in a Criegee intermediate. When the thermal rate is evaluated from a Boltzmann average over our new microcanonical formalism, it can overcome some problems of conventional instanton theory. In particular, it predicts a smooth transition at the crossover temperature and is able to describe bimolecular reactions with pre-reactive complexes such as CH₃OH + OH.

A Discrete Chloride Monohydrate: A Solid-State Structural and Spectroscopic Characterization

The solid-state structure of a discrete chloride monohydrate species, [Cl(H₂O)]^-, is reported for the first time. It was isolated as a salt of the tris(dipropylamino)cyclopropenium cation and has been structurally characterized by X-ray and neutron diffraction. Infrared (IR), far-infrared, and Raman spectroscopic studies were also carried out. Additionally, the D₂O and HDO isotopomers were investigated. Of the six fundamental vibrational modes, only the out-of-plane bend ν₃ was not observed as it forms an IR- and Raman-inactive local mode phonon.

Mapping the temperature-dependent and network site-specific onset of spectral diffusion at the surface of a water cluster cage

We explore the kinetic processes that sustain equilibrium in a microscopic, finite system. This is accomplished by monitoring the spontaneous, time-dependent frequency evolution (the frequency autocorrelation) of a single OH oscillator, embedded in a water cluster held in a temperature-controlled ion trap. The measurements are carried out by applying two-color, infrared-infrared photodissociation mass spectrometry to the D₃O⁺·(HDO)(D₂O)₁₉ isotopologue of the "magic number" protonated water cluster, H⁺·(H₂O)₂₁ The OH group can occupy any one of the five spectroscopically distinct sites in the distorted pentagonal dodecahedron cage structure. The OH frequency is observed to evolve over tens of milliseconds in the temperature range (90 to 120 K). Starting at 100 K, large "jumps" are observed between two OH frequencies separated by ∼300 cm^-1, indicating migration of the OH group from the bound OH site at 3,350 cm^-1 to the free position at 3,686 cm^-1 Increasing the temperature to 110 K leads to partial interconversion among many sites. All sites are observed to interconvert at 120 K such that the distribution of the unique OH group among them adopts the form one would expect for a canonical ensemble. The spectral dynamics displayed by the clusters thus offer an unprecedented view into the molecular-level processes that drive spectral diffusion in an extended network of water molecules.

Hydrogen-bonded host-guest systems are stable in ionic liquids

We show that H-bonded host–guest systems associate in ionic liquids (ILs), pure salts with melting point below room temperature, in which dipole–dipole electrostatic interactions should be negligible in comparison with dipole-charge interactions. Binding constants (K_a) obtained from titrations of four H-bonded host–guest systems in two organic solvents and two ionic liquids yield smaller yet comparable K_a values in ionic liquids than in organic solvents. We also detect the association event using force spectroscopy, which confirms that the binding is not solely due to (de)solvation processes. Our results indicate that classic H-bonded host–guest supramolecular chemistry takes place in ILs. This implies that strong H-bonds are only moderately affected by surroundings composed entirely of charges, which can be interpreted as an indication that the balance of Coulombic to covalent forces in strong H-bonds is not tipped towards the former.

Instanton theory of tunneling in molecules with asymmetric isotopic substitutions

We consider quantum tunneling in asymmetric double-well systems for which the local minima in the two wells have the same energy, but the frequencies differ slightly. In a molecular context, this situation can arise if the symmetry is broken by isotopic substitutions. We derive a generalization of instanton theory for these asymmetric systems, leading to a semiclassical expression for the tunneling matrix element and hence the energy-level splitting. We benchmark the method using a set of one- and two-dimensional models, for which the results compare favorably with numerically exact quantum calculations. Using the ring-polymer instanton approach, we apply the method to compute the level splittings in various isotopomers of malonaldehyde in full dimensionality and analyze the relative contributions from the zero-point energy difference and tunneling effects.

Calculations of quantum tunnelling rates for muonium reactions with methane, ethane and propane

Thermal rate constants for Mu + CH₄, Mu + C₂H₆ and Mu + C₃H₈ and their equivalent reactions with H were evaluated with ab initio instanton rate theory. The potential-energy surfaces are fitted using Gaussian process regression to high-level electronic-structure calculations evaluated around the tunnelling pathway. This method was able to successfully reproduce various experimental measurements for the rate constant of these reactions. However, it was not able to reproduce the faster-than-expected rate of Mu + C₃H₈ at 300 K reported by Fleming et al. [Phys. Chem. Chem. Phys., 2015, 17, 19901 and Phys. Chem. Chem. Phys., 2020, 22, 6326]. Analysis of our results indicates that the kinetic isotope effect at this temperature is not significantly influenced by quantum tunnelling. We consider many possible factors for the discrepancy between theory and experiment but conclude that in each case, the instanton approximation is unlikely to be the cause of the error. This is in part based on the good agreement we find between the instanton predictions and new multiconfigurational time-dependent Hartree (MCTDH) calculations for Mu + CH₄ using the same potential-energy surface. Further experiments will therefore be needed to resolve this issue.

Ionization of Water as an Effect of Quantum Delocalization at Aqueous Electrode Interfaces

The enhanced probability of water dissociation at the aqueous electrode interfaces is predicted by path-integral ab initio molecular dynamics. The ionization process is observed at the aqueous platinum interface when nuclear quantum effects are introduced in the statistical sampling, while minor effects have been observed at the gold interface. We characterize the dissociation mechanism of the formed water ions. In spite of the fact that the concentration and lifetime of the ions might be challenging to experimentally detect, they may serve as a guide to future experiments. Our observation might have a significant impact on the understanding of electrochemical processes occurring at the metal electrode surface.

Specific Ion Effects on Hydrogen-Bond Rearrangements in the Halide-Dihydrate Complexes

Small aqueous ionic clusters represent ideal systems to investigate the microscopic hydrogen-bonding structure and dynamics in ion hydration shells. In this context, halide-dihydrate complexes are the smallest systems where the interplay between halide-water and water-water interactions can be studied simultaneously. Here, quantum molecular dynamics simulations unravel specific ion effects on the temperature-dependent structural transition in X^-(H₂O)₂ complexes (X = Cl, Br, and I), which is induced by the breaking of the water-water hydrogen bond. A systematic analysis of the hydrogen-bonding rearrangements at low temperature provides fundamental insights into the competition between halide-water and water-water interactions depending on the properties of the halide ion. While the halide-water hydrogen-bond strength decreases going from Cl^-(H₂O)₂ to I^-(H₂O)₂, the opposite trend in observed in the strength of the water-water hydrogen-bond, suggesting that nontrivial many-body effects may also be at play in the hydration shells of halide ions in solution, especially in frustrated systems (e.g., interfaces) where the water molecules can have dangling OH bonds.