Extension of the element parameter set for ultra-fast excitation spectra calculation (sTDA-xTB)

Photoinduced Nonadiabatic Dynamics of a Single-Walled Carbon Nanotube-Porphyrin Complex

Single-walled carbon nanotubes (SWCNTs) have gained a lot of attention in the past few decades due to their promising optoelectronic properties. In addition, SWCNTs can form complexes that have good chemical stability and transport properties with other optical functional materials through noncovalent interactions. Elucidating the detailed mechanism of these complexes is of great significance for improving their optoelectronic properties. Nevertheless, simulating the photoinduced dynamics of these complexes accurately is rather challenging since they usually contain hundreds of atoms. To save computational efforts, most of the previous works have ignored the excitonic effects by employing nonadiabatic carrier (electron and hole) dynamics simulations. To properly consider the influence of excitonic effects on the photoinduced ultrafast processes of the SWCNT-tetraphenyl porphyrin (H2TPP) complex and to further improve the computational efficiency, we developed the nonadiabatic molecular dynamics (NAMD) method based on the extended tight binding-based simplified Tamm-Dancoff approximation (sTDA-xTB), which is applied to study the ultrafast photoinduced dynamics of the noncovalent SWCNT-porphyrin complex. In combination with statically electronic structure calculations, the present work successfully reveals the detailed microscopic mechanism of the ultrafast excitation energy transfer process of the complex. Upon local excitation on the H2TPP molecule, an ultrafast energy transfer process occurs from H2TPP (SWCNT-H2TPP*) to SWCNT (SWCNT*-H2TPP) within 10 fs. Then, two slower processes corresponding to the energy transfer from H2TPP to SWCNT and hole transfer from H2TPP to SWCNT take place in the 1 ps time scale. The sTDA-xTB-based electronic structure calculation and NAMD simulation results not only match the previous experimental observations from static and transient spectra but also provide more insights into the detailed information on the complex's photoinduced dynamics. Therefore, the sTDA-xTB-based NAMD method is a powerful theoretical tool for studying the ultrafast photoinduced dynamics in large extended systems with a large number of electronically excited states, which could be helpful for the subsequent design of SWCNT-based functional materials.

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.

Excited-State Properties for Extended Systems: Efficient Hybrid Density Functional Methods

Time-dependent density functional theory has become state-of-the-art for describing photophysical and photochemical processes in extended materials because of its affordable cost. The inclusion of exact exchange was shown to be essential for the correct description of the long-range asymptotics of electronic interactions and thus a well-balanced description of valence, Rydberg, and charge-transfer excitations. Several approaches for an efficient treatment of exact exchange have been established for the ground state, while implementations for excited-state properties are rare. Furthermore, the high computational costs required for excited-state properties in comparison to ground-state computations often hinder large-scale applications on periodic systems with hybrid functional accuracy. We therefore propose two approximate schemes for improving computational efficiency for the treatment of exact exchange. Within the auxiliary density matrix method (ADMM), exact exchange is estimated using a relatively small auxiliary basis and the introduced basis set incompleteness error is compensated by an exchange density functional correction term. Benchmark results for a test set of 35 molecules demonstrate that the mean absolute error introduced by ADMM is smaller than 0.3 pm for excited-state bond lengths and in the range of 0.02-0.04 eV for vertical excitation, adiabatic excitation, and fluorescence energies. Computational timings for a series of covalent-organic frameworks demonstrate that a speed-up of at least 1 order of magnitude can be achieved for excited-state geometry optimizations in comparison to conventional hybrid functionals. The second method is to use a semiempirical tight binding approximation for both Coulomb and exchange contributions to the excited-state kernel. This simplified Tamm-Dancoff approximation (sTDA) achieves an accuracy comparable to approximated hybrid density functional theory when referring to highly accurate coupled-cluster reference data. We find that excited-state bond lengths deviate by 1.1 pm on average and mean absolute errors in vertical excitation, adiabatic excitation, and fluorescence energies are in the range of 0.2-0.5 eV. In comparison to ADMM-approximated hybrid functional theory, sTDA accelerates the computation of broad-band excitation spectra by 1 order of magnitude, suggesting its potential use for large-scale screening purposes.

Perspective on Simplified Quantum Chemistry Methods for Excited States and Response Properties

We review recent developments in the framework of simplified quantum chemistry for excited state and optical response properties (sTD-DFT) and present future challenges for new method developments to improve accuracy and extend the range of application. In recent years, the scope of sTD-DFT was extended to molecular response calculations of the polarizability, optical rotation, first hyperpolarizability, two-photon absorption (2PA), and excited-state absorption for large systems with hundreds to thousands of atoms. The recently introduced spin-flip simplified time-dependent density functional theory (SF-sTD-DFT) variant enables an ultrafast treatment for diradicals and related strongly correlated systems. A few drawbacks were also identified, specifically for the computation of 2PA cross sections. We propose solutions to this problem and how to generally improve the accuracy of simplified schemes. New possible simplified schemes are also introduced for strongly correlated systems, e.g., with a second-order perturbative correlation correction. Interpretation tools that can extract chemical structure-property relationships from excited state or response calculations are also discussed. In particular, the recently introduced method-agnostic RespA approach based on natural response orbitals (NROs) as the key concept is employed.