Ab initio calculations of anharmonic vibrational spectroscopy for hydrogen fluoride (HF)n (n = 3, 4) and mixed hydrogen fluoride/water (HF)n(H2O)n (n = 1, 2, 4) clusters

Anharmonic vibrational frequencies and intensities are computed for hydrogen fluoride clusters (HF)n, with n = 3, 4 and mixed clusters of hydrogen fluoride with water (HF)n(H2O)n where n = 1, 2. For the (HF)4(H2O)4 complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Potential energy surfaces for these systems are obtained at the MP2/TZP level of electronic structure theory. Vibrational states are calculated from the potential surface points using the correlation-corrected vibrational self-consistent field method. The method accounts for the anharmonicities and couplings between all vibrational modes and provides fairly accurate anharmonic vibrational spectra that can be directly compared with experimental results without a need for empirical scaling. For (HF)n, good agreement is found with experimental data. This agreement shows that the M?ller-Plesset (MP2) potential surfaces for these systems are reasonably reliable. The accuracy is best for the stiff intramolecular modes, which indicates the validity of MP2 in describing coupling between intramolecular and intermolecular degrees of freedom. For (HF)n(H2O)n experimental results are unavailable. The computed intramolecular frequencies show a strong dependence on cluster size. Intensity features are predicted for future experiments.     

Study of NH stretching vibrations in small ammonia clusters by infrared spectroscopy in He droplets and ab initio calculations

Infrared spectra of the NH stretching vibrations of (NH3)n clusters (n = 2-4) have been obtained using the helium droplet isolation technique and first principles electronic structure anharmonic calculations. The measured spectra exhibit well-resolved bands, which have been assigned to the nu1, nu3, and 2nu4 modes of the ammonia fragments in the clusters. The formation of a hydrogen bond in ammonia dimers leads to an increase of the infrared intensity by about a factor of 4. In the larger clusters the infrared intensity per hydrogen bond is close to that found in dimers and approaches the value in the NH3 crystal. The intensity of the 2nu4 overtone band in the trimer and tetramer increases by a factor of 10 relative to that in the monomer and dimer, and is comparable to the intensity of the nu1 and nu3 fundamental bands in larger clusters. This indicates the onset of the strong anharmonic coupling of the 2nu4 and nu1 modes in larger clusters. The experimental assignments are compared to the ones obtained from first principles electronic structure anharmonic calculations for the dimer and trimer clusters. The anharmonic calculations were performed at the M?ller-Plesset (MP2) level of electronic structure theory and were based on a second-order perturbative evaluation of rovibrational parameters and their effects on the vibrational spectra and average structures. In general, there is excellent (<20 cm(-1)) agreement between the experimentally measured band origins for the N-H stretching frequencies and the calculated anharmonic vibrational frequencies. However, the calculations were found to overestimate the infrared intensities in clusters by about a factor of 4.     

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.     

Vibrational spectroscopy of (SO4(2-)).(H2O)n clusters, n=1-5: harmonic and anharmonic calculations and experiment

The vibrational spectroscopy of (SO4(2-)).(H2O)n is studied by theoretical calculations for n=1-5, and the results are compared with experiments for n=3-5. The calculations use both ab initio MP2 and DFT/B3LYP potential energy surfaces. Both harmonic and anharmonic calculations are reported, the latter with the CC-VSCF method. The main findings are the following: (1) With one exception (H2O bending mode), the anharmonicity of the observed transitions, all in the experimental window of 540-1850 cm(-1), is negligible. The computed anharmonic coupling suggests that intramolecular vibrational redistribution does not play any role for the observed linewidths. (2) Comparison with experiment at the harmonic level of computed fundamental frequencies indicates that MP2 is significantly more accurate than DFT/B3LYP for these systems. (3) Strong anharmonic effects are, however, calculated for numerous transitions of these systems, which are outside the present observation window. These include fundamentals as well as combination modes. (4) Combination modes for the n=1 and n=2 clusters are computed. Several relatively strong combination transitions are predicted. These show strong anharmonic effects. (5) An interesting effect of the zero point energy (ZPE) on structure is found for (SO4(2-)).(H2O)(5): The global minimum of the potential energy corresponds to a C(s) structure, but with incorporation of ZPE the lowest energy structure is C2v, in accordance with experiment. (6) No stable structures were found for (OH-).(HSO4-).(H2O)n, for n相似文献   

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.     

Ab initio vibrational calculations for H2SO4 and H2SO4 x H2O: spectroscopy and the nature of the anharmonic couplings

Vibrational frequencies for fundamental, overtone, and combination excitations of sulfuric acid (H2SO4) and of sulfuric acid monohydrate cluster (H2SO4 x H2O) are computed directly from ab initio MP2/TZP potential surface points using the correlation-corrected vibrational self-consistent field (CC-VSCF) method, which includes anharmonic effects. The results are compared with experiment. The computed transitions show in nearly all cases good agreement with experimental data and consistent improvement over the harmonic approximation. The CC-VSCF improvements over the harmonic approximation are largest for the overtone and combination excitations and for the OH stretching fundamental. The agreement between the calculations and experiment also supports the validity of the MP2/TZP potential surfaces. Anharmonic coupling between different vibrational modes is found to significantly affect the vibrational frequencies. Analysis of the mean magnitude of the anharmonic coupling interactions between different pairs of normal modes is carried out. The results suggest possible mechanisms for the internal flow of vibrational energy in H2SO4 and H2SO4 x H2O.     

Fourier transform infrared spectroscopy and ab initio theory of acid-hydrogen sulfide clusters: H2S-HCl, D2S-DCl and H2S-(HCl)(2)

The rotationally resolved Fourier transform infrared (FTIR) spectrum of the nu(s) HCl and DCl stretching bands for the hydrogen bonded complex H2S-HCl and its isotopomer D2S-DCl have been observed in a supersonic jet at 0.02 cm(-1) resolution. In the same experimental conditions, two additional bands observed without rotational structure in the HCl range of the dimer have been assigned to the cyclic trimer H2S-(HCl)(2). The multidimensional coupling picture involving the donor stretch mode nu(s) and low frequency intermolecular modes already evidenced in several medium strength hydrogen bonded complexes is beautifully confirmed by the observation of completely separated hot band progressions in the 198 K cell spectrum of both dimers. Based on our anharmonic adiabatic approach for the treatment of the coupled vibrations, absolute vibrational frequencies, diagonal and off-diagonal anharmonicities as well as rovibrational coupling constants obtained from analyses of several 2-D subspaces at MP2 and CCSD(T) level are in excellent agreement with spectroscopic results. In the case of small light complexes, the combination of elevated rotational constants and a negligible contribution of intramolecular vibrational redistribution (IVR) improve the reliability of predissociation lifetime measurements, estimated to 180 ps for H2S-HCl and above 200 ps for D2S-DCl.     

Pseudorotation in the D2O trimer

Far-infrared spectra of Ar/D2O supersonic expansions have been recorded between 39 and 44 cm-1, in which a c-type band centered at 41.1 cm-1 was observed. This new band, attributed to the cyclic water trimer, has the same ground state rotational constants as the (D20)3 98 cm-1 a-type band reported by Liu et al. [J. Am. Chem. Soc. 116 (1994) 3507]. Fits of the observed far-infrared band positions to an effective one-dimensional potential require a substantially larger pseudorotation barrier or effective moment of inertia than predicted theoretically. Estimated frequencies of unobserved transitions in the pseudorotation manifolds of (D2O)3 and (H20)3 are presented as a guide to future experimental searches.     

Vibrational spectroscopy for glycine adsorbed on silicon clusters: Harmonic and anharmonic calculations for models of the Si(1 0 0)-2 × 1 surface

The vibrational spectroscopy of a glycine molecule adsorbed on a silicon surface is studied computationally, using different clusters as models for the surface. Harmonic frequencies are computed using density functional theory (DFT) with the B3LYP functional. Anharmonic frequency calculations are carried out using vibrational self-consistent field (VSCF) algorithms on an improved PM3 potential energy surface. The results are compared with experiments on Glycine@Si(1 0 0)-2 × 1.

The main findings are: (1) Agreement of the computed frequencies with experiment improves with cluster size. (2) The anharmonic calculations are generally in better agreement with experiment than the harmonic ones. The improvements due to anharmonicity are most significant for hydrogenic stretching. (3) An important part of the anharmonic effects is due to anharmonic coupling between different normal modes of the system. (4) The anharmonic coupling between glycine vibrational modes is much larger than the anharmonic coupling between glycine and “phonon” (cluster) modes.

Implications of the results for surface vibrational spectroscopy are discussed.     

Dissociation chemistry of hydrogen halides in water

To understand the mechanism of aqueous acid dissociation, which plays a fundamental role in aqueous chemistry, the ionic dissociation of HX acids (X=F, Cl, Br, and I) in water clusters up to hexamer is examined using density functional theory and M?ller-Plesset second-order perturbation methods (MP2). Further accurate analysis based on the coupled clusters theory with singles and doubles excitations agrees with the MP2 results. The equilibrium structures, binding energies, electronic properties, stretching frequencies, and rotational constants of HX(H(2)O)(n) and X(-)(H(3)O)(+)(H(2)O)(n-1) are calculated. The dissociated structures of HF and HCl can be formed for n>/=4, while those of HBr and HI can be formed for n>/=3. Among these, the dissociated structures of HX (X=Cl, Br, and I) are more stable than the undissociated ones for n>/=4, while such cases for HF would require much more than six water molecules, in agreement with previous reports. The IR spectra of stable clusters including anharmonic frequencies are predicted to facilitate IR experimental studies. Undissociated systems have X-H stretching modes which are highly redshifted by hydration. Dissociated hydrogen halides show three characteristic OH stretching modes of hydronium moiety, which are redshifted from the OH stretching modes of water molecules.     

Ab initio and DFT studies on van der Waals trimers: the OCS.(CO2)2 complexes

Ab initio calculations [MP2, MP4SDTQ, and QCISD(T)] using different basis sets [6-31G(d,p), cc-pVXZ (X = D, T, Q), and aug-cc-pVDZ] and density functional theory [B3LYP/6-31G(d,p)] calculations were carried out to study the OCS.(CO2)2 van der Waals trimer. The DFT has proved inappropriate to the study of this type of systems where the dispersion forces are expected to play a relevant role. Three minima isomers (two noncyclic and one cyclic) were located and characterized. The most stable isomer exhibits a noncyclic barrel-like structure whose bond lengths, angles, rotational constants, and dipole moment agree quite well with the corresponding experimental values of the only structure observed in recent microwave spectroscopic studies. The energetic proximity of the three isomers, with stabilization energies of 1442, 1371, and 1307 cm-1, respectively, at the CBS-MP2/cc-pVXZ (X = D, T, Q) level, strongly suggests that the two unobserved structures should also be detected as in the case of the (CO2)3 trimer where both noncyclic and cyclic isomers have been reported to exist. The many-body symmetry-adapted perturbation theory is employed to analyze the nature of the interactions leading to the formation of the different structures. The three-body contributions are small and stabilizing for the two most stable structures and almost negligible for the cyclic isomer.     

Vibrational overtone spectrum of matrix isolated cis, cis-HOONO

Cis, cis-peroxynitrous acid is known to be an intermediate in atmospheric reactions between OH and NO2 as well as HOO and NO. The infrared absorption spectra of matrix-isolated cc-HOONO and cc-DOONO in argon have been observed in the range of 500-8000 cm-1. Besides the seven fundamental vibrational modes that have been assigned earlier for this molecule [Zhang et al., J. Chem. Phys. 124, 084305 (2006)], more than 50 of the overtone and combination bands have been observed for cc-HOONO and cc-DOONO. Ab initio CCSD(T)/atomic natural orbital anharmonic force field calculations were used to help guide the assignments. Based on this study of the vibrational overtone transitions of cis, cis-HOONO that go as high as 8000 cm-1 and the earlier paper on the vibrational fundamentals, we conclude that the CCSD(T)/ANO anharmonic frequencies seem to correct to +/-35 cm-1. The success of the theoretically predicted anharmonic frequencies {upsilon} in assigning overtone spectra of HOONO up to 8000 cm-1 suggests that the CCSD(T)/ANO method is producing a reliable potential energy surface for this reactive molecule.     

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.     

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.     

Theoretical Study of Small Water Clusters of Sulfur Dioxide

The intermolecular clusters of sulfur dioxide with water, SO2(H2O)n (n = 2～5), are studied by using B3LYP density functional theory and MP2 ab initio methods along with the large basis sets (6-311++G(d,p) and aug-cc-pVDZ). The equilibrium geometries, intermolecular binding energies, and anharmonic frequencies of the clusters are calculated and compared with those of pure water clusters and available experiments. SO2 tends to form cyclic hydrogen-bonded complexes with two or three water molecules. In the larger clusters, however, water molecules begin to retain the structure of pure water clusters and segregate from SO2. Infrared absorption assignments for the small clusters are discussed to resolve a possible incorrect assignment in a recent spectroscopic experiment on the clusters.     

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.     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Dimethylnitroamine Clusters

Ab initio SCF and Mφller-Plesset correlation correction methods in combination with counterpose procedure for BSSE correction have been applied to the theroetical studying of dimethylnitroamine and its dimers and trimers.Three optimized stable dimers and two trimers have been obtained.The corrected binding energies of the most stable dimer and trimer were predicted to be -24.68kJ/mol and -47.27kJ/mol,respectively at the MP2/6-31G^*//HF/6-31G^* level.The proportion of correlated interation energies to their total interaction energies for all clusters was at least 29.3 percent,and the BSSE of ΔE(MP2) was at least 10.0kJ/mol.Dispersion and/or electrostatic force were dominant in all clusters.There exist cooperative effects in both the chain and the cyclic trimers.The vibrational frequencies associated with N-O stretches or wags exhibit slight red shifts,but the modes associated with the motion of hydrogen atoms of the methyl group show somewhat blue shifts with respect to those of monomer.Thermodynamic properties of dimethylnitroamine and its clusters at different temperatures have been calculated on the basis of vibrational analyses.The changes of the Gibbs free energies for the aggregation from monomer to the most stable dimer and trimer were predicted to be 14.37kJ/mol and 30.40kJ/mol,respectively,at 1 atm and 298.15K.     

Anharmonic analysis of the vibrational spectra of some cyanides and related molecules of astrophysical importance

A detailed analysis of the vibrational spectra of carbonyl cyanide, diethynyl ketone and acetyl cyanide has been conducted in harmonic and anharmonic approximations. RHF, MP2 and density functional theory (DFT) methods with 6-311++G(2df,2p) basis sets and B3LYP functionals have been employed. Spectroscopic constants such as anharmonicity constants, rotational and centrifugal distortion constants, rotation-vibration coupling constants and Coriolis coupling coefficients have been calculated for each molecule and compared with the experimental data, where available. A close agreement between the calculated and experimental values of the spectroscopic constants has been obtained. Complete assignments have been provided to the fundamental bands, overtones and combination tones of the molecules. Density functional theory based anharmonic frequencies compare well with the experimental frequencies within +/-18 cm(-1) on an average. RHF and MP2 methods, however, give much higher values for the frequencies that need scaling even in the anharmonic approximation.     

Infrared spectra and ab initio calculations for the F- -(CH4)n (n = 1-8) anion clusters

Infrared spectra of mass-selected F- -(CH4)n (n = 1-8) clusters are recorded in the CH stretching region (2500-3100 cm-1). Spectra for the n = 1-3 clusters are interpreted with the aid of ab initio calculations at the MP2/6-311++G(2df 2p) level, which suggest that the CH4 ligands bind to F- by equivalent, linear hydrogen bonds. Anharmonic frequencies for CH4 and F--CH4 are determined using the vibrational self-consistent field method with second-order perturbation theory correction. The n = 1 complex is predicted to have a C3v structure with a single CH group hydrogen bonded to F-. Its spectrum exhibits a parallel band associated with a stretching vibration of the hydrogen-bonded CH group that is red-shifted by 380 cm-1 from the nu1 band of free CH4 and a perpendicular band associated with the asymmetric stretching motion of the nonbonded CH groups, slightly red-shifted from the nu3 band of free CH4. As n increases, additional vibrational bands appear as a result of Fermi resonances between the hydrogen-bonded CH stretching vibrational mode and the 2nu4 overtone and nu2+nu4 combination levels of the methane solvent molecules. For clusters with n < or = 8, it appears that the CH4 molecules are accommodated in the first solvation shell, each being attached to the F- anion by equivalent hydrogen bonds.     

Ab initio investigation of vibrational spectra of water-(CO2)n complexes (n = 1, 2)

In this paper, we have calculated using the ab initio method the IR vibrational spectra of complexes of CO2 formed with water (sp3 O-donating atom). Binding energies and structures of the CO2-H2O and water-(CO2)2 complexes have been determined at the second-order level of the Moller-Plesset perturbation theory (MP2) using Dunning's basis sets. The results are presented and critically discussed in terms of the nature of the water-CO2 interactions, electron donor acceptor (EDA) and weak O...H-O interactions. For water-(CO2)2 trimer, it is also shown that the contribution to the interaction energy of the irreducible three-bodies remains relatively negligible. We have analyzed the frequency shifts and the IR and Raman intensity variations under the complex formation. We have particularly emphasized the splitting of the 2 bending mode of CO2 and stretching modes of water, which have been revealed as the most pertinent probes to assess the nature of the forces involved in the different complexes. Finally, because water can play the role of Lewis base and acid as well, we found that weak O...H-O interactions can cooperate with EDA interactions in trimer, leading to very specific spectral signatures that are further discussed.