Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

Abstract

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low-energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. Intensity borrowing by these spurious states can affect intensities of the valence excitations, altering electronic bandshapes. At higher excitation energies, it is difficult to distinguish spurious charge-transfer states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces ``natural charge-transfer orbitals'' that provide a means to isolate orbitals that are most likely to participate in a CTTS excitation. Projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation and be used as an automated subspace selection algorithm for projection-based embedding of a high-level description of the CTTS state in a lower-level description of its environment. We apply this method to an ab initio molecular dynamics trajectory of I -(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with experiment and only one-third of the water molecules in the I -(H 2 O) 96 simulation cell need to be described with LR-TDDFT in order to obtain excitation energies that are converged to <0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments.File list (2)download file view on ChemRxiv NCTA.pdf (10.09 MiB) download file view on ChemRxiv NCTA-SI.pdf (3.58 MiB)

Hidden Hemibonding in the Aqueous Hydroxyl Radical

“…The use of LRC functionals in aqueous-phase TD-DFT calculations is necessary to avoid contamination of the spectrum by a plethora of spurious, low-energy charge-transfer states. − The presence of such states adds considerably to the cost of the calculation but can also impact line shapes, via spurious intensity borrowing . Functionals such as BH&HLYP with a larger-than-normal fraction of exact exchange have also been used to address this problem, but they tend to overestimate valence excitation energies, pushing them closer to Hartree–Fock values that are typically overestimated by 0.8 eV or more .…”

The Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory

“…Recently, Herbert and co-workers have used PbE in combination with TDDFT for the determination of the excited state in solution including selective solvent molecules by determining an effective CT acceptor space by projection . Significantly contributing solvent molecules within this subspace are then included in the active subsystem in the subsequent PbE-sTDDFT calculation.…”

Comprehensive Basis-Set Testing of Extended Symmetry-Adapted Perturbation Theory and Assessment of Mixed-Basis Combinations to Reduce Cost

“…(The GDD-tuned values of ω for each of the systems considered in this work can be found in Tables S1–S8 of the Supporting Information.) Previous work has shown that tuned values of ω are sensitive to the fraction of short-range exact exchange, but for a given density functional, these values are rather insensitive to the choice of basis set . The SG-1 quadrature grid is used for LRC-ωPBE.…”