Electronic excitation spectrum of thiophene studied by symmetry-adapted cluster configuration interaction method

A multiphoton ionisation photoelectron imaging study of thiophene

Thiophene is a prototype for the excited state photophysics that lies at the heart of many technologies within the field of organic electronics. Here, we report a multiphoton ionisation photoelectron imaging study of gas-phase thiophene using a range of photon energies to excite transitions from the ground electronic state to the first two electronically excited singlet states, from the onset of absorption to the absorption maximum. Analysis of the photoelectron spectra and angular distributions reveal features arising from direct photoionisation from the ground electronic state, and resonance-enhanced photoionisation via the electronically excited singlet states. The first two ionisation energies from the ground electronic state were confirmed to be 8.8 eV (adiabatic) and 9.6 eV (vertical). The ionisation energies from the first two electronically excited singlet states were found to be 3.7 eV (adiabatic) and 4.4 eV (vertical).

A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Medium Sized Molecules

Following our previous work focusing on compounds containing up to 3 non-hydrogen atoms [J. Chem. Theory Comput. 2018, 14, 4360-4379], we present here highly accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms: acetone, acrolein, benzene, butadiene, cyanoacetylene, cyanoformaldehyde, cyanogen, cyclopentadiene, cyclopropenone, cyclopropenethione, diacetylene, furan, glyoxal, imidazole, isobutene, methylenecyclopropene, propynal, pyrazine, pyridazine, pyridine, pyrimidine, pyrrole, tetrazine, thioacetone, thiophene, thiopropynal, and triazine. To obtain these energies, we use equation-of-motion/linear-response coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ) and selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational n-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/aug-cc-pVQZ vertical excitation energies as well as CC3/aug-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of ±0.04 eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse, and accurate set of more than 200 highly accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, n → π*, π → π*, ...). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, and NEVPT2. The results of these benchmarks are compared to the available literature data.

Qualitatively Incorrect Features in the TDDFT Spectrum of Thiophene-Based Compounds

Ab initio molecular electronic structure computations of thiophene-based compounds constitute an active field of research prompted by the growing interest in low-cost materials for organic electronic devices. In particular, the modeling of electronically excited states and other time-dependent phenomena has moved toward the description of more realistic albeit challenging systems. We demonstrate that due to its underlying approximations, time-dependent density functional theory predicts results that are qualitatively incorrect for thiophene and thienoacenes, although not for oligothiophene chains. The failure includes spurious state inversion and excitation characters, wrong distribution of oscillator strengths and erroneous potential energy surfaces. We briefly analyze possible origins of this behavior and identify alternative methods that alleviate these problems.

Charge transfer or biradicaloid character: assessing TD-DFT and SAC-CI for squarylium dye derivatives

Le Bahers's diagnostic indexes have been applied to study the electronic transition character of symmetrical squarylium dye derivatives through TD-DFT and SAC-CI methods.

A theoretical exploration on electronically excited states of protonated furan and thiophene

MPE profiles of the ground and S1 excited states of protonated furan along C–O bond-stretching reaction coordinate, computed at the RI-CC2/aug-cc-pVDZ level of theory.

Continuous process of the vacuum ultraviolet- (VUV-) photochemical oxidation of thiophene in the gas phase

Based on kinetic data and intermediate products of oxidation, reaction pathways of the oxidative degradation of gaseous thiophene are discussed.

Origin of Near-Infrared Absorption and Large Second Hyperpolarizability in Oxyallyl Diradicaloids: A Three-State Model Approach

Bis(benzofuranonyl)methanolate (BM4i4i) dye and croconate dyes (derivatives of oxyallyl molecules) in general are known to have intense transitions in the near-infrared (NIR) region, indicating small transition energies and large transition dipole moments. These molecules have been reported in the literature to have very large resonant third-order nonlinear optical (NLO) susceptibilities and molecular hyperpolarizabilities (gamma). In this work we investigate using density functional theory (DFT)/ab initio/symmetry adapted cluster-configuration interaction (SAC-CI) techniques the oxyallyl substructure and attribute the NIR transition and the NLO activity to this substructure, which is common in all these molecules. Using valence bond (VB) theory, an analysis of a three-state model of this substructure is carried out. It is seen that the mixture of an intermediate diradical character and some zwitterionic character in the molecule and a large coupling between these two VB resonance forms is responsible for large gamma values. This can be used as a design principle for increasing NLO activity in oxyallyl derivatives.

Theoretical study of the vertical excited states of benzene, pyrimidine, and pyrazine by the symmetry adapted cluster—Configuration interaction method

The ground state and the excited states of benzene, pyrimidine, and pyrazine have been examined by using the symmetry adapted cluster-configuration interaction (SAC-CI) method. Detailed characterizations and the structures of the absorption peaks in the vacuum ultraviolet (VUV), low energy electron impact (LEEI), and electron energy loss (EEL) spectra were theoretically clarified by calculating the excitation energy and the oscillator strength for each excited state. We show that SAC-CI has the power to well reproduce the electronic excitation spectra (VUV, LEEI, and EEL) simultaneously to an accuracy for both the singlet and the triplet excited states originated from the low-lying pi --> pi*, n --> pi*, pi --> sigma* and n --> sigma* excited states of the titled compounds. The present results are compared with those of the previous theoretical studies by methods, such as EOM-CCSD(T), STEOM-CCSD, CASPT2 and TD-B3LYP, etc.

Structure, Bonding, and Lowest Energy Transitions in Unsymmetrical Squaraines: A Computational Study

Natural resonance theory (NRT) and natural bond orbital (NBO) analysis have been carried out on a simple symmetrical and an unsymmetrical substituted squaraine with a view of understanding the structure of the latter type of squaraines. It is found that there are some fundamental differences in the structure and bonding between these two types of squaraines particularly in the resonance weights and delocalization energies. These differences are expected to reflect in the low energy transitions and charge transfer in these squaraines. To investigate this, the nature of the lowest energy transitions occurring on excitation in unsymmetrical squaraines has been studied using high-level symmetry adapted cluster-configuration interaction method (SAC/SAC-CI) and compared with reported experimental observations. In general the agreement with the experimental data is very good. The transition dipole moment always lies on the pi-backbone and is quite large in magnitude. The ground state dipole moment in some cases does not change in the excited state upon excitation while in some other cases there is a large reduction/enhancement in the magnitude indicative of some charge rearrangement in this direction. Inclusion of the solvent using the IEFPCM model, a slightly better agreement with the experiment is found in some cases. Studies are carried out with a different basis set and it is found that the change in basis set has very little effect on the transition energies. In the case of weak side donor groups attached to the central ring the larger charge transfer to the central acceptor ring in general takes place from the O- atoms of the squarylium moiety while in the case of strong donors the charge transfer from the O- atoms to the central rings drop down. We have not observed any correlation between the charge transfer in the excited state to the central ring from the side donor groups and the lowest energy excitation in the molecules. Reduction of the HOMO-LUMO gap (an indication of increase of the diradicaloid character) always leads to a bathochromic shift.

Singly and doubly excited states of butadiene, acrolein, and glyoxal: Geometries and electronic spectra

Excited-state geometries and electronic spectra of butadiene, acrolein, and glyoxal have been investigated by the symmetry adapted cluster configuration interaction (SAC-CI) method in their s-trans conformation. Valence and Rydberg states below the ionization threshold have been precisely calculated with sufficiently flexible basis sets. Vertical and adiabatic excitation energies were well reproduced and the detailed assignments were given taking account of the second moments. The deviations of the vertical excitation energies from the experiment were less than 0.3 eV for all cases. The SAC-CI geometry optimization has been applied to some valence and Rydberg excited states of these molecules in the planar structure. The optimized ground- and excited-state geometries agree well with the available experimental values; deviations lie within 0.03 A and 0.7 degrees for the bond lengths and angles, respectively. The force acting on the nuclei caused by the excitations has been discussed in detail by calculating the SAC-CI electron density difference between the ground and excited states; the geometry relaxation was well interpreted with the electrostatic force theory. In Rydberg excitations, geometry changes were also noticed. Doubly excited states (so-called 2 (1)A(g) states) were investigated by the SAC-CI general-R method considering up to quadruple excitations. The characteristic geometrical changes and large energetic relaxations were predicted for these states.

Near-Infrared Absorption in Symmetric Squarylium and Croconate Dyes: A Comparative Study Using Symmetry-Adapted Cluster-Configuration Interaction Methods

Symmetric croconate (CR) and squarylium dyes (SQ) are well-known near-infrared (NIR) dyes and, in general, are considered to be donor-acceptor-donor type molecules. It is established in the literature that CR dyes absorb in a longer wavelength region than the corresponding SQ dyes. This has been attributed to the CR ring being a better acceptor than the SQ ring. Thus increasing the donor capacity should lead to a bathochromic shift in both SQ and CR. On the other hand, some experiments reported in the literature have revealed that increasing the conjugation in the donor part of the SQ molecule leads first to red shift, which upon a further increase of the conjugation changes to a blue shift. Hence, to understand the role of the central ring and the substitutions in the absorption of these dyes, we carried out high-level symmetry-adapted cluster-configuration interaction (SAC-CI) calculations of some substituted SQ and CR dyes and compare the absorption energy with the existing experimental data. We found that there is very good agreement. We also carried out SAC-CI calculations of some smaller model molecules, which contain the main oxyallyl substructure. We varied the geometry (angle) of the oxyallyl subgroup and the substitution in these model molecules to establish a correlation with the bathochromic shift. We found that the charge transfer is very small and does not play the key role in the red shift, but on the other hand, the perturbation of the HOMO-LUMO gap (HLG) from both the geometry and substitution seems to be responsible for this shift. We suggest as a design principle that increasing the donor capacity of the groups may not help in the red shift, but introducing groups which perturb the HLG and decrease it without changing the MO character should lead to a larger bathochromic shift.

Theoretical fine spectroscopy with symmetry-adapted-cluster configuration-interaction method: Outer- and inner-valence ionization spectra of furan, pyrrole, and thiophene

Theoretical fine spectroscopy has been performed for the valence ionization spectra of furan, pyrrole, and thiophene with the symmetry-adapted-cluster configuration-interaction general-R method. The present method described that the pi(1) state interacts with the pi(3) (-2)pi*, pi(2) (-2)pi*, and pi(2) (-1)pi(3) (-1)pi* shake-up states providing the split peaks and the outer-valence satellites, both of which are in agreement with the experiments. The intensity distributions were analyzed in detail for the inner-valence region. In particular, for furan, theoretical intensities were successfully compared with the intensity measured by the electron momentum spectroscopy. The interactions of the 3b(2) and 5a(1) states with the shake-up states were remarkable for furan and pyrrole, while the 4b(2) state of thiophene had relatively large intensity.

Theoretical Study of Singlet and Triplet Excitation Energies in Oligothiophenes

We have analyzed singlet and triplet excitation energies in oligothiophenes (up to five rings) using time-dependent density-functional theory (TD-DFT) with different exchange-correlation functionals and compared them with results from the approximate coupled-cluster singles and doubles model (CC2) and experimental data. The excitation energies have been calculated in geometries obtained by TD-DFT optimization of the lowest excited singlet state and in the ground-state geometries of the neutral and anionic systems. TD-DFT methods underestimate photoluminescence energies but the energy difference between singlet and triplet states shows trends with the chain-length similar to CC2. We find that the second triplet excited state is below the first singlet excited state for long oligomers in contrast with the previous assignment of Rentsch et al. (Phys.Chem. Chem. Phys. 1999, 1, 1707). Their photodetachment photoelectron spectroscopy measurements are better described by considering higher triplet excited states.

Role of the Oxyallyl Substructure in the Near Infrared (NIR) Absorption in Symmetrical Dye Derivatives: A Computational Study

It is well-known from experimental studies that the oxyallyl-substructure-based squarylium and croconium dyes absorb in the NIR region of the spectrum. Recently, another dye has been reported (J. Am. Chem. Soc. 2003, 125, 348) which contains the same basic chromophore, but the absorption is red-shifted by at least 300 nm compared to the former dyes and is observed near 1100 nm. To analyze the reasons behind the large red shift, in this work we have carried out symmetry-adapted cluster-configuration interaction (SAC-CI) studies on some of these NIR dyes which contain the oxyallyl substructure. From this study, contrary to the earlier reports, it is seen that the donor groups do not seem to play a major role in the red-shift of the absorption. On the other hand, on the basis of the results of the high-level calculations carried out here and using qualitative molecular orbital theory, it is observed that the orbital interactions play a key role in the red shift. Finally, design principles for the oxyallyl-substructure-based NIR dyes are suggested.

Spin-orbit coupling of DFT/MRCI wavefunctions: method, test calculations, and application to thiophene

During the past decade the one-center mean-field approximation has proven to be a very appropriate framework for the accurate description of spin-orbit effects at the correlated all-electron level. Here, a new efficient code, SPOCK, is introduced that calculates spin-orbit matrix elements in the one-center mean-field approximation for multireference CI wave functions. For the first time, the computation of spin-dependent interactions within a Kohn-Sham orbital based CI (DFT/MRCI) scheme1 is made possible. The latter approach is suitable for large scale systems with up to 100-200 valence electrons. Test calculations are performed on well-known diatomic molecules and the thiocarbonyl pyranthione. Spin-orbit matrix elements show good agreement with their Hartree-Fock orbital based counterparts but are obtained at considerably lower expense, thus demonstrating the power of the method. As an application singlet-triplet couplings in thiophene are investigated that are important for the photophysics and photochemistry. Spin-orbit matrix elements between all pi --> pi* excited states are found to be small. Considerably larger spin-orbit matrix elements are observed only for cases in which pi --> sigma* excited configurations are involved.