Networks of Hydrogen Bond Networks in Water Clusters

Water clusters play a prominent role in atmospheric and solution chemistry. Numerous arrangements of protons, H-bond configurations or networks, shape the cluster properties. Studies of small water clusters by cryogenic scanning tunneling microscopy and high-resolution rovibrational spectroscopy have established proton rearrangement mechanisms forming pathways between H-bond networks. The mechanisms, concerted tunneling in particular, describe the local processes connecting pairs of configurations. Here, proton rearrangement networks mapping these transformations are defined and explored to provide a global view of the H-bond configurations in clusters. The networks are constructed for clusters of different sizes and structures. Their analysis reveals an odd-even effect with respect to the number of water molecules, exponential growth of the small-world character, bimodality of the degree distributions, and gapped assortativity of the networks. The last two properties signify the unexpected division of H-bond configurations into two classes according to their network connectivity. The results demonstrate qualitative differences between proton rearrangement mechanisms, suggest a strong influence of the cluster structure. The generated networks are of interest as real-world models for network rewiring; they establish an alternative platform for studies of proton rearrangements in H-bonded 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.

Collective Proton Transfers in Cyclic Water-Ammonia Tetramers: A Path Integral Machine-Learning Study

We present results from machine-learning-based path integral molecular dynamics simulations that describe isomerization paths articulated via collective proton transfers along mixed, cyclic tetramers combining water and ammonia at cryogenic conditions. The net result of such isomerizations is a reverse of the chirality of the global hydrogen-bonding architecture along the different cyclic moieties. In monocomponent tetramers, the classical free energy profiles associated with these isomerizations present the usual symmetric double-well characteristics whereas the reactive paths exhibit full concertedness among the different intermolecular transfer processes. Contrastingly, in mixed water/ammonia tetramers, the incorporation of a second component introduces imbalances in the strengths of the different hydrogen bonds leading to a partial loss of concertedness, most notably at the vicinity of the transition state. As such, the highest and lowest degrees of progression are registered along OH···N and O···HN coordinations, respectively. These characteristics lead to polarized transition state scenarios akin to solvent-separated ion-pair configurations. The explicit incorporation of nuclear quantum effects promotes drastic depletions in the activation free energies and modifications in the overall shape of the profiles which include central plateau-like stages, indicating the prevalence of deep tunneling regimes. On the other hand, the quantum treatment of the nuclei partially restores the degree of concertedness among the evolutions of the individual transfers.

Equilibrium and Dynamical Characteristics of Hydrogen Bond Bifurcations in Water-Water and Water-Ammonia Dimers: A Path Integral Molecular Dynamics Study

We present path integral molecular dynamics results that describe the effects of nuclear quantum fluctuations on equilibrium and dynamical characteristics pertaining to bifurcation pathways in hydrogen bonded dimers combining water and ammonia, at cryogenic temperatures of the order of 20 K. Along these isomerizations, the hydrogen atoms in the molecules acting as hydrogen-bond donors interchange their original dangling/connective characters. Our results reveal that the resulting quantum transition paths comprise three stages: the initial and final ones involve overall rotations during which the two protons retain their classical-like characteristics. Effects from quantum fluctuation are clearly manifested in the changes operated at the intermediate passages over transition states, as the spatial extents of the protons stretch over typical lengths comparable to the distances between connective and dangling basins of attractions. Consequently, the classical over-the-hill path is replaced by a tunneling controlled mechanism which, within the path integral perspective, can be cast in terms of concerted inter-basin migrations of polymer beads from dangling-to-connective and from connective-to-dangling, at practically no energy costs. We also estimated the characteristic timescales describing such interconversions within the approximate ring polymer rate theory. Effects derived from full and partial deuteration are also discussed.

Determination of concerted or stepwise mechanism of hydrogen tunneling from isotope effects: Departure between experiment and theory

Isotope substitution is an important experimental technique that offers deep insight into reaction mechanisms, as the measured kinetic isotope effects (KIEs) can be directly compared with theory. For multiple proton transfer processes, there are two types of mechanisms: stepwise transfer and concerted transfer. The Bell-Limbach model provides a simple theory to determine whether the proton transfer mechanism is stepwise or concerted from KIEs. Recent STM experiments have studied the proton switching process in water tetramers on NaCl(001). Theoretical studies predict that this process occurs via a concerted mechanism, however, the experimental KIEs resemble the Bell-Limbach model for stepwise tunneling, raising question on the underlying mechanism or the validity of the model. We study this system using ab initio instanton theory, and in addition to thermal rates, we also considered microcanonical rates, as well as tunneling splittings. Instanton theory predicts a concerted mechanism, and the KIEs for tunneling rates (both thermal and microcanonical) upon deuteration are consistent with the Bell-Limbach model for concerted tunneling, but could not explain the experiments. For tunneling splittings, partial and full deuteration changes the size of it in a similar fashion to how it changes the rates. We further examined the Bell-Limbach model in another system, porphycene, which has both stepwise and concerted tunneling pathways. The KIEs predicted by instanton theory are again consistent with the Bell-Limbach model. This study highlights differences between KIEs in stepwise and concerted tunneling, and the discrepancy between theory and recent STM experiments. New theory/experiments are desired to settle this problem.

Nanoconfinement Raises the Energy Barrier to Hydrogen Atom Exchange between Water and Glucose

In bulk aqueous environments, the exchange of protons between labile hydroxyl groups typically occurs easily and quickly. Nanoconfinement can dramatically change this normally facile process. Through exchange spectroscopy (EXSY) NMR measurements, we observe that nanoconfinement of glucose and water within AOT (sodium bis(2-ethylhexyl) sulfosuccinate) reverse micelles raises the energy barrier to labile hydrogen exchange, which suggests a disruption of the hydrogen bond network. Near room temperature, we measure barriers high enough to slow the process by as much as 2 orders of magnitude. Although exchange rates slow with decreasing temperatures in these nanoconfined environments, the barrier we measure below ∼285 K is 3-5 times lower than the barrier measured at room temperature, indicating a change in mechanism for the process. These findings suggest the possibility of hydrogen tunneling at a surprisingly high-temperature threshold. Furthermore, differences in exchange rates depend on the hydroxyl group position on the glucose pyranose ring and suggest a net orientation of glucose at the reverse micelle interface.

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.

The Case of Formic Acid on Anatase TiO2 (101): Where is the Acid Proton?

Carboxylic-acid adsorption on anatase TiO₂ is a relevant process in many technological applications. Yet, despite several decades of investigations, the acid-proton localization-either on the molecule or on the surface-is still an open issue. By modeling the adsorption of formic acid on top of anatase(101) surfaces, we highlight the formation of a short strong hydrogen bond. In the 0 K limit, the acid-proton behavior is ruled by quantum delocalization effects in a single potential well, while at ambient conditions, the proton undergoes a rapid classical shuttling in a shallow two-well free-energy profile. This picture, supported by agreement with available experiments, shows that the anatase surface acts like a protecting group for the carboxylic acid functionality. Such a new conceptual insight might help rationalize chemical processes involving carboxylic acids on oxide surfaces.

Cross-term Splittings Due to the Orientational Inequivalence of Proton Magnetic Shielding Tensors: Do Water Molecules Trapped in Crystals Hop or Tunnel?

Water molecules trapped in crystals of barium chlorate monohydrate have been investigated by magic-angle spinning (MAS) proton NMR spectroscopy in the temperature range 110-300 K. At high temperatures, a single spinning sideband pattern is observed. Below 150 K, however, two interleaved spinning sideband manifolds appear, with distinct centerbands that do not coincide with the average isotropic chemical shift seen at high temperatures. This hitherto unknown "cross-term splitting" results from the interplay of the homonuclear dipole-dipole coupling and two anisotropic proton shielding tensors that have identical principal components but nonequivalent orientations. The resulting cross terms cannot be averaged out by rotation about the magic angle. The analysis of the exchange-induced broadening, coalescence, and narrowing of the cross-term splitting in MAS spectra allows one to estimate the rate of exchange of the two protons between 140 and 190 K. The experimental data is compared with ²H and ¹H NMR studies of the same sample reported in the literature. Density functional theory methods are utilized to estimate the thermal activation energy for a 2-fold hopping process of proton exchange about the H-O-H bisector. The Bell-Limbach model allows one to take into account contributions due to incoherent quantum tunneling in the low-temperature regime.

UVC photon-induced denaturing of DNA: A possible dissipative route to Archean enzyme-less replication

Non-equilibrium thermodynamics is a relevant framework from within which to address formidable difficulties encountered in explaining the origin of life; from molecular synthesis and complexation, enzyme-less proliferation, to evolution (including the acquisition of homochirality and information). From within this framework we have proposed that the origin of life was the origin of the dissipative structuring of organic pigments which became the fundamental molecules of life (e.g. RNA and DNA) proliferated through autocatalytic photochemical reactions under the thermodynamic imperative of dissipating the imposed UVC solar photon flux available at the Archean surface. Here we present experimental evidence demonstrating that the absorption and dissipation of UVC light by synthetic DNA of 25 base pairs (and also natural salmon sperm DNA) over a range of temperatures, including below their melting temperature, leads to denaturing. Since denaturing is a non-trivial step on route to enzyme-less replication, our data suggest the possibility of a dissipative route to DNA replication at the origin of life. Such a dissipation-replication relation provides a simple mechanism for the early accumulation of both homochirality and information. Possible mechanisms of UVC photon-induced denaturing of DNA are discussed.

Isotope effect on hydrogen bond symmetrization in hydrogen and deuterium fluoride crystals by molecular dynamics simulation

The isotope effect on the collective proton/deuteron transfer in hydrogen and deuterium fluoride crystals has been investigated at 100 K by ab initio quantum-thermal-bath path-integral molecular dynamics (QTB-PIMD) simulation. The deuterons within a planar zigzag chain of the orthorhombic structure simultaneously flip between covalent and hydrogen bonds due to the barrier crossing through tunnelling. The height of the corresponding static barrier normalized for one deuteron is 29.2 meV. In the HF crystal, all the protons are located at the center of the heavy-atom distance. This evidences the symmetrization of the H-bonds, and indicates that the proton zero-point energy is above the barrier top. The decrease of the heavy-atom distance due to quantum fluctuations in both HF and DF crystals corresponds to a large decrease and an increase of the hydrogen and covalent bond lengths, respectively. Upon deuteration, the increase of the heavy-atom distance (Ubbelohde effect) is in agreement with experimental data.

Unprecedented O:⇔:O compression and H↔H fragilization in Lewis solutions

Charge injection in terms of protons, lone pairs, cations and anions by acid and base solvation mediates the HB network and properties of Lewis solutions through H↔H fragilization, O:⇔:O compression and polarization, ionic polarization and hydrating H₂O dipolar screen shielding, anion–anion repulsion, compressed solvent H–O bond elongation and undercoordinated solute H–O bond contraction.

The collective and quantum nature of proton transfer in the cyclic water tetramer on NaCl(001)

Proton tunneling is an elementary process in the dynamics of hydrogen-bonded systems. Collective tunneling is known to exist for a long time. Atomistic investigations of this mechanism in realistic systems, however, are scarce. Using a combination of ab initio theoretical and high-resolution experimental methods, we investigate the role played by the protons on the chirality switching of a water tetramer on NaCl(001). Our scanning tunneling spectroscopies show that partial deuteration of the H₂O tetramer with only one D₂O leads to a significant suppression of the chirality switching rate at a cryogenic temperature (T), indicating that the chirality switches by tunneling in a concerted manner. Theoretical simulations, in the meantime, support this picture by presenting a much smaller free-energy barrier for the translational collective proton tunneling mode than other chirality switching modes at low T. During this analysis, the virial energy provides a reasonable estimator for the description of the nuclear quantum effects when a traditional thermodynamic integration method cannot be used, which could be employed in future studies of similar problems. Given the high-dimensional nature of realistic systems and the topology of the hydrogen-bonded network, collective proton tunneling may exist more ubiquitously than expected. Systems of this kind can serve as ideal platforms for studies of this mechanism, easily accessible to high-resolution experimental measurements.