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.