Excitation Energies of Canonical Nucleobases Computed by Multiconfigurational Perturbation Theories

In this computational work, we assessed the performance of ab initio multireference (MR) methods for the calculation of vertical excitation energies of five nucleobases: adenine, guanine, cytosine, thymine and uracil. In total, we have studied 38 singlet and 30 triplet excited states. Where possible we used the multireference configuration interaction (MRCI) method as a reference for various flavors of multireference perturbation theory to second order. In particular, we have benchmarked CASPT2, NEVPT2 and XMCQDPT2. For CASPT2, we have analyzed the single‐state, multistate (MS) and extended MS variants. In addition, we have assessed the effect of the ionization potential electron affinity (IPEA) shift. For NEVPT2, we have used the partially and the strongly contracted variants. Further, we have tested the commonly used RI‐CC2, RI‐ADC2 and EOM‐CCSD methods. Generally, we observe the following trends for singlet excited states: NEVPT2 is the closest MR method to MRCISD+Q, closely followed by CASPT2 with the default IPEA shift. The same trend is observed for triplet states, although NEVPT2 and CASPT2‐IPEA are getting closer. Interestingly, the n, π* singlet excited states were described more accurately than π, π* excited states, while for triplet states the trend is inverted except for NEVPT2. This work is an important benchmark for future photochemical investigations.     

Triplet states of the amide group. Trapped electron spectra of formamide and related molecules

Trapped electron (TE) spectra are obtained using ion cyclotron resonance detection of scavenged electrons. The lowest singlet-triplet transitions, ³(n→π^*), in formamide (HCONH₂) and N,N-dimethyl formamide (HCONMe₂) are found at vertical energies of 5.30 and 5.00 eV, respectively. An unresolved band containing the ³(π→π^*) and ³(n→3s) states appears at higher energies, centered at 6.60 and 6.00 eV, respectively. The TE spectra of formaldehyde (HCHO), acetaldehyde (MeCHO) and acetone (Me₂CO) are obtained for comparison and are used along with results from ab initio theoretical calculations in establishing assignments. Singlet-triplet transitions dominate the spectra of all of these carbonyl containing molecules, to the exclusion of low lying singlet-singlet transitions. This is in agreement with other TE spectra and the expectation that (dσ/dE) will be higher near threshold for singlet-triplet as compared to singlet-singlet transitions.     

Semiempirical molecular orbital studies of phthalocyanines

Theπ andσ lone pair electron system of the phthalocyanine molecule has been studied by a semiempirical SCF-MO method. Electronic transitions of bothπ-π ^* andn-π ^* types have been considered. The excited states have been calculated by means of the method of superposition of configurations where all singly excited states are included. Assignments for the electronic spectrum of phthalocyanine could be made in good agreement with experiment. The position of the lowest electronically allowedn-π ^* transition is predicted to be found in the region of the strong Soret band.     

A MS-CASPT2 study of the low-lying electronic excited states of CH2BrCl

The low-lying singlet excited states of CH₂BrCl have been calculated using multiconfigurational CASSCF, second-order perturbation theory CASPT2 and its multistate extension MS-CASPT2. The CASSCF method shows spurious valence–Rydberg mixing and a wrong order of states. Inclusion of dynamical correlation by single root CASPT2 lowers dramatically the energy of the valences states but does not lead to a complete separation between valence and Rydberg states. This situation is improved by the MS-CASPT2 calculations, which gives two valence states for both A^′ and A″ symmetries below the lowest Rydberg state, corresponding to n(Br)→σ^*(C–Br) and n(Cl)→σ^*(C–Cl) transitions at 6.1 eV (203 nm) and 7.2 eV (173 nm), and being repulsive along C–Br and C–Cl coordinates.     

Quasirelativistic valence ab initio calculation of the potential-energy curves for Cd–rare gas atom pairs

The results of large-scale valence ab initio calculations of the potential-energy curves for the ground states and several excited states of Cd–rare gas (RG) van der Waals molecules are reported. In the calculations, Cd²⁰⁺ and RG⁸⁺ cores are simulated by energy-consistent pseudopotentials, which also account for scalar-relativistic effects and spin-orbit interaction within the valence shell. The potential energies of the Cd–RG species in the ΛS coupling scheme have been evaluated by means of ab initio complete-active-space multiconfiguration self-consistent-field (CASSCF)/CAS multireference second-order perturbation theory (CASPT2) calculations with a total 28 valence electrons, but the spin-orbit matrix has been computed in a reduced configuration interaction space restricted to the CASSCF level. Finally, the Ω potential curves are obtained by diagonalization of the modified spin-orbit matrix (its diagonal elements before diagonalization substituted by the corresponding CASPT2 eigenenergies). The calculated potential curves, especially the spectroscopic parameters derived for the ground states and several excited states of the Cd–RG species are presented and discussed in the context of available experimental data. The theoretical results exhibit very good agreement with experiment. Received: 20 April 2000 / Accepted: 1 September 2000 / Published online: 21 December 2000     

Non-empirical SCF and Cl Study of the ground and excited states of thioformaldehyde

Ab initio SCF and Cl calculations are reported for ground and various low-lying Rydberg and valence excited states of thioformaldehyde H₂CS. A double-zeta basis of near Hartree-Fock quality is employed in this work and the importance of polarization functions is also assessed. The calculations indicate uniformly larger CX bond lengths in this system than for H₂CO in the corresponding electronic states; they also lind potential minima for H₂CS non-planar nuclear conformations in the (n,π^*) and (π,π^*) excited states but in each case the calculated inversion barriers are seen to be smaller than those encountered in formaldehyde. The vertical transition energies to the various excited states studied are also found to be significantly smaller in H₂CS than in H₂CO but the order of electronic states is concluded to be virtually identical for the two systems. The lowest-lying excited states are the ^3,1(n,π^*) species calculated at 1.84 and 2.17 eV respectively; the first two allowed transitions are indicated to be the Rydberg species (n,s_R) and (n,px_R) at 5.83 and 6.62 eV. These are followed by the two allowed transitions σ → π^* and π → π^* at 7.51 and 7.92 eV respectively, both well below the first ionization limit in H₂CS. The much smaller splitting between the ^3,1(π,π^*) species in H₂CS than in H₂CO is attributed to the relatively diffuse charge distribution of the sulfur atom compared to that of oxygen.     

On the lower excited electronic states of biacetyl

An ab initio SCF and CI study has been carried out for the ground and electronically excited states of biacetyl (CH₃COCOCH₃). The second absorption band in the 4.40 eV region has been assigned to a ¹A_g → ¹B_g nπ^* transition The character of the lower-lying states has been analyzed in terms of the CI wavefunctions.     

Activite optique de la methyl-3 cyclopentenanone; un calcul ab initio

Ab initio calculations for two conformers of 3-methyl cyclopentanone, with and without 3d orbitals, show that inclusion of d orbitals on the carbonyl group lower energies and is thus necessary, as follows from the variation principle; it also markedly affects rotatory strengths of the n→π^* and π→π^* transitions, but conformational equilibria of the cyclic skeleton and the methyl group are almost independent of the basis chosen. Internal rotation of the methyl group exerts a negligible influence on the rotatory strength of the n→π^* transition, so that this group can be considered as a sphere in this respect.     

Ab initio excited states calculations of Kr 3 + , probing semi-empirical modelling

The accuracy of the diatomics-in-molecules (DIM) model for the krypton ionic trimer is examined in a series of ab initio calculations. In the C_2v symmetry, the ground states of irreducible representations B₂ and A₁ were calculated using partially spin restricted open-shell coupled cluster method with perturbative triple connections (RHF-RCCSD-T), the relativistic effective core potential (RECP) and an extended basis set of atomic orbitals. Internally contracted multireference configuration interaction method (icMRCI) with the extended and restricted basis set was used to generate the potential energy surfaces (PESs) of the nine electronic states of Kr ₃ ⁺ corresponding to Kr(¹S) + Kr(¹S) + Kr⁺(²P) dissociation limit in a wide interval of nuclear geometries. The overall agreement of the accurate ab initio PESs and the diatomics-in-molecules PESs confirms the quality of the DIM Hamiltonian for the Kr ₃ ⁺ clusters and justifies its use in dynamical and spectroscopic studies of the Kr _n ⁺ clusters. Inclusion of the spin–orbit coupling into the ab initio PESs through a semi-empirical scheme is proposed.     

Electronic structure of pyrazine: a study of the geometries in the lowest excited singlet π *← n and π*←π states

Electronic, excitation energies, charge distributions and geometries of pyrazine in the lowest excited singlet π^*←π and π^*←n states have been studied by the VE—PPP, CNDO/2 and CNDO/s-CI molecular orbital methods. Study of the change of geometry in the π^*←n excitation requires localization of the density matrices in the ground and excited states, and with the help of these σ-bond orders are defined. Charge distributions and bond orders in the lowest excited singlet π^*←π and π^*←n states are compared. Whereas in the lowest singlet π*←π excitation the pyrazine ring expands uniformly, in the case of the π^*←n excitation C-C bonds contract whereas C-N bonds elongate. The predictions of theory are in agreement with experimental results, showing that the method used can be employed to obtain reliably the trends of geometry changes following a π^*←π excitation of a molecule before a more complete theoretical or experimental study is performed.     

On the lower excited states of the formyl radical HCO

Extensive ab initio Cl calculations were carried out on the lower four doublet states of HCO. The results are used in assigning hydrocarbon flame bands, the upper state of the A and B bands being the B²A′ and C²A″ states. These two states are characterized as bonding of H with the π-to-π^* state of CO.     

Pseudopotential approach of the electronic structure in clusters: application to Alkali Halides and Rare Gases

We examine here the use of pseudopotentials in limit cases where only a very small number of electrons (much less than the usual number of valence electrons), eventually excited, can be singled out and considered as active to determine the electronic structure. Two applications are considered. The first one concerns the family of nonstoechiometric ionic clusters for which only the excess electrons are active. The relevance and accuracy of such ab initio determined pseudopotentials is illustrated on Na_nF_n-1 clusters which are treated as one-electron systems. The second application concerns excited states of neutral rare gas clusters Rg _n ^* . In this case, the excited electron may be coupled to different core states and a resonant hole-particle treatment involving e-Rg and e-Rg ⁺ pseudopotentials is presented.     

Theoretical studies on the electronic structures and spectroscopic properties for a series of Osmium(II)-2,2′,6′,2′′-terpyridine complexes

A series of Osmium(II) complexes [Os (trpy-R)₂]²⁺(trpy=2,2′,6′,2′′-terpyridine and R=H (1), OH (2), and C₆H₅(3)) have been investigated by the density functional (DF) and ab initio calculations. The structures of 1–3 in the ground and excited states were fully optimized at the B3LYP and CIS level, respectively, and their absorption and emission spectra in the acetonitrile solution were obtained using the TD-DFT (B3LYP) method associated with the PCM model. The calculations indicated that, for 1–3, the variation of the substituents on the terpyridine ligand only slightly changes their geometrical structures in the ground and excited states but leads to a sizable difference in the electronic structures. The results show that the low-lying MLCT/ILCT transitions (at 446 (1), 465 (2), and 499 nm (3)) are red-shifted according to the electron-donating ability of substituents on the terpyridine ligand, but blue-shift trend of the high-lying ILCT transitions (at 301 (1), 297 (2), and 272 nm (3)). It also reveals that the lowest energy emissions of 1–3 at 649 nm, 656 nm, and 676 nm have the character of mixing ³[π*(trpy) → d(Os)] and ³ ππ* (³MLCT/³ILCT) transitions localized on the terpyridine ligand, which are identical to the transition properties of the lowest-energy absorptions.     

Spektroskopische Untersuchung einiger α,β-ungesättigter cyclischer Ketone

The results of the spectroscopic investigation of the steroidal enones 1–6 can be summarized as follows:

• 1. Direct absorption and phosphorescence excitation techniques have been used to locate the ³(n,π*) states, and in each case it has been found to be the second triplet state.
• 2. The lowest excited state in each case is assigned as ³(π,π*) state.
• 3. The diffuseness in the phosphorescence emission from the ³(π,π*) states is attributed to a large change in the molecular geometry upon excitation (probably to a non-planar configuration).
• 4. The diffuseness in the S→ T_n,π* absorption is correspondingly attributed to interaction between the ³(n,π*) and ³(π,π*) states. A summary of the energy levels for these compounds is given in Fig. 4.

     

Potential energy curves of several excited states of the Ne2* excimer: Assignment of the transient absorption spectra of the excimer

Most of the excited states of Ne₂, which are correlated with the Rydberg state transitions 2p → 3s, 3p, and 4s of Ne, are studied by ab initio CI calculations. Two transient absorption spectra from the lowest excimer state Σ_u⁺ recently observed by Arai et al., are discussed on the basis of calculated potential energy curves. Possible assignments are presented. The calculated transition energies are in good agreement with the observed ones.     

Computational insight into excited states of the ring-opening radicals from the pyrolysis of furan biofuels

The low-lying valence excited states and Rydberg states of the radical species from the ring-opening reactions in pyrolysis of furan biofuels have been determined by extensive density functional theory and sophisticated wave function theory calculations. The radicals 1-C₄H₅O-2, 2-furylCH₂, and 4-C₆H₇O with the delocalized π-type single electron are predicted to be most stable among the reactive species here for furan, 2-methyfuran, and 2,5-dimethylfuran, respectively. Predicted vertical transition energies by TD-CAM-B3LYP show good agreement with those by CASPT2. Some among the electronic excitations to low-lying states can take place in the visible light region, and they may be involved in the combustion process. Further surface hopping dynamics simulations on the excited states of the most stable ring-opening radical 1-C₄H₅O-2 of furan as an example reveal that 89.9% sampling trajectories at the initial excited state of 2²A”(π¹π*²) decay to the 1²A’(n¹π*²) state within an average of 384 fs, and then 81.2% trajectories at the 1²A’ state go to the ground state within an average of 114 fs. At the end of the simulation for 1000 fs, 18.8% trajectories still stay on the excited states of 2²A” and 1²A’, suggesting that the reactive radicals in the ground state are mainly responsible for the combustion chemistry of furan biofuels. © 2018 Wiley Periodicals, Inc.     

Effect of nonproximate atomic substitution on excited state intramolecular proton transfer

The central C atom of the OCCCO skeleton of the malonaldehyde molecule is replaced by N, and the effects upon the intramolecular H-bond and the proton transfer are monitored by ab initio calculations in the ground and excited electronic states. The H-bond is weakened in the singlet and triplet states arising from n→π* excitation in both molecules, which is accompanied by a heightened barrier to proton transfer.³ππ* behaves in the same manner, but the singlet ππ* state has a stronger H-bond and lower barrier. Replacement of the central C atom by N strengthens the intramolecular H-bond. Although the proton transfer barrier in the ground state of formimidol is lower than in malonaldehyde, the barriers in all four excited states are higher in the N-analog. The latter substitution also dampens the effect of the n→π* excitation upon the H-bond and increases the excitation energies of the various states, particularly ππ*. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 129–138, 1998     

Theoretical studies on structures and electronic spectra of linear carbon chains C2nH+ (n = 1−5)

The density functional theory (DFT) and the complete active space self‐consistent‐field (CASSCF) method have been used for full geometry optimization of carbon chains C_2nH⁺ (n = 1–5) in their ground states and selected excited states, respectively. Calculations show that C_2nH⁺ (n = 1–5) have stable linear structures with the ground state of X³Π for C₂H⁺ or X³Σ^? for other species. The excited‐state properties of C_2nH⁺ have been investigated by the multiconfigurational second‐order perturbation theory (CASPT2), and predicted vertical excitation energies show good agreement with the available experimental values. On the basis of our calculations, the unsolved observed bands in previous experiments have been interpreted. CASSCF/CASPT2 calculations also have been used to explore the vertical emission energy of selected low‐lying states in C_2nH⁺ (n = 1–5). Present results indicate that the predicted vertical excitation and emission energies of C_2nH⁺ have similar size dependences, and they gradually decrease as the chain size increases. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009