A comparative electron correlation treatment in H(2)S-benzene dimer with DFT and wavefunction-based ab initio methods

The S₀ (X¹A′), T₁ (a³A″), S₁ (A¹A″), T₂ (b³A′), and S₂ (B¹A′) states of the (trans-)HONO molecule were studied by using the CASSCF and CASPT2 methods. The CASPT2(//CASPT2) adiabatic and vertical excitation energy values are in good agreement with available experimental data. The CASPT2//CASSCF potential energy curves (PECs) calculations indicate that: (i) all the five states correlate with the products of OH (X²Π) + NO (X²Π); (ii) along each of the T₁, S₁, and T₂ PECs there is a minimum followed by a transition state (barrier); and (iii) the repulsive S₂ PEC crosses the T₂, S₁, and T₁ PECs. The geometries and relative energies for the stationary points along these PECs were calculated at the CASPT2(//CASPT2) level, and the calculations predict that the barrier height value for S₁ is negligibly small (0.0018 eV).     

A theoretical estimate for quasihomopolar singlet levels in long polyenes with a spin model hamiltonian

A spin model hamiltonian is employed to estimate the energy levels of quasihomopolar singlet states in long polyenes. To this end, the singlet excitation energies of the alternating spin chain are calculated as functions of γ and N using the Hartree-Fock method developed based on the spinless fermion formalism. The result indicates that the excitation energies decrease monotonically with increasing chain length in accord with the known experimental results. Using the parameters J and γ empirically determined from the data on short polyenes, an estimation is given for the energy levels of the lowest-lying 2¹A_g states as well as those of the second ¹B_u quasihomopolar states in long polyenes. In particular, the 2¹A_g level in polyacetylene is predicted to be at 1.1 eV in good agreement with a recent photoluminescence result (1.19 eV).     

DNA环境中硒代胸腺嘧啶和腺嘌呤碱基对的激发态性质和光物理机理的理论研究

应用高精度的多态完全活化自洽场二级微扰理论方法,在量子力学/分子力学组合方法的理论框架QM(MS-CASPT2//CASSCF)/MM下,系统研究了DNA环境中2-硒和4-硒取代胸腺嘧啶和腺嘌呤碱基对(2SeT-A和4SeT-A)的最低5个电子态(S₀, S₁, S₂, T₂和T₁)的结构、性质和光物理过程. QM(MS-CASPT2//CASSCF)/MM计算揭示了DNA环境中2SeT-A和4SeT-A碱基对激发态性质和光物理过程差异性的来源,提出的机理将有助于理解DNA类似物的光物理过程,在光动力学治疗中具有潜在的应用.     

The potential energy surfaces and the radiationless dynamics of excited states of benzene and pyrazine

The potential energy surfaces of the lowest excited states of benzene and pyrazine are investigated as a function of some of the symmetry-adapted internal coordinates by means of the INDO/S method. A large stabilization of the T₂ (ππ*) state of pyrazine (≈ 0.5 eV) along the S_8b vibrational coordinate is found. The calculated potential energy in some excited states (T₁ in benzene, T₂ and S₂ in pyrazine) is a very flat function of the S_16b vibrational coordinate, leading to a crossing with the potential energy of the ground state at relatively small excess of vibrational energy (≈ 1 eV). Thus the ν_16b vibrational mode is postulated to play an important role in the radiationless relaxation to the ground states of these systems. No such crossing has been found near the “channel three” threshold of benzene.     

亚硝酸甲酯光解反应的从头算和Cl研究

用从头算和MP2方法求得亚硝酸甲酯的基态、第一和第二激发态解离为CH₃O和NO自由基的解离能分別为238.14、68.99和-183.97kJ/mol,而CH₃O和NO易于生成甲醛和硝酰。由CI方法求出的亚硝酸甲酯直接生成甲醛和硝酰的基态和激发态反应曲线表明,该反应难以按这种机理进行。因此,以上计算支持了实验提出的亚硝酸甲酯光反应生成甲醛和硝酰的两种机理中的光解离机制。     

Modeling the electronic structures of the ground and excited states of the ytterbium atom and the ytterbium dimer: A modern quantum chemistry perspective

We present a comprehensive theoretical study of the electronic structures of the Yb atom and the Yb₂ molecule, respectively, focusing on their ground and lowest-lying electronically excited states. Our study includes various state-of-the-art quantum chemistry methods such as CCSD, CCSD(T), CASPT2 (including spin-orbit coupling), and EOM-CCSD as well as some recently developed pCCD-based approaches and their extensions to target excited states. Specifically, we scan the lowest-lying potential energy surfaces of the Yb₂ dimer and provide a reliable benchmark set of spectroscopic parameters including optimal bond lengths, vibrational frequencies, potential energy depths, and adiabatic excitation energies. Our in-depth analysis unravels the complex nature of the electronic spectrum of Yb₂, which is difficult to model accurately by any conventional quantum chemistry method. Finally, we scrutinize the bi-excited character of the first excited state and its evolution along the potential energy surface.     

AuX(X=O,S)低电子态的理论研究

通常要用多参考态方法才能合理处理需考虑旋轨耦合(SOC)效应的开壳层分子如AuO和AuS的低电子态.事实上,通过选取合适的参考态,采用运动方程耦合簇方法(EOM-CC)也能计算这些分子的一些低电子态,而且EOM-CC方法是单参考态方法,使用起来比多参考态方法更加简单.本文采用最近发展的含旋轨耦合的EOM-CC计算电离能的方法(EOMIP-CC),选取对应的负离子为参考态,在CCSD级别上计算了AuO和AuS低电子态的性质.在不考虑旋轨耦合时,通过比较EOMIP-CCSD和EOMIP-CCSDT的结果考察EOMIPCCSD的精度.此外,与EOMIP-CCSDT结果相比,如果自旋污染较为显著而且T1的模较大时,UCCSD(T)方法对能量最低的某一特定对称性的电子态的所对应的电离能误差约为0.1-0.15 eV.在考虑了旋轨耦合效应后,我们的方法得到的键长和振动频率与实验值吻合较好.另一方面,虽然EOMIP-SOC-CCSD高估了能量较高的2Δ3/2态、2Σ+1/2态和2Π1/2态的能量,但是对于其它能量更低的电子态,它们的能量与已有实验值误差在0.2 eV左右.这显示我们所用的含SOC的EOMIP-CCSD方法对原本需要用多参考态方法才能处理的AuO和AuS低电子态能给出可靠的结果.     

A TDDFT study of the low-lying excitation energies of polycyclic cinnolines and their carbocyclic analogues

The low-lying excited states of 21 compounds of polycyclic cinnoline, monoaza-hydrocarbon and their corresponding polycyclic aromatic hydrocarbon (PAH) analogues have been investigated in the framework of time-dependent density functional theory (TDDFT). The gradient corrected BLYP, hybrid B3LYP and B3P86 functionals together with 6-31G (d) basis set have been used. The hybrid-type B3LYP and B3P86 systematically overestimate the excitation energies for states with dominating ionic character (corresponding to group III band), with a mean absolute deviation (MAD) of 0.33 eV (B3LYP) and 0.34 eV (B3P86), respectively. However, they can accurately predict the excitation energies for states with covalent character (corresponding to group II and I bands). The MAD for group II (B3LYP: 0.05 eV; B3P86: 0.05 eV) and I bands (B3LYP: 0.12 eV; B3P86: 0.12 eV) are significantly smaller. The BLYP outperforms B3LYP and B3P86 for group III band (MAD: 0.09 eV), but has a worse performance for group II (MAD: 0.15 eV) and I (MAD: 0.13 eV) bands. Comparison of the lowest-lying excited states for polycyclic cinnolines with those of the corresponding PAH analogues, the first excited states of polycyclic cinnolines mainly result from n→π₀* transitions. Therefore, in non-polar solvents, the spectra of some polycyclic cinnolines exhibit an additional absorption band at longer wavelengths.     

A CASSCF/CASPT2 study on the low-lying excited states of HSiCN, HSiNC and their ions

Equilibrium geometries of low-lying electronic states of cyanosilylene (HSiCN), isocyanosilylene (HSiNC), and their ions have been investigated using the complete active space self-consistent field (CASSCF) approach. The harmonic vibrational frequencies on the optimized geometries were calculated using the multiconfiguration linear response (MCLR) method. Taking the further correlation effects into account, the complete active space perturbation theory of second-order (CASPT2) was carried out for the energetic correction. The CASPT2 calculations have been performed to obtain the vertical excitation energies of selected low-lying excited states of HSiCN and HSiNC. Computed results show that the singlet-triplet splittings are calculated to be 0.99 and 1.30 eV for HSiCN and HSiNC, respectively. The vertical excitation energies of the lowest singlet and triplet excited states in HSiCN are lower than those in HSiNC. The first vertical ionization energy of HSiCN (10.04 eV) is higher than that of HSiNC (9.97 eV). The ground-state adiabatic electron affinities are found to be rather high, and the value of HSiCN (1.85 eV) higher than that of HSiNC (1.52 eV). The existences of dipole-bound excited negative ion states have been discovered within HSiCN and HSiNC.     

Equation-of-motion coupled-cluster study on exciton states of polyethylene with periodic boundary condition

Equation-of-motion coupled cluster with singles and doubles (EOM-CCSD) method has been applied to exciton states of polyethylene using ab initio crystal Hartree-Fock method with one-dimensional periodic boundary condition. Full transformation of two-electron integrals from atomic-orbital basis to crystal-orbital basis has been performed for EOM-CCSD calculations. In order to make transformed integrals to have correct properties of translational symmetry, a lattice summation scheme has been proposed. The EOM-CCSD excitation energies have been obtained for the lowest singlet and triplet exciton states of polyethylene. The excitation energies converge with system size much faster than oligomer calculations using n-alkanes. Quasiparticle energy-level calculations by second-order many-body perturbation theory and by solving the inverse Dyson equation have also been performed to obtain exciton binding energies. Basis set dependencies on excitation energy, quasiparticle band gap, and exciton binding energy have been investigated. At the 6-31+G level, the excitation energy of the lowest singlet-exciton state and its binding energy are calculated to be 8.1 and 3.2 eV, respectively. The calculated excitation energy is well comparable with the corresponding experimental value, 7.6 eV.     

A theoretical study of the excited electronic states of the SiCN system

Ab initio methods have been used to study the lowest-lying electronic states of the SiCN radical, which has two stable linear isomers in its electronic ground state, SiCN and SiNC. Vertical excitation energies and oscillator strengths have been computed for a number of states lying up to 8 eV. The geometries of the lowest-lying doublet and quartet states have been determined. The lowest-lying excited doublet state of SiNC (1²Σ⁺, 4.0 eV) arises from a HOMO–LUMO excitation (3π → 10σ), although the 1²Δ state (9σ → 3π) is very close in energy. In the case of the SiCN isomer the lowest excited state is 1²Δ, which arises from an excitation from the highest occupied σ orbital into the HOMO (9σ → 3π) and lies 3.6 eV above the ground state. SiCN should present very strong absorptions at 4.9 and 6.1 eV whereas SiNC should have relatively strong absorptions in the region of 5.7–5.9 eV. The smallest adiabatic energy gaps with respect to the ground state of SiNC and SiCN are very close (about 2.8 eV) and the excited state is the same 1²A′, which has angular equilibrium geometries for both isomers. We have determined accurate values for enthalpies of formation of the two linear doublet forms and .     

Dissociation reactions of the vinyl radical in the AA″ state

In the present theoretical work we have explored mechanisms of dissociation reactions of the vinyl radical in the A²A″ state (C₂H₃ (A²A″)) and examined possible pathways for nonadiabatic dissociation of C₂H₃ (A²A″) into C₂H₂ (X¹Σ_g⁺). In the calculations we used the complete active space self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) methods in conjunction with the cc-pVDZ and cc-pVTZ basis sets. Mechanisms for the following three dissociation channels of C₂H₃ in the A²A″ state were explored: (1) C₂H₃ (A²A″) → C₂H₂ (trans, ³A_u) + H, (2) C₂H₃ (A²A″) → C₂H₂ (cis, ³A₂) + H, and (3) C₂H₃ (A²A″) → H₂CC (³A₂) + H. The CASSCF and CASPT2 potential energy curve calculations for the C₂H₃ (A²A″) dissociation channels (1)–(3) indicate that there is neither transition state nor intermediate for each of the channels. At the CASPT2//CASSCF/cc-pVTZ level, the dissociation energies for channels (1)–(3) are predicted to be 84.3, 91.1, and 86.9 kcal/mol, respectively. For a recently observed nonadiabatic dissociation of C₂H₃ (A²A″) into C₂H₂ (X¹Σ_g⁺) + H [J. Chem. Phys. 111 (1999) 3783], two previously suggested internal conversion (IC) pathways were examined based on our CASSCF and CASPT2 calculations. Our preliminary CASSCF and CASPT2 calculations indicate that the assumed IC pathway via the twisted C₂H₃ (A²A^′) structure might be feasible. The CASSCF/cc-pVTZ geometry optimization and frequency analysis calculations were performed for the four C_2v bridge structures in the ²B₂, ²A₂, ²B₁, and ²A₁ states along the pathways of the 1²A^′ (X²A^′), 1²A″ (A²A″), 2²A″, and 2²A^′ states of C₂H₃, respectively, and the CASPT2//CASSCF/cc-pVTZ energetic results indicate that the assumed IC pathway, via a C_2v (²A₂) structure and then ²A₂/²A₁ surface crossing, be not feasible since at their excitation wavelengths (327.4 and 366.2 nm) the C_2v (²A₂) structure could not be accessed.     

Locally correlated equation-of-motion coupled cluster theory for the excited states of large molecules

We report an extension of the local correlation concept to electronically excited states via the equation-of-motion coupled cluster singles and doubles (EOM-CCSD) method. We apply the same orbital domain structure used successfully for ground-state CCSD by Werner and co-workers and find that the resulting localized excitation energies are in error generally by less than 0.2 eV relative to their canonical EOM-CCSD counterparts, provided the basis set is flexible and includes Rydberg-like functions. In addition, we account for weak-pair contributions efficiently using a correction to local-EOM-CCSD transition energies based on the perturbative (D) correction used with configuration interaction singles (CIS).     

Predictions of the geometries and fluorescence emission energies of oxyluciferins

The complete active space self-consistent field (CASSCF) method and multiconfigurational second-order perturbation theory (CASPT2) have been used to study the structures and spectra of oxyluciferins (OxyLH2). The ground and lowest-lying singlet excited states geometries have been optimized using CASSCF. CASPT2 has been used to predict relaxed emission energies. The focus is on the lowest-lying singlet excited states of the anionic keto and enol forms of OxyLH2(-1) at the optimized excited-state geometries. The planar keto and enol forms of OxyLH2(-1) are minima on both the S0 and the S1 potential energy surfaces. The twisted keto and enol forms of OxyLH2(-1) are transition states on the S0 and S1 potential energy surfaces. The S1 --> S0 fluorescence emission energies are in the range of 54.2-58.4 kcal/mol for the anionic planar keto forms of OxyLH2, and in the range of 55.7-63.2 kcal/mol for the anionic enol forms of OxyLH2. S0 and S1 potential energy surfaces and thus are not implicated in the emission spectra in the gas phase.     

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.     

The triplet state of cytosine and its derivatives: electron impact and quantum chemical study

The excitation of the lowest electronic states and vibrational excitation of cytosine (C) have been studied using electron energy loss spectroscopy (EELS, 0-100 eV) with angular analysis. The singlet states have been found to be in good agreement with UV-VIS absorption results on sublimed films, slightly blueshifted by about 0.1 eV. The EEL spectra recorded at residual energy below 2 eV show clear shoulders at energy losses of 3.50 and 4.25 eV (+/-0.1 eV). They are assigned to the lowest triplet electronic states of cytosine. Energies and molecular structures of the lowest-lying triplet state of C and its methylated and halogenated 5-X-C, 6-X-C, and 5-X, 6-X-C substituted derivatives (X=CH3, F, Cl, and Br) have been studied using quantum chemical calculations with both molecular orbital and density functional methods, in conjunction with the 6-311++G(d,p), 6-311++G(3df,2p), and aug-cc-pVTZ basis sets. The triplet-singlet energy gap obtained using coupled-cluster theory [CCSD(T)] and density functional theory (DFT) methods agrees well with those derived from EELS study. The first C's vertical triplet state is located at 3.6 eV, in good agreement with experiment. The weak band observed at 4.25 eV is tentatively assigned to the second C's vertical triplet excitation. For the substituted cytosines considered, the vertical triplet state is consistently centered at 3.0-3.2 eV above the corresponding singlet ground state but about 1.0 eV below the first excited singlet state. Geometrical relaxation involving out-of-plane distortions of hydrogen atoms leads to a stabilization of 0.6-1.0 eV in favor of the equilibrium triplet. The lowest-lying adiabatic triplet states are located at 2.3-3.0 eV. Halogen substitution at both C(5) and C(6) positions tends to reduce the triplet-singlet separations whereas methylation tends to enlarge it. The vibrational modes of triplet cytosine and the ionization energies of substituted derivatives were also evaluated.     

Triplet states of ketene by trapped electron spectroscopy

The threshold electron impact excitation spectrum of ketene is reported. The spectrum is interpreted in comparison with results of an ab initio frozen core calculation. Triplet states (³A₂, ³A₁, ³B₁) are observed at 3.8, 5.0 and 5.8 eV excitation energy.     

Theoretical prediction of the electronic excited states and resonance Raman intensities in formamide from coupled cluster calculations

Electronic excitations and the resonance Raman spectrum of formamide were obtained from ab initio electron correlation calculations using the equation of motion coupled cluster (EOM-CCSD) method. Interpretation of the UV spectrum on the basis of calculated vertical excitation energies and oscillator strengths accounts for all experimental bands previously assigned. Our assignment, however, suggests an additional Rydberg band at about 7.4 eV which may be hidden under the main absorption. We also show that the Rydberg states appear pairwise, corresponding to n and π hole states, respectively. Using analytic derivative techniques, derivatives of the excited state energies with respect to normal coordinates of the ground state were calculated. Approximate resonance Raman intensities have been determined.