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.     

Potential curves for proton motion in H-bonded systems: HF and H2O polymers

The presence of long range coupling between hydrogen atoms is shown for the HF and H₂O hydrogen bonded systems. The coupling of H atoms critically depends on the spatial orientation of the H atoms being considered. Explicit calculations of the potential curves of the protons are performed using as a model a ring of six HF, or H₂O, molecules. The method of calculation is the CNDO/2. The strong similarities of the results for H₂O and HF polymers supports the conclusion that the coupling is essentially due to factors such as the asymmetric equilibrium position of the H atoms, the high electronic polarizability of the system, etc.     

Theoretical Investigations on the Interaction Nature of Br2 with HF, H2O and NH3

The interactions of HF, H2O and NH3 with Br2 are investigated at the MP2(full)/ aug-cc-pVDZ level. It is found that two kinds of stable complexes, halogen-bonded and hydrogen- bonded complexes, exist between Br2 and HF and between Br2 and H2O. The interaction energy analysis and natural population analysis (NPA) are conducted to these two kinds of complexes, indicating the halogen-bonded complexes are more stable than the corresponding hydrogen-bonded ones, and the binding energies of the former increase in the order HFH2O for the latter.     

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.     

Theoretical study of the structure and vibrational spectra of the (H2O)2…HF and H2O…(HF)2 molecular complexes

Using the Hartree–Fock and MP 2 methods with bases of up to 6-31++G (2d, 2p) quality, the optimum geometry of the 1:2 and 2:1 (H₂O)_n… (HF)_m complexes of water and hydrogen fluoride is searched in a systematic way. Two minimum-energy conformations are found for the 1:2 complex connected through a low-energy transition state. For the 2:1 complex, only one minimum-energy structure is obtained. The analysis of the geometries of the minima and their vibrational frequencies shows that none of them can be used to explain the existence of the H …F? H reverse complex detected experimentally. © 1994 John Wiley & Sons, Inc.     

Theoretical study of formation of ion pairs in (NH3��HCl)(H2O)6 and (NH3��HF)(H2O)6

We have performed theoretical studies on sixteen molecular cubes for both (NH₃·HCl)(H₂O)₆ and (NH₃·HF)(H₂O)₆. We use an empirical gauge, based upon the N?CH and H?CX bond lengths, to categorize the degree to which the cubes are neutral adduct or ion pair in character. On this basis, we describe all sixteen cubes of the former as highly ionized, but only five of the latter as greater than 85% ionic in character. Addition of one or two bridging water molecules to form (NH₃·HF)(H₂O)₇ or (NH₃·HF)(H₂O)₈ raises the percent ionic character to greater than 85% for these systems. The relative energy of the cubes can be categorized based on simple chemical principles. The computed vibrational frequency corresponding to the proton stretch in the N?CH?CF framework shows the highest degree of redshifting for systems near 50% ion-pair character. Molecular cubes close to neutral adduct or to ion-pair character show less redshifting of this vibrational motion.     

Femtosecond dynamics of Cu(H(2)O)(2)

The ultrafast relaxation dynamics of Cu(H(2)O)(2) is investigated using femtosecond photodetachment-photoionization spectroscopy. In addition, stationary points on the Cu(H(2)O)(2) anion, neutral, and cation potential energy surfaces are characterized by ab initio electronic structure calculations. Electron photodetachment from Cu(-)(H(2)O)(2) initiates the dynamics on the ground-state potential energy surface of neutral Cu(H(2)O)(2). The resulting Cu(H(2)O)(2) complexes experience large-amplitude H(2)O reorientation and dissociation. The time evolution of the Cu(H(2)O)(2) fragmentation products is monitored by time-resolved resonant multiphoton ionization. The parent ion, Cu(+)(H(2)O)(2), is not detected above background levels. The rise to a maximum of the Cu(+) signal from Cu(-)(H(2)O)(2), and the decay of the Cu(+)(H(2)O) signal from Cu(-)(H(2)O)(2) have similar tau approximately 10 ps time dependences to the corresponding signals from Cu(-)(H(2)O), but display clear differences at very short and long times. The experimental observations can be understood in terms of the following picture. Prompt dissociation of H(2)O from nascent Cu(H(2)O)(2) gives rise to a vibrationally excited Cu(H(2)O) complex, which dissociates to Cu+H(2)O due to coupling of H(2)O internal rotation to the dissociation coordinate. This prompt dissociation removes all intra-H(2)O vibrational excitation from the intermediate Cu(H(2)O) fragment, which quenches the long time vibrational predissociation to Cu+H(2)O previously observed in analogous experiments on Cu(-)(H(2)O).     

Imaging H2O photofragments in the predissociation of the HCl-H2O hydrogen-bonded dimer

The state-to-state vibrational predissociation (VP) dynamics of the hydrogen-bonded HCl-H(2)O dimer was studied following excitation of the dimer's HCl stretch by detecting the H(2)O fragment. Velocity map imaging (VMI) and resonance-enhanced multiphoton ionization (REMPI) were used to determine pair-correlated product energy distributions. Following vibrational excitation of the HCl stretch of the dimer, H(2)O fragments were detected by 2 + 1 REMPI via the C (1)B(1) (000) ← X (1)A(1) (000) transition. REMPI spectra clearly show H(2)O from dissociation produced in the ground vibrational state. The fragments' center-of-mass (c.m.) translational energy distributions were determined from images of selected rotational states of H(2)O and were converted to rotational state distributions of the HCl cofragment. The distributions were consistent with the previously measured dissociation energy of D(0) = 1334 ± 10 cm(-1) and show a clear preference for rotational levels in the HCl fragment that minimize translational energy release. The usefulness of 2 + 1 REMPI detection of water fragments is discussed.     

Study of the H-F stretching band in the absorption spectrum of (CH3)2O...HF in the gas phase

The absorption spectra of the (CH3)2O...HF complex in the range of 4200-2800 cm(-1) were recorded in the gas phase at a resolutions of 0.1 cm(-1) at T = 190-340 K. The spectra obtained were used to analyze their structure and to determine the temperature dependencies of the first and second spectral moments. The band shape of the (CH3)2O...HF complex in the region of the nu1(HF) stretching mode was reconstructed nonempirically. The nu1 and nu3 stretching vibrations and four bending vibrations responsible for the formation of the band shape were considered. The equilibrium geometry and the 1D-4D potential energy surfaces were calculated at the MP2 6-311++G(2d,2p) level with the basis set superposition error taken into account. On the basis of these surfaces, a number of one- and multidimensional anharmonic vibrational problems were solved by the variational method. Solutions of auxiliary 1D and 2D vibrational problems showed the strong coupling between the modes. The energy levels, transition frequencies and intensities, and the rotational constants for the combining vibrational states necessary to reconstruct the spectrum were obtained from solutions of the 4D problem (nu1, nu3, nu5(B2), nu6(B2)) and the 2D problem (nu5(B1), nu6(B1)). The theoretical spectra reconstructed for different temperatures as a superposition of rovibrational bands associated with the fundamental, hot, sum, and difference transitions reproduce the shape and separate spectral features of the experimental spectra. The calculated value of the nu1 frequency is 3424 cm(-1). Along with the frequencies and absolute intensities, the calculation yields the vibrationally averaged values of the separation between the centers of mass of the monomers Rc.-of-m., R(O...F), and r(HF) for different states. In particular, upon excitation of the nu1 mode, Rc.-of-m. becomes shorter by 0.0861 A, and r(HF) becomes longer by 0.0474 A.     

Reactive scattering dynamics in atom+polyatomic systems: F+C2H6-->HF(v,J)+C2H5

State-to-state scattering dynamics of F+C2H6-->HF(v,J)+C2H5 have been investigated at Ecom=3.2(6) kcalmol under single-collision conditions, via detection of nascent rovibrationally resolved HF(v,J) product states with high-resolution infrared laser absorption methods. State-resolved Doppler absorption profiles are recorded for multiple HF(v,J) transitions originating in the v=0,1,2,3 manifold, analyzed to yield absolute column-integrated densities via known HF transition moments, and converted into nascent probabilities via density-to-flux analysis. The spectral resolution of the probe laser also permits Doppler study of translational energy release into quantum-state-resolved HF fragments, which reveals a remarkable linear correlation between (i) HF(v,J) translational recoil and (ii) the remaining energy available, Eavail=Etot-E(HF(v,J)). The dynamics are interpreted in the context of a simple impulsive model based on conservation of linearangular momentum that yields predictions in good agreement with experiment. Deviations from the model indicate only minor excitation of ethyl vibrations, in contrast with a picture of extensive intramolecular vibrational energy flow but consistent with Franck-Condon excitation of the methylene CH2 bending mode. The results suggest a relatively simple dynamical picture for exothermic atom+polyatomic scattering, i.e., that of early barrier dynamics in atom+diatom systems but modified by impulsive recoil coupling at the transition state between translationalrotational degrees of freedom.     

Model systems for probing metal cation hydration: the V+(H2O) and ArV+(H2O) complexes

In support of mass-selected infrared photodissociation (IRPD) spectroscopy experiments, coupled-cluster methods including all single and double excitations (CCSD) and a perturbative contribution from connected triple excitations [CCSD(T)] have been used to study the V+(H2O) and ArV+(H2O) complexes. Equilibrium geometries, harmonic vibrational frequencies, and dissociation energies were computed for the four lowest-lying quintet states (5A1, 5A2, 5B1, and 5B2), all of which appear within a 6 kcal mol(-1) energy range. Moreover, anharmonic vibrational analyses with complete quartic force fields were executed for the 5A1 states of V+(H2O) and ArV+(H2O). Two different basis sets were used: a Wachters+f V[8s6p4d1f] basis with triple-zeta plus polarization (TZP) for O, H, and Ar; and an Ahlrichs QZVPP V[11s6p5d3f2g] and Ar[9s6p4d2f1g] basis with aug-cc-pVQZ for O and H. The ground state is predicted to be 5A1 for V+(H2O), but argon tagging changes the lowest-lying state to 5B1 for ArV+(H2O). Our computations show an opening of 2 degrees -3 degrees in the equilibrium bond angle of H2O due to its interaction with the metal ion. Zero-point vibrational averaging increases the effective bond angle further by 2.0 degrees -2.5 degrees, mostly because of off-axis motion of the heavy vanadium atom rather than changes in the water bending potential. The total theoretical shift in the bond angle of about +4 degrees is significantly less than the widening near 9 degrees deduced from IRPD experiments. The binding energies (D0) for the successive addition of H2O and Ar to the vanadium cation are 36.2 and 9.4 kcal mol(-1), respectively.     

Association patterns in (HF)(m)(H2O)(n) (m + n = 2-8) clusters

In an attempt to understand the phase behavior of aqueous hydrogen fluoride, the clustering in the mixture is investigated at the molecular level. The study is performed at the mPW1B95/6-31+G(d,p) level of theory. Several previous studies attempted to describe the dissociation of HF in water, but in this investigation, the focus is only on the association patterns that are present in this binary mixture. A total of 214 optimized geometries of (HF)n(H2O)m clusters, with m + n as high as 8, were investigated. For each cluster combination, several different conformations are investigated, and the preferred conformations are presented. Using multiple linear regressions, the average strengths of the four possible H-bonding interactions are obtained. The strongest H-bond interaction is reported to be the H2O...H-F interaction. The most probable distributions of mixed clusters as a function of composition are also deduced. It is found that the larger (HF)n(H2O)m clusters are favored both energetically and entropically compared to the ones that are of size m + n < or = 3. Also, the clusters with equimolar contributions of HF and H2O are found to have the strongest interactions.     

The dynamics of formation of vibrationally excited HF in reactions of O (21 D2) atoms with partially fluorinated alkanes

The nascent vibrational energy distributions of the HF? formed in the reactions of a series of partially fluorinated alkanes (R_FH; R_F = CH₂F, CHF₂, CF₃, C₂F₅, C₃F₇, and C₇F₁₅) with electronically excited oxygen atoms O(2¹D₂) have been determined by measuring the appearance times of stimulated emissions from various vibration–rotation transitions in a grating-tuned optical cavity. The vibrational energy contents of the HF formed in these reactions were found to be considerably greater than statistically expected. These reactions are believed to occur via vibrationally excited short-lived α;-fluorinated alcohols (R_FOH?), formed by insertion of the O(2¹D₂) atoms into C? H bonds. The observation of nonstatistical energy partitioning in the above reactions is in clear contrast to the result obtained from the O(2¹D₂) + CF₃CH₃ reaction that produces the β-fluorinated alcohol CF₃CH₂OH, from which the HF product carries a near statistical vibrational energy distribution. A mechanism for HF? formation in these very exothermic reactions is presented.     

Vibrational deactivation of HF(υ = 1) in the H2O + HF dimer: HF(υ = 1) + H2O(000) → HF(υ = 0) + H2O(001) + ΔE

The vibration-vibration energy transfer in the near-resonant collision HF(υ = 1) + H₂O(000) → HF(υ = 0) + H₂O(001) + ΔE = 205 cm^?1 has been investigated on the basis of the model of the nonrigid H₂O-HF dimer formation for temperatures not greatly higher than room temperature. The energy mismatch ΔE is considered to be removed by the slow translational motion of two molecules in the complex about their equilibrium separation. A strong negative temperature dependence of the energy exchange rate is shown between 300 and 500 K.     

Infrared absorptions of the H(2)O ... H(2) complex trapped in solid neon

When a sample of neon to which have been added less than 1% each of H(2) and H(2)O is deposited at 4.3 K, the infrared spectrum of the resulting solid includes an absorption by the vibrational fundamental of H(2), which is normally infrared inactive. New absorptions are also associated with the vibrational fundamentals of the H(2)O in the sample. Similar results are obtained for deuterium-enriched samples. The new peaks are assigned to the van der Waals complex of H(2)O with H(2). As has been found in earlier theoretical, gas-phase, and solid-state studies of this and closely related systems, the infrared absorptions arise principally from complexes involving ortho-H(2), for which J=1.     

Ab initio direct dynamics trajectory simulation of C2H5F-->C2H4 + HF product energy partitioning

Direct dynamics classical trajectory simulations were performed to study product energy partitioning in C(2)H(5)F-->C(2)H(4)+HF dissociation. The intrinsic reaction coordinate potential energy curve, reaction energetics, and transition state (TS) properties were calculated for this reaction at different levels of electronic structure theory, and MP2/6-31G( *) was chosen as a meaningful and practical method for performing the direct dynamics. The trajectories show that the HF bond, uncoupled from the other degrees of freedom, is formed within the first 10 fs as the system moves from the TS towards products. The populations of the HF vibration states, determined from the simulations, decrease monotonically as found from experiments. However, the simulation's populations for the low and high energy vibration states are larger and smaller, respectively, than the experimental results. The HF rotational temperature found from the simulations is in agreement with experiment. Increasing the TS's excess energy gives higher rotational temperatures for both C(2)H(4) and HF. Energy is partitioned to the products from both the excess energy in the TS and the potential energy release in the exit channel. Partitioning from these two energy sources is distinguished by varying the TS's excess energy. An analysis of the simulation's energy disposal shows that the fractions of the excess energy partitioned to relative translation, C(2)H(4) vibration, C(2)H(4) rotation, HF vibration, and HF rotation, are 0.17, 0.64, 0.076, 0.067, and 0.046, respectively, and are in good agreement with previous simulations on empirical potentials and experiments. The partitioning found for the potential energy release is 81%, <0.05%, 5%, 11%, and 3% to relative translation, C(2)H(4) vibration, C(2)H(4) rotation, HF vibration, and HF rotation. This result is substantially different than the deduction from experiments, which summarizes the partitioning as 20%, 45%, 24%, and <12% to relative translation, C(2)H(4) vibration+rotation, HF vibration, and HF rotation. Possible origins of the difference between the simulations and experiments in the release of the potential energy is discussed.     

Relaxation of H2O from its /04>- vibrational state in collisions with H2O, Ar, H2, N2, and O2

We report rate coefficients at 293 K for the collisional relaxation of H2O molecules from the highly excited /04>(+/-) vibrational states in collisions with H2O, Ar, H2, N2, and O2. In our experiments, the mid R:04(-) state is populated by direct absorption of radiation from a pulsed dye laser tuned to approximately 719 nm. Evolution of the population in the (/04>(+/-)) levels is observed using the combination of a frequency-quadrupled Nd:YAG laser, which selectively photolyses H2O(/04>(+/-)), and a frequency-doubled dye laser, which observes the OH(v=0) produced by photodissociation via laser-induced fluorescence. The delay between the pulse from the pump laser and those from the photolysis and probe lasers was systematically varied to generate kinetic decays. The rate coefficients for relaxation of H2O(/04>(+/-)) obtained from these experiments, in units of cm3 molecule(-1) s(-1), are: k(H2O)=(4.1+/-1.2) x 10(-10), k(Ar)=(4.9+/-1.1) x 10(-12), k(H2)=(6.8+/-1.1) x 10(-12), k(N2)=(7.7+/-1.5) x 10(-12), k(O2)=(6.7+/-1.4) x 10(-12). The implications of these results for our previous reports of rate constants for the removal of H2O molecules in selected vibrational states by collisions with H atoms (P. W. Barnes et al., Faraday Discuss. Chem. Soc. 113, 167 (1999) and P. W. Barnes et al., J. Chem. Phys. 115, 4586 (2001).) are fully discussed.     

Ab initio potential energy surfaces and nonadiabatic collision dynamics in H(+)+O(2) system

The adiabatic potential energy surfaces for the lowest five electronic states of (3)A" symmetry for the H(+)+O(2) collision system have been obtained at the multireference configuration interaction level of accuracy using Dunning's correlation consistent polarized valence triple zeta basis set. The radial nonadiabatic coupling terms and the mixing angle between the lowest two electronic states (1 (3)A" and 2 (3)A"), which adiabatically correlate in the asymptotic limit to H((2)S)+O(2) (+)(X (2)Pi(g)) and H(+)+O(2)(X (3)Sigma(g)(-)), respectively, have been computed using ab initio procedures at the same level of accuracy to yield the corresponding quasidiabatic potential energy matrix. The computed strengths of the vibrational coupling matrix elements reflect the trend observed for inelastic vibrational excitations of O(2) in the experiments at collision energy of 9.5 eV. The quantum dynamics has been preformed on the newly obtained coupled quasidiabatic potential energy surfaces under the vibrational close-coupling rotational infinite-order sudden framework at the experimental collision energy of 9.5 eV. The present theoretical results for vibrational elastic/inelastic excitations of O(2) are in overall good agreement with the available experimental data obtained from the proton energy-loss spectra in molecular beam experiments [F. A. Gianturco et al., J. Phys. B 14, 667 (1981)]. The results for the complementary charge transfer processes are also presented at this collision energy.     

Skin Bond Electron Relaxation Dynamics of Germanium Manipulated by Interactions with H2, O2, H2O,H2O2, HF,and Au

Although germanium performs amazingly well at sites surrounding hetero‐coordinated impurities and under‐coordinated defects or skins with unusual properties, having important impact on electronic and optical devices, understanding the behavior of the local bonds and electrons at such sites remains a great challenge. Here we show that a combination of density functional theory calculations, zone‐resolved X‐ray photoelectron spectroscopy, and bond order length strength correlation mechanism has enabled us to clarify the physical origin of the Ge 3d core‐level shift for the under‐coordinated (111) and (100) skin with and without hetero‐coordinated H₂, O₂, H₂O, H₂O₂, HF impurities. The Ge 3d level shifts from 27.579 (for an isolated atom) by 1.381 to 28.960 eV upon bulk formation. Atomic under‐coordination shifts the binding energy further to 29.823 eV for the (001) and to 29.713 eV for the (111) monolayer skin. Addition of O₂, HF, H₂O, H₂O₂ and Au impurities results in quantum entrapment by different amounts, but H adsorption leads to polarization.     

Dissociation of hydrogen fluoride in HF(H(2)O)(7)

We have previously demonstrated that H-bond arrangement has a significant influence on the energetics, structure and chemistry of water clusters. In this work, the effect of H-bond orientation on the dissociation of hydrogen fluoride with seven water molecules is studied by means of graph theory and high level ab initio methods. It is found that cubic structures of HF(H(2)O)(7) are more stable than structures of other topologies reported in the literature. Electronic calculations on all possible H-bond orientations of cubie-HF(H(2)O)(7) show that ionized structures are energetically more favorable than nonionized ones. This is an indication that seven water molecules might be capable of ionizing hydrogen fluoride.