FCclasses3: Vibrationally-resolved spectra simulated at the edge of the harmonic approximation

Vibronic Coupling and Multiple Electronic States Effect in ABS and ECD Spectra: Three [7]Helicene Derivatives

Vibronic coupling and multiple electronic states effect play a pivotal role in the molecular spectroscopy of large systems. Herein, we present a detailed theoretical study on the absorption (ABS) and electronic circular dichroism (ECD) spectra of three [7]helicene derivatives in chloroform, with a particular emphasis on the significance of vibronic coupling and the multiple electronic states effect in spectral simulations. The vertical gradient (VG) and vertical Hessian (VH) models, incorporating the Franck-Condon (FC) effect and Herzberg-Teller (HT) contribution, are considered in the vibronic calculations. The results indicate that the simulated vibrationally resolved spectra obtained by the VG model combined with the FC effect are more reliable, showing advantages in the rationality of spectral shapes, the accuracy of the relative heights of absorption peaks, and the correctness of positive and negative signals. Reliable predictions of the three [7]helicene derivatives allowed us to further explore the importance of the multiple electronic states effect in the vibrationally resolved spectra, demonstrating that the high-energy electronic excited states, particularly the fifth (S5) and sixth (S6) excited states, are essential for accurately capturing the fine structures observed in the experimental spectra. Our study predicts reliable theoretical reference spectra for the family of [7]helicene derivatives and provides a fundamental understanding of the vibronic coupling of chiral organic molecules with multiple electronic states.

Localized and Excimer Triplet Electronic States of Naphthalene Dimers: A Computational Study

We perform DFT calculations with different hybrid (ωB97X-D and M05-2X) and double hybrid (B2PLYP-D3 and ωB2PLYP) functionals to characterize the lowest energy triplet excited states of naphthalene monomer and dimers in different stacking arrangements and to simulate their absorption spectra. We show that both excimer and localized triplet minima exist. In the former, the spin density is delocalized over the two monomers, adopting a face-to-face arrangement with a short inter-molecular distance. In the latter, the spin density is localized on a single naphthalene molecule, and different minima or pseudo-minima are possible, the most stable one corresponding to a slipped parallel arrangement. According to B2PLYP-D3 calculations, excimer minima are the most stable, in line with the indications of ADC(2) studies. However, the relative stability of the minima is reverted when including thermal and vibrational effects. Excimer minima exhibit a very intense absorption spectrum, peaking above 500 nm. The computed absorption spectra of localized minima significantly depend on the stacking geometry and do not coincide with that of isolated naphthalene. Hybrid functionals provide very accurate vibronic absorption spectra for naphthalene monomer, both in the singlet and in the triplet state, but underestimate the stability of the excimer triplet.

Charge Transfer-Induced Weakening of Vibronic Coupling for Single Terrylene Molecules Adsorbed onto Hexagonal Boron Nitride

Fluorescence spectra of single terrylene molecules adsorbed on hexagonal boron nitride flakes were recorded at cryogenic temperatures. The pure electronic transitions of terrylene molecules are spread over a broad energy scale from 570 to 610 nm. Surprisingly, peaks in the vibrationally resolved fluorescence spectrum show intensity variations of ≤20-fold between molecules. We find an extreme case in which the Debye-Waller-Franck-Condon factor of the zero-phonon line exceeds 0.8. The vibronic intensity correlates with both the spectral position of the electronic transition and the frequency of the longitudinal stretch mode, which varies between 243 and 257 cm-1. Using density functional theory calculations, we show that these observations can be explained by terrylene chemisorption on charge-donating defect sites. The electronic states of molecules at such chemisorption sites would be very attractive for the efficient emission of single photons with narrow lines and for the generation of indistinguishable photons.

Vibronic Effects Analysis of the Substituent Effect on the Spectral Properties of the EMI and CPL of Three [7]Helicene Derivatives

The substituent effect has a significant influence on the optical properties of spectral shape, width, and wavelength, and the intensities of the maximum peaks of emission (EMI) and circularly polarized luminescence (CPL). In this work, we conducted a systematic theoretical study to investigate how substituents alter the optical response in the EMI and CPL spectra of three [7]helicene derivatives at the vibronic level. To incorporate the vibronic effect, a state-of-the-art time-dependent (TD) method was used to achieve the fully converged spectra. In the meantime, a time-independent (TI) approach also provided a way to show the progression of the spectra, serving as a complementary strategy and creating reliable documentation for the experiment. The experimental spectra of EMI and CPL are nicely reproduced, which validates the reliability of the Adiabatic Hessian (AH) model in simulating experimental data. This allowed us to analyze in detail the effect of substituents, particularly on the optical responses. The introduction of cyano and methoxy groups is highlighted, as they altered the transition dipole moments and led to a 1000-fold increase in the intensity of EMI and CPL. Moreover, substituents can also rationally alter the spectral shape of EMI and CPL by affecting the responsible normal modes. The unique CN stretching and the MeO rotation in the substituted [7]helicene are highlighted as key factors contributing to the different behaviors of EMI and CPL. This sheds light on the design and synthesis of helicenes that can serve as ideal full-color EMI and CPL emitters.

Structures, Chiroptical Properties, and Unexpectedly Facile Helical Inversion of Highly Elongated Anthracene-Fused Expanded Helicenes

Helical fused anthracenes were elongated by fusing additional aromatic units at both ends to yield novel expanded helicenes. Compounds [5]HA2N and [7]HA consisting of 19 and 21 benzene rings, respectively, were synthesized by fourfold cycloisomerization of the corresponding terminal alkyne precursors. The helical structures were confirmed by X-ray crystallographic analysis, where the aromatic frameworks stacked effectively with the helical turn numbers exceeding two. The enantiomers of the two compounds were resolved by chiral HPLC. Whereas [5]HA2N readily underwent enantiomerization at room temperature at the barrier to enantiomerization of 91 kJ mol-1, the barrier was enhanced to 99 kJ mol-1 for the long analog [7]HA. The enantiomers of [7]HA exhibited strong responses in the circular dichroism (CD) and circularly polarized luminescence (CPL) spectra, as scaled by dissymmetry factors |gabs|=0.034 and |glum|=0.012. Theoretical calculations by the r2SCAN-3c method suggested stepwise mechanisms for the enantiomerization via helical inversion with acceptable barrier heights. The unexpectedly flexible nature of the aromatic frameworks of [5]HA2N and [7]HA was discussed on the basis of the proposed mechanism.

The Golden Age of Thermally Activated Delayed Fluorescence Materials: Design and Exploitation

Since the seminal report by Adachi and co-workers in 2012, there has been a veritable explosion of interest in the design of thermally activated delayed fluorescence (TADF) compounds, particularly as emitters for organic light-emitting diodes (OLEDs). With rapid advancements and innovation in materials design, the efficiencies of TADF OLEDs for each of the primary color points as well as for white devices now rival those of state-of-the-art phosphorescent emitters. Beyond electroluminescent devices, TADF compounds have also found increasing utility and applications in numerous related fields, from photocatalysis, to sensing, to imaging and beyond. Following from our previous review in 2017 ( Adv. Mater. 2017, 1605444), we here comprehensively document subsequent advances made in TADF materials design and their uses from 2017-2022. Correlations highlighted between structure and properties as well as detailed comparisons and analyses should assist future TADF materials development. The necessarily broadened breadth and scope of this review attests to the bustling activity in this field. We note that the rapidly expanding and accelerating research activity in TADF material development is indicative of a field that has reached adolescence, with an exciting maturity still yet to come.

Calibration of several first excited state properties for organic molecules through systematic comparison of TDDFT with experimental spectra

Time-dependent density functional theory (TDDFT) is a powerful computational tool for investigating excitation properties in organic electronics, and it holds significant potential for high-throughput virtual screening (HTVS) in this field. While most benchmarks focus on excitation energies, less attention has been paid to evaluating the accuracy of computed oscillator strengths and exciton reorganization energies against experimental data. In this work, we provide a systematic approach to evaluate in parallel the accuracy of these three quantities on the basis of a suitable fitting of the experimental absorption spectra of 71 molecules in solution. After considering 18 computational methodologies, the results from the M06-2X/def2-TZVP/PCM method demonstrate the strongest correlation with experimental data across the desired properties. For HTVS, the M06-2X/6-31G(d)/PCM method appears to be a particularly convenient choice among all methodologies due to its balance of computational efficiency and accuracy. Our results provide an additional benchmark needed before employing TDDFT methods for the discovery and design of organic electronic molecules.

Influence of symmetry breaking on the absorption spectrum of crystal violet: from isolated cations at 5 K to room temperature solutions

We report the resolution of a long-standing puzzle in molecular spectroscopy: the origin of the shoulder in the room temperature solution absorption spectrum of crystal violet (CV) - an archetypal cationic triphenylmethane dye. This was achieved by comparing experimental and theoretical results for CV in solution at room temperature and as an isolated cation in gas-phase at 5 K. The two lowest energy electronically excited states involved in the visible region absorption are degenerate and coupled via a Jahn-Teller (JT) mechanism involving phenyl torsions, making CV particularly sensitive to environmental perturbations. The shoulder is absent in the low-temperature isolated cation spectrum, and vibronic simulations based on time dependent density functional theory (TD-DFT) indicate negligible JT effects under these conditions. Combining vibronic simulations with molecular dynamics, demonstrates that in water and toluene solution at room temperature the shoulder arises mainly from an intermolecular, Jahn-Teller-like symmetry-breaking effect induced by the fluctuating electrostatic potential of the disordered solvent environment, rather than from molecular distortions.

Toward Efficient Modeling of Nonradiative Decay in Extended INVEST: Overcoming Computational Challenges in Quantum Dynamics Simulations

In recent years, an increasing number of fully organic molecules capable of thermally activated delayed fluorescence (TADF) have been reported, often with very small or even inverted singlet-triplet (INVEST) energy gaps. These molecules typically exhibit complex photophysics due to the close energy levels of multiple singlet and triplet states, which create various transition pathways toward emission. A predictive model for the rates of these transitions is thus essential for assessing the suitability of new materials for light-emitting devices. Quantum Dynamics (QD) calculations are ideal for this purpose, as they include quantum effects, without the limitations of first-order perturbative approaches, also allowing taking into account more than two electronic states at once. However, the huge computational demands of QD methodologies, especially for large molecules, currently limit their use as a standard tool. To address this problem, we here employ a strategy that allows us to include almost the whole set of the vibrational coordinates by selecting the key elements of the Hilbert space that significantly impact dynamics, thereby hugely reducing the computational burden. Application of this protocol to two relatively large INVEST molecules reveals that internal conversion in these systems is very fast, making indirect emissive pathways a possible channel for the population of the S1 state. More importantly, this study demonstrates that the dynamics can be accurately described even with a significantly reduced vibrational space, thus allowing quantum dynamics calculations that yield accurate transition rates in a few minutes of computational time.

Rational Control of Maximum EMI/CPL Intensity and Wavelength of Bora[6]helicene via Polarity and Vibronic Effects

Solvent polarity control as an efficient methodology to regulate the chiroptical properties, including spectral shape, width, intensity, wavelength, etc., has emerged as a novel frontier in optical materials design. However, the underling relationship connecting polarity to the optical property remains unclear. Herein, using state-of-the-art computations and the FC|VG model, the solvent effect on the chiroptical properties of bora[6]helicene was accurately and systematically computed to shed light on this issue. It is found that the vibronic coupling is crucial in explaining the spectral shape, width, and relative intensity of different peaks. Moreover, the intensity and position of the emission (EMI) and circularly polarized luminescence (CPL) are closely related to the polarity of the solvent. Intriguingly, we got a series of good linear relationships between polarity and EMI|CPL (|r| ≥ 0.95). Thus, this parameter can be used as a potential descriptor to estimate the intensity and position of EMI|CPL, leading to new strategies for designing fully colored fluorescent materials.

Study of the Energy Crossing Between Excited States Affected by the Electronegativity of Substituents for Three 4-Azido-1,8-naphthalimide Derivatives

Rapid detection of H2S is crucial for human physiological health and natural ecosystems. In this study, the fluorescent sensing mechanisms of three 4-azido-1,8-naphthalimide-based fluorescent probes to monitor H2S were theoretically investigated by density functional theory and time-dependent density functional theory. The potential energy curve of the charge transfer (CT) state has a crossover with that of the locally excited (LE) state proved by the constructed linear interpolating internal coordinate pathway. Thus, the transform takes place from the LE state to the CT state causing the fluorescence quenching of the probes from a nonradiative transition process of the CT state. The distance between the Franck-Condon point and the minimal energy conical intersection becomes larger with the increase of the electronegativity of substituents on the 1,8-naphthalimide fluorophore. In addition, the degree of charge separation is closely related to the energy difference between the CT and the LE states which are also essentially affected by the electronegativity of the substituents. Since the electronegativity of the substituents has proved important for the probes, our work lays a certain theoretical foundation for the design and synthesis of more sensitive 4-azido-1,8-naphthalimide-based fluorescent probes.

Ultra-Narrowband Circularly Polarized Luminescence from Multiple 1,4-Azaborine-Embedded Helical Nanographenes

In this manuscript we present a strategy to achieve ultranarrowband circularly polarized luminescence (CPL) from multiple 1,4-azaborine-embedded helical nanographenes. The impact of number and position of boron and nitrogen atoms in the rigid core of the molecule on optical properties─including absorption and emission maxima, photoluminescence quantum yield, Stokes shift, excited singlet-triplet energy gap and full width at half-maximum (fwhm) for CPL and fluorescence─was investigated. The molecules reported here exhibits ultranarrowband fluorescence (fwhm 16-17.5 nm in toluene) and CPL (fwhm 18-19 nm in toluene). To the best of our knowledge, this is among the narrowest CPL for any organic molecule reported to date. Quantum chemical calculations, including computed CPL spectra involving vibronic contributions, provide valuable insights for future molecular design aimed at achieving narrowband CPL.

AtmoSpec-A Tool to Calculate Photoabsorption Cross-Sections for Atmospheric Volatile Organic Compounds

Characterizing the photolysis processes undergone by transient volatile organic compounds (VOCs) in the troposphere requires the knowledge of their photoabsorption cross-section-quantities often challenging to determine experimentally, particularly due to the reactivity of these molecules. We present a computational tool coined AtmoSpec, which can predict a quantitative photoabsorption cross-section for volatile organic compounds by using computational photochemistry. The user enters the molecule of interest as a SMILES code and, after selecting a level of theory for the electronic structure (and waiting for the calculations to take place), is presented with a photoabsorption cross-section for the low-energy conformers and an estimate of the photolysis rate coefficient for different standardized actinic fluxes. More specifically, AtmoSpec is an automated workflow for the nuclear ensemble approach, an efficient technique to approximate the absolute intensities and excitation wavelengths of a photoabsorption cross-section for a molecule in the gas phase of interest in atmospheric chemistry and astrochemistry. This work provides background information on the nuclear ensemble approach, a guided example of a typical AtmoSpec calculation, details about the architecture of the code, and the current limitations and future developments of this tool.

Toward the improvement of vibronic spectra and non-radiative rate constants using the vertical Hessian method

For the computation of vibrationally resolved electronic spectra, various approaches can be employed. Adiabatic approaches simulate vibronic transitions using harmonic potentials of the initial and final states, while vertical approaches extrapolate the final state potential from the gradients and Hessian at the Franck-Condon point, avoiding a full exploration of the potential energy surface of the final state. Our implementation of the vertical Hessian (VH) method has been validated with a benchmark set of four small molecules, each presenting unique challenges, such as complex topologies, problematic low-frequency vibrations, or significant geometrical changes upon electronic excitation. We assess the quality of both adiabatic and vertical approaches for simulating vibronic transitions. For two types of donor-acceptor compounds with promising thermally activated delayed fluorescence properties, our computations confirm that the vertical approaches outperform the adiabatic ones. The VH method significantly reduces computational costs and yields meaningful emission spectra, where adiabatic models fail. More importantly, we pioneer the use of the VH method for the computation of rate constants for non-radiative processes, such as intersystem crossing and reverse intersystem crossing along a relaxed interpolated pathway of a donor-acceptor compound. This study highlights the potential of the VH method to advance computational vibronic spectroscopy by providing meaningful simulations of intricate decay pathway mechanisms in complex molecular systems.

The Excited State Dynamics of a Mutagenic Guanosine Etheno Adduct Investigated by Femtosecond Fluorescence Spectroscopy and Quantum Mechanical Calculations

Femtosecond fluorescence upconversion experiments were combined with CASPT2 and time dependent DFT calculations to characterize the excited state dynamics of the mutagenic etheno adduct 1,N2-etheno-2'-deoxyguanosine (ϵdG). This endogenously formed lesion is attracting great interest because of its ubiquity in human tissues and its highly mutagenic properties. The ϵdG fluorescence is strongly modified with respect to that of the canonical nucleoside dG, notably by an about 6-fold increase in fluorescence lifetime and quantum yield at neutral pH. In addition, femtosecond fluorescence upconversion experiments reveal the presence of two emission bands with maxima at 335 nm for the shorter-lived and 425 nm for the longer-lived. Quantum mechanical calculations rationalize these findings and provide absorption and fluorescence spectral shapes similar to the experimental ones. Two different bright minima are located on the potential energy surface of the lowest energy singlet excited state. One planar minimum, slightly more stable, is associated with the emission at 335 nm, whereas the other one, with a bent etheno ring, is associated with the red-shifted emission.

Rational Design of Color-Pure Blue Organic Emitters by Poly-Heteroaromatic Omni-Delocalization

Current research on organic light emitters which utilize multiple resonance-induced thermally activated delayed fluorescence (MR-TADF) materials is gaining significant interest because of the materials' ability to efficiently generate color-pure blue emission. However, the underlying reasons for high color purity remain unclear. It is shown here that these emitters share a common electronic basis, which is deduced from resonance structure considerations following Clar's rule, and which is termed as "poly-heteroaromatic omni-delocalization" (PHOD). The simple and clear design rules derived from the PHOD concept allow extending the known chemical space by new structural motifs. Based on PHOD, a set of novel high-efficiency color-pure emitters with brilliant deep-blue hue is specifically designed.

Theoretical Study on Vibrationally Resolved Electronic Spectra of Chiral Nanographenes

Nanographenes are of increasing importance owing to their potential applications in the photoelectronic field. Meanwhile, recent studies have primarily focused on the pure electronic spectra of nanographenes, which have been found to be inadequate for describing the experimental spectra that contain vibronic progressions. In this study, we focused on the vibronic effect on the electronic transition of a range of chiral nanographenes, especially in the low-energy regions with distinct vibronic progressions, using theoretical calculations. All the calculations were performed at the PBE0-D3(BJ)/def2-TZVP level of theory, adopting both time-dependent and time-independent approaches with Franck-Condon approximation. The resulting calculated curves exhibited good alignment with the experimental data. Notably, for the nanographenes incorporating helicene units, owing to the increasing π-extension, the major vibronic modes in the vibrationally resolved spectra differed significantly from those of the primitive helicenes. This investigation suggests that calculations that account for the vibronic effect could have better reproducibility compared with calculations based solely on pure electronic transitions. We anticipate that this study could pave the way for further investigations into optical and chiroptical properties, with a deeper understanding of the vibronic effect, thereby providing theoretical explanations with higher precision on more sophisticated nanographenes.

Explicit Modelling of Spectral Bandshapes by a Mixed Quantum-Classical Approach: Solvent Order and Temperature Effects in the Optical Spectra of Distyrylbenzene

The absorption and emission spectral shapes of a flexible organic probe, the distyrylbenzene (DSB) dye, are simulated accounting for the effect of different environments of increasing complexity, ranging from a homogeneous, low-molecular- weight solvent, to a long-chain alkane, and, eventually, a channel-forming organic matrix. Each embedding is treated explicitly, adopting a mixed quantum-classical approach, the Adiabatic Molecular Dynamics - generalized vertical Hessian (Ad-MD|gVH) model, which allows a direct simulation of the environment-induced constraining effects on the vibronic spectral shapes. In such a theoretical framework, the stiff modes of the dye are described at a quantum level within the harmonic approximation, including Duschinsky mixing effects, while flexible degrees of freedom of the solute (e. g. torsions) and those of the solvent are treated classically by means of molecular dynamics sampling. Such a setup is shown to reproduce the distinct effects exerted by the different environments in varied thermodynamic conditions. Besides allowing for a first-principles rationale on the supramolecular mechanism leading to the experimental spectral features, this result represents the first successful application of the Ad-MD|gVH method to complex embeddings and supports its potential application to other heterogeneous environments, such as for instance, pigment-protein complexes or organic dyes adsorbed into metal-organic frameworks.

Time-dependent surface-enhanced Raman scattering: A theoretical approach

A new procedure for computing the time-dependent Raman scattering of molecules in the proximity of plasmonic nanoparticles (NPs) is proposed, drawing inspiration from the pioneering Lee and Heller's theory. This strategy is based on a preliminary simulation of the molecular vibronic wavefunction in the presence of a plasmonic nanostructure and an incident light pulse. Subsequently, the Raman signal is evaluated through an inverse Fourier Transform of the coefficients' dynamics. Employing a multiscale approach, the system is treated by coupling the quantum mechanical description of the molecule with the polarizable continuum model for the NP. This method offers a unique advantage by providing insights into the time evolution of the plasmon-enhanced Raman signal, tracking the dynamics of the incident electric field. It not only provides for the total Raman signal at the process's conclusion but also gives transient information. Importantly, the flexibility of this approach allows for the simulation of various incident electric field profiles, enabling a closer alignment with experimental setups. This adaptability ensures that the method is relevant and applicable to diverse real-world scenarios.

Full Electron Delocalization across the Cluster in 1,12-bisBMes2-p-carborane Radical Anion

Conjugation between three-dimensional (3D) carboranes and the attached substituents is commonly believed to be very weak. In this paper, we report that reducing 1,12-bis(BMes2)-p-carborane (B2pCab) with one electron gives a radical anion with a centrosymmetric semiquinoidal structure. This radical anion shows extensive electron delocalization between the two boron centers over the p-carborane bridge due to the overlap of carborane lowest unoccupied molecular orbital (LUMO) and the BMes2 LUMO. Unlike dianions of other C2B10H12 carboranes, which rearrange to a nido-form, two-electron reduction of B2pCab leads to a rearrangement into a basket-shaped intermediate.

Newly Designed Quasi-intrinsic Photosensitizers for Fluorescence Image-Guided Two-Photon Photodynamic Therapy with Type I/II Photoreactions

In this work, a set of quasi-intrinsic photosensitizers are theoretically proposed based on the 2-amino-8-(1'-β-d-2'-deoxyribofuranosyl)-imidazo[1,2-α]-1,3,5-triazin-4(8H)-one (P), which could pair with the 6-amino-5-nitro-3-(1'-β-d-2'-deoxyribofuranosyl)-2(1H)-pyridone (Z) and keep the essential structural characters of nucleic acid. It is revealed that the ring expansion and electron-donating/electron-withdrawing substitution bring enhanced two-photon absorption and bright photoluminescence of these monomers, thereby facilitating the selective excitation and tumor localization through fluorescence imaging. However, instead of undergoing radiative transition (S1 → S0), the base pairing induced fluorescence quenching and rapid intersystem crossing (S1 → Tn) are observed and characterized by the reduced singlet-triplet energy gaps and large spin-orbit coupling values. To ensure the phototherapeutic properties of the considered base pairs in long-lived T1 state, we examined the vertical electron affinity as well as vertical ionization potential for production of superoxide anions via Type I photoreaction, and their required T1 energy (0.98 eV) to generate singlet oxygen 1O2 via Type II mechanism.

Structure and Chiroptical Properties of Anthra[1,2-a]anthracene-1-yl Dimers as New Biaryls

Dimers of anthra[1,2-a]anthracene-1-yl units and its mesityl derivative were synthesized by Ni(0)-mediated coupling of the corresponding chloro derivatives as new biaryls. The X-ray analysis and DFT calculations revealed that two polycyclic aromatic units with nonplanar deformations took a twisted conformation about the single bond as a chiral axis. Enantiomers of the nonsubstituted compound were resolved by chiral HPLC, and the enantiopure samples showed intense Cotton effects at 321 nm in the circular dichroism (CD) spectra and emission bands at 449 nm in the circularly polarized luminescence (CPL) spectra with dissymmetry factor of |glum| 3.6×10-3. The absolute stereochemistry of this biaryl was determined by the theoretical calculation of CD spectrum by the time-dependent DFT method. The barrier to enantiomerization was determined to be 108 kJ mol-1 at 298 K. The dynamic process proceeded via a stepwise mechanism involving the helical inversion of each aromatic unit and the rotation about the biaryl axis as analyzed by the DFT calculations.

Deconvolution of the X-ray absorption spectrum of trans-1,3-butadiene with resonant Auger spectroscopy

High-resolution carbon K-edge X-ray photoelectron, X-ray absorption, non-resonant and resonant Auger spectra are presented of gas phase trans-1,3-butadiene alongside a detailed theoretical analysis utilising nuclear ensemble approaches and vibronic models to simulate the spectroscopic observables. The resonant Auger spectra recorded across the first pre-edge band reveal a complex evolution of different electronic states which remain relatively well-localised on the edge or central carbon sites. The results demonstrate the sensitivity of the resonant Auger observables to the weighted contributions from multiple electronic states. The gradually evolving spectral features can be accurately and feasibly simulated within nuclear ensemble methods and interpreted with the population analysis.

Linear Vibronic Coupling Approach for Surface-Enhanced Raman Scattering: Quantifying the Charge-Transfer Enhancement Mechanism

The outstanding amplification observed in surface-enhanced Raman scattering (SERS) is due to several enhancement mechanisms, and standing out among them are the plasmonic (PL) and charge-transfer (CT) mechanisms. The theoretical estimation of the enhancement factors of the CT mechanism is challenging because the excited-state coupling between bright plasmons and dark CT states must be properly introduced into the model to obtain reliable intensities. In this work, we aim at simulating electrochemical SERS spectra, considering models of pyridine on silver clusters subjected to an external electric field E⃗ that represents the effect of an electrode potential Vel. The method adopts quantum dynamical propagations of nuclear wavepackets on the coupled PL and CT states described with linear vibronic coupling models parametrized for each E⃗ through a fragment-based maximum-overlap diabatization. By presenting results at different values of E⃗, we show that indeed there is a relation between the population transferred to the CT states and the total scattered intensity. The tuning and detuning processes of the CT states with the bright PLs as a function of the electric field are in good agreement with those observed in experiments. Finally, our estimations for the CT enhancement factors predict values in the order of 105 to 106, meaning that when the CT and PL states are both in resonance with the excitation wavelength, the CT and PL enhancements are comparable, and vibrational bands whose intensity is amplified by different mechanisms can be observed together, in agreement with what was measured by typical experiments on silver electrodes.

Charge Transfer-Promoted Excited State of a Heavy-Atom-Free Photosensitizer for Efficient Application of Mitochondria-Targeted Fluorescence Imaging and Hypoxia Photodynamic Therapy

Conventional photosensitizers (PSs) used in photodynamic therapy (PDT) have shown preliminary success; however, they are often associated with several limitations including potential dark toxicity in healthy tissues, limited efficacy under acidic and hypoxic conditions, suboptimal fluorescence imaging capabilities, and nonspecific targeting during treatment. In response to these challenges, we developed a heavy-atom-free PS, denoted as Cz-SB, by incorporating ethyl carbazole into a thiophene-fused BODIPY core. A comprehensive investigation into the photophysical properties of Cz-SB was conducted through a synergistic approach involving experimental and computational investigations. The enhancement of intersystem crossing (kISC) and fluorescence emission (kfl) rate constants was achieved through a donor-acceptor pair-mediated charge transfer mechanism. Consequently, Cz-SB demonstrated remarkable efficiency in generating reactive oxygen species (ROS) under acidic and low-oxygen conditions, making it particularly effective for hypoxic cancer PDT. Furthermore, Cz-SB exhibited good biocompatibility, fluorescence imaging capabilities, and a high degree of localization within the mitochondria of living cells. We posit that Cz-SB holds substantial prospects as a versatile PS with innovative molecular design, representing a potential "one-for-all" solution in the realm of cancer phototheranostics.

One- and two-photon absorption spectra of organoboron complexes: vibronic and environmental effects

We synthesized a series of four parent aza-β-ketoiminate organoboron complexes and performed spectroscopic studies using both experimental and computational techniques. We studied how benzannulation influences the vibronic structure of the UV/Vis absorption bands with a focus on the bright lowest-energy π → π* electronic excitation. Theoretical simulations, accounting for inhomogeneous broadening effects using different embedding schemes, allowed gaining in-depth insights into the observed differences in band shapes induced by structural modifications. We observed huge variations in the distributions of vibronic transitions depending on the position of benzannulation. By and large, the harmonic approximation combined with the adiabatic hessian model delivers qualitatively correct band shapes for the one-photon absorption spectra, except in one case. We also assessed the importance of non-Condon effects (accounted for by the linear term in Herzberg-Teller expansion of the dipole moment) for S0 → S1 band shapes. It turned out that non-Condon contributions have no effect on the band shape in one-photon absorption spectra. In contrast, these effects significantly change the Franck-Condon band shapes of the two-photon absorption spectra. For one of the studied organoboron complexes we also performed a preliminary exploration of mechanical anharmonicity, resulting in an increase of the intensity of the 0-0 transition, which improves the agreement with the experimental data compared to the harmonic model.

Blue-Shifted and Broadened Fluorescence Enhancement by Visible and Mode-Selective Infrared Double Excitations

Mode-selective vibrational excitations to modify the electronic states of fluorescein dianion in methanol solutions are carried out with a femtosecond visible pulse synchronized with a tunable high-power, narrow-band picosecond infrared (IR) pulse. In this work, simultaneous intensity enhancement, peak blueshift, and line width broadening of fluorescence are observed in the visible/IR double resonance experiments. Comprehensive investigations on the modulation mechanism with scanning the vibrational excitation frequencies, tuning the time delay between the two excitation pulses, theoretical calculations, and nonlinear and linear spectroscopic measurements suggest that the fluorescence intensity enhancement is caused by the increase of the Franck-Condon factor induced by the vibrational excitations at the electronic ground state. Various enhancement effects are observed as vibrations initially excited by the IR photons relax to populate the vibrational modes of lower frequencies. The peak blueshift and line width broadening are caused by both increasing the Franck-Condon factors among different subensembles because of IR pre-excitation and the long-lived processes induced by the initial IR excitation. The results demonstrate that the fluorescence from the visible/IR double resonance experiments is not a simple sum frequency effect, and vibrational relaxations can produce profound effects modifying luminescence.

Efficient near-infrared emission benefits from slowing down the internal conversion process

Organic deep-red (DR) and near-infrared (NIR) emitters with high photoluminescence quantum yield (PLQY) are rare due to the strong non-radiative (knr) decay. Here, we report two DR/NIR emitters with high PLQY, TPANZPyPI and TPANZ3PI. Interestingly, the TPANZPyPI film exhibits 46.5% PLQY at 699 nm. Theoretical calculations indicate that TPANZPyPI can achieve this high PLQY in the near-infrared emission region due to its small S1 to S0 internal conversion (IC) rate. Meanwhile, research has found that, compared to TPANZ3PI, TPANZPyPI with a more rigid structure can effectively suppress the T2 to T1 IC process, which is conducive to higher exciton utilization efficiency (EUE). TPANZPyPI's non-doped OLED shows NIR emission with 4.6% @ 684 nm maximum external quantum efficiency (EQEmax). Its doped OLEDs radiate DR with an EQEmax of 6.9% @ 666 nm. These EQEs are among the highest values for hybridized local charge transfer state materials emitting more than 640 nm. This work demonstrates for the first time, based on a combination of theory and experiment, that increasing the molecular rigidity can inhibit the excited state IC process in addition to the S1 to S0 IC, realizing efficient electroluminescence.

Significance of Vibronic Coupling that Shapes Circularly Polarized Luminescence of Double Helicenes

The circularly polarized luminescence (CPL) spectra of S- and X-shaped double helicenes exhibit distinct vibrational structures and overall shape variations. In this study, we conducted an in-depth investigation into the vibronic effects influencing the CPL spectra of two double helicenes, namely DPC and DNH. Employing state-of-the-art computations utilizing the FC-HT1|VH model at the CAM-B3LYP/def2-TZVP level, we unveiled the paramount impact of Franck-Condon (FC), Herzberg-Teller (HT), and Duschinsky effects on their chiroptical responses. Our research underscores the pivotal role of structural deformations associated with the S1-to-S0 electronic transition in molding CPL spectra and wavelength-dependent dissymmetry (g) factor values, as well as the significance of HT effects in shaping and enhancing CPL responses. This extensive investigation not only advances our comprehension of the vibronic characteristics in configurationally distinct double helicenes but also offers valuable insights for the design of chiral molecules featuring controllable or finely-tunable CPL responses.

Tetracene Diacid Aggregates for Directing Energy Flow toward Triplet Pairs

A comprehensive investigation of the solution-phase photophysics of tetracene bis-carboxylic acid [5,12-tetracenepropiolic acid (Tc-DA)] and its related methyl ester [5,12-tetracenepropynoate (Tc-DE)], a non-hydrogen-bonding counterpart, reveals the role of the carboxylic acid moiety in driving molecular aggregation and concomitant excited-state behavior. Low-concentration solutions of Tc-DA exhibit similar properties to the popular 5,12-bis((triisopropylsilyl)ethynl)tetracene, but as the concentration increases, evidence for aggregates that form excimers and a new mixed-state species with charge-transfer (CT) and correlated triplet pair (TT) character is revealed by transient absorption and fluorescence experiments. Aggregates of Tc-DA evolve further with concentration toward an additional phase that is dominated by the mixed CT/TT state which is the only state present in Tc-DE aggregates and can be modulated with the solvent polarity. Computational modeling finds that cofacial arrangement of Tc-DA and Tc-DE subunits is the most stable aggregate structure and this agrees with results from 1H NMR spectroscopy. The calculated spectra of these cofacial dimers replicate the observed broadening in ground-state absorption as well as accurately predict the formation of a near-UV transition associated with a CT between molecular subunits that is unique to the specific aggregate structure. Taken together, the results suggest that the hydrogen bonding between Tc-DA molecules and the associated disruption of hydrogen bonding with solvent produce a regime of dimer-like behavior, absent in Tc-DE, that favors excimers rather than CT/TT mixed states. The control of aggregate size and structure using distinct functional groups, solute concentration, and solvent in tetracene promises new avenues for its use in light-harvesting schemes.

A CASSCF/MRCI trajectory surface hopping simulation of the photochemical dynamics and the gas phase ultrafast electron diffraction patterns of cyclobutanone

We present the simulation of the photochemical dynamics of cyclobutanone induced by the excitation of the 3 s Rydberg state. For this purpose, we apply the complete active space self-consistent field method together with the spin-orbit multireference configuration interaction singles treatment, combined with the trajectory surface hopping for the inclusion of nonadiabatic effects. The simulations were performed in the spin-adiabatic representation, including nine electronic states derived from three singlet and two triplet spin-diabatic states. Our simulations reproduce the two previously observed primary dissociation channels: the C2 pathway yielding C2H4 + CH2CO and the C3 pathway producing c-C3H6 + CO. In addition, two secondary products, CH2 + CO from the C2 pathway and C3H6 from the C3 pathway, both of them previously reported, are also observed in our simulation. We determine the ratio of the C3:C2 products to be about 2.8. Our findings show that most of the trajectories reach their electronic ground state within 200 fs, with dissociation events finished after 300 fs. We also identify the minimum energy conical intersections that are responsible for the relaxation and provide an analysis of the photochemical reaction mechanism based on multidimensional scaling. Furthermore, we demonstrate a minimal impact of triplet states on the photodissociation mechanism within the observed timescale. In order to provide a direct link to experiments, we simulate the gas phase ultrafast electron diffraction patterns and connect their features to the underlying structural dynamics.

Threading a Linear Molecule Through a Macrocycle Thanks to Boron: Optical Properties of the Threaded Species and Synthesis of a Rotaxane

Two BODIPYs and two boron β-diketonates were threaded through a macrocycle bearing a 2,2'-biphenol unit, showing thus the ability of boron to act as a gathering atom. The new threaded species were characterized by 1D and 2D NMR spectroscopy as well as by X-ray crystallography for one of them and their properties rationalized with quantum chemistry to unravel the vibronic contributions. The BODIPYs exhibited interesting fluorescence features with quantum yields up to 91 % and enhanced photostability compared to their non-threaded homologues. A rotaxane was synthesized using this threading strategy after stoppering and removing the boron with potassium hydroxide.

Diradical-Based Strategy in Designing Narrowband Thermally Activated Delayed Fluorescence Molecules with Tunable Emission Wavelengths

In the design of thermally activated delayed fluorescence (TADF) materials, narrow-band emission is of particular importance for the development of organic light-emitting diodes (OLEDs). In this work, we proposed a new strategy for designing TADF molecules utilizing degenerate nonbonding (NB) orbitals of diradical parent molecules, and these designed molecules are termed NB-TADF molecules. Based on this strategy, a series of NB-TADF molecules is finely designed and systematically studied by theoretical calculations. Taking advantage of the nonbonding properties, these NB-TADF molecules exhibit desirable narrowband emissions and high quantum yields. More importantly, the emission bands can be easily tuned from blue to near-infrared by changing the conjugate length of the parent group in the NB-TADF molecules. We hope that this new strategy can open a new door for the design of novel TADF materials.

The Electronic Structures of Azaphenanthrenes and Their Dimers

Insertion of a nitrogen atom modifies the electronic structures and photochemistry of polycyclic aromatic hydrocarbons by introducing nπ* states into the molecules. To better understand the electronic structures of isolated polycyclic aromatic nitrogen-containing hydrocarbons (PANHs) and their dimers as well as the influence of the position of the nitrogen atom in the molecule, we investigate three different azaphenanthrenes, benzo[f]quinoline, benzo[h]quinoline, and phenanthridine, in a joint experimental and computational study. Experimentally, resonance-enhanced multiphoton ionization (REMPI) spectroscopy is applied to characterize the excited electronic states. The REMPI spectra of the azaphenanthrene monomers have a rather similar appearance, with origins between 3.645 and 3.670 eV for the 1ππ* ← S0 transition. In contrast to the phenanthrene parent, 2ππ* ← S0 is broad and unstructured even at the band origin. The experiments are accompanied by density functional theory computation, and vibrationally resolved spectra are simulated using a time-independent approach. The differences between phenanthrene and the azaphenanthrenes are assigned to perturbations due to the low-lying 1(nπ*) state, which accelerates nonradiative deactivation. For the dimers, it is found that two π-stacked isomers with two electronic transitions each contribute to the electronic spectrum, leading to overlapping bands that are difficult to assign.

Accidental triplet harvesting in donor-acceptor dyads with low spin-orbit coupling

We present an accidental mechanism for efficient intersystem crossing (ISC) between singlet and triplet states with low spin-orbit coupling (SOC) in molecules having donor-acceptor (D-A) moieties separated by a Sigma bond. Our study shows that SOC between the lowest singlet excited state and the higher-lying triplet states, together with nuclear motion-driven coupling of this triplet state with lower-lying triplet states during the free rotation about a Sigma bond, is one of the possible ways to achieve the experimentally observed ISC rate for a class of D-A type photoredox catalysts. This mechanism is found to be the dominant contributor to the ISC process with the corresponding rate reaching a maximum at a dihedral angle in the range of 72°-78° between the D-A moieties of 10-(naphthalen-1-yl)-3,7-diphenyl-10H-phenoxazine and other molecules included in the study. We have further demonstrated that the same mechanism is operative in a specific spirobis[anthracene]dione molecule, where the D and A moieties are interlocked near to the optimal dihedral angle, indicating the plausible effectiveness of the proposed mechanism. The present finding is expected to have implications in strategies for the synthesis of new generations of triplet-harvesting organic molecules.

Intramolecular rate-constant calculations based on the correlation function using temperature dependent quantum Green's functions

A theoretical method for calculating rate constants for internal conversion (IC), intersystem crossing (ISC) and radiative (R) electronic transitions is presented. The employed method uses temperature-dependent quantum Green's functions, which give the opportunity to consider almost any nth-order polynomial perturbation operator and the influence of external electromagnetic fields on the rate constants. The rate constants of the IC, ISC and R processes are calculated for two important indocyanine molecules namely indocyanine green (ICG) and heptamethine cyanine (IR808) at the Franck-Condon level using the temperature-dependent quantum Green's function approach. Calculations at the time-dependent density functional theory level with the MN15 functional show that ICG and IR808 have only one triplet state below the S1 state. The main deactivation channel of the S1 state is the IC process with a large (kIC(S1 → S0)) rate constant of ∼109-1011 s-1. The estimated quantum yield of fluorescence (φfl) is ∼0.001-0.24 for the two studied molecules, which agrees rather well with experimental values. Thus, the present approach enables calculations of the three kinds of rate constants and the quantum yield of fluorescence using the same computational methodology.

Modeling the Fluorescence Quantum Yields of Aromatic Compounds: Benchmarking the Machinery to Compute Intersystem Crossing Rates

The from-first-principles calculation of fluorescence quantum yields (FQYs) and lifetimes of organic dyes remains very challenging. In this article, we extensively test the machinery to calculate FQYs. Specifically, we perform an extensive analysis on the parameters influencing the intersystem crossing (ISC), internal conversion (IC), and fluorescence rate constants calculations. The impact of (i) the electronic structure (chosen exchange-correlation functional and spin-orbit Hamiltonian), (ii) the vibronic parameters (coordinate system, broadening function, and dipole expansion), and (iii) the excited-state kinetic models are systematically assessed for a series of seven rigid aromatic molecules. Our studies provide more insights into the choice of parameters and the expected accuracy for the computational protocols aiming to deliver FQY values. Some challenges are highlighted, such as, on the one hand, the difficulty to benchmark against the experimental nonradiative rate constants, for which the separation between the IC and ISC contributions is often not provided in the literature and, on the other hand, the need to go beyond the harmonic approximation for the calculation of the IC rates.

Prediction of fluorescence quantum yields using the extended thawed Gaussian approximation

Spontaneous emission and internal conversion rates are calculated within harmonic approximations and compared to the results obtained within the semi-classical extended thawed Gaussian approximation (ETGA). This is the first application of the ETGA in the calculation of internal conversion and emission rates for real molecular systems, namely, formaldehyde, fluorobenzene, azulene, and a dicyano-squaraine dye. The viability of the models as black-box tools for prediction of spontaneous emission and internal conversion rates is assessed. All calculations were done using a consistent protocol in order to investigate how different methods perform without previous experimental knowledge using density functional theory (DFT) and time-dependent DFT (TD-DFT) with B3LYP, PBE0, ωB97XD, and CAM-B3LYP functionals. Contrasting the results with experimental data shows that there are further improvements required before theoretical predictions of emission and internal conversion rates can be used as reliable indicators for the photo-luminescence properties of molecules. We find that the ETGA performs rather similar to the vertical harmonical model. Including anharmonicities in the calculation of internal conversion rates has a moderate effect on the quantitative results in the studied systems. The emission rates are fairly stable with respect to computational parameters, but the internal conversion rate reveals itself to be highly dependent on the choice of the spectral line shape function, particularly the width of the Lorentzian function, associated with homogeneous broadening.

First-Principles Calculations of Excited-State Decay Rate Constants in Organic Fluorophores

In this Perspective, we discuss recent advances made to evaluate from first-principles the excited-state decay rate constants of organic fluorophores, focusing on the so-called static strategy. In this strategy, one essentially takes advantage of Fermi's golden rule (FGR) to evaluate rate constants at key points of the potential energy surfaces, a procedure that can be refined in a variety of ways. In this way, the radiative rate constant can be straightforwardly obtained by integrating the fluorescence line shape, itself determined from vibronic calculations. Likewise, FGR allows for a consistent calculation of the internal conversion (related to the non-adiabatic couplings) in the weak-coupling regime and intersystem crossing rates, therefore giving access to estimates of the emission yields when no complex photophysical phenomenon is at play. Beyond outlining the underlying theories, we summarize here the results of benchmarks performed for various types of rates, highlighting that both the quality of the vibronic calculations and the accuracy of the relative energies are crucial to reaching semiquantitative estimates. Finally, we illustrate the successes and challenges in determining the fluorescence quantum yields using a series of organic fluorophores.

Gold-Catalyzed 1,2-Aryl Shift and Double Alkyne Benzannulation

The tandem intramolecular hydroarylation of alkynes accompanied by a 1,2-aryl shift is described. Harnessing the unique electron-rich character of 1,4-dihydropyrrolo[3,2-b]pyrrole scaffold, we demonstrate that the hydroarylation of alkynes proceeds at the already occupied positions 2 and 5 leading to a 1,2-aryl shift. Remarkably, the reaction proceeds only in the presence of cationic gold catalyst, and it leads to heretofore unknown π-expanded, centrosymmetric pyrrolo[3,2-b]pyrroles. The utility is verified in the preparation of 13 products that bear six conjugated rings. The observed compatibility with various functional groups allows for increased tunability with regard to the photophysical properties as well as providing sites for further functionalization. Computational studies of the reaction mechanism revealed that the formation of the six-membered rings accompanied with a 1,2-aryl shift is both kinetically and thermodynamically favourable over plausible formation of products containing 7-membered rings. Steady-state UV/Visible spectroscopy reveals that upon photoexcitation, the prepared S-shaped N-doped nanographenes undergo mostly radiative relaxation leading to large fluorescence quantum yields. Their optical properties are rationalized through time-dependent density functional theory calculations. We anticipate that this chemistry will empower the creation of new materials with various functionalities.

Theoretical Insight into the Effect of Phosphorus Oxygenation on Nonradiative Decays: Comparative Analysis of P-Bridged Stilbene Analogs

Incorporation of the phosphorus element into a π-conjugated skeleton offers valuable prospects for adjusting the electronic structure of the resulting functional π-electron systems. Trivalent phosphorus has the potential to decrease the LUMO level through σ*-π* interaction, which is further enhanced by its oxygenation to the pentavalent P center. This study shows that utilizing our computational analysis to examine excited-state dynamics based on radiative/nonradiative rate constants and fluorescence quantum yield (ΦF) is effective for analyzing the photophysical properties of P-containing organic dyes. We theoretically investigate how the trivalent phosphanyl group and pentavalent phosphine oxide moieties affect radiative and nonradiative decay processes. We evaluate four variations of P-bridged stilbene analogs. Our analysis reveals that the primary decay pathway for photoexcited bis-phosphanyl-bridged stilbene is the intersystem crossing (ISC) to the triplet state and nonradiative. The oxidation of the phosphine moiety, however, suppresses the ISC due to the relative destabilization of the triplet states. The calculated rate constants match an increase in experimental ΦF from 0.07 to 0.98, as simulated from 0.23 to 0.94. The reduced HOMO-LUMO gap supports a red shift in the fluorescence spectra relative to the phosphine analog. The thiophene-fused variant with the nonoxidized trivalent P center exhibits intense emission with a high ΦF, 0.95. Our prediction indicates that the ISC transfer is obstructed owing to the relatively destabilized triplet state induced by the thiophene substitution. Conversely, the thiophene-fused analog with the phosphine oxide moieties triggers a high-rate internal conversion mediated by conical intersection, leading to a decreased ΦF.

Benchmarks for transition metal spin-state energetics: why and how to employ experimental reference data?

Accurate prediction of energy differences between alternative spin states of transition metal complexes is essential in computational (bio)inorganic chemistry-for example, in characterization of spin crossover materials and in the theoretical modeling of open-shell reaction mechanisms-but it remains one of the most compelling problems for quantum chemistry methods. A part of this challenge is to obtain reliable reference data for benchmark studies, as even the highest-level applicable methods are known to give divergent results. This Perspective discusses two possible approaches to method benchmarking for spin-state energetics: using either theoretically computed or experiment-derived reference data. With the focus on the latter approach, an extensive general review is provided for the available experimental data of spin-state energetics and their interpretations in the context of benchmark studies, targeting the possibility of back-correcting the vibrational effects and the influence of solvents or crystalline environments. With a growing amount of experience, these effects can be now not only qualitatively understood, but also quantitatively modeled, providing the way to derive nearly chemically accurate estimates of the electronic spin-state gaps to be used as benchmarks and advancing our understanding of the phenomena related to spin states in condensed phases.

A Single Organic Fluorescent Probe for the Discrimination of Dual Spontaneous ROS in Living Organisms: Theoretical Approach

Single-organic-molecule fluorescent probes with double-lock or even multi-lock response modes have attracted the attention of a wide range of researchers. The number of corresponding reports has rapidly increased in recent years. The effective application of the multi-lock response mode single-molecule fluorescent probe has improved the comprehensive understanding of the related targets' functions or influences in pathologic processes. Building a highly efficient functional single-molecule fluorescent probe would benefit the diagnosis and treatment of corresponding diseases. Here, we conducted a theoretical analysis of the synthesizing and sensing mechanism of this kind of functional single-molecule fluorescent probe, thereby guiding the design and building of new efficient probes. In this work, we discuss in detail the electronic structure, electron excitation, and fluorescent character of a recently developed single-molecule fluorescent probe, which could achieve the discrimination and profiling of spontaneous reactive oxygen species (ROS, •OH, and HClO) simultaneously. The theoretical results provide insights that will help develop new tools for fluorescent diagnosis in biological and medical fields.

Quantum-Classical Protocol for Efficient Characterization of Absorption Lineshape and Fluorescence Quenching upon Aggregation: The Case of Zinc Phthalocyanine Dyes

A quantum-classical protocol that incorporates Jahn-Teller vibronic coupling effects and cluster analysis of molecular dynamics simulations is reported, providing a tool for simulations of absorption spectra and ultrafast nonadiabatic dynamics in large molecular photosystems undergoing aggregation in solution. Employing zinc phthalocyanine dyes as target systems, we demonstrated that the proposed protocol provided fundamental information on vibronic, electronic couplings and thermal dynamical effects that mostly contribute to the absorption spectra lineshape and the fluorescence quenching processes upon dye aggregation. Decomposing the various effects arising upon dimer formation, the structure-property relations associated with their optical responses have been deciphered at atomistic resolution.

Curved Nanographenes: Multiple Emission, Thermally Activated Delayed Fluorescence, and Non-Radiative Decay

The intriguing and rich photophysical properties of three curved nanographenes (CNG 6, 7, and 8) are investigated by time-resolved and temperature-dependent photoluminescence (PL) spectroscopy. CNG 7 and 8 exhibit dual fluorescence, as well as dual phosphorescence at low temperature in the main PL bands. In addition, hot bands are detected in fluorescence as well as phosphorescence, and, in the narrow temperature range of 100-140 K, thermally activated delayed fluorescence (TADF) with lifetimes on the millisecond time-scale is observed. These findings are rationalized by quantum-chemical simulations, which predict a single minimum of the S1 potential of CNG 6, but two S1 minima for CNG 7 and CNG 8, with considerable geometric reorganization between them, in agreement with the experimental findings. Additionally, a higher-lying S2 minimum close to S1 is optimized for the three CNG, from where emission is also possible due to thermal activation and, hence, non-Kasha behavior. The presence of higher-lying dark triplet states close to the S1 minima provides mechanistic evidence for the TADF phenomena observed. Non-radiative decay of the T1 state appears to be thermally activated with activation energies of roughly 100 meV and leads to disappearance of phosphorescence and TADF at T > 140 K.

Extensive Analysis of the Parameters Influencing Radiative Rates Obtained through Vibronic Calculations

Defining a theoretical model systematically delivering accurate ab initio predictions of the fluorescence quantum yields of organic dyes is highly desirable for designing improved fluorophores in a systematic rather than trial-and-error way. To this end, the first required step is to obtain reliable radiative rates (kr), as low kr typically precludes effective emission. In the present contribution, using a series of 10 substituted phenyls with known experimental kr, we analyze the impact of the computational protocol on the kr determined through the thermal vibration correlation function (TVCF) approach on the basis of time-dependent density functional theory (TD-DFT) calculations of the energies, structures, and vibrational parameters. Both the electronic structure (selected exchange-correlation functional, application or not of the Tamm-Dancoff approximation) and the vibronic parameters (line-shape formalism, coordinate system, potential energy surface model, and dipole expansion) are tackled. Considering all possible combinations yields more than 3500 cases, allowing to extract statistically-relevant information regarding the impact of each computational parameter on the magnitude of the estimated kr. It turns out that the selected vibronic model can have a significant impact on the computed kr, especially the potential energy surface model. This effect is of the same order of magnitude as the difference noted between B3LYP and CAM-B3LYP estimates. For the treated compounds, all evaluated functionals do deliver reasonable trends, fitting the experimental values.

Multiconfigurational Calculations and Photodynamics Describe Norbornadiene Photochemistry

Storing solar energy is a vital component of using renewable energy sources to meet the growing demands of the global energy economy. Molecular solar thermal (MOST) energy storage is a promising means to store solar energy with on-demand energy release. The light-induced isomerization reaction of norbornadiene (NBD) to quadricyclane (QC) is of great interest because of the generally high energy storage density (0.97 MJ kg-1) and long thermal reversion lifetime (t1/2,300K = 8346 years). However, the mechanistic details of the ultrafast excited-state [2 + 2]-cycloaddition are largely unknown due to the limitations of experimental techniques in resolving accurate excited-state molecular structures. We now present a full computational study on the excited-state deactivation mechanism of NBD and its dimethyl dicyano derivative (DMDCNBD) in the gas phase. Our multiconfigurational calculations and nonadiabatic molecular dynamics simulations have enumerated the possible pathways with 557 S2 trajectories of NBD for 500 fs and 492 S1 trajectories of DMDCNBD for 800 fs. The simulations predicted the S2 and S1 lifetimes of NBD (62 and 221 fs, respectively) and the S1 lifetime of DMDCNBD (190 fs). The predicted quantum yields of QC and DCQC are 10 and 43%, respectively. Our simulations also show the mechanisms of forming other possible reaction products and their quantum yields.

Non-Phenomenological Description of the Time-Resolved Emission in Solution with Quantum-Classical Vibronic Approaches-Application to Coumarin C153 in Methanol

We report a joint experimental and theoretical work on the steady-state spectroscopy and time-resolved emission of the coumarin C153 dye in methanol. The lowest energy excited state of this molecule is characterized by an intramolecular charge transfer thus leading to remarkable shifts of the time-resolved emission spectra, dictated by the methanol reorganization dynamics. We selected this system as a prototypical test case for the first application of a novel computational protocol aimed at the prediction of transient emission spectral shapes, including both vibronic and solvent effects, without applying any phenomenological broadening. It combines a recently developed quantum-classical approach, the adiabatic molecular dynamics generalized vertical Hessian method (Ad-MD|gVH), with nonequilibrium molecular dynamics simulations. For the steady-state spectra we show that the Ad-MD|gVH approach is able to reproduce quite accurately the spectral shapes and the Stokes shift, while a ∼0.15 eV error is found on the prediction of the solvent shift going from gas phase to methanol. The spectral shape of the time-resolved emission signals is, overall, well reproduced, although the simulated spectra are slightly too broad and asymmetric at low energies with respect to experiments. As far as the spectral shift is concerned, the calculated spectra from 4 ps to 100 ps are in excellent agreement with experiments, correctly predicting the end of the solvent reorganization after about 20 ps. On the other hand, before 4 ps solvent dynamics is predicted to be too fast in the simulations and, in the sub-ps timescale, the uncertainty due to the experimental time resolution (300 fs) makes the comparison less straightforward. Finally, analysis of the reorganization of the first solvation shell surrounding the excited solute, based on atomic radial distribution functions and orientational correlations, indicates a fast solvent response (≈100 fs) characterized by the strengthening of the carbonyl-methanol hydrogen bond interactions, followed by the solvent reorientation, occurring on the ps timescale, to maximize local dipolar interactions.

A direct observation of up-converted room-temperature phosphorescence in an anti-Kasha dopant-matrix system

It is common sense that emission maxima of phosphorescence spectra (λ_P) are longer than those of fluorescence spectra (λ_F). Here we report a serendipitous finding of up-converted room-temperature phosphorescence (RTP) with λ_P < λ_F and phosphorescence lifetime > 0.1 s upon doping benzophenone-containing difluoroboron β-diketonate (BPBF₂) into phenyl benzoate matrices. The up-converted RTP is originated from BPBF₂'s T_n (n ≥ 2) states which show typical ³n-π* characters from benzophenone moieties. Detailed studies reveal that, upon intersystem crossing from BPBF₂'s S₁ states of charge transfer characters, the resultant T₁ and T_n states build T₁-to-T_n equilibrium. Because of their ³n-π* characters, the T_n states possess large phosphorescence rates that can strongly compete RTP(T₁) to directly emit RTP(T_n) which violates Kasha's rule. The direct observation of up-converted RTP provides deep understanding of triplet excited state dynamics and opens an intriguing pathway to devise visible-light-excitable deep-blue afterglow emitters, as well as stimuli-responsive afterglow materials.

Cost-Effective Simulations of Vibrationally-Resolved Absorption Spectra of Fluorophores with Machine-Learning-Based Inhomogeneous Broadening

The results of electronic and vibrational structure simulations are an invaluable support for interpreting experimental absorption/emission spectra, which stimulates the development of reliable and cost-effective computational protocols. In this work, we contribute to these efforts and propose an efficient first-principle protocol for simulating vibrationally-resolved absorption spectra, including nonempirical estimations of the inhomogeneous broadening. To this end, we analyze three key aspects: (i) a metric-based selection of density functional approximation (DFA) so to benefit from the computational efficiency of time-dependent density function theory (TD-DFT) while safeguarding the accuracy of the vibrationally-resolved spectra, (ii) an assessment of two vibrational structure schemes (vertical gradient and adiabatic Hessian) to compute the Franck-Condon factors, and (iii) the use of machine learning to speed up nonempirical estimations of the inhomogeneous broadening. In more detail, we predict the absorption band shapes for a set of 20 medium-sized fluorescent dyes, focusing on the bright ππ^★ S₀ → S₁ transition and using experimental results as references. We demonstrate that, for the studied 20-dye set which includes structures with large structural variability, the preselection of DFAs based on an easily accessible metric ensures accurate band shapes with respect to the reference approach and that range-separated functionals show the best performance when combined with the vertical gradient model. As far as band widths are concerned, we propose a new machine-learning-based approach for determining the inhomogeneous broadening induced by the solvent microenvironment. This approach is shown to be very robust offering inhomogeneous broadenings with errors as small as 2 cm^-1 with respect to genuine electronic-structure calculations, with a total CPU time reduced by 98%.