On the strength of hydrogen bonding within water clusters on the coordination limit

Thermodynamic and optical properties of HCOOH(H2O)n and HCOOH(NH3)(H2O)(n-1) clusters at various temperatures and pressures: a computational study

Density functional theory has been used to compute the gas-phase geometries, binding energies, ZPE-corrected binding energies, BSSE-corrected binding energies, binding enthalpies, and binding free energies of HCOOH(H₂O)_n and HCOOH(NH₃)(H₂O)_(n-1) clusters with n = 1-8, 10, 12, 14, 16, 18, and 20. Enthalpies and free energies are calculated for a range of atmospherically relevant temperatures (T) and pressures (P) (from T = 298.15 K, P = 1013.25 hPa to T = 216.65 K, P = 226.32 hPa). The optical properties of those clusters have been studied at the CAM-B3LYP/aug-cc-pVDZ level of theory.

The IQA Energy Partition in a Drug Design Setting: A Hepatitis C Virus RNA-Dependent RNA Polymerase (NS5B) Case Study

The interaction of the thumb site II of the NS5B protein of hepatitis C virus and a pair of drug candidates was studied using a topological energy decomposition method called interacting quantum atoms (IQA). The atomic energies were then processed by the relative energy gradient (REG) method, which extracts chemical insight by computation based on minimal assumptions. REG reveals the most important IQA energy contributions, by atom and energy type (electrostatics, sterics, and exchange–correlation), that are responsible for the behaviour of the whole system, systematically from a short-range ligand–pocket interaction until a distance of approximately 22 Å. The degree of covalency in various key interatomic interactions can be quantified. No exchange–correlation contribution is responsible for the changes in the energy profile of both pocket–ligand systems investigated in the ligand–pocket distances equal to or greater than that of the global minimum. Regarding the hydrogen bonds in the system, a “neighbour effect” was observed thanks to the REG method, which states that a carbon atom would rather not have its covalent neighbour oxygen form a hydrogen bond. The combination of IQA and REG enables the automatic identification of the pharmacophore in the ligands. The coarser Interacting Quantum Fragments (IQF) enables the determination of which amino acids of the pocket contribute most to the binding and the type of energy of said binding. This work is an example of the contribution topological energy decomposition methods can make to fragment-based drug design.

A QCT View of the Interplay between Hydrogen Bonds and Aromaticity in Small CHON Derivatives

The somewhat elusive concept of aromaticity plays an undeniable role in the chemical narrative, often being considered the principal cause of the unusual properties and stability exhibited by certain π skeletons. More recently, the concept of aromaticity has also been utilised to explain the modulation of the strength of non-covalent interactions (NCIs), such as hydrogen bonding (HB), paving the way towards the in silico prediction and design of tailor-made interacting systems. In this work, we try to shed light on this area by exploiting real space techniques, such as the Quantum Theory of Atoms in Molecules (QTAIM), the Interacting Quantum Atoms (IQA) approaches along with the electron delocalisation indicators Aromatic Fluctuation (FLU) and Multicenter (MCI) indices. The QTAIM and IQA methods have been proven capable of providing an unbiased and rigorous picture of NCIs in a wide variety of scenarios, whereas the FLU and MCI descriptors have been successfully exploited in the study of diverse aromatic and antiaromatic systems. We used a collection of simple archetypal examples of aromatic, non-aromatic and antiaromatic moieties within organic molecules to examine the changes in π delocalisation and aromaticity induced by the Aromaticity and Antiaromaticity Modulated Hydrogen Bonds (AMHB). We observed fundamental differences in the behaviour of systems containing the HB acceptor within and outside the ring, e.g., a destabilisation of the rings in the former as opposed to a stabilisation of the latter upon the formation of the corresponding molecular clusters. The results of this work provide a physically sound basis to rationalise the strengthening and weakening of AMHBs with respect to suitable non-cyclic non-aromatic references. We also found significant differences in the chemical bonding scenarios of aromatic and antiaromatic systems in the formation of AMHB. Altogether, our investigation provide novel, valuable insights about the complex mutual influence between hydrogen bonds and π systems.

Partition of the electronic energy of the PM7 method via the interacting quantum atoms approach

Partitions of the electronic energy such as that provided by the Interacting Quantum Atoms (IQA) approach have given valuable insights for numerous chemical systems and processes. Unfortunately, this kind of analysis may involve the integration of scalar fields over very irregular volumes, a condition which leads to a large and often prohibitive computational effort. These circumstances have limited the use of these energy partitions to systems comprising a few tens of atoms at most. On the other hand, semiempirical methods have proved useful in the study of systems of several thousands of atoms. Therefore, the goal of this work is to carry out partitions of the semiempirical method PM7 in compliance with the IQA approach. For this purpose, we computed one- and two-atomic energetic contributions whose sum equals the PM7 electronic energy. We illustrate how one might exploit the partition of electronic energies computed via the PM7 method by considering small organic and inorganic molecules and the energetics of individual hydrogen bond interactions within several water clusters which include (H₂O)₃₀, (H₂O)₅₀ and (H₂O)₁₀₀. We also considered the solvation of the amphiphilic caprylate anion to exemplify how to exploit the energy partition proposed in this paper. Overall, this investigation shows how the approach put forward herein might give further insights of the interactions occurring within complex systems in physical and biological chemistry.

High-Density “Windowpane” Coordination Patterns of Water Clusters and Their NBO/NRT Characterization

Cluster mixture models for liquid water at higher pressures suggest the need for water clusters of higher coordination and density than those commonly based on tetrahedral H-bonding motifs. We show here how proton-ordered water clusters of increased coordination and density can assemble from a starting cyclic tetramer or twisted bicyclic (Möbius-like) heptamer to form extended Aufbau sequences of stable two-, three-, and four-coordinate “windowpane” motifs. Such windowpane clusters exhibit sharply reduced (~90°) bond angles that differ appreciably from the tetrahedral angles of idealized crystalline ice I_h. Computed free energy and natural resonance theory (NRT) bond orders provide quantitative descriptors for the relative stabilities of clusters and strengths of individual coordinative linkages. The unity and consistency of NRT description is demonstrated to extend from familiar supra-integer bonds of the molecular regime to the near-zero bond orders of the weakest linkages in the present H-bond clusters. Our results serve to confirm that H-bonding exemplifies resonance–covalent (fractional) bonding in the sub-integer range and to further discount the dichotomous conceptions of “electrostatics” for intermolecular bonding vs. “covalency” for intramolecular bonding that still pervade much of freshman-level pedagogy and force-field methodology.

Implementation of the interacting quantum atom energy decomposition using the CASPT2 method

We present an implementation of the interacting quantum atom (IQA) energy decomposition scheme using the complete active space second-order perturbation theory (CASPT2). This combination yields a real-space interpretation tool with a proper account of the static and dynamic correlation that is particularly relevant for the description of processes in electronic excited states. The IQA/CASPT2 approach allows determination of the energy redistribution that takes place along a photophysical/photochemical deactivation path in terms of self- and interatomic contributions. The applicability of the method is illustrated by the description of representative processes spanning different bonding regimes: noble gas excimer and exciplex formation, the reaction of ozone with a chlorine atom, and the photodissociations of formaldehyde and cyclobutane. These examples show the versatility of using CASPT2 with the significant information provided by the IQA partition to describe chemical processes with a large multiconfigurational character.

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.

Hydration of selenolate moiety: Ab initio investigation of properties of O-H⋯Se(-) hydrogen bonds in CH3 Se(-)⋯(H2 O)n clusters

This work is devoted to investigations of the influence of O-H···Se(-) hydrogen bonds on the electronic shells of selenolate R-Se(-) fragment (R═CH₃ ). The geometric, energetic and nuclear magnetic resonance (NMR) spectral parameters for various conformers of CH₃ Se(-)⋯(H₂ O)_n clusters with n = 0-6 were calculated at CCSD/aug-cc-pVDZ level of theory. For selenolate anion CH₃ Se(-) solvation free energy was calculated, and for water media it is equal to -71.41 kcal/mol. For O-H···Se(-) hydrogen bonds the proportionality coefficients between QTAIM parameters at (3; -1) bond critical point and the strength of an individual hydrogen bond ∆E were proposed. It was shown, that despite a relative weakness of O-H···Se(-) hydrogen bonds, the outer electronic shell of the selenium atom changes significantly upon formation of each hydrogen bond. This, in turn, cause the dramatic change of NMR parameters of selenium nuclei.

Bonding analysis of water clusters using quasi-atomic orbitals

The quasi-atomic orbital (QUAO) bonding analysis introduced by Ruedenberg and co-workers is used to develop an understanding of the hydrogen bonds in small water clusters, from the dimer through the hexamer (bag, boat, book, cyclic, prism and cage conformers). Using kinetic bond orders as a metric, it is demonstrated that as the number of waters in simple cyclic clusters increases, the hydrogen bonds strengthen, from the dimer through the cyclic hexamer. However, for the more complex hexamer isomers, the strength of the hydrogen bonds varies, depending on whether the cluster contains double acceptors and/or double donors. The QUAO analysis also reveals the three-center bonding nature of hydrogen bonds in water clusters.

Appraisal of individual hydrogen bond strengths and cooperativity in ammonia clusters via a molecular tailoring approach

In this work, we propose and test a method, based on the molecular tailoring approach (MTA), for the evaluation of individual hydrogen bond (HB) energies in ammonia (NH₃)_n clusters. This methodology was tested, in our earlier work, on water clusters. Liquid ammonia being a universal, non-aqueous ionizing solvent, such information of individual HB strength is indispensable in many studies. The estimated HB energies by an MTA-based method, in (NH₃)_n for n = 3-8, were calculated to be in the range of 0.65 to 5.54 kcal mol^-1 with the cooperativity contribution falling between -0.54 and 1.88 kcal mol^-1 both calculated at the MP2(full)/aug-cc-pVTZ level of theory. It is seen that the strong HBs in (NH₃)_n clusters were additionally strengthened by the large contribution of HB cooperativity. The accuracy of these estimated HB energies was validated by approximately estimating the molecular energy of a given cluster by adding the sum of HB energies to the sum of monomer energies. This approximately estimated molecular energy of a given cluster was found to be in excellent agreement with the actual calculated values. The negligibly small difference (less than 5.6 kcal mol^-1) in these two values suggests that the estimated individual HB energies in ammonia clusters are quite reliable. Furthermore, these estimated HB energies by MTA are in excellent qualitative agreement with the other indirect measures of HB strength, such as HB bond distances and angles, N-H stretching frequency and the electron density values at the (3,-1) bond critical points.

On the Relationship between Hydrogen Bond Strength and the Formation Energy in Resonance-Assisted Hydrogen Bonds

Resonance-assisted hydrogen bonds (RAHB) are intramolecular contacts that are characterised by being particularly energetic. This fact is often attributed to the delocalisation of π electrons in the system. In the present article, we assess this thesis via the examination of the effect of electron-withdrawing and electron-donating groups, namely -F, -Cl, -Br, -CF3, -N(CH3)2, -OCH3, -NHCOCH3 on the strength of the RAHB in malondialdehyde by using the Quantum Theory of Atoms in Molecules (QTAIM) and the Interacting Quantum Atoms (IQA) analyses. We show that the influence of the investigated substituents on the strength of the investigated RAHBs depends largely on its position within the π skeleton. We also examine the relationship between the formation energy of the RAHB and the hydrogen bond interaction energy as defined by the IQA method of wave function analysis. We demonstrate that these substituents can have different effects on the formation and interaction energies, casting doubts regarding the use of different parameters as indicators of the RAHB formation energies. Finally, we also demonstrate how the energy density can offer an estimation of the IQA interaction energy, and therefore of the HB strength, at a reduced computational cost for these important interactions. We expected that the results reported herein will provide a valuable understanding in the assessment of the energetics of RAHB and other intramolecular interactions.

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.