面向有机光电材料的电子激发态结构与过程的计算方法

共轭聚合物与有机分子材料中的电子激发结构与过程决定了材料的光电功能：根据Kasha规则,低能级激发态的排序决定能否发光;最低激发态至基态的辐射跃迁与无辐射跃迁之间的竞争决定了发光效率,后者主要由非绝热耦合（声子作用）决定;电荷激发态载体的传输由电子分布与振动耦合或杂质和无序的散射弛豫过程决定．本文针对有机功能材料的发光性能,介绍两种理论方法的研究进展,即可用于计算共轭聚合物激发态结构的量子化学密度矩阵重整化群方法和计算发光效率的多模耦合无辐射跃迁速率方法．这些方法被应用于有机功能材料的性能预测和分子设计中．     

An intermolecular perturbation approach to hydrogen bondings in systems in the ground and excited electronic states

The intermolecular interaction energy for binary systems in the ground and excited electronic states was partitioned into the Coulomb, exchange-repulsion, induction, dispersion and charge-transfer interaction terms by the perturbation expansion method. The various interaction terms were evaluated for the hydrogen bondings in (HF)₂, (H₂O)₂, (CH₃OH)₂, (RCOOH)₂, and HF·H₂O in various geometrical configurations. It has been found that the Coulombic interaction plays a dominant role in the stability of these hydrogen bonded systems. The method was further applied to the HCOOH·H₂O codimer in both the ground and excited singlet electronic states. The results were in accord with the well-known water solvent effects on the shifts of absorption spectral bands.     

Conformational stability, optimized geometries, vibrational and electronic spectra of methacryloyl bromide in ground and excited electronic states

In order to understand conformational isomerism in methacryloyl bromide (MABR) in the ground (S(0)) and the first excited (S(1)) electronic states and to interpret the vibrational and electronic spectra of its conformers in the S(0) state, quantum mechanical calculations using Density Functional Theory (DFT) and RHF methods with extended basis sets 6-31G, 6-31G** and 6-311+G(d,p) have been conducted. In RHF calculations, electron correlation effects have been included at the M?ller-Plesset MP2 level. It is inferred that in both the electronic states the molecule may exist in two isomeric forms-s-trans and s-cis; the former being more stable than the later by about 1.629 kcal mol(-1) in the S(0) state and by about 2.218 kcal mol(-1) in the S(1) state. Electronic transition tends to increase the s-trans/s-cis and s-cis/s-trans, rotational barriers from 7.059 kcal mol(-1) (2468.1 cm(-1)) and 5.428 kcal mol(-1) (1897.8 cm(-1)) in S(0) state to 23.594 kcal mol(-1) (8249.4 cm(-1)) and 21.376 kcal mol(-1) (7473.9 cm(-1)) in the S(1) state. Completely optimized geometries of the two conformers in S(0) state reveal that while there is no significant difference in their bond lengths, some of the bond angles associated with COBr group are appreciably different. Electronic excitation tends to change both the bond lengths and bond angles. Based on suitably scaled DFT and RHF results obtained from the use of 6-31G** and 6-311+G(d,p) basis sets, a complete assignment is provided to the fundamental vibrational bands of both the s-trans and s-cis conformers in terms of frequency, form and intensity of vibrations and potential distribution across the symmetry coordinates in the S(0) state and a comparison has been made with experimental assignments. A theoretical prediction of the electronic transitions in the near UV-region in the two conformers and their tentative assignment has been provided on the basis of CI level calculations using 6-31G basis set.     

Ab initio studies of ground and excited electronic states of MgAr,CdAr, and BeAr

Summary The ground state (X ¹⁺) and several excited state (A ³,c ³+,C ¹,D ¹⁺, andE ³⁺) potential energy surfaces for the diatomic molecules MgAr, CdAr, and BeAr have been computed using complete active space self-consistent field (CASSCF) wavefunctions and valence double- and triplezeta quality basis sets augmented with polarization and diffuse functions. Pump-and-probe laser experiments have examined the quenching, of excited singlet states of metal-rare gas complexes such as CdXe to produce triplets that dissociate to³ P _Jmetal atoms. This quenching, which is detected for CdXe but not for CdAr or MgAr, is thought to occur via a crossing or strong coupling of a repulsive triplet curve correlating to the underlying³ P state of the metal, with an attractive singlet curve that correlates to the higher¹ P state of the metal. The present work indicates that the attractiveC ¹ and repulsivec ³⁺ curves of MgAr and CdArdo not intersect in the energetically accessible region of theC ¹ surface, unlike the corresponding curves for the CdXe diatom. These data are consistent with the absence of³ P _J Cd atoms in the MgAr and CdAr experiments, respectively. However, an alternative quenching mechanism involving vibronic coupling between theC ¹ vibrational eigenstates and the continuum eigenstates of the underlying repulsive³⁺ surface may be operative; this possibility is examined qualitatively and predicted to be unlikely for MgAr (due to small spin-orbit coupling) and CdAr (due to unfavorable vibronic factors). BeAr, which has yet to be probed experimentally, is predicted to be bound by 770 and 900 cm^–1 in theD ¹⁺ state (which has metal 2s2p character) and theE ³⁺ state (which has Rydberg metal 2s3s character), respectively, and to display interesting potential curve intersections.Dedicated to Prof. Klaus Ruedenberg     

Characterization of singlet ground and low-lying electronic excited states of phosphaethyne and isophosphaethyne

The singlet ground ((approximate)X(1)Sigma1+) and excited (1Sigma-,1Delta) states of HCP and HPC have been systematically investigated using ab initio molecular electronic structure theory. For the ground state, geometries of the two linear stationary points have been optimized and physical properties have been predicted utilizing restricted self-consistent field theory, coupled cluster theory with single and double excitations (CCSD), CCSD with perturbative triple corrections [CCSD(T)], and CCSD with partial iterative triple excitations (CCSDT-3 and CC3). Physical properties computed for the global minimum ((approximate)X(1)Sigma+HCP) include harmonic vibrational frequencies with the cc-pV5Z CCSD(T) method of omega1=3344 cm(-1), omega2=689 cm(-1), and omega3=1298 cm(-1). Linear HPC, a stationary point of Hessian index 2, is predicted to lie 75.2 kcal mol(-1) above the global minimum HCP. The dissociation energy D0[HCP((approximate)X(1)Sigma+)-->H(2S)+CP(X2Sigma+)] of HCP is predicted to be 119.0 kcal mol(-1), which is very close to the experimental lower limit of 119.1 kcal mol(-1). Eight singlet excited states were examined and their physical properties were determined employing three equation-of-motion coupled cluster methods (EOM-CCSD, EOM-CCSDT-3, and EOM-CC3). Four stationary points were located on the lowest-lying excited state potential energy surface, 1Sigma- -->1A", with excitation energies Te of 101.4 kcal mol(-1) (1A"HCP), 104.6 kcal mol(-1)(1Sigma-HCP), 122.3 kcal mol(-1)(1A" HPC), and 171.6 kcal mol(-1)(1Sigma-HPC) at the cc-pVQZ EOM-CCSDT-3 level of theory. The physical properties of the 1A" state with a predicted bond angle of 129.5 degrees compare well with the experimentally reported first singlet state ((approximate)A1A"). The excitation energy predicted for this excitation is T0=99.4 kcal mol(-1) (34 800 cm(-1),4.31 eV), in essentially perfect agreement with the experimental value of T0=99.3 kcal mol(-1)(34 746 cm(-1),4.308 eV). For the second lowest-lying excited singlet surface, 1Delta-->1A', four stationary points were found with Te values of 111.2 kcal mol(-1) (2(1)A' HCP), 112.4 kcal mol(-1) (1Delta HPC), 125.6 kcal mol(-1)(2(1)A' HCP), and 177.8 kcal mol(-1)(1Delta HPC). The predicted CP bond length and frequencies of the 2(1)A' state with a bond angle of 89.8 degrees (1.707 A, 666 and 979 cm(-1)) compare reasonably well with those for the experimentally reported (approximate)C(1)A' state (1.69 A, 615 and 969 cm(-1)). However, the excitation energy and bond angle do not agree well: theoretical values of 108.7 kcal mol(-1) and 89.8 degrees versus experimental values of 115.1 kcal mol(-1) and 113 degrees. of 115.1 kcal mol(-1) and 113 degrees.     

Linear-scaling localized-density-matrix method for the ground and excited states of one-dimensional molecular systems

A linear-scaling localized-density-matrix (LDM) method is developed to evaluate the ground-state reduced single-electron density matrices of one-dimensional molecular systems. The new method may be combined with the existing linear-scaling LDM method for the excited states (Yokojima and Chen, Chem. Phys. Lett. 292 (1998) 379), and thus leads to a linear-scaling calculation method for the properties of both the ground and excited states. The combined method is applied to the polyacetylene oligomers and the linear-scaling of the total computational time is clearly demonstrated.     

Alternated quinoid/aromatic units in terthiophenes building blocks for electroactive narrow band gap polymers. Extended spectroscopic, solid state, electrochemical, and theoretical study

We report here the synthesis of three novel pi-conjugated heterocyclic mixed trimers that contain two electron-donating 3,4-ethylenedioxy-2-thienyl (EDOT) units covalently attached to a central proquinoid electron-accepting thienopyrazine moiety (two of these narrow-HOMO-LUMO gap D-A-D compounds also bear hexyl side chains attached either to the outermost alpha positions of the EDOT end rings or to the beta positions of the pyrazine fused ring). The modification of the terthiophene structure upon EDO, pyrazine, and hexyl substitutions has been treated in detail with spectroscopic and theoretical arguments. Solid-state properties reveal the occurrence of short intramolecular contacts between heteroatoms of adjacent rings. The analysis of the structure of the pi-conjugated backbone of each molecule is consistent with a partial quinoid-like pattern which partially reverts to be subtly more aromatic depending on the topology of the positive inductive effect of the hexyl chains. This quinoidization is a consequence of the appearance of a D(EDOT)-->A(PyT)<--D(EDOT) intramolecular charge transfer which further polarizes the structure. The same chemical concepts have been applied to address their electrochemical behavior. The three mixed trimers exhibit amphoteric properties due to the combination of electron acceptor and donor groups. Given their relative low HOMO-LUMO energy gap, these trimers promise to be good candidates for obtaining polymers with significant low energy gap combining electroactivity.     

Vibrational analysis of excited and ground electronic states of all-trans retinal protonated Schiff-bases

We report on vibrational coherence dynamics in excited and ground electronic states of all-trans retinal protonated Schiff-bases (RPSB), investigated by time-resolved Degenerate Four-Wave-Mixing (DFWM). The results show that wave packet dynamics in the excited state of RPSB consist of only low-frequency (<800 cm(-1)) modes. Such low-frequency wave packet motion is observed over a broad range of detection wavelengths ranging from excited state absorption (～500 nm) to stimulated emission (>600 nm). Our results indicate that low-frequency coherences in the excited state are not activated directly by laser excitation but rather by internal vibrational energy redistribution. This is supported by the observation that similar coherence dynamics are not observed in the electronic ground state. Challenging previous experimental results, we show that the formation of low-frequency coherence dynamics in RPSB does not require significant excess vibrational energy deposition in the excited state vibrational manifolds. Concerning ground state wave packet dynamics, we observe a set of high-frequency (>800 cm(-1)) modes, reflecting mainly single and double bond stretching motion in the retinal polyene-chain. Dephasing of these high-frequency coherences is mode-dependent and partially differs from analogous vibrational dephasing of the all-trans retinal chromophore in a protein environment (bacteriorhodopsin).     

Geometric and electronic structure of ground and excited states of groupVA diatomics. A theoretical LCGTO-MP-LSD study

LCGTO-MP-LSD calculation was performed for the ground and several low-lying excited states of homo- (N₂, P₂, As₂, and Sb₂) and hetero-nuclear (PN, AsN, AsP, AsSb, SbN, and SbP) groupVA diatomics. For all the systems the ground state is found to be¹Σ⁺. For N₂ and P₂, the¹Σ _g ⁺ ground state is followed by the³Σ _u ⁺ ,³Π _g ,³Δ _u ,¹Π _g , and¹Δ _u low-lying exited states while for As₂ the order is found to be³Σ _u ⁺ ,³Δ _u ,³Π _g ,¹Δ _u ,¹Π _g . Finally for Sb₂ the relative stability of excited states is³Σ _u ⁺ ,³Δ _u ,¹Δ _u ,³Π _g ,¹Π _g . For the hetero-nuclear diatomics the¹Σ⁺ ground state is, in the case of PN, AsN, AsP, SbN, and SbP, followed by the³Σ⁺,³Δ,³Π,¹Π and¹Δ low-lying excited states while for the AsSb diatomic an inversion of stability of the two last singlets occurs. The calculated spectroscopic parameters (Re, ωe, andDe) are in good agreement with all the available experimental results while, theTe values are overestimated by about 0.5 eV. Mulliken population analysis shows that both homo- and hetero-nuclear groupVA diatomics are essentially triple bonded systems.     

Reactivity and electronic structure of the ground and excited states of oxazole, 2-phenyloxazole, and 2-phenylthiazole

The relationship of the electronic structure, of oxazole, 2-phenyloxazole, and 2-phenylthiazole to the properties of the electronic excited states and transitions is examined. Spectral properties of these compounds in the free state (no effect from external perturbations) and in complexes with the proton and aprotic acids are measured and calculated by quantum-chemical PPP/S (-approximation), and INDO/S (sp-basis) methods. Features of the electronic excitation of the atoms and the vibronic interaction of bonds in the singlet and triplet states are examined for a change of the various structural forms of azoles, which determine their spectral fluorescence properties and reactivity. Possible direction control of reactions and optimized syntheses of new compounds with given properties are discussed based on a study of the properties of the ground and excited states.Russian People's Friendship University, Moscow 117302 N. D. Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, Moscow 117913 Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 7, pp. 948–974, July, 1999.     

Decomposition of nitramine energetic materials in excited electronic states: RDX and HMX

Ultraviolet excitation (8-ns duration) is employed to study the decomposition of RDX (1,3,5-trinitro-1,3,5-triazacyclohexane) and HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane) from their first excited electronic states. Isolated RDX and HMX are generated in the gas phase utilizing a combination of matrix-assisted laser desorption and supersonic jet expansion techniques. The NO molecule is observed as one of the initial dissociation products by both time-of-flight mass spectroscopy and laser-induced fluorescence spectroscopy. Four different vibronic transitions of NO are observed: A (2)Sigma(v(') = 0)<--X (2)Pi(v(") = 0,1,2,3). Simulations of the NO rovibronic intensities for the A<--X transitions show that dissociated NO from RDX and HMX is rotationally cold (approximately 20 K) and vibrationally hot (approximately 1800 K). Another potential initial product of RDX and HMX excited state dissociation could be OH, generated along with NO, perhaps from a HONO intermediate species. The OH radical is not observed in fluorescence even though its transition intensity is calculated to be 1.5 times that found for NO per radical generated. The HONO intermediate is thereby found not to be an important pathway for the excited electronic state decomposition of these cyclic nitramines.     

Vibrational co-assignment of the calculated frequencies in the ground and two lowest excited electronic states

The feasibility of co-assignments of the vibrational frequencies in the ground and T₁ and S₁ excited electronic states has been demonstrated for trans-C₂O₂F₂. Matrices analogous to the Duschinsky matrix were used to juxtapose the vibrational frequencies of this molecule calculated at the CASPT2/cc-pVTZ level in the ground S₀ and excited triplet T₁ and singlet S₁ electronic states. The calculations suggest that the calculated CC and CF stretching frequencies of trans-C₂O₂F₂ in these three electronic states should be mutually reassigned in comparison with the previous interpretation.     

Pump-probe spectroscopy of molecules driven by infrared field in both ground and excited electronic states

Pump-probe spectra of molecules driven by strong infrared (IR) field in both ground and excited states are studied theoretically. The role of the final state interaction becomes important when pump and probe pulses overlap, and the Rabi frequency is comparable with the lifetime broadening of the excited state and the duration of the pump pulse. Our theoretical approach is applied to x-ray absorption of nitrogen molecule and valence photoionization of carbon monoxide. It is shown that IR-x-ray pump-probe spectroscopy can directly evidence the delocalization of core hole.     

The ground and two lowest-lying singlet excited electronic states of copper hydroxide (CuOH)

Various ab initio methods, including self-consistent field (SCF), configuration interaction, coupled cluster (CC), and complete-active-space SCF (CASSCF), have been employed to study the electronic structure of copper hydroxide (CuOH). Geometries, total energies, dipole moments, harmonic vibrational frequencies, and zero-point vibrational energies are reported for the linear 1Sigma+ and 1Pi stationary points, and for the bent ground-state X 1A', and excited-states 2 1A' and 1 1A". Six different basis sets have been used in the study, Wachters/DZP being the smallest and QZVPP being the largest. The ground- and excited-state bending modes present imaginary frequencies for the linear stationary points, indicating that bent structures are more favorable. The effects of relativity for CuOH are important and have been considered using the Douglas-Kroll approach with cc-pVTZ/cc-pVTZ_DK and cc-pVQZ/cc-pVQZ_DK basis sets. The bent ground and two lowest-lying singlet excited states of the CuOH molecule are indeed energetically more stable than the corresponding linear structures. The optimized geometrical parameters for the X 1A' and 1 1A" states agree fairly well with available experimental values. However, the 2 1A' structure and rotational constants are in poor agreement with experiment, and we suggest that the latter are in error. The predicted adiabatic excitation energies are also inconsistent with the experimental values of 45.5 kcal mol(-1) for the 2 1A' state and 52.6 kcal mol(-1) for the 1 1A" state. The theoretical CC and CASSCF methods show lower adiabatic excitation energies for the 1 1A" state (53.1 kcal mol(-1)) than those for the corresponding 2 1A' state (57.6 kcal mol(-1)), suggesting that the 1 1A" state might be the first singlet excited state while the 2 1A' state might be the second singlet excited state.     

Internal rotation potential functions of an acryloyl fluoride molecule in the ground (S 0) and excited (S 1) electronic states

Journal of Structural Chemistry - Structural parameters of trans- and cis-isomers of an acryloyl fluoride molecule in the ground (S 0) and excited (S 1) electronic states are determined. The...     

Theoretical study of isomerization and dissociation of acetylene dication in the ground and excited electronic states

Ab initio calculations employing the configuration interaction method including Davidson's corrections for quadruple excitations have been carried out to unravel the dissociation mechanism of acetylene dication in various electronic states and to elucidate ultrafast acetylene-vinylidene isomerization recently observed experimentally. Both in the ground triplet and the lowest singlet electronic states of C2H2(2+) the proton migration barrier is shown to remain high, in the range of 50 kcal/mol. On the other hand, the barrier in the excited 2 3A" and 1 3A' states decreases to about 15 and 34 kcal/mol, respectively, indicating that the ultrafast proton migration is possible in these states, especially, in 2 3A", even at relatively low available vibrational energies. Rice-Ramsperger-Kassel-Marcus calculations of individual reaction-rate constants and product branching ratios indicate that if C2H(2)2+ dissociates from the ground triplet state, the major reaction products should be CCH+(3Sigma-)+H+ followed by CH+(3Pi)+CH+(1Sigma+) and with a minor contribution (approximately 1%) of C2H+(2A1)+C+(2P). In the lowest singlet state, C2H+(2A1)+C+(2P) are the major dissociation products at low available energies when the other channels are closed, whereas at Eint>5 eV, the CCH+(1A')+H+ products have the largest branching ratio, up to 70% and higher, that of CH+(1Sigma+)+CH+(1Sigma+) is in the range of 25%-27%, and the yield of C2H++C+ is only 2%-3%. The calculated product branching ratios at Eint approximately 17 eV are in qualitative agreement with the available experimental data. The appearance thresholds calculated for the CCH++H+, CH++CH+, and C2H++C+ products are 34.25, 35.12, and 34.55 eV. The results of calculations in the presence of strong electric field show that the field can make the vinylidene isomer unstable and the proton elimination spontaneous, but is unlikely to significantly reduce the barrier for the acetylene-vinylidene isomerization and to render the acetylene configuration unstable or metastable with respect to proton migration.     

Ab initio potential energy surface and excited vibrational states for the electronic ground state of Li2H

A 285-point multi-reference configuration-interaction involving single and double excitations (MRS-DCI) potential energy surface for the electronic ground state of Li₂H is determined by using 6-311G (2df, 2pd) basis set. A Simons-Parr-Finlan polynomial expansion is used to fit the discrete surface with a X² of 4.64 × 10^-6. The equilibrium geometry occurs at R_e =0.172 nm and <LiHLi =94.10. The dissociation energy for reaction Li₂H(²A)⇑ Li₂(¹⌆_g)+H(²S) is 243.910 kJ/mol. and that for reaction Li₂H(²A)⇑HLi(¹be)+Li(²S) is 106.445 kJ/mol. The inversion barrier height is 50.388 kJ/mol. The vibrational energy levels are calculated using the discrete variable representation (DVR) method. Project supported by the National Natural Science Foundation of China (grant No. 29673029) and by the Special Doctoral Research Foundation of the State Education Commission of China.