Thermodynamics of forming water clusters at various temperatures and pressures by Gaussian-2, Gaussian-3, complete basis set-QB3, and complete basis set-APNO model chemistries; implications for atmospheric chemistry

The Gaussian-2, Gaussian-3, complete basis set- (CBS-) QB3, and CBS-APNO methods have been used to calculate Delta H degrees and Delta G degrees values for neutral clusters of water, (H(2)O)(n), where n = 2-6. The structures are similar to those determined from experiment and from previous high-level calculations. The thermodynamic calculations by the G2, G3, and CBS-APNO methods compare well against the estimated MP2(CBS) limit. The cyclic pentamer and hexamer structures release the most heat per hydrogen bond formed of any of the clusters. While the cage and prism forms of the hexamer are the lowest energy structures at very low temperatures, as temperature is increased the cyclic structure is favored. The free energies of cluster formation at different temperatures reveal interesting insights, the most striking being that the cyclic trimer, cyclic tetramer, and cyclic pentamer, like the dimer, should be detectable in the lower troposphere. We predict water dimer concentrations of 9 x 10(14) molecules/cm(3), water trimer concentrations of 2.6 x 10(12) molecules/cm(3), tetramer concentrations of approximately 5.8 x 10(11) molecules/cm(3), and pentamer concentrations of approximately 3.5 x 10(10) molecules/cm(3) in saturated air at 298 K. These results have important implications for understanding the gas-phase chemistry of the lower troposphere.     

In search of CS2(H2O)(n=1-4) clusters

Gaussian-3 and MP2/aug-cc-pVnZ methods have been used to calculate geometries and thermochemistry of CS(2)(H2O)n, where n=1-4. An extensive molecular dynamics search followed by optimization using these two methods located two dimers, six trimers, six tetramers, and two pentamers. The MP2/aug-cc-pVDZ structure matched best with the experimental result for the CS(2)(H2O) dimer, showing that diffuse functions are necessary to model the interactions found in this complex. For larger CS(2)(H2O)n clusters, the MP2/aug-cc-pVDZ minima are significantly different from the MP2(full)6-31G* structures, revealing that the G3 model chemistry is not suitable for investigation of sulfur containing van der Waals complexes. Based on the MP2/aug-cc-pVTZ free energies, the concentration of saturated water in the atmosphere and the average amount of CS(2) in the atmosphere, the concentrations of these clusters are predicted to be on the order of 10(5) CS(2)(H2O) clusters.cm(-3) and 10(2) CS(2)(H2O)(2) clusters.cm(-3) at 298.15 K. The MP2/aug-cc-pVDZ scaled harmonic and anharmonic frequencies of the most abundant dimer cluster at 298 K are presented, along with the MP2/aug-cc-pVDZ scaled harmonic frequencies for the CS(2)(H(2)O)(n) structures predicted to be present in a low-temperature molecular beam experiment.     

Electric field effects on water clusters (n = 3-5): systematic ab initio study of structures, energetics, and transition states

The structures, energetics, and transition states of water clusters (trimer to pentamer, n = 3-5) are investigated as a function of electric field by using ab initio calculations. With an increasing strength of the field, the most stable cyclic structures of trimer, tetramer, and pentamer open up to align their dipole moments along the direction of the field. For the lower strength (below 0.3 V/angstroms) of the electric field, the dipole moment of each water monomer is along the same direction with the field, while it retains the cyclic structure. For the higher strength of the field, to have a higher dipole moment for the cluster along the field direction, each cyclic structure opens up to form a linear chain or "water wire." We have investigated the transition state structures between the cyclic and linear forms for the field strengths of 0.3-0.4 V/angstroms where both cyclic and linear forms are energetically comparable.     

Prediction of accurate anharmonic experimental vibrational frequencies for water clusters, (H2O)n, n=2-5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2-5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.     

Hydration of OCS with one to four water molecules in atmospheric and laboratory conditions

Carbonyl sulfide is the most abundant sulfur gas in the atmosphere. We have used MP2 and CCSD(T) theory to study the structures and thermochemistries of carbonyl sulfide interacting with one to four water molecules. We have completed an extensive search for clusters of OCS(H2O)n , where n = 1-4. We located three dimers, two trimers, five tetramers, and four pentamers with the MP2/aug-cc-pVDZ method. In each of the complexes with two or more waters, OCS preferentially interacts with low-energy water clusters. Our results match current theoretical and experimental literature, showing correlation with available geometries and frequencies for the OCS(H2O) species. The CCSD(T)/aug-cc-pVTZ thermochemical values combined with the average amount of OCS and the saturated concentration of H2O in the troposphere, lead to the prediction of 10(6) OCS(H2O) clusters x cm(-3) and 10(2) OCS(H2O)2 clusters x cm(-3) at 298 K. We predict the structures of OCS(H2O)n , n = 1-4 that should predominate in a low-temperature molecular beam and identify specific infrared vibrations that can be used to identify these different clusters.     

Theoretical studies on the bonding and thermodynamic properties of Ge n Si m (m+n=5) clusters: the precursors of germanium/silicon nanomaterials

Theoretical studies on the Ge n Si m clusters have been carried out using advanced ab initio approaches. The lowest energy isomers were determined for the clusters with compositions n+m=2-5. All possible isomers arising due to permutations of Ge and Si atoms were investigated. The L-shaped structure for the trimers, tetragonal with diagonal bond for tetramers, and a trigonal bipyramid for pentamers represent the energy optimized ground state geometries. The bonding analyses revealed that the trimers and tetramers are stabilized through multicenter pi bonding. In pentamers, this stabilizing factor is eliminated due to the further cluster growth. The ionization of clusters does not change their geometrical characteristics. The agreement of the calculated ionization and atomization energies with those obtained from the mass spectrometric studies (through estimated appearance potential) validated the reported structures of the clusters. The bonding properties of these species are discussed using their molecular orbital characteristics and analysis of natural bond orbital population data.     

Density Functional Theory Study of Water Diffusion and Clustering on Pd（111）

1 INTRODUCTION The interaction of water molecules with metal sur- faces plays a vital role in a number of important pro- cesses, such as corrosion, heterogeneous catalysis, electrochemical processes in aqueous solutions, hydrogen production, etc.[1] The structure and pro- perties of water adsorbed on well-defined metal sur- faces have been the subject of numerous experi- mental and theoretical investigations. There have been a number of experimental studies of water on metal surfaces throu…     

An ab initio study of water clusters in gas phase and bulk aqueous media: (H2O)n,n=1–12

Geometries of several clusters of water molecules including single minimum energy structures of n‐mers (n=1–5), several hexamers and two structures of each of heptamer to decamer derived from hexamer cage and hexamer prism were optimized. One structural form of each of 11‐mer and 12‐mer were also studied. The geometry optimization calculations were performed at the RHF/6‐311G* level for all the cases and at the MP2/6‐311++G** level for some selected cases. The optimized cluster geometries were used to calculate total energies of the clusters in gas phase employing the B3LYP density functional method and the 6‐311G* basis set. Frequency analysis was carried out in all the cases to ensure that the optimized geometries corresponded to total energy minima. Zero‐point and thermal free energy corrections were applied for comparison of energies of certain hexamers. The optimized cluster geometries were used to solvate the clusters in bulk water using the polarized continuum model (PCM) of the self‐consistent reaction field (SCRF) theory, the 6‐311G* basis set, and the B3LYP density functional method. For the cases for which MP2/6‐311++G** geometry optimization was performed, solvation calculations in water were also carried out using the B3LYP density functional method, the 6‐311++G** basis set, and the PCM model of SCRF theory, besides the corresponding gas‐phase calculations. It is found that the cage form of water hexamer cluster is most stable in gas phase among the different hexamers, which is in agreement with the earlier theoretical and experimental results. Further, use of a newly defined relative population index (RPI) in terms of successive total energy differences per water molecule for different cluster sizes suggests that stabilities of trimers, hexamers, and nonamers in gas phase and those of hexamers and nonamers in bulk water would be favored while those of pentamer and decamer in both the phases would be relatively disfavored. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 90–104, 2001     

Monte Carlo simulations of equilibrium solubilities and structure of water in n-alkanes and polyethylene

Gibbs ensemble Monte Carlo methods based on a force field that combines the simple point charge [Berendsen et al., in Intermolecular Forces, edited by Pullman (Reidel, Dordrecht, 1981), p. 331] and transferable potentials for phase equilibria [Martin and Siepmann, J. Phys. Chem. B 102, 2659 (1998)] models were used to study the equilibrium properties of binary systems consisting of water and n-alkanes with chain lengths from hexane to hexadecane. In addition, systems where extended linear alkane chains (up to 300 carbon units long) were used to represent amorphous polyethylene were simulated in the presence of water using a connectivity altering osmotic Gibbs ensemble. In these simulations the equilibrium between a liquid water phase and a polymer phase into which water was inserted was studied. The predicted solubilities, which were determined between 350 and 550 K, are in good agreement with experiment, where experimental results are available, and the density of water molecules in the hydrocarbons is approximately 63% as high as in saturated water vapor under the same conditions. At the lower temperatures most of the water exists as monomers; increasing the temperature leads to an increase in the density of water in the alkane phase and hence in the fraction of molecules that participate in clusters. Dimers are the most prevalent clusters in all hydrocarbons and at all temperatures studied, and the fraction of clusters of given size decrease with increasing cluster size. A large fraction of trimers, tetramers, and pentamers, which are the cluster sizes for which topologies have been studied, are cyclic at low temperatures, but at higher temperatures linear structures predominate. The same properties are observed for pure water vapor clusters in equilibrium with the liquid phase, showing that the cluster topologies are not significantly affected by the surrounding hydrocarbon.     

Clusters of hydrated methane sulfonic acid CH3SO3H.(H2O)n (n = 1-5): a theoretical study

Ab initio and density functional methods have been used to examine the structures and energetics of the hydrated clusters of methane sulfonic acid (MSA), CH3SO3H.(H2O)n (n = 1-5). For small clusters with one or two water molecules, the most stable clusters have strong cyclic hydrogen bonds between the proton of OH group in MSA and the water molecules. With three or more water molecules, the proton transfer from MSA to water becomes possible, forming ion-pair structures between CH3SO3- and H3O+ moieties. For MSA.(H2O)3, the energy difference between the most stable ion pair and neutral structures are less than 1 kJ/mol, thus coexistence of neutral and ion-pair isomers are expected. For larger clusters with four and five water molecules, the ion-pair isomers are more stable (>10 kJ/mol) than the neutral ones; thus, proton transfer takes place. The ion-pair clusters can have direct hydrogen bond between CH3SO3- and H3O+ or indirect one through water molecule. For MSA.(H2O)5, the energy difference between ion pairs with direct and indirect hydrogen bonds are less than 1 kJ/mol; namely, the charge separation and acid ionization is energetically possible. The calculated IR spectra of stable isomers of MSA.(H2O)n clusters clearly demonstrate the significant red shift of OH stretching of MSA and hydrogen-bonded OH stretching of water molecules as the size of cluster increases.     

Theoretical studies on photoelectron and IR spectral properties of Br2.-(H2O)n clusters

We report vertical detachment energy (VDE) and IR spectra of Br2.-.(H2O)n clusters (n=1-8) based on first principles electronic structure calculations. Cluster structures and IR spectra are calculated at Becke's half-and-half hybrid exchange-correlation functional (BHHLYP) with a triple split valence basis function, 6-311++G(d,p). VDE for the hydrated clusters is calculated based on second order Moller-Plesset perturbation (MP2) theory with the same set of basis function. On full geometry optimization, it is observed that conformers having interwater hydrogen bonding among solvent water molecules are more stable than the structures having double or single hydrogen bonded structures between the anionic solute, Br2.-, and solvent water molecules. Moreover, a conformer having cyclic interwater hydrogen bonded network is predicted to be more stable for each size hydrated cluster. It is also noticed that up to four solvent H2O units can reside around the solute in a cyclic interwater hydrogen bonded network. The excess electron in these hydrated clusters is localized over the solute atoms. Weighted average VDE is calculated for each size (n) cluster based on statistical population of the conformers at 150 K. A linear relationship is obtained for VDE versus (n+3)(-1/3) and bulk VDE of Br2.- aqueous solution is calculated as 10.01 eV at MP2 level of theory. BHHLYP density functional is seen to make a systematic overestimation in VDE values by approximately 0.5 eV compared to MP2 data in all the hydrated clusters. It is observed that hydration increases VDE of bromine dimer anion system by approximately 6.4 eV. Calculated IR spectra show that the formation of Br2.--water clusters induces large shifts from the normal O-H stretching bands of isolated water keeping bending modes rather insensitive. Hydrated clusters, Br2.-.(H2O)n, show characteristic sharp features of O-H stretching bands of water in the small size clusters.     

Electronic structure and stability of AlnPn (n = 2-4) clusters

The geometry, electronic configurations, harmonic vibrational frequencies, and stability of the structural isomers of aluminum phosphide clusters have been investigated using the density functional theory. For dimers and trimers, the lowest energy structures are cyclic (IIs, IIIs) with D(nh) symmetry. The caged structure with Td symmetry (Xs) lie lowest in energy among the tetramers. The Al--P bond dominates the structures for many isomers so that one preferred dissociation channel is loss of the AlP monomer. The hybridization and chemical bonding in the different structures are also discussed. Comparisons with silicon and boron nitride clusters, the ground state structures of Al(n)P(n) clusters are analogous to those of their corresponding Si(2n) counterparts. This similarity follows the isoelectronic principle.     

Dihydrogen trioxide clusters, (HOOOH)n (n = 2-4), and the hydrogen-bonded complexes of HOOOH with acetone and dimethyl ether: implications for the decomposition of HOOOH

Hydrogen-bonded gas-phase molecular clusters of dihydrogen trioxide (HOOOH) have been investigated using DFT (B3LYP/6-311++G(3df,3pd)) and MP2/6-311++G(3df,3pd) methods. The binding energies, vibrational frequencies, and dipole moments for the various dimer, trimer, and tetramer structures, in which HOOOH acts as a proton donor as well as an acceptor, are reported. The stronger binding interaction in the HOOOH dimer, as compared to that in the analogous cyclic structure of the HOOH dimer, indicates that dihydrogen trioxide is a stronger acid than hydrogen peroxide. A new decomposition pathway for HOOOH was explored. Decomposition occurs via an eight-membered ring transition state for the intermolecular (slightly asynchronous) transfer of two protons between the HOOOH molecules, which form a cyclic dimer, to produce water and singlet oxygen (Delta (1)O 2). This autocatalytic decomposition appears to explain a relatively fast decomposition (Delta H a(298K) = 19.9 kcal/mol, B3LYP/6-311+G(d,p)) of HOOOH in nonpolar (inert) solvents, which might even compete with the water-assisted decomposition of this simplest of polyoxides (Delta H a(298K) = 18.8 kcal/mol for (H 2O) 2-assisted decomposition) in more polar solvents. The formation of relatively strongly hydrogen-bonded complexes between HOOOH and organic oxygen bases, HOOOH-B (B = acetone and dimethyl ether), strongly retards the decomposition in these bases as solvents, most likely by preventing such a proton transfer.     

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.     

氟化氢团簇(HF)_n(n=2～8)拓扑性区别氢键构型的图论列举法和量子化学计算法研究

利用提出的图论程序列举出氟化氢团簇(HF)_n(n=2~8)所有可能存在的拓扑性区别氢键构型,通过精密调查获得有可能存在的拓扑性区别构型,发现了满足HF团簇稳定性的若干条件,在这些条件的基础上编写FORTRAN程序和Python语言执行程序,再用画图软件包Graph Viz2.37自动画出对应的有向图或条件性有向图.以对应的有向图作理论框架,分别利用从头算法Moller-Plesset(MP2)二级微扰方法和密度泛函理论(DFT)方法 B3LYP计算水平的6-31G**(d,p)基组对氟化氢团簇(HF)_n(n=3~7)所有拓扑性区别条件性有向图对应的初始结构进行结构优化并作振动频率分析,获得氟化氢团簇(HF)_n(n=2~7)的最稳定构型,发现了氟化氢团簇的五聚体(HF)_5、六聚体(HF)_6和七聚体(HF)_7等一些新的稳定结构.     

Density functional calculation of the structure and electronic properties of Cu(n)O(n) (n = 1-8) clusters

Ab initio simulations and calculations were used to study the structures and stabilities of copper oxide clusters, Cu(n)O(n) (n = 1-8). The lowest energy structures of neutral and charged copper oxide clusters were determined using primarily the B3LYP/LANL2DZ model chemistry. For n ≥ 4, the clusters are nonplanar. Selected electronic properties including atomization energies, ionization energies, electron affinities, and Bader charges were calculated and examined as a function of n.     

叠氮二氢硼多聚体结构和性质的理论研究(英文)

本文采用DFT-B3LYP方法,以不同基组对叠氮二氢硼多聚体(H2BN3)n (n=1-4)进行计算研究.二聚体(H2BN3)2(C2h对称性)中含B2N2平面四元环结构.船式(Cs对称性)和椅式(C3v对称性)三聚体(H2BN3)3的结合能相近(-122 和 -126 kJ·mol-1),其中均含B3N3六元环结构.拥有B4N4八元环结构的四个四聚体的结合能只有稍微差别.与单体相比,簇合物的结构参数变化较大.由ΔG0T可知,298.2 K下单体形成二聚体在热力学上是不利的,而形成三聚体和四聚体是有利的.     

固态配位聚合物{[Ag2（bpp）2（H2O）2]（chd）.9H2O}n中之新颖“三叉星”笼状（H2O）5水簇的发现

固态配位聚合物{[Ag2（bpp）2（H2O）2]（chd）.9H2O}n（1）（bpp=1,3-二（4-吡啶基）丙烷、H2chd=1,4-环己烷二羧酸）是通过溶液相超声波合成技术获得的,建立在X射线结构分析技术基础上的结构（单斜晶系C2/c空间群,Z=4;晶胞参数a=30.37（2）？,b=9.271（5）？,c=18.89（1）？;β=128.47o;V=4164（4）？3）展现出该化合物捕捉了一条有趣的、由＂三叉星＂笼状五聚体及D2h构型四聚体水簇交替构成的一维水带,该化合物之所以新奇在于其代表了被认为是不太稳定的笼状水五聚体团簇的晶体学范例.第一性原理的密度涵理论（DFT）计算表明,笼状五聚体及D2h构型四聚体水簇的相互连接对于此类欠稳定水簇的稳定性是有益的.     

羰基硼化合物(BCO)n(n==1～12)的理论研究

运用HF/3-21G方法和密度泛函理论(DFT)的B3LYP/6-31G*方法, 对羰基硼化合物(BCO)n(n=1~12)的各种可能结构进行了优化, 对在B3LYP/6-31G*水平上得到的几何构型、电子态、结合能、振动频率、核独立化学位移(NICS)、能量二次差分和热力学性质进行了理论研究, 得到了(BCO)n(n=1~12)结构的稳定性信息. 十二种基态结构都是端配位(μ1-CO), (BCO)n(n=1~3, 5, 6) 的基态是线型或平面结构, (BCO)n(n=4, 7~12)的基态是笼状结构; B—C平均键能呈现奇偶交替现象, 偶数的结构比奇数稳定; 能量二次差分得到同样的结论；羰基的振动频率与实验值非常吻合; 热力学性质的研究对实验具有重要的指导意义.     

Observation of the novel “three-pointed star” cage-like （H2O）5cluster in a polymeric solid {[Ag2（bpp）2（H2O）2]（chd）·9H2O}n

Coordination polymeric solid, {[Ag 2 (bpp) 2 (H 2 O) 2 ](chd)·9H 2 O} n (1) (bpp = 1,3-bis(4-pyridyl) propane, H 2 chd = 1,4-cyclo-hexanedicarboxylic acid), has been obtained by the solution phase ultrasonic synthesis techniques. The structure established through X-ray structural analysis shows that the compound traps an interesting 1D water tape built by the alternating "three-pointed star" cage-like pentameric and D 2h tetrameric clusters (C2/c, Z = 4; a = 30.37(2) , b = 9.271(5) , and c = 18.89(1) ; β = 128.47°; V = 4164(4) 3 ). The novelty of the present complex is the rarely crystallographic example of the cage-shaped water pentamer, which is usually ascribed to a less-stable species. A first-principle density functional theory (DFT) calculation demonstrates that the interconnectivity between cage-like pentamers and D 2h tetramers is beneficial for contribution to the structural stabilization of these less-stable water cluster species.