Assessment of Functionals for TDDFT Calculations of One- and Two-Photon Absorption Properties of Neutral and Anionic Fluorescent Proteins Chromophores

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.

Toward Accurate Two-Photon Absorption Spectrum Simulations: Exploring the Landscape beyond the Generalized Gradient Approximation

In this Letter, we present a pioneering analysis of the density functional approximations (DFAs) beyond the generalized gradient approximation (GGA) for predicting two-photon absorption (2PA) strengths of a set of push-pull π-conjugated molecules. In more detail, we have employed a variety of meta-generalized gradient approximation (meta-GGA) functionals, including SCAN, MN15, and M06-2X, to assess their accuracy in describing the 2PA properties of a chosen set of 48 organic molecules. Analytic quadratic response theory is employed for these functionals, and their performance is compared against the previously studied DFAs and reference data obtained at the coupled-cluster CC2 level combined with the resolution-of-identity approximation (RI-CC2). A detailed analysis of the meta-GGA functional performance is provided, demonstrating that they improve upon their predecessors in capturing the key electronic features of the π-conjugated two-photon absorbers. In particular, the Minnesota functional MN15 shows very promising results as it delivers pleasingly accurate chemical rankings for two-photon transition strengths and excited-state dipole moments.

Molecular modelling of fullerene C60 functionalized by nitric oxide for use in biological environment

The unique potential of fullerene C60 for various biological applications has ignited significant interest. However, its inherent non-polarity poses a critical challenge for its effective integration within biological systems. This study delves into the intricate physicochemical characteristics of the innovative [C60 + NO] complex using density functional theory and time-dependent density functional theory. The computational analyses encompass molecular charge, surface electrostatic potential, and dipole moment evaluations. Impressively, the dipole moment of the [C60 + NO] complex significantly increases to 12.92 D. Meticulous surface analysis reveals a subtle interplay between molecular structures, indicating weak interactions. The analysis of the absorption spectrum unveils a noteworthy red-shift of 200 nm subsequent to complex formation. To elucidate the electron transfer mechanisms, we explore photo-induced electron transfer through CAM-B3LYP. This exploration elucidates intricate pathways governing electron transfer, with complementary insights gleaned from Marcus theory's outputs, especially the Gibbs free energy of electron transfer. Changes in the physicochemical properties of approaching C60 and NO molecules reveal interesting results compared to separate molecules. These findings resonate profoundly in the context of potential biological and pharmaceutical utilization. With implications for the biomedical area, the outcomes linked to the [C60 + NO] complex kindle optimism for pioneering biomedical applications.

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.

Benchmarking two-photon absorption strengths of rhodopsin chromophore models with CC3 and CCSD methodologies: An assessment of popular density functional approximations

This work presents the investigations of the impact of an increasing electron correlation in the hierarchy of coupled-cluster methods, i.e., CC2, CCSD, and CC3, on two-photon absorption (2PA) strengths for the lowest excited state of the minimal rhodopsin's chromophore model-cis-penta-2,4-dieniminium cation (PSB3). For a larger chromophore's model [4-cis-hepta-2,4,6-trieniminium cation (PSB4)], CC2 and CCSD calculations of 2PA strengths were performed. Additionally, 2PA strengths predicted by some popular density functional theory (DFT) functionals differing in HF exchange contribution were assessed against the reference CC3/CCSD data. For PSB3, the accuracy of 2PA strengths increases in the following order: CC2 < CCSD < CC3, with the CC2 deviation from both higher-level methods exceeding 10% at 6-31+G* basis sets and 2% at aug-cc-pVDZ basis set. However, for PSB4, this trend is reversed and CC2-based 2PA strength is larger than the corresponding CCSD value. Among the DFT functionals investigated, CAM-B3LYP and BHandHLYP provide 2PA strengths in best compliance with reference data, however, with the error approaching an order of magnitude.

Ultrafast Evaluation of Two-Photon Absorption with Simplified Time-Dependent Density Functional Theory

This work presents the theoretical background to evaluate two-photon absorption (2PA) cross-sections in the framework of simplified time-dependent density functional theory (sTD-DFT). Our new implementation allows the ultrafast evaluation of 2PA cross-sections for large molecules based on a regular DFT ground-state determinant as well as a variant employing our tight-binding sTD-DFT-xTX flavor for very large systems. The method is benchmarked against higher-level calculations for trans-stilbene and typical fluorescent protein chromophores. For eGFP, a quadrupolar chromophore and its branched version, the flavine mono-nucleotide, and the iLOV protein, we compare sTD-DFT 2PA spectra to experimental ones. This includes extension and testing of our all-atom quantum chemistry methodology for the evaluation of 2PA for a system of ∼2000 atoms, providing striking agreement with the experimental spectrum.

Quantum chemistry study of the multiphoton absorption in enhanced green fluorescent protein at the single amino acid residue level

The chromophore (CRO) of fluorescent proteins (FPs) is embedded in a complex environment that is a source of specific interactions with the CRO. Understanding how these interactions influence FPs spectral properties is important for a directed design of novel markers with desired characteristics. In this work, we apply computational chemistry methods to gain insight into one-, two- and three-photon absorption (1PA, 2PA, 3PA) tuning in enhanced green fluorescent protein (EGFP). To achieve this goal, we built EGFP models differing in: i) number and position of hydrogen-bonds (h-bonds) donors to the CRO and ii) the electric field, as approximated by polarizable force  field, acting on the CRO. We  find that h-bonding to the CRO's phenolate oxygen results in stronger one- and multiphoton absorption. The brighter absorption can be also achieved by creating more positive electric field near the CRO's phenolate moiety. Interestingly, while individual CRO-environment h-bonds usually enhance 1PA and 2PA, it takes a few h-bond donors to enhance 3PA. Clearly, response of the absorption intensity to many-body effects depends on the excitation mechanism.  We further employ symmetry-adapted perturbation theory (SAPT) to reveal excellent (2PA) and good (3PA) correlation of multiphoton intensity with electrostatic and induction interaction energies. This points to importance of accounting for mutual CRO-environment polarization in quantitative calculations of absorption spectra in FPs.

Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges

Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin–orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm–Dancoff approximation to TDDFT, spin–flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine–Ag20 complex), we emphasize that caution is needed to properly account for system–environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.

The role of hydrogen bonds and electrostatic interactions in enhancing two-photon absorption in green and yellow fluorescent proteins

The spectral properties of fluorescent proteins (FPs) depend on the protein environment of the chromophore (CRO). A deeper understanding of the CRO - environment interactions in terms of FPs spectral characteristics will allow for a rational design of novel markers with desired properties. Here, we are taking a step towards achieving this important goal. With the time-dependent density functional theory (TDDFT), we calculate one- and two-photon absorption (OPA and TPA) spectra for 5 green FPs (GFPs) and 3 yellow FPs (YFPs) differing in amino acid sequence. The goal is to reveal a role of: (i) electrostatic interactions, (ii) hydrogen-bonds (h-bonds), and (iii) h-bonds together with distant electrostatic field in absorption spectra tuning. Our results point to design hypothesis towards FPs optimised for TPA-based applications. Both h-bonds and electrostatic interactions co-operate in enhancing TPA cross-section (σ TPA ) for the S 0 ->S 1   transition in GFPs. Furthermore, it seems that details of h-bonds network in the CRO's vicinity influences σ TPA   response to CRO - environment electrostatic interactions in YFPs. We postulate that engineering FPs with more hydrophilic CRO's environment can lead to greater σ TPA . We also find that removing h-bonds formed with the CRO's phenolate leads to TPA enhancement for transition to higher excited states than S 1 . Particularly Y145 and T203 residues are important in this regard.

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.

E/Z Molecular Photoswitches Activated by Two-Photon Absorption: Comparison between Different Families

Nonlinear optical techniques as two-photon absorption (TPA) have raised relevant interest within the last years due to the capability to excite chromophores with photons of wavelength equal to only half of the corresponding one-photon absorption energy. At the same time, its probability being proportional to the square of the light source intensity, it allows a better spatial control of the light-induced phenomenon. Although a consistent number of experimental studies focus on increasing the TPA cross section, very few of them are devoted to the study of photochemical phenomena induced by TPA. Here, we show a design strategy to find suitable E/Z photoswitches that can be activated by TPA. A theoretical approach is followed to predict the TPA cross sections related to different excited states of various photoswitches’ families, finally concluding that protonated Schiff-bases (retinal)-like photoswitches outperform compared to the others. The donor-acceptor substitution effect is therefore rationalized for the successful TPA activatable photoswitch, in order to maximize its properties, finally also forecasting a possible application in optogenetics. Some experimental measurements are also carried out to support our conclusions.

Removing artifacts in polarizable embedding calculations of one- and two-photon absorption spectra of fluorescent proteins

The multiscale calculations involving excited states may suffer from the electron spill-out (ESO) problem. This seems to be especially the case when the environment of the core region, described with the electronic structure method, is approximated by a polarizable force field. The ESO effect often leads to incorrect physical character of electronic excitations, spreading outside the quantum region, which, in turn, results in erroneous absorption spectra. In this work, we investigate means to remove the artifacts in one-photon absorption (OPA) and two-photon absorption (TPA) spectra of green and yellow fluorescent protein representatives. This includes (i) using different basis sets, (ii) extending the core subsystem beyond the chromophore, (iii) modification of polarization interaction between the core region and its environment, and (iv) including the Pauli repulsion through effective core potentials (ECPs). Our results clearly show that ESO is observed when diffuse functions are used to assemble the multielectron wave function regardless of the exchange-correlation functional used. Furthermore, extending the core region, thus accounting for exchange interactions between the chromophore and its environment, leads to even more spurious excited states. Also, damping the interactions between the core subsystem and the polarizable force field is hardly helpful. In contrast, placing ECPs in the position of sites creating the embedding potential leads to the removal of artificious excited states that presumably should not be observed in the OPA and TPA spectra. We prove that it is a reliable and cost-effective approach for systems where the covalent bond(s) between the core region and its environment must be cut.

Benchmarking Magnetizabilities with Recent Density Functionals

We have assessed the accuracy of the magnetic properties of a set of 51 density functional approximations, including both recently published and already established functionals. The accuracy assessment considers a series of 27 small molecules and is based on comparing the predicted magnetizabilities to literature reference values calculated using coupled-cluster theory with full singles and doubles and perturbative triples [CCSD(T)] employing large basis sets. The most accurate magnetizabilities, defined as the smallest mean absolute error, are obtained with the BHandHLYP functional. Three of the six studied Berkeley functionals and the three range-separated Florida functionals also yield accurate magnetizabilities. Also, some older functionals like CAM-B3LYP, KT1, BHLYP (BHandH), B3LYP, and PBE0 perform rather well. In contrast, unsatisfactory performance is generally obtained with Minnesota functionals, which are therefore not recommended for calculations of magnetically induced current density susceptibilities and related magnetic properties such as magnetizabilities and nuclear magnetic shieldings. We also demonstrate that magnetizabilities can be calculated by numerical integration of magnetizability density; we have implemented this approach as a new feature in the gauge-including magnetically induced current (GIMIC) method. Magnetizabilities can be calculated from magnetically induced current density susceptibilities within this approach even when analytical approaches for magnetizabilities as the second derivative of the energy have not been implemented. The magnetizability density can also be visualized, providing additional information that is not otherwise easily accessible on the spatial origin of magnetizabilities.

Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy

Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore-environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore-environment configurations with Franck-Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Employing an average Franck-Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.

Two-photon conversion of a bacterial phytochrome

In nature, sensory photoreceptors underlie diverse spatiotemporally precise and generally reversible biological responses to light. Photoreceptors also serve as genetically encoded agents in optogenetics to control by light organismal state and behavior. Phytochromes represent a superfamily of photoreceptors that transition between states absorbing red light (Pr) and far-red light (Pfr), thus expanding the spectral range of optogenetics to the near-infrared range. Although light of these colors exhibits superior penetration of soft tissue, the transmission through bone and skull is poor. To overcome this fundamental challenge, we explore the activation of a bacterial phytochrome by a femtosecond laser emitting in the 1 μm wavelength range. Quantum chemical calculations predict that bacterial phytochromes possess substantial two-photon absorption cross sections. In line with this notion, we demonstrate that the photoreversible Pr ↔ Pfr conversion is driven by two-photon absorption at wavelengths between 1170 and 1450 nm. The Pfr yield was highest for wavelengths between 1170 and 1280 nm and rapidly plummeted beyond 1300 nm. By combining two-photon activation with bacterial phytochromes, we lay the foundation for enhanced spatial resolution in optogenetics and unprecedented penetration through bone, skull, and soft tissue.

What is the Optimal Size of the Quantum Region in Embedding Calculations of Two-Photon Absorption Spectra of Fluorescent Proteins?

We systematically investigate an impact of the size and content of a quantum (QM) region, treated at the density functional theory level, in embedding calculations on one- (OPA) and two-photon absorption (TPA) spectra of the following fluorescent proteins (FPs) models: Aequorea victoria green FP (avGFP) with neutral (avGFP-n) and anionic (avGFP-a) chromophore as well as Citrine FP. We find that amino acid (a.a.) residues as well as water molecules hydrogen-bonded (h-bonded) to the chromophore usually boost both OPA and TPA processes intensity. The presence of hydrophobic a.a. residues in the quantum region also non-negligibly affects both absorption spectra but decreases absorption intensity. We conclude that to reach a quantitative description of OPA and TPA spectra in multiscale modeling of FPs, the quantum region should consist of a chromophore and most of a.a. residues and water molecules in a radius of 0.30–0.35 nm (ca. 200–230 atoms) when the remaining part of the system is approximated by the electrostatic point-charges. The optimal size of the QM region can be reduced to 80–100 atoms by utilizing a more advanced polarizable embedding model. We also find components of the QM region that are specific to a FP under study. We propose that the F165 a.a. residue is important in tuning the TPA spectrum of avGFP-n but not other investigated FPs. In the case of Citrine, Y203 and M69 a.a. residues must definitely be part of the QM subsystem. Furthermore, we find that long-range electrostatic interactions between the QM region and the rest of the protein cannot be neglected even for the most extensive QM regions (ca. 350 atoms).