Energetic and fragmentation stability of water clusters (H2O)n, n=2–30

Ambiguous Faces of Water-Based Inclusion Compounds: L4(4)8(8) Intercalato-Clathrate Hydrate of Pt(II) Complex

Clathrate hydrates are among the most intensively studied H-bond inclusion compounds. Despite the broad definition for this class of compounds, their meaning commonly refers to closed polyhedral nanocages that encapsulate small guest molecules. On the other hand, larger solutes enforce another type of encapsulation because of the solute size effect. Herein, we report a series of structures containing various molecules encapsulated by intercalated water layers constructed of polycyclic moieties of L4(4)8(8) topology. We parametrized the corrugation of individual layers and characterized interactions governing their formation. We suggested that these could be categorized as two-dimensional clathrates based on the character of intra-layer interactions and the effects observed between entrapped molecules and water-based intercalators.

A computational study of H-bonded networks in cyclic water clusters, (H2O)n (n = 3-12)

CONTEXT

We have performed a detailed MM and DFT investigation of neutral water clusters (H2O)n (n = 3-12). Our results show the trend of interaction energies in these clusters as a function of the size of the cluster. They show that the H-bond strength increases with cluster size and that the model of water is better described if two different partial charges are used on the hydrogen, depending on whether hydrogen is H-bonded or not. The average binding enthalpy change due to the formation of H-bonds between water molecules is found to be - 25.9 kJ mol-1 at B3LYP/aug-cc-pVDZ level of theory. We observe the formation of cyclic H-bonded networks through the analysis of frontier orbitals and IR vibrational frequencies spectra. For the water cluster with n = 11, we observe an unusual reduction of the bandgap indicative of a cyclic H-bonded network.

METHODS

Calculations were performed with the MMFF94 force field and the B3LYP method using various large basis sets. Molecular orbital diagrams and population analysis were done using standard tools in Gaussian.

Magic Numbers and Stabilities of Photoionized Water Clusters: Computational and Experimental Characterization of the Nanosolvated Hydronium Ion

The stability and distributions of small water clusters generated in a supersonic beam expansion are interrogated by tunable vacuum ultraviolet (VUV) radiation generated at a synchrotron. Time-of-flight mass spectrometry reveals enhanced population of various protonated water clusters (H⁺(H₂O)_n) based upon ionization energy and photoionization distance from source, suggesting there are "magic" numbers below the traditional n = 21 that predominates in the literature. These intensity distributions suggest that VUV threshold photoionization (11.0-11.5 eV) of neutral water clusters close to the nozzle exit leads to a different nonequilibrium state compared to a skimmed molecular beam. This results in the appearance of a new magic number at 14. Metadynamics conformer searches coupled with modern density functional calculations are used to identify the global minimum energy structures of protonated water clusters between n = 2 and 21, as well as the manifold of low-lying metastable minima. New lowest energy structures are reported for the cases of n = 5, 6, 11, 12, 16, and 18, and special stability is identified by several measures. These theoretical results are in agreement with the experiments performed in this work in that n = 14 is shown to exhibit additional stability, based on the computed second-order stabilization energy relative to most cluster sizes, though not to the extent of the well-known n = 21 cluster. Other cluster sizes that show some additional energetic stability are n = 7, 9, 12, 17, and 19. To gain insight into the balance between ion-water and water-water interactions as a function of the cluster size, an analysis of the effective two-body interactions (which sum exactly to the total interaction energy) was performed. This analysis reveals a crossover as a function of cluster size between a water-hydronium-dominated regime for small clusters and a water-water-dominated regime for larger clusters around n = 17.

Adsorption free energy of phenol onto coronene: Solvent and temperature effects

Molecular modeling can considerably speed up the discovery of materials with high adsorption capacity for wastewater treatment. Despite considerable efforts in computational studies, the molecular modeling of adsorption processes has several limitations in reproducing experimental conditions. Handling the environmental effects (solvent effects) and the temperature effects are part of the important limitations in the literature. In this work, we address these two limitations using the adsorption of phenol onto coronene as case study. In the proposed model, for the solvent effects, we used a hybrid solvation model, with n explicit water molecules and implicit solvation. We increasingly used n=1 to n=12 explicit water molecules. To account for the temperature effects, we evaluated the adsorption efficiency using the adsorption free energy for temperatures varying from 200 to 400K. We generated initial configurations using classical molecular dynamics, before further optimisation at the ωB97XD/aug-cc-pVDZ level of theory. Polarisable continuum solvation model (PCM) is used for the implicit solvation. The adsorption free energy is evaluated to be -1.3kcal/mol at room temperature. It has been found that the adsorption free energy is more negative at low temperatures. Above 360K, the adsorption free energy is found to be positive.

Understanding the Underlying Field Evaporation Mechanism of Pure Water Tips in High Electrical Fields

We investigated the field evaporation process of frozen water in atom probe tomography (APT) by density functional simulations. In previous experiments, a strong tailing effect was observed for peaks caused by the molecular structure (H₂O)_nH⁺, in contrast to other peaks. In purely field-induced and thermally assisted evaporation simulations, we found that chains of protonated water molecules were pulled out of the dielectric surface by up to 6 Å, which are stable over a wide range of field strengths. Therefore, the resulting water clusters experience only part of the acceleration after evaporation compared to molecules evaporating directly from the surface and, thus, exhibit an energy deficit, which explains the tailing effect. Our simulations provide new insight into the complex evaporation behavior of water in high electrical fields and reveal possibilities for adapting the existing reconstruction algorithms.

Hierarchical clustering analysis of hydrogen bond networks in aqueous solutions

To understand the relation between the macroscopic properties and microscopic structure of hydrogen bond networks in solutions, we introduced a hierarchical clustering method to analyze the typical configurations of water clusters in this type of network. Regarding hydrogen bonds as frames, the rings, fragments and clusters are defined and analyzed to provide a comprehensive perspective for the distributional and dynamic characteristics of the hydrogen-bonding network in NaCl solution at different concentrations. The properties of the radial distribution function and hydrogen bonds are first analyzed. Destruction and shorter lifetimes of hydrogen bonds are observed in solutions. In further analysis of the two-dimensional configuration, i.e., ring, and three-dimensional configuration, i.e., fragment, the average number, size and lifetime of these structures consistently decrease as the concentration increases. Ionic effects on disrupting rings and fragments are significant in the first hydration shell, especially with sodium cations, and these effects weaken beyond the first hydration shell. Regarding the clusters obtained using the Louvain algorithm, our results indicate that clusters break and become smaller as the NaCl concentration increases. The presence of ions also leads to the isolation of clusters and therefore the inhibition of changes. The lifetime of clusters increases with NaCl concentration, indicating the slowed breakage and reformation of clusters in NaCl solutions. This method can be further applied to quantitatively characterize hydrogen bond networks to elucidate more properties of aqueous solutions.

Comparative Study of Proton Exchange in Tri- and Hexatitanates: Correlations between Stability and Electronic Properties

A hydrothermal method is considered to be convenient and is extensively used in preparing titanate architectures, but the intermediate and final products are complicated and variable. To date, it is accepted that intermediates are tri- and hexatitanates. Here, atomic structures, energetics, and correlations between stability and electronic properties of proton exchange in tri- and hexatitanates, i.e., Na_2-xH_xTi₃O₇ and Na_2-xH_xTi₆O₁₃, are investigated by first-principles calculations. We found that the bond length of Na-O bonds plays a significant role in determining the activity of tunnel oxygen atoms, while the proton substitution sites are closely related to the activity of tunnel O atom in titanates. As H⁺ concentration increases, the formation energy of Na_2-xH_xTi₃O₇ and Na_2-xH_xTi₆O₁₃ decreases first and then increases, suggesting that completely protonated titanates, i.e., H₂Ti₃O₇ and H₂Ti₆O₁₃, are unstable. However, we found that H⁺ substitution would take place even in an alkaline solution both for Na₂Ti₃O₇ and Na₂Ti₆O₁₃. With a decrease in the pH, the process of H⁺ exchange becomes more energetically favorable. Compared to Na₂Ti₃O₇, Na⁺ ions are more easily exchanged by H⁺ ions in Na₂Ti₆O₁₃ at the same pH value. We found that there is a strong correlation between stability and electronic properties during the Na⁺-H⁺ exchange process. Finally, hydrogen bonds are observed in H₂Ti₃O₇ and Na_2-xH_xTi₆O₁₃ complexes, which make them more stable than Na_2-xH_xTi₃O₇ complexes without H-bonds. All of these findings provide insight into understanding the geometry of possible intermediates in the preparation of titanates and suitable conditions for the synthesis of titanates.

Effect of the co-adsorption of small molecules from air on the properties of penta-graphene and their proton transfer calculation

Penta-graphene has attracted considerable attention due to its unique structure and novel properties. Herein, we studied the effect of the co-adsorption of small molecules from air on the properties of penta-graphene using first-principles calculations. Our results show that oxygen molecules can be self-decomposed on the surface of penta-graphene and the process of O₂ decomposition is an exothermic reaction. On the contrary, the adsorption of H₂O or N₂ molecule on penta-graphene exhibits weak interaction characteristic. For co-adsorption systems, the adsorption of N₂ molecule has no effect on the electronic properties of penta-graphene because the N₂ molecule is more inert than other molecules. Hydrogen bonds (H-bonds) have been observed in the co-adsorption of H₂O and O₂ on penta-graphene. We find that shorter H-bonds lead to higher stability of the systems. We also explore the proton transfer process between H₂O and oxidized penta-graphene. Our results show that the proton transfer process is relatively difficult due to the high energy barrier. However, double-proton transfer is an exothermic process since the energy of the final state is 0.11 eV lower than that of the initial state. These results indicate that the configuration of oxidized penta-graphene is complicated. Our research provides a theoretical basis and important guidance for the experimental synthesis and functionalization of penta-graphene.

Increase of network hydrophilicity from sql to lvt supramolecular isomers of Cu-MOFs with the bifunctional 4-(3,5-dimethyl-1H-pyrazol-4-yl)benzoate linker

A slight difference in the H-bonding of the linker pyrazole-NH group changes the framework hydrophilicity drastically.

Ring-Stacking Water Clusters: Morphology and Stabilities

The structures and interaction energies of water clusters with ring stacking motifs are studied by using ab initio calculations. The structures of the water clusters are constructed by stacking either single rings or multi-rings of tetramer, pentamer, and hexamer. We found that, in the single-ring-stacking motif, the most stable isomers exhibit an alternative clockwise-anticlockwise stacking pattern. We also show that four-layer single-ring-stacking isomers are not energetically favorable in comparison with those of two-layer multi-ring-stacking isomers. The relative stability of the isomers is also analyzed in terms of H-bond strength and elastic distortions of the water molecules.

Structures, relative stability and binding energies of neutral water clusters, (H2O)2–30

We have revised the structures of neutral water clusters, (H₂O)_n=2–30, with the affordable M06-2X functional, presenting up to 25 isomers for each cluster size.

Exploring water radiolysis in proton cancer therapy: Time-dependent, non-adiabatic simulations of H+ + (H2O)1-6

To elucidate microscopic details of proton cancer therapy (PCT), we apply the simplest-level electron nuclear dynamics (SLEND) method to H⁺ + (H₂O)_1-6 at E_Lab = 100 keV. These systems are computationally tractable prototypes to simulate water radiolysis reactions—i.e. the PCT processes that generate the DNA-damaging species against cancerous cells. To capture incipient bulk-water effects, ten (H₂O)_1-6 isomers are considered, ranging from quasi-planar/multiplanar (H₂O)_1-6 to “smallest-drop” prism and cage (H₂O)₆ structures. SLEND is a time-dependent, variational, non-adiabatic and direct method that adopts a nuclear classical-mechanics description and an electronic single-determinantal wavefunction in the Thouless representation. Short-time SLEND/6-31G* (n = 1–6) and /6-31G** (n = 1–5) simulations render cluster-to-projectile 1-electron-transfer (1-ET) total integral cross sections (ICSs) and 1-ET probabilities. In absolute quantitative terms, SLEND/6-31G* 1-ET ICS compares satisfactorily with alternative experimental and theoretical results only available for n = 1 and exhibits almost the same accuracy of the best alternative theoretical result. SLEND/6-31G** overestimates 1-ET ICS for n = 1, but a comparable overestimation is also observed with another theoretical method. An investigation on H⁺ + H indicates that electron direct ionization (DI) becomes significant with the large virtual-space quasi-continuum in large basis sets; thus, SLEND/6-31G** 1-ET ICS is overestimated by DI contributions. The solution to this problem is discussed. In relative quantitative terms, both SLEND/6-31* and /6-31G** 1-ET ICSs precisely fit into physically justified scaling formulae as a function of the cluster size; this indicates SLEND’s suitability for predicting properties of water clusters with varying size. Long-time SLEND/6-31G* (n = 1–4) simulations predict the formation of the DNA-damaging radicals H, OH, O and H₃O. While “smallest-drop” isomers are included, no early manifestations of bulk water PCT properties are observed and simulations with larger water clusters will be needed to capture those effects. This study is the largest SLEND investigation on water radiolysis to date.

Understanding the structure and hydrogen bonding network of (H2O)32 and (H2O)33: an improved Monte Carlo temperature basin paving (MCTBP) method and quantum theory of atoms in molecules (QTAIM) analysis

Large water clusters are of particular interest because of their connection to liquid water and the intricate hydrogen bonding networks they possess.

Internal electric fields in small water clusters [(H2O)n; n = 2-6]

The electric field experienced by a water molecule within a water cluster depends on its position relative to the rest of the water molecules. The stabilization energies and the red-shifts in the donor O-H stretching vibrations in the water clusters increase with the cluster size concomitant with the increase in the electric field experienced by the donor O-H of a particular water molecule due to the hydrogen bonding network. The red-shifts in O-H stretching frequencies show a spread of about ±100 cm(-1) against the corresponding electric fields. Deviations from linearity were marked in the region of 100-160 MV cm(-1), which can be attributed to the strain in the hydrogen bonding network, especially for structures with DDAA and DDA motifs. The linear Stark effect holds up to 200 MV cm(-1) of internal electric field for the average red-shifts in the O-H stretching frequencies, with a Stark tuning rate of 2.4 cm(-1) (MV cm(-1))(-1), suggesting the validity of the classical model in small water clusters.

Delayed luminescence induced by complex domains in water and in TEOS aqueous solutions

Many recent studies on water have conjectured a complex structure composed of hydrogen bonded low- and high-density domains. In this work the structure of pure water and aqueous solutions of silica gel (TEOS) has been investigated by using delayed luminescence, which has previously shown a significant increase in aqueous salt solutions where low-density domain formation is expected. Photon emission shows an Arrhenius trend with an activation energy in water-TEOS solutions larger than in pure water and salt-water solutions. Moreover, delayed photon emission decay shows an intrinsic lifetime of about 5 μs both in solutions and in pure water that, along with secondary lifetimes induced by the presence of TEOS, could be related to the formation of different domains.

Density functional studies of fused dodecahedral and irregular-dodecahedral water cages

The fused cages of dodecahedral and irregular-dodecahedral water cages with the maximum number of t1d hydrogen bonds were studied using the dispersion corrected density functional method.

Structures and relative stabilities of ammonia clusters at different temperatures: DFT vs. ab initio

The global minimum energy structures of (NH₃)_n=2–10are pointed out for the first time at a given temperature.

Facilitating Minima Search for Large Water Clusters at the MP2 Level via Molecular Tailoring

Water clusters (H2O)20 and (H2O)25 are explored at the Møller-Plesset second-order perturbation (MP2) level of theory. Geometry optimization is carried out on favorable structures, initially generated by the temperature basin paving (TBP) method, utilizing the fragment-based molecular tailoring approach (MTA). MTA-based stabilization energies at the complete basis set limit are accurately estimated by grafting the energy correction using a smaller basis set. For prototypical cases, the minima are established via MTA-based vibrational frequency calculations at the MP2/aug-cc-pVDZ level. The potential of MTA in tackling large clusters is further demonstrated by performing geometry optimization at MP2/aug-cc-pVDZ starting with the global minimum of (H2O)30 reported by Monte Carlo (MC) and molecular dynamics (MD) investigations. The present study brings out the efficacy of MTA in performing computationally expensive ab initio calculations with minimal off-the-shelf hardware without significant loss of accuracy.