Global Search for Minimum Energy (H2O)n Clusters, n = 3−5

Revisiting the IR Spectra of H2S in an Argon Matrix: Identification of (H2S)n (n ≤ 4) Clusters via the Spectral Assignment of νS-H Transitions

Small-molecular clusters of H2S up to tetramer have been experimentally identified in a cold and solid argon matrix by the spectral assignment of νS-H fundamental transitions for different H2S:argon mixing ratios. Normal-mode frequency calculations at the MP2-CP/aug-cc-pV(Q + d)Z level have been used to support the spectral assignments. In addition, modulations in relative populations of different clusters due to the annealing of the deposited matrix and the preheating of the H2S-argon gas mixture before deposition reinforced the spectral assignments. Variations in mixing ratio, annealing of the matrix, and preheating of the gas mixture have also been used, in a combined manner, to unambiguously identify the ν1 and ν3 bands of the H2S monomer, which has been a matter of dispute for a long period. The two bands have been identified at 2634.4 and 2648.0 cm-1, respectively, while three bands at 2581.5, 2568.4, and 2547.6 cm-1 have been assigned to H-bonded dimers, cyclic trimers, and cyclic tetramers, respectively. Multiple bands within 2550-2580 cm-1 have been assigned to caged tetramers. The cooperative strengthening of S-H···S H bonds in cyclic H2S clusters was evident from the linear increment in νS-H spectral shifts with cluster size.

Water Arrangements upon Interaction with a Rigid Solute: Multiconfigurational Fenchone-(H2O)4-7 Hydrates

Insight into the arrangements of water molecules around solutes is important to understand how solvation proceeds and to build reliable models to describe water-solute interactions. We report the stepwise solvation of fenchone, a biogenic ketone, with 4-7 water molecules. Multiple hydrates were observed using broadband rotational spectroscopy, and the configurations of four fenchone-(H2O)4, three fenchone-(H2O)5, two fenchone-(H2O)6, and one fenchone-(H2O)7 complexes were characterized from the analysis of their rotational spectra in combination with quantum-chemical calculations. Interactions with fenchone deeply perturb water configurations compared with the pure water tetramer and pentamer. In two fenchone-(H2O)4 complexes, the water tetramer adopts completely new arrangements, and in fenchone-(H2O)5, the water pentamer is no longer close to being planar. The water hexamer interacts with fenchone as the least abundant book isomer, while the water heptamer adopts a distorted prism structure, which forms a water cube when including the fenchone oxygen in the hydrogen bonding network. Differences in hydrogen bonding networks compared with those of pure water clusters show the influence of fenchone's topology. Specifically, all observed hydrates except one show two water molecules binding to fenchone through each oxygen lone pair. The observation of several water arrangements for fenchone-(H2O)4-7 complexes highlights water adaptability and provides insight into the solvation process.

Insight into the Binding of Argon to Cyclic Water Clusters from Symmetry-Adapted Perturbation Theory

This work systematically examines the interactions between a single argon atom and the edges and faces of cyclic H2O clusters containing three-five water molecules (Ar(H2O)n=3-5). Full geometry optimizations and subsequent harmonic vibrational frequency computations were performed using MP2 with a triple-ζ correlation consistent basis set augmented with diffuse functions on the heavy atoms (cc-pVTZ for H and aug-cc-pVTZ for O and Ar; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the two-body-many-body (2b:Mb) and three-body-many-body (3b:Mb) techniques; here, high-level CCSD(T) computations capture up through the two-body or three-body contributions from the many-body expansion, respectively, while less demanding MP2 computations recover all higher-order contributions. Five unique stationary points have been identified in which Ar binds to the cyclic water trimer, along with four for (H2O)4 and three for (H2O)5. To the best of our knowledge, eleven of these twelve structures have been characterized here for the first time. Ar consistently binds more strongly to the faces than the edges of the cyclic (H2O)n clusters, by as much as a factor of two. The 3b:Mb electronic energies computed with the haTZ basis set indicate that Ar binds to the faces of the water clusters by at least 3 kJ mol-1 and by nearly 6 kJ mol-1 for one Ar(H2O)5 complex. An analysis of the interaction energies for the different binding motifs based on symmetry-adapted perturbation theory (SAPT) indicates that dispersion interactions are primarily responsible for the observed trends. The binding of a single Ar atom to a face of these cyclic water clusters can induce perturbations to the harmonic vibrational frequencies on the order of 5 cm-1 for some hydrogen-bonded OH stretching frequencies.

Prebiotic dimer and trimer peptide formation in gas-phase atmospheric nanoclusters of water

Insight into the origin of prebiotic molecules is key to our understanding of how living systems evolved into the complex network of biological processes on Earth. By modelling diglycine and triglycine peptide formation in the prebiotic atmosphere, we provide a plausible pathway for peptide growth. By examining different transition states (TSs), we conclude that the formation of diglycine and triglycine in atmospheric nanoclusters of water in the prebiotic atmosphere kinetically favors peptide growth by an N-to-C synthesis of glycines through a trans conformation. Addition of water stabilizes the TS structures and lowers the Gibbs free activation energies. At temperatures that model the prebiotic atmosphere, the free energies of activation with a six water nanocluster as part of the TS are predicted to be 16 kcal mol-1 relative to the prereactive complex. Examination of the trans vs. cis six water transition states reveals that a homodromic water network that maximizes the acceptor/donor nature of the six waters is responsible for enhanced kinetic favorability of the trans N-to-C pathway. Compared to the non-hydrated trans TS, the trans six-water TS accelerates the reaction of diglycine and glycine to form triglycine by 13 orders of magnitude at 217 K. Nature uses the trans N-to-C pathway to synthesize proteins in the ribosome, and we note the similarities in hydrogen bond stabilization between the transition state for peptide synthesis in the ribosome and the transition states formed in nanoclusters of water in the same pathway. These results support the hypothesis that small oligomers formed in the prebiotic atmosphere and rained onto earth's surface.

Amino Acids Compete with Ammonia in Sulfuric Acid-Based Atmospheric Aerosol Prenucleation: The Case of Glycine and Serine

We present a computational investigation of the sulfuric acid, glycine, serine, ammonia, and water system to understand if this system can form prenucleation clusters, which are precursors to larger aerosols in the atmosphere. We have performed a comprehensive configurational search of all possible clusters in this system, starting with the four different monomers and zero to five waters. Accurate Gibbs free energies of formation have been calculated with the DLPNO-CCSD(T)/complete basis set (CBS) method on ωb97xd/6-31++G** geometries. For the dry dimers of sulfuric acid, the weakest base, serine, is found to form the most stable complex, which is a consequence of the strong di-ionic complex formed between the bisulfate ion and the protonated serine cation. For the dry dimers without sulfuric acid, the glycine-serine complex is more stable than the glycine-ammonia or serine-ammonia complexes, stemming from the detailed structure and not related to base strength. For the larger complexes, sulfuric acid deprotonates and the proton is shifted to glycine, serine, or ammonia. The two amino acids and ammonia are almost interchangeable and there is no easy way to predict which molecule will be protonated without the calculated results. Assuming reasonable starting concentrations and a closed system of sulfuric acid, glycine, serine, ammonia, and five waters, we predict the concentrations of all possible complexes at two temperatures spanning the troposphere. The most negative ΔG° values are a function of the detailed molecular interactions of these clusters. These details are more important than the base strength of ammonia, glycine, and serine.

Raman study of water deposited in solid argon matrix

New Raman data are presented concerning H₂O and D₂O water aggregation in argon matrix having the ratio of number of argon atoms to water molecules close to 40:1. Experiments were conducted at temperatures from 8 K to 34 K allowing observation of OH and OD stretching vibrations of water monomers, dimers, trimers and higher multimers, as well as broad bands corresponding to solid amorphous water. Molecular dynamics simulations were performed for thirteen or sometimes fourteen water molecules dispersed among 500 argon atoms. Resulting final configurations included dimers, trimers, tetramers and pentamers, all in open chain configurations which upon optimization resulted in mostly cyclic conformations. Observed OH stretching vibrations were assigned by comparing calculated normal modes in harmonic approximation at the B3LYP/aug-cc-pVDZ and PBEPBE1/aug-cc-pVDZ level of theory with our data and previously observed bands from infrared matrix isolation studies and Raman jet cooled experiments. Raman bands assigned to water multimers in argon matrix are shifted 20 to 25 cm^-1 towards lower wavenumbers with respect to the positions of OH stretching vibrations of almost free water clusters.

A unified cost-effective method for the construction of reliable potential energy surfaces for H2S and H2O clusters

A DFT-based methodology has been used to construct the potential energy surface of H2S clusters up to pentamers. Geometrical parameters and energetics show very good agreement with the existing experimental and high-level theoretical results. Distinct stable conformers of three dimers, six trimers, eleven tetramers and twenty-three pentamers have been identified. Both S-HS H-bond and SS interactions are identified in dimers, trimers and pentamers, while no SS interactions could be found in any of the 11 tetramer conformers. The binding energies of the most stable dimer, trimer, tetramer and pentamer are -1.66, -5.21, -8.57 and -12.54 kcal mol-1, respectively. The PES has been found to be exceedingly flat and the energy gap between the most and the least stable conformers was found to be only 0.09, 2.13, 1.65 and 1.13 kcal mol-1, from the dimer to the pentamer, respectively. The proposed method has also been used for water clusters up to the pentamer. The results obtained were found to agree closely with the existing results. Only one conformer was found for the water dimer, whereas four, five and fifteen conformers were obtained for the trimer, tetramer and pentamer, respectively. Atoms in molecular calculations were found to corroborate with the geometric and energetic results for both clusters.

Molecular Dynamics Simulation of the Excited-State Proton Transfer Mechanism in 3-Hydroxyflavone Using Explicit Hydration Models

3-Hydroxyflavon (3-HF) represents an interesting paradigmatic compound to study excited-state intramolecular proton transfer (ESIPT) and intermolecular (ESInterPT) processes to explain the experimentally observed dual fluorescence in solvents containing protic contamination (water) as opposed to single fluorescence in highly purified nonpolar solvents. In this work, adiabatic on-the-fly molecular dynamics simulations have been performed for isolated 3-HF in an aqueous solution using a polarizable continuum model and including explicit water molecules to represent adequately hydrogen bonding. For the calculation of the excited state, time-dependent density functional theory and the Becke-3-Lee-Yang-Parr (B3LYP) functional have been used. For the isolated 3-HF, ultrafast ESIPT from the enol group to the neighboring keto group has been observed. The calculated PT time of 48 fs agrees well with the experimental value of 39 fs. Addition of one water molecule quenches this ESIPT process but shows an intermolecular concerted or stepwise tautomerization process via the bridging water molecule. Adding a second or more water molecules inhibits this ESInterPT process to a large degree. Most of the trajectories do not show any PT, preserving the initial excited-state enol structure, which is the origin of the violet-blue fluorescence appearing in the solvents contaminated with protic components.

Stabilization of Near Identical Hydrogen Bonded Octameric Water Clusters in Crystal Structures of Three Distinct Non-Charged Polyamide Macrocyclic Host Molecules

In this paper, we present a comparative analysis of the solid state structures of three well-resolved hydrates of macrocyclic host molecules 1a, 1b, and 2 containing an intrannular amide-aryl substituent (lariat arm) connected to a fixed 26-membered ring in a normal (-NHCOAr, hosts 1a and 1b) or reverse manner (-CONHAr, host 2). Despite different chemical structures, these hosts crystallize as isostructural tetrahydrates in the same P-1 space group. Moreover, their crystals exhibit identical hydrogen bond motifs resulting in a stabilization of an almost identical unusual octameric water cluster built from the cyclic tetramer core and four water molecules, attached sequentially in an "up-and-down" manner. Further analysis reveals that, among the series, the structure of host 2 provides the most suitable environment for the accommodation of this type of water cluster.

Temperature Dependence of Intramolecular Vibrational Bands in Small Water Clusters

Cyclic water clusters are pivotal for understanding atmospheric reactions as well as liquid water, yet the temperature (T) dependence of their dynamics and spectroscopy is poorly studied. The development of highly accurate water potentials, such as MB-pol, partly rectifies this. It remains to account for the quantum nuclear effects (NQE), because quantum nuclear dynamics become increasingly inaccurate at low temperatures. From a practical point of view, we find that NQE can be accounted for simply by subtracting a constant from the frequencies obtained from the velocity autocorrelation functions (VACF) of classical nuclear dynamics, resulting in unprecedented agreement with experiment, mostly within 5 cm^-1. We have performed classical simulations of (H₂O)_n clusters (n = 2-5) from 20 K and up to their melting temperature, calculating both all-atom and partial VACF, thus generating the temperature dependence of the vibrational frequencies (IR and Raman bands). Focusing on the hydrogen-bonded (HBed) OH stretch and HOH bend, we find opposing T dependencies. The HBed OH modes blue shift linearly with T, attributed to ring expansion rather than any specific conformational change. The lowest-frequency Raman concerted mode is predicted to show the largest such shift. In contrast, the HOH bend undergoes a red-shift, with the highest frequency concerted band undergoing the largest red-shift. These results can be explained by a coupled-oscillator model for n hydrogen atoms on a ring, constrained to move either tangentially (stretch) or perpendicularly (bend) to the ring. With increasing temperature and weakening of HBs, the intrinsic force constant increases (stretch) or remains constant (bend), while the nearest-neighbor coupling constant decreases, and this results in the interesting behavior revealed herein. T-dependent Raman studies are required for testing some of these predictions.

Ab Initio Molecular Dynamics Investigation of the Electronic and Structural Stability of Anionic O2-(H2O)n, n = 1-16 Clusters

We report an ab initio molecular dynamics investigation of the electronic and structural stability of water molecules binding to a nucleation O₂^- particle, O₂^-(H₂O)_n with n = 1-16, to ascertain the factors that create particularly stable species. Our results compare well with previous experimental and theoretical reports for clusters with less water content, find three new geometries for species with 7, 9, and 10 water molecules, and determine that 8, 11, 13, and 15 water molecules form remarkably stable structures around O₂^-. These special clusters correspond to well-defined compact structures formed by cubes and four-member rings made of water's hydrogen bonds interacting with a negative kernel formed by O₂^- with five water molecules, O₂^-(H₂O)₅, in which the negative charge is localized in the first four water molecules, while the fifth molecule provides geometrical stability. We assess the clusters' energetic stability based on dissociation energies, analyze electron detachment energies to understand its geometrical evolution, and investigate its charge distribution based upon isosurfaces of the highest occupied molecular orbital (HOMO). This research can help provide theoretical insight into the starting steps of nucleation of water clusters around ionic particles.

Hydrogen-bond pattern to characterize water network

Hydrogen-bond (HB) patterns correspond to topologically distinct isomers of water clusters, and can be expressed by digraphs. The HB pattern is used to divide the configuration space of water cluster at a finite temperature. The populations of the HB patterns are transformed into the relative Helmholtz energies. The method is based on the combination of molecular simulation with graph theory. At a finite temperature it can be observed that other isomers than local minimum structures on the potential energy surface are highly populated. The dipole moment of a constituent molecule in a water cluster is enhanced depending on the local HB network around the water molecule. Rooted digraph is used to represent topologically distinct isomers of protonated water (PW) clusters. O–H bonds of PW clusters are classified into 10 topological types based on the combination of the local HB types of the contributing water molecules to the O–H bond. If the topological type is the same, vibrational frequencies of those O–H bonds of PW clusters are similar even in different isomers; i.e. vibrational frequency of O–H bond is transferable, and can be used as a vibrational spectral signature of PW clusters.

Stimulus-responsive reversible thermochromism and exciplex emission of a Zn(ii) complex and selective sensing of NH3 gas

A Zn(ii) complex exhibits reversible thermochromism and emission switch, and the resulting free radical species responds fast and selectively toward NH₃.

Energetics competition in centrally four-coordinated water clusters and Raman spectroscopic signature for hydrogen bonding

Viaseparating the H-bonded neighbour molecules of centrally four-coordinated water molecules from other molecules in outer cages, the calculations discover these two regions interact competitively with the central molecule.

Water inhibits CO oxidation on gold cations in the gas phase. Structures and binding energies of the sequential addition of CO, H2O, O2, and N2 onto Au+

We report a detailed experimental and theoretical study of the gas phase reactivity of Au⁺ with CO, O₂, N₂ and their mixtures in the presence of a trace amount of water impurity.

Distribution of topologically distinct isomers of water clusters and dipole moments of constituent water molecules at finite atmospheric temperatures

Hydrogen-bonding (HB) patterns correspond to topologically distinct isomers of a water cluster and can be expressed by digraphs. We make use of the HB pattern to divide the configuration space of a water cluster ((H2O)n, n = 3-8) at a finite temperature. Each configuration of a water cluster generated in Monte Carlo (MC) simulation is classified into an HB pattern. The number of observed HB patterns increases exponentially with the cluster size, whereas the population of the most abundant HB pattern decreases. The populations of the HB patterns are transformed into the relative Helmholtz energies. At a finite temperature, it can be observed that isomers other than local minimum structures on the potential energy surface are highly populated. The dipole moment of a constituent molecule in a water cluster is enhanced, depending on the molecular circumstance. The change is rationalized by the difference in the local HB type of the water molecule in the HB network.

EVOLUTIONARY DISCOVERY OF TRANSITION STATES IN WATER CLUSTERS

As a basic Aristotle element, water is the most abundant and more importantly crucial substance on earth. Without water, there would not be any form of life as we know. Understanding many phenomena in water such as water evaporation and ice melting and formation requires a deep understanding of hydrogen bond breaking and formation. In particular transition states play a key role in the understanding of such hydrogen bond behavior. Transition states, unlike other metastable states, are energy maxima along the minimum energy path connecting two isomers of molecular clusters. Geometry optimization of transition state structures, however, is a difficult task, and becomes even more arduous, especially when dealing with complex biochemical systems using first-principles calculations. In this paper, a novel molecular memetic algorithm (MOL-MA) composing of specially designed molecular-based water evolutionary operators coupled with a transition-state-local search solver and valley adaptive clearing scheme for the discovery of multiple precise transition states structures is proposed. The transition states of water clusters up to four water molecules uncovered using MOL-MA are reported. MOL-MA is shown not only to reproduce previously found transition states in water clusters, but also established newly discovered transition states for sizes 2–4 water molecules. The search performance of MOL-MA is also shown to outperform its compeers when pitted against those reported in the literature for finding transition states as well as recent advances in niching algorithms in terms of solution precision, computational effort, and number of transition states uncovered.

Anion Cluster: Assembly of Dihydrogen Phosphates for the Formation of a Cyclic Anion Octamer

Structural characterization of a dihydrogen phosphate complex of triprotonated tris[2-(2-thienylmethylamino)ethyl] amine shows that eight dihydrogen phosphate anions are assembled around the host by strong interactions of H-bond donors and acceptors to form a new type of cyclic anion octamer as (H(2)PO(4) (-))(8), an analogy of cyclic water octamer. The presence of an anion cluster has also been identified by electrospray ionization mass spectrometry and (31)P NMR experiments.

Accurate predictions of water cluster formation, (H₂O)(n=2-10)

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2-10 water molecules. This approach reveals new low energy conformers for (H(2)O)(n=7,9,10). The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Reversal of the relative stability of the isomeric radicals HSO and HOS upon hydration and their reactions with ozone

The radical HSO is an oxidation product of pollutants such as H(2)S and CH(3)SH in Earth's atmosphere. For the first time, the interaction of HSO and its tautomer HOS with single water molecules to yield the hydrates HSO.nH(2)O and HOS.nH(2)O was studied for n = 1-3, applying the high-level G3X(MP2) theory. A large number of structures corresponding to local minima on the potential energy surfaces has been identified. While gaseous HSO is more stable than HOS, the enthalpy diffference between HSO.nH(2)O and HOS.nH(2)O decreases with increasing degree of hydration and becomes practically zero for n = 3. Thus, in aqueous solution as well as in fog and rain droplets, HOS is expected to compete with HSO. The barrier for the tautomerization of HSO to HOS is dramatically lowered by the presence of water molecules since a cyclic transition state allows a concerted proton shift within the system of neighboring hydrogen bonds. The corresponding activation enthalpy of only 73.5 kJ mol(-1) predicted for the transformation of HSO.2H(2)O into HOS.2H(2)O may be compared to the 202 kJ mol(-1) reported for the tautomerization of the unhydrated gaseous HSO/HOS molecules. The impact of water of hydration on the fundamental vibrational modes of HSO and HOS has also been studied. Furthermore, HOS is predicted to dimerize at low temperatures to give two van der Waals molecules with singlet (symmetry C(2)) or triplet configuration (symmetry C(2h)), the latter being more stable than the singlet isomer. The disproportionation of 2HSO to H(2)S and SO(2) is predicted to be exothermic by -263.5 kJ mol(-1). The reaction of HSO with ozone to HSO(2) and O(2) is also strongly exothermic by -274.0 kJ mol(-1) and seems to proceed without any barrier. HOS forms a 1:1 van der Waals complex with O(3); the redox reaction of its two components is calculated as exothermic by -410.9 kJ mol(-1) and results in a rather stable adduct between HOSO and O(2) with the structure of a peroxo isomer of HOSO(3). This unprecedented hydrogen peroxosulfite radical might open a novel route to atmospheric sulfate without the intermediate formation of SO(2) and SO(3).