Ab initio potential energy surfaces and nonadiabatic interactions in the H+ +NO collision system

Ab initio calculations on the H(+)+NO system have been carried out in Jacobi coordinates at the multireference configuration interaction level employing Dunning's correlation-consistent polarized valence triple zeta basis set to analyze the role of low-lying electronic excited states in influencing the collision dynamics relevant to the experimental collision energy range of 9.5-30 eV. The lowest two adiabatic potential energy surfaces, asymptotically correlating to H(+)+NO(X (2)Pi) and H((2)S)+NO(+)(X (1)Sigma(+)), have been obtained. Using ab initio procedures, the (radial) nonadiabatic couplings and the mixing angle between the lowest two electronic states (1 (2)A' and 2 (2)A') have been obtained to yield the corresponding quasidiabatic potential energy matrix. The strengths of the computed vibrational coupling matrix elements reflect a similar trend, as has been observed experimentally in the magnitudes of the state-to-state transition probability for the inelastic vibrational excitations [J. Krutein and F. Linder, J. Chem. Phys. 71, 559 (1979); F. A. Gianturco et al., J. Phys. B 14, 667 (1981)].     

Cross sections and rate constants for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a quasiclassical trajectory method

First quasiclassical trajectory calculations have been carried out for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a recent ab initio potential energy surface for the ground electronic state, X 2A', of HCO/COH. Total and state-specific integral cross sections have been determined for a wide range of collision energies (0.001-1 eV). Then, thermal and state-specific rate constants have been calculated in the 1-500 K temperature range. The thermal rate constant varies from 1.78x10(-10) cm3 s-1 at 1 K down to 5.96x10(-11) cm3 s-1 at 500 K with a maximum value of 3.39x10(-10) cm3 s-1 obtained at 7 K. Cross sections and rate constants are found to be almost independent of the rovibrational state of OH.     

Quantum chemical and statistical rate investigation of the CF2(a3B1)+NO(X2Pi) reaction: a fast chemical quenching process

The reaction of CF2(a3B1) with NO(X2Pi) was theoretically investigated using the B3LYP, MP2, CCSD(T), G2M, CASSCF, and CASPT2 quantum chemical methods with various basis sets including 6-31G(d), 6-311G(d), 6-311+G(3df), cc-pVDZ, and cc-pVTZ. In agreement with the experimental kinetic data, the CF2(a3B1)+NO(X2Pi) reaction is found to proceed via a fast, barrier-free combination. This process, occurring on the doublet potential energy surface, leads to the electronically excited adduct F2C-NO(22A'), which readily undergoes a surface hopping to the 12A' electronic surface, with a Landau-Zener transition probability estimated to be close to 90% per C-N vibration. The metastable adduct F2C-NO(12A') can then either spontaneously decompose into CF2(X1A1)+NO(X2Pi) in a direct chemical quenching mechanism or relax to its ground-state equilibrium structure F2CNO(X2A'). The product distribution resulting from the latter, chemically activated intermediate was evaluated by solution of the master equation (ME), under different reaction conditions, using the exact stochastic simulation method; microcanonical rate constants were computed using Rice-Ramsperger-Kassel-Marcus (RRKM) theory, based on the potential energy surfaces (PESs) constructed using both G2M and CASPT2 methods. The RRKM/ME analysis reveals that the hot F2CNO(X2A') rapidly fragments almost exclusively to the same products as above, CF2(X1A1)+NO(X2Pi), which amounts to an indirect chemical quenching mechanism. The reaction on the quartet PES is unlikely to be significant except at very high temperatures. The high crossing probability (up to 90%) between the two "avoided" doublet PESs points out the inherent difficulty in treating chemically activated reactions with fast-moving nuclei within the Born-Oppenheimer approximation.     

Ab initio potential energy surface and spectrum of the B(3Pi) state of the HeI2 complex

The three-dimensional interaction potential for I2(B 3Pi0u+)+He is computed using accurate ab initio methods and a large basis set. Scalar relativistic effects are accounted for by large-core relativistic pseudopotentials for the iodine atoms. Using multireference configuration interaction calculations with subsequent treatment of spin-orbit coupling, it is shown for linear and perpendicular structures of the complex that the interaction potential for I2(B 3Pi0u+)+He is very well approximated by the average of the 3A' and 3A" interaction potentials obtained without spin-orbit coupling. The three-dimensional 3A' and 3A" interaction potentials are computed at the unrestricted open-shell coupled-cluster level of theory using large basis sets. Bound state calculations based on the averaged surface are carried out and binding energies, vibrationally averaged structures, and frequencies are determined. These results are found to be in excellent accord with recent experimental measurements from laser-induced fluorescence and action spectra of HeI2. Furthermore, in combination with a recent X-state potential, the spectral blueshift is obtained and compared with available experimental values.     

Ab initio study on the decomposition of first excited state HOOO radicals

The HOOO radical plays a crucial role in atmospheric processes involving the OH radical and O(2) molecule. We present an ab initio molecular orbital theory study on the decomposition reaction of the first excited state HOOO((2)A') with respect to OH and O(2). The geometries and harmonic vibrational frequencies of all stationary points are calculated at the CASSCF and MRCI levels of theory in conjunction with the 6-31+G(d,p) basis set. The potential energy profile of the decomposition reaction is studied at the CASSCF/6-31+G(d,p) level of theory, in which the complete valence orbitals and electrons are included in the active space. The energies of the potential energy profile are further refined at the CASPT2 and MRCI levels of the theory. Additionally, we have determined the interesting reaction process: the HOOO((2)A') radical with C(s) symmetry does not dissociate to OH((2)Pi) and O(2)((3)Sigma(-)(g)) directly as this is forbidden by orbital symmetry, but dissociates to OH((2)Pi) and O(2)((3)Sigma(-)(g)) via the change in symmetry from C(s) to C(infinity v) symmetry with a low barrier.     

A study of the reaction of N+ with O2: experimental quantification of NO+(a 3Sigma+) production (298-500 K) and computational study of the overall reaction pathways

The product branching ratios for NO+(X 1Sigma+) and NO+(a 3Sigma+) produced from the reaction of N+ with O2 have been measured at 298 and 500 K in a selected ion flow tube. Approximately 0.5% of the total products are in NO+(a) at both temperatures, despite the fact that the reaction to form NO+(a) is 0.3 eV exothermic. High-level ab initio calculations of the potential energy surfaces for the N+ + O2 reaction show that the reaction from N+(3P) + O2(3Sigma(g)) reactants starts with an efficient early stage charge transfer to the N(2D) + O2+(X 2Pi) channel, which gives rise to the O2+(X 2Pi) product and, at the same time, serves as the starting point for all of the reaction channels leading to NO+ and O+ products. Pathways to produce NO+(a 3Sigma+) are found to be less favorable than pathways leading to the major product NO+(X 1Sigma+). Production of N(2D) has implications for the concentration of NO in the mesosphere.     

Collision energy dependence of the O(1D) + HCl --> OH + Cl(2P) reaction studied by crossed beam scattering and quasiclassical trajectory calculations on ab initio potential energy surfaces

The dynamics of the O(1D) + HCl --> OH + Cl(2P) reaction are investigated by a crossed molecular beam ion-imaging method and quasiclassical trajectory calculations on the three ab initio potential energy surfaces, the ground 1(1)A' and two excited (1(1)A' and 2(1)A') states. The scattering experiment was carried out at collision energies of 4.2, 4.5, and 6.4 kcal/mol. The observed doubly differential cross sections (DCSs) for the Cl(2P) product exhibit almost no collision energy dependence over this inspected energy range. The nearly forward-backward symmetric DCS indicates that the reaction proceeds predominantly on the ground-state potential energy surface at these energies. Variation of the forward-backward asymmetry with collision energy is interpreted using an osculating complex model. Although the potential energy surfaces obtained by CASSCF-MRCI ab initio calculations exhibit relatively low potential barriers of 1.6 and 6.5 kcal/mol for 1(1)A' and 2(1)A', respectively, the dynamics calculations indicate that contributions of these excited states are small at the collision energies lower than 15.0 kcal/mol. Theoretical DCSs calculated for the ground-state reaction pathway agree well with the observed ones. These experimental and theoretical results suggest that the titled reaction at collision energies less than 6.5 kcal/mol is predominantly via the ground electronic state.     

Dynamics of the O(3P) + CH4 → OH + CH3 reaction is similar to that of a triatomic reaction

The O((3)P) + CH(4) reaction has been investigated using the quasi-classical trajectory (QCT) method and an ab initio pseudotriatomic potential energy surface (PES). This has been mainly motivated by very recent experiments which support the reliability of the triatomic modeling even at high collision energy ( = 64 kcal mol(-1)). The QCT results agree rather well with the experiments (translational and angular distributions of products); i.e., the ab initio pseudotriatomic modeling "captures" the essence of the reaction dynamics, although the PES was not optimized for high E(col). Furthermore, similar experiments on the O((3)P) + CD(4) reaction at moderate E(col) (12.49 kcal mol(-1)) have also been of a large interest here and, under these softer reaction conditions, the QCT method leads to results which are almost in quantitative agreement with experiments. The utility of the ab initio pseudotriatomic modeling has also been recognized for other analogous systems (X + CH(4)) but with very different PESs.     

Experimental and theoretical studies of the CN-Ar van der Waals complex

The CN-Ar van der Waals complex has been observed using the B (2)Sigma(+)-X (2)Sigma(+) and A (2)Pi-X (2)Sigma(+) electronic transitions. The spectra yield a dissociation energy of D(0")=102+/-2 cm(-1) and a zero-point rotational constant of B(0")=0.067+/-0.005 cm(-1) for CN(X)-Ar. The dissociation energy for CN(A)-Ar was found to be D(0')=125+/-2 cm(-1). Transitions to vibrationally excited levels of CN(B)-Ar dominated the B-X spectrum, indicative of substantial differences in the intermolecular potential energy surfaces (PESs) for the X and B states. Ab initio PESs were calculated for the X and B states. These were used to predict rovibrational energy levels and van der Waals bond energies (D(0")=115 and D(0')=183 cm(-1)). The results for the X state were in reasonably good agreement with the experimental data. Spectral simulations based on the ab initio potentials yielded qualitative insights concerning the B-X spectrum, but the level of agreement was not sufficient to permit vibronic assignment. Electronic predissociation was observed for both CN(A)-Ar and CN(B)-Ar. The process leading to the production of CN(A,nu=8,9) fragments from the predissociation of CN(B,nu=0)-Ar was characterized using time-resolved fluorescence and optical-optical double resonance measurements.     

Ab initio prediction of the infrared-absorption spectrum of the C2Cl radical

The three lowest (1(2)A', 2(2)A', and 1(2)A") potential-energy surfaces of the C2Cl radical, correlating at linear geometries with 2Sigma+ and 2Pi states, have been studied ab initio using a large basis set and multireference configuration-interaction techniques. The electronic ground state is confirmed to be bent with a very low barrier to linearity, due to the strong nonadiabatic electronic interactions taking place in this system. The rovibronic energy levels of the 12C12C35Cl isotopomer and the absolute absorption intensities at a temperature of 5 K have been calculated, to an upper limit of 2000 cm(-1), using diabatic potential-energy and dipole moment surfaces and a recently developed variational method. The resulting vibronic states arise from a strong mixture of all the three electronic components and their assignments are intrinsically ambiguous.     

Interaction of lead atom with atmospheric hydroxyl radical. An ab initio and density functional theory study of the resulting complexes PbOH and HPbO

The two potential hypersurfaces 2A' (ground state) and 2A" (excited state) have been studied through ab initio and density functional theory (DFT) methods for the Pb(OH) complex. Two processes have been identified. The first one concerns the hydrogen inversion process in the coordination of PbOH and the second one the isomerization of PbOH into HPbO. Eight stationary points have been found; four of them correspond to the stable structures with symmetries PbOH(2A'), PbOH(2A"), HPbO(2A'), and HPbO(2Pi), and four correspond to transition states [TS] with the symmetries 2Pi, 2A', 2Sigma+, and 2A". The hydrogen inversion process in PbOH exhibits the so-called Renner-Teller effect with a rather low barrier, whereas the isomerization process PbOH-->HPbO exhibits a rather high barrier. The energetic, structural, spectroscopy, and thermodynamics results obtained at various levels through, e.g., DFT with BLYP, B3LYP exchange-correlation functionals, coupled clusters methods, namely CCSD (single and double excitations) and CCSD(T) (with triple excitations, by perturbation) are presented for the whole sets of the stationary points and their dissociation products. The relativistic effects, as well as spin-orbit interaction, taken into account in the case of the BLYP exchange-correlation functional, have been estimated and discussed in order to measure their importance in the case of system including heavy metals such as Pb. Reactions of lead (Pb) with oxidizing atmospheric molecules (OH, HO2, O2, and O3) have been studied at various levels of approximation in order to study the possible existence of PbOH in the atmosphere.     

A new reaction path for the C + NO reaction: dynamics on the 4A' potential-energy surface

We present a new reaction path without significant barriers for the C + NO reaction, forming ground state N((4)S) and CO. Electronic structure (CASPT2) calculations have been performed for the two lowest (4)A' states of the CNO system. The lowest of these states shows no significant barriers against reaction in the C + NO or O + CN channels. This surface has been fitted to an analytical function using a many-body expansion. Using this surface, and the previously published (2)A' and (2)A' surfaces [Andersson et al., Phys. Chem. Chem. Phys., 2000, 2, 613; Andersson et al., Chem. Phys., 2000, 259, 99], we have performed quasiclassical trajectory (QCT) calculations, obtaining rate coefficients for the C((3)P) + NO(X(2)Pi) --> CO(X(1)Sigma(+)) + N((4)S,(2)D) and C((3)P) + NO(X(2)Pi) --> O((3)P) + CN(X(2)Sigma(+)) reactions. We have also simulated the crossed molecular beam experiments of Naulin et al. [Chem. Phys., 1991, 153, 519] The inclusion of the (4)A' surface in the QCT calculations gives excellent agreement with experiments. This is the first time an adiabatic pathway from C((3)P) + NO(X(2)Pi) to CO(X(1)Sigma(+))+N((4)S) has been reported.     

Ab initio treatment of the chemical reaction precursor complex Cl(2P)-HF. 1. Three-dimensional diabatic potential energy surfaces

The three adiabatic potential surfaces of the Cl(2P)-HF complex that correlate with the 2P ground state of the Cl atom were calculated with the ab initio RCCSD(T) method (partially spin-restricted coupled cluster theory including single and double excitations and perturbative correction for the triples). With the aid of a geometry-dependent diabatic mixing angle, calculated by the complete active space self-consistent field (CASSCF) and multireference configuration-interaction (MRCI) methods, these adiabatic potential surfaces were converted to a set of four distinct diabatic potential surfaces required to define the full 3 x 3 matrix of diabatic potentials. Each of these diabatic potential surfaces was expanded in terms of the appropriate spherical harmonics in the angle theta between the HF bond axis r and the Cl-HF intermolecular axis R. The dependence of the expansion coefficients on the Cl-HF distance R and the HF bond length r(HF) was fit to an analytic form. The strongest binding occurs for the hydrogen-bonded linear Cl-HF geometry, with D(e) = 676.5 cm(-1) and R(e) = 6.217 a0 when r(HF) = r(e) = 1.7328 a0. This binding energy D(e) depends strongly on r(HF), with larger r(HF) causing stronger binding. An important contribution to the binding energy is provided by the interaction between the quadrupole moment of the Cl(2P) atom and the dipole of HF. In agreement with this electrostatic picture, the ground state of linear Cl-HF is a 2-fold degenerate electronic Pi state. For the linear Cl-FH geometry the states are in opposite order, i.e., the Sigma state is lower in energy than the Pi state. The following paper in this issue describes full three-dimensional computations of the bound states of the Cl-HF complex, based on the ab initio diabatic potentials of this paper.     

Ab initio study of the low lying electronic states of ZnF and ZnF-

Highly correlated ab initio calculations have been performed for an accurate determination of the electronic structure and of the spectroscopy of the low lying electronic states of the ZnF system. Using effective core pseudopotentials and aug-cc-pVQZ basis sets for both atoms, the potential curves, the dipole moment functions, and the transition dipole moments between relevant electronic states have been calculated at the multireference-configuration-interaction level. The spectroscopic constants calculated for the X(2)Sigma(+) ground state are in good agreement with the most recent theoretical and experimental values. It is shown that, besides the X(2)Sigma(+) ground state, the B(2)Sigma(+), the C(2)Pi, and the D(2)Sigma(+) states are bound. The A(2)Pi state, which has been mentioned in previous works, is not bound but its potential presents a shoulder in the Franck-Condon region of the X(2)Sigma(+) ground state. All of the low lying quartet states are found to be repulsive. The absorption transitions from the v=0 level of the X(2)Sigma(+) ground state toward the three bound states have been evaluated and the spectra are presented. The potential energy of the ZnF(-) molecular anion has been determined in the vicinity of its equilibrium geometry and the electronic affinity of ZnF (EA=1.843 eV with the zero energy point correction) has been calculated in agreement with the photoelectron spectroscopy experiments.     

Interaction potential and transport properties of NeO(+)

The results of both experimental and theoretical studies of the mobility of O(+) in Ne are reported. Errors in the experiments have been carefully assessed, allowing the obtained data to serve as stringent tests of the ab initio potentials. These potentials were calculated using the RCCSD(T) method, employing basis sets of quintuple-zeta quality. Curves were calculated for the lowest (4)Sigma(-) state [arising from O(+)((4)S) interacting with Ne] and for the (2)Pi state [arising from O(+)((2)D) interacting with Ne]. Then, the effects of spin-orbit coupling were incorporated by using the Breit-Pauli operator. The resulting ground state (Omega=32) of NeO(+) gives mobility values in good agreement with experiment at all field strengths. Values of spectroscopic quantities for the ground electronic state of NeO(+) are also presented.     

Theoretical study on the electronic states of NaLi

Configuration interaction calculations have been carried out on electronic states of the NaLi molecule and the cation NaLi(+). Potential energy curves are presented for the lowest nine (1)Sigma(+), seven (1)Pi, four (1)Delta, eight (3)Sigma(+), seven (3)Pi, and four (3)Delta states of NaLi as well as for the lowest ten (2)Sigma(+), six (2)Pi, and two (2)Delta states of NaLi(+). The results of the present many-electron configuration interaction calculations on the cation are in support of previous core-polarization effective potential calculations. The present calculations on the NaLi molecule are complementary to previous theoretical work on this system, including recently observed electronic states that had not been calculated previously as well as an investigation of nonadiabatic effects leading to spectral perturbations. Furthermore, ab initio potential energy curves of the neutral and the ground state of the cation are employed to determine quantum defect that may be employed to generate potential energy curves for nd and (n+1)p (for n>3) Rydberg states of NaLi. The present results on the 3 (1)Pi and 4 (1)Pi states are in good agreement with recent experimental work, whereas on the basis of theoretical data, the recently observed state 5 (1)Pi is better described as 6 (1)Pi.     

Photodissociation of vibrationally excited SH and SD radicals at 288 and 291 nm: the S(1D2) channel

Ultraviolet photodissociation of SH (X 2Pi, upsilon"=2-7) and SD (X 2Pi, upsilon"=3-7) has been studied at 288 and 291 nm, using the velocity map imaging technique to probe the angular and speed distributions of the S(1D2) products. Photodissociation cross sections for the A 2Sigma+<--X 2Pi(upsilon") and 2Delta<--X 2Pi(upsilon") transitions have been obtained by ab initio calculations at the CASSCF-MRSDCI/aug-cc-pV5Z level of theory. Both the experimental and theoretical results show that SH/SD photodissociation from X 2Pi (upsilon"相似文献   

Heavier alkali-metal monosulfides (KS, RbS, CsS, and FrS) and their cations

The heavier alkali-metal monosulfides (KS, RbS, CsS, and FrS) have been studied by high-level ab initio calculations. The RCCSD(T) method has been employed, combined with large flexible valence basis sets. All-electron basis sets are used for potassium and sulfur, with effective core potentials being used for the other metals, describing the core electrons. Potential-energy curves are calculated for the lowest two neutral and cationic states: all neutral monosulfide species have a (2)Pi ground state, in contrast with the alkali-metal monoxide species, which undergo a change in the electronic ground state from (2)Pi to (2)Sigma(+) as the group is descended. In the cases of KS, RbS, and CsS, spin-orbit curves are also calculated. We also calculate potential-energy curves for the lowest (3)Sigma(-) and (3)Pi states of the cations. From the potential-energy curves, spectroscopic constants are derived, and for KS the spectroscopic results are compared to experimental spectroscopic values. Ionization energies, dissociation energies, and heats of formation are also calculated; for KS, we explore the effects of relativity and basis set extrapolation on these values.