Taking up the cyanine challenge with quantum tools

Exploring the Te(II)/Te(IV) Redox Couple of a Tellurorosamine Chromophore: Photophysical, Photochemical, and Electrochemical Studies

A tellurorosamine dye [Te(II)] undergoes aerobic photooxidation. Although Te(IV) species have been used in a number of oxidations, key Te(IV)-oxo and Te(IV)-bis(hydroxy) intermediates are challenging to study. Under aerobic irradiation with visible light, Te(II) (λmax = 600 nm) transforms into a Te(IV) species (λmax = 669 nm). The resultant Te(IV) species is not stable in the dark or at -20 °C, decomposing back to Te(II) and other byproducts over many hours. To eliminate the structural ambiguity of the Te(IV) photoproduct, we used spectroelectrochemistry, wherein the bis(hydroxy) Te(IV)-(OH)2 was electrochemically generated under anaerobic conditions. The absorption of Te(IV)-(OH)2 matches that of the Te(IV) photoproduct. Because isosbestic points are maintained both photochemically and electrochemically, the oxo core formed photochemically must rapidly equilibrate with Te(IV)-(OH)2. Calculations on the bis(hydroxy) versus oxo species further corroborate that the equilibration is rapid and the spectra of the two species are similar. To further explore Te(IV) cores, two novel compounds, Te(IV)-Cl2 and Te(IV)-Br2, were synthesized. Characterization of Te(IV)-X2 was simplified because these cores have no analogue to the Te(IV)-(O)/Te(IV)-(OH)2 equilibrium. This work provides insights into the photophysical and electrochemical behavior of Te analogues of chalcogenoxanthylium dyes, which are relevant for a broad range of photochemical applications.

Conformational and environmental effects on the electronic and vibrational properties of dyes for solar cell devices

The main challenge for solar cell devices is harvesting photons beyond the visible by reaching the red-edge (650-780 nm). Dye-sensitized solar cell (DSSC) devices combine the optical absorption and the charge separation processes by the association of a sensitizer as a light-absorbing material (dye molecules, whose absorption can be tuned and designed) with a wide band gap nanostructured semiconductor. Conformational and environmental effects (i.e., solvent, pH) can drastically influence the photophysical properties of molecular dyes. This study proposes a combined experimental and computational approach for the comprehensive investigation of the electronic and vibrational properties of a unique class of organic dye compounds belonging to the family of red-absorbing dyes, known as squaraines. Our focus lies on elucidating the intricate interplay between the molecular structure, vibrational dynamics, and optical properties of squaraines using state-of-the-art density functional theory calculations and spectroscopic techniques. Through systematic vibrational and optical analyses, we show that (i) the main absorption peak in the visible range is influenced by the conformational and protonation equilibria, (ii) the solvent polarity tunes the position of the UV-vis absorption, and (iii) the vibrational spectroscopy techniques (infrared and Raman) can be used as informative tools to distinguish between different conformations and protonation states. This comprehensive understanding offers valuable insights into the design and optimization of squaraine-based DSSCs for enhanced solar energy conversion efficiency.

Synthesis and Regioselective Functionalization of Tetrafluorobenzo-[α]-Fused BOPYPY Dyes

The synthesis of a new bis-BF2 tetrafluorobenzo-[α]-fused BOPYPY dye from 4,5,6,7-tetrafluoroisoindole and 2-hydrazinopyrazine is reported. The regioselectivity of nucleophilic substitution reactions at the periphery of the tetrafluorinated BOPYPY and its α-bromo derivative were investigated using N-, O-, S-, and C-based nucleophiles. Among the aromatic fluorine atoms, the F2 atom is consistently regioselectively substituted, except when the α-position contains a thiophenol group; in this case, F4 is substituted instead due to stabilizing π-π-stacking between the two aromatic groups. The α-bromo BOPYPY derivative also reacts under Stille cross-coupling reaction conditions to produce the corresponding α-substituted product. The spectroscopic properties of these new fluorinated BOPYPYs were investigated and compared with nonfluorinated analogs.

Synthesis, Structures and Photophysical Properties of Asymmetric Fulvene-[b]-fused BODIPYs

Herein, we developed a one-pot procedure to synthesize novel fulvene-[b]-fused BODIPYs from α-(2-alkynylphenyl)-pyrrole and acylpyrrole, using 5-exo cyclization as the key transformation. Compared to benzene-[b]-fused BODIPYs, although they have similar chemical compositions, their structures and properties significantly differ from each other, which can be attributed to the less aromaticity of the fulvene linker than benzene. Notably, fulvene-[b]-fused BODIPY 1 b exhibits helical-twisted core skeleton, intensified red-shifted absorption, and peak fluorescence. In addition, the pathway of this one-pot reaction and the mechanism of POCl3 mediated 5-exo cyclization have been proposed by a combining experimental and computational study.

Electrochemical Diselenation of BODIPY Fluorophores for Bioimaging Applications and Sensitization of 1 O2

We report a rapid, efficient, and scope-extensive approach for the late-stage electrochemical diselenation of BODIPYs. Photophysical analyses reveal red-shifted absorption - corroborated by TD-DFT and DLPNO-STEOM-CCSD computations - and color-tunable emission with large Stokes shifts in the selenium-containing derivatives compared to their precursors. In addition, due to the presence of the heavy Se atoms, competitive ISC generates triplet states which sensitize 1 O2 and display phosphorescence in PMMA films at RT and in a frozen glass matrix at 77 K. Importantly, the selenium-containing BODIPYs demonstrate the ability to selectively stain lipid droplets, exhibiting distinct fluorescence in both green and red channels. This work highlights the potential of electrochemistry as an efficient method for synthesizing unique emission-tunable fluorophores with broad-ranging applications in bioimaging and related fields.

Cyclopentene ring effects in cyanine dyes: a handle to fine-tune photophysical properties

The aim of this study is to investigate the photophysical properties of a cyanine dye analogue by performing first-principles calculations based on density functional theory (DFT) and time dependent-DFT. Cationic cyanine dyes are the subject of great importance due to their versatile applications and the tunability of their photophysical properties, such as by modifying their end groups and chain length. An example of this is the vinylene shift, which is experimentally known for these molecules, and it consists of a bathochromic (red) shift of approximately 100 nm of the 0-0 vibronic transition when a vinyl group is added to the polymethine chain. Our study shows that when the saturated moiety C2H4 of the cyclopentene ring is added to the chain, it interacts with the conjugated π-system, resulting in a smaller HOMO-LUMO gap. Here, we demonstrate the origin of this interaction and how it can be used to fine tune the absorption energies of this class of dyes.

First-Principles Investigation of the Optical Properties of Eumelanin Protomolecules

Using first-principles calculations, we investigate the absorption spectra (in the near-infrared, visible, and first UV range) of the two most probable eumelanin tetrameric molecules exhibiting either a linear open-chain or a cyclic porphyrine-like configuration. In order to simulate a realistic molecular system, an implicit solvent model is used in our calculations to mimic the effect of the solvated environment around the eumelanin molecule. Although the presence of solvent is found not to significantly affect the absorption pattern of both molecules, the onset of the spectra are shifted toward higher energies, especially for the linear tetramer. Interestingly, the absorption spectra and optical onsets of the two molecules differ significantly both in a vacuum and in ethanol. However, the two predicted spectra do not allow us to definitely discriminate between the two configurations when comparing the theoretical predictions with the available experimental spectrum. In addition, a mix of the two eumelanin configurations (close to fifty-fifty) leads to a maximum overlap between theoretical and experimental spectra. Consequently, this theoretical research shows that deeper insight can be gained using beyond DFT techniques on the real form of eumelanin protomolecules present in living systems as well as on their possible use in hybrid solar cells.

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.

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.

Pitfall in simulations of vibronic TD-DFT spectra: diagnosis and assessment

In this Communication, we study the effect of spurious oscillations in the profiles of energy derivatives with respect to nuclear coordinates calculated with density functional approximations (DFAs) for formaldehyde, pyridine, and furan in their ground and electronic excited states. These spurious oscillations, which can only be removed using extensive integration grids that increase enormously the CPU cost of DFA calculations, are significant in the case of third- and fourth-order energy derivatives of the ground and excited states computed by M06-2X and ωB97X functionals. The errors in question propagate to anharmonic vibronic spectra computed under the Franck-Condon approximation, i.e., positions and intensities of vibronic transitions are affected to a large extent (shifts as significant as hundreds of cm-1 were observed). On the other hand, the LC-BLYP and CAM-B3LYP functionals show a much less pronounced effect due to spurious oscillations. Based on the results presented herein, we recommend either LC-BLYP or CAM-B3LYP with integration grids (250, 974) (or larger) for numerically stable simulations of vibronic spectra including anharmonic effects.

Tris-benzo[cd]indole Cyanine Enables the NIR-photosensitized Radical and Thiol-ene Polymerizations at 940 nm

A near-infrared-absorbing heptamethine (HM+ ) incorporating three bulky benzo[cd]indole heterocycles was designed to efficiently prevent self-aggregation of the dye, which results in a strong enhancement of its photoinitiating reactivity as compared to a parent bis-benzo[cd]indole heptamethine (HMCl+ ) used as a reference system. In this context, we highlight an efficient free-radical NIR-polymerization up to a 100 % acrylates C=C bonds conversion even under air conditions. Such an important initiating performance was obtained by incorporating our NIR-sensitizer into a three-component system leading to its self-regeneration. This original photoredox cycle was thoroughly investigated through the identification of each intermediary species using EPR spectroscopy.

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.

Near-infrared voltage-sensitive dyes based on chromene donor

Voltage-sensitive dyes (VSDs) are used to image electrical activity in cells and tissues with submillisecond time resolution. Most of these fast sensors are constructed from push-pull chromophores whose fluorescence spectra are modulated by the electric field across the cell membrane. It was found that the substitution of naphthalene with chromene produces a 60 to 80 nm red-shift in absorption and emission spectra while maintaining fluorescence quantum efficiency and voltage sensitivity. One dye was applied to ex vivo murine heart with excitation at 730 nm, by far the longest wavelength reported in voltage imaging. This VSD resolves cardiac action potentials in single trials with 12% ΔF/F per action potential. The well-separated excitation spectra between these long-wavelength VSDs and channelrhodopsin (ChR2) enabled monitoring of action potential propagation in ChR2 hearts without any perturbation of electrical dynamics. Importantly, by employing spatially localized optogenetic manipulation, action potential dynamics can be assessed in an all-optical fashion with no artifact related to optical cross-talk between the reporter and actuator. These new environmentally sensitive chromene-based chromophores are also likely to have applications outside voltage imaging.

Aggregates of Cyanine Dyes: When Molecular Vibrations and Electrostatic Screening Make the Difference

Aggregates of cyanine dyes are currently investigated as promising materials for advanced electronic and photonic applications. The spectral properties of aggregates of cyanine dyes can be tuned by altering the supramolecular packing, which is affected by the length of the dye, the presence of alkyl chains, or the nature of the counterions. In this work, we present a joint experimental and theoretical study of a family of cyanine dyes forming aggregates of different types according to the length of the polymethinic chain. Linear and nonlinear optical spectra of aggregates are rationalized here in terms of an essential-state model accounting for intermolecular interactions together with the molecular polarizability and vibronic coupling. A strategy is implemented to properly account for screening effects, distinguishing between electrostatic intermolecular interactions relevant to the ground state (mean-field effect) and the interactions relevant to the excited states (excitonic effects). To the best of our knowledge, this is the first attempt to simulate nonlinear spectral properties of aggregates of symmetric dyes accounting for molecular vibrations.

Spin-Vibronic Control of Intersystem Crossing in Iodine-Substituted Heptamethine Cyanines

Spin-orbit coupling between electronic states of different multiplicity can be strongly coupled to molecular vibrations, and this interaction is becoming recognized as an important mechanism for controlling the course of photochemical reactions. Here, we show that the involvement of spin-vibronic coupling is essential for understanding the photophysics and photochemistry of heptamethine cyanines (Cy7), bearing iodine as a heavy atom in the C3' position of the chain and/or a 3H-indolium core, as potential triplet sensitizers and singlet oxygen producers in methanol and aqueous solutions. The sensitization efficiency was found to be an order of magnitude higher for the chain-substituted than the 3H-indolium core-substituted derivatives. Our ab initio calculations demonstrate that while all optimal structures of Cy7 are characterized by negligible spin-orbit coupling (tenths of cm^-1) with no dependence on the position of the substituent, molecular vibrations lead to its significant increase (tens of cm^-1 for the chain-substituted cyanines), which allowed us to interpret the observed position dependence.

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

Probing DNA structural heterogeneity by identifying conformational subensembles of a bicovalently bound cyanine dye

DNA is a re-configurable, biological information-storage unit, and much remains to be learned about its heterogeneous structural dynamics. For example, while it is known that molecular dyes templated onto DNA exhibit increased photostability, the mechanism by which the structural dynamics of DNA affect the dye photophysics remains unknown. Here, we use femtosecond, two-dimensional electronic spectroscopy measurements of a cyanine dye, Cy5, to probe local conformations in samples of single-stranded DNA (ssDNA-Cy5), double-stranded DNA (dsDNA-Cy5), and Holliday junction DNA (HJ-DNA-Cy5). A line shape analysis of the 2D spectra reveals a strong excitation-emission correlation present in only the dsDNA-Cy5 complex, which is a signature of inhomogeneous broadening. Molecular dynamics simulations support the conclusion that this inhomogeneous broadening arises from a nearly degenerate conformer found only in the dsDNA-Cy5 complex. These insights will support future studies on DNA's structural heterogeneity.

Double Hybrid Density Functionals for the Electronic Excitation Energies of Linear Cyanines

The lowest bright electronic excitations of seven model linear cyanines (CN3-CN15) using 28 double-hybrid (DH) density functionals are benchmarked against accurate and recent CC3 results. Some of these DH functionals are recently designed specifically for excited electronic state calculations. In addition, CIS, CIS(D), SCS-CIS(D), and SOS-CIS(D) were also tested. Four different basis sets were used for the vertical electronic excitation calculations: cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis. Augmented basis sets (e.g. aug-cc-pVDZ and aug-cc-pVTZ) are found to be required for accurate and consistent results using DH functionals. The DH functionals tested in this work are classified into four main groups: global double-hybrids (GDH), range-separated double-hybrids (RSDH), spin-component and spin-opposite scaling global double-hybrids (SCS/SOS-GDH), and spin-component and spin-opposite scaling range-separated double-hybrids (SCS/SOS-RSDH). Within these groups, the SCS/SOS-RSDH group of functionals is found to provide the lowest mean absolute error (MAE) values (in the range 0.020-0.148 eV) in comparison with the GDH group (0.195-0.441 eV), the RSDH group (0.186-0.511 eV), and the SCS/SOS-GDH group (0.079-0.461 eV). Of all the DH functionals and ab initio methods investigated in the present contribution, the following functionals are found to be the most accurate and consistent: SCS-ωB2GPPLYP (MAE = 0.036 eV), SOS-ωB2GPPLYP (MAE = 0.020 eV), SOS-ωB88PP86 (MAE = 0.035 eV), and SOS-ωPBEPP86 (MAE = 0.037 eV). In general, the ab initio methods tested here show mediocre performance as compared to many DH functionals.

Al(III) and Ga(III) Bisphenolate Azadipyrromethene-Based “N2O2” Complexes as Efficient NIR-Fluorophores

Aza-boron-dipyrromethenes (Aza-BODIPYs) are an increasingly studied class of fluorophores. They can be seen as an azadipyrromethene ("aza-DIPY") ligand rigidified by a metalloid, a boron atom. Based on this idea, a series of complexes of group 13 metals (aluminum and gallium) have been synthesized and characterized. The impact of the metal and of the nature of the substituents of aza-DIPY core were investigated. The photophysical and electrochemical properties were determined, and an X-ray structure of an azaGaDIPY was obtained. These data reveal that azaGaDIPY and azaAlDIPY exhibit significant red-shifted fluorescence compared to their analogue aza-BODIPY. Their emission can go up to 800 nm for the maximum emission length and up to NIR-II for the emission tail. This, associated with their electrochemical stability (no metal release whether oxidized or reduced) makes them a promising class of fluorophores for optical medical imaging. Moreover, X-ray structure and molecular modeling studies have shown that this redshift seems to be more due to the geometry around the boron/metal than to the nature of the metal.

Accurate Vertical Excitation Energies of BODIPY/Aza-BODIPY Derivatives from Excited-State Mean-Field Calculations

We report a benchmark study of vertical excitation energies and oscillator strengths for the HOMO → LUMO transitions of 17 boron–dipyrromethene (BODIPY) structures, showing a large variety of ring sizes and substituents. Results obtained at the time-dependent density functional theory (TDDFT) and at the delta-self-consistent-field (ΔSCF) by using 13 different exchange correlation kernels (within LDA, GGA, hybrid, and range-separated approximations) are benchmarked against the experimental excitation energies when available. It is found that the time-independent ΔSCF DFT method, when used in combination with hybrid PBE0 and B3LYP functionals, largely outperforms TDDFT and can be quite competitive, in terms of accuracy, with computationally more costly wave function based methods such as CC2 and CASPT2.

Rational Design of a Near-Infrared Ratiometric Probe with a Large Stokes Shift: Visualization of Polarity Abnormalities in Non-Alcoholic Fatty Liver Model Mice

Tracking liver polarity with noninvasive and dynamic imaging techniques is helpful to better understand the non-alcoholic fatty liver (NAFL). Herein, a novel near-infrared (NIR) fluorescent probe Cy-Mp is constructed using a "symmetry collapse" strategy. The structure modification leads to the conversion of locally excited state fluorescence to charge transfer state fluorescence. Cy-Mp emits at near-infrared (NIR) wavelengths with high photostability as well as a large Stokes shift. Cy-Mp exhibits a ratiometric response to polarity, providing more accurate analysis of intracellular polarity via the built-in internal reference correction. Most importantly, the in vivo studies indicate that Cy-Mp can accumulate in the liver and the decreased polarity in the liver of mice with NAFL is verified by the ratiometric imaging, implying the great potential of Cy-Mp in the diagnosis of NAFL.

An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules

The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.

Modeling the Electronic Absorption Spectra of the Indocarbocyanine Cy3

Accurate modeling of optical spectra requires careful treatment of the molecular structures and vibronic, environmental, and thermal contributions. The accuracy of the computational methods used to simulate absorption spectra is limited by their ability to account for all the factors that affect the spectral shapes and energetics. The ensemble-based approaches are widely used to model the absorption spectra of molecules in the condensed-phase, and their performance is system dependent. The Franck–Condon approach is suitable for simulating high resolution spectra of rigid systems, and its accuracy is limited mainly by the harmonic approximation. In this work, the absorption spectrum of the widely used cyanine Cy3 is simulated using the ensemble approach via classical and quantum sampling, as well as, the Franck–Condon approach. The factors limiting the ensemble approaches, including the sampling and force field effects, are tested, while the vertical and adiabatic harmonic approximations of the Franck–Condon approach are also systematically examined. Our results show that all the vertical methods, including the ensemble approach, are not suitable to model the absorption spectrum of Cy3, and recommend the adiabatic methods as suitable approaches for the modeling of spectra with strong vibronic contributions. We find that the thermal effects, the low frequency modes, and the simultaneous vibrational excitations have prominent contributions to the Cy3 spectrum. The inclusion of the solvent stabilizes the energetics significantly, while its negligible effect on the spectral shapes aligns well with the experimental observations.

Sensitization of Nd3+ Luminescence by Simultaneous Two-Photon Excitation through a Coordinating Polymethinic Antenna

We have designed and synthesized two new cyaninic Nd³⁺ complexes where the lanthanide emission can be induced from simultaneous two-photon absorption followed by energy migration. These complexes correspond to a molecular design that uses an antenna ligand formed by the functionalization of a heptamethine dye with 5-ol-phenanthroline or 4-phenyl-terpyridine derivatives. These complexes employ the important nonlinear optical properties of symmetric polymethines to sensitize the lanthanide ion. We verified that simultaneous biphotonic excitation indirectly induces the ⁴F_3/2 → ⁴I_11/2 Nd³⁺ emission using femtosecond laser pulses tuned below the first electronic transition of the antenna. The simultaneous two-photon excitation events initially form the nonlinear-active second excited singlet of the polymethine antenna, which rapidly evolves into its first excited singlet. This state in turn induces the formation of the emissive Nd³⁺ states through energy transfer. The role of the first excited singlet of the antenna as the donor state in this process was verified through time resolution of the antenna's fluorescence. These measurements also provided the rates for antenna-lanthanide energy transfer, which indicate that the phenanthroline-type ligand is approximately five times more efficient for energy transfer than the phenyl-terpyridine derivative due to their relative donor-acceptor distances. The simultaneous two-photon excitation of this polymethine antenna allows for high spatial localization of the Nd³⁺excitation events.

Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies

Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events but also to provide reference values for vertical transition energies, hence allowing benchmarking of lower-order models. In this category, both the complete-active-space second-order perturbation theory (CASPT2) and the N-electron valence state second-order perturbation theory (NEVPT2) are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches, the present work relies on highly accurate (±0.03 eV) aug-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 n → π*, 119 π → π*, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially contracted (PC) and strongly contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of 0.11 and 0.13 eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focused on single-reference wave function methods ( J. Chem. Theory Comput. 2018, 14, 4360; J. Chem. Theory Comput. 2020, 16, 1711), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.

Femtosecond Spectroscopy and Nonlinear Optical Properties of aza-BODIPY Derivatives in Solution

The fast relaxation processes in the excited electronic states of functionalized aza-boron-dipyrromethene (aza-BODIPY) derivatives (1-4) were investigated in liquid media at room temperature, including the linear photophysical, photochemical, and nonlinear optical (NLO) properties. Optical gain was revealed for nonfluorescent derivatives 3 and 4 in the near infrared (NIR) spectral range under femtosecond excitation. The values of two-photon absorption (2PA) and excited-state absorption (ESA) cross-sections were obtained for 1-4 in dichloromethane using femtosecond Z-scans, and the role of bromine substituents in the molecular structures of 2 and 4 is discussed. The nature of the excited states involved in electronic transitions of these dyes was investigated using quantum-chemical TD-DFT calculations, and the obtained spectral parameters are in reasonable agreement with the experimental data. Significant 2PA (maxima cross-sections ∼2000 GM), and large ESA cross-sections ∼10^-20  m² of these new aza-BODIPY derivatives 1-4 along with their measured high photostability reveal their potential for photonic applications in general and optical limiting in particular.

Human Serum Albumin Dimerization Enhances the S2 Emission of Bound Cyanine IR806

Cyanine molecules are important phototheranostic compounds given their high fluorescence yield in the near-infrared region of the spectrum. We report on the frequency and time-resolved spectroscopy of the S₂ state of IR806, which demonstrates enhanced emission upon binding to the hydrophobic pocket of human serum albumin (HSA). From excitation-emission matrix spectra and electronic structure calculations, we identify the emission as one associated with a state having the polymethine chain twisted out of plane by 103°. In addition, we find that this configuration is significantly stabilized as the concentration of HSA increases. Spectroscopic changes associated with the S₁ and S₂ states of IR806 as a function of HSA concentration, as well as anisotropy measurements, confirm the formation of HSA dimers at concentrations greater than 10 μM. These findings imply that the longer-lived S₂ state configuration can lead to more efficient phototherapy agents, and cyanine S₂ spectroscopy may be a useful tool to determine the oligomerization state of HSA.

Direct transformation of coumarins into orange-red emitting rhodols

The lactone carbonyl group of coumarin derivatives has been shown to participate in intramolecular Knoevenagel condensations, enabling the unprecedented direct transformation of coumarins into rhodols. The resulting rhodols, possessing two ester groups, have very intense orange-red fluorescence.

Choosing Bad versus Worse: Predictions of Two-Photon-Absorption Strengths Based on Popular Density Functional Approximations

We present a benchmark study of density functional approximation (DFA) performances in predicting the two-photon-absorption strengths in π-conjugated molecules containing electron-donating/-accepting moieties. A set of 48 organic molecules is chosen for this purpose, for which the two-photon-absorption (2PA) parameters are evaluated using different DFAs, including BLYP, PBE, B3LYP, PBE0, CAM-B3LYP, LC-BLYP, and optimally tuned LC-BLYP. Minnesota functionals and ωB97X-D are also used, applying the two-state approximation, for a subset of molecules. The efficient resolution-of-identity implementation of the coupled-cluster CC2 model (RI-CC2) is used as a reference for the assessment of the DFAs. Two-state models within the framework of both DFAs and RI-CC2 are used to gain a deeper insight into the performance of different DFAs. Our results give a clear picture of the performance of the density functionals in describing the two-photon activity in dipolar π-conjugated systems. The results show that global hybrids are best suited to reproduce the absolute values of 2PA strengths of donor–acceptor molecules. The range-separated functionals CAM-B3LYP and optimally tuned LC-BLYP, however, show the highest linear correlations with the reference RI-CC2 results. Hence, we recommend the latter DFAs for structure–property studies across large series of dipolar compounds.

Reference Excitation Energies of Increasingly Large Molecules: A QMC Study of Cyanine Dyes

We revisit here the lowest vertical excitations of cyanine dyes using quantum Monte Carlo and leverage recent developments to systematically improve on previous results. In particular, we employ a protocol for the construction of compact and accurate multideterminant Jastrow-Slater wave functions for multiple states, which we have recently validated on the excited-state properties of several small prototypical molecules. Here, we obtain quantum Monte Carlo excitation energies in excellent agreement with high-level coupled cluster for all the cyanines where the coupled cluster method is applicable. Furthermore, we push our protocol to longer chains, demonstrating that quantum Monte Carlo is a viable methodology to establish reference data at system sizes which are hard to reach with other high-end approaches of similar accuracy. Finally, we determine which ingredients are key to an accurate treatment of these challenging systems and rationalize why a description of the excitation based on only active π orbitals lacks the desired accuracy for the shorter chains.

Accurate predictions of the electronic excited states of BODIPY based dye sensitizers using spin-component-scaled double-hybrid functionals: a TD-DFT benchmark study

The vertical excitation energies of 13 BODIPY based dye sensitizers are benchmarked by means of TD-DFT, using 36 functionals from different DFT rungs. Most TD-DFT results were found to overestimate the excitation energies, and show mean absolute error (MAE) values in the range 0.2–0.5 eV. The dispersion-corrected, spin-component-scaled, double-hybrid (DSD) functionals DSD-BLYP and DSD-PBEP86 were found to have the smallest MAE values of 0.083 eV and 0.106 eV, respectively, which is close to the range of average errors found in the more expensive coupled-cluster methods. Moreover, DSD-BLYP and DSD-PBEP86 functionals show excellent consistency and quality of results (standard deviation = 0.048 eV and 0.069 eV respectively). However, the range separated hybrid (RSH) and the range separated double hybrid (RSDH) functionals were found to provide the best predictability (linear determination coefficient R² > 0.97 eV).

The excitation energies of 13 BODIPY dye sensitizers are benchmarked by means of TD-DFT, using 36 functionals. Spin-component-scaled double-hybrid (DSD) functionals are found to show the best performance.

Give or Take: Effects of Electron-Accepting/-Withdrawing Groups in Red-Fluorescent BODIPY Molecular Rotors

Mapping microviscosity, temperature, and polarity in biosystems is an important capability that can aid in disease detection. This can be achieved using fluorescent sensors based on a green-emitting BODIPY group. However, red fluorescent sensors are desired for convenient imaging of biological samples. It is known that phenyl substituents in the β position of the BODIPY core can shift the fluorescence spectra to longer wavelengths. In this research, we report how electron-withdrawing (EWG) and -donating (EDG) groups can change the spectral and sensory properties of β-phenyl-substituted BODIPYs. We present a trifluoromethyl-substituted (EWG) conjugate with moderate temperature sensing properties and a methoxy-substituted (EDG) molecule that could be used as a lifetime-based polarity probe. In this study, we utilise experimental results of steady-state and time-resolved fluorescence, as well as quantum chemical calculations using density functional theory (DFT). We also explain how the energy barrier height (Ea) for non-radiative relaxation affects the probe’s sensitivity to temperature and viscosity and provide appropriate Ea ranges for the best possible sensitivity to viscosity and temperature.

Going beyond the borders: pyrrolo[3,2-b]pyrroles with deep red emission

A two-step route to strongly absorbing and efficiently orange to deep red fluorescent, doubly B/N-doped, ladder-type pyrrolo[3,2-b]pyrroles has been developed. We synthesize and study a series of derivatives of these four-coordinate boron-containing, nominally quadrupolar materials, which mostly exhibit one-photon absorption in the 500-600 nm range with the peak molar extinction coefficients reaching 150 000, and emission in the 520-670 nm range with the fluorescence quantum yields reaching 0.90. Within the family of these ultrastable dyes even small structural changes lead to significant variations of the photophysical properties, in some cases attributed to reversal of energy ordering of alternate-parity excited electronic states. Effective preservation of ground-state inversion symmetry was evidenced by very weak two-photon absorption (2PA) at excitation wavelengths corresponding to the lowest-energy, strongly one-photon allowed purely electronic transition. π-Expanded derivatives and those possessing electron-donating groups showed the most red-shifted absorption- and emission spectra, while displaying remarkably high peak 2PA cross-section (σ _2PA) values reaching ∼2400 GM at around 760 nm, corresponding to a two-photon allowed higher-energy excited state. At the same time, derivatives lacking π-expansion were found to have a relatively weak 2PA peak centered at ca. 800-900 nm with the maximum σ _2PA ∼50-250 GM. Our findings are augmented by theoretical calculations performed using TD-DFT method, which reproduce the main experimental trends, including the 2PA, in a nearly quantitative manner. Electrochemical studies revealed that the HOMO of the new dyes is located at ca. -5.35 eV making them relatively electron rich in spite of the presence of two B^--N⁺ dative bonds. These dyes undergo a fully reversible first oxidation, located on the diphenylpyrrolo[3,2-b]pyrrole core, directly to the di(radical cation) stage.

A Peierls Transition in Long Polymethine Molecular Wires: Evolution of Molecular Geometry and Single-Molecule Conductance

Molecules capable of mediating charge transport over several nanometers with minimal decay in conductance have fundamental and technological implications. Polymethine cyanine dyes are fascinating molecular wires because up to a critical length, they have no bond-length alternation (BLA) and their electronic structure resembles a one-dimensional free-electron gas. Beyond this threshold, they undergo a symmetry-breaking Peierls transition, which increases the HOMO-LUMO gap. We have investigated cationic cyanines with central polymethine chains of 5-13 carbon atoms (Cy3⁺-Cy11⁺). The absorption spectra and crystal structures show that symmetry breaking is sensitive to the polarity of the medium and the size of the counterion. X-ray crystallography reveals that Cy9·PF₆ and Cy11·B(C₆F₅)₄ are Peierls distorted, with high BLA at one end of the π-system, away from the partially delocalized positive charge. This pattern of BLA distribution resembles that of solitons in polyacetylene. The single-molecule conductance is essentially independent of molecular length for the polymethine salts of Cy3⁺-Cy11⁺ with the large B(C₆F₅)₄^- counterion, but with the PF₆^- counterion, the conductance decreases for the longer molecules, Cy7⁺-Cy11⁺, because this smaller anion polarizes the π-system, inducing a symmetry-breaking transition. At higher bias (0.9 V), the conductance of the shorter chains, Cy3⁺-Cy7⁺, increases with length (negative attenuation factor, β = -1.6 nm^-1), but the conductance still drops in Cy9⁺ and Cy11⁺ with the small polarizing PF₆^- counteranion.

Radiative lifetime of a BODIPY dye as calculated by TDDFT and EOM-CCSD methods: solvent and vibronic effects

The radiative emission lifetime and associated S₁ excited state properties of a BODIPY dye are investigated with TDDFT and EOM-CCSD calculations. The effects of a solvent are described with the polarizable continuum model using the linear response (LR) approach as well as state-specific methods. The Franck-Condon (FC), Herzberg-Teller (HT) and Duschinsky vibronic effects are evaluated for the absorption and emission spectra, and for the radiative lifetime. The transition energies, spectra shapes and radiative lifetime are assessed with respect to experimental results. It is found that the TDDFT transition energies are overestimated by about 0.4-0.5 eV, whereas EOM-CCSD improves the vertical emission energy by about 0.1 eV in comparison to TDDFT. The solvatochromic and Stokes shifts are better reproduced by the state-specific solvation methods, which show that these methods are more suited than the LR model to describe the solvent effects on the BODIPY dye. The vibronic effects lead to an increase of the radiative lifetime of about 0.4 to 1.0 ns depending on the theoretical approach, which highlights the importance of such effects. Moreover, the HT effects are negligible on both the spectra and lifetime, which demonstrates that the FC approximation is accurate for the BODIPY dye. Finally, the comparison with experimental data shows that the radiative lifetimes predicted by EOM-CCSD and TDDFT have comparable accuracy.

Linear and Nonlinear Optical Processes Controlling S2 and S1 Dual Fluorescence in Cyanine Dyes

We report on the changes in the dual fluorescence of two cyanine dyes IR144 and IR140 as a function of viscosity and probe their internal conversion dynamics from S₂ to S₁ via their dependence on a femtosecond laser pulse chirp. Steady-state and time-resolved measurements performed in methanol, ethanol, propanol, ethylene glycol, and glycerol solutions are presented. Quantum calculations reveal the presence of three excited states responsible for the experimental observations. Above the first excited state, we find an excited state, which we designate as S₁', that relaxes to the S₁ minimum, and we find that the S₂ state has two stable configurations. Chirp-dependence measurements, aided by numerical simulations, reveal how internal conversion from S₂ to S₁ depends on solvent viscosity and pulse duration. By combining solvent viscosity, transform-limited pulses, and chirped pulses, we obtain an overall change in the S₂/S₁ population ratio of a factor of 86 and 55 for IR144 and IR140, respectively. The increase in the S₂/S₁ ratio is explained by a two-photon transition to a higher excited state. The ability to maximize the population of higher excited states by delaying or bypassing nonradiative relaxation may lead to the increased efficiency of photochemical processes.

Assessment of hybrid functionals for singlet and triplet excitations: Why do some local hybrid functionals perform so well for triplet excitation energies?

The performance of various hybrid density functionals is assessed for 105 singlet and 105 corresponding triplet vertical excitation energies from the QUEST database. The overall lowest mean absolute error is obtained with the local hybrid (LH) functional LH12ct-SsirPW92 with individual errors of 0.11 eV (0.11 eV) for singlet (triplet) n → π* excitations and 0.29 eV (0.17 eV) for π → π* excitations. This is slightly better than with the overall best performing global hybrid M06-2X [n → π*: 0.13 eV (0.17 eV), π → π*: 0.30 eV (0.20 eV)], while most other global and range-separated hybrids and some LHs suffer from the "triplet problem" of time-dependent density functional theory. This is exemplified by correlating the errors for singlet and triplet excitations on a state-by-state basis. The excellent performance of LHs based on a common local mixing function, i.e., an LMF constructed from the spin-summed rather than the spin-resolved semilocal quantities, is systematically investigated by the introduction of a spin-channel interpolation scheme that allows us to continuously modulate the fraction of opposite-spin terms used in the LMF. The correlation of triplet and singlet errors is systematically improved for the n → π* excitations when larger fractions of the opposite-spin-channel are used in the LMF, whereas this effect is limited for the π → π* excitations. This strongly supports a previously made hypothesis that attributes the excellent performance of LHs based on a common LMF to cross-spin-channel nondynamical correlation terms.

Designing a Red-Emitting Viscosity-Sensitive BODIPY Fluorophore for Intracellular Viscosity Imaging

Viscosity imaging at a microscopic scale can provide important information about biosystems, including the development of serious illnesses. Microviscosity imaging is achievable with viscosity-sensitive fluorophores, the most popular of which are based on the BODIPY group. However, most of the BODIPY probes fluoresce green light, whereas the red luminescence is desired for the imaging of biological samples. Designing a new viscosity probe with suitable spectroscopic properties is a challenging task because it is difficult to preserve viscosity sensitivity after modifying the molecular structure. Here we describe how we developed a new red-emitting, viscosity-sensitive, BODIPY fluorophore BP-PH-2M-NO₂ that is suitable for reliable intracellular viscosity imaging of lipid droplets in MCF-7 breast cancer cells. The design of BP-PH-2M-NO₂ was aided by DFT calculations that allowed a successful prediction of the viscosity sensitivity of fluorophores before synthesis. In summary, we report a new red viscosity probe possessing monoexponential fluorescence decay that makes it attractive for lifetime-based viscosity imaging.

Accurate Modeling of Excitonic Coupling in Cyanine Dye Cy3

Accurate modeling of excitonic coupling in molecules is of great importance for inferring the structures and dynamics of coupled systems. Cy3 is a cyanine dye that is widely used in molecular spectroscopy. Its well-separated excitation bands, high sensitivity to the surroundings, and the high energy transfer efficiency make it a perfect choice for excitonic coupling experiments. Many methods have been used to model the excitonic coupling in molecules with varying degrees of accuracy. The atomic transition charge model offers a high-accuracy and cost-effective way to calculating the excitonic coupling. The main focus of this work is to generate high-quality atomic transition charges that can accurately model the Cy3 dye's transition density. The transition density of the excitation of the ground to first excited state is calculated using configuration-interaction singles and time-dependent density functional theory and is benchmarked against the algebraic diagrammatic construction method. Using the transition density we derived the atomic transition charges using two approaches: Mulliken population analysis and charges fitted to the transition electrostatic potential. The quality of the charges is examined, and their ability to accurately calculate the excitonic coupling is assessed via comparison to experimental data of an artificial biscyanine construct. Theoretical comparisons to the supermolecule ab initio couplings and the widely used point-dipole approximation are also made. Results show that using the transition electrostatic potential is a reliable approach for generating the transition atomic charges. A high-quality set of charges, that can be used to model the Cy3 dye dimer excitonic coupling with high-accuracy and a reasonable computational cost, is obtained.

Mild Open-Shell Character of BODIPY and Its Impact on Singlet and Triplet Excitation Energies

This study describes and rationalizes the electronic structure of BODIPY combining a large variety of quantum chemistry methods and computational tools. Examination of the obtained results using state-of-the-art electronic structure analyses provides a new and complete interpretation of the nature of low-lying electronic states in BODIPY and elucidates the limitations of excited-state methods in the computation of T₁ and S₁ energies, that is, systematic under- and overestimation of time-dependent density functional theory energies, respectively, and a large overestimation of the T₁/S₁ energy gap. Our analysis identifies the important role and physical origin of the mild open-shell character in the BODIPY ground state, that is, strong highest occupied and lowest unoccupied molecular orbital exchange interactions. The study provides guidelines for the accurate quantification of the T₁/S₁ gap, which is extremely relevant for the computational investigation of the photophysical properties of BODIPY and its derivatives. These conclusions should be taken into consideration in order to predict and interpret conspicuous photoactivated phenomena such as intersystem crossing, singlet fission, and triplet-triplet annihilation. Moreover, we believe that our study might provide new ideas and strategies for the analysis of other molecular chromophores.