Two-Photon Absorption Strengths of Small Molecules: Reference CC3 Values and Benchmarks.

We present a large dataset of highly accurate two-photon transition strengths (δTPA) determined for standard small molecules. Our reference values have been calculated using the quadratic response implementation of the third-order coupled cluster method including iterative triples (Q-CC3). The aug-cc-pVTZ atomic basis set is used for molecules with up to five non-hydrogen atoms, while larger molecules are assessed with aug-cc-pVDZ; the differences due to the basis sets are discussed. This dataset, encompassing 82 singlet transitions of various characters (Rydberg, valence, and double excitations), enables a comprehensive benchmark of smaller basis sets and alternative wavefunction methods when Q-CC3 calculations become beyond reach as well as time-dependent density functional theory (TD-DFT) approaches. The evaluated wavefunction methods include quadratic response and equation-of-motion CCSD approximations, Q-CC2, and second-order algebraic diagrammatic construction in its intermediate state representation (I-ADC2). In the TD-DFT framework, a set of five commonly used exchange-correlation functionals are evaluted. This extensive analysis provides a quantitative assessment of these methods, revealing how different system sizes, response intensities, and types of transitions affect their performances.

Dual luminescence and infrared circularly polarized luminescence up to 900 nm with platinum complexes bearing a helical donor-acceptor ligand.

Chiral molecular materials able to emit circularly polarized luminescence (CPL) have attracted considerable interest in the last few decades, due to the potential of CP-light in a wide range of applications. While CP luminescent molecules with blue, green, and yellow emissions are now well-reported, NIR CPL from organic and organometallic compounds lags behind due to the dual challenge of promoting radiative deexcitation of the excited state in this low energy region while assuring a significant magnetic dipole transition moment, a prerequisite for generating CPL. Based on a versatile axially chiral arylisoquinoline ligand, we report the synthesis and chiroptical properties of chiral donor-acceptor platinum(ii) complexes displaying CPL that extends up to almost 900 nm. Interestingly, these emitters show both fluorescence and phosphorescence emissions in solution, with intensities depending on the charge-transfer character of the organic ligand. Experimental and theoretical investigations show that this feature strongly impacts the intersystem crossing event between the singlet and triplet excited states of these complexes and the related phosphorescence lifetime. The effect is less important regarding the CPL, and most complexes show luminescence dissymmetry factors with values up to ca. 2 × 10-3 around 800 nm.

Light- and Temperature-Controlled Hybridization, Chiral Induction and Handedness of Helical Foldamers.

Helical foldamers have attracted much attention over the last decades given their resemblance to certain biomacromolecules and their potential in domains as different as pharmaceutics, catalysis and photonics. Various research groups have successfully controlled the right- or left- handedness of these oligomers by introducing stereogenic centers through covalent or non-covalent chemistry. However, developing helical structures whose handedness can be reversibly switched remains a major challenge for chemists. To date, such an achievement has been reported with light-responsive single-stranded foldamers only. Herein, we demonstrate that grafting a unidirectional motor onto foldamer strands constitutes a relevant strategy to i) control the single or double helical state of a foldamer, ii) switch on the chiral induction process from the motor to the helical strands and iii) select the handedness of double helical structures through photochemical and thermal stimulations.

Reference CC3 Excitation Energies for Organic Chromophores: Benchmarking TD-DFT, BSE/GW, and Wave Function Methods.

To expand the QUEST database of highly accurate vertical transition energies, we consider a series of large organic chromogens ubiquitous in dye chemistry, such as anthraquinone, azobenzene, BODIPY, and naphthalimide. We compute, at the CC3 level of theory, the singlet and triplet vertical transition energies associated with the low-lying excited states. This leads to a collection of more than 120 new highly accurate excitation energies. For several singlet transitions, we have been able to determine CCSDT transition energies with a compact basis set, finding minimal deviations from the CC3 values for most states. Subsequently, we employ these reference values to benchmark a series of lower-order wave function approaches, including the popular ADC(2) and CC2 schemes, as well as time-dependent density-functional theory (TD-DFT), both with and without applying the Tamm-Dancoff approximation (TDA). At the TD-DFT level, we evaluate a large panel of global, range-separated, local, and double hybrid functionals. Additionally, we assess the performance of the Bethe-Salpeter equation (BSE) formalism relying on both G0W0 and evGW quasiparticle energies evaluated from various starting points. It turns out that CC2 and ADC(2.5) are the most accurate models among those with respective O(N5) and O(N6) scalings with system size. In contrast, CCSD does not outperform CC2. The best performing exchange-correlation functionals include BMK, M06-2X, M06-SX, CAM-B3LYP, ωB97X-D, and LH20t, with average deviations of approximately 0.20 eV or slightly below. Errors on vertical excitation energies can be further reduced by considering double hybrids. Both SOS-ωB88PP86 and SOS-ωPBEPP86 exhibit particularly attractive performances with overall quality on par with CC2, whereas PBE0-DH and PBE-QIDH are only slightly less efficient. BSE/evGW calculations based on Kohn-Sham starting points have been found to be particularly effective for singlet transitions, but much less for their triplet counterparts.

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.

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.

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.

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.

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.

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.

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.

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.