Computing Excited States of Very Large Systems with Range-Separated Hybrid Functionals and the Exact Integral Simplified Time-Dependent Density Functional Theory (XsTD-DFT)

Simplified quantum chemistry (sQC) methods can routinely compute excited states for very large systems in an "all-atom" fashion. They are viable alternatives to regular multiscale schemes. sQC methods have the advantage of accounting explicitly for all of the environment at a quantum mechanical (QM) level. The treatment of charge-transfer states is now improved by the native implementation of range-separated hybrid (RSH) exchange-correlation functionals into the eXact integral simplified time-dependent density functional theory (XsTD-DFT). After the RSH XsTD-DFT/TDA scheme was benchmarked, XsTD-DFT(/TDA) ultraviolet/visible absorption, circular dichroism (CD), and/or two-photon absorption (2PA) spectra were directly compared to the results of experiments for four challenging and increasingly large systems: eYFP model system, Λ-shaped multimodular D-π-A-π-D'-π-A-π-D chromophore, mixed donor/acceptor ligand Pd(II) double cage [3BF4@Pd4DmA8-m], and the photoactive yellow protein (PYP). Among the results, this study shows that an "all-atom" approach is unavoidable for reproducing absorption and CD spectra of PYP because one of the main transitions involves a local excitation on a tryptophan.

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

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

Non-Degenerate Two-Photon Absorption of Fluorescent Protein Chromophores

Two-photon absorption (2PA), where a pair of photons are absorbed simultaneously, is recognized as a potent bioimaging technique, which depends on the quantified 2PA probability, defined as cross-section (σ2PA). The absorbed photons either have equivalent (ω1 = ω2) or different frequencies (ω1 ≠ ω2), where the former is degenerate 2PA (D-2PA) and the latter is nondegenerate 2PA (ND-2PA). ND-2PA is of particular interest since it is a promising imaging technology with flexibility of photon frequencies and enhanced cross sections, however, it remains a relatively unexplored area compared to D-2PA. This work utilizes time-dependent density functional theory (TD-DFT) and second-order approximate coupled-cluster with the resolution-of-identity approximation (RI-CC2), for the excitation from S0 to S1, to investigate σD-2PA and σND-2PA of FP chromophore models. Interestingly, comparing CAM-B3LYP with the RI-CC2 computations shows qualitative and, in fact, near quantitative agreement in the computed improvements of σND-2PA for comparable (relative) frequency detunings, despite the known underestimations of 2PA cross sections, for TD-DFT results relative to RI-CC2 values. As expected from the 2-state model, the computed values of σND-2PA are quantitatively larger than σD-2PA, where chromophores with the largest values of σD-2PA show greater potential for σND-2PA improvement. Anionic chromophores demonstrated improvements up to 14%, while substantial enhancements were observed in neutral chromophores with some achieving a 30% increase. This work investigates the ND-2PA photophysical characteristics of FP chromophores and identifies qualitative patterns in the computed properties of ND-2PA relative to D-2PA.

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.

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.

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.

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

Multiscale Modeling of Two-Photon Probes for Parkinson's Diagnostics Based on Monoamine Oxidase B Biomarker

Monoamine oxidase B (MAO-B) is a potential biomarker for Parkinson's disease (PD), a neurodegenerative disease associated with the loss of motor activities in human subjects. The disease state is associated with dopamine deprival, and so the inhibitors of MAO-B can serve as therapeutic drugs for PD. Since the expression level of MAO-B directly correlates to the disease progress, the distribution and population of this enzyme can be employed to monitor disease development. One of the approaches available for estimating the population is two-photon imaging. The ligands used for two-photon imaging should have high binding affinity and binding specificity toward MAO-B along with significant two-photon absorption cross sections when they are bound to the target. In this article, we study using a multiscale modeling approach, the binding affinity and spectroscopic properties (one- and two-photon absorption) of three (Flu1, Flu2, Flu3) of the currently available probes for monitoring the MAO-B level. We report that the binding affinity of the probes can be explained using the molecular size and binding cavity volume. The experimentally determined one-photon absorption spectrum is well reproduced by the employed QM/MM approaches, and the most accurate spectral shifts, on passing from one probe to another, are obtained at the coupled-cluster (CC2) level of theory. An important conclusion from this study is also the demonstration that intrinsic molecular two-photon absorption strengths (δ^2PA) increase in the order δ^2PA (Flu1) > δ^2PA (Flu2) > δ^2PA (Flu3). This is in contrast with experimental data, which predict similar values of two-photon absorption cross sections for Flu1 and Flu3. We demontrate, based on the results of electronic-structure calculations for Flu1 that this discrepancy cannot be explained by an explicit account for neighboring residues (which could lead to charge transfer between a probe and neighboring aromatic amino acids thus boosting δ^2PA). In summary, we show that the employed multiscale approach not only can optimize two-photon absorption properties and verify binding affinity, but it can also help in detailed analyses of experimental data.

Two-Photon Absorption Cross-Sections in Fluorescent Proteins Containing Non-canonical Chromophores Using Polarizable QM/MM

Multi-photon absorption properties, particularly two-photon absorption (2PA), of fluorescent proteins (FPs) have made them attractive tools in deep-tissue clinical imaging. Although the diversity of photophysical properties for FPs is wide, there are some caveats predominant among the existing FP variants that need to be overcome, such as low quantum yields and small 2PA cross-sections. From a computational perspective, Salem et al. (2016) suggested the inclusion of non-canonical amino acids in the chromophore of the red fluorescent protein DsRed, through the replacement of the tyrosine amino acid. The 2PA properties of these new non-canonical chromophores (nCCs) were determined in vacuum, i.e., without taking into account the protein environment. However, in the computation of response properties, such as 2PA cross-sections, the environment plays an important role. To account for environment and protein-chromophore coupling effects, quantum mechanical/molecular mechanical (QM/MM) schemes can be useful. In this work, the polarizable embedding (PE) model is employed along with time-dependent density functional theory to describe the 2PA properties of a selected set of chromophores made from non-canonical amino acids as they are embedded in the DsRed protein matrix. The objective is to provide insights to determine whether or not the nCCs could be developed and, thus, generate a new class of FPs. Results from this investigation show that within the DsRed environment, the nCC 2PA cross-sections are diminished relative to their values in vacuum. However, further studies toward understanding the 2PA limit of these nCCs using different protein environments are needed.

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

Performance of DFT functionals with different percentages of exact Hartree-Fock exchange energy (EX) is assessed for recovery of the CC2 reference one- (OPA) and two-photon absorption (TPA) spectra of fluorescent proteins chromophores in vacuo. The investigated DFT functionals, together with their EX contributions are BLYP (0%), B3LYP (20%), B1LYP (25%), BHandHLYP (50%), and CAM-B3LYP (19% at short range and 65% at long range). Our test set consists of anionic and neutral chromophores as naturally occurring in the fluorescent proteins. For the first time, we compare TDDFT and CC2 methods for higher excited states than the S₁ state, exhibiting relatively large TPA intensity. Our TDDFT results for neutral chromophores reveal an increase in excitation energies as well as TPA and OPA intensities errors, compared to CC2-derived results, as the DFT functional contains less exact exchange. The long-range-corrected CAM-B3LYP functional performs the best, closely followed by BHandHLYP, while BLYP usually significantly underestimates all investigated spectral properties, hence being the worst in reproducing the reference CC2 results. The hybrid B3LYP and B1LYP functionals can be roughly placed in between. We propose that TDDFT may underestimate the TPA intensities for neutral chromophores of fluorescent proteins due to underestimated oscillator strengths between some excited states. In the case of anionic chromophores, we find that B3LYP and B1LYP functionals overcome others in terms of reproducing CC2 excitation energies. On the other hand, however, TPA intensity is usually significantly underestimated, and in this respect, CAM-B3LYP functional seems to be again superior. In contrast to the case of neutral chromophores, it seems that a large magnitude of excited-state dipole moments or changes in dipole moments upon excitation may be the driving force behind high TPA transition moments.

Facile Synthesis of a Terephthalic Acid-Based Organic Fluorophore with Strong and Color-Tunable Emission in Both Solution and Solid States for LED Applications

Dual-state emission (DSE) fluorophores with strong photoluminescence in both solution and solid states have wide applications in fluorescent probes and photoelectric devices. However, most of the existing DSE fluorophores involve complex synthetic strategies or only one or two emission colors in the solid state, which may hinder their practical applications. Herein, we report a facile and effective strategy to fabricate a kind of novel adjustable DSE fluorophore TPAA-Cu based on the small molecule terephthalic acid (TPA) and ascorbic acid (AA) with CuCl₂ as the catalyst. Not only can the TPAA-Cu solution emit blue, blue-green, and yellow colors depending upon the pH of solution, but also it shows reversible pH-tunable fluorescence colors. More importantly, corresponding TPAA-Cu solid fluorophores with blue, blue-green, and yellow colors can also be obtained with the quantum yields of 16.3, 26.7, and 29.6%, respectively.

Two-photon absorption spectroscopy of trans-stilbene, cis-stilbene, and phenanthrene: Theory and experiment

Two-photon absorption (2PA) spectroscopy provides complementary, and sometimes more detailed, information about the electronic structure of a molecule relative to one-photon absorption (1PA) spectroscopy. However, our understanding of the 2PA processes is rather limited due to technical difficulties in measuring experimental 2PA spectra and theoretical challenges in computing higher-order molecular properties. This paper examines the 2PA spectroscopy of trans-stilbene, cis-stilbene, and phenanthrene by a combined experimental and theoretical approach. The broadband 2PA spectra of all three compounds are measured under identical conditions in order to facilitate a direct comparison of the absolute 2PA cross sections in the range 3.5-6.0 eV. For comparison, the theoretical 2PA cross sections are computed using the equation-of-motion coupled-cluster method with single and double substitutions. Simulated 2PA spectra based on the calculations reproduce the main features of the experimental spectra in solution, although the quantitative comparison is complicated by a number of uncertainties, including limitations of the theoretical model, vibronic structure, broadening of the experimental spectra, and solvent effects. The systematic comparison of experimental and theoretical spectra for this series of structurally similar compounds provides valuable insight into the nature of 2PA transitions in conjugated molecules. Notably, the orbital character and symmetry-based selection rules provide a foundation for interpreting the features of the experimental 2PA spectra in unprecedented detail.

Steady-State Linear and Non-linear Optical Spectroscopy of Organic Chromophores and Bio-macromolecules

Bio-macromolecules as DNA, lipid membranes and (poly)peptides are essential compounds at the core of biological systems. The development of techniques and methodologies for their characterization is therefore necessary and of utmost interest, even though difficulties can be experienced due to their intrinsic complex nature. Among these methods, spectroscopies, relying on optical properties are especially important to determine their macromolecular structures and behaviors, as well as the possible interactions and reactivity with external dyes—often drugs or pollutants—that can (photo)sensitize the bio-macromolecule leading to eventual chemical modifications, thus damages. In this review, we will focus on the theoretical simulation of electronic spectroscopies of bio-macromolecules, considering their secondary structure and including their interaction with different kind of (photo)sensitizers. Namely, absorption, emission and electronic circular dichroism (CD) spectra are calculated and compared with the available experimental data. Non-linear properties will be also taken into account by two-photon absorption, a highly promising technique (i) to enhance absorption in the red and infra-red windows and (ii) to enhance spatial resolution. Methodologically, the implications of using implicit and explicit solvent, coupled to quantum and thermal samplings of the phase space, will be addressed. Especially, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are explored for a comparison with solely QM methods, in order to address the necessity to consider an accurate description of environmental effects on spectroscopic properties of biological systems.

Molecular cocrystals: design, charge-transfer and optoelectronic functionality

This perspective article primarily focuses on the research work related to optoelectronic properties of organic charge transfer cocrystals.

Protein charge transfer absorption spectra: an intrinsic probe to monitor structural and oligomeric transitions in proteins

The utility of ProCharTS as an intrinsic spectral probe to track protein aggregation and monitor conformational changes is reported.

Prediction of two-photon absorption enhancement in red fluorescent protein chromophores made from non-canonical amino acids

Two-photon spectroscopy of fluorescent proteins is a powerful bio-imaging tool known for deep tissue penetration and little cellular damage. Being less sensitive than the one-photon microscopy alternatives, a protein with a large two-photon absorption (TPA) cross-section is needed. Here, we use time-dependent density functional theory (TD-DFT) at the B3LYP and CAM-B3LYP/6-31+G(d,p) levels of theory to screen twenty-two possible chromophores that can be formed upon replacing the amino-acid Tyr66 that forms the red fluorescent protein (RFP) chromophore with a non-canonical amino acid. The two-level model for TPA was used to assess the properties (i.e., transition dipole moment, permanent dipole moment difference, and the angle between them) leading to the TPA cross-sections determined via response theory. Computing TPA cross-sections with B3LYP and CAM-B3LYP yields similar overall trends. Results using both functionals agree that the RFP-derived model of the Gold Fluorescent Protein chromophore (Model 20) has the largest intrinsic TPA cross-section at the optimized geometry. TPA was further computed for selected chromophores following conformational changes: variation of both the dihedral angle of the acylimine moiety and the tilt and twist angles between the rings of the chromophore. The TPA cross-section assumed an oscillatory trend with the rotation of the acylimine dihedral, and the TPA is maximized in the planar conformation for almost all models. Model 21 (a hydroxyquinoline derivative) is shown to be comparable to Model 20 in terms of TPA cross-section. The conformational study on Model 21 shows that the acylimine angle has a much stronger effect on the TPA than its tilt and twist angles. Having an intrinsic TPA ability that is more than 7 times that of the native RFP chromophore, Models 20 and 21 appear to be very promising for future experimental investigation.

Resolving Electronic Transitions in Synthetic Fluorescent Protein Chromophores by Magnetic Circular Dichroism

The detailed electronic structures of fluorescent chromophores are important for their use in imaging of living cells. A series of green fluorescent protein chromophore derivatives is examined by magnetic circular dichroism (MCD) spectroscopy, which allows the resolution of more bands than plain absorption and fluorescence. Observed spectral patterns are rationalized with the aid of time-dependent density functional theory (TDDFT) computations and the sum-over-state (SOS) formalism, which also reveals a significant dependence of MCD intensities on chromophore conformation. The combination of organic and theoretical chemistry with spectroscopic techniques also appears useful in the rational design of fluorescence labels and understanding of the chromophore's properties. For example, the absorption threshold can be heavily affected by substitution on the phenyl ring but not much on the five-member ring, and methoxy groups can be used to further tune the electronic levels.

Two-photon absorption of fluorescent protein chromophores incorporating non-canonical amino acids: TD-DFT screening and classical dynamics

Two-photon spectroscopy of fluorescent proteins is a powerful bio-imaging tool characterized by deep tissue penetration and little damage. However, two-photon spectroscopy has lower sensitivity than one-photon microscopy alternatives and hence a protein with a large two-photon absorption cross-section is needed. We use time-dependent density functional theory (TD-DFT) at the B3LYP/6-31+G(d,p) level of theory to screen twenty-two possible chromophores that can be formed upon replacing the amino-acid Tyr66 that forms the green fluorescent protein (GFP) chromophore with a non-canonical amino acid. A proposed chromophore with a nitro substituent was found to have a large two-photon absorption cross-section (29 GM) compared to other fluorescent protein chromophores as determined at the same level of theory. Classical molecular dynamics are then performed on a nitro-modified fluorescent protein to test its stability and study the effect of the conformational flexibility of the chromophore on its two-photon absorption cross-section. The theoretical results show that the large cross-section is primarily due to the difference between the permanent dipole moments of the excited and ground states of the nitro-modified chromophore. This large difference is maintained through the various conformations assumed by the chromophore in the protein cavity. The nitro-derived protein appears to be very promising as a two-photon absorption probe.

Open-Ended Recursive Approach for the Calculation of Multiphoton Absorption Matrix Elements

We present an implementation of single residues for response functions to arbitrary order using a recursive approach. Explicit expressions in terms of density-matrix-based response theory for the single residues of the linear, quadratic, cubic, and quartic response functions are also presented. These residues correspond to one-, two-, three- and four-photon transition matrix elements. The newly developed code is used to calculate the one-, two-, three- and four-photon absorption cross sections of para-nitroaniline and para-nitroaminostilbene, making this the first treatment of four-photon absorption in the framework of response theory. We find that the calculated multiphoton absorption cross sections are not very sensitive to the size of the basis set as long as a reasonably large basis set with diffuse functions is used. The choice of exchange–correlation functional, however, significantly affects the calculated cross sections of both charge-transfer transitions and other transitions, in particular, for the larger para-nitroaminostilbene molecule. We therefore recommend the use of a range-separated exchange–correlation functional in combination with the augmented correlation-consistent double-ζ basis set aug-cc-pVDZ for the calculation of multiphoton absorption properties.

Temperature and pressure effects on GFP mutants: explaining spectral changes by molecular dynamics simulations and TD-DFT calculations

By combining spectroscopic measurements under high pressure with molecular dynamics simulations and quantum mechanics calculations we investigate how sub-angstrom structural perturbations are able to tune protein function. We monitored the variations in fluorescence output of two green fluorescent protein mutants (termed Mut2 and Mut2Y, the latter containing the key T203Y mutation) subjected to pressures up to 600 MPa, at various temperatures in the 280-320 K range. By performing 150 ns molecular dynamics simulations of the protein structures at various pressures, we evidenced subtle changes in conformation and dynamics around the light-absorbing chromophore. Such changes explain the measured spectral tuning in the case of the sizable 120 cm(-1) red-shift observed for pressurized Mut2Y, but absent in Mut2. Previous work [Barstow et al., Proc. Natl. Acad. Sci. U. S. A., 2008, 105, 13362] on pressure effects on GFP also involved a T203Y mutant. On the basis of cryocooling X-ray crystallography, the pressure-induced fluorescence blue shift at low temperature (77 K) was attributed to key changes in relative conformation of the chromophore and Tyr203 phenol ring. At room temperature, however, a red shift was observed at high pressure, analogous to the one we observe in Mut2Y. Our investigation of structural variations in compressed Mut2Y also explains their result, bridging the gap between low-temperature and room-temperature high-pressure effects.

Benchmarking two-photon absorption cross sections: performance of CC2 and CAM-B3LYP

We investigate the performance of CC2 and TDDFT/CAM-B3LYP for the calculation of two-photon absorption (TPA) strengths and cross sections and contrast our results to a recent coupled cluster equation-of-motion (EOM-EE-CCSD) benchmark study [K. D. Nanda and A. I. Krylov, J. Chem. Phys., 2015, 142, 064118]. In particular, we investigate whether CC2 TPA strengths are significantly overestimated compared to higher-level coupled-cluster calculations for fluorescent protein chromophores. Our conclusion is that CC2 TPA strengths are only slightly overestimated compared to the reference EOM-EE-CCSD results and that previously published overestimated cross sections are a result of inconsistencies in the conversion of the TPA strengths to macroscopic units. TDDFT/CAM-B3LYP TPA strengths, on the other hand, are found to be 1.5 to 3 times smaller than the coupled-cluster reference for the molecular systems considered. The unsatisfactory performance of TDDFT/CAM-B3LYP can be linked to an underestimation of excited-state dipole moments predicted by TDDFT/CAM-B3LYP.

Two-photon absorption cross sections within equation-of-motion coupled-cluster formalism using resolution-of-the-identity and Cholesky decomposition representations: Theory, implementation, and benchmarks

The equation-of-motion coupled-cluster (EOM-CC) methods provide a robust description of electronically excited states and their properties. Here, we present a formalism for two-photon absorption (2PA) cross sections for the equation-of-motion for excitation energies CC with single and double substitutions (EOM-CC for electronically excited states with single and double substitutions) wave functions. Rather than the response theory formulation, we employ the expectation-value approach which is commonly used within EOM-CC, configuration interaction, and algebraic diagrammatic construction frameworks. In addition to canonical implementation, we also exploit resolution-of-the-identity (RI) and Cholesky decomposition (CD) for the electron-repulsion integrals to reduce memory requirements and to increase parallel efficiency. The new methods are benchmarked against the CCSD and CC3 response theories for several small molecules. We found that the expectation-value 2PA cross sections are within 5% from the quadratic response CCSD values. The RI and CD approximations lead to small errors relative to the canonical implementation (less than 4%) while affording computational savings. RI/CD successfully address the well-known issue of large basis set requirements for 2PA cross sections calculations. The capabilities of the new code are illustrated by calculations of the 2PA cross sections for model chromophores of the photoactive yellow and green fluorescent proteins.

Two-Photon Absorption in Fluorescent Protein Chromophores: TDDFT and CC2 Results

Two-photon spectroscopy of fluorescent proteins is a powerful bioimaging tool. Considerable effort has been made to measure absolute two-photon absorption (TPA) for the available fluorescent proteins. Being a technically involved procedure, there is significant variation in the published experimental measurements even for the same protein. In this work, we present a time-dependent density functional theory (TDDFT) study on isolated chromophores comparing the ability of four functionals (PBE0, B3LYP, CAM-B3LYP, and LC-BLYP) combined with the 6-31+G(d,p) basis set to reproduce averaged experimental TPA energies and cross sections. The TDDFT energies and TPA cross sections are also compared to corresponding CC2/6-31+G(d,p) results for excitation to S1 for the five smallest chromophores. In general, the computed TPA energies are less functional dependent than the TPA cross sections. The variation between functionals is more pronounced when higher-energy transitions are studied. Changes to the conformation of a chromophore are shown to change the TPA cross-section considerably. This adds to the difficulty of comparing an isolated chromophore to the one embedded in the protein environment. All functionals considered give moderate agreement with the corresponding CC2 results; in general, the TPA cross sections determined by TDDFT are 1.5-10 times smaller than the corresponding CC2 values for excitation to S1. LC-BLYP and CAM-B3LYP give erroneously large TPA cross sections in the higher-energy regions. On the other hand, B3LYP and PBE0 yield values that are of the same order of magnitude and in some cases very close to the averaged experimental data. Thus, based on the results reported here, B3LYP and PBE0 are the preferred functionals for screening chromphores for TPA. However, at best, TDDFT can be used to semiquantitatively scan chromophores for potential TPA probes and highlight spectroscopic peaks that could be present in the mature protein.

Chemical control of channel interference in two-photon absorption processes

The two-photon absorption (TPA) process is the simplest and hence the most studied nonlinear optical phenomenon, and various aspects of this process have been explored in the past few decades, experimentally as well as theoretically. Previous investigations have shown that the two-photon (TP) activity of a molecular system can be tuned, and at present, performance-tailored TP active materials are easy to develop by monitoring factors such as length of conjugation, dimensionality of charge-transfer network, strength of donor-acceptor groups, polarity of solvents, self-aggregation, H-bonding, and micellar encapsulation to mention but a few. One of the most intriguing phenomena affecting the TP activity of a molecule is channel interference. The phrase "channel interference" implies that if the TP transition from one electronic state to another involves more than one optical pathway or channel, characterized by the corresponding transition dipole moment (TDM) vectors, the channels may interfere with each other depending upon the angles between the TDM vectors and hence can either increase (constructive interference) or decrease (destructive interference) the overall TP activity of a system to a significant extent. This phenomenon was first pointed out by Cronstrand, Luo, and Ågren [Chem. Phys. Lett. 2002, 352, 262-269] in two-dimensional systems (i.e., only involving two components of the transition moment vectors). For three-dimensional molecules, an extended version of this idea was required. In order to fill this gap, we developed a generalized model for describing and exploring channel interference, valid for systems of any dimensionality. We have in particular applied it to through-bond (TB) and through-space (TS) charge-transfer systems both in gas phase and in solvents with different polarities. In this Account, we will, in addition to briefly describing the concept of channel interference, discuss two key findings of our recent work: (1) how to control the channel interference by chemical means, and (2) the role of channel interference in the anomalous solvent dependence of certain TP chromophores. For example, we will show that simple structurally induced changes in certain dihedral angles of the well-known betaine dye (TB type) will help fine-tune the constructive channel interference and hence increase the overall TP activity of molecules with this general TP channel structure. Another intriguing result we will discuss is observed for a tweezer-trinitrofluorinone complex (TS type) where, on moving from polar to essentially nonpolar solvents, the nature of the channel interference switches from destructive to constructive, leading to a net abnormal solvent dependence of the TP activity of the system. The present Account highlights the usefulness of the channel interference effect and establishes it as a new and unique way of controlling the TP transition probability in different types of three-dimensional molecules.