Insights into dynamics of the S2 state of thiophosgene from ab initio calculations

The S2 potential energy surface for Cl2CS dissociation has been characterized with a combined complete active space self-consistent field and multireference configuration interaction method. The S3/S2 minimum-energy intersection has been determined with the state-averaged complete active space self-consistent field method. The S2 direct dissociation was found to have a barrier of 6.0 kcal/mol, leading to formation of Cl(X2P)+ClCS(A2A") in the excited electronic state. Dynamics of the S2 state of Cl2CS can be summarized as follows: (1) The S2-S0 fluorescence occurs with high quantum yield at low excess energies; (2) Both the S(2) dissociation and the S2-->S3 internal conversion cause the loss of the S2-S0 fluorescence upon photoexcitation at 235-253 nm; (3) The S2-->S3 internal conversion (IC) followed by the direct IC to the ground electronic state results in the fragments produced in the ground state, while the S2 dissociation leads to formation of the fragments in excited electronic states.     

Insights into photodissociation dynamics of benzamide and formanilide from ab initio calculations

In the present study, the five lowest electronic states that control the UV photodissociation of formanilide and benzamide have been characterized using the complete active space self-consistent field theory. The mechanisms for the initial relaxation and subsequent dissociation processes have been determined on the basis of the calculated potential energy surfaces and their intersections. It was found that the S(1)/T(1)/T(2) three-surface intersection plays an important role in the photodissociation processes of benzamide. However, the dissociation behavior of formanilide and benzamide was found to be quite different from that for aliphatic amides. The present study provides several insights into the photodissociation dynamics of formanilide and benzamide.     

The lowest-lying electronic singlet and triplet potential energy surfaces for the HNO-NOH system: energetics, unimolecular rate constants, tunneling and kinetic isotope effects for the isomerization and dissociation reactions

The lowest-lying electronic singlet and triplet potential energy surfaces (PES) for the HNO-NOH system have been investigated employing high level ab initio quantum chemical methods. The reaction energies and barriers have been predicted for two isomerization and four dissociation reactions. Total energies are extrapolated to the complete basis set limit applying focal point analyses. Anharmonic zero-point vibrational energies, diagonal Born-Oppenheimer corrections, relativistic effects, and core correlation corrections are also taken into account. On the singlet PES, the (1)HNO → (1)NOH endothermicity including all corrections is predicted to be 42.23 ± 0.2 kcal mol(-1). For the barrierless decomposition of (1)HNO to H + NO, the dissociation energy is estimated to be 47.48 ± 0.2 kcal mol(-1). For (1)NOH → H + NO, the reaction endothermicity and barrier are 5.25 ± 0.2 and 7.88 ± 0.2 kcal mol(-1). On the triplet PES the reaction energy and barrier including all corrections are predicted to be 7.73 ± 0.2 and 39.31 ± 0.2 kcal mol(-1) for the isomerization reaction (3)HNO → (3)NOH. For the triplet dissociation reaction (to H + NO) the corresponding results are 29.03 ± 0.2 and 32.41 ± 0.2 kcal mol(-1). Analogous results are 21.30 ± 0.2 and 33.67 ± 0.2 kcal mol(-1) for the dissociation reaction of (3)NOH (to H + NO). Unimolecular rate constants for the isomerization and dissociation reactions were obtained utilizing kinetic modeling methods. The tunneling and kinetic isotope effects are also investigated for these reactions. The adiabatic singlet-triplet energy splittings are predicted to be 18.45 ± 0.2 and 16.05 ± 0.2 kcal mol(-1) for HNO and NOH, respectively. Kinetic analyses based on solution of simultaneous first-order ordinary-differential rate equations demonstrate that the singlet NOH molecule will be difficult to prepare at room temperature, while the triplet NOH molecule is viable with respect to isomerization and dissociation reactions up to 400 K. Hence, our theoretical findings clearly explain why (1)NOH has not yet been observed experimentally.     

Study of the C(3P) + OH(X2Pi) --> CO(X1Sigma(g)+) + H(2S) reaction: a fully global ab initio potential energy surface of the X2A' state

The C((3)P) + OH(X (2)Pi) --> CO(X (1)Sigma(g)(+)) + H((2)S) reaction has been investigated by ab initio electronic structure calculations of the X(2)A' state based on the multireference (MR) internally contracted single and double configuration interaction (SDCI) method plus Davidson correction (+Q) using Dunning aug-cc-pVQZ basis sets. In particular, the multireference space is taken to be a complete active space (CAS). Improvement over previously proposed potential energy surfaces for HCO/COH is obtained in the sense that present surface describes also the potential part where the CO interatomic distance is large. A large number of geometries (around 2000) have been calculated and analytically fitted using the reproducing kernel Hilbert space (RKHS) method of Ho and Rabitz both for the two-body and three-body terms following the many-body decomposition of the total electronic energies. Results show that the global reaction is highly exothermic ( approximately 6.4 eV) and barrierless (relative to the reactant channel), while five potential barriers are located on this surface. The three minima and five saddle points observed are characterized and found to be in good agreement with previous work. The three minima correspond to the formation of HCO and COH complexes and to the CO + H products, with the COH complex being a metastable minimum relative to the product channel. The five saddle points correspond to potential barriers for both the dissociation/formation of HCO and COH into/from CO + H, to barriers for the isomerization of HCO into COH and to barriers for the inversion of HCO and COH through their respective linear configuration.     

Approach to the atmospheric chemistry of methyl nitrate and methylperoxy nitrite. Chemical mechanisms of their formation and decomposition reactions in the gas phase

Potential energy surfaces, minimum energy reaction paths, minima, transition states, reaction barriers, and conical intersections for the most important atmospheric reactions of methyl nitrate (CH(3)ONO(2)) and methylperoxy nitrite (C(3)HOONO) on the electronic ground state have been studied (i) with the second-order multiconfigurational perturbation theory (CASPT2) by computation of numerical energy gradients for stationary points and (ii) with the density functional theory (DFT). The proposed mechanism explains the conversion of unreactive alkyl peroxy radicals into alkoxy radicals: CH(3)O(2) + NO <=> CH(3)OONO <=> CH(3)O + NO(2) left arrow over right arrow CH(3)ONO(2). Additionally, several discrepancies found in the comparison of the results obtained from the two employed approaches are analyzed. CASPT2 predicts that all dissociation reactions into radicals occur without an extra exit energy barrier. In contrast, DFT finds transition states for the dissociations of cis- and trans-methylperoxy nitrite into CH(3)O + NO(2). Furthermore, multiconfigurational methods [CASPT2 and complete active space SCF (CAS-SCF)] predict the isomerization of CH3ONO2 to CH3OONO to occur in a two-step mechanism: (i) CH(3)ONO(2) --> CH(3)O + NO(2); and (ii) CH(3)O + NO(2) --> CH(3)OONO. The reason for this has to do with the coupling of the ground electronic state with the first excited state. Therefore, it is demonstrated that DFT methods based on single determinantal wave functions give an incorrect picture of the aforementioned reaction mechanisms.     

Methane dehydrogenation on Cu and Ni surfaces with low and moderate oxygen coverage

First-principles density functional theory calculations are carried out to evaluate energy barriers and mechanisms for the dehydrogenation reactions of CH₄ on clean and oxygen-covered surfaces of Cu (111) and Ni (111) with low and moderate oxygen coverage. In the presence of oxygen, two possible pathways have been evaluated. The more likely pathway, which is further analyzed, is that CH₄ loses an H to the surface O. Results from this pathway agree with previous findings showing that oxygen promotes CH₄ dissociation on Cu (111) and hinders that on Ni (111). In addition, our results show lower energy barriers on Cu with higher oxygen coverages up to 0.38 monolayer. However, such an increase in oxygen coverage did not show any favorable effect for CH₄ dissociation on Ni (111). The findings are analyzed through electronic factors revealed by charge analysis and density of states.     

A CASSCF and CASPT2 study on the excited states of s-trans-formaldazine

The low-lying excited states of s-trans-formaldazine (H2CN-NCH2) have been investigated using the complete active space self-consistent field (CASSCF) and the multiconfigurational second-order perturbation (CASPT2) methods. The vertical excitation energies have been calculated at the state-average CASSCF and multistate CASPT2 levels employing the cc-pVTZ basis set. The photodissociation mechanisms starting from the S1 state have been determined. The lowest energy points along the seams of surface intersections have been located in both the Franck-Condon region and the N-N dissociation pathway in the S1 state. Once the system populates the S1 state, in the viewpoint of energy, the radiationless decay via S1/S0(3) conical intersection followed by the N-N bond fission in the ground-state is more favorable in comparison with the N-N dissociation process in the S1 state. A three-surface crossing region (S1/T1/T2), where the S1, T1, and T2 states intersect, was also found. However, the intersystem crossing process via S1/T1/T2 is not energetically competitive with the internal conversion via S1/S0(3).     

Striving to understand the photophysics and photochemistry of thiophosgene: a combined CASSCF and MR-CI study

The potential energy surfaces for Cl(2)CS dissociation into ClCS + Cl in the five lowest electronic states have been determined with the combined complete active space self-consistent field (CASSCF) and MR-CI method. The wavelength-dependent photodissociation dynamics of Cl(2)CS have been characterized through computed potential energy surfaces, surface crossing points, and CASSCF molecular dynamics calculations. Irradiation of the Cl(2)CS molecules at 360-450 nm does not provide sufficient internal energy to overcome the barrier on S(1) dissociation, and the S(1)/T(2) intersection region is energetically inaccessible at this wavelength region; therefore, S(1) --> T(1) intersystem crossing is the dominant process, which is the main reason S(1)-S(0) fluorescence breaks off at excess energies of 3484-9284 cm(-1). Also, the S(1) --> T(2) intersystem crossing process can take place via the S(1)-T(2) vibronic interaction in this range of excess energies, which is mainly responsible for the quantum beats observed in the S(1) emission. Both S(2) direct dissociation and S(2) --> S(3) internal conversion are responsible for the abrupt breakoff of S(2)-S(0) fluorescence at higher excess energies. S(2) direct dissociation leads to the formation of the fragments of Cl(X(2)P) + ClCS(A(2)A' ') in excited electronic states, while S(2) --> S(3) internal conversion followed by direct internal conversion to the ground electronic state results in the fragments produced in the ground state.     

Accurate ab initio structure, dissociation energy, and vibrational spectroscopy of the F(-)-CH4 anion complex

Accurate equilibrium structure, dissociation energy, global potential energy surface (PES), dipole moment surface (DMS), and the infrared vibrational spectrum in the 0-3000 cm(-1) range of the F(-)-CH4 anion complex have been obtained. The equilibrium electronic structure calculations employed second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster (CC) method up to single, double, triple, and perturbative quadruple excitations using the aug-cc-p(C)VXZ [X = 2(D), 3(T), 4(Q), and 5] correlation-consistent basis sets. The best equilibrium geometry has been obtained at the all-electron CCSD(T)/aug-cc-pCVQZ level of theory. The dissociation energy has been determined based on basis set extrapolation techniques within the focal-point analysis (FPA) approach considering (a) electron correlation beyond the all-electron CCSD(T) level, (b) relativistic effects, (c) diagonal Born-Oppenheimer corrections (DBOC), and (d) variationally computed zero-point vibrational energies. The final D(e) and D0 values are 2398 +/- 12 and 2280 +/- 20 cm(-1), respectively. The global PES and DMS have been computed at the frozen-core CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ levels of theory, respectively. Variational vibrational calculations have been performed for CH4 and F(-)-CH4 employing the vibrational configuration interaction (VCI) method as implemented in Multimode.     

Excited electronic state decomposition of furazan based energetic materials: 3,3'-diamino-4,4'-azoxyfurazan and its model systems, diaminofurazan and furazan

We report the first experimental and theoretical study of gas phase excited electronic state decomposition of a furazan based, high nitrogen content energetic material, 3,3'-diamino-4,4'-azoxyfurazan (DAAF), and its model systems, diaminofurazan (DAF) and furazan (C2H2N2O). DAAF has received major attention as an insensitive high energy explosive; however, the mechanism and dynamics of the decomposition of this material are not clear yet. In order to understand the initial decomposition mechanism of DAAF and those of its model systems, nanosecond energy resolved and femtosecond time resolved spectroscopies and complete active space self-consistent field (CASSCF) calculations have been employed to investigate the excited electronic state decomposition of these materials. The NO molecule is observed as an initial decomposition product from DAAF and its model systems at three UV excitation wavelengths (226, 236, and 248 nm) with a pulse duration of 8 ns. Energies of the three excitation wavelengths coincide with the (0-0), (0-1), and (0-2) vibronic bands of the NO A 2Sigma+<--X 2Pi electronic transition, respectively. A unique excitation wavelength independent dissociation channel is observed for DAAF, which generates the NO product with a rotationally cold (20 K) and a vibrationally hot (1265 K) distribution. On the contrary, excitation wavelength dependent dissociation channels are observed for the model systems, which generate the NO product with both rotationally cold and hot distributions depending on the excitation wavelengths. Potential energy surface calculations at the CASSCF level of theory illustrates that two conical intersections between the excited and ground electronic states are involved in two different excitation wavelength dependent dissociation channels for the model systems. Femtosecond pump-probe experiments at 226 nm reveal that the NO molecule is still the main observed decomposition product from the materials of interest and that the formation dynamics of the NO product is faster than 180 fs. Two additional fragments are observed from furazan with mass of 40 amu (C2H2N) and 28 amu (CH2N) employing femtosecond laser ionization. This observation suggests a five-membered heterocyclic furazan ring opening mechanism with rupture of a CN and a NO bond, yielding NO as a major decomposition product. NH2 is not observed as a secondary decomposition product of DAAF and DAF.     

Theoretical study of the electronic states of Nb4, Nb5 clusters and their anions (Nb4-,Nb5-)

Geometries and energy separations of the various low-lying electronic states of Nb(n) and Nb(n) (-) (n=4,5) clusters with various structural arrangements have been investigated. The complete active space multiconfiguration self-consistent field method followed by multireference singles and doubles configuration interaction (MRSDCI) calculations that included up to 52x10(6) configuration spin functions have been used to compute several electronic states of these clusters. The ground states of both Nb(4) ((1)A('), pyramidal) and Nb(4) (-) ((2)B(3g), rhombus) are low-spin states at the MRSDCI level. The ground state of Nb(5) cluster is a doublet with a distorted trigonal bipyramid (DTB) structure. The anionic cluster of Nb(5) has two competitive ground states with singlet and triplet multiplicities (DTB). The low-lying electronic states of these clusters have been found to be distorted due to Jahn-Teller effect. On the basis of the energy separations of our computed electronic states of Nb(4) and Nb(5), we have assigned the observed photoelectron spectrum of Nb(n) (-) (n=4,5) clusters. We have also compared our MRSDCI results with density functional calculations. The electron affinity, ionization potential, dissociation and atomization energies of Nb(4) and Nb(5) have been calculated and the results have been found to be in excellent agreement with the experiment.     

Relativistic multireference calculation of photodissociation of o-, m-, and p-bromofluorobenzene

Quantum chemical calculations with relativistic effects were performed on the photodissociation of o-, m-, and p-bromofluorobenzene (o-, m-, and p-BrFPh) at 266 nm. The method of multistate second-order multiconfigurational perturbation theory in conjunction with spin-orbit interaction through complete active space state interaction was employed to calculate the potential energy curves for the ground and low-lying excited states of o-, m-, and p-BrFPh along their photodissociation reaction coordinates. The dissociation mechanisms with products of Br((2)P(3∕2)) and Br(?)((2)P(1∕2)) states were clarified with the computed potential energy curves and the surface crossings. The current calculations augmented previous theoretical investigations by including relativistic effects and resolved some differences of experimental assignment regarding the dissociation channels of o-, m-, and p-BrFPh.     

The HCSHSC and HCS+HSC+ systems: molecular properties, isomerization, and energetics

The isomer pair HCSHSC and related cations have been studied by means of a highly accurate level of theory. For all the species investigated the near-equilibrium potential energy surface has been calculated using the coupled cluster method in conjunction with correlation consistent basis sets ranging in size from quadruple to sextuple zeta. After extrapolation to the complete basis set limit, additional corrections due to core-valence correlation and scalar relativistic effects have also been included. Consequently, the molecular and spectroscopic properties as well as the ionization potentials and dissociation energies have been predicted to high accuracy. Isomerization path and energy for both radical and cationic species have also been investigated. Finally, the anharmonic vibrational frequencies have been employed in order to obtain zero-point corrections to ionization potentials, dissociation energies, and isomerization barriers: IP0(HCS) = 7.57(4) eV and IP0(HSC) = 9.00(5) eV; D0(C-H) = 49.29(55) kcal/mol and D0(S-H) = 9.99(37) kcal/mol; deltaE0(HCSHSC) = 39.29(49) kcal/mol, and deltaE0(HCS+HSC+) = 72.24(75) kcal/mol.     

Dissociation dynamics of fluorinated ethene cations: from time bombs on a molecular level to double-regime dissociators

The dissociative photoionization mechanism of internal energy selected C(2)H(3)F(+), 1,1-C(2)H(2)F(2)(+), C(2)HF(3)(+) and C(2)F(4)(+) cations has been studied in the 13-20 eV photon energy range using imaging photoelectron photoion coincidence spectroscopy. Five predominant channels have been found; HF loss, statistical and non-statistical F loss, cleavage of the C-C bond post H or F-atom migration, and cleavage of the C=C bond. By modelling the breakdown diagrams and ion time-of-flight distributions using statistical theory, experimental 0 K appearance energies, E(0), of the daughter ions have been determined. Both C(2)H(3)F(+) and 1,1-C(2)H(2)F(2)(+) are veritable time bombs with respect to dissociation via HF loss, where slow dissociation over a reverse barrier is followed by an explosion with large kinetic energy release. The first dissociative ionization pathway for C(2)HF(3) and C(2)F(4) involves an atom migration across the C=C bond, giving CF-CHF(2)(+) and CF-CF(3)(+), respectively, which then dissociate to form CHF(2)(+), CF(+) and CF(3)(+). The nature of the F-loss pathway has been found to be bimodal for C(2)H(3)F and 1,1-C(2)H(2)F(2), switching from statistical to non-statistical behaviour as the photon energy increases. The dissociative ionization of C(2)F(4) is found to be comprised of two regimes. At low internal energies, CF(+), CF(3)(+) and CF(2)(+) are formed in statistical processes. At high internal energies, a long-lived excited electronic state is formed, which loses an F atom in a non-statistical process and undergoes statistical redistribution of energy among the nuclear degrees of freedom. This is followed by a subsequent dissociation. In other words only the ground electronic state phase space stays inaccessible. The accurate E(0) of CF(3)(+) and CF(+) formation from C(2)F(4) together with the now well established Δ(f)H(o) of C(2)F(4) yield self-consistent enthalpies of formation for the CF(3), CF, CF(3)(+) and CF(+) species.     

First-principles study of C adsorption, O adsorption, and CO dissociation on flat and stepped Ni surfaces

The adsorption of atomic oxygen and carbon was studied with plane wave density functional theory on four Ni surfaces, Ni(110), Ni(111), Ni(210), and Ni(531). Various adsorption sites on these surfaces are examined in order to identify the most favorable adsorption site for each atomic species. The dependence of surface bonding on adsorbate coverage is also investigated. Adsorption energies and structural information are obtained and compared with existing experimental results for Ni(110) and Ni(111). In addition, activation barriers to CO dissociation have been determined on Ni(111) and Ni(531) by locating the transition states for these processes. Our results indicate that the binding energies of C are comparatively stronger on stepped surfaces than on flat surfaces, and the energy barriers associated with CO dissociation strongly favor reactions occurring near surface steps.     

X, A, B, C, and D states of the C6H5F+ ion studied using multiconfiguration wave functions

Electronic states of the C6H5F+ ion have been studied within C2v symmetry by using the complete active space self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) methods in conjunction with an atomic natural orbital basis. Vertical excitation energies (Tv) and relative energies (Tv') at the ground-state geometry of the C6H5F molecule were calculated for 12 states. For the five lowest-lying states, 1(2)B1, 1(2)A2, 2(2)B1, 1(2)B2, and 1(2)A1, geometries and vibrational frequencies were calculated at the CASSCF level, and adiabatic excitation energies (T0) and potential energy curves (PEC) for F-loss dissociations were calculated at the CASPT2//CASSCF level. On the basis of the CASPT2 T0 calculations, we assign the X, A, B, C, and D states of the ion to 1(2)B1, 1(2)A2, 2(2)B1, 1(2)B2, and 1(2)A1, respectively, which supports the suggested assignment of the B state to (2)(2)B1 by Anand et al. based on their experiments. Our CASPT2 Tv and Tv' calculations and our MRCI T0, Tv, and Tv' calculations all indicate that the 2(2)B1 state of C6H5F+ lies below 1(2)B2. By checking the relative energies of the asymptote products and checking the fragmental geometries and the charge and spin density populations in the asymptote products along the CASPT2//CASSCF PECs, we conclude that the 1(2)B1, 1(2)B2, and 1(2)A1 states of C6H5F+ correlate with C6H5+ (1(1)A1) + F (2P) (the first dissociation limit). The energy increases monotonically along the 1(2)B1 PEC, and there are barriers and minima along the 1(2)B2 and 1(2)A1 PECs. The predicted appearance potential value for C6H5+ (1(1)A1) is very close to the average of the experimental values. Our CASPT2//CASSCF PEC calculations have led to the conclusion that the 1(2)A2 state of C6H5F+ correlates with the third dissociation limit of C6H5+ (1(1)A2) + F (2P), and a preliminary discussion is presented.     

Theoretical investigation on o-, m-, and p-chlorotoluene photodissociations at 193 and 266 nm

Photodissociations of the o-, m-, and p-chlorotoluene at 193 and 266 nm were investigated by ab initio calculations with and without spin-orbit interaction. The experimentally observed photodissociation channels were clearly assigned by multistate second order multiconfigurational perturbation theory (MS-CASPT2) calculated potential energy curves. The dissociation products with spin-orbit-coupled states of Cl*(2P1/2) and Cl(2P3/2) were identified by MS-CASPT2 in conjunction with spin-orbit interaction through complete active space state interaction (MS-CASPT2/CASSI-SO) calculations. The effects of methyl rotation and substituent on the photodissociation mechanism were discussed in detail.     

Full six-dimensional nonadiabatic quantum dynamics calculation for the energy pooling reaction O(2)(a (1)Delta)+O(2)(a (1)Delta)-->O(2)(b (1)Sigma)+O(2)(X (3)Sigma)

Six new potential energy surfaces of four singlet states and two triplet states for the title oxygen molecule reaction along with the spin-orbit coupling among them have been constructed from the complete active space second-order perturbation theory with a 6-311+G(d) basis. Accurate integral cross sections are calculated with a full six-dimensional nonadiabatic time-dependent quantum wave packet method. The thermal rate constant based on the integral cross sections agrees well with the result of the experimental measurements, and the intersystem crossing effects are also discussed in this electronic energy-transfer process.     

CASPT2 study of the decomposition of nitrosomethane and its tautomerization reactions in the ground and low-lying excited states

The dissociation reaction of nitrosomethane into methyl and nitric oxide and the tautomerization reactions to formaldehyde oxime, nitrone, and methoxy nitrene have been studied with the second-order multiconfigurational perturbation theory (CASPT2) by the computation of numerical energy gradients. The prevailing reactions in both the ground and the excited states are dissociations. The structures of the ground and excited states are compared with the corresponding complete active space SCF (CAS-SCF) geometries. It is found that changes in the individual bond lengths are rather large (0.01-0.02 A), while the character and energetics of the CASPT2 optimizations remain similar to the CAS-SCF values.