Combining Local Range Separation and Local Hybrids: A Step in the Quest for Obtaining Good Energies and Eigenvalues from One Functional.

Some of the most successful exchange-correlation approximations in density functional theory are "hybrids", i.e., they rely on combining semilocal density functionals with exact nonlocal Fock exchange. In recent years, two classes of hybrid functionals have emerged as particularly promising: range-separated hybrids on the one hand, and local hybrids on the other hand. These functionals offer the hope to overcome a long-standing "observable dilemma", i.e., the fact that density functionals typically yield either a good description of binding energies, as obtained, e.g., in global and local hybrids, or physically interpretable eigenvalues, as obtained, e.g., in optimally tuned range-separated hybrids. Obtaining both of these characteristics from one and the same functional with the same set of parameters has been a long-standing challenge. We here discuss combining the concepts of local range separation and local hybrids as part of a constraint-guided quest for functionals that overcome the observable dilemma.

Predicting fundamental gaps accurately from density functional theory with non-empirical local range separation.

We present an exchange-correlation approximation in which the Coulomb interaction is split into long- and short-range components and the range separation is determined by a non-empirical density functional. The functional respects important constraints, such as the homogeneous and slowly varying density limits, leads to the correct long-range potential, and eliminates one-electron self-interaction. Our approach is designed for spectroscopic purposes and closely approximates the piecewise linearity of the energy as a function of the particle number. The functional's accuracy for predicting the fundamental gap in generalized Kohn-Sham theory is demonstrated for a large number of systems, including organic semiconductors with a notoriously difficult electronic structure.

Meta-generalized gradient approximations in time dependent generalized Kohn-Sham theory: Importance of the current density correction.

We revisit the use of Meta-Generalized Gradient Approximations (mGGAs) in time-dependent density functional theory, reviewing conceptual questions and solving the generalized Kohn-Sham equations by real-time propagation. After discussing the technical aspects of using mGGAs in combination with pseudopotentials and comparing real-space and basis set results, we focus on investigating the importance of the current-density based gauge invariance correction. For the two modern mGGAs that we investigate in this work, TASK and r2SCAN, we observe that for some systems, the current density correction leads to negligible changes, but for others, it changes excitation energies by up to 40% and more than 0.8 eV. In the cases that we study, the agreement with the reference data is improved by the current density correction.

Exact exchange-like electric response from a meta-generalized gradient approximation: A semilocal realization of ultranonlocality.

We review the concept of ultranonlocality in density functional theory and the relation between ultranonlocality, the derivative discontinuity of the exchange energy, and the static electric response in extended molecular systems. We present the construction of a new meta-generalized gradient approximation for exchange that captures the ultranonlocal response to a static electric field in very close correspondence to exact exchange, yet at a fraction of its computational cost. This functional, in particular, also captures the dependence of the response on the system size. The static electric polarizabilities of hydrogen chains and oligo-acetylene molecules calculated with this meta-GGA are quantitatively close to the ones obtained with exact exchange. The chances and challenges associated with the construction of meta-GGAs that are intended to combine a substantial derivative discontinuity and ultranonlocality with an accurate description of electronic binding are discussed.

Delocalized electronic excitations and their role in directional charge transfer in the reaction center of Rhodobacter sphaeroides.

In purple bacteria, the fundamental charge-separation step that drives the conversion of radiation energy into chemical energy proceeds along one branch-the A branch-of a heterodimeric pigment-protein complex, the reaction center. Here, we use first principles time-dependent density functional theory (TDDFT) with an optimally-tuned range-separated hybrid functional to investigate the electronic and excited-state structure of the six primary pigments in the reaction center of Rhodobacter sphaeroides. By explicitly including amino-acid residues surrounding these six pigments in our TDDFT calculations, we systematically study the effect of the protein environment on energy and charge-transfer excitations. Our calculations show that a forward charge transfer into the A branch is significantly lower in energy than the first charge transfer into the B branch, in agreement with the unidirectional charge transfer observed experimentally. We further show that the inclusion of the protein environment redshifts this excitation significantly, allowing for energy transfer from the coupled Qx excitations. Through analysis of transition and difference densities, we demonstrate that most of the Q-band excitations are strongly delocalized over several pigments and that both their spatial delocalization and charge-transfer character determine how strongly affected they are by thermally-activated molecular vibrations. Our results suggest a mechanism for charge-transfer in this bacterial reaction center and pave the way for further first-principles investigations of the interplay between delocalized excited states, vibronic coupling, and the role of the protein environment in this and other complex light-harvesting systems.

Understanding Primary Charge Separation in the Heliobacterial Reaction Center.

The homodimeric reaction center of heliobacteria retains features of the ancestral reaction center and can thus provide insights into the evolution of photosynthesis. Primary charge separation is expected to proceed in a two-step mechanism along either of the two reaction center branches. We reveal the first charge-separation step from first-principles calculations based on time-dependent density functional theory with an optimally tuned range-separated hybrid and ab initio Born-Oppenheimer molecular dynamics: the electron is most likely localized on the electron transfer cofactor 3 (EC3, OH-chlorophyll a), and the hole on the adjacent EC2. Including substantial parts of the surrounding protein environment into the calculations shows that a distinct structural mechanism is decisive for the relative energetic positioning of the electronic excitations: specific charged amino acids in the vicinity of EC3 lower the energy of charge-transfer excitations and thus facilitate efficient charge separation. These results are discussed considering recent experimental insights.

Hybrid functionals with local range separation: Accurate atomization energies and reaction barrier heights.

Range-separated hybrid approximations to the exchange-correlation density functional mix exact and semi-local exchange in a position-dependent manner. In their conventional form, the range separation is controlled by a constant parameter. Turning this constant into a density functional leads to a locally space-dependent range-separation function and thus a more powerful and flexible range-separation approach. In this work, we explore the self-consistent implementation of a local range-separated hybrid, taking into account a one-electron self-interaction correction and the behavior under uniform density scaling. We discuss different forms of the local range-separation function that depend on the electron density, its gradient, and the kinetic energy density. For test sets of atomization energies, reaction barrier heights, and total energies of atoms, we demonstrate that our best model is a clear improvement over common global range-separated hybrid functionals and can compete with density functionals that contain multiple empirical parameters. Promising results for equilibrium bond lengths, harmonic vibrational frequencies, and vertical ionization potentials further underline the potential and flexibility of our approach.

Investigating Primary Charge Separation in the Reaction Center of Heliobacterium modesticaldum.

We compute the primary charge separation step in the homodimeric reaction center (RC) of Heliobacterium modesticaldum from first principles. Using time-dependent density functional theory with the optimally tuned range-separated hybrid functional ωPBE, we calculate the excitations of a system comprising the special pair, the adjacent accessory bacteriochlorophylls, and the most relevant parts of the surrounding protein environment. The structure of the excitation spectrum can be rationalized from coupling of the individual bacteriochlorophyll pigments similar to molecular J- and H-aggregates. We find excited states corresponding to forward-charge transfer along the individual branches of the RC of H. modesticaldum. In the spectrum, these are located at an energy between the coupled Qy and Qx transitions. With ab initio Born-Oppenheimer molecular dynamics simulations, we reveal the influence of thermal vibrations on the excited states. The results show that the energy gap between the coupled Qy and the forward-charge transfer excitations is ∼0.4 eV, which we consider to conflict with the concept of a direct transfer mechanism. Our calculations, however, reveal a certain spectral overlap of the forward-charge transfer and the coupled Qx excitations. The reliability and robustness of the results are demonstrated by several numerical tests.

Self-interaction correction, electrostatic, and structural influences on time-dependent density functional theory excitations of bacteriochlorophylls from the light-harvesting complex 2.

First-principles calculations offer the chance to obtain a microscopic understanding of light-harvesting processes. Time-dependent density functional theory can have the computational efficiency to allow for such calculations. However, the (semi-)local exchange-correlation approximations that are computationally most efficient fail to describe charge-transfer excitations reliably. We here investigate whether the inexpensive average density self-interaction correction (ADSIC) remedies the problem. For the systems that we study, ADSIC is even more prone to the charge-transfer problem than the local density approximation. We further explore the recently reported finding that the electrostatic potential associated with the chromophores' protein environment in the light-harvesting complex 2 beneficially shifts spurious excitations. We find a great sensitivity on the chromophores' atomistic structure in this problem. Geometries obtained from classical molecular dynamics are more strongly affected by the spurious charge-transfer problem than the ones obtained from crystallography or density functional theory. For crystal structure geometries and density-functional theory optimized ones, our calculations confirm that the electrostatic potential shifts the spurious excitations out of the energetic range that is most relevant for electronic coupling.

Molecular excitations from meta-generalized gradient approximations in the Kohn-Sham scheme.

Meta-Generalized Gradient Approximations (meta-GGAs) can, in principle, include spatial and temporal nonlocality in time-dependent density functional theory at a much lower computational cost than functionals that use exact exchange. We here test whether a meta-GGA that has recently been developed with a focus on capturing nonlocal response properties and the particle number discontinuity can realize such features in practice. To this end, we extended the frequency-dependent Sternheimer formalism to the meta-GGA case. Using the Krieger-Li-Iafrate (KLI) approximation, we calculate the optical response for the selected paradigm molecular systems and compare the meta-GGA Kohn-Sham response to the one found with exact exchange and conventional (semi-)local functionals. We find that the new meta-GGA captures important properties of the nonlocal exchange response. The KLI approximation, however, emerges as a limiting factor in the evaluation of charge-transfer excitations.

Piecewise linearity, freedom from self-interaction, and a Coulomb asymptotic potential: three related yet inequivalent properties of the exact density functional.

The exact energy functional of density functional theory (DFT) is well known to obey various constraints. Three conditions that must be obeyed by the exact energy functional, but may or may not be obeyed by approximate ones, are often pointed out as important in general and for accurate computation of spectroscopic observables in particular. These are: (1) piecewise linearity as a function of the fractional particle number, (2) freedom from one-electron self-interaction, and (3) for a finite system, the functional derivative with respect to the density results in an asymptotic -1/r potential (in Hartree atomic units), where r is the distance from the system center. In this overview, we explain what these conditions are, what they address, and why each one is of importance for spectroscopy. We then show, using specific examples from the literature, that these three properties are related, but are not equivalent and need to be assessed individually.

Exploring local range separation: The role of spin scaling and one-electron self-interaction.

Range-separated hybrid functionals with a fitted or tuned global range-separation parameter are frequently used in density functional theory. We here explore the concept of local range separation, i.e., of turning the range-separation parameter into an explicit semilocal density functional. We impose three simple constraints on the local range-separation parameter that are frequently used in density functional construction: uniform density scaling, the homogeneous electron gas limit, and freedom from one-electron self-interaction. We further discuss different ways of how to model the spin dependence in combination with local range separation. We evaluate our local range-separation energy functionals exactly for closed-shell atoms using the previously suggested hypergeneralized gradient approximation for molecules and assess the quality of this approximation. We find a local range-separated hybrid functional that yields accurate binding energies for a set of small molecules.

Assessing density functional theory in real-time and real-space as a tool for studying bacteriochlorophylls and the light-harvesting complex 2.

We use real-time density functional theory on a real-space grid to calculate electronic excitations of bacteriochlorophyll chromophores of the light-harvesting complex 2 (LH2). Comparison with Gaussian basis set calculations allows us to assess the numerical trust range for computing electron dynamics in coupled chromophores with both types of techniques. Tuned range-separated hybrid calculations for one bacteriochlorophyll as well as two coupled ones are used as a reference against which we compare results from the adiabatic time-dependent local density approximation (TDLDA). The tuned range-separated hybrid calculations lead to a qualitatively correct description of the electronic excitations and couplings. They allow us to identify spurious charge-transfer excitations that are obtained with the TDLDA. When we take into account the environment that the LH2 protein complex forms for the bacteriochlorophylls, we find that it substantially shifts the energy of the spurious charge-transfer excitations, restoring a qualitatively correct electronic coupling of the dominant excitations also for TDLDA.

Linear response time-dependent density functional theory without unoccupied states: The Kohn-Sham-Sternheimer scheme revisited.

The Sternheimer approach to time-dependent density functional theory in the linear response regime is attractive because of its computational efficiency. The latter results from avoiding the explicit calculation of unoccupied orbitals and from the basic structure of the Sternheimer equations, which naturally lend themselves to far-reaching parallelization. In this article, we take a fresh look at the frequency-dependent Sternheimer equations. We first give a complete, self-contained derivation of the equations that complements previous derivations. We then discuss several aspects of an efficient numerical realization. As a worked example, we compute the photoabsorption spectra of small hydrogenated silicon clusters and confirm that for these the quality of the Kohn-Sham eigenvalues is more important than the effects of the exchange-correlation kernel. Finally, we demonstrate how triplet excitations can readily be computed from the Sternheimer approach.

Dielectric Screening Meets Optimally Tuned Density Functionals.

A short overview of recent attempts at merging two independently developed methods is presented. These are the optimal tuning of a range-separated hybrid (OT-RSH) functional, developed to provide an accurate first-principles description of the electronic structure and optical properties of gas-phase molecules, and the polarizable continuum model (PCM), developed to provide an approximate but computationally tractable description of a solvent in terms of an effective dielectric medium. After a brief overview of the OT-RSH approach, its combination with the PCM as a potentially accurate yet low-cost approach to the study of molecular assemblies and solids, particularly in the context of photocatalysis and photovoltaics, is discussed. First, solvated molecules are considered, with an emphasis on the challenge of balancing eigenvalue and total energy trends. Then, it is shown that the same merging of methods can also be used to study the electronic and optical properties of molecular solids, with a similar discussion of the pros and cons. Tuning of the effective scalar dielectric constant as one recent approach that mitigates some of the difficulties in merging the two approaches is considered.

Accurate Evaluation of Real-Time Density Functional Theory Providing Access to Challenging Electron Dynamics.

We demonstrate that electronic excitations and their transition densities can be obtained from real-time density functional theory calculations with great accuracy by relating the data from numerical propagation to the analytical form of the electronic response after a boost excitation. The latter is derived in this article. This approach facilitates quantitatively obtaining oscillator strengths, identifying excitations that carry very small oscillator strengths, and assessing electronic couplings from transition densities based on comparatively short propagation times. These features are of interest in particular for studying light-harvesting systems. For demonstration purposes, we study the Q band excitations of bacteriochlorophyll a (BChl a) and calculate coupling strengths between two BChl a's to check the validity of the dipole-dipole and pure Coulomb coupling mechanisms. For further illustration, we investigate the paradigm test system Na4 and the coupling between two Na2 dimers.

Investigating the electronic structure of a supported metal nanoparticle: Pd in SiCN.

We investigate the electronic structure of a Palladium nanoparticle that is partially embedded in a matrix of silicon carbonitride. From classical molecular dynamics simulations we first obtain a representative atomic structure. This geometry then serves as input to density-functional theory calculations that allow us to access the electronic structure of the combined system of particle and matrix. In order to make the computations feasible, we devise a subsystem strategy for calculating the relevant electronic properties. We analyze the Kohn-Sham density of states and pay particular attention to d-states which are prone to be affected by electronic self-interaction. We find that the density of states close to the Fermi level is dominated by states that originate from the Palladium nanoparticle. The matrix has little direct effect on the electronic structure of the metal. Our results contribute to explaining why silicon carbonitride does not have detrimental effects on the catalytic properties of palladium particles and can serve positively as a stabilizing mechanical support.

Hydrogen binding energies and electronic structure of Ni-Pd particles: a clue to their special catalytic properties.

We investigate the electronic structure of nickel-palladium systems with first-principles density functional theory (DFT). Our study is motivated by the experimental observation that nickel-palladium nanoalloys containing approximately equal amounts of nickel and palladium show a higher catalytic activity in hydrogenation reactions than pure particles of either metal. Our calculations show that the energy with which hydrogen is bound to the metal surface, a decisive factor in catalytic activity, strongly depends on the nickel-palladium ratio. The weakest binding and Gibbs free energies of hydrogen adsorption close to zero are found for systems with roughly equal amounts of nickel and palladium. While for extended Ni-Pd surfaces this observation can be explained by the well-established d-band model, for Ni-Pd clusters analysis of the electronic density of states shows a complex interplay of s- and d-orbital contributions. Studying extended surfaces further reveals that the formation of nanostructures on the surface can influence hydrogen adsorption pronouncedly. We discuss the important role that the exchange-correlation functional approximation plays in DFT calculations for nickel-palladium systems.

Effect of ensemble generalization on the highest-occupied Kohn-Sham eigenvalue.

There are several approximations to the exchange-correlation functional in density-functional theory, which accurately predict total energy-related properties of many-electron systems, such as binding energies, bond lengths, and crystal structures. Other approximations are designed to describe potential-related processes, such as charge transfer and photoemission. However, the development of a functional which can serve the two purposes simultaneously is a long-standing challenge. Trying to address it, we employ in the current work the ensemble generalization procedure proposed by Kraisler and Kronik [Phys. Rev. Lett. 110, 126403 (2013)]. Focusing on the prediction of the ionization potential via the highest occupied Kohn-Sham eigenvalue, we examine a variety of exchange-correlation approximations: the local spin-density approximation, semi-local generalized gradient approximations, and global and local hybrid functionals. Results for a test set of 26 diatomic molecules and single atoms are presented. We find that the aforementioned ensemble generalization systematically improves the prediction of the ionization potential, for various systems and exchange-correlation functionals, without compromising the accuracy of total energy-related properties. We specifically examine hybrid functionals. These depend on a parameter controlling the ratio of semi-local to non-local functional components. The ionization potential obtained with ensemble-generalized functionals is found to depend only weakly on the parameter value, contrary to common experience with non-generalized hybrids, thus eliminating one aspect of the so-called "parameter dilemma" of hybrid functionals.