Intersystem crossing in a dibenzofuran-based room temperature phosphorescent luminophore investigated by non-adiabatic dynamics

The use of phosphorescent luminophores is highly beneficial in diverse high-technological and biological applications. Yet, because of the formally forbidden character of intersystem crossing, the use of heavy metals or atoms is usually necessary to achieve high quantum yields. This choice imposes serious constraints in terms of high device cost and inherent toxicity. In this contribution we resort to density functional based surface hopping non-adiabatic dynamics of a potential organic luminophore intended for room-temperature applications. We confirm that intersystem crossing is operative in a ps time-scale without requiring the activation of large-scale movements, thus confirming the suitability of the El Sayed-based strategy for the rational design of fully organic phosphorescent emitters.

Exploring the ultrafast and isomer-dependent photodissociation of iodothiophenes via site-selective ionization

C-I bond extension and fission following ultraviolet (UV, 262 nm) photoexcitation of 2- and 3-iodothiophene is studied using ultrafast time-resolved extreme ultraviolet (XUV) ionization in conjunction with velocity map ion imaging. The photoexcited molecules and eventual I atom products are probed by site-selective ionization at the I 4d edge using intense XUV pulses, which induce multiple charges initially localized to the iodine atom. At C-I separations below the critical distance for charge transfer (CT), charge can redistribute around the molecule leading to Coulomb explosion and charged fragments with high kinetic energy. At greater C-I separations, beyond the critical distance, CT is no longer possible and the measured kinetic energies of the charged iodine atoms report on the neutral dissociation process. The time and momentum resolved measurements allow determination of the timescales and the respective product momentum and kinetic energy distributions for both isomers, which are interpreted in terms of rival 'direct' and 'indirect' dissociation pathways. The measurements are compared with a classical over the barrier model, which reveals that the onset of the indirect dissociation process is delayed by ∼1 ps relative to the direct process. The kinetics of the two processes show no discernible difference between the two parent isomers, but the branching between the direct and indirect dissociation channels and the respective product momentum distributions show isomer dependencies. The greater relative yield of indirect dissociation products from 262 nm photolysis of 3-iodothiophene (cf. 2-iodothiophene) is attributed to the different partial cross-sections for (ring-centred) π∗ ← π and (C-I bond localized) σ∗ ← (n/π) excitation in the respective parent isomers.

The structure of liquid thiophene from total neutron scattering

The structure of pure liquid thiophene is revealed by using a combination of total neutron scattering experiments with isotopic substitution and molecular simulations via the next generation empirical potential refinement software, Dissolve. In the liquid, thiophene presents three principle local structural motifs within the first solvation shell, in plane and out of the plane of the thiophene ring. Firstly, above/below the ring plane thiophenes present a single H towards the π cloud, due to a combination of electrostatic and dispersion interactions. Secondly, around the ring plane, perpendicular thiophene molecules find 5 preferred sites driven by bifurcated C-H⋯S interactions, showing that hydrogen-sulfur bonding prevails over the charge asymmetry created by the heteroatom. Finally, parallel thiophenes sit above and below the ring, excluded from directly above the ring center and above the sulfur.

Probing C-I bond fission in the UV photochemistry of 2-iodothiophene with core-to-valence transient absorption spectroscopy

The UV photochemistry of small heteroaromatic molecules serves as a testbed for understanding fundamental photo-induced chemical transformations in moderately complex compounds, including isomerization, ring-opening, and molecular dissociation. Here, a combined experimental-theoretical study of 268 nm UV light-induced dynamics in 2-iodothiophene (C4H3IS) is performed. The dynamics are experimentally monitored with a femtosecond extreme ultraviolet (XUV) probe that measures iodine N-edge 4d core-to-valence transitions. Experiments are complemented by density functional theory calculations of both the pump-pulse induced valence excitations and the XUV probe-induced core-to-valence transitions. Possible intramolecular relaxation dynamics are investigated by ab initio molecular dynamics simulations. Gradual absorption changes up to ∼0.5 to 1 ps after excitation are observed for both the parent molecular species and emerging iodine fragments, with the latter appearing with a characteristic rise time of 160 ± 30 fs. Comparison of spectral intensities and energies with the calculations identifies an iodine dissociation pathway initiated by a predominant π → π* excitation. In contrast, initial excitation to a nearby n⟂ → σ* state appears unlikely based on a significantly smaller oscillator strength and the absence of any corresponding XUV absorption signatures. Excitation to the π → π* state is followed by contraction of the C-I bond, enabling a nonadiabatic transition to a dissociative π→σC-I* state. For the subsequent fragmentation, a relatively narrow bond-length region along the C-I stretch coordinate between 230 and 280 pm is identified, where the transition between the parent molecule and the thienyl radical + iodine atom products becomes prominent in the XUV spectrum due to rapid localization of two singly occupied molecular orbitals on the two fragments.

Energy-dependent timescales in the dissociation of diiodothiophene dication

Photodissociation molecular dynamics of gas-phase 2,5-diiodothiophene molecules was studied in an electron-energy-resolved electron-multi-ion coincidence experiment performed at the FinEstBeAMS beamline of MAX IV synchrotron. Following the photoionization of the iodine 4d subshell and the Auger decay, the dissociation landscape of the molecular dication was investigated as a function of the Auger electron energy. Concentrating on an major dissociation pathway, C₄H₂I₂S²⁺ → C₄H₂S⁺ + I⁺ + I, and accessing the timescales of the process via ion momentum correlation analysis, it was revealed how this three-body process changes depending on the available internal energy. Using a generalized secondary dissociation model, the process was shown to evolve from secondary dissociation regime towards concerted dissociation as the available energy increased, with the secondary dissociation time constant changing from 1.5 ps to 129 fs. The experimental results were compared with simulations using a stochastic charge-hopping molecular mechanics model. It represented the observed trend and also gave a fair quantitative agreement with the experiment.

Nonradiative relaxation mechanisms of the elusive silole molecule

Silole derivatives have been extensively employed for developing organic optoelectronics, but few studies focused on the photophysical properties of the silole molecule. In this work, we investigate these properties by computing the absorption spectra and performing nonadiabatic molecular dynamics of silole employing the algebraic diagrammatic construction [ADC(2)] and extended multi-state XMS-CASPT2 ab initio electronic structure methods. For vertical excitations and excited state optimizations, the equation of motion coupled-cluster singles and doubles (EOM-CCSD) was also used. The nuclear ensemble and the fewest-switches surface hopping molecular dynamics methods coupled with the first two high-level electronic structure methods were applied to probe the relaxation mechanisms of silole. We could reproduce the experimental first absorption maximum value and found an ultrafast relaxation process occurring exclusively through ring-puckering distortions without breaking ring bonds or hydrogen elimination. Minimum energy conical intersection optimizations were carried out and potential energy curves, including triplet states, were calculated to further elucidate the relaxation process of silole.

Surface hopping dynamics reveal ultrafast triplet generation promoted by S1-T2-T1 spin-vibronic coupling in 2-mercaptobenzothiazole

We investigate the T₁ formation upon populating the optically "bright" S₂ in 2-mercaptobenzothiazole to interpret the underlying relaxation pathways associated with the experimental decay constants reported by D. Koyama and A. J. Orr-Ewing, Phys. Chem. Chem. Phys., 2016, 18, 26224-26235. Energetics, electronic populations and geometries of various stationary points of low-lying electronic states are computed using the semi-classical ab initio surface hopping dynamics simulations. Estimated decay constants of S₂-S₁ internal conversion (IC) and S₁-T₂ intersystem crossing (ISC) are in excellent agreement with the experiment. The observed ultrafast ISC is analyzed based on the S₁-T₂-T₁ spin-vibronic coupling mechanism. In contrast to the previous assignment of 6 ps to the T₂-T₁ IC, our findings enable us to attribute this decay constant to the combined events of T₂-T₁ IC followed by relaxation of vibrationally hot T₁.

Enabling a Unified Description of Both Internal Conversion and Intersystem Crossing in Formaldehyde: A Global Coupled Quasi-Diabatic Hamiltonian for Its S0, S1, and T1 States

In our recent work, a diabatic Hamiltonian that couples the S₀ and S₁ states of formaldehyde was constructed using a robust fitting-and-diabatizing procedure with artificial neural networks, which is capable of representing adiabatic energies, energy gradients, and derivative couplings over a wide range of geometries including seams of conical intersection. In this work, based on the diabatization of S₀ and S₁, the spin-orbit couplings between singlet states (S₀, S₁) and triplet state T₁ are also determined in the same diabatic representation. The diabatized spin-orbit couplings are then fit with a symmetrized neural-network functional form. The ab initio spin-orbit couplings are well reproduced in large configuration space. Together with the neural-network-based potential energy surface for T₁, the full quasi-diabatic Hamiltonian for the S₀, S₁, and T₁ states is completed, enabling a unified description of both internal conversion and intersystem crossing in formaldehyde. The vibrational levels on the three adiabatic states are found to be in good agreement with known experimental band origins.

Tailoring CIPSI Expansions for QMC Calculations of Electronic Excitations: The Case Study of Thiophene

The perturbatively selected configuration interaction scheme (CIPSI) is particularly effective in constructing determinantal expansions for quantum Monte Carlo (QMC) simulations with Jastrow–Slater wave functions: fast and smooth convergence of ground-state properties and balanced descriptions of ground and excited states of different symmetries have been reported. In particular, accurate excitation energies have been obtained by the pivotal requirement of using CIPSI expansions with similar second-order perturbation corrections for each state, that is, a similar estimated distance to the full configuration interaction limit. Here, we elaborate on the CIPSI selection criterion for excited states of the same symmetry as the ground state, generating expansions from a common orbital set. Using these expansions in QMC as determinantal components of Jastrow–Slater wave functions, we compute the lowest, bright excited state of thiophene, which is challenging due to its significant multireference character. The resulting vertical excitation energies are within 0.05 eV of the best theoretical estimates, already with expansions of only a few thousand determinants. Furthermore, we relax the ground- and excited-state structures following the corresponding root in variational Monte Carlo and obtain bond lengths that are accurate to better than 0.01 Å. Therefore, while the full treatment at the CIPSI level of this system is quite demanding, in QMC, we can compute high-quality excitation energies and excited-state structural parameters building on affordable CIPSI expansions with relatively few, well-chosen determinants.

Exploring the effects of quantum decoherence on the excited-state dynamics of molecular systems

Fewest-switches surface hopping (FSSH) is employed in order to investigate the nonadiabatic excited-state dynamics of thiophene and related compounds and hence to establish a connection between the electronic system, the critical points in configuration space and the deactivation dynamics. The potential-energy surfaces of the studied molecules were calculated with complete active space self-consistent field and time-dependent density-functional theory. They are analyzed thoroughly to locate and optimize minimum-energy conical intersections, which are essential to the dynamics of the system. The influence of decoherence on the dynamics is examined by employing different decoherence schemes. We find that irrespective of the employed decoherence algorithm, the population dynamics of thiophene give results which are sound with the expectations grounded on the analysis of the potential-energy surface. A more detailed look at single trajectories as well as on the excited-state lifetimes, however, reveals a substantial dependence on how decoherence is accounted for. In order to connect these findings, we describe how ensemble averaging cures some of the overcoherence problems of uncorrected FSSH. Eventually, we identify carbon–sulfur bond cleavage as a common feature accompanying electronic transitions between different states in the simulations of all thiophene-related compounds studied in this work, which is of interest due to their relevance in organic photovoltaics.

Nonadiabatic dynamics simulation of photoinduced ring-opening reaction of 2(5H)-thiophenone with internal conversion and intersystem crossing

In the present work, the quantum trajectory mean-field approach, which is able to overcome the overcoherence problem, was generalized to simulate internal conversion and intersystem crossing processes simultaneously. The photoinduced ring-opening and subsequent rearrangement reactions of isolated 2(5H)-thiophenone were studied based on geometry optimizations on critical structures and nonadiabatic dynamics simulations using this method. Upon 267 nm irradiation, the molecule is initially populated in the 1ππ* state. After a sudden rupture of one C-S bond within 100 fs in this state, the lowest two singlet excited states and the lowest two triplet excited states become quasi-degenerated, and then the intersystem crossing processes between singlet and triplet states accompanied by rearrangement reactions can be observed several times. Compared with our previous nonadiabatic simulations in the absence of intersystem crossing (ChemPhotoChem, 2019, 3, 897-906), some new nonadiabatic relaxation pathways involving triplet states and different ring-opening products were identified. The present work provides new mechanistic insights into the photoinduced ring-opening of thio-substituted heterocyclic molecules and reveals the importance of nonadiabatic dynamics simulation that is able to deal with multiple electronic states with different spin multiplicities.

Formation of excited triplet states in naphthalene diimide and perylene diimide derivatives: A detailed theoretical study

Mechanistic details of the excited triplet state formation upon photoexcitation to the low-lying singlet manifold in naphthalene diimide and perylene diimide derivatives are explored theoretically. Static and dynamic aspects of two singlets (S₁ and S₂) and six triplets (T₁-T₆) of these molecules are investigated. Suitable vibronic Hamiltonians are constructed to investigate the internal conversion dynamics in both the singlet and triplet manifolds. Computed singlet-triplet energetics, spin-orbit coupling matrix elements, and intersystem crossing rates strongly suggest an efficient intersystem crossing process involving higher triplet states (T₆, T₅, and T₄). Separate full dimensional quantum wavepacket simulations of singlet and triplet manifolds in the approximate linear vibronic model by assuming initial Franck-Condon conditions are carried out to unravel the internal conversion decay dynamics in the respective manifolds. The obtained diabatic electronic populations and nuclear densities are analyzed to illustrate the triplet generation pathways involving higher triplet states in these molecules.

Experimental and theoretical studies on thymine photodimerization mediated by oxidatively generated DNA lesions and epigenetic intermediates

Combined spectroscopic and computational studies reveal that, in spite of their structural similarities, 5-formyluracil and 5-formylcytosine photosensitize cyclobutane thymine dimers through two different types of mechanisms.

Ground-state dissociation pathways for the molecular cation of 2-chlorothiophene: A time-of-flight mass spectrometry and computational study

RATIONALE

Halogenated thiophenes are an important class of compounds mostly used in the synthesis of various materials, showing unusual electronic and optical properties. The Thiophene Ring Fragmentation (TRF) process is widely used in synthetic chemistry. In this study, the fragmentation pattern of the molecular cation of halogenated thiophene, namely, 2-chlorothiophene, was monitored to establish its dissociation mechanism.

METHODS

The molecular cation of 2-chlorothiophene was prepared using multiphoton excitation using a laser at 235 nm. Various product ions upon fragmentation of the molecular ion were mass analyzed using time-of-flight mass spectrometry. Laser power dependence studies were also conducted for various product ions to arrive at the dissociation mechanism. Theoretical calculations were carried out to estimate the reaction enthalpies for various reactions and compared with the experimental data available in the literature.

RESULTS

The most abundant product ion was observed as the HCS⁺ radical cation followed by the C₃ H₃ ⁺ ion and the H₂ CCCCS⁺ radical cation. Other product ions such as SCCl⁺ , ClHCCS⁺ radical cations were also observed to a lesser extent in the fragmentation pattern of the parent molecular ion. Various dissociation channels were identified and supported with ab initio calculation. It has been inferred that the TRF process is usually initiated by the H/Cl atom transfer process. The appearance energies of the various fragment ions were also estimated theoretically and compared with literature values.

CONCLUSIONS

In conclusion, the fragmentation pattern of the molecular cation of 2-chlorothiophene was studied and the formation mechanisms of various product ions have been assigned. The appearance energies of the various fragment ions were also calculated. Finally, it is inferred that a TRF process is initiated by the H/Cl atom migration and subsequent ring opening either by C-C or C-S bond cleavage leading to the various isomers and their subsequent fragmentation. The ionization energies were accurately predicted for various species using ab initio calculation.

A combined multi-reference pump-probe simulation method with application to XUV signatures of ultrafast methyl iodide photodissociation

UV pump-XUV/X-ray probe measurements have been successfully applied in the study of photo-induced chemical reactions. Although rich element-specific electronic structure information is accessible within XUV/X-ray (inner-shell) absorption spectra, it can be difficult to interpret the chemistry directly from the spectrum without supporting theoretical simulations. A multireference method to completely simulate UV pump-XUV/X-ray probe measurement has been developed and applied to study the methyl iodide photodissociation process. Multireference, fewest-switches surface hopping (FSSH) trajectories were used to explore the coupled electronic and ionic dynamics upon photoexcitation of methyl iodide. Interpretation of previous measurements is provided by associated multireference, restricted active space, inner-shell spectral simulations. This combination of multireference FSSH trajectories and XUV spectra provides an interpretation of transient features appearing in previous measurements within the first 100 fs after photoexcitation and validates the significant branching ratio in the final excited-state population. This methodology should prove useful for interpretation of the increasing number of inner-shell probe studies of molecular excited states or for directing new experiments toward interesting regions of the potential energy landscape.

Observation of Ultrafast Intersystem Crossing in Thymine by Extreme Ultraviolet Time-Resolved Photoelectron Spectroscopy

We studied the photoinduced ultrafast relaxation dynamics of the nucleobase thymine using gas-phase time-resolved photoelectron spectroscopy. By employing extreme ultraviolet pulses from high harmonic generation for photoionization, we substantially extend our spectral observation window with respect to previous studies. This enables us to follow relaxation of the excited state population all the way to low-lying electronic states including the ground state. In thymine, we observe relaxation from the optically bright ¹ππ* state of thymine to a dark ¹nπ* state within 80 ± 30 fs. The ¹nπ* state relaxes further within 3.5 ± 0.3 ps to a low-lying electronic state. By comparison with quantum chemical simulations, we can unambiguously assign its spectroscopic signature to the ³ππ* state. Hence, our study draws a comprehensive picture of the relaxation mechanism of thymine including ultrafast intersystem crossing to the triplet manifold.

Large Excited-State Conformational Displacements Expedite Triplet Formation in a Small Conjugated Oligomer

Intersystem crossing in conjugated organic molecules is most conveniently viewed from pure electronic perspectives; yet, vibrational displacements may often drive these transitions. We investigate an alkyl-substituted thienylene-vinylene dimer (dTV) displaying efficient triplet formation. Steady-state electronic and Raman spectra display large Stokes shifts (∼4000 cm^-1) involving high-frequency skeletal symmetric stretching modes (∼900-1600 cm^-1) in addition to large displacements of low-frequency torsional motions (∼300-340 cm^-1). Transient absorption spectroscopy reveals the emergence of distorted singlet (S₁) and triplet signatures following initial vibrational relaxation dynamics that dominate spectral dynamics on time scales > 100 ps, with the latter persisting on time scales up to ca. 7 μs. Potential energy surfaces calculated along the dominant displaced out-of-plane torsional mode reveal shallow energy barriers for entering the triplet manifold from S₁. We propose that dTV is a good model system for understanding vibrational contributions to intersystem crossing events in related polymer systems.

Conical intersections and intersystem crossings explain product formation in photochemical reactions of aryl azides

Photochemistry of substituted aryl azides is governed by surface crossings. Internal conversion and intersystem crossing govern photodecomposition of 3-methoxyphenyl azide and 4-methoxyphenyl azide.

Nonadiabatic dynamics: The SHARC approach

We review the Surface Hopping including ARbitrary Couplings (SHARC) approach for excited-state nonadiabatic dynamics simulations. As a generalization of the popular surface hopping method, SHARC allows simulating the full-dimensional dynamics of molecules including any type of coupling terms beyond nonadiabatic couplings. Examples of these arbitrary couplings include spin-orbit couplings or dipole moment-laser field couplings, such that SHARC can describe ultrafast internal conversion, intersystem crossing, and radiative processes. The key step of the SHARC approach consists of a diagonalization of the Hamiltonian including these couplings, such that the nuclear dynamics is carried out on potential energy surfaces including the effects of the couplings-this is critical in any applications considering, for example, transition metal complexes or strong laser fields. We also give an overview over the new SHARC2.0 dynamics software package, released under the GNU General Public License, which implements the SHARC approach and several analysis tools. The review closes with a brief survey of applications where SHARC was employed to study the nonadiabatic dynamics of a wide range of molecular systems. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and KineticsSoftware > Simulation MethodsSoftware > Quantum Chemistry.

Theoretical Investigations on Mechanisms and Pathways of C₂H₅O₂ with BrO Reaction in the Atmosphere

In this work, feasible mechanisms and pathways of the C₂H₅O₂ + BrO reaction in the atmosphere were investigated using quantum chemistry methods, i.e., QCISD(T)/6-311++G(2df,2p)//B3LYP/6-311++G(2df,2p) levels of theory. Our result indicates that the title reaction occurs on both the singlet and triplet potential energy surfaces (PESs). Kinetically, singlet C₂H₅O₃Br and C₂H₅O₂BrO were dominant products under the atmospheric conditions below 300 K. CH₃CHO₂ + HOBr, CH₃CHO + HOBrO, and CH₃CHO + HBrO₂ are feasible to a certain extent thermodynamically. Because of high energy barriers, all products formed on the triplet PES are negligible. Moreover, time-dependent density functional theory (TDDFT) calculation implies that C₂H₅O₃Br and C₂H₅O₂BrO will photolyze under the sunlight.