skip to main content

Title: How Good Is the Density-Corrected SCAN Functional for Neutral and Ionic Aqueous Systems, and What Is So Right about the Hartree–Fock Density?
Award ID(s):
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Journal of Chemical Theory and Computation
Page Range / eLocation ID:
4745 to 4761
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    A major challenge in density functional theory (DFT) is the development of density functional approximations (DFAs) to overcome errors in existing DFAs, leading to more complex functionals. For such functionals, we consider roles of the noninteracting reference systems. The electron density of the Kohn–Sham (KS) reference with a local potential has been traditionally defined as being equal to the electron density of the physical system. This key idea has been applied in two ways: the inverse calculation of such a local KS potential for the reference from a given density and the direct calculation of density and energy based on given DFAs. By construction, the inverse calculation can yield a KS reference with the density equal to the input density of the physical system. In application of DFT, however, it is the direct calculation of density and energy from a DFA that plays a central role. For direct calculations, we find that the self-consistent density of the KS reference defined by the optimized effective potential (OEP), is not the density of the physical system, when the DFA is dependent on the external potential. This inequality holds also for the density of generalized KS (GKS) or generalized OEP reference, which allows a nonlocal potential, when the DFA is dependent on the external potential. Instead, the density of the physical system, consistent with a given DFA, is given by the linear response of the total energy with respect to the variation of the external potential. This is a paradigm shift in DFT on the use of noninteracting references: the noninteracting KS or GKS references represent the explicit computational variables for energy minimization, but not the density of the physical system for external potential-dependent DFAs. We develop the expressions for the electron density so defined through the linear response for general DFAs, demonstrate the results for orbital functionals and for many-body perturbation theory within the second-order and the random-phase approximation, and explore the connections to developments in DFT. 
    more » « less
  2. ABSTRACT We study the evolutionary trend of the total density profile of early-type galaxies (ETGs) in IllustrisTNG. To this end, we trace ETGs from z = 0 to 4 and measure the power-law slope γ′ of the total density profile for their main progenitors. We find that their slopes γ′ steepen on average during z ∼ 4–2, then becoming shallower until z = 1, after which they remain almost constant, aside from a residual trend of becoming shallower towards z = 0. We also compare to a statistical sample of ETGs at different redshifts, selected based on their luminosity profiles and stellar masses. Due to different selection effects, the average slopes of the statistical samples follow a modified evolutionary trend. They monotonically decrease since z = 3, and after z ≈ 1, they remain nearly invariant with a mild increase towards z = 0. These evolutionary trends are mass dependent for both samples, with low-mass galaxies having in general steeper slopes than their more massive counterparts. Galaxies that transitioned to ETGs more recently have steeper mean slopes as they tend to be smaller and more compact at any given redshift. By analysing the impact of mergers and AGN feedback on the progenitors’ evolution, we conjecture a multiphase path leading to isothermality in ETGs: dissipation associated with rapid wet mergers tends to steepen γ′ from z = 4 to 2, whereas subsequent AGN feedback (especially in the kinetic mode) makes γ′ shallower again from z = 2 to 1. Afterwards, passive evolution from z = 1 to 0, mainly through gas-poor mergers, mildly decreases γ′ and maintains the overall mass distribution close to isothermal. 
    more » « less