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)     

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.     

The role of Ag(I) ions in the electronic spectroscopy of adenine–cytosine mispairs: A MS-CASPT2 theoretical study

Adenine–cytosine (AC) mispairs have been theoretically studied with MS-CASPT2//CASSCF methods in the presence and absence of Ag ions. The electronically excited states of the most stable AC mispair in the reverse-Wobble (RW) conformation have been compared with those of different Ag(I)–AC complexes, including (i) metalated RW conformations, and (ii) the most stable structures in gas phase which contain the Ag ion bridging A and C. The spectra of these complexes are characterized by charge-transfer (CT) and strong locally excited (LE) states. The metal-to-metal, metal-to-ligand, and Rydberg transitions are very weak in comparison to the nucleobase transitions. Attending to the LE and CT states, and except for the shifts induced by the presence of the Ag, the electronic spectrum of metalated AC mispairs resembles the one of the RW, showing two intense LE bands around 4.5 and 5.5 eV, corresponding to transitions within the adenine and cytosine π-system, respectively. Additionally TD-DFT results obtained with the B3LYP functional are compared with MS-CASPT2//CASSCF calculations. The results clearly evidence the weakness of TD-DFT to describe long range exchange interactions leading to strongly underestimated CT states.     

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 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.     

Exploring the Mechanism of a Chiral N-Alkyl Imine-Based Light-Driven Molecular Rotary Motor at MS-CASPT2//CASSCF and MS-CASPT2//(TD) DFT Levels

The working mechanism including the photoisomerization and thermal isomerization steps of a chiral N-alkyl imine-based motor synthesized by Lehn et al. are revealed by MS-CASPT2//CASSCF and MS-CASPT2//(TD-)DFT methods. For the photoisomerization process of the imine-based motor, it involves both the bright (π,π*) state and the dark (n,π*) state. In addition, the MECI has similar geometry and energy to the minimum of the S₁ state, which shows that the process is barrierless and keeps the unidirectionality of rotation well; the result confirms the imine-based motor is a good candidate for a light-driven molecular rotary motor. For the thermal isomerization process of the imine-based motor, there are two even isomerization paths: one with the mechanism of the in-plane N inversion, the energy barriers of which are 29.6 kcal mol⁻¹ at MS3-CASPT2//CAM-B3LYP level and 29.2 kcal mol⁻¹ at MS3-CASPT2//CASSCF level; the other with the mechanism of ring inversion of the cycloheptatriene moiety, with energy barriers of 28.1 kcal mol⁻¹ at MS3-CASPT2//CAM-B3LYP level and 18.1 kcal mol⁻¹ at MS3-CASPT2//CASSCF level. According to the structural feature of the stator moiety, the imine molecule can be used as a two-step or a four-step light-driven rotary motor.     

Low-lying electronic states and Cl-loss photodissociation of the C2H3Cl+ ion studied using multiconfiguration second-order perturbation theory

We studied the 1(2)A' '(X2A' '), 1(2)A' (A2A'), 2(2)A' ' (B2A' '), and 2(2)A' (C2A') states of the C2H3Cl+ ion using the complete active space self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) methods. For the four ionic states, we calculated the equilibrium geometries, adiabatic (T0) and vertical (Tv) excitation energies, and relative energies (Tv') at the geometry of the molecule at the CASPT2 level and the Cl-loss dissociation potential energy curves (PECs) at the CASPT2//CASSCF level. The computed oscillator strength f value for the X2A' ' <-- A2A' transition is very small, which is in line with the experimental fact that the A state has a long lifetime. The CASPT2 geometry and T0 value for the A2A' state are in good agreement with experiment. The CASPT2 Tv' values for the A2A', B2A' ', and C2A' states are in good agreement with experiment. The Cl-loss PEC calculations predict that the X2A' ', A2A', and C2A' states correlate to C2H3+ (XA1) and the BA' ' state to C2H3+ (1A' ') (the B2A' ' and C2A' PECs cross at R(C-Cl) approximately 2.24 A). Our calculations indicate that at 357 nm the X2A' ' state can undergo a transition to B2A' ' followed by a predissociation of B2A' ' by the repulsive C2A' state (via the B/C crossing), leading to C2H3+ (X1A1), and therefore confirm the experimentally proposed pathway for the photodissociation of X2A' ' at 357 nm. Our CASPT2 D0 calculations support the experimental fact that the X state does not undergo dissociation in the visible spectral region and imply that a direct dissociation of the A state to C2H3+ (X1A1) is energetically feasible.     

Photoexcitation of dinucleoside radical cations: a time-dependent density functional study

The excited states of dinucleoside phosphates (dGpdG, dApdA, dApdT, TpdA, and dGpdT) in their cationic radical states were studied with time-dependent density functional theory (TD-DFT). The ground-state geometries of all the dinucleoside phosphate cation radicals considered, in their base stacked conformation, were optimized with the B3LYP/6-31G(d) method. Further, to take into account the effect of the aqueous environment surrounding the dinucleoside phosphates, the polarized continuum model (PCM) was considered and the excitation energies were computed by using the TD-B3LYP/6-31G(d) method. From this study, we find that the first transition in all the dinucleoside molecules involves hole transfer from base to base. dG*+pdG and dApdA*+ were found to have substantially lower first transition energies than others with two different DNA bases. Higher energy transitions involve base to sugar as well as base to base hole transfer. The calculated TD-B3LYP/6-31G(d) transition energies are in good agreement with previous calculations with CASSCF/CAS-PT2 level of theory. This TD-DFT work supports the experimental findings that sugar radicals formed upon photoexcitation of G*+ in gamma-irradiated DNA and suggests an explanation for the wavelength dependence found.     

OClO里德堡态激发能的准确预测及其阴离子低能激发态的从头算研究

采用全活化空间自洽场方法(CASSCF)研究了OClO阴离子7个低能电子态及其自由基的基态. 为了进一步考虑动态电子相关效应, 采用二级多组态微扰理论(CASPT2)获得更加可靠的能量值. 此外, 在ANO-L基组的基础上, 在OClO自由基的电荷中心增加了为研究里德堡态所建立的1s1p1d的波函数, 并应用多组态二级微扰理论(MS-CASPT2)方法获得了里德堡态的准确电子激发能.     

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.     

A theoretical study on 1-(thiophen-2-yl-methyl)-2-(thiophen-2-yl)-1H-benzimidazole

The molecular geometry and vibrational frequencies of 1-(thiophen-2-yl-methyl)-2-(thiophen-2-yl)-1H-benzimidazole (C(16)H(12)N(2)S(2)) in the ground state has been calculated using the Hartree-Fock (HF) and density functional method (B3LYP) with 6-31G(d) basis set. The optimized geometric bond lengths and bond angles obtained by using HF and DFT (B3LYP) show the best agreement with the experimental data. Comparison of the observed fundamental vibrational frequencies of 1-(thiophen-2-yl-methyl)-2-(thiophen-2-yl)-1H-benzimidazole (C(16)H(12)N(2)S(2)) and calculated results by density functional B3LYP and Hartree-Fock methods indicate that B3LYP is superior to the scaled Hartree-Fock approach for molecular vibrational problems.     

Theoretical studies of molecular structure and vibrational spectra of 2-amino-5-phenyl-1,3,4-thiadiazole

The molecular geometry and vibrational frequencies of 2-amino-5-phenyl-1,3,4-thiadiazole (C8H7N3S) in the ground state has been calculated using the Hartree-Fock and density functional method (B3LYP) with 6-31G(d) basis set. The optimized geometric bond lengths and bond angles obtained by using HF and DFT (B3LYP) show the best agreement with the experimental data. Comparison of the observed fundamental vibrational frequencies of 2-amino-5-phenyl-1,3,4-thiadiazole (C8H7N3S) and calculated results by density functional B3LYP and Hartree-Fock methods indicate that B3LYP is superior to the scaled Hartree-Fock approach for molecular vibrational problems.     

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.     

[M(N)X2]-(M=Ru, Os; X=S2C6H4, mnt)的电子结构和光谱性质的理论研究

利用密度泛函理论(DFT)中的B3LYP方法优化了氮化钌和氮化锇配合物[M(N)X2]-[M=Ru, Os; X=S2C6H4, mnt(maleonitriledithiolate)]的基态几何结构, 得到的几何参数与实验结果吻合得很好. 采用TD-DFT方法, 得到了配合物在CH3CN溶液中的激发态电子结构和电子吸收光谱. 利用SCRF方法中的CPCM模型来模拟溶剂化效应. 研究结果表明, 配合物1~4在CH3CN溶液中的吸收跃迁性质相似, 低能吸收均被指认为LMCT和LLCT的混合跃迁, 高能吸收均被指认为ILCT/LLCT跃迁.     

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

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

TD-DFT Study on the Electronic Spectrum Properties of 2,7＇-（Ethylene）- bis-8-hydroxyquinoline and Its Derivatives

The structures of 2,7＇-（ethylene）-bis-8-hydroxyquinoline and its derivatives were optimized at the ground states using ab initio HF and B3LYP methods. At the same time, the molecular structures of the first singlet excited state for 2,7＇-（ethylene）-bis-8-hydroxyqulnoline and its derivatives were optimized by CIS/6-31G（d）. The absorption and emission spectra based on the above structures were obtained by the time-dependent density functional theory （TD-DFT） by the B3LYP method with the 6-31G（d） basis set. The calculated results of luminescence originate from the electronic transition from the hydroxphenol ring of 8-hydroxyquinoline A to the pyridine ring of 8-hydroxyquinoline B. Their luminescence wave bands can be tuned by different substituents on the ligand of 8-hydroxyquinoline.     

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.     

联吡啶Ir(Ⅲ)配合物电子结构及光谱性质的理论研究

采用密度泛函理论(DFT)对配合物Ir(ppy)2(N^N)+ [ppy=2-phenylpyrine, N^N=bpy= 2,2’-bipyridine(1); N^N=H2dcbpy=4.4’-dicarboxy-2,2’-bipyridine(2), N^N=Hcmbpy=4-carboxy-4’-methyl-2,2’-bipyridine(3)] 的基态和激发态几何构型进行优化, 通过TDDFT/B3LYP方法得到这些化合物在乙腈溶液中的吸收光谱和磷光发射光谱及其跃迁性质. 研究结果表明, 化合物1 (384 nm), 2(433 nm)和3 (413 nm) 最低的吸收谱被指认为MLCT/LLCT[dIr+π(ppy)→π*(N^N)]电荷跃迁. 化合物1(486 nm), 2(576 nm)和3 (567 nm)最低的磷光发射可以描述为[dIr+π(ppy)]→[π*(N^N)]跃迁. 这是由于联吡啶配体上吸电子基团的引入, 稳定了相应的空轨道, 导致了化合物2和3的吸收和发射光谱红移. 同时, 化合物非线性光学性质的计算结果表明, 三种化合物均具有较大的一阶超极化率(β), 联吡啶配体中吸电子基团的增加, 使得分子内电子转移增强, 导致一阶超极化率增大.     

二(2-苯基-8-羟基喹啉)锌及其衍生物的电子光谱性质的含时密度泛函理论研究

采用从头算(ab initio)和密度泛函理论(DFT B3LYP)方法。对二(2-苯基-8-羟基喹啉)锌(Zn(qPh)₂)及其衍生物的基态结构进行优化，同时用ab initio HF单激发组态相互作用(CIS)法在6-31G基组上优化其最低激发单重态几何结构，用含时密度泛函理论(TD-DFT/B3LYP)及6-31G基组计算吸收和发射光谱。计算表明，该类物质电子在基态与激发态间的跃迁，主要是电子云分布由定域化向离域化的转变。吸收及发射光谱的计算值与实验值基本符合。该类化合物的电子亲和能较大，都是优良的电子传输材料，改变中心金属原子对配合物光谱性质影响不大。而羟基氧被硫原子取代后，化合物的吸收光谱产生明显红移。     

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.