Search

Calibration-quality ab initio adiabatic potential energy surfaces (PES) have been determined for all isotopologues of the molecular ion H3+. The underlying Born-Oppenheimer electronic structure computations used optimized explicitly correlated shifted Gaussian functions. The surfaces include diagonal Born-Oppenheimer corrections computed from the accurate electronic wave functions. A fit to the 41 655 ab initio points is presented which gives a standard deviation better than 0.1 cm-1 when restricted to the points up to 6000 cm-1 above the first dissociation asymptote. Nuclear motion calculations utilizing this PES, called GLH3P, and an exact kinetic energy operator given in orthogonal internal coordinates are presented. The ro-vibrational transition frequencies for H3+, H2D+, and HD2+ are compared with high resolution measurements. The most sophisticated and complete procedure employed to compute ro-vibrational energy levels, which makes explicit allowance for the inclusion of non-adiabatic effects, reproduces all the known ro-vibrational levels of the H3+ isotopologues considered to better than 0.2 cm-1. This represents a significant (order-of-magnitude) improvement compared to previous studies of transitions in the visible. Careful treatment of linear geometries is important for high frequency transitions and leads to new assignments for some of the previously observed lines. Prospects for further investigations of non-adiabatic effects in the H3+ isotopologues are discussed. In short, the paper presents (a) an extremely accurate global potential energy surface of H3+ resulting from high accuracy ab initio computations and global fit, (b) very accurate nuclear motion calculations of all available experimental line data up to 16 000 cm-1, and (c) results suggest that we can predict accurately the lines of H3+ towards dissociation and thus facilitate their experimental observation. [ABSTRACT FROM AUTHOR]

The abstraction and addition reactions of H with trans-N2H2 are studied by high-level ab initio methods and density functional theory. Rate constants were calculated for these two reactions by multistructural variational transition state theory with multidimensional tunneling and including torsional anharmonicity by the multistructural torsion method. Rate constants of the abstraction reaction show large variational effects, that is, the variational transition state yields a smaller rate constant than the conventional transition state; this results from the fact that the variational transition state has a higher zero-point vibrational energy than the conventional transition state. The addition reaction has a classical barrier height that is about 1 kcal/mol lower than that of the abstraction reaction, but the addition rates are lower than the abstraction rates due to vibrational adiabaticity. The calculated branching ratio of abstraction to addition is 3.5 at 200 K and decreases to 1.2 at 1000 K and 1.06 at 1500 K. [ABSTRACT FROM AUTHOR]

Limiting high pressure rate constants for the recombination reaction H+O2→HO2 are modeled between 0 and 5000 K on an ab initio potential. Quantum capture theory is employed for the temperature range from 0 to about 1 K, while classical trajectory calculations are suitable for covering temperatures above about 200 K. The intermediate temperature range is analyzed by adiabatic channel capture theory. The system is characterized by transition-state switching from outer transition states in the long-range-C6/R6 potential to inner transition states in the range of a “shoulder” of the potential. The limiting high pressure rate constants from the trajectory calculations are sufficient for comparison with the experimental data which are available over the range from 300 to 900 K. Specific rate constants k(E,J) for HO2 dissociation are also given and analyzed with respect to internal consistency with capture cross sections. [ABSTRACT FROM AUTHOR]

In open quantum systems, decoherence occurs through interaction of a quantum subsystem with its environment. The computation of expectation values requires a knowledge of the quantum dynamics of operators and sampling from initial states of the density matrix describing the subsystem and bath. We consider situations where the quantum evolution can be approximated by quantum-classical Liouville dynamics and examine the circumstances under which the evolution can be reduced to surface-hopping dynamics, where the evolution consists of trajectory segments exclusively evolving on single adiabatic surfaces, with probabilistic hops between these surfaces. The justification for the reduction depends on the validity of a Markovian approximation on a bath averaged memory kernel that accounts for quantum coherence in the system. We show that such a reduction is often possible when initial sampling is from either the quantum or classical bath initial distributions. If the average is taken only over the quantum dispersion that broadens the classical distribution, then such a reduction is not always possible. [ABSTRACT FROM AUTHOR]

The adiabatic potential energy surfaces for the lowest five electronic states of 3A″ symmetry for the H++O2 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 3A″ and 2 3A″), which adiabatically correlate in the asymptotic limit to H(2S)+O2+(X 2Πg) and H++O2(X 3Σ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 O2 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 O2 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. [ABSTRACT FROM AUTHOR]

The state-to-state dynamics of the H+D2 reaction is studied by the reactant-product decoupling method using the double many-body expansion potential energy surface. Two approaches are compared: one uses only the lowest adiabatic sheet while the other employs both coupled diabatic sheets. Rotational distributions for the reaction H+D2 (υ=0,j=0)→HD(υ′=3,j′)+D are obtained at eight different collision energies between 1.49 and 1.85 eV; no significant difference are found between the two approaches. Initial state-selected total reaction probabilities and integral cross sections are also given for energies ranging from 0.25 up to 2.0 eV with extremely small differences being observed between the two sets of results, thus showing that the nonadiabatic effects in the title reaction are negligible at least for small energies below 2.0 eV. [ABSTRACT FROM AUTHOR]

The Drude model for treating the interaction of excess electrons with polar molecules is extended to calculate continuum functions and to evaluate photodetachment cross sections. The approach is applied to calculate the cross sections for photodetachment of dipole-bound electrons from HCN- and HNC-. In addition, an adiabatic model separating the angular and radial degrees of freedom of the excess electron is introduced and shown to account in a qualitative manner for the cross sections. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

A study of equilibrium states of a thermodynamic system whose evolution is governed not only by the temperature, but also by the ordering field is carried out. It is found that an adiabatically insulated system may have a new type of nonuniform state of equilibrium which is inhomogeneous in temperature. The comparison is made of the stability conditions in isothermal and adiabatic systems. The steady motion of an interface boundary during a first-order phase transition is investigated. It is shown, that depending upon the values of the diffusion coefficients, different regimes can exist. For small thermal diffusivity, the temperature of the final phase after the exothermal transition can be above the equilibrium point. The kinetic problem is reformulated to a dynamical system, and a numerical procedure to solve the latter is presented. Numerical results are discussed in comparison with analytic ones. [ABSTRACT FROM AUTHOR]

One of us has previously developed a technique for obtaining the leading terms in a perturbation expansion of the Born–Oppenheimer energies and derivative couplings for the X3 system near the conical intersection at a C3v configuration. In a preceding paper, these leading terms have been utilized to study various aspects of the adiabatic approximation for both model and real systems of this type. In the present article, this technique is generalized and extended to all orders, yielding rigorous functional forms for both energies and derivative couplings. In particular, the ‘‘nonremovable’’ part of the derivative coupling, which cannot be transformed away by going to a diabatic basis, is explicitly exhibited. Convenient approximations for both removable and nonremovable couplings are also obtained. These should facilitate the estimation of the effect of the different couplings in various situations. [ABSTRACT FROM AUTHOR]