Search

Quasiclassical trajectories have been run for N 2 A 3 Σ u + + H 2 S collisions allowing for electronic transitions through surface hopping. An accurate multi-sheeted representation of the HN2 potential energy surface has been employed. For the rate constant at room temperature, we obtain 9.85 ± 0.16 × 10 - 10 cm 3 s - 1 , a value fivefold larger than the recommended experimental one. Unprecedently, we investigate the temperature dependence of the rate constant, and obtained an essentially flat Arrehenius plot. We estimate that the NH yield is at least one order of magnitude lower than the experimental upper limit.

The multiple energy switching scheme [J. Chem. Phys. 119 (2003) 2596] has been used to improve the double many-body expansion (DMBE II) potential energy surface of C 3 near its linear global minima by morphing it with an accurate Taylor-series expansion [J. Chem. Phys. 144 (2016) 044307]. The final ES form attains the accuracy of the local form in reproducing the rovibrational spectrum of C 3 while keeping unaltered all key attributes of the original DMBE II, namely conical intersection seams and dissociative channels. The ES form is therefore commended for adiabatic spectroscopic and reaction dynamics studies.

A potential energy surface is reported at the complete basis set limit for the ground electronic state of HO 2 + via a four dimensional interpolation in relaxed hyperspherical coordinates. The dissociated products OH + and O 2 + have also been studied at the complete basis set level and subsequently analytically modelled. In the case of O 2 + , the vibrational states predicted from the model potential have also been directly included in the fit. The potential form so obtained matches the RKR levels with a root mean square deviation below 1 cm - 1 . Spin orbit + relativistic effects have also been taken into account for the O 2 + ion.

Despite its paradigmatic role in early molecular beam experiments and theoretical molecular reaction dynamics, an attempt is made in the present work to elucidate basic aspects of the title reaction. Both thermal rate constants and its rotation specificity are studied using two recently reported ab initio potential energy surfaces by running quasiclassical trajectories. Some inconsistencies in previous studies are pinpointed, but the results show that they do not alter significantly the agreement between experiment and theory, thus suggesting that further work at both levels is necessary for a more complete understanding.

We report quantum dynamics calculations of the O + OH → H + O2 reaction on the CHIPR (Varandas, 2013) and DMBE IV (Pastrana et al., 1990) potential energy surfaces (PESs) for ground-state HO2 using the 3D time-dependent wavepacket formalism based on hyperspherical coordinates. Reaction probabilities for J = 0 are calculated for several initial rovibrational states of the OH radical (v = 0; j = 0–5). The J-shifting approximation is used to obtain initial state selected rate-coefficients in the range 0–400 K. Total and state-to-state rate-coefficients are predicted in good agreement with recent experimental studies and theoretical calculations.

A general program to fit global adiabatic potential energy surfaces of up to tetratomic molecules to ab initio points and available spectroscopic data for simple diatomics is reported. It is based on the Combined-Hyperbolic-Inverse-Power-Representation (CHIPR) method. The final form describes all dissociating fragments and long-range/valence interactions, while obeying the system permutational symmetry. The code yields as output a Fortran 90 subroutine that readily evaluates the potential and gradient at any arbitrary geometry. Program summary Program title: CHIPR-4.0 CPC Library link to program files: http://dx.doi.org/10.17632/8wdv87gt5x.2 Licensing provisions: GPLv3 Programming language: Fortran 90 Journal reference of previous version: C. M. R. Rocha, A. J. C. Varandas, Comput. Phys. Commun. 247 (2020) 106913 Does the new version supersede the previous version?: Yes Reasons for the new version: CHIPR-4.0 is a major extension of the previous CHIPR-3.0 code [1] which, besides diatomics and triatomics, also includes the possibility of fitting global adiabatic potential energy surfaces (PESs) of A 4 -, A B 3 -, A 2 B 2 -, A B C 2 - and A B C D -type tetratomic molecules [3]. Additionally, a new feature is that it allows the user to fine-tune ab initio diatomic curves using available spectroscopic data. Summary of revisions: • Implementation of CHIPR’s polynomial form and permutation operators for tetratomic molecules • Implementation of a direct-fit module to fine-tune ab initio potential energy curves using spectroscopic data • Automatic global minimum search on the diatomic, triatomic or tetratomic final PESs • Automatic harmonic vibrational analysis at the optimum molecular geometries Nature of problem: This version of the CHIPR code, CHIPR-4.0, provides a set of subroutines to fit global adiabatic potential energy surfaces of up to tetratomic molecules using ab initio and (for diatomics) optional fine-tuning experimental data. Solution method: CHIPR-4.0 uses the Combined-Hyperbolic-Inverse-Power-Representation (CHIPR) [2] method to interpolate and extrapolate ab initio data points. Additional comments including restrictions and unusual features: For triatomic and tetratomic fits, the user must supply subroutines containing the analytic forms defining the sum of two- and two-plus-three-bodies, respectively. These are utilized internally to calculate (from ab initio data) the corresponding three- and four-body energies to be fitted. The current version of this code is strictly based on single-sheeted analytical PESs. High-accuracy in diatomic refinements is currently limited to the 1 Σ single-curve case. References [1] C. M. R. Rocha, A. J. C. Varandas, Comput. Phys. Commun. 247 (2020) 106913 DOI: 10.1016/j.cpc.2019.106913 [2] A. J. C. Varandas, J. Chem. Phys. 138 (2013) 054120 DOI: 10.1063/1.4788912 [3] C. M. R. Rocha, A. J. C. Varandas, Phys. Chem. Chem. Phys. 21 (2019) 24406 DOI: 10.1039/C9CP04890A

A recent generalization of the Jahn-Teller coordinate proposed for modeling cusps on single-sheeted adiabatic potential energy surfaces for triatomics is extended to any type of crossing seam. The model is applied to the 2 A ″ state of NO 2 and the more complicated HN 2 ( 2 A ′ ) system, which shows permutationally equivalent and strongly curved seams.

Sub-femtosecond nuclear dynamics and high-order harmonic generation (HHG) studies are reported for the X 2 B 1 and A 2 A 1 states of Mu 2 O + . The photoelectron spectra and autocorrelation functions are calculated by solving the time-dependent Schrodinger equation, and the HHG signals from the autocorrelation functions for the two cationic states. Good agreement is observed with our earlier studies, with the autocorrelation function ratios revealing maxima as a function of time. Expectation values of bond lengths and bond angle show quasiperiodic oscillations that reflect repeated passages of the wavepacket at minima of the potential surfaces, thence being responsible for the HHG peaks.

The performance of various hierarchized basis sets has been assessed by extrapolating the correlation energy to the estimated one-electron complete basis set limit. New hierarchical numbers are obtained which are found to lie close to previously reported ones for the Dunning correlation-consistent basis. Universal hierarchical numbers are also suggested, with correlation extrapolations showing a significantly better performance than when based on traditional protocols. It is observed that the accuracy of the extrapolation is enhanced when the underlying plot approaches linearity.

Time-independent quantum reactive scattering studies are performed for the X + O2 → OX + O reactions, where X = Mu, H, D and T are atomic hydrogen isotopes, using the accurate CHIPR potential energy surface (Varandas, 2013) and the popular ABC code. J-shifting thermal rate coefficients have been computed, thence extending our recently reported (Teixidor and Varandas, 2015) values for the reaction H ( 2 S ) + O 2 ( 3 Σ g − ) → OH ( 2 Π ) + O ( 3 P ) . The results are compared with each other and the best experimental and theoretical data available in the literature, with the effect of isotopic substitution on state-specific rate constants being highlighted.

A general program to fit global potential energy surfaces of ABC, AB 2, and A3 triatomic molecules to ab initio points is reported. It is based on the Combined-Hyperbolic-Inverse-Power-Representation (CHIPR) method. The final form describes all dissociating fragments, long-range and valence interactions, while warranting the system permutational symmetry. The code yields as output a Fortran 90 subroutine that readily evaluates the potential and gradient at any arbitrary geometry. Program summary Program Title: CHIPR-3.0 Program Files doi: http://dx.doi.org/10.17632/8wdv87gt5x.1 Licensing provisions: GPLv3 Programming language: Fortran 90 Supplementary material: Directory “EXAMPLES” contains example files and subroutines. Nature of problem: The version zero of the CHIPR code for up to triatomics, CHIPR-3.0, provides a set of subroutines to fit global potential energy surfaces of any triatomic molecule using ab initio data points. Solution method: CHIPR-3.0 uses the Combined-Hyperbolic-Inverse-Power-Representation (CHIPR) method to interpolate and extrapolate ab initio points. Additional comments including restrictions and unusual features: Besides ab initio data, the program generates a subroutine containing two-body terms whose analytic form is subsequently utilized as program input for the triatomic. The first step is dispensable in case the user supplies his own diatomic curves.

A simple scheme is suggested to extrapolate the valence electron correlation to the one-electron complete-basis-set limit in species formed by atoms of hydrogen through neon. The performance of the novel scheme is benchmarked for the correlation energy with a data set of 106 systems and, for the polarizability, with 8 molecules. For the former, the results are excellent, most often better than when extrapolated with the best available two-parameter protocols. For the polarizability, they show an enhancement relative to the raw ab initio values, and good agreement with the experimental ones.

Through the numerical solution of the 3D time-dependent Schrodinger equation, the two-dimensional photoelectron energy spectra (PES) of argon have been calculated, showing good accordance with the experimental measurements. Further analysis on the TDSE calculated data and comparison of PES under linearly and circularly polarized laser pulses demonstrated that the forward rescattering of the recolliding electron is mainly responsible for the appearance of the low energy structure (LES) in the PES. It is further revealed that the LES is correlated with the high angular momenta ( l = 8 – 13 ) of the photoelectron.

We report a theoretical investigation of the minimum energy path for isomerization of HO3 with density functional theory (DFT). Specifically, we search for a functional that can reproduce the energy difference between the cis- and trans-isomers of HO3 which has been accurately determined in previous work. By envisaging a full-dimensional map of the isomerization path, the calculations are restricted to a cost-effective model chemistry with a medium-sized cc-pVTZ basis, with the fraction of exact exchange in one-parameter hybrids used to minimize the differences between the ab initio and DFT calculations.

The stability of dipositronium is discussed. Our purpose is not to compete with highly accurate calculations that are available but point out that its dissociation energy can be understood on the basis of simplified models if satisfying the important requirement of consistency in the description of the ‘atoms’ and molecules of positronium. The importance of correlations in positronium and dipositronium is analyzed. Also studied is the replacement of positrons by heavier Fermions in positronium and dipositronium. It is found that the Born–Oppenheimer approximation is valid provided the mass of the heavier Fermions exceeds 200 times the mass of the electron.

The picture of an insertion reaction with formation of a fully lithiated carbamate is presented inspired in recent experimental information and based on accurate electronic structure calculations. A study is performed of the Li 3 N + CO 2 reaction, which was recently suggested as a candidate to capture carbon dioxide and yield useful chemicals. Investigated in the gas phase, the results predict Li 3 N to leave its planar equilibrium form when capturing CO 2 that bends itself during the capture process. A species with the sum formula Li 3 NCO 2 is formed which, jointly with its fragments, is characterized in detail by Moller–Plesset perturbation theory and the coupled cluster method. The practical interest of the title reaction and work ahead are also assessed.

Full geometry optimizations have been carried out in order to determine the structure of the global minimum for the quintet state of the isolated α -center [FeO] 2+ using the neutral cluster OFe(OH) 2 as a model. The intrinsic reaction coordinate and a potential energy cut were calculated aiming at the examination of other low-energy structures. The absolute minimum here reported for the title species differs from the structures reported in other molecular orbital-based studies. Such a structure has been analyzed in order to determine the oxidation state of the Fe atom in the ferryl catalyst.

Rate coefficients for the N ( 2 D ) + N 2 collisions were calculated employing quasiclassical trajectories and the first available set of potential energy surfaces for such excited nitrogen interactions. The details of the vibrational energy transfer are discussed, such as the contributions from reactive and non-reactive trajectories as well as the contribution of each electronic symmetry. The calculated state-to-state and state-to-all rate coefficients show that deactivation is far more probable than excitation, and multi-quanta deactivation play an important role.

We report a theoretical investigation on the role of water in the reaction between the hydroperoxyl radical and ozone. Due to the encouraging results obtained in our previous work on the HO2 + O3 reaction, and also because of the low computational cost involved, all calculations have employed the Kohn–Sham DFT formalism. In particular, functionals with a high percentage of exact exchange are utilized in order to attain a high accuracy at the saddle points for reaction. It is found that ozone reacts with the HO2 · H2O complex and that the classical barrier height of the oxygen-abstraction mechanism is now lower than the barrier to the hydrogen-abstraction mechanism. The implications of the results on the mechanism of the HO2 + O3 reaction are discussed.

Potential energy curve for the X 2 Σ u + state of H 2 - has been computed at the full CI/cc-pVXZ (X = D − Q) level. An analytic function is fitted to the resulting potential energy curve and is then used to calculate the vibrational bound states of H 2 - and compared with those for H2.

The Longuet–Higgins theorem has recently been extended to Jahn–Teller type electronic degeneracy manifolds of arbitrary dimension and suggested to obtain generalized Born–Oppenheimer equations including the geometrical phase effect [A.J.C. Varandas, Chem. Phys. Lett. 487 (2010) 139]. By replacing the traditional sign-change diagnostic of the adiabatic electronic wave vectors (upon adiabatic transportation in configuration space along a loop that encircles the point of degeneracy) with a sign-change criterion based upon their rotation by π in the N -dimensional electronic wave-vector space, it has been shown that only one pair of adiabatic electronic states changes sign upon encircling the degeneracy in configuration space. Here, we examine in detail the four-fold problem and report relevant dynamics studies. The results are compared with those obtained from the simple Born–Oppenheimer approximation (geometrical phase ignored) as well as two- and three-fold single-surface generalized Born–Oppenheimer equations. A rationalization is provided based on the symmetry-allowed reactive probabilities.

The single-sheeted double many-body expansion potential energy surface for the ground-state of NH 2 is refined to attain near spectroscopic accuracy from a multiproperty fit to accurate ab initio energies and experimental vibrational levels, including the Renner–Teller effect for total angular momentum (excluding spin) of N = 0. Quasiclassical trajectory calculations on both the original and newly reported potential energy surfaces suggest that the dynamical properties of the original form remain essentially unaltered as aimed.

The present work is concerned with the reaction of the hydroperoxyl radical with ozone which is key in the atmosphere. Emphasis is given to the theoretical work developed at the authors’ Group, in perspective with developments from other laboratories. Although written in the form of a short review, the present manuscript contains also novel, yet exploratory, results in particular relating to the role of adding a chaperon water molecule to reactants.

We present quantum mechanical calculations for the title reaction employing an adiabatic time-dependent wavepacket method. The initial state specified rate constants calculated from the fitted integral cross sections (ICSs) cover temperatures between 10 and 4000 K. ICSs obtained from the capture model method are also given. This approximation should give reliable estimates of the exact ICSs if denser partial waves are used for the low-energy region.

Using a highly accurate potential energy surface reported by Poveda and Varandas (2005) [12] for the quartet state of NH 2 , we have applied the quasi-classical trajectory and quantum mechanical methods to calculate reaction probabilities for the title reaction at zero total angular momentum. The total reaction cross section and product rotational alignment have also been calculated for the ground vibrational–rotational state of the reactants. The specific vibrational–rotational rate constant shows an Arrhenius-type behavior over the investigated temperature range of 200–2500 K, in unprecedentedly good agreement with the experimental data available.

The question of how many states are subject to geometrical phase in a N -fold electronic degeneracy is examined for Jahn–Teller systems. It is shown that it can be diagnosed from the sign change of the adiabatic electronic vector states when a rotation operator is applied to them in functional space defined in N -dimensional spherical polar coordinates, rather than by adiabatic transportation along a closed path in the higher-dimensional configuration space. It is also demonstrated that the Longuet-Higgins theorem for twofold conical intersections (namely that the electronic wave function changes sign on completing a closed path around a doubly-degenerate intersection) can be used to rationalize the results, irrespectively of N .

We report a simple, yet reliable, scheme for the direct diabatization of two adiabatic states showing an avoided intersection. The method requires only discrete ab initio points for the adiabatic curves, but can provide also a robust tool for analytic modeling purposes. The approach is illustrated for the dicarbon molecule.

using time-dependent nonadiabatic and adiabatic wave packet methods, we investigate nonadiabatic effects in collisions of d+(h+) with h-2(d-2) for collision energies from 1.7 to 2.5 ev. a new accurate diabatic potential energy surface [l.p. viegas, a. alijah, a.j.c. varandas, j.chem. phys. 126 (2007) 74309] of the h-3(+) system fitted from full multireference configuration interaction calculations with a cc-pv5z basis set is incorporated into the quantum calculations of initial state resolved integral cross sections. comparison with previous calculations and experimental measurements are used to give insight on the nonadiabatic charge transfer process and isotope effect. (c) 2009 elsevier b. v. all rights reserved.

The potential energy surface of exohedral HeC 36 is analyzed in detail and extrapolated to the complete basis set limit using Moller-Plesset perturbation calculations. A limited number of coupled cluster singles and doubles calculations has also been performed. Exohedral complexes above the hexagonal faces toward the major axis are found to be the most stable, but the helium atom can move almost freely from the equator to the poles of C 36 without dissociating the complex due to the disparity of masses. The effective 1D interaction potential has van der Waals attributes in good agreement with those estimated from He 2 and a Girifalco function for the C 36 dimer, supporting three bound states.

Three cuts of the H 2 S ( 1 A ′ ) potential energy surface, which correspond to the dissociation of a single S–H bond [cut (i)], the simultaneous, C 2 v -symmetric, dissociation of both S–H bonds [cut (ii)], and the C 2 v -symmetric minimum energy dissociation pathway leading to H 2 ( X 1 Σ g + ) and S ( 3 p 4 1 D ) [cut (iii)], are examined with the conventional and completely renormalized (CR) coupled-cluster (CC) methods, and the multi-reference configuration interaction approach [MRCI(Q)]. The size extensive CR-CC method with singles, doubles, and non-iterative triples, termed CR-CC(2,3), provides the results of the MRCI(Q) quality for cuts (i) and (iii). To obtain a similar quality for cut (ii), the CR-CC(2,3) energy must be corrected for the effect of quadruply excited clusters.

A recently proposed scheme that enables a potential energy surface corresponding to a calculation with a large target basis set to be obtained from small basis set calculations via scaling of the electron correlation at a single-pivotal geometry is generalized to include an arbitrary number of such geometries. If the correlation is extrapolated to the complete basis set limit at the pivotal geometries and used for the scaling, the method can yield accurate potentials at costs up to factors nearly fifty times cheaper than required otherwise. The approach, free from parameters alien to the ab initio methods, is tested on N2, O2, and F2.

Accurate time-dependent wave packet calculations have been performed for the isotopic N( 2 D) + D 2 reaction on the recent reported 1 2 A″ DMBE (double many body expansion) potential energy surface of the NH 2 system. The total reaction probabilities for total angular momentum J up to 75 have been calculated to get the converged integral cross-section over a collision energy range of 0.0–1.0 eV. Then, the ground-state thermal rate constant has been derived from the integral cross-section and the intermolecular isotope effect is further discussed.

By expressing an unknown state in terms of a complete set, a simple scheme for approximate quantization of the continuous vibrational-rotational energy distributions that are obtained from quasi-classical trajectory calculations is suggested. The problem of zero-point energy leakage is also revisited, and the new method tested on the prototype O + OH and H + D2 reactions. http://www.sciencedirect.com/science/article/B6TFN-4NCJCX1-3/1/cba1ca88eec1685b1502de957a8cc175

We report a quasi-classical trajectory study of the title reaction using a double many-body expansion potential energy surface for ground state HSO 2 . Calculations for specific low ro-vibrational states of the reactants, as well as thermalized reactants over the range of temperatures 1400 ⩽ T /K ⩽ 2200, have been performed. The calculated cross sections indicate a barrier-type mechanism, and nearly thermalized product distributions. Both HOSO and HSO 2 isomers are formed in quite different instances during the collisional process. The rate constant is found to be enhanced by adding internal energy to the reactants, with the thermalized result being in good agreement with existing experimental data.

We report a quasiclassical trajectory study of the title reaction using a recently reported double many-body expansion potential energy surface for ground state HSO2. Two methods that aim to account for zero point energy leakage are used. A capture-type behavior is predicted, with recrossing effects manifesting mainly at high translational energies. Reaction cross sections and a comparison of the calculated rate constant with available literature data are also reported.

It is shown that the molecular potential energy surface corresponding to a high level of ab initio theory can be accurately predicted by performing calculations with smaller basis sets and then scaling the electron correlation at a single point calculated with the larger target basis set.

By direct-fitting the parameters in the physically motivated extended Hartree–Fock approximate correlation energy model to spectroscopic data, potential curves have been obtained for a variety of diatomic systems that mimic the ro-vibrational data often within much less than 1 cm−1. The approach offers a simple, yet accurate, scheme for modeling the two-body terms that are the leading contributions in a cluster expansion of potential energy surfaces of larger dimensionality.

We report accurate close coupling quantum wave packet calculations of the rate constant for the reaction N( 2 D) + H2 ! NH + H using a recently reported single-sheeted double many-body expansion potential energy surface for NH2(1 2 A 00 ) that has been calibrated from high level ab initio MRCI energies. Results from quasi-classical trajectories run on the same potential energy surface are also given, and used to rationalize the quantum ones. The calculated rotationally averaged rate constant is found to be in good agreement with experimental data. A comparison with other theoretical results is also presented. � 2006 Elsevier B.V. All rights reserved.

The mechanism of the H + ClONO2 reaction is examined by performing QCISD calculations at geometries optimized at the MP2 level. Each of the six reaction channels involves stereoisomeric transition states that have identical energy barriers. The lowest energy barrier is 24.2 kcal mol(-1) for the indirect metathetical pathway leading to OH + cis-CIONO, being the corresponding rate constant calculated employing TST theory. The NO2-elimination channel and the indirect metathetical pathway leading to OH + trans-ClONO should compete with each other as they have barriers of 24.8 and 25.1 kcal mol(-1). For Cl-substitution, Cl-abstraction, and N-attack, the barriers are 27.4, 35.1, and 41.3 kcal mol(-1). (c) 2005 Elsevier B.V. All rights reserved.

We investigate the relaxation of highly vibrationally excited HO 2 in collisions with vibrationally cold O 2 . The calculations employ the quasiclassical trajectory method and a global double many-body expansion potential energy surface for HO 4 ( 2 A). Both deactivation and activation processes are observed, with the rate constants being found to be strongly influenced by the donor internal energy content. It is also found that the title process may not be ignorable when modelling the ozone budget in the middle atmosphere.

Standard single-valued double many-body expansion (DMBE) theory has been extended to allow an accurate representation of a double-valued potential energy surface. With this new theory, the degeneracy of the two sheets along the conical intersection line can be guaranteed by construction. The method has been applied to the lowest triplet state of H3+ ( a 3 E ′ ) for which the lowest vibrational levels have also been calculated.

We report a dynamics study of both reactive and non-reactive O + OH( v ′, j ′) collisions under local thermodynamic disequilibrium. Besides the total rate constant for temperatures over the range 110 ⩽ T /K ⩽ 450, vibrationally specific ones are reported for quantum numbers v ′ ⩽ 9. Unknown thus far, such kinetic parameters may be of relevance in modelling the OH nightglow in the middle atmosphere.

The title reaction is studied for the first time using a global double many-body expansion potential energy surface. The major differences between the calculated and experimental results are: (a) the calculated activation energy is larger than the experimental one; (b) the reaction is predicted to occur exclusively via oxygen abstraction from ozone while it has been suggested from experimental work that hydrogen abstraction accounts for 88–95% of the reactive encounters. This discrepancy in mechanism is tentatively attributed to scrambling reactions, with further theoretical and experimental work being required to clarify pending issues.

We have performed extensive calculations of the ro-vibrational states of triplet H 3 + , using the method of hyperspherical harmonics and our recently reported double many-body expansion potential energy surface. The rotational term values of the lowest 19 states are presented here for a total angular momentum of J ⩽10.

A single-valued double many-body expansion potential energy surface (DMBE I) recently obtained for the ground electronic state of the sulfur dioxide molecule by fitting correlated ab initio energies suitably corrected by scaling the dynamical correlation energy is now refined by fitting simultaneously available spectroscopic levels up to 6886 cm −1 above the minimum. The topographical features of the novel potential energy surface (DMBE II) are examined in detail, and the method is emphasized as a robust route to fit together state-of-the-art theoretical calculations and spectroscopic measurements using a single fully dimensional potential form.

We report new potential energy curves for the X1Σ+ and A1Π states of the CO molecule by employing the realistic extended Hartree–Fock approximate correlation energy model for diatomics with inclusion of the united-atom limit (EHFACE2U). Special care has been taken to describe the A1Π potential maximum, since this affects the vibrational wave functions of the highest vibrational levels. The corresponding vibronic spectroscopy and Franck–Condon factors (v″=0, 1) have also been calculated using accurate ab initio transition dipole moments. Comparison has been made with recent theoretical and experimental data.

We discuss the title atmospheric reaction for vibrational excited states of OH over the range 0⩽v⩽9. All calculations have employed the quasi-classical trajectory (QCT) method and a realistic double many-body expansion (DMBE) potential energy surface for HO 4 ( 2 A ) . For v=9, the calculated thermal rate coefficient is found to agree with measurements of OH quenching reported by others, which has been attributed to both chemical reaction and non-reactive collisions. The temperature dependence of the rate constant is described by a simple model for all studied vibrational states of OH.

We report a comparative dynamics study of the title reaction using a realistic double many-body expansion potential energy surface (DMBE I) for ground-state HO3 [Mol. Phys. 91 (1997) 301], and a recently reported DMBE II [Chem. Phys. Lett. 334 (2001) 173] surface obtained by modifying DMBE I to fit extensive ab initio calculations corrected for the complete basis set limit. The comparison provides strong evidence that the title reaction can be a significant source of ozone under stratospheric conditions.

We report a single-valued potential energy surface for HO4(2A) from the double many-body expansion method. All n-body (n=2-4) energy terms are taken from published studies on the relevant fragments, with a five-body energy term of Gaussian form added to mimic the experimental activation energy for the OH(v=0)+O3 reaction. A detailed dynamics study of this reaction is also reported using classical trajectories. Good agreement with existing experimental data is obtained. http://www.sciencedirect.com/science/article/B6TFN-41WBCGK-M/1/f846c63bce20f7a4c43e3cdd2338c852

We discuss the singularities which arise in the Hamiltonian operator at a crossing seam involving two potential energy surfaces of the same global symmetry. The Mead–Truhlar and our own equations are discussed and found to differ from each other, although leading to identical phases up to a constant factor and sign in the vicinity of the crossing seam. Also established are the relations which link the vector and scalar gauge potentials with the mixing angle.