Stimuli-Induced Fluorescence Switching in Azine-Containing Fluorophores Displaying Resonance-Stabilized ESIPT Emission.

This article reports the synthesis, along with structural and photophysical characterization of 2-(2'-hydroxyphenyl)benzazole derivatives functionalized with various azaheterocycles (pyridine, pyrimidine, terpyridine). These compounds show dual-state emission properties, that is intense fluorescence both in solution and in the solid-state with a range of fluorescent color going from blue to orange. Moreover, the nature of their excited state can be tuned by the presence of external stimuli such as protons or metal cations. In the absence of stimuli, these dyes show emission stemming from anionic species obtained after deprotonation (D* transition), whereas upon protonation or metal chelation, ESIPT process occurs leading to a stabilized and highly emissive K* transition. With the help of extensive ab initio calculations, we confirm that external stimuli can switch the nature of the transitions, making this series of dyes attractive candidates for the development of stimuli-responsive fluorescent ratiometric probes.

Reference Energies for Double Excitations: Improvement and Extension.

In the realm of photochemistry, the significance of double excitations (also known as doubly excited states), where two electrons are concurrently elevated to higher energy levels, lies in their involvement in key electronic transitions essential in light-induced chemical reactions as well as their challenging nature from the computational theoretical chemistry point of view. Based on state-of-the-art electronic structure methods (such as high-order coupled-cluster, selected configuration interaction, and multiconfigurational methods), we improve and expand our prior set of accurate reference excitation energies for electronic states exhibiting a substantial amount of double excitations [Loos et al. J. Chem. Theory Comput. 2019, 15, 1939]. This extended collection encompasses 47 electronic transitions across 26 molecular systems that we separate into two distinct subsets: (i) 28 "genuine" doubly excited states where the transitions almost exclusively involve doubly excited configurations and (ii) 19 "partial" doubly excited states which exhibit a more balanced character between singly and doubly excited configurations. For each subset, we assess the performance of high-order coupled-cluster (CC3, CCSDT, CC4, and CCSDTQ) and multiconfigurational methods (CASPT2, CASPT3, PC-NEVPT2, and SC-NEVPT2). Using as a probe the percentage of single excitations involved in a given transition (%T1) computed at the CC3 level, we also propose a simple correction that reduces the errors of CC3 by a factor of 3, for both sets of excitations. We hope that this more complete and diverse compilation of double excitations will help future developments of electronic excited-state methodologies.

Post-functionalization of triamino-phenazinium dyes to reach near-infrared emission.

This study presents the synthesis and characterization of phenazinium dyes with absorption ranging from red to far-red, as well as emission extending into the far-red to near-infrared (NIR) region. The procedure involves the post-functionalization of a triamino-phenazinium that was recently reported as a theranostic agent. The introduction of electron-withdrawing moieties is accomplished through acylation or aromatic nucleophilic substitution. For one of the obtained products, a further substitution step could be achieved with primary amines to tune the electron density of the phenazinium core. The isolated dyes exhibit promising features that hold potential for future applications as biological markers or therapeutic agents.

State-Specific Coupled-Cluster Methods for Excited States.

We reexamine ΔCCSD, a state-specific coupled-cluster (CC) with single and double excitations (CCSD) approach that targets excited states through the utilization of non-Aufbau determinants. This methodology is particularly efficient when dealing with doubly excited states, a domain in which the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here to evaluate the effectiveness of ΔCCSD when applied to other types of excited states, comparing its consistency and accuracy with EOM-CCSD. To this end, we report a benchmark on excitation energies computed with the ΔCCSD and EOM-CCSD methods for a set of molecular excited-state energies that encompasses not only doubly excited states but also doublet-doublet transitions and (singlet and triplet) singly excited states of closed-shell systems. In the latter case, we rely on a minimalist version of multireference CC known as the two-determinant CCSD method to compute the excited states. Our data set, consisting of 276 excited states stemming from the quest database [Véril et al., WIREs Comput. Mol. Sci. 2021, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of ΔCCSD. Except for the doubly excited states, we found that ΔCCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of 0.10 and 0.07 eV) is less pronounced than that obtained for singly excited states of closed-shell systems (MAEs of 0.15 and 0.08 eV). This discrepancy is largely attributed to a greater number of excited states in the latter set exhibiting multiconfigurational characters, which are more challenging for ΔCCSD. We also found typically small improvements by employing state-specific optimized orbitals.

Experimental and theoretical comprehension of ESIPT fluorophores based on a 2-(2'-hydroxyphenyl)-3,3'-dimethylindole (HDMI) scaffold.

Excited-State Intramolecular Proton Transfer (ESIPT) emission is associated with intense single or multiple fluorescence in the solid-state, along with enhanced photostability and sensitivity to the close environment. As a result, ESIPT probes are attractive candidates for ratiometric sensing of a variety of substrates. A new family of ESIPT fluorophores is described herein, inspired by the well-known 2-(2'hydroxyphenyl)benzazole (HBX) organic scaffold. The connection of 3,3'-dimethylindole (or 3H-indole) derivatives with phenol rings triggers the formation of novel 2-(2'-hydroxyphenyl)-3,3'-dimethylindole (HDMI) fluorophores, capable of stimuli-responsive ESIPT emission. This brand new family of dyes displays redshifted emission, as compared to HBX, along with an unprecedented acid/base-mediated stabilization of different rotamers, owing to supramolecular interactions with methyl groups. These compounds are therefore highly sensitive to external stimuli, such as the presence of acid or base, where protonated and deprotonated species have specific optical signatures. Moreover, a new pyridine-functionalized HDMI dye displays acid-sensitive AIE properties. The photophysical properties of all compounds have also been studied using ab initio calculations to support experiments in deciphering the nature of the various radiative transitions observed and the related excited rotameric species.

Quadrupolar, Highly Polarized Dyes: Emission Dependence on Viscosity and Selective Mitochondria Staining.

Quadrupolar A-D-A-type 1,4-dihydropyrrolo[3,2-b]pyrroles (DHPPs) bearing pyridinium and quinolinium substituents emit in the 500-600 nm region. The enhancement of electronic communication between the electron-rich heterocyclic core and electron-deficient peripheral substituents turned out to be crucial for achieving emission enhancement in viscous media. DHPP bearing two 4-pyridinium substituents has optical brightness 34,000 in glycerol and only 700 in MeOH, as evidenced by measurements of the emission intensity and fluorescence lifetimes in a series of polar solvents. Such behavior makes it an excellent candidate for viscosity probes in fluorescence microscopy, as demonstrated by the fluorescence imaging of H9C2 cardiomyocytes.

Assessing the accuracy of TD-DFT excited-state geometries through optimal tuning with GW energy levels.

We study the accuracy of excited state (ES) geometries using optimally tuned LC-PBE functionals with tuning based on GW quasiparticle energies. We compare the results obtained with the PBE, PBE0, non-tuned, and tuned LC-PBE functionals with available high-level CC reference values as well as experimental data. First, we compare ES geometrical parameters obtained for three different types of systems: molecules composed of a few atoms, 4-(dimethylamino)benzonitrile (DMABN), and conjugated dyes. To this end, we used wave-function results as benchmarks. Next, we evaluate the accuracy of the theoretically simulated spectra as compared to the experimental ones for five large dyes. Our results show that, besides small compact molecules for which tuning LC-PBE does not allow obtaining geometries more accurate than those computed with standard functionals, tuned range-separated functionals are clearly to be favored, not only for ES geometries but also for 0-0 energies, band shapes, and intensities for absorption and emission spectra. In particular, the results indicate that GW-tuned LC-PBE functionals provide improved matching with experimental spectra as compared to conventionally tuned functionals. It is an open question whether TD-DFT with GW-tuned functionals can qualitatively mimic the actual many-body Bethe-Salpeter (BSE/GW) formalism for which analytic ionic gradients remain to be developed.

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.

A mountaineering strategy to excited states: Accurate vertical transition energies and benchmarks for substituted benzenes.

In an effort to expand the existing QUEST database of accurate vertical transition energies [Véril et al. WIREs Comput. Mol. Sci. 2021, 11, e1517], we have modeled more than 100 electronic excited states of different natures (local, charge-transfer, Rydberg, singlet, and triplet) in a dozen of mono- and di-substituted benzenes, including aniline, benzonitrile, chlorobenzene, fluorobenzene, nitrobenzene, among others. To establish theoretical best estimates for these vertical excitation energies, we have employed advanced coupled-cluster methods including iterative triples (CC3 and CCSDT) and, when technically possible, iterative quadruples (CC4). These high-level computational approaches provide a robust foundation for benchmarking a series of popular wave function methods. The evaluated methods all include contributions from double excitations (ADC(2), CC2, CCSD, CIS(D), EOM-MP2, STEOM-CCSD), along with schemes that also incorporate perturbative or iterative triples (ADC(3), CCSDR(3), CCSD(T)(a) ⋆ , and CCSDT-3). This systematic exploration not only broadens the scope of the QUEST database but also facilitates a rigorous assessment of different theoretical approaches in the framework of a homologous chemical series, offering valuable insights into the accuracy and reliability of these methods in such cases. We found that both ADC(2.5) and CCSDT-3 can provide very consistent estimates, whereas among less expensive methods SCS-CC2 is likely the most effective approach. Importantly, we show that some lower order methods may offer reasonable trends in the homologous series while providing quite large average errors, and vice versa. Consequently, benchmarking the accuracy of a model based solely on absolute transition energies may not be meaningful for applications involving a series of similar compounds.

Difluoro Dipyridomethene Boron Complexes: Synthesis, Characterization, and Ab Initio Calculations.

Molecular engineering studies on the meso-cyano difluoro dipyridomethene boron complexes are presented and two series (a and b) of novel fluorophores are extensively studied. Halogenated derivatives were reacted under Suzuki-Miyaura or Sonogashira cross coupling reactions to introduce electron-donating or electron-withdrawing functional groups on positions 1 and 2 of the aromatic ligand. All derivatives were obtained in 14-90% yields and studied in detail by structural, photophysical, and computational analyses. Both series display excellent emissive properties in solution with blue to orange fluorescence emission upon blue light absorption and promising features as solid emitters. All the spectroscopic measurements are supported and confirmed by first-principles theoretical calculations combining TD-DFT and CC2. Series b, featuring an aryl substituent onto position 1 of the aromatic core, showed significantly large Stokes shifts values.

Azaborahelicene fluorophores derived from four-coordinate N,C-boron chelates: synthesis, photophysical and chiroptical properties.

A series of six azaborahelicenes with varying electron-donor substitution at the 4-position of the aryl residue (i.e., naphthyl) or with variable π-extension of the aryl residue (thianthrenyl, anthryl, pyrenyl) was prepared with an efficient and flexible synthetic protocol. These different types of functionalization afforded notably pronounced intramolecular charge-transfer (ICT) character for the dyes with the strongest electron donor substitution (NMe2) or easiest to oxidize aryl residues, as evidenced by photophysical investigations. These effects also impact the corresponding chiroptical properties of the separated M- and P-enantiomers, which notably display circularly polarized luminescence (CPL) with dissymmetry factors in the order of magnitude of 10-4 to 10-3. Theoretical calculations confirm the optical spectroscopy data and are in agreement with the proposed involvement of ICT processes.

Insights into extended coupled polymethines through the investigation of dual UV-to-NIR acidochromic switches based on heptamethine-oxonol dyes.

A series of heptamethine-oxonol dyes featuring different heterocyclic end groups were designed with the aim to explore structure-property relationships in π-extended coupled polymethines. These dyes can be stabilised under three different protonation states, affording dicationic derivatives with an aromatic core, cationic heptamethines, and zwitterionic bis-cyanine forms. The variation of the end groups directly impacts the absorption and emission properties and mostly controls reaching either a colourless neutral dispirocyclic species or near-infrared zwitterions. The acidochromic switching between the three states involves profound electronic rearrangements leading to notable shifts of their optical properties that were investigated using a parallel experiment-theory approach, providing a comprehensive description of these unique systems.

Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study.

Nonadiabatic dynamics simulations have become a standard approach to explore photochemical reactions. Such simulations require underlying potential energy surfaces and couplings between them, calculated at a chosen level of theory, yet this aspect is rarely assessed. Here, in combination with the popular trajectory surface hopping dynamics method, we use a high-accuracy XMS-CASPT2 electronic structure level as a benchmark for assessing the performances of various post-Hartree-Fock methods (namely, CIS, ADC(2), CC2, and CASSCF) and exchange-correlation functionals (PBE, PBE0, and CAM-B3LYP) in a TD-DFT/TDA context, using the isomerization around a double bond as test case. Different relaxation pathways are identified, and the ability of the different methods to reproduce their relative importance and time scale is discussed. The results show that multireference electronic structure methods should be preferred, when studying nonadiabatic decay between excited and ground states. If not affordable, TD-DFT with TDA and hybrid functionals and ADC(2) are efficient alternatives but overestimate the nonradiative decay yield and thus may miss deexcitation pathways.

Non-uniform sampling to enhance the performance of compact NMR for characterizing new psychoactive substances.

Efficient and robust analytical methods are needed to improve the identification and subsequent regulation of new psychoactive substances (NPS). NMR spectroscopy is a unique method able to determine the structure of small molecules such as NPS even in mixtures. However, high-field NMR analysis is associated with expensive purchase and maintenance costs. For more than a decade, compact NMR spectrometers have changed this paradigm. It was recently shown that a dedicated analytical workflow combining compact NMR and databases could identify the molecular structure of NPS, in spite of the lower spectral dispersion and sensitivity of compact spectrometers. This approach relies on 1H-13C HSQC to both recognize NPS and elucidate the structure of unknown substances. Still, its performance is limited by the need to compromise between resolution and experiment time. Here, we show that this strategy can be significantly improved by implementing non-uniform sampling (NUS) to improve spectral resolution in the 13C dimension of HSQC at no cost in terms of experiment time. Gains in the range of 3 to 4 in resolution are achieved for pure NPS and for a mixture. Finally, 2D HSQC with NUS was applied to improve the identification of NPS with the assistance of databases. The resulting method appears as a useful tool for the characterization of NPS in mixtures, which is essential for forensic laboratories.

Heptazine, Cyclazine, and Related Compounds: Chemically-Accurate Estimates of the Inverted Singlet-Triplet Gap.

Molecules that violate Hund's rule and exhibit an inverted gap between the lowest singlet S1 and triplet T1 excited states have attracted considerable attention due to their potential applications in optoelectronics. Among these molecules, the triangular-shaped heptazine, and its derivatives, have been in the limelight. However, conflicting reports have arisen regarding the relative energies of S1 and T1. Here, we employ highly accurate levels of theory, such as CC3, to not only resolve the debate concerning the sign but also quantify the magnitude of the S1-T1 gap. We also determined the 0-0 energies to evaluate the significance of the vertical approximation. In addition, we compute reference S1-T1 gaps for a series of 10 related molecules. This enables us to benchmark lower-order methods for future applications in larger systems within the same family of compounds. This contribution can serve as a foundation for the design of triangular-shaped molecules with enhanced photophysical properties.

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.

Phospha-cyanines in Their Ideal Polymethine State: Synthesis and Structure-Property Relationships.

We report the synthesis and full characterization of a family of phosphorus-containing polymethine cyanines (phospha-cyanines). The compounds are easily prepared in two steps, starting from readily available phosphanes. The impact of the P-substituents and the counterions on the structural and optical properties is investigated through a joint experimental/theoretical approach. Based on the study of the single-crystal X-ray diffraction structures, all phospha-cyanines present a bond length alternation close to zero, independently of the substituents and counterions, which indicates an ideal polymethine state. All these compounds display the typical cyanine-like UV-vis absorption with an intense and sharp transition with a vibronic shoulder despite possessing a reverse electronic configuration compared to "classical" cyanines. Time-dependent density-functional theory calculations allowed us to fully rationalize the optical properties (absorption/emission wavelengths, luminescence quantum yields). Interestingly, due to the tetrahedral shape of the P atom, the optical properties are independent of the counterion, which is in marked contrast with N-analogues, which enables predictive engineering of the phospha-cyanines regardless of the medium in which they are used.

Reference Vertical Excitation Energies for Transition Metal Compounds.

To enrich and enhance the diversity of the quest database of highly accurate excitation energies [Véril, M.; et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], we report vertical transition energies in transition metal compounds. Eleven diatomic molecules with a singlet or doublet ground state containing a fourth-row transition metal (CuCl, CuF, CuH, ScF, ScH, ScO, ScS, TiN, ZnH, ZnO, and ZnS) are considered, and the corresponding excitation energies are computed using high-level coupled-cluster (CC) methods, namely, CC3, CCSDT, CC4, and CCSDTQ, as well as multiconfigurational methods such as CASPT2 and NEVPT2. In many cases, to provide more comprehensive benchmark data, we also provide full configuration interaction estimates computed with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. Based on these calculations, theoretical best estimates of the transition energies are established in both the aug-cc-pVDZ and aug-cc-pVTZ basis sets. This allows us to accurately assess the performance of the CC and multiconfigurational methods for this specific set of challenging transitions. Furthermore, comparisons with experimental data and previous theoretical results are also reported.

The magic of biaryl linkers: the electronic coupling through them defines the propensity for excited-state symmetry breaking in quadrupolar acceptor-donor-acceptor fluorophores.

Charge transfer (CT) is key for molecular photonics, governing the optical properties of chromophores comprising electron-rich and electron-deficient components. In photoexcited dyes with an acceptor-donor-acceptor or donor-acceptor-donor architecture, CT breaks their quadrupolar symmetry and yields dipolar structures manifesting pronounced solvatochromism. Herein, we explore the effects of electronic coupling through biaryl linkers on the excited-state symmetry breaking of such hybrid dyes composed of an electron-rich core, i.e., 1,4-dihydropyrrolo[3,2-b]pyrrole (DHPP), and pyrene substituents that can act as electron acceptors. Experimental and theoretical studies reveal that strengthening the donor-acceptor electronic coupling decreases the CT rates and the propensity for symmetry breaking. We ascribe this unexpected result to effects of electronic coupling on the CT thermodynamics, which in its turn affects the CT kinetics. In cases of intermediate electronic coupling, the pyrene-DHPP conjugates produce fluorescence spectra, spreading over the whole visible range, that in addition to the broad CT emission, show bands from the radiative deactivation of the locally excited states of the donor and the acceptors. Because the radiative deactivation of the low-lying CT states is distinctly slow, fluorescence from upper locally excited states emerge leading to the observed anti-Kasha behaviour. As a result, these dyes exhibit white fluorescence. In addition to demonstrating the multifaceted nature of the effects of electronic coupling on CT dynamics, these chromophores can act as broad-band light sources with practical importance for imaging and photonics.

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.