Theoretical study of isomerization and dissociation of acetylene dication in the ground and excited electronic states

Ab initio calculations employing the configuration interaction method including Davidson's corrections for quadruple excitations have been carried out to unravel the dissociation mechanism of acetylene dication in various electronic states and to elucidate ultrafast acetylene-vinylidene isomerization recently observed experimentally. Both in the ground triplet and the lowest singlet electronic states of C2H2(2+) the proton migration barrier is shown to remain high, in the range of 50 kcal/mol. On the other hand, the barrier in the excited 2 3A" and 1 3A' states decreases to about 15 and 34 kcal/mol, respectively, indicating that the ultrafast proton migration is possible in these states, especially, in 2 3A", even at relatively low available vibrational energies. Rice-Ramsperger-Kassel-Marcus calculations of individual reaction-rate constants and product branching ratios indicate that if C2H(2)2+ dissociates from the ground triplet state, the major reaction products should be CCH+(3Sigma-)+H+ followed by CH+(3Pi)+CH+(1Sigma+) and with a minor contribution (approximately 1%) of C2H+(2A1)+C+(2P). In the lowest singlet state, C2H+(2A1)+C+(2P) are the major dissociation products at low available energies when the other channels are closed, whereas at Eint>5 eV, the CCH+(1A')+H+ products have the largest branching ratio, up to 70% and higher, that of CH+(1Sigma+)+CH+(1Sigma+) is in the range of 25%-27%, and the yield of C2H++C+ is only 2%-3%. The calculated product branching ratios at Eint approximately 17 eV are in qualitative agreement with the available experimental data. The appearance thresholds calculated for the CCH++H+, CH++CH+, and C2H++C+ products are 34.25, 35.12, and 34.55 eV. The results of calculations in the presence of strong electric field show that the field can make the vinylidene isomer unstable and the proton elimination spontaneous, but is unlikely to significantly reduce the barrier for the acetylene-vinylidene isomerization and to render the acetylene configuration unstable or metastable with respect to proton migration.     

Photoelectron imaging of I2- at 5.826 eV

We report the anion photoelectron spectrum of I2- taken at 5.826 eV detachment energy using velocity mapped imaging. The photoelectron spectrum exhibits bands resulting from transitions to the bound regions of the X 1Sigmag+(0g+), A' 3Piu(2u), A 3Piu(1u), and B 3Piu(0u+) electronic states as well as bands resulting from transitions to the repulsive regions of several I2 electronic states: the B' 3Piu(0u-), B" 1Piu(1u), 3Pig(2g), a 3Pig(1g), 3Pig(0g-), and C 3Sigmau+(1u) states. We simulate the photoelectron spectrum using literature parameters for the I2- and I2 ground and excited states. The photoelectron spectrum includes bands resulting from transitions to several high-lying excited states of I2 that have not been seen experimentally: 3Pig(0g-), 1Pig3(1g), 1 3Sigmag-3(0g+), and the 1Sigmag-3(0u-) states of I2. Finally, the photoelectron spectrum at 5.826 eV allows for the correction of a previous misassignment for the vertical detachment energy of the I2 B 3Piu(0u+) state.     

Photodissociation of N2O: triplet states and triplet channel

The role of triplet states in the UV photodissociation of N(2)O is investigated by means of quantum mechanical wave packet calculations. Global potential energy surfaces are calculated for the lowest two (3)A' and the lowest two (3)A' states at the multi-reference configuration interaction level of electronic structure theory using the augmented valence quadruple zeta atomic basis set. Because of extremely small transition dipole moments with the ground electronic state, excitation of the triplet states has only a marginal effect on the far red tail of the absorption cross section. The calculations do not show any hint of an increased absorption around 280 nm as claimed by early experimental studies. The peak observed in several electron energy loss spectra at 5.4 eV is unambiguously attributed to the lowest triplet state 1(3)A'. Excitation of the 2(1)A' state and subsequent transition to the repulsive branch of the 2(3)A' state at intermediate NN-O separations, promoted by spin-orbit coupling, is identified as the main pathway to the N(2)((1)Σ(g)(+))+O((3)P) triplet channel. The yield, determined in two-state wave packet calculations employing calculated spin-orbit matrix elements, is 0.002 as compared to 0.005 ± 0.002 measured by Nishida et al. [J. Phys. Chem. A 108, 2451 (2004)].     

HO2自由基电子激发态的理论研究

采用CASPT2/CASSCF方法对HO2自由基进行统计算, 优化了三个电子态的稳定点几何构型, 得到详细的频率数据. 利用垂直激发计算确定了3个里德堡态、11个价电子态的电子结构以及在三种理论水平上(CASSCF, SS-CASPT2和MS-CASPT2)的能量信息. 计算中使用了ANO-L和ANO-L+基组, 验证了已知实验数据的同时, 通过与其它理论计算结果的对比, 揭示了应用弥散轨道系数对于该体系激发态研究的重要性.     

A study of the ground state of manganese dimer using quasidegenerate perturbation theory

We study the electronic structure of the ground state of the manganese dimer using the state-averaged complete active space self-consistent field method, followed by second-order quasidegenerate perturbation theory. Overall potential energy curves are calculated for the 1Sigmag+, 11Sigmau+, and 11Piu states, which are candidates for the ground state. Of these states, the 1Sigmag+ state has the lowest energy and we therefore identify it as the ground state. We find values of 3.29 A, 0.14 eV, and 53.46 cm(-1) for the bond length, dissociation energy, and vibrational frequency, in good agreement with the observed values of 3.4 A, 0.1 eV, and 68.1 cm(-1) in rare-gas matrices. These values show that the manganese dimer is a van der Waals molecule with antiferromagnetic coupling.     

Pulsed field ionization electron spectroscopy and molecular structure of aluminum uracil

Al-uracil (Al-C4H4N2O2) was synthesized in a laser-vaporization supersonic molecular beam source and studied with pulsed field ionization-zero electron kinetic energy (ZEKE) photoelectron spectroscopy and density functional theory (DFT). The DFT calculations predicted several low-energy Al-uracil isomers with Al binding to the diketo, keto-enol, and dienol tautomers of uracil. The ZEKE spectroscopic measurements of Al-uracil determined the ionization energy of 43 064(5) cm-1 [or 5.340(6) eV] and a vibrational mode of 51 cm-1 for the neutral complex and several vibrational modes of 51, 303, 614, and 739 cm-1 for the ionized species. Combination of the ZEEK spectrum with the DFT and Franck-Condon factor calculations determined the preferred isomeric structure and electronic states of the Al-uracil complex. This isomer is formed by Al binding to the O4 atom of the diketo tautomer of uracil and has a planar Cs symmetry. The ground electronic states of the neutral and ionized species are 2A' ' and 1A', respectively. The 2A' ' neutral state has a slightly shorter Al-O4 distance than the 1A' ion state. However, the 1A' ion state has stronger metal-ligand binding compared to the 2A' ' state. The increased Al-O4 distance from the 2A' ' state to the 1A' state is attributed to the loss of the pi binding interaction between Al and O4 in the singlet ion state, whereas the increased metal-ligand binding strength is due to the additional charge-dipole interaction in the ion that surpasses the loss of the pi orbital interaction.     

Interaction of triatomic germanium with lithium atoms: electronic structure and stability of Ge3Lin clusters

Quantum chemical calculations were applied to investigate the electronic structure of mono-, di-, and tri- lithiated triatomic germanium (Ge3Lin) and their cations (n = 0-3). Computations using a multiconfigurational quasi-degenerate perturbation approach (MCQDPT2) based on complete active space CASSCF wavefunctions, MRMP2 and density functional theory reveal that Ge3Li has a 2A' ground state with a doublet-quartet gap of 24 kcal/mol. Ge3Li2 has a singlet ground state with a singlet-triplet (3A' '-1A1) gap of 30 kcal/mol, and Ge3Li3 a doublet ground state with a doublet-quartet (4A' '-2A') separation of 16 kcal/mol. The cation Ge3Li+ has a 1A' ground state, being 18 kcal/mol below the 3A' state. The computed electron affinities for triatomic germanium are EA(1) = 2.2 eV (experimental value is 2.23 eV), EA(2) = -2.5 eV, and EA(3) = -5.9 eV, for Ge3-, Ge32-, and Ge33-, respectively, indicating that only the monoanion is stable with respect to electron detachment, in such a way that Ge3Li is composed of Ge3-Li+ ions. An atoms in molecules (AIM) analysis shows the absence of a Ge-Ge-Li ring critical point in Ge3Li. An electron localization function (ELF) map of Ge3Li supports the view that the Ge-Li bond is predominantly ionic; however, a small covalent character could be anticipated from the Laplacian at the Ge-Li bond critical point. The ionic picture of the Ge-Li bond is further supported by the natural bond orbital (NBO) results. The calculated Li affinity value for Ge3 is 2.17 eV, and the Li+ cation affinity value for Ge3- amounts to 5.43 eV. The larger Li+ cation affinity of Ge3- favors an electron transfer, resulting in a Ge3-Li+ interaction.     

Excited state electronic structures and dynamics of NOCl: a new potential function set, absorption spectrum, and photodissociation mechanism

A set of analytical potential energy surfaces (PESs) for six singlet excited states of NOCl are constructed based on multireference configuration interaction calculations. The total absorption cross section at the energy range of 2-7 eV is calculated by quantum dynamics calculations with the present PESs and transition dipole moments. The calculated absorption spectrum agrees well with the experiment. It is also found that the A band with the absorption maximum at 6.3 eV is attributed to the transition to the 4 1A' state, though the excitations to the 3 1A' and 3 1A" states contribute to the spectrum at the energy range between 4 and 5 eV. The spin-forbidden transitions are concluded to be negligibly weak. The mechanism of photodissociation reaction at the energy region corresponding to the A band is examined. The nonadiabatic transition rates from the 4 1A' state to lower singlet and triplet states are estimated by Fermi's golden rule, and the transitions to the 1 1A' and 3 1A' states induced by vibronic coupling are found to be the predominant dissociation pathways. The experimentally observed energy dependence of the recoil anisotropy of the fragments is discussed based on the calculated nonadiabatic transition rates.     

Quantum dynamical study of the O(1D)+HCl reaction employing three electronic state potential energy surfaces

Quantum dynamical calculations are reported for the title reaction, for both product arrangement channels and using potential energy surfaces corresponding to the three electronic states, 1 1A', 2 1A', and 1 1A", which correlate with both reactants and products. The calculations have been performed for J=0 using the time-dependent real wavepacket approach by Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. Reaction probabilities for both product arrangement channels on all three potential energy surfaces are presented for total energies between 0.1 and 1.1 eV. Product vibrational state distributions at two total energies, 0.522 and 0.722 eV, are also presented for both channels and all three electronic states. Product rotational quantum state distributions are presented for both product arrangement channels and all three electronic states for the first six product vibrational states.     

Photoelectron spectroscopic study of the oxyallyl diradical

The photoelectron spectrum of the oxyallyl (OXA) radical anion has been measured. The radical anion has been generated in the reaction of the atomic oxygen radical anion (O(?-)) with acetone. Three low-lying electronic states of OXA have been observed in the spectrum. Electronic structure calculations have been performed for the triplet states ((3)B(2) and (3)B(1)) of OXA and the ground doublet state ((2)A(2)) of the radical anion using density functional theory (DFT). Spectral simulations have been carried out for the triplet states based on the results of the DFT calculations. The simulation identifies a vibrational progression of the CCC bending mode of the (3)B(2) state of OXA in the lower electron binding energy (eBE) portion of the spectrum. On top of the (3)B(2) feature, however, the experimental spectrum exhibits additional photoelectron peaks whose angular distribution is distinct from that for the vibronic peaks of the (3)B(2) state. Complete active space self-consistent field (CASSCF) method and second-order perturbation theory based on the CASSCF wave function (CASPT2) have been employed to study the lowest singlet state ((1)A(1)) of OXA. The simulation based on the results of these electronic structure calculations establishes that the overlapping peaks represent the vibrational ground level of the (1)A(1) state and its vibrational progression of the CO stretching mode. The (1)A(1) state is the lowest electronic state of OXA, and the electron affinity (EA) of OXA is 1.940 ± 0.010 eV. The (3)B(2) state is the first excited state with an electronic term energy of 55 ± 2 meV. The widths of the vibronic peaks of the X? (1)A(1) state are much broader than those of the a? (3)B(2) state, implying that the (1)A(1) state is indeed a transition state. The CASSCF and CASPT2 calculations suggest that the (1)A(1) state is at a potential maximum along the nuclear coordinate representing disrotatory motion of the two methylene groups, which leads to three-membered-ring formation, i.e., cyclopropanone. The simulation of b? (3)B(1) OXA reproduces the higher eBE portion of the spectrum very well. The term energy of the (3)B(1) state is 0.883 ± 0.012 eV. Photoelectron spectroscopic measurements have also been conducted for the other ion products of the O(?-) reaction with acetone. The photoelectron imaging spectrum of the acetylcarbene (AC) radical anion exhibits a broad, structureless feature, which is assigned to the X? (3)A' state of AC. The ground ((2)A') and first excited ((2)A') states of the 1-methylvinoxy (1-MVO) radical have been observed in the photoelectron spectrum of the 1-MVO ion, and their vibronic structure has been analyzed.     

Density Functional Theory (DFT) studies on the ground state of NO_3(~2A''''_2) radical and the first triplet state of NO_3~ cation

Density Functional Theory (DFT) studies on the ground states (2A'2) of NO3 radical and on the ground state (1A1') and the first triplet state (3E") of NO3 cation provide an unambiguous prediction about their geometrical structure-, the ground states of both NO3 radical and NO3 cation have D3h symmetry and the geometrical configuration of the first triplet state 3E" of NO3 cation has C2v symmetry. It is shown that as far as the ionization energy calculations on NO, radical are concerned, the results are only slightly different, no mater that gradient corrections of the exchange-correlation energy are included during self-consistent iterations or they are included as perturbations after the self-consistent iterations.     

Theoretical study of low-lying triplet states of aniline

Multireference complete active space self-consistent-field CASSCF(10,12)/ANO and second-order perturbation theory MS-CASPT2 calculations were performed to determine the vertical low-lying singlet and triplet states of aniline. The sequence of the seven lower lying triplet states is T1(1(3)A'), T2(1(3)A' '), T3(2(3)A'), T4(3(3)A'), T5(2(3)A' '), T6(4(3)A'), and T7(3(3)A' '). The 3(3)A', 4(3)A', and 3(3)A' ' states are assigned as 3s, 3py, and 3pz Rydberg states, respectively, while other states correspond to pi <-- pi excitations. Both the T1 and T2 states are found to be below at the lowest-lying singlet S1 (1(1)A' ') state. Geometry, vibrational modes, and electron distribution of the lowest lying T1 state were determined using UB3LYP calculations. The vertical and adiabatic singlet-triplet energy gaps DeltaE(S0-T1) amount to 3.7 and 3.5 +/- 0.2 eV, respectively. In clear contrast with the S0 state, the triplet aniline is no longer aromatic, and its protonation occurs preferentially at the ring meta-carbon site, with a proton affinity PA = 243 +/- 3 kcal/mol.     

Spiro[4.4]nonatetraene and its Positive and Negative Radical Ions: Molectronic Structure Investigations

The electronic structure of spiro[4.4]nonatetraene 1 as well as that of its radical anion and cation were studied by different spectroscopies. The electron‐energy‐loss spectrum in the gas phase revealed the lowest triplet state at 2.98 eV and a group of three overlapping triplet states in the 4.5 – 5.0 eV range, as well as a number of valence and Rydberg singlet excited states. Electron‐impact excitation functions of pure vibrational and triplet states identified various states of the negative ion, in particular the ground state with an attachment energy of 0.8 eV, an excited state corresponding to a temporary electron attachment to the 2b₁ MO at an attachment energy of 2.7 eV, and a core excited state at 4.0 eV. Electronic‐absorption spectroscopy in cryogenic matrices revealed several states of the positive ion, in particular a richly structured first band at 1.27 eV, and the first electronic transition of the radical anion. Vibrations of the ground state of the cation were probed by IR spectroscopy in a cryogenic matrix. The results are discussed on the basis of density‐functional and CASSCF/CASPT2 quantum‐chemical calculations. In their various forms, the calculations successfully rationalized the triplet and the singlet (valence and Rydberg) excitation energies of the neutral molecule, the excitation energies of the radical cation, its IR spectrum, the vibrations excited in the first electronic absorption band, and the energies of the ground and the first excited states of the anion. The difference of the anion excitation energies in the gas and condensed phases was rationalized by a calculation of the Jahn‐Teller distortion of the anion ground state. Contrary to expectations based on a single‐configuration model for the electronic states of 1 , it is found that the gap between the first two excited states is different in the singlet and the triplet manifold. This finding can be traced to the different importance of configuration interaction in the two multiplicity manifolds.     

Vibrationally resolved photoelectron spectroscopy of BO- and BO2-: a joint experimental and theoretical study

We report a photoelectron spectroscopy and computational study of two simple boron oxide species: BO- and BO2-. Vibrationally resolved photoelectron spectra are obtained at several photon energies (355, 266, 193, and 157 nm) for the 10B isotopomers, 10BO- and 10BO2-. In the spectra of 10BO-, we observe transitions to the 2Sigma+ ground state and the 2Pi excited state of 10BO at an excitation energy of 2.96 eV. The electron affinity of 10BO is measured to be 2.510+/-0.015 eV. The vibrational frequencies of the ground states of 10BO- and 10BO and the 2Pi excited state are measured to be 1725+/-40, 1935+/-30, and 1320+/-40 cm-1, respectively. For 10BO2-, we observe transitions to the 2Pig ground state and two excited states of 10BO2, 2Piu, and 2Sigmau+, at excitation energies of 2.26 and 3.04 eV, respectively. The electron affinity of 10BO2 is measured to be 4.46+/-0.03 eV and the symmetrical stretching vibrational frequency of the 2Piu excited state of 10BO2 is measured to be 980+/-30 cm-1. Both density functional and ab initio calculations are performed to elucidate the electronic structure and chemical bonding of the two boron oxide molecules. Comparisons with the isoelectronic AlO- and AlO2- species and the closely related molecules CO, N2, CN-, and CO2 are also discussed.     

Low-energy states of manganese-oxo corrole and corrolazine: multiconfiguration reference ab initio calculations

Manganese(V)-oxo corrole and corrolazine have been studied with ab initio multiconfiguration reference methods (CASPT2 and RASPT2) and large atomic natural orbital (ANO) basis sets. The calculations confirm the expected singlet d(δ)(2) ground states for both complexes and rule out excited states within 0.5 eV of the ground states. The lowest excited states are a pair of Mn(V) triplet states with d(δ)(1)d(π)(1) configurations 0.5-0.75 eV above the ground state. Manganese(IV)-oxo macrocycle radical states are much higher in energy, ≥1.0 eV relative to the ground state. The macrocyclic ligands in the ground states of the complexes are thus unambiguously 'innocent'. The approximate similarity of the spin state energetics of the corrole and corrolazine complexes suggests that the latter macrocycle on its own does not afford any special stabilization for the Mn(V)O center. The remarkable stability of an Mn(V)O octaarylcorrolazine thus appears to be ascribable to the steric protection afforded by the β-aryl groups.     

Calculating energy levels of isomerizing tetra-atomic molecules. II. The vibrational states of acetylene and vinylidene

A general, full-dimensional computational method for the accurate calculation of rotationally and vibrationally excited states of tetra-atomic molecules is further developed. The resulting computer program may be run in serial and parallel modes and is particularly appropriate for molecules executing wide-amplitude motions and isomerizations. An application to the isomerizing acetylene/vinylidene system is presented. Large-scale calculations using a coordinate system based on orthogonal satellite vectors have been performed in six dimensions and vibrational term values and wave functions for acetylene and vinylidene states up to approximately 23 000 cm(-1) above the potential minimum have been determined. This has permitted the characterization of acetylene and vinylidene states at and above the isomerization barrier. These calculations employ more extensive vibrational basis sets and hence consider a much higher density of states than in any variational calculations reported hitherto for this system. Comparison of the calculated density of states with that determined empirically suggests that our calculations are the most realistic achieved for this system to date. Indeed more states have been converged than in any previous study of this system. Calculations on lower lying excited states of acetylene based on HC-CH diatom-diatom coordinates give nearly identical results to those based on orthogonal satellite vectors. Comparisons are also made with calculations based on HH-CC diatom-diatom coordinates.     

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.     

Photofragmentation of dichloromethanol Cl2CHOH along C-O and C-Cl cleavages: a theoretical study

Multireference configuration interaction calculations are carried out for ground and excited states of dichloromethanol, Cl2CHOH, to investigate two important photofragmentation processes relevant to atmospheric chemistry. Five low-lying excited states (1(1)A", 2(1)A', 1(3)A", 2(3)A" and 1(3)A') in the energy range between 6.4 and 7.5 eV are found to be highly repulsive for C-Cl elongation, leading to ClCHOH (X2A) and Cl (X2P). Photodissociation along the C-O bond resulting in CHCl2 (X2A') and OH (X2II) has to overcome a barrier of about 0.5 eV because the low-lying excited states 1(1)A", 1(3)A' and 1(3)A" become repulsive only after the C-O bond is elongated by about 0.3 A.     

F-loss and H-loss dissociations in low-lying electronic states of the CH3F+ ion studied using multiconfiguration second-order perturbation theory

Complete active space self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) calculations with an atomic natural orbital basis were performed for the 1(2)A', 1(2)A', 2(2)A', 2(2)A', and 3(2)A' (X2E, A2A1, and B2E) states of the CH3F+ ion. The 1(2)A' state is predicted to be the ground state, and the C(s)-state energy levels are different from those of the CH3Cl+ ion. The 2(2)A' (A2A1) state is predicted to be repulsive, and the calculated adiabatic excitation energies for 2(2)A' and 3(2)A' are very close to the experimental value for the B state. The CASPT2//CASSCF potential energy curves (PECs) were calculated for F-loss dissociation from the five C(s) states and H-loss dissociation from the 1(2)A', 1(2)A', and 2(2)A' states. The electronic states of the CH3+ and CH2F+ ions as the dissociation products were carefully determined by checking the energies and geometries of the asymptote products, and appearance potentials for the two ions in different states are predicted. The F-loss PEC calculations for CH3F+ indicate that F-loss dissociation occurs from the 1(2)A', 1(2)A', and 2(2)A' states [all correlating with CH3+(X1A1')], which supports the experimental observations of direct dissociation from the X and A states, and that direct F-loss dissociation can occur from the two Jahn-Teller component states of B2E, 2(2)A' and 3(2)A' [correlating with CH3+(1(3)A') and CH3+(1(3)A'), respectively]. Some aspects of the 3(2)A' Cl-loss PEC of the CH3Cl+ ion are inferred on the basis of the calculation results for CH3F+. The H-loss PEC calculations for CH3F+ indicate that H-loss dissociation occurs from the 1(2)A', 1(2)A', and 2(2)A' states [correlating with CH2F+(1(3)A'), CH2F+(X1A1), and CH2F+(1(1)A'), respectively], which supports the observations of direct dissociation from the X and B states. As the 2(2)A' H-loss PEC of CH3Cl+, the 2(2)A' H-loss PEC of CH3F+ does not lead to H + CH2X+, but the PECs of the two ions represent different types of reactions.