Theoretical prediction on HAlS+ and HSAl+ cations using multiconfiguration second‐order perturbation theory

Some low‐lying states of the HAlS⁺ and HSAl⁺ cations have been studied for the first time by large‐scale theoretical calculations using three methods: complete active space self‐consistent field (CASSCF), complete active second‐order perturbation theory (CASPT2), and density functional theory Becke's three‐parameter hybrid function with the nonlocal correlation of Lee–Yang–Parr (B3LYP) with the contracted atomic natural orbital (ANO‐L) and cc‐pVTZ basis sets. The geometries of all stationary points along the potential energy surfaces (PESs) were optimized at the CASSCF/ANO‐L and B3LYP/cc‐pVTZ levels. The ground and the first excited states of linear HAlS⁺ are predicted to be X²Π and A²Σ⁺ states, respectively. For the linear HSAl⁺ structure, the first excited state is A²Σ⁺. The X²Π state of linear HSAl⁺ is a second‐order saddle point, because it has two imaginary frequencies. Two bent global minima M1 and M2 were found along the 1²A′ and 1²A″ PESs, respectively. The CASPT2/ANO‐L potential energy curves of isomerization reactions were calculated as a function of HAlS bond angle. According to our calculations, the ground‐state HAlS⁺ is linear, whereas the ground‐state HSAl⁺ is bent. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Theoretical Study on the Formation of H‐ and O‐Atoms,HONO, OH,NO, and NO2 from the Lowest Lying Singlet and Triplet States in Ortho‐Nitrophenol Photolysis

The photolysis of nitrophenols was proposed as a source of reactive radicals and NOx compounds in polluted air. The S₀ singlet ground state and T₁ first excited triplet state of nitrophenol were investigated to assess the energy dependence of the photofragmentation product distribution as a function of the reaction conditions, based on quantum chemical calculations at the G3SX//M06–2X/aug‐cc‐pVTZ level of theory combined with RRKM master equation calculations. On both potential energy surfaces, we find rapid isomerization with the aci‐nitrophenol isomer, as well as pathways forming NO, NO₂, OH, HONO, and H‐, and O‐atoms, extending earlier studies on the T₁ state and in agreement with available work on other nitroaromatics. We find that accessing the lowest photofragmentation channel from the S₀ ground state requires only 268 kJ/mol of activation energy, but at a pressure of 1 atm collisional energy loss dominates such that significant fragmentation only occurs at internal energies exceeding 550 kJ/mol, making this surface unimportant for atmospheric photolysis. Intersystem crossing to the T₁ triplet state leads more readily to fragmentation, with dissociation occurring at energies of ～450 kJ/mol above the singlet ground state even at 1 atm. The main product is found to be OH + nitrosophenoxy, followed by formation of hydroxyphenoxy + NO and phenyloxyl + HONO. The predictions are compared against available experimental data.     

Complete TDA and RPA Calculations on the Electronic Transitions of Fullerene‐C60 in the CNDO/S and INDO/S Approximations*

Complete single‐excitation mixing calculations on the electronic transitions of the icosahedral C₆₀ molecule have been carried out with the Tamm–Dancoff approximation (TDA) and random‐phase approximation (RPA) schemes in the CNDO/S and INDO/S approximations. The complete space of 14,400 (1p–1h) pairs is partitioned into subspaces classified according to the irreducible representations of the Ih group. For this purpose, matrix representations of the group generators are obtained on a fixed set of basis functions and are used to construct the projection operators. Degenerate molecular orbitals in each energy level are symmetry‐adapted to these projection operators. Degenerate (1p–1h) pairs or singly excited configuration wave functions are similarly symmetrized. In addition, the Clebsch–Gordan coefficients are obtained and listed in an Appendix. The TDA and RPA equations are then solved for each irreducible representation separately. Both schemes with the projection operators and with the Clebsch–Gordan coefficients gave the same results as expected, indicating that the calculations were correctly done. The transition energies from the ground state 1¹A_g to low‐lying singlet and triplet excited states and the oscillator strengths for the allowed transitions (n¹T_1u–1¹A_g) are given in tables. A proper way to normalize is discussed for the eigenvectors of the RPA‐type matrix equation. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Low‐lying singlet excited states of isocyanogen

Recent photofragment fluorescence excitation (PHOFEX) spectroscopy experiments have observed the ùA″ singlet excited state of isocyanogen (CNCN) for the first time. The observed spectrum is not completely assigned and significant questions remain about the excited states of this system. To provide insight into the energetically accessible excited states of CNCN, optimized geometries, harmonic vibrational frequencies, and excitation energies for the first three singlet excited states are determined using equation‐of‐motion coupled‐cluster theory with singles and doubles (EOM‐CCSD) and correlation‐consistent basis sets. Additionally, excited state coupled‐cluster methods which approximate the contributions from triples (CC3) are utilized to estimate the effect of higher‐order correlation on the energy of each excited state. For the ùA″ state, our best estimate for T₀ is about 42,200 cm^?1, in agreement with the experimentally estimated upper limit for the zero‐point level of 42,523 cm^?1. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

CASPT2 study on low‐lying states of HMgO and HOMg

The HMgO and magnesium monohydroxide (HOMg) have been reinvestigated using the complete active space self‐consistent field (CASSCF) and multiconfiguration second‐order perturbation theory (CASPT2) methods with the contracted atomic natural orbital (ANO) basis sets. The geometries of all stationary points along the potential energy surfaces (PESs) were optimized at the CASSCF/ANO levels. The ground and the first excited states of HMgO are predicted to be X²Π and A²Σ⁺ states, respectively. It was predicted that the ground state of HOMg is X²Σ⁺ state. The A²Π state of HOMg has unique imaginary frequency. A bent local minimum M1 was found for the first time along the 1²A″ PES and the A²Π state of HOMg should be the transition state of the isomerization reactions for M1 ? M1. The CASPT2/ANO potential energy curves of isomerization reactions were calculated as a function of HMgO bond angle. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

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.     

Theoretical study on the vertical electronic spectrum of carbonyl fluoride,F2CO

Ab initio self-consistent-field (SCF ) and configuration interaction (CI ) calculations on the ground and excited states of carbonyl fluoride (F₂CO) were carried out at its experimental ground-state equilibrium geometry. Vertical transition energies deduced from the CI results provide assignments for the electronic systems I–IV, experimentally observed by Workman and Duncan. The singlet excited state, ¹A₁ (π→π*), is found to be a mixed valence–Rydberg state and to he 1 to 1.2 eV above the suggested experimental value, irrespective of the choice of the basis used for the CI calculations.     

The Role of the nπ* 1Au State in the Photoabsorption and Relaxation of Pyrazine

The geometric, energetic, and spectroscopic properties of the ground state and the lowest four singlet excited states of pyrazine have been studied by using DFT/TD‐DFT, CASSCF, CASPT2, and related quantum chemical calculations. The second singlet nπ* state, ¹A_u, which is conventionally regarded dark due to the dipole‐forbidden ¹A_u←¹A_g transition, has been investigated in detail. Our new simulation has shown that the state could be visible in the absorption spectrum by intensity borrowing from neighboring nπ* ¹B_3u and ππ* ¹B_2u states through vibronic coupling. The scans on potential‐energy surfaces further indicated that the ¹A_u state intersects with the ¹B_2u states near the equilibrium of the latter, thus implying its participation in the ultrafast relaxation process.     

Ab initio calculations on urea

Multiconfiguration wave functions constructed from contracted Gaussian-lobe functions have been found for the ground and valence-excited states of urea. ICSCF molecular orbitals of the excited states were used as the parent configurations for the CI calculations except for the ¹A₁(π → π*) state. The ¹A₁(π → π*) state used as its parent configuration an orthogonal linear combination of natural orbitals obtained from the second root of a three-configuration SCF calculation. The lowest excited states are predicted to be the n π → π* and π → π* triplet states. The lowest singlet state is predicted to be the n π → π* state with an energy in good agreement with the one known UV band at 7.2 eV. The π → π* singlet state is predicted to be about 1.9 eV higher, contrary to several previous assignments which assumed the lowest band was a π → π* amide resonance band. The predicted ionization energy of 9.0 eV makes this and higher states autoionizing.     

An AB initio SCF MO CI study of the n → π* transition of methylenimine

The electron correlation energies of both the ground and n → π* excited states of methylenimine (CH₂NH) are investigated by means of ab initio SCF MO CI calculations. Then n → π* singlet and triplet state energies of methylenimine are obtained through 3461-dimensional CI including the singly, doubly and triply excited configurations. the excitation energy from the ground state to the ¹(n → π*) state nearly coincides with that obtained in the framework of the singly excited configuration interaction (SECI) procedure. This result suggests that there is good cancellation of the correlation energy between the ground and the excited singlet sates, proving the usefulness of the SECI method for the excitation energies.     

High‐level ab initio study of the N+(3P)+SH2 reactions in the gas phase: Role of spin‐forbidden pathways

The singlet and triplet potential energy surfaces involved in N⁺+SH₂ reactions have been explored using high‐level ab initio techniques. The geometries of the stationary points were optimized at the QCISD/6‐311G(df,p) level. The final energies were obtained in CCSD(T)/6‐311+G(3df,2p) single‐point calculations. The results obtained show that, although the N⁺(¹D)+SH₂ entrance channel is higher in energy than the N⁺(³P)+SH₂ one, most of the [H₂, S, N]⁺ singlet state cations are lower in energy than the corresponding triplets, due to their different bonding characteristics. Both singlet and triplet potential energy surfaces are quite close each other, and crossover between them can occur. The minimum energy crossing points were located by means of CASSCF(6,5) calculations. The spin‐orbit couplings show that the transition probability from the triplet to the singlet potential energy surface is significantly large. One of the most important consequences is that some of the products of the reaction, such as SH⁺, can be formed in typical spin‐forbidden processes. Since all the relevant structures along these pathways are much lower in energy than the reactants, this mechanism should be accessible even at low impact energies and therefore could be important in processes taking place in interstellar media. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Theoretical studies of cyclopropylsilylenes: The structures and stability of cyclopropylsilylene C3H5SiH

Ab initio molecular orbital calculations at the G2(MP2) level have been carried out on cyclopropylsilylene C₃H₅SiH. Four equilibrium structures were located. Like H₂Si, the ground state of C₃H₅SiH is singlet and the triplet is the low‐lying excited state. The singlet–triplet separation energy is 127.9 kJ/mol. The cis‐trans isomerization path of singlet cyclopropylsilylene was investigated by intrinsic reaction coordinate (IRC) calculations. The calculations show that no gauche conformers exist along the potential energy curve of the cis‐trans isomerization and the isomerization happens with a barrier of 30.1 kJ/mol. Changes (ΔH and ΔG) in thermodynamic functions, equilibrium constant K(T), and A factor and reaction rate constant k(T) in Eyring transition state theory of the cis‐trans isomerization were also calculated. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

CASPT2 study on low‐lying states of HBN and HNB radicals

In this study, some low‐lying states of the HBN and HNB radicals have been studied using multiconfiguration second‐order perturbation theory. The geometries of all stationary points along the potential energy surfaces (PESs) were optimized at the CASPT2/cc‐pVQZ level. The ground and the first excited states of HBN were predicted to be X²Π and A²Σ⁺ states, respectively. It was predicted that the ground state of HNB is X²Σ⁺ state. The A²Π state of HNB has unique imaginary frequency, which was different from the previously published results. A bending local minimum M1 was found for the first time along the 1²A″ PES, and the A²Π state of HNB should be the transition state of the isomerization reactions for M1 ? M1. The CASPT2/ANO potential energy curves (PECs) of isomerization reactions for HBN ? HNB were calculated as a function of HBN bond angle. By comparing the CASPT2 and CASSCF calculated results, we concluded that the influence of the dynamic electron correlation on HBN ? HNB system is not large. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Nonadiabatic coupling and diabatic electronic population dynamics on 11A2 and 11B1 states of ozone molecule

In this work, we examine nonadiabatic population dynamics for 1¹B₁ and 1¹A₂ states of ozone molecule (O₃). In O₃, two lowest singlet excited states, ¹A₂ and ¹B₁, can be coupled. Thus, population transfer between them occurs through the seam involving these two states. At any point of the seam (conical intersection), the Born-Oppenheimer approximation breaks down, and it is necessary to investigate nonadiabatic dynamics. We consider a linear vibronic coupling Hamiltonian model and evaluate vibronic coupling constant, diabatic frequencies for three modes of O₃, bilinear and quadratic coupling constants for diabatic potentials, displacements, and Huang-Rhys coupling constants using ab initio calculations. The electronic structure calculations have been performed at the multireference configuration interaction and complete active space with second-order perturbation theory with a full-valence complete active space self-consistent field methods and augmented Dunning's standard correlation-consistent-polarized quadruple zeta basis set to determine ab initio potential energy surfaces for the ground state and first two excited states of O₃, respectively. We have chosen active space comprising 18 electrons distributed over 12 active orbitals. Our calculations predict the linear vibronic coupling constant 0.123 eV. We have obtained the population on the 1¹B₁ and 1¹A₂ excited electronic states for the first 500 fs after photoexcitation.     

The Low‐Lying Excited States of 2,2′‐Bithiophene: A Theoretical Analysis

The low‐energy regions of the singlet→singlet, singlet→triplet, and triplet→triplet electronic spectra of 2,2′‐bithiophene are studied using multiconfigurational second‐order perturbation theory (CASPT2) and extended atomic natural orbitals (ANO) basis sets. The computed vertical, adiabatic, and emission transition energies are in agreement with the available experimental data. The two lowest singlet excited states, 1¹B_u and 2¹B_u, are computed to be degenerate, a novel feature of the system to be borne in mind during the rationalization of its photophysics. As regards the observed high triplet quantum yield of the molecule, it is concluded that the triplet states 2³A_g and 2³B_u, separated about 0.4 eV from the two lowest singlet excited states, can be populated by intersystem crossing from nonplanar singlet states.     

Theoretical investigation of the cyclic GaO2 and GaS2 molecules at DFT and correlated wave function levels

The geometries and the bonding properties have been predicted for cyclic GaO₂ and GaS₂ species at density functional theory (DFT), MPn (n=2,3,4 with different substitutions), QCISD(T), and CCSD(T) all‐electron correlation levels with 6‐311+G* basis set. The geometrical optimizations and the harmonic vibrational frequency analysis are performed using DFT and second‐order Møller–Plesset (MP2) methods. The relevant energy quantities are also calibrated at the high‐order electron correlation levels [MP3, MP4, quadratic configuration interaction (QCI), and coupled cluster (CC)]. Each species possesses a ²A₂ ground state with a higher energy level ²A₁ state. The corresponding state–state separations are about 32 kcal/mol for GaO₂ species and about 20 kcal/mol for GaS₂ species at the QCISD(T)/6‐311+G* level. The QCISD(T) and CCSD(T) calculations yield dissociation energies of 42.0 and 59.0 kcal/mol for two species, respectively, and other methods yield dissociation energies within ∼5 kcal/mol. Result analysis has indicated that the cyclic GaO₂ should be classified as superoxide and the GaS₂ species should be classified as supersulfide in their ground state, and those in the excited state (²A₁) should not be. However, the cyclic GaS₂ (²A₂) is less ionic than the GaO₂ (²A₂) and they are far less ionic than NaO₂. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 222–231, 2001     

Construction and applications of symmetrized valence bond wave functions

A method constructing symmetry-adapted bonded Young tableau bases is proposed, based on the symmetry properties of bonded tableaus and the projection operator associated with a point group. Several examples including the ground states and π excited states of O₃⁻, O₃, O₃⁺, and C₃⁻ are shown for instruction to construct the symmetrized valence bond (VB) wave function. Excitation energies of transitions from the ground states to π excited states of O₃⁻, C₃H₅, and C₃⁻ are calculated with an optimized symmetrized valence bond wave function in the σ–π separation approximation. Good agreement between the VB and experimental excitation energies is observed. The bonding features of the ground state and the first π excited singlet and triplet states for S₃ are discussed according to bonding populations from VB calculations. Both the singlet-biradical and the dipole structures have significant contributions to the ground state X ¹A₁ of S₃, while the excited state 1 ¹B₂ is essentially composed of the dipole structures, and the 1 ³B₂ excited state is comprised from a triplet-biradical structure. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 1–7, 1998     

Energy and angular differential cross sections for K + NO2 charge transfer reactions

Doubly differential cross sections, in energy and angle, are reported for the electron transfer reaction between potassium and nitrogen dioxide in a crossed beam apparatus at relative collision energies between 2.7 and 30.8 eV. The formation of NO^?₂ in its ground ¹A₁ and excited ³B₁ state has been observed. Theoretical consideration of these processes indicates that bond bending during the collision has a stronger influence on ion-pair formation than bond stretching. At the lower collision energies most of the excess energy is converted into internal energy of NO^?₂.     

Are They Linear,Bent, or Cyclic? Quantum Chemical Investigation of the Heavier Group 14 and Group 15 Homologues of HCN and HNC

The singlet potential‐energy surface (PES) of the system involving the atoms H, X, and E (the (H, X, E) system) in which X=N–Bi and E=C–Pb has been explored at the CCSD(T)/TZVPP and BP86/TZ2P+ levels of theory. The nature of the X? E bonding has been analyzed with charge‐ and energy‐partitioning methods. The calculations show that the linear isomers of the nitrogen systems lin ‐HEN and lin ‐HNE are minima on the singlet PES. The carbon compound lin ‐HCN (HCN=hydrogen cyanide) is 14.9 kcal mol^?1 lower in energy than lin ‐HNC but the heavier group 14 homologues lin ‐HEN (E=Si–Pb) are between 64.8 and 71.5 kcal mol^?1 less stable than the lin ‐HNE isomers. The phosphorous system (H, P, E) exhibits significant differences concerning the geometry and stability of the equilibrium structures compared with the nitrogen system. The linear form lin ‐HEP of the former system is much more stable than lin ‐HPE . The molecule lin ‐HCP is the only minimum on the singlet PES. It is 78.5 kcal mol^?1 lower in energy than lin ‐HPC , which is a second‐order saddle point. The heavier homologues lin ‐HPE , in which E=Si–Pb, are also second‐order saddle points, whereas the bent ‐HPE structures are the global minima on the PES. They are between 10.3 (E=Si) and 36.5 kcal mol^?1 (E=Pb) lower in energy than lin ‐HEP . The bent ‐HPE structures possess rather acute bending angles H‐P‐E between 60.1 (E=Si) and 79.7° (E=Pb). The energy differences between the heavier group 15 isomers lin ‐HEX (X=P–Bi) and the bent structures bent ‐HXE become continuously smaller. The silicon species lin ‐HSiBi is even 3.1 kcal mol^?1 lower in energy than bent ‐HBiSi . The bending angle H‐X‐E becomes more acute when X becomes heavier. The drastic energy differences between the isomers of the system (H, X, E) are explained with three factors that determine the relative stabilities of the energy minima: 1) The different bond strength between the hydrogen bonds H? X and H? E. 2) The electronic excitation energy of the fragment HE from the X ²Π ground state to the ⁴Σ^? excited state, which is required to establish a E≡X triple bond in the molecules lin ‐HEX . 3) The strength of the intrinsic X? E interactions in the molecules. The trends of the geometries and relative energies of the linear, bent, and cyclic isomers are explained with an energy‐decomposition analysis that provides deep insight into the nature of the bonding situation.     

Multireference Calculation of the Photodissociation of Benzyl Chloride

The photodissociation mechanism of benzyl chloride （BzCl） under 248 nm has been investigated by the complete active space SCF （CASSCF） method by calculating the geometries of the ground （S0） and lower excited states, the vertical （Tv） and adiabatic （T0） excitation energies of the lower states, and the dissociation reaction pathways on the potential energy surfaces （PES） of SI, TI and T2 states. The calculated results clearly elucidated the photodissociation mechanism of BzCl, and indicated that the photodissociation on the PES of T1 state is the most favorable.