Density functionals examined include SCAN, , SCAN0 (with 25% exact exchange), 90 B3LYP, ωB97X-V, Becke's "half-and-half" functional (BH&HLYP) containing 50% exact exchange, CAM-B3LYP with range-separation parameter ω = 0.33 bohr -1 , and the long-range corrected (LRC) functionals LRC-ωPBE (with ω = 0.3 bohr -1 ) and LRC-ωPBEh (ω = 0.2 bohr -1 ), the latter of which includes 20% exact exchange at short range. Finally, we examine the short-range corrected (SRC) functional SRC1-r1, which was parameterized for TD-DFT calculations of K-edge XAS transition energies for the elements C, O, N, and F. The def2-QZVP basis set is used for all calculations, in an effort to separate basis-set errors from methodological errors. Previous results for core-level binding energies indicate that DFT/def2-QZVP values are converged, such that uncontracting the core functions makes negligible difference. (For def2-TZVP, uncontracting the basis set changes electron binding energies by ∼ 0.4 eV. ) For generalized gradient approximations (GGAs) we use the SG-1 quadrature grid, 100 whereas SG-2 is used for meta-GGAs. For symmetry-equivalent atoms, Boys localization 102 is applied prior to the fractional-electron and ∆SCF calculations.