Unusually strong hydrogen bond cooperativity in particular (H2O)20 clusters

Insights into the Effect of Charges on Hydrogen Bonds

Previous computational and experimental studies showed that charges located at the surroundings of hydrogen bonds can exert two opposite effects on them: rupture or strengthening of the hydrogen bond. This work aims to generalize the effect of charges in different hydrogen-bonded systems and to propose a coherent explanation of this effect. For these purposes, 19 systems with intra- and intermolecular hydrogen bonds were studied computationally with DFT. The FT-IR spectra of the systems were simulated, and two energy components of the hydrogen bond were studied separately to determine their variation upon the presence of a charge: charge transfer and molecular overlap. It was determined that either the breaking or strengthening of the hydrogen bond can be favored one over the other, for instance, depending on the heteroatom involved in the hydrogen bond. In addition, it is showed that the strengthening of the hydrogen bond by the presence of a charge is directly related to the decrease in charge transfer between the monomers, which is explained by an increase in molecular overlapping, suggesting a more covalent character of the interaction. The understanding of how hydrogen bonds are affected by charges is important, as it is a key towards a strategy to manipulate hydrogen bonds at convenience.

Pressure-Induced Aggregation of Associating Liquids as a Driving Force Enhancing Hydrogen Bond Cooperativity

The behavior of hydrogen bonds under extreme pressure is still not well understood. Until now, the shift of the stretching vibration band of the X-H group (X = the donor atom) in infrared spectra has been attributed to the variation in the length of the covalent X-H bond. Herein, we combined infrared spectroscopy and X-ray diffraction experimental studies of two H-bonded liquid hexane derivatives, i.e., 2-ethyl-1-hexanol and 2-ethyl-1-hexylamine, in diamond anvil cells at pressures up to the GPa level, with molecular dynamics simulations covering similar thermodynamic conditions. Our findings revealed that the observed changes in the X-H stretching vibration bands under compression are not primarily due to H-bond shortening resulting from increased density but mainly due to cooperative enhancement of H-bonds caused by intensified molecular clustering. This sheds new light on the nature of H-bond interactions and the structure of liquid molecular systems under compression.

Exploring the Origins of the Intensity of the OH Stretch-HOH Bend Combination Band in Water

While the intensity of the OH stretching fundamental transition is strongly correlated to hydrogen-bond strength, the intensity of the corresponding transition to the state with one quantum of excitation in both the OH stretching and HOH bending vibrations in the same water molecule shows a much weaker sensitivity to the hydrogen-bonding environment. The origins of this difference are explored through analyses of the contributions of terms in the expansion of the dipole moment to the calculated intensity. It is found that the leading contribution to the stretch-bend intensity involves the second derivative of the dipole moment with respect to the OH bond length and HOH angle. While this is not surprising, the insensitivity of this derivative to the hydrogen-bonding environment is unexpected. Possible contributions of mode mixing are also explored. While mode mixing leads to splittings of the energies of nearly degenerate excited states, it does not result in significant changes in the sum of the intensities of these transitions. Analysis of changes in the partial charges on the hydrogen atoms upon displacement of the HOH angles shows that these charges generally increase with increasing HOH angle. This effect is partially canceled by a decrease in the charge of the hydrogen atom when a hydrogen bond is broken. The extent of this cancellation increases with the hydrogen bond strength, which is reflected in the observed insensitivity of the intensity of the stretch-bend transition to hydrogen-bond strength.

Enhanced Hydrogen Bonds of the (H2O)n (n = 4-8) Clusters Confined in Uranyl Peroxide Cluster Na20(UO2)20(O2)30

Water is a basic resource and an essential component of living organisms. It often exhibits some novel properties under confinement. The water clusters (H₂O)_n (n = 4-8) confined in the cavity of uranyl peroxide cluster Na₂₀(UO₂)₂₀(O₂)₃₀ (U₂₀) have been computationally investigated by using ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) calculations in this study. The results show that the confined water clusters can form hydrogen bonds with the internal oxygen atoms (O_uranyl) of U₂₀, and their conformations changed significantly. The average lengths (2.553-2.645 Å) of hydrogen bonds in confined (H₂O)_n are shorter than those (2.731-2.841 Å) in the corresponding free water clusters. Moreover, these confined hydrogen bonds show better hydrogen bond patterns according to the quantified indices. The natural bond orbital (NBO) calculations determine that there is electron transferring from the U₂₀ to its interior (H₂O)_n. It is the main reason for enhancing hydrogen bond interactions among the confined water molecules because their oxygen atoms are more negatively charged and their hydrogen atoms are more positively charged. The quantum theory of atoms in molecules (QTAIM) and interacting quantum atoms (IQA) analyses indicate that the confined hydrogen bonds are more covalent, based on the significant electron density ρ(r) and local energy density H(r) at the bond critical points (BCPs), and the stronger energies of interatomic exchange interactions (V_xc). These findings may help to promote the communication of confined water clusters and enrich the understating of confined hydrogen bonds.

Statistical Equivalence of Quantum Chemical Methods for Energy Distribution in Water Clusters

Methods of quantum chemistry are instrumental in understanding molecular structures and properties. However, the results demonstrate significant variability, which is difficult to predict and rationalize. The fundamental question is whether some molecular systems exhibit properties invariant with respect to the computational method. The idea explored here is that collective properties of statistical ensembles should be more robust than characteristics of individual molecules and their arbitrary sets. This effect is demonstrated for the complete set of hydrogen-bond topologies of the dodecahedral water cluster (H₂ O)₂₀ . Non-Gaussian energy distributions produced by various methods have the same functional form despite strong differences in mean values and standard deviations. The conclusion is tested on methods of different complexity and origin employing a number of criteria. A linear mapping between the energies produced by different methods is discussed. The significance of the results is in establishing a collective equivalence property of quantum chemical methods.

Perturbative vibration of the coupled hydrogen-bond (O:H-O) in water

Perturbation Raman spectroscopy has underscored the hydrogen bond (O:H-O or HB) cooperativity and polarizability (HBCP) for water, which offers a proper parameter space for the performance of the HB and electrons in the energy-space-time domains. The OO repulsive coupling drives the O:H-O segmental length and energy to relax cooperatively upon perturbation. Mechanical compression shortens and stiffens the O:H nonbond while lengthens and softens the HO bond associated with polarization. However, electrification by an electric field or charge injection, or molecular undercoordination at a surface, relaxes the O:H-O in a contrasting way to the compression with derivation of the supersolid phase that is viscoelastic, less dense, thermally diffusive, and mechanically and thermally more stable. The HO bond exhibits negative thermal expansivity in the liquid and the ice-I phase while its length responds in proportional to temperature in the quasisolid phase. The O:H-O relaxation modifies the mass densities, phase boundaries, critical temperatures and the polarization endows the slipperiness of ice and superfluidity of water at the nanometer scale. Protons injection by acid solvation creates the H↔H anti-HB and introduction of electron lone pairs derives the O:⇔:O super-HB into the solutions of base or H₂O₂ hydrogen-peroxide. The repulsive H↔H and O:⇔:O interactions lengthen the solvent HO bond while the solute HO bond contracts because its bond order loss. Differential phonon spectroscopy quantifies the abundance, structure order, and stiffness of the bonds transiting from the mode of pristine water to the perturbed states. The HBCP and the perturbative spectroscopy have enabled the dynamic potentials for the relaxing O:H-O bond. Findings not only amplified the power of the Raman spectroscopy but also substantiated the understanding of anomalies of water subjecting to perturbation.

Anisotropy, segmental dynamics and polymorphism of crystalline biogenic carboxylic acids

Carboxylic acids of the Krebs cycle possess invaluable biochemical significance. Still, there are severe gaps in the availability of thermodynamic and crystallographic data, as well as ambiguities prevailing in the literature on the thermodynamic characterization and polymorph ranking. Providing an unambiguous description of the structure, thermodynamics and polymorphism of their neat crystalline phases requires a complex multidisciplinary approach. This work presents results of an extensive investigation of the structural anisotropy of the thermal expansion and local dynamics within these crystals, obtained from a beneficial cooperation of NMR crystallography and ab initio calculations of non-covalent interactions. The observed structural anisotropy and spin-lattice relaxation times are traced to large spatial variations in the strength of molecular interactions in the crystal lattice, especially in the orientation of the hydrogen bonds. A completely resolved crystal structure for oxaloacetic acid is reported for the first time. Thanks to multi-instrumental calorimetric effort, this work clarifies phase behavior, determines third-law entropies of the crystals, and states definitive polymorph ranking for succinic and fumaric acids. These thermodynamic observations are then interpreted in terms of first-principles quasi-harmonic calculations of cohesive properties. A sophisticated model capturing electronic, thermal, and configurational-entropic effects on the crystal structure approaches captures the subtle Gibbs energy differences governing polymorph ranking for succinic and fumaric acids, representing another success story of computational chemistry.

Character of the OH Bend-Stretch Combination Band in the Vibrational Spectra of the "Magic" Number H3O+(H2O)20 and D3O+(D2O)20 Cluster Ions

The fundamental transitions that contribute to the diffuse OH stretching spectrum of water are known to increase in width and intensity with increasing red shift from the free OH frequency. In contrast, the profile of the higher-energy combination band involving the OH stretching and the intramolecular HOH bending modes displays a qualitatively different spectral shape with a much faster falloff on the lower-energy side. We elucidate the molecular origin of this difference by analyzing the shapes of the combination bands in the IR spectra of cryogenically cooled H₃O⁺(H₂O)₂₀ and D₃O⁺(D₂O)₂₀ clusters. The difference in the shapes of the bands is traced to differences in the dependence of their transition dipole matrix elements on the hydrogen-bonding environment. The fact that individual transitions across the combination band envelope have similar intensities makes it a useful way to determine the participation of various sites in extended H-bonding networks.

Correlated Local Fluctuations in the Hydrogen Bond Network of Liquid Water

The hypothesis that liquid water can separate into two phases in the supercooled state has been supported by recent experimental and theoretical studies. However, whether such structural inhomogeneity extends to ambient conditions is under intense debate. Due to the dynamic nature of the hydrogen bond network of liquid water, exploring its structure requires detailed insight into the collective motion of neighboring water molecules, a missing link that has not been examined so far. Here, highly sensitive quantum mechanical calculations detect that the time evolution of nearby hydrogen bonds is strongly correlated, revealing a direct mechanism for the appearance of short-range structural fluctuations in the hydrogen bond network of liquid water for the first time. This correlated dynamics is found to be closely connected to the static structural picture. The distortions from the tetrahedral structure do not occur independently but are correlated due to the preference of nearby donors and acceptors to be in similar environments. The existence of such cooperative fluctuations is further supported by the temperature dependence of the local structural evolution and explained by conventional analysis of localized orbitals. It was found that such correlated structural fluctuations are only observed on a short length scale in simulations at ambient conditions. The correlations of the nearby hydrogen bond pairs of liquid water unveiled here are expected to offer a new insight into connecting the dynamics of individual water molecules and the local structure of the hydrogen bond network.

Decay of hydrogen bonding in mixtures of aliphatic heptanols and bistriflimide ionic liquids

Hydrogen bonding in liquids of the constitution isomers of heptan-1-ol mixed with 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic liquids (ILs), [C_nmim][NTf₂], is investigated using both computational and experimental techniques. All-atom non-polarizable molecular-dynamics (MD) simulations predict that the hydrogen bonds gradually decay with increasing temperature. This decay is more pronounced for the branched alcohols and in the presence of the ionic liquids. The primary and linear isomer, heptan-1-ol, and its tertiary and bulky analogue 3-ethylpentan-3-ol are identified as the opposite extremes of the spectrum of hydrogen bonding stability in the bulk liquid. While neat heptan-1-ol exhibits strong hydrogen bonding at 350 K, 3-ethylpentan-3-ol is prone to hydrogen bonding decay already at 300 K. The presence of ionic liquids is found to affect the hydrogen bonding comparably as a 50 K temperature increase. Since the heat capacities of the associating liquids are very sensitive to any variation in hydrogen bonding strength and to the character of the hydrogen-bonded clusters in the bulk liquid, the calorimetric effort provides useful experimental data to confirm the results predicted by MD simulations. In this work, excess heat capacity is measured for equimolar single-phase mixtures of alcohols and ILs, and it differs largely in its sign and magnitude for individual heptanol isomers. Temperature trends of the excess heat capacities suggest that the stability of hydrogen bonding for individual heptanol isomers is temperature-shifted, based on their capability of hydrogen bonding. The predicted hierarchy of hydrogen bonding in individual alcohols and its impact on the excess heat capacity trends are qualitatively confirmed via thermodynamic modelling of the associative contribution to the excess heat capacities. These terms are found to predetermine the observed non-monotonous excess heat capacity trends.

Advances in the study of supercooled water

In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.

Cooperativity and Anticooperativity in Ion-Water Interactions: Implications for the Aqueous Solvation of Ions

Non-additive effects in hydrogen bonds (HB) take place as a consequence of electronic charge transfers. Therefore, it is natural to expect cooperativity and anticooperativity in ion-water interactions. Nevertheless, investigations on this matter are scarce. This paper addresses the interactions of (i) the cations Li⁺ , Na⁺ , K⁺ , Be²⁺ , Mg²⁺ , and Ca²⁺ together with (ii) the anions F^- , Cl^- , Br^- , NO₃ ^- and SO₄ ^2- with water clusters (H₂ O)_n , n=1-8, and the effects of these ions on the HBs within the complete molecular adducts. We used quantum chemical topology tools, specifically the quantum theory of atoms in molecules and the interacting quantum atoms energy partition to investigate non-additive effects among the interactions studied herein. Our results show a decrease on the interaction energy between ions and the first neighbouring water molecules with an increment of the coordination number. We also found strong cooperative effects in the interplay between HBs and ion-dipole interactions within the studied systems. Such cooperativity affects considerably the interactions among ions with their first and second solvation shells in aqueous environments. Overall, we believe this article provides valuable information about how ion-dipole contacts interact with each other and how they relate to other interactions, such as HBs, in the framework of non-additive effects in aqueous media.

DFT study of phenol alkylation with propylene on H-BEA in the absence and presence of water

Water reduces the activation barrier of the rate-limiting step of phenol alkylation with propylene in H-BEA. This, in turn, increases the transition-state theory rate coefficient by two orders-of-magnitude, suggesting much faster alkylation.

Isolating the Contributions of Specific Network Sites to the Diffuse Vibrational Spectrum of Interfacial Water with Isotopomer-Selective Spectroscopy of Cold Clusters

Decoding the structural information contained in the interfacial vibrational spectrum of water requires understanding how the spectral signatures of individual water molecules respond to their local hydrogen bonding environments. In this study, we isolated the contributions for the five classes of sites that differ according to the number of donor (D) and acceptor (A) hydrogen bonds that characterize each site. These patterns were measured by exploiting the unique properties of the water cluster cage structures formed in the gas phase upon hydration of a series of cations M⁺·(H₂O)_n (M = Li, Na, Cs, NH₄, CH₃NH₃, H₃O, and n = 5, 20-22). This selection of ions was chosen to systematically express the A, AD, AAD, ADD, and AADD hydrogen bonding motifs. The spectral signatures of each site were measured using two-color, IR-IR isotopomer-selective photofragmentation vibrational spectroscopy of the cryogenically cooled, mass selected cluster ions in which a single intact H₂O is introduced without isotopic scrambling, an important advantage afforded by the cluster regime. The resulting patterns provide an unprecedented picture of the intrinsic line shapes and spectral complexities associated with excitation of the individual OH groups, as well as the correlation between the frequencies of the two OH groups on the same water molecule, as a function of network site. The properties of the surrounding water network that govern this frequency map are evaluated by dissecting electronic structure calculations that explore how changes in the nearby network structures, both within and beyond the first hydration shell, affect the local frequency of an OH oscillator. The qualitative trends are recovered with a simple model that correlates the OH frequency with the network-modulated local electron density in the center of the OH bond.