Nonlinear response time-dependent density functional theory combined with the effective fragment potential method

The eXact integral simplified time-dependent density functional theory (XsTD-DFT)

In the framework of simplified quantum chemistry methods, we introduce the eXact integral simplified time-dependent density functional theory (XsTD-DFT). This method is based on the simplified time-dependent density functional theory (sTD-DFT), where all semi-empirical two-electron integrals are replaced by exact one- and two-center two-electron integrals, while other approximations from sTD-DFT are kept. The performance of this new parameter-free XsTD-DFT method was benchmarked on excited state and (non)linear response properties, including ultra-violet/visible absorption, first hyperpolarizability, and two-photon absorption (2PA). For a set of 77 molecules, the results from the XsTDA approach were compared to the TDA data. XsTDA/B3LYP excitation energies only deviate on average by 0.14 eV from TDA while drastically cutting computational costs by a factor of 20 or more depending on the energy threshold chosen. The absolute deviations of excitation energies with respect to the full scheme are decreasing with increasing system size, showing the suitability of XsTDA/XsTD-DFT to treat large systems. Comparing XsTDA and its predecessor sTDA, the new scheme generally improves excitation energies and oscillator strengths, in particular, for charge transfer states. TD-DFT first hyperpolarizability frequency dispersions for a set of push-pull π-conjugated molecules are faithfully reproduced by XsTD-DFT, while the previous sTD-DFT method provides redshifted resonance energy positions. Excellent performance with respect to the experiment is observed for the 2PA spectrum of the enhanced green fluorescent protein. The obtained robust accuracy similar to TD-DFT at a fraction of the computational cost opens the way for a plethora of applications for large systems and in high throughput screening studies.

One- and two-photon absorption spectra of organoboron complexes: vibronic and environmental effects

We synthesized a series of four parent aza-β-ketoiminate organoboron complexes and performed spectroscopic studies using both experimental and computational techniques. We studied how benzannulation influences the vibronic structure of the UV/Vis absorption bands with a focus on the bright lowest-energy π → π* electronic excitation. Theoretical simulations, accounting for inhomogeneous broadening effects using different embedding schemes, allowed gaining in-depth insights into the observed differences in band shapes induced by structural modifications. We observed huge variations in the distributions of vibronic transitions depending on the position of benzannulation. By and large, the harmonic approximation combined with the adiabatic hessian model delivers qualitatively correct band shapes for the one-photon absorption spectra, except in one case. We also assessed the importance of non-Condon effects (accounted for by the linear term in Herzberg-Teller expansion of the dipole moment) for S0 → S1 band shapes. It turned out that non-Condon contributions have no effect on the band shape in one-photon absorption spectra. In contrast, these effects significantly change the Franck-Condon band shapes of the two-photon absorption spectra. For one of the studied organoboron complexes we also performed a preliminary exploration of mechanical anharmonicity, resulting in an increase of the intensity of the 0-0 transition, which improves the agreement with the experimental data compared to the harmonic model.

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.

Photoswitchable Molecular Units with Tunable Nonlinear Optical Activity: A Theoretical Investigation

The first-, second-, and third-order molecular nonlinear optical properties, including two-photon absorption of a series of derivatives, involving two dithienylethene (DTE) groups connected by several molecular linkers (bis(ethylene-1,2-dithiolato)Ni- (NiBDT), naphthalene, quasilinear oligothiophene chains), are investigated by employing density functional theory (DFT). These properties can be efficiently controlled by DTE switches, in connection with light of appropriate frequency. NiBDT, as a linker, is associated with a greater contrast, in comparison to naphthalene, between the first and second hyperpolarizabilities of the "open-open" and the "closed-closed" isomers. This is explained by invoking the low-lying excited states of NiBDT. It is shown that the second hyperpolarizability can be used as an index, which follows the structural changes induced by photochromism. Assuming a Förster type transfer mechanism, the intramolecular excited-state energy transfer (EET) mechanism is studied. Two important parameters related to this are computed: the electronic coupling (VDA) between the donor and acceptor fragments as well as the overlap between the absorption and emission spectra of the donor and acceptor groups. NiBDT as a linker is associated with a low electronic coupling, VDA, value. We found that VDA is affected by molecular geometry. Our results predict that the linker strongly influences the communication between the open-closed DTE groups. The sensitivity of the molecular nonlinear optical properties could assist with identification of molecular isomers.

Impact of the Protein Environment on Two-Photon Absorption Cross-Sections of the GFP Chromophore Anion Resolved at the XMCQDPT2 Level of Theory

The search for fluorescent proteins with large two-photon absorption (TPA) cross-sections and improved brightness is required for their efficient use in bioimaging. Here, we explored the impact of a single-point mutation close to the anionic form of the GFP chromophore on its TPA activity. We considered the lowest-energy transition of EGFP and its modification EGFP T203I. We focused on a methodology for obtaining reliable TPA cross-sections for mutated proteins, based on conformational sampling using molecular dynamics simulations and a high-level XMCQDPT2-based QM/MM approach. We also studied the numerical convergence of the sum-over-states formalism and provide direct evidence for the applicability of the two-level model for calculating TPA cross-sections in EGFP. The calculated values were found to be very sensitive to changes in the permanent dipole moments between the ground and excited states and highly tunable by internal electric field of the protein environment. In the case of the GFP chromophore anion, even a single hydrogen bond was shown to be capable of drastically increasing the TPA cross-section. Such high tunability of the nonlinear photophysical properties of the chromophore anions can be used for the rational design of brighter fluorescent proteins for bioimaging using two-photon laser scanning microscopy.

Curing the Divergence in Time-Dependent Density Functional Quadratic Response Theory

The adiabatic approximation in time-dependent density functional theory is known to give an incorrect pole structure in the quadratic response function, leading to unphysical divergences in excited state-to-state transition probabilities and hyperpolarizabilties. We find the form of the exact quadratic response kernel and derive a practical and accurate approximation that cures the divergence. We demonstrate our results on excited state-to-state transition probabilities of a model system and of the LiH molecule.

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

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.

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.

Ultrafast relaxation investigated by photoelectron circular dichroism: an isomeric comparison of camphor and fenchone

We study the isomeric effects using time resolved photoelectron circular dichroism (TR-PECD). Using a (1 + 1') pump-probe ionisation scheme with photoelectrons collected by the velocity map imaging technique, we compare the relaxation dynamics from the 3s-Rydberg state in 1R,4R-(+)-camphor with the one in its chiral isomer, 1R,4S-(-)-fenchone [Comby et al., 2016, JPCL, 7, 4514]. Our measurements revealed a similar lifetime for both isomers. However, the circular dichroism in the photoelectron angular distribution decays exponentially in ∼730 fs from a +9% forward amplitude during the first hundreds of femtoseconds to reach an asymptotic -2% backward amplitude. This time-scale is drastically shorter than in fenchone. Our analysis allows us to evaluate the impact of the anisotropy of excitation; the relaxation dynamics, following photoexcitation by the linearly polarized pump, is then compared to that induced by a circularly polarized pump pulse (CPL). With such a CPL pump, we then retrieve time constants of our chiral observables similar to the ones recorded in fenchone. Quantum and classical simulations are developed and used to decipher the dependence of the PECD on the anisotropy of excitation and the spatial distribution of the 3s-Rydberg electron wavefunction. Our experimental investigations, supported by our simulations, suggest that varying the pump ellipticity enables us to reveal the breakdown of the Franck-Condon approximation.

All-Atom Quantum Mechanical Calculation of the Second-Harmonic Generation of Fluorescent Proteins

Fluorescent proteins (FPs) are biotags of choice for second-harmonic imaging microscopy (SHIM). Because of their large size, computing their second-harmonic generation (SHG) response represents a great challenge for quantum chemistry. In this contribution, we propose a new all-atom quantum mechanics methodology to compute SHG of large systems. This is now possible because of two recent implementations: the tight-binding GFN2-xTB method to optimize geometries and a related version of the simplified time-dependent density functional theory (sTD-DFT-xTB) to evaluate quadratic response functions. In addition, a new dual-threshold configuration selection scheme is introduced to reduce the computational costs while retaining overall similar accuracy. This methodology was tested to evaluate the SHG of the proteins iLOV and bacteriorhodopsin (bR). In the case of bR, quantitative agreement with respect to experiment was reached for the out-of-resonance low-energy part of the β_HRS frequency dispersion. This work paves the way toward an accurate prediction of the SHG of large structures-a requirement for the design of new and improved SHIM biotags.

Two-Photon Absorption Activity of BOPHY Derivatives: Insights from Theory

We present a theoretical study of a two-photon absorption (2PA) process in dipolar and quadrupolar systems containing two BF₂ units. For this purpose, we considered 13 systems studied by Ponce-Vargas et al. [J. Phys. Chem. B2017, 121, 10850−1085829136383] and performed linear and quadratic response theory calculations based on the RI-CC2 method to obtain the 2PA parameters. Furthermore, using the recently developed generalized few-state model, we provided an in-depth view of the changes in 2PA properties in the molecules considered. Our results clearly indicate that suitable electron-donating group substitution to the core BF₂ units results in a large red-shift of the two-photon absorption wavelength, thereby entering into the desired biological window. Furthermore, the corresponding 2PA strength also increases significantly (up to 30-fold). This makes the substituted systems a potential candidate for biological imaging.

Optimization of Three State Conical Intersections by Adaptive Penalty Function Algorithm in Connection with the Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory Method (MRSF-TDDFT)

A new adaptive algorithm for penalty function optimization for minimum-energy three-states conical intersections (ME3CI) is suggested. The new algorithm differs from the original penalty function algorithm by (a) removing the redundancy in the target function, (b) using an adaptive increment for the penalty function weighting factor, and (c) using tighter convergence criteria for the energy gap. The latter was introduced to guarantee convergence to a true conical intersection rather than to a narrowly avoided crossing geometry. The new algorithm was tested in the optimization of the ME3CI geometries in butadiene and malonaldehyde, where all of the previously found true ME3CI geometries were recovered. The previously found butadiene's CI_3/2/1 turned out to be a narrowly avoided crossing. For butadiene, seven new ME3CI geometries have been located. Because of the removal of the redundancy and the use of the adaptive weighting factor, the convergence rate of the new algorithm is noticeably improved as compared to that of the previously proposed penalty function algorithm. The application to malonaldehyde and butadiene demonstrates that the three-state conical intersections may be more abundant and hence more involved in the photochemistry than previously thought. The recently developed mixed-reference spin flip (MRSF)-TDDFT method yields ME3CI geometries and relative energies quantitatively consistent with the previously reported calculations at a much reduced computational cost.

Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory (MRSF-TDDFT) as a Simple yet Accurate Method for Diradicals and Diradicaloids

Due to their multiconfigurational nature featuring strong electron correlation, accurate description of diradicals and diradicaloids is a challenge for quantum chemical methods. The recently developed mixed-reference spin-flip (MRSF)-TDDFT method is capable of describing the multiconfigurational electronic states of these systems while avoiding the spin-contamination pitfalls of SF-TDDFT. Here, we apply MRSF-TDDFT to study the adiabatic singlet-triplet (ST) gaps in a series of well-known diradicals and diradicaloids. On average, MRSF displays a very high prediction accuracy of the adiabatic ST gaps with the mean absolute error (MAE) amounting to 0.14 eV. In addition, MRSF is capable of accurately describing the effect of the Jahn-Teller distortion occurring in the trimethylenemethane diradical, the violation of the Hund rule in a series of the didehydrotoluene diradicals, and the potential energy surfaces of the didehydrobenzene (benzyne) diradicals. A convenient criterion for distinguishing diradicals and diradicaloids is suggested on the basis of the easily obtainable quantities. In all of these cases, which are difficult for the conventional methods of density functional theory (DFT), MRSF shows results consistent with the experiment and the high-level ab initio computations. Hence, the present study documents the reliability and accuracy of MRSF and lays out the guidelines for its application to strongly correlated molecular systems.

Predictions of High-Order Electric Properties of Molecules: Can We Benefit from Machine Learning?

There is an exigency of adopting machine learning techniques to screen and discover new materials which could address many societal and technological challenges. In this work, we follow this trend and employ machine learning to study (high-order) electric properties of organic compounds. The results of quantum-chemistry calculations of polarizability and first hyperpolarizability, obtained for more than 50,000 compounds, served as targets for machine learning-based predictions. The studied set of molecular structures encompasses organic push-pull molecules with variable linker lengths. Moreover, the diversified set of linkers, composed of alternating single/double and single/triple carbon-carbon bonds, was considered. This study demonstrates that the applied machine learning strategy allows us to obtain the correlation coefficients, between predicted and reference values of (hyper)polarizabilities, exceeding 0.9 on training, validation, and test set. However, in order to achieve such satisfactory predictive power, one needs to choose the training set appropriately, as the machine learning methods are very sensitive to the linker-type diversity in the training set, yielding catastrophic predictions in certain cases. Furthermore, the dependence of (hyper)polarizability on the length of spacers was studied in detail, allowing for explanation of the appreciably high accuracy of employed approaches.

Nonlinear-response properties in a simplified time-dependent density functional theory (sTD-DFT) framework: Evaluation of the first hyperpolarizability

Recent developments in nonlinear imaging microscopy show the need to implement new theoretical tools, which are able to characterize nonlinear optical properties in an efficient way. For second-harmonic imaging microscopy (SHIM), quantum chemistry could play an important role to design new exogenous dyes with enhanced first hyperpolarizabilities or to characterize the response origin in large endogenous biological systems. Such methods should be able to screen a large number of compounds while reproducing their trends and to treat large systems in reasonable computation times. To fulfill these requirements, we present a new simplified time-dependent density functional theory (sTD-DFT) implementation to evaluate the first hyperpolarizability where the Coulomb and exchange integrals are approximated by short-range damped Coulomb interactions of transition density monopoles. For an ultra-fast computation of the first hyperpolarizability, a tight-binding version (sTD-DFT-xTB) is also proposed. In our implementation, a sTD-DFT calculation is more than 600 time faster with respect to a regular TD-DFT treatment, while the xTB version speeds up the entire calculation further by at least two orders of magnitude. We challenge our implementation on three test cases: typical push-pull π-conjugated compounds, fluorescent proteins, and a collagen model, which were selected to model requirements for SHIM applications.