An official website of the United States government
In this work, we have implemented the time-dependent density functional theory approximate auxiliary s function (TDDFT-aas) method, which is an approximate TDDFT method. Instead of calculating the exact two-center electron integrals in the K coupling matrix when solving the Casida equation, we approximate the integrals, thereby reducing the computational cost. In contrast to the related TDDFT plus tight-binding (TDDFT+TB) method, a new type of gamma function is used in the coupling matrix that does not depend on the tight-binding parameters. The calculated absorption spectra of silver and gold nanoparticles using TDDFT-aas show good agreement with TDDFT and TDDFT+TB results. In addition, we have implemented the analytical excited-state gradients for the TDDFT-aas method, which makes it possible to calculate the emission energy of molecular systems.
more »
« less
Analytical excited state gradients for time-dependent density functional theory plus tight binding (TDDFT + TB)
Understanding photoluminescent mechanisms has become essential for photocatalytic, biological, and electronic applications. Unfortunately, analyzing excited state potential energy surfaces (PESs) in large systems is computationally expensive, and hence limited with electronic structure methods such as time-dependent density functional theory (TDDFT). Inspired by the sTDDFT and sTDA methods, time-dependent density functional theory plus tight binding (TDDFT + TB) has been shown to reproduce linear response TDDFT results much faster than TDDFT, particularly in large nanoparticles. For photochemical processes, however, methods must go beyond the calculation of excitation energies. Herein, this work outlines an analytical approach to obtain the derivative of the vertical excitation energy in TDDFT + TB for more efficient excited state PES exploration. The gradient derivation is based on the Z vector method, which utilizes an auxiliary Lagrangian to characterize the excitation energy. The gradient is obtained when the derivatives of the Fock matrix, the coupling matrix, and the overlap matrix are all plugged into the auxiliary Lagrangian, and the Lagrange multipliers are solved. This article outlines the derivation of the analytical gradient, discusses the implementation in Amsterdam Modeling Suite, and provides proof of concept by analyzing the emission energy and optimized excited state geometry calculated by TDDFT and TDDFT + TB for small organic molecules and noble metal nanoclusters.
more »
« less
- Award ID(s):
- 1955823
- PAR ID:
- 10469966
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 22
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The Hohenberg-Kohn theorem of density-functional theory establishes the existence of a bijection between the ground-state electron density and the external potential of a many-body system. This guarantees a one-to-one map from the electron density to all observables of interest including electronic excited-state energies. Time-Dependent Density-Functional Theory (TDDFT) provides one framework to resolve this map; however, the approximations inherent in practical TDDFT calculations, together with their computational expense, motivate finding a cheaper, more direct map for electronic excitations. Here, we show that determining density and energy functionals via machine learning allows the equations of TDDFT to be bypassed. The framework we introduce is used to perform the first excited-state molecular dynamics simulations with a machine-learned functional on malonaldehyde and correctly capture the kinetics of its excited-state intramolecular proton transfer, allowing insight into how mechanical constraints can be used to control the proton transfer reaction in this molecule. This development opens the door to using machine-learned functionals for highly efficient excited-state dynamics simulations.more » « less
-
Following the formulation of cavity quantum-electrodynamical time-dependent density functional theory (cQED-TDDFT) models [Flick et al., ACS Photonics 6, 2757–2778 (2019) and Yang et al., J. Chem. Phys. 155, 064107 (2021)], here, we report the derivation and implementation of the analytic energy gradient for polaritonic states of a single photochrome within the cQED-TDDFT models. Such gradient evaluation is also applicable to a complex of explicitly specified photochromes or, with proper scaling, a set of parallel-oriented, identical-geometry, and non-interacting molecules in the microcavity.more » « less
-
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–Ag 20 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.more » « less
-
The QMol-grid package provides a suite of routines for performing quantum-mechanical simulations in atomic and molecular systems, currently implemented in one spatial dimension. It supports ground- and excited-state calculations for the Schrödinger equation, density-functional theory, and Hartree–Fock levels of theory as well as propagators for field-free and field-driven time-dependent Schrödinger equation (TDSE) and real-time time-dependent density-functional theory (TDDFT), using symplectic-split schemes. The package is written using MATLAB’s object-oriented features and handle classes. It is designed to facilitate access to the wave function(s) (TDSE) and the Kohn–Sham orbitals (TDDFT) within MATLAB’s environment.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0142240
