A Hierarchical Approach to Study the Thermal Behavior of Protonated Water Clusters H+(H2O)n

Structure, stability, infrared spectra, and bonding of OHm(H2O)7 (m = 0, ±1) clusters: ab initio study combining the particle swarm optimization algorithm

The various structural candidates of anionic, neutral, and cationic water clusters OHm(H2O)7 (m = 0, ±1) have been globally predicted by combining the particle swarm optimization method and quantum chemical calculations. Geometry optimization and vibrational analysis for the optimal structures were performed with the MP2/aug-cc-pVDZ method, and the energy profile was further refined at the CCSD(T)/CBS level. Special attention was paid to the relationships between configurations and energies, particularly the first solvation shell coordination number of OH- and OH. For OH-(H2O)7, OH(H2O)7, and OH+(H2O)7 clusters, the most stable species at room temperature are predicted to be the tetra-solvated multi-ring structure A6, the tri-solvated hemibond cage structure N1, and the single five-membered ring structure C2, respectively. The temperature effects on the stability of these three systems were also explored via Gibbs free energies. Furthermore, for the OH-(H2O)7 clusters, the assignments of vibrational transitions in the OH stretching region are in good agreement with the studies of small hydroxide ion-water clusters, and the IR spectra of two isomers (tetra-solvated multi-ring A6 and penta-solvated cage A3) may match future experimental observation well. By topological analysis and reduced density gradient analysis, the structural characteristics and bonding strengths of the studied clusters were investigated. This work indicates the excellent performance of the PSO search algorithm and CALYPSO on water clusters, and may further provide extensive insights into the chemical behavior such as the transport mechanism of OH- ions and OH radicals in the aqueous phase.

A liquid crucible model for aggregation of phenylacetylene in the gas phase

Structural transformation from a π-stacked dimer to an aromatic C–H⋯π trimer and a tetramer.

Revisit the landscape of protonated water clusters H+(H2O)n with n = 10-17: An ab initio global search

Using a genetic algorithm incorporated with density functional theory, we explore the ground state structures of protonated water clusters H⁺(H₂O)_n with n = 10-17. Then we re-optimize the isomers at B97-D/aug-cc-pVDZ level of theory. The extra proton connects with a H₂O molecule to form a H₃O⁺ ion in all H⁺(H₂O)_10-17 clusters. The lowest-energy structures adopt a monocage form at n = 10-16 and core-shell structure at n = 17 based on the MP2/aug-cc-pVTZ//B97-D/aug-cc-pVDZ+ZPE single-point-energy calculation. Using second-order vibrational perturbation theory, we further calculate the infrared spectra with anharmonic correction for the ground state structures of H⁺(H₂O)_10-17 clusters at the PBE0/aug-cc-pVDZ level. The anharmonic correction to the spectra is crucial since it reproduces the experimental results quite well. The extra proton weakens the O-H bond strength in the H₃O⁺ ion since the Wiberg bond order of the O-H bond in the H₃O⁺ ion is smaller than that in H₂O molecules, which causes a red shift of the O-H stretching mode in the H₃O⁺ ion.

Hydrogen bond network structures of protonated short-chain alcohol clusters

Protonated alcohol clusters enable extraction of the physical essence of the nature of hydrogen bond networks.

Solvation energies of the proton in methanol revisited and temperature effects

Various functionals assessing solvation free energies and enthalpies of the proton in methanol.

Exploration of hydrogen bond networks and potential energy surfaces of methanol clusters using a two-stage clustering algorithm

A two-stage algorithm based both on the similarity in shape and hydrogen bond network is developed to explore the potential energy surface of methanol clusters.

Ab Initio Study of Ionized Water Radical Cation (H2O)8+ in Combination with the Particle Swarm Optimization Method

The structures of cationic water clusters (H₂O)₈⁺ have been globally explored by the particle swarm optimization method in combination with quantum chemical calculations. Geometry optimization and vibrational analysis for the 15 most interesting clusters were computed at the MP2/aug-cc-pVDZ level and infrared spectrum calculation at MPW1K/6-311++G** level. Special attention was paid to the relationships between their configurations and energies. Both MP2 and B3LYP-D3 calculations revealed that the cage-like structure is the most stable, which is different from a five-membered ring lowest energy structure but agrees well with a cage-like structure in the literature. Furthermore, our obtained cage-like structure is more stable by 0.87 and 1.23 kcal/mol than the previously reported structures at MP2 and B3LYP-D3 levels, respectively. Interestingly, on the basis of their relative Gibbs free energies and the temperature dependence of populations, the cage-like structure predominates only at very low temperatures, and the most dominating species transforms into a newfound four-membered ring structure from 100 to 400 K, which can contribute greatly to the experimental infrared spectrum. By topological analysis and reduced density gradient analysis, we also investigated the structural characteristics and bonding strengths of these water cluster radical cations.

DFT Study of Hydrogen-Bonding Interaction, Solvation Effect, and Electric-Field Effect on Raman Spectra of Hydrated Proton

Strong hydrogen-bonding interaction and Raman spectra of hydrated proton have been investigated using hybrid density functional theory method B3LYP. The solvation model of density (SMD) approach is employed in the present calculation to simulate hydrated protons in aqueous solution. Focusing on the different hydrogen-bonded Eigen-water and Zundel-water interactions, we present a better assignment of Raman signals of hydrated proton on the basis of vibrational analysis in different environments. Our results showed that B3LYP calculations could give a good prediction for characteristic vibrational frequencies of Eigen and Zundel isomers in liquid phase. The O-H stretching vibrational frequencies from Eigen and Zundel units are very sensitive to hydrogen-bonding interaction with solvent water molecules. Moreover, the solvation effect and the external electric-field effect lead to the proton deviating from the central position of Zundel structure and finally resulting in a transition to Eigen one in aqueous solution. Furthermore, by combining theoretical prediction and Raman scattering theory, we calculate absolute Raman intensities of characteristic signals based on the polarizability tensor derivatives of hydrated proton clusters. This is very helpful to infer the microstructure of hydrated protons in aqueous solution by using Raman measurements.

Ab initio investigation of structure, stability, thermal behavior, bonding, and infrared spectra of ionized water cluster (H2O)6. J Chem Phys 2016;145:154307. [PMID: 27782468 DOI: 10.1063/1.4964860]

The low-lying isomers of cationic water cluster (H₂O)₆⁺ have been globally explored by using particle swarm optimization algorithm in conjunction with quantum chemical calculations. Compared with previous results, our searching method covers a wide range of structural isomers of (H₂O)₆⁺ and therefore turns out to be more effective. With these local minima, geometry optimization and vibrational analysis are performed for the most interesting clusters at second-order Møller-Plesset (MP2)/aug-cc-pVDZ level, and their energies are further refined at MP2/aug-cc-pVTZ and coupled-cluster theory with single, double, and perturbative triple excitations/aug-cc-pVDZ level. The interaction energies using the complete basis set limits at MP2 level are also reported. The relationships between their structure arrangement and their energies are discussed. Based on the results of thermal simulation, structural change from a four-numbered ring to a tree-like structure occurs at T ≈ 45 K, and the relative population of six lowest-free-energy isomers is found to exceed 4% at some point within the studied temperature range. Studies reveal that, among these six isomers, two new-found isomers constitute 10% of isomer population at 180 K, and the experimental spectra can be better explained with inclusions of the two isomers. The molecular orbitals for six representative cationic water clusters are also studied. Through topological and reduced density gradient analysis, we investigated the structural characteristics and the bonding strengths of these water cluster radical cations.

Characterizing the local solvation environment of OH(-) in water clusters with AIMD

In this work, we use ab initio molecular dynamics coupled with metadynamics to explore and characterize the glassy potential energy landscape of the OH(-) in a 20 and 48 water cluster. The structural, energetic, and topological properties of OH(-) are characterized for both clusters and the molecular origins of the IR signatures are examined. We find that in both the small and large clusters, the OH(-) can donate or accept a varying number of hydrogen bonds confirming that the amphiphilic character does not depend on cluster size. However, we highlight some important differences found between the energetic and topological properties of both families of clusters which may have implications on understanding the changes in the solvation structure of OH(-) between bulk and interfacial environments. By studying the IR spectra of smaller subsets of molecules within the 20 water molecule cluster, we find that the IR spectrum of the bare OH(-) as well as the water molecule donating a strong hydrogen bond to it exhibits characteristic absorption along the amphiphilic band between 1500 and 3000 cm(-1) at positions very similar to those found for the entire hydroxide cluster. The results presented here will be useful in the calibration and improvement of both ab initio and semi-empirical methods to model this complex anion.

Stepwise formation of H3O(+)(H2O)n in an ion drift tube: Empirical effective temperature of association/dissociation reaction equilibrium in an electric field

We measured equilibrium constants for H3O(+)(H2O)n-1 + H2O↔H3O(+)(H2O)n (n = 4-9) reactions taking place in an ion drift tube with various applied electric fields at gas temperatures of 238-330 K. The zero-field reaction equilibrium constants were determined by extrapolation of those obtained at non-zero electric fields. From the zero-field reaction equilibrium constants, the standard enthalpy and entropy changes, ΔHn,n-1 (0) and ΔSn,n-1 (0), of stepwise association for n = 4-8 were derived and were in reasonable agreement with those measured in previous studies. We also examined the electric field dependence of the reaction equilibrium constants at non-zero electric fields for n = 4-8. An effective temperature for the reaction equilibrium constants at non-zero electric field was empirically obtained using a parameter describing the electric field dependence of the reaction equilibrium constants. Furthermore, the size dependence of the parameter was thought to reflect the evolution of the hydrogen-bond structure of H3O(+)(H2O)n with the cluster size. The reflection of structural information in the electric field dependence of the reaction equilibria is particularly noteworthy.

Structures and spectroscopy of protonated ammonia clusters at different temperatures

Protonated ammonia clusters are all Eigen structures and the first solvation shell of the related ammonium ion core is saturated by four ammonia molecules.

Disentangling the Contribution of Multiple Isomers to the Infrared Spectrum of the Protonated Water Heptamer

We use infrared/infrared double-resonance population labeling (IR(2)MS(2)) spectroscopy in the spectral region of the free and hydrogen-bonded OH stretching fundamentals (2880-3850 cm(-1)) to identify the number and to isolate the vibrational signatures of individual isomers contributing to the gas-phase IR spectra of the cryogenically cooled protonated water clusters H(+)(H2O)n·H2/D2 with n = 7-10. For n = 7, four isomers are identified and assigned. Surprisingly, the IR(2)MS(2) spectra of the protonated water octa-, nona-, and decamer show no evidence for multiple isomers. The present spectra support the prediction that the quasi-2D to 3D structural transition occurs in between n = 8 and 9 in the cold cluster regime. However, the same models have difficulty explaining the remarkable size dependence of the isomer population reported here.

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.

An infrared spectroscopic and theoretical study on (CH3)3N–H+–(H2O)n, n = 1–22: highly polarized hydrogen bond networks of hydrated clusters

Highly polarized water networks are found in the micro hydaration of protonated trimethylamine.

Temperature dependent structural variations of OH−(H2O)n, n = 4–7: effects on vibrational and photoelectron spectra

In this work, we identified a large number of structurally distinct isomers of midsized deprotonated water clusters using first-principles methods.

Hydrogen-bonded ring closing and opening of protonated methanol clusters H+(CH3OH)n (n = 4–8) with the inert gas tagging

Temperature dependence of hydrogen bond network structures of protonated methanol clusters is explored by IR spectroscopy and DFT simulations.

Structural evolution and solvation of the OH radical in ionized water radical cations (H2O)n+, n = 5–8

Structural evolution of ionized water radical cations (H₂O)_n⁺, n = 5–8, is studied by ab intio methods.

Characterization of a solvent-separated ion-radical pair in cationized water networks: infrared photodissociation and Ar-attachment experiments for water cluster radical cations (H2O)n+(n = 3-8)

We present infrared spectra of nominal water cluster radical cations (H(2)O)(n)(+) (n = 3-8), or to be precise, ion-radical complexes H(+)(H(2)O)(n-1)(OH), with and without an Ar tag. These clusters are closely related to the ionizing radiation-induced processes in water and are a good model to characterize solvation structures of the ion-radical pair. The spectra of Ar-tagged species show narrower bandwidths relative to those of the bare clusters due to the reduced internal energy via an Ar-attachment. The observed spectra are analyzed by comparing with those of the similar system, H(+)(H(2)O)(n), and calculated ones. We find that the observed spectra are attributable to ion-radical-separated motifs in n ≥ 5, as reported in the previous study (Mizuse et al. Chem. Sci.2011, 2, 868-876). Beyond the structural trends found in the previous study, we characterize isomeric structures and determine the number of water molecules between the charged site and the OH radical in each cluster size. In all of the characterized cluster structures (n = 5-8), the most favorable position of OH radical is the next neighbor of the charged site (H(3)O(+) or H(5)O(2)(+)). The positions and cluster structures are governed by the balance among the hydrogen-bonding abilities of the charged site, H(2)O, and OH radical. These findings on the ionized water networks lead to understanding of the detailed processes of ionizing radiation-initiated reactions in liquid water and aqueous solutions.

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.