The geometries, absorption and fluorescence wavelengths of solvated fluorescent coumarins: A CIS and TD-DFT comparative study

Beyond Explored Functionals: A Computational Journey of Two-Photon Absorption

We present a thorough investigation into the efficacy of 19 density functional theory (DFT) functionals, relative to RI-CC2 results, for computing two-photon absorption (2PA) cross sections (σ2PA) and key dipole moments (|μ00|, |μ11|, |Δμ|, |μ01|) for a series of coumarin dyes in the gas-phase. The functionals include different categories, including local density approximation (LDA), generalized gradient approximation (GGA), hybrid-GGA (H-GGA), range-separated hybrid-GGA (RSH-GGA), meta-GGA (M-GGA), and hybrid M-GGA (HM-GGA), with 14 of them being subjected to analysis for the first time with respect to predicting σ2PA values. Analysis reveals that functionals integrating both short-range (SR) and long-range (LR) corrections, particularly those within the RSH-GGA and HM-GGA classes, outperform the others. Furthermore, the range-separation approach was found more impactful compared to the varying percentages of Hartree-Fock exchange (HF Ex) within different functionals. The functionals traditionally recommended for 2PA do not appear among the top 9 in our study, which is particularly interesting, as these top-performing functionals have not been previously investigated in this context. This list is dominated by M11, QTP variants, ωB97X, ωB97X-V, and M06-2X, surpassing the performance of other functionals, including the commonly used CAM-B3LYP.

Degenerate and non-degenerate two-photon absorption of coumarin dyes

Two-photon absorption (2PA) spectroscopy is a robust bioimaging tool that depends on the determined cross-sections (σ^2PA). The absorption of both photons occurs simultaneously with equivalent (degenerate) or different (non-degenerate) photon energies, D-2PA and ND-2PA, respectively. The former has been investigated experimentally and computationally for many systems, while the latter remains relatively unexplored computationally and limited experimentally. In this study, response theory using time-dependent density functional theory (TD-DFT) and the 2-state model (2SM) have been utilized to investigate σ^D-2PA and σ^ND-2PA for the excitation to the lowest energy singlet state (S₁) of coumarin, coumarin 6, coumarin 120, coumarin 307, and coumarin 343. Solvents involved were methanol (MeOH), chloroform (ClForm), and dimethylsulfoxide (DMSO), where the latter leads to the largest σ^2PA. Values of σ^2PA are largest for coumarin 6 and lowest for coumarin, which illustrates the effect of substituents. The 2SM clarifies how the largest cross-sections correspond to molecules with the largest transition dipole moments, μ₀₁. In general, σ^D-2SM computations agree with σ^D-2PA. Moreover, σ^ND-2SM are in qualitative agreement with σ^ND-2PA with comparable enhancement relative to σ^D-2PA. Overall, σ^ND-2PA are larger than σ^D-2PA where the increase is in the range of 22% to 49%, depending on the coumarin as well as the relative energies of the two photons. This work aids in future investigations into various fluorophores to understand their photophysical properties for ND-2PA.

The Curious Case of Aqueous Warfarin: Structural Isomers or Distinct Excited States?

Warfarin is a potent anti-coagulant drug and is on the World Health Organization's List of Essential Medicines. Additionally, it displays fluorescence enhancement upon binding to human serum albumin, making warfarin a prototype fluorescent probe in biology. Despite its biological significance, the current structural assignment of warfarin in aqueous solution is based on indirect evidence in organic solvents. Warfarin is known to exist in different isomeric forms-open-chain, hemiketal, and anionic forms-based on the solvent and pH. Moreover, warfarin displays a dual absorption feature in several solvents, which has been employed to study the ring-chain isomerism between its open-chain and hemiketal isomers. In this study, our pH-dependent experiments on warfarin and structurally constrained warfarin derivatives in aqueous solution demonstrate that the structural assignment of warfarin solely on the basis of its absorption spectrum is erroneous. Using a combination of steady-state and time-resolved spectroscopic experiments, along with quantum chemical calculations, we assign the observed dual absorption to two distinct π → π* transitions in the 4-hydroxycoumarin moiety of warfarin. Furthermore, we unambiguously identify the isomeric form of warfarin that binds to human serum albumin in aqueous buffer.

Modeling the Electron Transfer Chain in an Artificial Photosynthetic Machine

The development of efficient artificial leaves relies on the subtle combination of molecular assemblies able to absorb sunlight, converting light energy into electrochemical potential energy and finally transducing it into accessible chemical energy. The electronic design of these charge transfer molecular machines is crucial to build a complex supramolecular architecture for the light energy conversion. Here, we present an ab initio simulation of the whole decay pathways of a recently proposed artificial molecular reaction center. A complete structural and energetic characterization has been carried out with methods based on density functional theory, its time-dependent version, and a broken-symmetry approach. On the basis of our findings we provide a revision of the pathway only indirectly postulated from an experimental point of view, along with unprecedented and significant insights on the electronic and nuclear structure of intramolecular charge-separated states, which are fundamental for the application of this molecular assembly in photoelectrochemical cells. Importantly, we unravel the molecular driving forces of the various charge transfer steps, in particular those leading to the proton-coupled electron transfer final product, highlighting key elements for the future design strategies of such molecular assays.

Evaluation and correction on quinones' quantification errors: Derived from the coexistence of different quinone species and pH-sensitive feature

Quinones are becoming an essential tool for refractory organics treatment, while their quantification may be not well-considered. In this paper, two kinds of potential errors in quantification were evaluated in multiple pH conditions. They were derived from the coexistence of oxidized/reduced quinone species (Type I) and pH-sensitive feature (Type II), respectively. These errors would remarkably influence the accuracy of quantification while they haven't been emphasized. Thus, to elaborate the relationship between the two types of errors and the absorbance or pH conditions, three typical quinones [Anthraquinone-1-sulfonate (α-AQS), anthraquinone-2,6-disulfonate (AQDS) and lawsone] were selected and their acid dissociation coefficients (pKa) as well as UV-Vis spectra were determined. Results revealed that, for Type I, the relative error (RE) of α-AQS concentration would exceed the limit (5%) when reduced α-AQS was below 48% of total α-AQS. Similar results were found for lawsone. However, the RE can be eliminated by the equation established in this paper. For Type II, the pH-sensitive feature was related to the pKa values of quinones. Absorbances of α-AQS and lawsone would change remarkably with pH variation. Therefore, a model for correction was established. Analog data showed high consistency with experimental data [r = 0.995 (n = 25, p < 0.01) and r = 0.997 (n = 36, p < 0.01), for lawsone and α-AQS respectively]. Especially, the determination of AQDS concentrations was noticed to be pH-independent at 437 nm under pH 4.00 to 9.18 conditions. Based on these features, a comprehensive data solution was proposed for handling these errors.

Photophysics of Derivatives of 3-Hydroxybenzo[c]coumarin

The photophysical studies of two phenols, derivatives of 3-hydroxybenzo[c]coumarin, were performed in n-nonane matrix at 5 K. Unstructured fluorescence spectrum of the derivative bearing a salicylaldehyde moiety, whose onset is shifted by ca. 3000 cm(-1) to lower energy in respect to that of absorption, and short decay time of this emission (0.75 ns) suggested the occurrence of excited-state intramolecular proton transfer (ESIPT). The experimental results were interpreted with the aid of quantum chemistry calculations performed with the DFT and TDDFT/B3LYP/6-31++G(d,p) methods.

Dye chemistry with time-dependent density functional theory

In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.

"Turn-on" fluorescence probe integrated polymer nanoparticles for sensing biological thiol molecules

A "turn-on" thiol-responsive fluorescence probe was synthesized and integrated into polymeric nanoparticles for sensing intracellular thiols. There is a photo-induced electron transfer process in the off state of the probe, and this process is terminated upon the reaction with thiol compounds. Configuration interaction singles (CIS) calculation was performed to confirm the mechanism of this process. A series of sensing studies were carried out, showing that the probe-integrated nanoparticles were highly selective towards biological thiol compounds over non-thiolated amino acids. Kinetic studies were also performed to investigate the relative reaction rate between the probe and the thiolated amino acids. Subsequently, the Gibbs free energy of the reactions was explored by means of the electrochemical method. Finally, the detection system was employed for sensing intracellular thiols in cancer cells, and the sensing selectivity could be further enhanced with the use of a cancer cell-targeting ligand in the nanoparticles. This development paves a path for the sensing and detection of biological thiols, serving as a potential diagnostic tool in the future.

Spectral properties of coumarin derivatives in various environments

The structural and spectral properties of coumarin derivatives in complex environments were investigated within the time-dependent density functional theory (TD DFT). Absorption spectra calculations were obtained at TD PBE0/6-31+G(d,p) level of theory for coumarin47 in the gas-phase and in various polar and non-polar organic solvents. The geometries of coumarins 6, 30, 47 and 522 in the gas phase and in inclusion complexes with the β-cyclodextrin (βCD) were determined by PM3 and DFT (HCTH/6-31G) calculations. Encapsulation of coumarin in βCD and associated changes in electronic structure produced either a red or blue shift in the absorption spectra of coumarins. A proposed cavity model for βCD-coumarin complex in water solution allowed identification of various contributions to the overall shift in the absorption spectra of coumarin upon complex formation in a solvent environment

TD-DFT Assessment of Functionals for Optical 0-0 Transitions in Solvated Dyes

Using TD-DFT, we performed simulations of the adiabatic energies of 40 fluorescent molecules for which the experimental 0-0 energies in condensed phase are available. We used six hybrid functionals (B3LYP, PBE0, M06, M06-2X, CAM-B3LYP, and LC-PBE) that have been shown to provide accurate transition energies in previous TD-DFT assessments, selected two diffuse-containing basis sets, and applied the most recent models for estimating bulk solvation effects. In each case, the correction arising from the difference of zero-point vibrational energies between the ground and the excited states has been consistently determined. Basis set effects have also been carefully studied. It turned out that PBE0 and M06 are the most effective functionals in terms of average deviation (mean absolute error of 0.22-0.23 eV). However, both the M06-2X global hybrid that contains more exact exchange and the CAM-B3LYP range-separated hybrid significantly improve the consistency of the prediction for a relatively negligible degradation of the average error. In addition, we assessed (1) the cross-structure/spectra relationships, (2) the importance of solvent effects, and (3) the differences between adiabatic and vertical energies.

Experimental and theoretical study of the inclusion complexes of 3-carboxycoumarin acid with β- and 2-hydroxypropyl-β-cyclodextrins

The association process of a host-guest system, cyclodextrins (CyD) - 3-carboxycoumarin acid (3CCA) was followed by means of UV-vis, circular dichroism and steady-state fluorescence spectroscopies in buffer solution at pH=1. The experimental data were analyzed in order to get information on the stoichiometry, the equilibrium constants and the geometry of the inclusion complexes. In the circular dichroism spectra, a positive induced signal was obtained reflecting that the guest penetrates the cavity in such a way that the transition moment of the electronic band is quasi parallel to the host main axis. The experimental data are supported by the DFT and TDDFT (B3LYP/6-31G) calculations performed on the isolated ligand and by TDDFT (ZINDO) calculations carried out on the supramolecular ligand-cyclodextrin system.

Spectroscopic study of solvent effects on the electronic absorption spectra of flavone and 7-hydroxyflavone in neat and binary solvent mixtures

The solvatochromic characteristics of flavone and 7-hydroxyflavone were investigated in neat and binary solvent mixtures. The spectral shifts of these solutes were correlated with the Kamlet and Taft parameters (α, β and π*) using linear solvation energy relationships. The multiparametric analysis indicates that both specific hydrogen bond donor ability and non-specific dipolar interactions of the solvents play an important role in absorption maxima of flavone in pure solvents. The hydrogen bond acceptor ability of the solvent was the main parameter affecting the absorption maxima of 7-hydroxyflavone. The simulated absorption spectra using a TD-DFT method were in good agreement with the experimental ones for both flavones. Index of preferential solvation was calculated as a function of solvent composition. Preferential solvation by ethanol was detected in cyclohexane-ethanol and acetonitrile-ethanol mixtures for flavone and in acetonitrile-ethanol mixtures for 7-hydroxyflavone. These results indicate that intermolecular hydrogen bonds between solute and solvent are responsible for the non-linear variation of the solvatochromic shifts on the mole fraction of ethanol in the analyzed binary mixtures.

Planar vs. twisted intramolecular charge transfer mechanism in Nile Red: new hints from theory

Using a time-dependent density functional theory approach and taking into account bulk solvent effects, we investigate the absorption and fluorescence spectra of Nile Red. In particular, we have assessed both the planar and twisted intramolecular charge transfer mechanism by using a panel of exchange correlation functionals including both global and range-separated hybrids, refined solvent models and the simulation of vibronic couplings. It turned out that the appropriate choice of the functional is of prime importance to obtain, not only quantitatively accurate values, but also qualitatively correct evolution of the spectral features with respect to the dihedral angles of the amino group. At the light of this study, the interpretation of the experimental data is critically re-examined and compared to typical dual-fluorescent molecules.

Accurate simulation of optical properties in dyes

Since Antiquity, humans have produced and commercialized dyes. To this day, extraction of natural dyes often requires lengthy and costly procedures. In the 19th century, global markets and new industrial products drove a significant effort to synthesize artificial dyes, characterized by low production costs, huge quantities, and new optical properties (colors). Dyes that encompass classes of molecules absorbing in the UV-visible part of the electromagnetic spectrum now have a wider range of applications, including coloring (textiles, food, paintings), energy production (photovoltaic cells, OLEDs), or pharmaceuticals (diagnostics, drugs). Parallel to the growth in dye applications, researchers have increased their efforts to design and synthesize new dyes to customize absorption and emission properties. In particular, dyes containing one or more metallic centers allow for the construction of fairly sophisticated systems capable of selectively reacting to light of a given wavelength and behaving as molecular devices (photochemical molecular devices, PMDs).Theoretical tools able to predict and interpret the excited-state properties of organic and inorganic dyes allow for an efficient screening of photochemical centers. In this Account, we report recent developments defining a quantitative ab initio protocol (based on time-dependent density functional theory) for modeling dye spectral properties. In particular, we discuss the importance of several parameters, such as the methods used for electronic structure calculations, solvent effects, and statistical treatments. In addition, we illustrate the performance of such simulation tools through case studies. We also comment on current weak points of these methods and ways to improve them.