Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H(5)O(2) (+) using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh et al. [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.     

Structures of HCN-Mgn (n=2-6) complexes from rotationally resolved vibrational spectroscopy and ab initio theory

High-resolution infrared laser spectroscopy has been used to determine the structures of HCN-Mgn complexes formed in helium nanodroplets. The magnesium atoms are first added to the droplets to ensure that the magnesium complexes are preformed before the HCN molecule is added. The vibrational frequencies, structures, and dipole moments of these complexes are found to vary dramatically with cluster size, illustrating the nonadditive nature of the HCN-magnesium interactions. All of the complexes discussed here have the nitrogen end of the HCN pointing towards the magnesium clusters. For Mg3, the HCN binds to the "threefold" site, yielding a symmetric top spectrum. Although the HCN-Mg4 complex also has C3v symmetry, the HCN sits "on-top" of a single magnesium atom. These structures are confirmed by both ab initio calculations and measurements of the dipole moments. Significant charge transfer is observed in the case of HCN-Mg4, indicative of charge donation from the lone pair on the nitrogen of HCN into the lowest unoccupied molecular orbital of the Mg4.     

Overtone vibrational spectroscopy in H(2)-H(2)O complexes: A combined high level theoretical ab initio, dynamical and experimental study

First results are reported on overtone (v(OH) = 2 ← 0) spectroscopy of weakly bound H(2)-H(2)O complexes in a slit supersonic jet, based on a novel combination of (i) vibrationally mediated predissociation of H(2)-H(2)O, followed by (ii) UV photodissociation of the resulting H(2)O, and (iii) UV laser induced fluorescence on the nascent OH radical. In addition, intermolecular dynamical calculations are performed in full 5D on the recent ab initio intermolecular potential of Valiron et al. [J. Chem. Phys. 129, 134306 (2008)] in order to further elucidate the identity of the infrared transitions detected. Excellent agreement is achieved between experimental and theoretical spectral predictions for the most strongly bound van der Waals complex consisting of ortho (I = 1) H(2) and ortho (I = 1) H(2)O (oH(2)-oH(2)O). Specifically, two distinct bands are seen in the oH(2)-oH(2)O spectrum, corresponding to internal rotor states in the upper vibrational manifold of Σ and Π rotational character. However, none of the three other possible nuclear spin modifications (pH(2)-oH(2)O, pH(2)-pH(2)O, or oH(2)-pH(2)O) are observed above current signal to noise level, which for the pH(2) complexes is argued to arise from displacement by oH(2) in the expansion mixture to preferentially form the more strongly bound species. Direct measurement of oH(2)-oH(2)O vibrational predissociation in the time domain reveals lifetimes of 15(2) ns and <5(2) ns for the Σ and Π states, respectively. Theoretical calculations permit the results to be interpreted in terms of near resonant energy levels and intermolecular alignment of the H(2) and H(2)O wavefunctions, providing insight into predissociation dynamical pathways from these metastable levels.     

Flexible H2O2 in water: electronic structure from photoelectron spectroscopy and ab initio calculations

The effect of hydration on the electronic structure of H(2)O(2) is investigated by liquid-jet photoelectron spectroscopy measurements and ab initio calculations. Experimental valence electron binding energies of the H(2)O(2) orbitals in water are, on average, 1.9 eV red-shifted with respect to the gas-phase molecule. A smaller width of the first peak was observed in the photoelectron spectrum from the solution. Our experiment is complemented by simulated photoelectron spectra, calculated at the ab initio level of theory (with EOM-IP-CCSD and DFT methods), and using path-integral sampling of the ground-state density. The observed shift in ionization energy upon solvation is attributed to a combination of nonspecific electrostatic effects (long-range polarization) and of the specific interactions between H(2)O(2) and H(2)O molecules in the first solvation shell. Changes in peak widths are found to result from merging of the two lowest ionized states of H(2)O(2) in water due to conformational changes upon solvation. Hydration effects on H(2)O(2) are stronger than on the H(2)O molecule. In addition to valence spectra, we report oxygen 1s core-level photoelectron spectra from H(2)O(2)(aq), and observed energies and spectral intensities are discussed qualitatively.     

Quantum calculations of vibrational energies of H3O2- on an ab initio potential

We report a full-dimensional potential energy surface for H3O2-, based on fitting 66,965 ab initio electronic energies. A major feature of this potential is a barrier of roughly 200 cm-1 to internal rotation of the two hydroxyl groups about a line connecting the two oxygen atoms and the bridging hydrogen atom. The potential is used in calculations of vibrational energies, performed with the "Reaction Path" version of the code "MULTIMODE". The results are compared to recent infrared messenger experiments and are used to propose interpretations of the experimental results.     

HF(H2O)n clusters with an excess electron: ab initio study

The structures of electron-bound and neutral clusters of HF(H2O)n (n=1-3) were optimized at the level of second-order Moller-Plesset perturbation theory (MP2). Then, the energies were studied using the coupled cluster singles, doubles, and perturbative triples correction [CCSD(T)] method. The vertical detachment energies of the electron-bound clusters for n=1-3 are 60, 180, and approximately 300 meV, respectively. In the case of the n=3, two structures are competing energetically. The electron-bound clusters for n=1 and 2 are 1.5 and 1.8 kcal/mol more stable than the neutral, while that for n=3 is 0.6-0.9 kcal/mol less stable. The excess electron is stabilized in the surface-bound state of the dipole oriented structures of the hydrated acid clusters. Vibrational spectra of the electron-bound clusters are discussed.     

High-level ab initio versus DFT calculations on (H2O2)2 and H2O2–H2O complexes as prototypes of multiple hydrogen bond systems

The performance of B-LYP, B-P86, B3-LYP, B3-P86, and B3-PW91 density functionals to describe multiple hydrogen bond systems was studied. For this purpose we have chosen the dimers of hydrogen peroxide and the hydrogen peroxide–water complexes. The geometries and vibrational frequencies obtained with a 6-311+G(d,p) basis set were compared with those obtained at the MP2 level using the same basis set expansion. The corresponding dimerization energies were obtained using a 6-311+G(3df,2p) basis set and compared with those obtained using the G2(MP2) theory. Red shiftings of the OH donor stretching frequencies were predicted by all approaches investigated; however, in all cases, the DFT values were sizably larger than the MP2 ones. Similarly, the blue shifting of the torsion of the hydrogen peroxide subunit was larger when evaluated at the DFT level. All functionals reproduced the G2(MP2) relative stabilities of the different local minima quite well. With the exception of the B-LYP and B3-PW91 approaches, all functionals yielded binding energies which deviated from the G2(MP2) values by less than 0.5 kcal/mol, provided that G2-type basis sets were used and that the corresponding BSSE corrections were included. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1124–1135     

Pulsed-field ionization electron spectroscopy and ab initio calculations of copper-diazine complexes

Copper complexes of pyrazine (1,4-C4H4N2), pyrimidine (1,3-C4H4N2), and pyridazine (1,2-C4H4N2) are produced in laser-vaporization supersonic molecular beams and studied by pulsed-field ionization zero electron kinetic energy (ZEKE) spectroscopy and second-order Moller-Plesset perturbation theory. Both sigma and pi complexes are considered by these ab initio calculations; only sigma structures are identified in these experiments. Adiabatic ionization energies and metal-ligand vibrational frequencies of the sigma complexes are measured from the ZEKE spectra. Metal-ligand bond dissociation energies of these complexes are obtained from a thermochemical cycle. The ionization energies follow the trend of Cu pyridazine (43,054 cm(-1)) < Cu pyrimidine (45,332 cm(-1)) < Cu pyrazine (46,038 cm(-1)); the bond energies are in the order of Cu pyridazine (56.2 kJ mol(-1)) > Cu pyrazine (48.5 kJ mol(-1)) approximately Cu pyrimidine (46.4 kJ mol(-1)). The stronger binding of pyridazine is due to its larger electric dipole moment and possibly bidentate binding.     

Photoionization and ab initio study of Ba(H2O)n (n = 1-4) clusters

An experimental and theoretical study of the photoionization energies (IE's) of Ba(H(2)O)(n) clusters containing up to n = 4 water molecules has been performed. The clusters were generated by a pick-up source combining laser vaporization with pulsed supersonic expansion, and then photoionized by radiation of 272.5-340 nm. The experimentally determined IE(e)'s for n = 1 to 4 are 4.56 ± 0.05, 4.26 ± 0.05, 3.90 ± 0.05 and 3.71 ± 0.05 eV. This cluster size dependence of IE is reproduced within ±0.06 eV employing the mPW1PW91 density-functional and CCSD(T, Full) quantum-chemical methods combined with the 6-311++G(d,p) basis set for the H and O atoms and three different relativistic effective core potentials for Ba atoms. The calculations indicate that the lowest energy hydration structures represent the most relevant contributions to both the vertical and adiabatic ionization energies. Experimental and theoretical evidence correlates with the progressive surface-delocalization of the electron from the hydration cavity around the Ba atom and suggests that the intra-cluster electron transfer is possible even for small aggregates.     

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.     

Tautomers and conformers of malonamide, NH2-C(O)-CH2-C(O)-NH2: vibrational analysis, NMR spectra and ab initio calculations

The conformational and tautomeric compositions of malonamide, NH2-C(O)-CH2-C(O)-NH2 were determined by vibrational spectroscopy and theoretical calculations (HF/6-31G*, B3PW91/6-31G*). Solid state Fourier transform infrared and Raman spectra were analysed. They reveal the existence of a diketo tautomer. Theoretical calculations predict a diketo structure belonging to the C1 symmetry group. No enol form is present in the molecule in the solid. 13C-NMR studies show only signals of a diketo tautomer.     

Vibrational spectra and structure of CH3Cl:(H2O)2 and CH3Cl:(D2O)2 complexes. IR matrix isolation and ab initio calculations

The infrared spectra of CH3Cl + H2O isolated in solid neon at low temperature have been investigated. High concentration studies of water (0.01%-4%) and subsequent annealing lead to the formation of the ternary CH3Cl:(H2O)2 complex. Detailed vibrational assignments were made on the observed spectra of water and deuterated water engaged in the complex. In parallel, structural, energetic, and vibrational properties of the complex have been studied at the second-order M?ller-Plesset perturbation theory using several basis sets. Anaharmonic correction to the vibrational frequencies has been done with the standard second-order perturbation approach. It was shown that the ground state of the complex has a cyclic form for which the nonadditive three-body contribution was found to be around 10% of the interaction energy.     

Structures, energetics, vibrational spectra of NH4+ (H2O)(n=4,6) clusters: Ab initio calculations and first principles molecular dynamics simulations

Important structural isomers of NH(4) (+)(H(2)O)(n=4,6) have been studied by using density functional theory, Moller-Plesset second order perturbation theory, and coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. The zero-point energy (ZPE) correction to the complete basis set limit of the CCSD(T) binding energies and free energies is necessary to identify the low energy structures for NH(4) (+)(H(2)O)(n=4,6) because otherwise wrong structures could be assigned for the most probable structures. For NH(4) (+)(H(2)O)(6), the cage-type structure, which is more stable than the previously reported open structure before the ZPE correction, turns out to be less stable after the ZPE correction. In first principles Car-Parrinello molecular dynamics simulations around 100 K, the combined power spectrum of three lowest energy isomers of NH(4) (+)(H(2)O)(4) and two lowest energy isomers of NH(4) (+)(H(2)O)(6) explains each experimental IR spectrum.     

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.     

A theoretical ab initio study on the H(2)NO + O(3) reaction

The deviation of the NH(2) pseudo-first-order decay Arrhenius plots of the NH(2) + O(3) reaction at high ozone pressures measured by experimentalists, has been attributed to the regeneration of NH(2) radicals due to the subsequent reactions of the products of this reaction with ozone. Although these products have not yet been characterized experimentally, the radical H(2)NO has been postulated, because it can regenerate NH(2) radicals through the reactions: H(2)NO + O(3) --> NH(2) + O(2) and H(2)NO + O(3) --> HNO + OH + O(2). With the purpose of providing a reasonable explanation from a theoretical point of view to the kinetic observed behaviour of the NH(2) + O(3) system, we have carried ab initio electronic structure calculations on both H(2)NO + O(3) possible reactions. The results obtained in this article, however, predict that of both reactions proposed, only the H(2)NO + O(3) --> NH(2) + O(2) reaction would regenerate indeed NH(2) radicals, explaining thus the deviation of the NH(2) pseudo-first-order decay observed experimentally.     

C4H5N-(H2O)n氢键团簇的多光子电离与从头计算研究

在355、532nm激光波长下用TOF质谱研究了C4H5N-(H2O)n氢键团簇体系的多光子电离。二波长下均得到一系列C4H5N-(H2O)n+及质子化产物C4H5N-(H2O)nH+。355nm下可能存在双光子共振电离过程,使得该波长下吡咯母体及团簇离子信号较532nm有明显增强。从头计算结果表明质子化产物的质子更可能连接于吡咯环的α-C原子,而不是N原子上,即光电离过程诱发了一个簇内的质子转移反应。在532nm下质子化产物的生成主要来自于一个发生于团簇内部的Penning电离或电荷转移过程。团簇的形成对吡咯光解产物的稳定化作用使得团簇系列C4H4N-(H2O)n+出现反常强度变化。     

Photodissociation of (SO2)m(H2O)n clusters employing femtosecond pump-probe spectroscopy

A femtosecond pump-probe technique was employed to study the photodissociation dynamics of (SO2)m(H2O)n clusters in real time for clusters, where m=1, 2 and n as large as 11. The pump (excitation) step occurs through a multiphoton process which populates the dissociative E state as well as a lower-lying bound state of the sulfur dioxide (SO2) chromophore. Dissociation of the SO2 monomer occurs through the E state and the decay is fit to a lifetime of 230 fs. The present study is in agreement with our previous investigations of homogeneous (SO2)m clusters that have shown that cluster formation inhibits the dissociation process owing to a steric effect induced by the cluster environment [K. L. Knappenberger, Jr. and A. W. Castleman, Jr., J. Chem. Phys. 121, 3540 (2004)]. The E state lifetime increases sequentially as a function of cluster size to as much as 668 fs when 11 water molecules solvate the chromophore. We have employed a method to compare the ratio of amplitude coefficients, which reflect a respective component of the mathematical fit, to determine the nature of the wave packet evolution in binary clusters. An increase of this ratio by as much as 440% was observed for large cluster sizes. A preferential ion state charge transfer, rather than dissociation, was observed in binary clusters. The significance of cluster size on evaporation processes has been investigated.     

High-level ab initio calculations on CH+2(2A1) + PO(2II) reactions

We present ab initio calculations carried out in the framework of the G 2 theory on the singlet and triplet potential energy surfaces corresponding to the gas-phase between CH⁺₂ and PO. The global minimum of both potential energy surfaces is a cyclic singlet-state cation. Oxygen attachment of PO to CH⁺₂ in a triplet configuration is accompanied by a P(SINGLEBOND)O bond fission, with the result that the corresponding global minimum is an ion-dipole complex between P⁺(³P) and formaldehyde. This is also consistent with the fact that our results predict the formation of formaldehyde to be highly exothermic, either as a neutral or as radical cation. Both charge-transfer processes yielding CH₂(³B₁) or CH₂(¹A₁) are also exothermic. The formation of other carbon and oxygen containing species are endothermic. © 1996 John Wiley & Sons, Inc.     

Lithium cluster anions: photoelectron spectroscopy and ab initio calculations

Structural and energetic properties of small, deceptively simple anionic clusters of lithium, Li(n)(-), n = 3-7, were determined using a combination of anion photoelectron spectroscopy and ab initio calculations. The most stable isomers of each of these anions, the ones most likely to contribute to the photoelectron spectra, were found using the gradient embedded genetic algorithm program. Subsequently, state-of-the-art ab initio techniques, including time-dependent density functional theory, coupled cluster, and multireference configurational interactions methods, were employed to interpret the experimental spectra.     

Analytical potential from ab initio calculations for the Fe+-H2O and Feo-H2O systems

Analytical intermolecular potentials for the Fe⁺?H₂O and Fe^o?H₂O systems have been determined from ab initio calculations. Interaction energies for a lot of points along the two potential energy surfaces were calculated using Huzinga's MINI ?2 basis set. The results obtained were fitted to an analytical function containing 11 adjustable parameters that we have already used with success for the Fe²⁺?H₂O system. The goodness of the generated intermolecular potentials is discussed.