Nuclear quantum effects on the structure and the dynamics of [H2O]8 at low temperatures

Impact of Nuclear Quantum Effects on the Structural Properties of Protonated Water Clusters

Nuclear quantum effects (NQEs) play a crucial role in hydrogen-bonded systems due to quantum tunneling and proton fluctuation. Our understanding of how NQEs affect microstructures mainly focuses on bulk phases of liquids and solids but remains deficient for water clusters, including their hydrogen nuclei, hydrogen-bonded configurations, and temperature dependence. Here, we conducted ab initio molecular dynamics (MD) and path integral MD simulations to investigate the influence of NQEs on the structural properties of protonated water clusters H+(H2O)n (n = 3, 6, 9, 12). The results reveal that the NQEs become less evident as the cluster size increases due to the competition between NQEs and electrostatic interactions. Simulations of several H+(H2O)6 isomers at different temperatures indicate that the effect of elevated temperature on proton transfer is related to the initial structure. Interestingly, the process of proton transfer also involves the interconversion between Zundel-type and Eigen-type isomers. These findings significantly deepen our understanding of ion-water and water-water interactions, opening new avenues for the study of hydrated ion clusters and related systems.

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.

Structure and dynamics of nanoconfined water and aqueous solutions

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Nuclear quantum effects on the hydrogen bond donor-acceptor exchange in water-water and water-methanol dimers

We present results from path integral molecular dynamics simulations that describe effects from the explicit incorporation of nuclear quantum fluctuations on the topology of the free energy associated with the geared exchange of hydrogen bonds in the water-water dimer. Compared to the classical treatment, our results reveal important reductions in the free energy barriers and changes at a qualitative level in the overall profile. Most notable are those manifested by a plateau behavior, ascribed to nuclear tunneling, which bridges reactant and product states, contrasting with the usual symmetric double-well profile. The characteristics of the proton localizations along the pathway are examined. An imaginary time analysis of the rotational degrees of freedom of the partners in the dimer at the vicinities of transition states shows a clear "anticorrelation" between intermolecular interactions coupling beads localized in connective and dangling basins of attractions. As such, the transfer is operated by gradual concerted inter-basin migrations in opposite directions, at practically no energy costs. Modifications operated by partial deuteration and by the asymmetries in the hydrogen bonding characteristics prevailing in water-methanol heterodimers are also examined.

Combined Effects from Solvation and Nuclear Quantum Fluctuations on Autoionization Mechanisms in Aqueous Clusters

Using path-integral molecular dynamics simulations, we examine isomerization paths involving collective proton transfers in [H₂O]₅ and [H₂O]₈ clusters under cryogenic conditions. We focused attention on combined effects derived from solvation and nuclear quantum fluctuations on the characteristics of free energy barriers and relative stabilities of reactants and products. In particular, we analyzed two different processes: the first one involves the exchange of donor-acceptor hydrogen bond roles along cyclic moieties, whereas the second one corresponds to charge separation leading to stable [H₃O]⁺[OH]^- ion pairs. In the first case, the explicit incorporation of quantum tunneling introduces important modifications in the classical free energy profile. The resulting quantum profile presents two main contributions, one corresponding to compressions of O-O distances and a second one ascribed to nuclear tunneling of the light protons. Solvation effects promote a moderate polarization of the cyclic structures and a partial loss of concertedness in the collective modes, most notably, at the onset of tunneling. Still, the latter effects are also sufficiently strong to promote the stabilization of ion pairs along the classical trajectories. In contrast, the explicit incorporation of nuclear quantum fluctuations leads to charge separated configurations that are marginally stable. As such, the latter states could also be regarded as short-lived intermediate states along the reactive exchange path.

Inclusion of nuclear quantum effects for simulations of nonlinear spectroscopy

The computation and interpretation of nonlinear vibrational spectroscopy is of vital importance for understanding a wide range of dynamical processes in molecular systems. Here, we introduce an approach to evaluate multi-time response functions in terms of multi-time double symmetrized Kubo transformed thermal correlation functions. Furthermore, we introduce a multi-time extension of ring polymer molecular dynamics to evaluate these Kubo transforms. Benchmark calculations show that the approximations are useful for short times even for nonlinear operators, providing a consistent improvement over classical simulations of multi-time correlation functions. The introduced methodology thus provides a practical way of including nuclear quantum effects in multi-time response functions of non-linear optical spectroscopy.

Surface Isotope Segregation as a Probe of Temperature in Water Nanoclusters

Using ring polymer molecular dynamics simulations, we examine equilibrium and dynamical characteristics of solid-like, aqueous clusters that combine isotopic mixtures of HDO dilute in H2O, at temperatures intermediate between 50 and 175 K. In particular, we focus attention on the relative thermodynamic stabilities of the two isotopes at dangling hydrogen bond sites. The water octamer is analyzed as a reference system. For this aggregate, decreasing temperature yields a gradual stabilization of the light isotope at dangling sites in molecules acting as single-donor-double-acceptors of hydrogen bonds. At T ∼ 50 K, the imbalance between the corresponding quantum kinetic energies leads to a free energy difference between dangling and hydrogen bonded sites of the order of ∼2kBT. Similar free energy differences were found at dangling sites in Nw = 50 water clusters. The extent of the H/D segregation can be adequately monitored by modifications in the peak intensity of the high frequency shoulder of the stretching band of the infrared spectrum. These signals, in turn, represent a potential experimental signature of the elusive temperature of clusters in molecular beams.