Earth Mover's Distance as a Metric to Evaluate the Extent of Charge Transfer in Excitations Using Discretized Real-Space Densities

Importance of Orbital Invariance in Quantifying Electron-Hole Separation and Exciton Size

A fundamental tenet of quantum mechanics is that properties should be independent of representation. In self-consistent field methods such as density functional theory, this manifests as a requirement that properties be invariant with respect to unitary transformations of the occupied molecular orbitals and (separately) the unoccupied molecular orbitals. Various ad hoc measures of excited-state charge separation that are commonly used to analyze time-dependent density-functional calculations violate this requirement, as they are based on incoherent averages of excitation amplitudes rather than expectation values involving coherent superpositions. As a result, these metrics afford markedly different values in various common representations, including canonical molecular orbitals, Boys-localized orbitals, and natural orbitals. Numerical values can be unstable with respect to basis-set expansion and may afford nonsensical results in the presence of extremely diffuse basis functions. In contrast, metrics based on well-defined expectation values are stable, representation-invariant, and physically interpretable. Use of natural transition orbitals improves the stability of the incoherent averages, but numerical values can only be interpreted as expectation value in the absence of superposition. To satisfy this condition, the particle and hole density matrices must each be dominated by a single eigenvector so that the transition density is well described by a single pair of natural transition orbitals. Counterexamples are readily found where this is not the case. Our results suggest that ad hoc charge-transfer diagnostics should be replaced by rigorous expectation values, which are no more expensive to compute.

Benchmarking DFT-based excited-state methods for intermolecular charge-transfer excitations

Intermolecular charge-transfer is a highly important process in biology and energy-conversion applications where generated charges need to be transported over several moieties. However, its theoretical description is challenging since the high accuracy required to describe these excited states must be accessible for calculations on large molecular systems. In this benchmark study, we identify reliable low-scaling computational methods for this task. Our reference results were obtained from highly accurate wavefunction calculations that restrict the size of the benchmark systems. However, the density-functional theory based methods that we identify as accurate can be applied to much larger systems. Since targeting charge-transfer states requires the unambiguous classification of an excited state, we first analyze several charge-transfer descriptors for their reliability concerning intermolecular charge-transfer and single out the charge-transfer distance calculated based on the variation of electron density upon excitation (DCT) as an optimal choice for our purposes. In general, best results are obtained for orbital-optimized methods and among those, the maximum overlap method proved to be the most numerically stable variant when using the initial MOs as reference orbitals. Favorable error cancellation with optimally-tuned range-separated hybrid functionals and a rather small basis set can provide an economical yet reasonable wavefunction when using time-dependent density functional theory, which provides relevant information about the excited-state character to be used in the orbital-optimized methods. The qualitative agreement makes these fast calculations attractive for high-throughput screening applications.

Optimal Transport Distances to Characterize Electronic Excitations

Understanding the character of electronic excitations is important in computational and reaction mechanistic studies, but their classification from simulations remains an open problem. Distances based on optimal transport have proven very useful in a plethora of classification problems and, therefore, seem a natural tool to try to tackle this challenge. We propose and investigate a new diagnostic Θ based on the Sinkhorn divergence from optimal transport. We evaluate a k-NN classification algorithm on Θ, the popular Λ diagnostic, and their combination, and assess their performance in labeling excitations, finding that (i) the combination only slightly improves the classification, (ii) Rydberg excitations are not separated well in any setting, and (iii) Θ breaks down for charge transfer in small molecules. We then define a length-scale-normalized version of Θ and show that the result correlates closely with Λ for results obtained with Gaussian basis functions. Finally, we discuss the orbital dependence of our approach and explore an orbital-independent version. Using an optimized combination of the optimal transport and overlap diagnostics together with a different metric is in our opinion the most promising for future classification studies.

Orbital-Optimized Versus Time-Dependent Density Functional Calculations of Intramolecular Charge Transfer Excited States

The performance of time-independent, orbital-optimized calculations of excited states is assessed with respect to charge transfer excitations in organic molecules in comparison to the linear-response time-dependent density functional theory (TD-DFT) approach. A direct optimization method to converge on saddle points of the electronic energy surface is used to carry out calculations with the local density approximation (LDA) and the generalized gradient approximation (GGA) functionals PBE and BLYP for a set of 27 excitations in 15 molecules. The time-independent approach is fully variational and provides a relaxed excited state electron density from which the extent of charge transfer is quantified. The TD-DFT calculations are generally found to provide larger charge transfer distances compared to the orbital-optimized calculations, even when including orbital relaxation effects with the Z-vector method. While the error on the excitation energy relative to theoretical best estimates is found to increase with the extent of charge transfer up to ca. -2 eV for TD-DFT, no correlation is observed for the orbital-optimized approach. The orbital-optimized calculations with the LDA and the GGA functionals provide a mean absolute error of ∼0.7 eV, outperforming TD-DFT with both local and global hybrid functionals for excitations with a long-range charge transfer character. Orbital-optimized calculations with the global hybrid functional B3LYP and the range-separated hybrid functional CAM-B3LYP on a selection of states with short- and long-range charge transfer indicate that inclusion of exact exchange has a small effect on the charge transfer distance, while it significantly improves the excitation energy, with the best-performing functional CAM-B3LYP providing an absolute error typically around 0.15 eV.

Visualizing and characterizing excited states from time-dependent density functional theory

Time-dependent density functional theory (TD-DFT) is the most widely-used electronic structure method for excited states, due to a favorable combination of low cost and semi-quantitative accuracy in many contexts, even if there are well recognized limitations. This Perspective describes various ways in which excited states from TD-DFT calculations can be visualized and analyzed, both qualitatively and quantitatively. This includes not just orbitals and densities but also well-defined statistical measures of electron-hole separation and of Frenkel-type exciton delocalization. Emphasis is placed on mathematical connections between methods that have often been discussed separately. Particular attention is paid to charge-transfer diagnostics, which provide indicators of when TD-DFT may not be trustworthy due to its categorical failure to describe long-range electron transfer. Measures of exciton size and charge separation that are directly connected to the underlying transition density are recommended over more ad hoc metrics for quantifying charge-transfer character.