Theoretical study of anharmonic vibrational spectra of HNO3, HNO3–H2O, HNO4: Fundamental, overtone and combination excitations

Gas phase dynamics, conformational transitions and spectroscopy of charged saccharides: the oxocarbenium ion, protonated anhydrogalactose and protonated methyl galactopyranoside

Protonated intermediates are postulated to be involved in the rate determining step of many sugar reactions. This paper presents a study of protonated sugar species, isolated in the gas phase, using a combination of infrared multiple photon dissociation (IRMPD) spectroscopy, classical ab initio molecular dynamics (AIMD) and quantum mechanical vibrational self-consistent field (VSCF) calculations. It provides a likely identification of the reactive intermediate oxocarbenium ion structure in a d-galactosyl system as well as the saccharide pyrolysis product anhydrogalactose (that suggests oxocarbenium ion stabilization), along with the spectrum of the protonated parent species: methyl d-galactopyranoside-H+. Its vibrational fingerprint indicates intramolecular proton sharing. Classical AIMD simulations for galactosyl oxocarbenium ions, conducted in the temperature range ∼300-350 K (using B3LYP potentials on-the-fly) reveal efficient transitions on the picosecond timescale. Multiple conformers are likely to exist under the experimental conditions and along with static VSCF calculations, they have facilitated the identification of the individual structural motifs of the galactosyl oxocarbenium ion and protonated anhydrogalactose ion conformers that contribute to the observed experimental spectra. These results demonstrate the power of experimental IRMPD spectroscopy combined with dynamics simulations and with computational spectroscopy at the anharmonic level to unravel conformer structures of protonated saccharides, and to provide information on their lifetimes.

Facet shapes and thermo-stabilities of H₂SO₄•HNO₃ hydrates involved in polar stratospheric clouds

The nucleation, ice crystal shapes and thermodynamic stability of polar stratospheric clouds particles are interesting concerns owing to their implication in the ozone layer destruction. Some of these particles are formed by conformers of H2O, HNO3, and H2SO4. We carried out calculations using density functional theory (DFT) to obtain optimized structures. Several stable trimers are achieved -divided in two groups, one with HNO3 moiety, second with H2SO4 moiety- after pre-optimization at B3LYP/6-31G and subsequently optimization at B3LYP/aug-cc-pVTZ level of theory. For both most stable conformers five H2O molecules are added to their optimized trimers to calculate hydrated geometries. The OH stretching harmonic frequencies are provided for all aggregates. The zero-point energy correction (ZEPC), relative electronic energies (∆E), relative reaction Gibbs free energies ∆(∆G)k-relative, and cooling constant (K cooling ) are reported at three temperatures: 188 K, 195 K, and 210 K. Shapes given in our calculations are compared with various experimental shapes as well as comparisons with their thermo-stabilities.

Largely reduced grid densities in a vibrational self-consistent field treatment do not significantly impact the resultingwavenumbers

Especially for larger molecules relevant to life sciences, vibrational self-consistent field (VSCF) calculations can become unmanageably demanding even when only first and second order potential coupling terms are considered. This paper investigates to what extent the grid density of the VSCF’s underlying potential energy surface can be reduced without sacrificing accuracy of the resulting wavenumbers. Including single-mode and pair contributions, a reduction to eight points per mode did not introduce a significant deviation but improved the computational efficiency by a factor of four. A mean unsigned deviation of 1.3% from the experiment could be maintained for the fifteen molecules under investigation and the approach was found to be applicable to rigid, semi-rigid and soft vibrational problems likewise. Deprotonated phosphoserine, stabilized by two intramolecular hydrogen bonds, was investigated as an exemplary application.

Quantum chemical and spectroscopic investigations of 3-methyladenine

FTIR, FT-Raman and UV-Vis spectra of 3-methyladenine have been recorded and investigated using quantum chemical calculations. The molecular geometry and vibrational spectra of 3-methyladenine in the ground state are computed by using HF and DFT methods with 6-311G(d,p) basis set. VSCF, CC-VSCF methods based on 2MR-QFF and PT2 (Barone method) have been utilized for computing anharmonic vibrational frequencies. These methods yield results that are in remarkable agreement with the experimental data. The magnitudes of coupling between pair of modes have been also computed. Vibrational modes are assigned with the help of visual inspection of atomic displacements. The electronic spectra, simulated at TD-B3LYP/6-311++G(d,p) level of theory, are compared to the experiment. The global quantities: electronic chemical potential, electrophilicity index, chemical hardness and softness based on HOMO and LUMO energy eigenvalues are also computed at B3LYP/6-311++G(d,p) level of theory.

Approximate first-principles anharmonic calculations of polyatomic spectra using MP2 and B3LYP potentials: comparisons with experiment

Anharmonic vibrational spectroscopy calculations using MP2 and B3LYP computed potential surfaces are carried out for a series of molecules, and frequencies and intensities are compared with those from experiment. The vibrational self-consistent field with second-order perturbation correction (VSCF-PT2) is used in computing the spectra. The test calculations have been performed for the molecules HNO3, C2H4, C2H4O, H2SO4, CH3COOH, glycine, and alanine. Both MP2 and B3LYP give results in good accord with experimental frequencies, though, on the whole, MP2 gives very slightly better agreement. A statistical analysis of deviations in frequencies from experiment is carried out that gives interesting insights. The most probable percentage deviation from experimental frequencies is about -2% (to the red of the experiment) for B3LYP and +2% (to the blue of the experiment) for MP2. There is a higher probability for relatively large percentage deviations when B3LYP is used. The calculated intensities are also found to be in good accord with experiment, but the percentage deviations are much larger than those for frequencies. The results show that both MP2 and B3LYP potentials, used in VSCF-PT2 calculations, account well for anharmonic effects in the spectroscopy of molecules of the types considered.

Molecular dynamic simulations of OH-stretching overtone induced photodissociation of fluorosulfonic and chlorosulfonic acid

We investigate the OH-stretching overtone dynamics of fluorosulfonic acid (FSO(3)H) and chlorosulfonic acid (ClSO(3)H) using classical trajectory simulations. The initial nuclear coordinates and momenta for the trajectories are sampled from anharmonic correlation-corrected vibrational self consistent field (CC-VSCF) wavefunctions. We consider both OH-stretching overtone states and combination states containing a mix of OH-stretching overtones and SOH-bending or OSOH-torsional modes. Our molecular dynamics simulations confirm that photodissociation of these sulfonic acids to form the corresponding hydrogen halides and sulfur trioxide (HF + SO(3) and HCl + SO(3)), is possible via highly vibrationally excited states on a picosecond timescale. Hydrogen-hopping events are also observed in the trajectories whereby the hydrogen atom of the excited OH group is found to migrate to one of the adjacent S=O groups. The transition states, activation energies and dissociation dynamics of FSO(3)H and ClSO(3)H are found to be similar to those of H(2)SO(4). We suggest that these halogenated sulfonic acids should be suitable proxy molecules for an experimental investigation of the OH-stretching overtone induced photodissociation of H(2)SO(4) thought to be important in Earths' atmosphere.

Spectroscopic identification and stability of the intermediate in the OH + HONO2 reaction

The reaction of nitric acid with the hydroxyl radical influences the residence time of HONO(2) in the lower atmosphere. Prior studies [Brown SS, Burkholder JB, Talukdar RK, Ravishankara AR (2001) J Phys Chem A 105:1605-1614] have revealed unusual kinetic behavior for this reaction, including a negative temperature dependence, a complex pressure dependence, and an overall reaction rate strongly affected by isotopic substitution. This behavior suggested that the reaction occurs through an intermediate, theoretically predicted to be a hydrogen-bonded OH-HONO(2) complex in a six-membered ring-like configuration. In this study, the intermediate is generated directly by the association of photolytically generated OH radicals with HONO(2) and stabilized in a pulsed supersonic expansion. Infrared action spectroscopy is used to identify the intermediate by the OH radical stretch (nu(1)) and OH stretch of nitric acid (nu(2)) in the OH-HONO(2) complex. Two vibrational features are attributed to OH-HONO(2): a rotationally structured nu(1) band at 3516.8 cm(-1) and an extensively broadened nu(2) feature at 3260 cm(-1), both shifted from their respective monomers. These same transitions are identified for OD-DONO(2). Assignments of the features are based on their vibrational frequencies, analysis of rotational band structure, and comparison with complementary high level ab initio calculations. In addition, the OH (v = 0) product state distributions resulting from nu(1) and nu(2) excitation are used to determine the binding energy of OH-HONO(2), D(0) <or= 5.3 kcal x mol(-1), which is in good accord with ab initio predictions.

Dynamics of proton recombination with NO3- anion in water clusters

Recombination events of a proton with NO3- at (H2O)8 clusters are studied by molecular dynamics, using "on-the-fly" reliable ab initio MP2 potentials. The main findings are: (1) the lifetime of the ions is less than 1.2 picoseconds; (2) the recombination step invariably involves H3O+, not H5O2+; and (3) an essentially unique transition-state structure of H3O+/NO3- for recombination is found in all cases. Proton migration involves both H3O+ and H5O2+ species: Grotthuss and other mechanisms contribute.

Dynamics of vibrational overtone excitations of H2SO4, H2SO4-H2O: hydrogen-hopping and photodissociation processes

Photochemical processes of sulfuric acid (H2SO4) and sulfuric acid monohydrate (H2SO4-H2O) following overtone excitation of the OH stretching mode are studied by classical trajectory simulations using the semiempirical PM3 potential suface in "on the fly" calculations. The main results are the following: (1) In the excitation of H2SO4 to the fifth OH-stretch overtone, hopping of the H atom between oxygen atoms is found to take place in 22% of the trajectories, only once during simulations of 400 ps. (2) All the trajectories for H2SO4 show a rapid cis-trans isomerization. (3) The photolysis of H2SO4 into SO3 + H2O takes place in 5% of the trajectories on a time scale of approximately 9 ps. (4) Only low overtone levels of H2SO4-H2O have sufficiently long lifetimes to be spectroscopically relevant. Excitation to these OH stretching overtones is found to result in the dissociation of the cluster. H hopping or dissociation of H2SO4 does not take place.

Vibrational spectroscopy of (SO42−)∙(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<or=5.

Photochemical processes induced by vibrational overtone excitations: dynamics simulations for cis-HONO, trans-HONO, HNO3, and HNO3-H2O

Photochemical processes in HNO3, HNO3-H2O, and cis- and trans-HONO following overtone excitation of the OH stretching mode are studied by classical trajectory simulations. Initial conditions for the trajectories are sampled according to the initially prepared vibrational wave function. Semiempirical potential energy surfaces are used in "on-the-fly" simulations. Several tests indicate at least semiquantitative validity of the potential surfaces employed. A number of interesting new processes and intermediate species are found. The main results include the following: (1) In excitation of HNO3 to the fifth and sixth OH-stretch overtone, hopping of the H atom between the oxygen atoms is found to take place in nearly all trajectories, and can persist for many picoseconds. H-atom hopping events have a higher yield and a faster time scale than the photodissociation of HNO3 into OH and NO2. (2) A fraction of the trajectories for HNO3 show isomerization into HOONO, which in a few cases dissociates into HOO and NO. (3) For high overtone excitation of HONO, isomerization into the weakly bound species HOON is seen in all trajectories, in part of the events as an intermediate step on the way to dissociation into OH + NO. This process has not been reported previously. Well-established processes for HONO, including cis-trans isomerization and H hopping are also observed. (4) Only low overtone levels of HNO3-H2O have sufficiently long liftimes to be spectrocopically relevant. Excitation of these OH stretching overtones is found to result in the dissociation of the cluster H hopping, or dissociation of HNO3 does not take place. The results demonstrate the richness of processes induced by overtone excitation of HNO(x) species, with evidence for new phenomena. Possible relevance of the results to atmospheric processes is discussed.

Theoretical calculations: Can Gibbs free energy for intermolecular complexes be predicted efficiently and accurately?

The theoretical study has been performed to refine the procedure for calculations of Gibbs free energy with a relative accuracy of less than 1 kcal/mol. Three benchmark intermolecular complexes are examined via several quantum-chemical methods, including the second-order Moller-Plesset perturbation (MP2), coupled cluster (CCSD(T)), and density functional (BLYP, B3LYP) theories augmented by Dunnings correlation-consistent basis sets. The effects of electron correlation, basis set size, and anharmonicity are systematically analyzed, and the results are compared with available experimental data. The results of the calculations suggest that experimental accuracy can be reached only by extrapolation of MP2 and CCSD(T) total energies to the complete basis set. The contribution of anharmonicity to the zero point energy and TDeltaSint values is fairly small. The new, economic way to reach chemical accuracy in the calculations of the thermodynamic parameters of intermolecular interactions is proposed. In addition, interaction energy (De) and free energy change (DeltaA) for considered species have been evaluated by Carr-Parrinello molecular dynamics (CPMD) simulations and static BLYP-plane wave calculations. The free energy change along the reaction paths were determined by the thermodynamic integration/"Blue Moon Ensemble" technique. Comparison between obtained values, and available experimental and conventional ab initio results has been made. We found that the accuracy of CPMD simulations is affected by several factors, including statistical uncertainty and convergence of constrained forces (TD integration), and the nature of DFT (density functional theory) functional. The results show that CPMD technique is capable of reproducing interaction and free energy with an accuracy of 1 kcal/mol and 2-3 kcal/mol respectively.

Vibrational structure theory: new vibrational wave function methods for calculation of anharmonic vibrational energies and vibrational contributions to molecular properties

A number of recently developed theoretical methods for the calculation of vibrational energies and wave functions are reviewed. Methods for constructing the appropriate quantum mechanical Hamilton operator are briefly described before reviewing a particular branch of theoretical methods for solving the nuclear Schrödinger equation. The main focus is on wave function methods using the vibrational self-consistent field (VSCF) as starting point, and includes vibrational configuration interaction (VCI), vibrational Møller-Plesset (VMP) theory, and vibrational coupled cluster (VCC) theory. The convergence of the different methods towards the full vibrational configuration interaction (FVCI) result is discussed. Finally, newly developed vibrational response methods for calculation of vibrational contributions to properties, energies, and transition probabilities are discussed.

Exploring the effect of anharmonicity of molecular vibrations on thermodynamic properties

Thermodynamic properties of selected small and medium size molecules were calculated using harmonic and anharmonic vibrational frequencies. Harmonic vibrational frequencies were obtained by normal mode analysis, whereas anharmonic ones were calculated using the vibrational self-consistent field (VSCF) method. The calculated and available experimental thermodynamic data for zero point energy, enthalpy, entropy, and heat capacity are compared. It is found that the anharmonicity and coupling of molecular vibrations can play a significant role in predicting accurate thermodynamic quantities. Limitations of the current VSCF method for low frequency modes have been partially removed by following normal mode displacements in internal, rather than Cartesian, coordinates.

A comparison of experimental and calculated spectra of HNO3 in the near-infrared using Fourier transform infrared spectroscopy and vibrational perturbation theory

This work combines new laboratory studies of the near-infrared vibrational spectra of HNO3 with theoretical predictions of these spectra as a means to understand the properties of this molecule at energies well above the fundamental region. Trends in overtone and combination band energy levels and intensities are compiled and examined. The theoretical calculations are in excellent agreement with the observed frequencies and intensities of the transitions in this spectral region. The calculations also serve as a valuable aid for assigning many of the transitions. This work validates the ab initio generated potential energy surface for HNO3 by comparing vibrational perturbation theory calculations to experimental spectra focused on combination band and overtone absorptions.

New Experimental and Theoretical Approach to the Heterogeneous Hydrolysis of NO2: Key Role of Molecular Nitric Acid and Its Complexes

Although heterogeneous chemistry on surfaces in the troposphere is known to be important, there are currently only a few techniques available for studying the nature of surface-adsorbed species as well as their chemistry and photochemistry under atmospheric conditions of 1 atm pressure and in the presence of water vapor. We report here a new laboratory approach using a combination of long path Fourier transform infrared spectroscopy (FTIR) and attenuated total reflectance (ATR) FTIR that allows the simultaneous observation and measurement of gases and surface species. Theory is used to identify the surface-adsorbed intermediates and products, and to estimate their relative concentrations. At intermediate relative humidities typical of the tropospheric boundary layer, the nitric acid formed during NO2 heterogeneous hydrolysis is shown to exist both as nitrate ions from the dissociation of nitric acid formed on the surface and as molecular nitric acid. In both cases, the ions and HNO3 are complexed to water molecules. Upon pumping, water is selectively removed, shifting the NO(3-)-HNO3(H2O)y equilibria toward more dehydrated forms of HNO3 and ultimately to nitric acid dimers. Irradiation of the nitric acid-water film using 300-400 nm radiation generates gaseous NO, while irradiation at 254 nm generates both NO and HONO, resulting in conversion of surface-adsorbed nitrogen oxides into photochemically active NO(x). These studies suggest that the assumption that deposition or formation of nitric acid provides a permanent removal mechanism from the atmosphere may not be correct. Furthermore, a potential role of surface-adsorbed nitric acid and other species formed during the heterogeneous hydrolysis of NO2 in the oxidation of organics on surfaces, and in the generation of gas-phase HONO on local to global scales, should be considered.

First-principles calculation of geometry and anharmonic vibrational spectra of thioformamide and thioformamide-d2

The equilibrium geometry of thioformamide HCSNH2 has been determined at the MP2 and CCSD(T) electron correlation levels under C(s) symmetry constraints using triple-zeta basis sets up to cc-pVTZ. All optimized planar structures are true minima on the potential-energy surface and are characterized by the C-N bond length within 1.353-1.343 A, C-S distances of 1.656-1.628 A, and NCS angle between 125.7 degrees and 125.9 degrees . The wave number of the NH2 out-of-plane wagging mode computed in the harmonic approximation shows stronger dependence on the basis set rather than the electron correlation level and varies from 85.9 cm(-1) at CCSD(T)cc-pVDZ level to 335 cm(-1) at MP2/aug-cc-pVTZ level. Anharmonic vibrational spectra of HCSNH2 and HCSND2 have been determined directly from the potential-energy surfaces computed at MP2 level in triple-zeta valence (TZV)(2df,2p) and TZV+(2df,2p) basis sets using vibrational self-consistent-field (VSCF) and correlation-corrected VSCF (CC-VSCF) methods. CC-VSCF wave numbers of fundamental, first overtone, and most intense combination transitions are reported for thioformamide and those of fundamentals for thioformamide-d2. The NH2 wagging (nu12) mode is strongly anharmonic and its fundamentals have been computed at 406.9 cm(-1) [TZV(2df,2p)] and 399.5 cm(-1) [TZV+(2df,2p)], which is remarkably close to the experimental energy of 393 cm(-1). Anharmonically computed fundamentals of this mode in thioformamide-d2, 299.7 cm(-1) [TZV(2df,2p)] and 299.6 cm(-1) [TZV+(2df,2p)], are only approximately 7 cm(-1) higher than the transition energy (293 cm(-1)) observed in the gas phase spectrum of HCSND2. The first overtone of the NH2 wagging mode of thioformamide (nu12 (02)) has been calculated by CC-VSCF procedure at 830.8 cm(-1) [TZV(2df,2p)] and 880.0 cm(-1) [TZV+(2df,2p)], which implies "negative" (nu12 (02)>2*nu12 (01)) anharmonicity of this mode.

Ab Initio Vibrational Calculations for H2SO4 and H2SO4·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.