Intersystem Crossing as a Key Component of the Nonadiabatic Relaxation Dynamics of Bithiophene and Terthiophene

Ab initio investigation of excited state charge transfer pathways in differently capped bithiophene cages

The electronic excitations of conformationally constrained bithiophene cage systems as previously investigated by Lewis et al. (J. Am. Chem. Soc. 143, 18548 (2021)) are revisited, employing the correlated ab initio Scaled Opposite-Spin Algebraic Diagrammatic Construction Second Order electronic structure method. Quantitative descriptors are determined to assess the extent of charge transfer between the bithiophene moieties and the capping domains, represented by either phenyl or triazine groups. The investigation substantiates intrinsic differences in the photophysical behavior of these two structural variants and reveals the presence of lower-energy excited states characterized by noteworthy charge transfer contributions in the triazine cage system. The manifestation of this charge transfer character is discernible even at the Franck-Condon geometry, persisting throughout the relaxation of the excited state. By examining isolated monomer building blocks, we confirm the existence of analogous charge transfer contributions in their excitations. Employing this methodological approach facilitates the prospective identification of potential wall/cap chromophore pairs, wherein charge transfer pathways can be accessed within the energetically favorable regime.

New Gradient Correction Scheme for Electronically Nonadiabatic Dynamics Involving Multiple Spin States

It has been recommended that the best representation to use for trajectory surface hopping (TSH) calculations is the fully adiabatic basis in which the Hamiltonian is diagonal. Simulations of intersystem crossing processes with conventional TSH methods require an explicit computation of nonadiabatic coupling vectors (NACs) in the molecular-Coulomb-Hamiltonian (MCH) basis, also called the spin-orbit-free basis, in order to compute the gradient in the fully adiabatic basis (also called the diagonal representation). This explicit requirement destroys some of the advantages of the overlap-based algorithms and curvature-driven algorithms that can be used for the most efficient TSH calculations. Therefore, although these algorithms allow one to perform NAC-free simulations for internal conversion processes, one still requires NACs for intersystem crossing. Here, we show that how the NAC requirement is circumvented by a new computation scheme called the time-derivative-matrix scheme.

Optical and electronic properties enhancement via chalcogenides: promising materials for DSSC applications

CONTEXT

Comparatively, metal-free sensitizers featuring the chalcogen family receive less attention despite known electronic properties for metal-chalcogenide materials. This work reports an array of optoelectronic properties using quantum chemical methods. Observed red-shifted bands within the UV/Vis to NIR regions with absorption maxima > 500 nm were consistent with increasing chalcogenide size. There is a monotonic down-shift in the LUMO and ESOP energy consistent with O 2p, S 3p, Se 4p, to Te 5p atomic orbital energies. Excited-state lifetime and charge injection free energies follow the decreasing order of chalcogenide electronegativity. Adsorption energies of dyes on TiO₂ anatase (101) range between - 0.08 and - 0.77 eV. Based on evaluated properties, selenium- and tellurium-based materials show potential use in DSSCs and futuristic device applications. Therefore, this work motivates continued investigation of the chalcogenide sensitizers and their application.

METHODS

The geometry optimization was performed at B3LYP/6-31 + G(d,p) and B3LYP/LANL2DZ level of theory for lighter and heavier atoms, respectively, using Gaussian 09. The equilibrium geometries were confirmed by the absence of imaginary frequencies. Electronic spectra were obtained at CAM-B3LYP/6-31G + (d,p)/LANL2DZ level of theory. Adsorption energies for dyes on a 4 × 5 supercell TiO₂ anatase (101) were obtained using VASP. The dye-TiO₂ optimizations were employed using GGA and PBE with the PAW pseudo-potentials. The energy cutoff was set at 400 eV and convergence threshold for self-consistent iteration was set to 10^-4, and van der Waals were accounted using DFT-D3 model and on-site Coulomb repulsion potential set at 8.5 eV for Ti.

Taking up the quest for novel molecular solar thermal systems: Pros and cons of storing energy with cubane and cubadiene

Molecular solar thermal (MOST) systems are working their way as a possible technology to store solar light and release it when necessary. Such systems could, in principle, constitute a solution to the energy storage problem characteristic of solar cells and are conceived, at a first instance, as simple molecular photoswitches. Nevertheless, the optimization of their different required properties is presently limiting their technological scale up. From the chemical perspective, we need to design a novel MOST system based on unconventional photoswitches. Here, by applying multi-configurational quantum chemistry methods, we unravel the potentialities of ad hoc-designed molecular photoswitches, which aim to photoproduce cubane or cubadiene as high-energy isomers that can be thermally (or eventually catalytically) reverted to the initial structure, releasing their stored energy. Specifically, while cubane can be photoproduced via different paths depending on the reactant tricycle diene conformation, an undesired bicyclic by-product limits its application to MOST systems. An evolution of this starting design toward cubadiene formation is therefore proposed, avoiding conformational equilibria and by-products, considerably red shifting the absorption to reach the visible portion of the solar spectrum and maintaining an estimated storage density that is expected to overcome the current MOST reference system (norbornadiene/quadricyclane), although consistently increasing the photoisomerization energy barrier.

Accurate Calculation of Excited-State Absorption for Small-to-Medium-Sized Conjugated Oligomers: Multiconfigurational Treatment vs Quadratic Response TD-DFT

Excited-state absorption (ESA) spectra of π-conjugated compounds are frequently calculated by (quadratic response) time-dependent density functional theory, (QR) TD-DFT, often giving a reasonable representation of the experimental results despite the (known) incomplete electronic description. To investigate whether this is inherent to the method, we calculate here the ESA spectra of small-to-medium-sized oligophenylenevinylenes (nPV) and oligothiophenes (nT) using QR TD-DFT as well as CASPT2 based on CASSCF geometries. CASPT2 gives indeed a reliable, theoretically correct description of the ESA features for all compounds; the computational effort can be reduced without significant loss of accuracy using TD-DFT geometries. QR TD-DFT, based on BHandHLYP and CAM-/B3LYP functionals, fails on short nTs but provides a reasonable description for spectral positions of nPVs and long nTs. The failure on short nTs is, however, only partly due to the incomplete configuration description but, in particular, related to an improper MO description, resulting in an asymmetric energy spacing of the occupied vs unoccupied MOs in the DFT scheme. Longer nTs, on the other side, adapt approximately the MO scheme for alternant hydrocarbons just like in nPVs, while contributions by two triplet excitations combined to a singlet (which inhibits an accurate treatment of polyenes with standard TD-DFT) do not play a relevant role in the current case. For such "well-behaved" systems, a reasonable representation of ESA spectra is found at the QR TD-DFT level due to the rather small energy shifts when including higher-order excitations.

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

Direct coherent switching with decay of mixing for intersystem crossing dynamics of thioformaldehyde: The effect of decoherence

We evaluate the effect of electronic decoherence on intersystem crossing in the photodynamics of thioformaldehyde. First, we show that the state-averaged complete-active-space self-consistent field electronic structure calculations with a properly chosen active space of 12 active electrons in 10 active orbitals can predict the potential energy surfaces and the singlet-triplet spin-orbit couplings quite well for CH₂S, and we use this method for direct dynamics by coherent switching with decay of mixing (CSDM). We obtain similar dynamical results with CSDM or by adding energy-based decoherence to trajectory surface hopping, with the population of triplet states tending to a small steady-state value over 500 fs. Without decoherence, the state populations calculated by the conventional trajectory surface hopping method or the semiclassical Ehrenfest method gradually increase. This difference shows that decoherence changes the nature of the results not just quantitatively but qualitatively.

The role of CO2 detachment in fungal bioluminescence: thermally vs. excited state induced pathways

Different fungi lineages are known to emit light on Earth, mainly in tropical climates. Although the preparation of bioluminescent cell-free extracts allowed one to characterize the enzymatic requirements, the molecular mechanism underlying luminescence is still largely unknown and is based on the experimental putative assumption that a high-energy intermediate should be formed by reaction with O2 and formation of an endoperoxide. Here, we aim at determining, through state-of-the-art multiconfigurational quantum chemistry, the full mechanistic landscape leading from the endoperoxide to the emitting species, envisaging different possible pathways and proposing their viability. Especially, thermal CO2 detachment followed by excited-state peroxide opening (thermal-chemiluminescence) can compete with a parallel pathway, i.e., first excited-state endoperoxide opening, followed by CO2 detachment on the same excited-state (excited state-chemiluminescence). Clear differences in the energy supplies, as well as the possibility to directly populate the emitting species from the intersection seam between ground and excited states, land credence to a kinetically efficient thermal-chemiluminescent pathway, establishing for the first time a detailed description of fungal bioluminescence.

Simulation of time-resolved x-ray absorption spectroscopy of ultrafast dynamics in particle-hole-excited 4-(2-thienyl)-2,1,3-benzothiadiazole

To date, alternating co-polymers based on electron-rich and electron-poor units are the most attractive materials to control functionality of organic semiconductor layers in which ultrafast excited-state processes play a key role. We present a computational study of the photoinduced excited-state dynamics of the 4-(2-thienyl)-2,1,3-benzothiadiazole (BT-1T) molecule, which is a common building block in the backbone of π-conjugated polymers used for organic electronics. In contrast to homo-polymer materials, such as oligothiophene, BT-1T has two non-identical units, namely, thiophene and benzothiadiazole, making it attractive for intramolecular charge transfer studies. To gain a thorough understanding of the coupling of excited-state dynamics with nuclear motion, we consider a scenario based on femtosecond time-resolved x-ray absorption spectroscopy using an x-ray free-electron laser in combination with a synchronized ultraviolet femtosecond laser. Using Tully's fewest switches surface hopping approach in combination with excited-state calculations at the level of configuration interaction singles, we calculate the gas-phase x-ray absorption spectrum at the carbon and nitrogen K edges as a function of time after excitation to the lowest electronically excited state. The results of our time-resolved calculations exhibit the charge transfer driven by non-Born-Oppenheimer physics from the benzothiadiazole to thiophene units during relaxation to the ground state. Furthermore, our ab initio molecular dynamics simulations indicate that the excited-state relaxation processes involve bond elongation in the benzothiadiazole unit as well as thiophene ring puckering at a time scale of 100 fs. We show that these dynamical trends can be identified from the time-dependent x-ray absorption spectrum.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Adiabatic-Molecular Dynamics Generalized Vertical Hessian Approach: A Mixed Quantum Classical Method To Compute Electronic Spectra of Flexible Molecules in the Condensed Phase

We present a general mixed quantum classical method that couples classical molecular dynamics (MD) and vibronic models to compute the shape of electronic spectra of flexible molecules in the condensed phase without, in principle, any phenomenological broadening. It is based on a partition of the nuclear motions of the solute + solvent system in "soft" and "stiff" vibrational modes and an adiabatic hypothesis that assumes that stiff modes are much faster than soft ones. In this framework, the spectrum is rigorously expressed as a conformational integral of quantum vibronic spectra along the stiff coordinates only. Soft modes enter at the classical level through the conformational distribution that is sampled with classical MD runs. In each configuration, reduced-dimensionality quadratic Hamiltonians are built in the space of the stiff coordinates only, thanks to a generalization of the Vertical Hessian harmonic model and an iterative application of projectors in internal coordinates to remove soft modes. Quantum vibronic spectra, specific for each sampled configuration of the soft coordinates, are then computed at the desired temperature with efficient time-dependent techniques, and the global spectrum simply arises from their average. For consistency of the whole procedure, classical MD runs are performed with quantum-mechanically derived force fields, parameterized at the same level of theory selected for generating the quadratic Hamiltonians along the stiff coordinates. Application to N-methyl-6-oxyquinolinium betaine in water, dithiophene in ethanol, and cyanidine in water is presented to show the performance of the method.

Luminescent lanthanide complexes with a pyridine-bis(carboxamide)-bithiophene sensitizer showing wavelength-dependent singlet oxygen generation

A new pyridine-bis(carboxamide)-based ligand with a bithiophene pendant, 2Tcbx, was synthesized.