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.

Robust Binding of Disulfide-Substituted Rhenium Bipyridyl Complexes for CO2 Reduction on Gold Electrodes.

Heterogenization of homogenous catalysts on electrode surfaces provides a valuable approach for characterization of catalytic processes in operando conditions using surface selective spectroelectrochemistry methods. Ligand design plays a central role in the attachment mode and the resulting functionality of the heterogenized catalyst as determined by the orientation of the catalyst relative to the surface and the nature of specific interactions that modulate the redox properties under the heterogeneous electrode conditions. Here, we introduce new [Re(L)(CO)3Cl] catalysts for CO2 reduction with sulfur-based anchoring groups on a bipyridyl ligand, where L = 3,3'-disulfide-2,2'-bipyridine (SSbpy) and 3,3'-thio-2,2'-bipyridine (Sbpy). Spectroscopic and electrochemical analysis complemented by computational modeling at the density functional theory level identify the complex [Re(SSbpy)(CO)3Cl] as a multi-electron acceptor that combines the redox properties of both the rhenium tricarbonyl core and the disulfide functional group on the bipyridyl ligand. The first reduction at -0.85 V (vs. SCE) involves a two-electron process that breaks the disulfide bond, activating it for surface attachment. The heterogenized complex exhibits robust anchoring on gold surfaces, as probed by vibrational sum-frequency generation (SFG) spectroscopy. The binding configuration is normal to the surface, exposing the active site to the CO2 substrate in solution. The attachment mode is thus particularly suitable for electrocatalytic CO2 reduction.

Isotopic equilibria in aqueous clusters at low temperatures: Insights from the MB-pol many-body potential.

By combining path-integrals molecular dynamics simulations with the accurate MB-pol potential energy surface, we investigate the role of alternative potential models on isotopic fractionation ratios between H and D atoms at dangling positions in water clusters at low temperatures. Our results show clear stabilizations of the lighter isotope at dangling sites, characterized by free energy differences ΔG that become comparable to or larger than kBT for temperatures below ∼75 K. The comparison between these results to those previously reported using the empirical q-TIP4P/F water model [P. E. Videla et al., J. Phys. Chem. Lett. 5, 2375 (2014)] reveals that the latter Hamiltonian overestimates the H stabilization by ∼25%. Moreover, predictions from the MB-pol model are in much better agreement with measured results reported for similar isotope equilibria at ice surfaces. The dissection of the quantum kinetic energies into orthogonal directions shows that the dominant differences between the two models are to be found in the anharmonic characteristics of the potential energy surfaces along OH bond directions involved in hydrogen bonds.

Isotope effects in aqueous solvation of simple halides.

We present a path-integral-molecular-dynamics study of the thermodynamic stabilities of DOH⋯ X- and HOD⋯ X- (X = F, Cl, Br, I) coordination in aqueous solutions at ambient conditions. In agreement with experimental evidence, our results for the F- case reveal a clear stabilization of the latter motif, whereas, in the rest of the halogen series, the former articulation prevails. The DOH⋯ X- preference becomes more marked the larger the size of the ionic solute. A physical interpretation of these tendencies is provided in terms of an analysis of the global quantum kinetic energies of the light atoms and their geometrical decomposition. The stabilization of the alternative ionic coordination geometries is the result of a delicate balance arising from quantum spatial dispersions along parallel and perpendicular directions with respect to the relevant O-H⋯X- axis, as the strength of the water-halide H-bond varies. This interpretation is corroborated by a complementary analysis performed on the different spectroscopic signals of the corresponding IR spectra.

Positional Isotope Exchange in HX·(H2O)n (X = F, I) Clusters at Low Temperatures.

We present molecular dynamics simulation results describing proton/deuteron exchange equilibria along hydrogen bonds at the vicinity of HX acids (X = F, I) in aqueous clusters at low temperatures. To allow for an adequate description of proton transfer processes, our simulation scheme resorted on the implementation of a multistate empirical valence bond hamiltonian coupled to a path integral scheme to account for effects derived from nuclear quantum fluctuations. We focused attention on clusters comprising a number of water molecules close to the threshold values necessary to stabilize contact-ion-pairs. For X = F, our results reveal a clear propensity of the heavy isotope to lie at the bond bridging the halide to the nearest water molecule. Contrasting, for X = I, the thermodynamic stability is reversed and the former connectivity is preferentially articulated via the light isotope. These trends remain valid for undissociated and ionic descriptions of the stable valence bond states. The preferences are rationalized in terms of differences in the quantum kinetic energies of the isotopes which, in turn, reflect the extent of the local spatial confinements prevailing along the different hydrogen bonds in the clusters. In most cases, these features are also clearly reflected in the characteristics of the corresponding stretching bands of the simulated infrared spectra. This opens interesting possibilities to gauge the extent of the isotopic thermodynamic stabilizations and the strengths of the different hydrogen bonds by following the magnitudes and shifts of the spectral signals in temperature-controlled experiments, performed on mixed clusters combining H2O and HOD.

Communication: Isotopic effects on tunneling motions in the water trimer.

We present results of ring polymer molecular dynamics simulations that shed light on the effects of nuclear quantum fluctuations on tunneling motions in cyclic [H2O]3 and [D2O]3, at the representative temperature of T = 75 K. In particular, we focus attention on free energies associated with two key isomerization processes: The first one corresponds to flipping transitions of dangling OH bonds, between up and down positions with respect to the O-O-O plane of the cluster; the second involves the interchange between connecting and dangling hydrogen bond character of the H-atoms in a tagged water molecule. Zero point energy and tunneling effects lead to sensible reductions of the free energy barriers. Due to the lighter nature of the H nuclei, these modifications are more marked in [H2O]3 than in [D2O]3. Estimates of the characteristic time scales describing the flipping transitions are consistent with those predicted based on standard transition-state-approximation arguments.

Isotopic Preferential Solvation of I(-) in Low-Temperature Water Nanoclusters.

We present results from ring polymer molecular dynamics experiments that provide microscopic insights into the characteristics of the isotopic stabilizations of H and D aqueous species in the first solvation shell of a halide I(-) ion in water nanoclusters at low temperatures. The analysis of the simplest I(-)·(HOD) dimer shows a clear propensity for the light isotope to lie at the non-hydrogen-bonded dangling position. Our results predict that, at T ∼ 50 K, I(-)·(DOH) isomers are three times more abundant than I(-)·(HOD) ones. The reasons for such stabilization can be traced back to differences in the nuclear kinetic energy projected along directions perpendicular to the plane of the water molecule. Dynamical implications of these imbalances are shown to be reflected in the characteristics of the corresponding bands of the infrared spectroscopic signals. A similar analysis performed in larger aggregates containing ∼20 water molecules reveals, in contrast, a stabilization of the light isotope along I(-)···HO hydrogen bonds. Effects derived from the consideration of smaller halide anions with larger electric fields at the surface are also examined.

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.

A quantum molecular dynamics study of aqueous solvation dynamics.

Ring polymer molecular dynamics experiments have been carried out to examine effects derived from nuclear quantum fluctuations at ambient conditions on equilibrium and non-equilibrium dynamical characteristics of charge solvation by a popular simple, rigid, water model, SPC/E , and for a more recent, and flexible, q-TIP4P/F model, to examine the generality of conclusions. In particular, we have recorded the relaxation of the solvent energy gap following instantaneous, ±e charge jumps in an initially uncharged Lennard-Jones-like solute. In both charge cases, quantum effects are reflected in sharper decays at the initial stages of the relaxation, which produce up to a ∼20% reduction in the characteristic timescales describing the solvation processes. For anionic solvation, the magnitude of polarization fluctuations controlling the extent of the water proton localization in the first solvation shell is somewhat more marked than for cations, bringing the quantum solvation process closer to the classical case. Effects on the solvation response from the explicit incorporation of flexibility in the water Hamiltonian are also examined. Predictions from linear response theories for the overall relaxation profile and for the corresponding characteristic timescales are reasonably accurate for the solvation of cations, whereas we find that they are much less satisfactory for the anionic case.

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

We use ring-polymer-molecular-dynamics (RPMD) techniques and the semi-empirical q-TIP4P/F water model to investigate the relationship between hydrogen bond connectivity and the characteristics of nuclear position fluctuations, including explicit incorporation of quantum effects, for the energetically low lying isomers of the prototype cluster [H2O]8 at T = 50 K and at 150 K. Our results reveal that tunneling and zero-point energy effects lead to sensible increments in the magnitudes of the fluctuations of intra and intermolecular distances. The degree of proton spatial delocalization is found to map logically with the hydrogen-bond connectivity pattern of the cluster. Dangling hydrogen bonds exhibit the largest extent of spatial delocalization and participate in shorter intramolecular O-H bonds. Combined effects from quantum and polarization fluctuations on the resulting individual dipole moments are also examined. From the dynamical side, we analyze the characteristics of the infrared absorption spectrum. The incorporation of nuclear quantum fluctuations promotes red shifts and sensible broadening relative to the classical profile, bringing the simulation results in much more satisfactory agreement with direct experimental information in the mid and high frequency range of the stretching band. While RPMD predictions overestimate the peak position of the low frequency shoulder, the overall agreement with that reported using an accurate, parameterized, many-body potential is reasonable, and far superior to that one obtains by implementing a partially adiabatic centroid molecular dynamics approach. Quantum effects on the collective dynamics, as reported by instantaneous normal modes, are also discussed.

Aqueous electrolytes confined within functionalized silica nanopores.

Molecular dynamics simulations have been carried out to investigate structural and dynamical characteristics of NaCl aqueous solutions confined within silica nanopores in contact with a "bulk-like" reservoir. Two types of pores, with diameters intermediate between 20 Å and 37.5 Å, were investigated: The first one corresponded to hydrophobic cavities, in which the prevailing wall-solution interactions were of the Lennard-Jones type. In addition, we also examined the behavior of solutions trapped within hydrophilic cavities, in which a set of unsaturated O-sites at the wall were transformed in polar silanol Si-OH groups. In all cases, the overall concentrations of the trapped electrolytes exhibited important reductions that, in the case of the narrowest pores, attained 50% of the bulk value. Local concentrations within the pores also showed important fluctuations. In hydrophobic cavities, the close vicinity of the pore wall was coated exclusively by the solvent, whereas in hydrophilic pores, selective adsorption of Na(+) ions was also observed. Mass and charge transport were also investigated. Individual diffusion coefficients did not present large modifications from what is perceived in the bulk; contrasting, the electrical conductivity exhibited important reductions. The qualitative differences are rationalized in terms of simple geometrical considerations.