Ab initio investigation on the structures and spectra of the firefly luciferin

The ground state (S0) geometry of the firefly luciferin (LH2) was optimized by both DFT B3LYP and CASSCF methods. The vertical excitation energies (Tv) of three low-lying states (S1, S2, and S3) were calculated by TD-DFT B3LYP//CASSCF method. The S1 geometry was optimized by CASSCF method. Its Tv and the transition energy (Te) were calculated by MS-CASPT2//CASSCF method. Both the TD-DFT and MS-CASPT2 calculated S1 state Tv values agree with the experimental one. The IPEA shift greatly affects the MS-CASPT2 calculated Tv values. Some important excited states of LH2 and oxyluciferin (oxyLH2) are charge-transfer states and have more than one dominant configuration, so for deeply researching the firefly bioluminescence, the multireference calculations are desired.     

Theoretical investigation of anthracene-9,10-endoperoxide vertical singlet and triplet excitation spectra

The electronic absorption spectrum of anthracene-9,10-endoperoxide (APO) has been investigated by means of multiconfigurational multi-state second order perturbation theory on complete active space self-consistent field wavefunctions (MS-CASPT2/CASSCF) and two single reference methods: time-dependent density functional theory (TD-DFT) and coupled cluster of second order (CC2). After testing several active spaces and basis sets, a CAS (14,12) active space together with an ANO-S basis set was found an appropriate choice to describe the vertical singlet and triplet electronic states of APO. Unfortunately, TD-DFT and CC2 methods cannot reproduce the MS-CASPT2 and experimental spectrum. Our MS-CASPT2//CASSCF(14,12)/ANO-S calculations predict a predominant pi*(OO)sigma*(OO) character for the lowest singlet excited state S(1) at 3.85 eV. Accordingly, the lowest singlet state of APO should be responsible for homolysis of the endoperoxide group. The next two absorbing excited states, experimentally proposed to be responsible for singlet oxygen production and therefore connected to the biological interest of APO, have been computed vertically at 4.34 and 4.59 eV and assigned to pi(CC)pi*(CC) and pi*(OO)pi*(CC) transitions, respectively. The vertical triplet electronic spectrum follows the singlet vertical spectrum ordering. The high density of triplet and singlet excited states of different nature within few eV points to the possibility of intersystem crossings between potential energy surfaces of different multiplicity.     

Ab initio investigation on the structures and spectra of the firefly luciferin

The ground state (S₀) geometry of the firefly luciferin (LH2) was optimized by both DFT B3LYP and CASSCF methods. The vertical excitation energies (T _v) of three low-lying states (S₁, S₂, and S₃) were calculated by TD-DFT B3LYP//CASSCF method. The S₁ geometry was optimized by CASSCF method. Its T _v and the transition energy (T _e) were calculated by MS-CASPT2//CASSCF method. Both the TD-DFT and MS-CASPT2 calculated S₁ state T _v values agree with the experimental one. The IPEA shift greatly affects the MS-CASPT2 calculated T _v values. Some important excited states of LH2 and oxyluciferin (oxyLH2) are charge-transfer states and have more than one dominant configuration, so for deeply researching the firefly bioluminescence, the multireference calculations are desired. Supported by the National Natural Science Foundation of China (Grant No. 20673012) and the Major State Basic Research Development Programs (Grant No. 2004CB719903)     

Influence of localized excited states on the transition moment directions of charge transfer complex absorptions

The influence of localized excited (LE) states on the spectroscopy of charge transfer (CT) complexes has been examined for a series of complexes formed between methyl-substituted benzene donors and 1,2,4,5-tetracyanobenzene as acceptor in 1,2-dichloroethane and octanenitrile solvents. A molecular orbital model was used to describe the appearance of multiple CT absorption bands that occur in the spectra of these complexes. The influence of LE states in these CT absorptions was explored using time-resolved linear dichroism spectroscopy where the direction of the CT transition moment vector (TMV) was used to probe the magnitude of intensity borrowing. The TMV directions for each of the observed CT transitions within the absorption spectra were determined for several complexes. In some cases, the observed CT transitions were interpreted as being pure CT transitions; in others the observed transitions are influenced significantly by a LE transition. The correlation between the TMV directions and the transition energy suggests that the magnitude of intensity borrowing is influenced not only by the energy difference between the CT and LE transitions but also by the specific character of the transitions under consideration.     

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.     

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

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

Ultraviolet absorption spectra of substituted phenols: a computational study

Vertical excitation energies for electronic transitions from the ground state to the first two excited states of phenol, mono- and disubstituted methoxyphenols and methyl-substituted phenols have been characterized with the Time-Dependent Density Functional Theory (TD-DFT), the Complete Active Space Self-Consistent Field method (CASSCF) and the Coupled Cluster with Single and Double Excitations Equation-of-Motion approach (CCSD-EOM) to simulate and interpret experimental ultraviolet absorption spectra. While CASSCF excitation energies for the first two transitions either are grossly overestimated or exhibit a weak correlation with experimental data, both TD-DFT and CCSD-EOM perform very well, reproducing the spectral shifts of both the primary band and secondary band observed upon substitution. The conformational dependence of the calculated excitation energies is generally smaller than the shifts caused by substitution.     

Electronic excitations of hypervalent configurations in amorphous selenium: Quantum-chemical modeling

The first excited states of the clusters modeling fragments of a continuous random network and the most stable local defects in solid selenium have been calculated by different ab initio methods (TD-DFT, EOM CCSD, CASSCF). A split-valence double-zeta basis set augmented with polarization and diffuse functions has been used. The effect of the cluster size and structure on the energies of transitions is considered, and the tendencies of geometric relaxation of excited states for the most typical structures are assessed.     

Understanding the differences in photochemical properties of substituted aminopyrimidines

The luminescent patterns of several members of the aminopyrimidine family are very different, showing not fluorescence at all, only a fluorescence band, normal or anomalous, or dual fluorescence, depending on the substituents and on the environment (gas phase vs. polar solvents). In this work, we study the lowest excited states of several members of this family that exhibit different fluorescence patterns to try to explain their photochemistry and to understand the effect of the substituents and the environment. We have found that several excited states (local excited (LE), charge transfer (CT) and n _N?C??* states) have minima on the lowest excited potential energy surface (S₁), being their relative energy the determinant factor of the luminescent behavior. If the more stable S₁ minima are of n _N?C??* character, a non-radiative deexcitation channel is the most efficient and the system shows no fluorescence. If the CT and/or LE states are the most stable, the non-radiative deactivation channel is not accessible and the system fluoresces. The relative energies of the CT and LE minima (affected by substituents and by the presence of a polar solvent) and the different magnitude of the oscillator strength for the radiative transition to the ground state determine which emission is more efficient, giving place to normal, anomalous or dual fluorescence. The study has been carried out by CASSCF/CASPT2 computations, including the solvent effect by means of the PCM model.     

Photochemistry of CH3Mn(CO)5: a multiconfigurational ab initio study

The electronic spectroscopy of CH3Mn(CO)5 has been investigated by means of ab initio multiconfigurational MS-CASPT2/CASSCF calculations. The absorption spectrum is characterized by a series of Metal-Centered (MC) excited states in the UV energy domain (below 290 nm) that could be responsible for the observed photoreactivity starting at 308 nm. The upper part of the spectrum is overcrowded between 264 and 206 nm and dominated by a high density of Metal-to-Ligand-Charge-Transfer (MLCT) states corresponding mainly to 3d(Mn) --> pi*(CO) excitations. A non-negligible contribution of Metal-to-sigma-Bond-Charge-Transfer (MSBCT) states corresponding to 3d(Mn) --> sigma*(Mn-CH3) excitations is also present in the theoretical spectrum of CH3Mn(CO)5. However, in contrast to other transition metal hydrides and methyl substituted (HMn(CO)5, HCo(CO)4, and CH3Co(CO)4) these MSBCT transitions do not participate to the lowest bands of the spectrum as main contributions. The photochemistry of CH3Mn(CO)5, namely the loss of a CO ligand vs. the metal-methyl bond homolysis, is investigated by means of MS-CASPT2 states correlation diagrams. This study illustrates the complexity of the photodissociation mechanism of this class of molecules, which involves a large number of nearly degenerate electronic states with several channels for fragmentation.     

MS-CASPT2 calculation of excess electron transfer in stacked DNA nucleobases

Calculations using the complete active space self-consistent field (CASSCF) and complete active space second-order perturbation (CASPT2) methods, and the multistate formulation of CASPT2 (MS-CASPT2), are performed for the ground and excited states of radical anions consisting of two pi-stacked nucleobases. The electronic couplings for excess electron transfer (EET) in the pi-stacks are estimated by using the generalized Mulliken-Hush approach. We compare results obtained within the different methods with data derived using Koopmans' theorem approximation at the Hartree-Fock level. The results suggest that although the one-electron scheme cannot be applied to calculate electron affinities of nucleobases, it provides reasonable estimates for EET energies. The electronic couplings calculated with KTA lie between the CASPT2 and the MS-CASPT2 based values in almost all cases.     

Can TD‐DFT calculations accurately describe the excited states behavior of stacked nucleobases? The cytosine dimer as a test case

By using calculations rooted in the time dependent density functional theory (TD-DFT) we have investigated how the lowest energy excited states of a face-to-face pi-stacked cytosine dimer vary with the intermonomer distance (R). The perfomances of different density functionals have been compared, focussing mainly on the lowest energy single excited state of the dimer (S(1))(2). TD-PBE0, TD-LC-omegaPBE, and TD-M05-2X provide a picture very similar to that obtained at the CASPT2 level by Merchan et al. (J Chem Phys 2006, 125, 231102), predicting that (S(1))(2) has a minimum for R approximately 3 A, with a binding energy of approximately 0.5 eV, whereas TD-B3LYP, TD-CAM-B3LYP, and TD-PBE understimate the binding energy. However, independently of the functional employed, no low-energy spurious charge transfer transitions are predicted by TD-DFT calculations, also when a nonsymmetric dimer is investigated, providing encouraging indications for the use of TD-DFT for studying the excited state of pi-stacked nucleobases.     

Theoretical study of the singlet electronically excited states of N4

Vertical excitation energies for the lowest eleven singlet states of T_d N₄ were calculated using the TD-DFT method with the B3LYP functional, and at the EOM-CCSD level of theory. The vertical excitation energies for the five lowest-lying excited states were also obtained using the state-averaged CASSCF, CASPT2, CASPT3, and MRCI + Q methods. Our results show that the five lowest-lying states are of valence character. EOM-CCSD/d-aug-cc-pVTZ calculations predict that there are two weakly allowed optical transitions of T₂ symmetry at 10.44 and 10.82 eV. The transition to the third T₂ state, which is predicted to be at 10.89 eV, has an oscillator strength about one order of magnitude higher.     

Ligand-to-diimine/metal-to-diimine charge-transfer excited states of [Re(NCS)(CO)3(alpha-diimine)] (alpha-diimine = 2,2'-bipyridine, di-iPr-N,N-1,4-diazabutadiene). A spectroscopic and computational study

Two new complexes fac-[Re(NCS)(CO)3(N,N)] (N,N = 2,2'-bipyridine (bpy), di-iPr-N,N-1,4-diazabutadiene (iPr-DAB)) were synthesized and their molecular structures determined by X-ray diffraction. UV-vis absorption, resonance Raman, emission, and picosecond time-resolved IR spectra were measured experimentally and calculated with TD-DFT. A good agreement between experimental and calculated ground- and excited-state spectra is obtained, but only if the solvent (MeCN) is included into calculations and excited state structures are fully optimized at the TD-DFT level. The lowest excited states of the bpy and iPr-DAB complexes are assigned by TD-DFT as 3aA' by comparison of calculated and experimental IR spectra. Excited-state lifetimes of 23 ns and ca. 625 ps were determined for the bpy and DAB complex, respectively, in a fluid solution at room temperature. Biexponential emission decay (1.3, 2.7 micros) observed for [Re(NCS)(CO)3(bpy)] in a 77 K glass indicates the presence of two unequilibrated emissive states. Low-lying electronic transitions and excited states of both complexes have a mixed NCS --> N,N ligand-to-ligand and Re --> N,N metal-to-ligand charge-transfer character (LLCT/MLCT). It originates in mixing between Re d(pi) and NCS pi characters in high-lying occupied MOs. Experimentally, the LLCT/MLCT mixing in the lowest excited state is manifested by shifting the nu(CO) and nu(NC) IR bands to higher and lower wavenumbers, respectively, upon excitation. Resonant enhancement of both nu(CO) and nu(NC) Raman bands indicates that the same LLCT/MLCT character mixing occurs in the lowest allowed electronic transition.     

Structure and bonding of Ag(I)-DNA base complexes and Ag(I)-adenine-cytosine mispairs: an ab initio study

High level ab initio calculations have been carried out to characterize the structure, bonding and energetics of Ag(I)-DNA base complexes, including adenine or cytosine, as well as Ag(I)-adenine-cytosine mispairs. The interactions of the Ag cation in all binding sites of all adenine and cytosine tautomers have been considered. The calculations show that in gas phase the canonical form of cytosine is stabilized upon metalation, whereas the lowest energy structure of Ag-adenine correspond to a rare tautomer. Interestingly, the theoretical inspection of metalated adenine-cytosine mispair reveals that the most stable structures are formed from the canonical cytosine and adenine tautomers. The lowest energy structure is planar with adenine and cytosine hydrogen-bonded. Within few kcal/mol nonplanar, conformationally very flexible structures are found, in which the Ag(I) crosslinks an endocyclic nitrogen of adenine and the oxygen of cytosine. Metalated reverse-Wobble type of structures, on the contrary, are predicted much higher in energy.     

TD-DFT Accuracy in Determining Excited-state Structures and Fluorescence Spectra of Firefly Emitter

We analyzed the excited-state structures and emission spectra of firefly emitter, the anionic keto form of firefly oxyluciferin(keto-1), determined by the time dependent-density functional theory(TD-DFT) approach. The analysis is based on a direct comparison with the highly correlated CASSCF(MS-CASPT2) ab initio approach. 49 DFT functionals were considered and applied to the study. Among the tested functionals, mPW3PBE, B3PW91 and B3P86 give the best performance for ground-state geometry, absorption spectrum, excited-state geometry and emission spectrum.     

Theoretical Understanding on the Facilitated Photoisomerization of a Carbonyl Supported Borane System

Boron compound BOMes₂ containing an internal B−O bond undergoes highly efficient photoisomerization, followed by sequential structural transformations, resulting in a rare eight-membered B, O-heterocycle (S. Wang, et al. Org. Lett. 2019 , 21, 5285–5289). In this work, the detailed reaction mechanisms of such a unique carbonyl-supported tetracoordinate boron system in the first excited singlet (S₁) state and the ground (S₀) state were investigated by using the complete active space self-consistent field and its second-order perturbation (MS-CASPT2//CASSCF) method combined with time-dependent density functional theory (TD-DFT). Moreover, an imine-substituted tetracoordinated organic boron system (BNMes₂) was selected for comparative study to explore the intrinsic reasons for the difference in reactivity between the two types of compounds. Steric factor was found to influence the photoisomerization activity of BNMes₂ and BOMes₂. These results rationalize the experimental observations and can provide helpful insights into understanding the excited-state dynamics of heteroatom-doped tetracoordinate organoboron compounds, which facilitates the rational design of boron-based materials with superior photoresponsive performances.     

Conformational effects in UV absorption spectra of tetrasilanes

UV absorption of cyclic carbosilanes (SiMe₂)₄(CH₂)_n, N = 1–4 (1–4), and Si₄Me₁₀ (5) provides an experimental counterpart to the singlet transition energy and intensity correlation diagrams for the syn-anti conformational transformation in tetrasilane. A new third transition is found between the two previously known singlet transitions. Transition energies are nearly independent of the dihedral angle, while intensities vary widely. All trends agree with CIS/3–21G^*//HF/3–21G^* calculations. The Sandorfy C and ladder C models of σ conjugation fail to describe electronically excited states of tetrasilane, since they do not consider the lateral bonds to substituents.