Symmetry Breaking with the SCAN Density Functional Describes Strong Correlation in the Singlet Carbon Dimer.

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C2 singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B2, C2, O2, and F2 in the local spin density approximation (LSDA), the Perdew-Burke-Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C2. This work shows that spin symmetry breaking in singlet C2, which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.

Spherical vs non-spherical and symmetry-preserving vs symmetry-breaking densities of open-shell atoms in density functional theory.

The atomization energies of molecules from first-principles density functional approximations improve from the local spin-density approximation to the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) to the strongly constrained and appropriately normed (SCAN) meta-GGA, and their sensitivities to non-spherical components of the density increase in the same order. Thus, these functional advances increase density sensitivity and imitate the exact constrained search over correlated wavefunctions better than that over ensembles. The diatomic molecules studied here, singlet C2 and F2 plus triplet B2 and O2, have cylindrically symmetric densities. Because the densities of the corresponding atoms are non-spherical, the approximate Kohn-Sham potentials for the atoms have a lower symmetry than that of the external (nuclear) potential so that the non-interacting wavefunctions are not eigenstates of the square of total orbital angular momentum, breaking a symmetry that yields a feature of the exact ground-state density. That spatial symmetry can be preserved by a non-self-consistent approach in which a self-consistent equilibrium-ensemble calculation is followed by integer re-occupation of the Kohn-Sham orbitals as the first of several steps. The symmetry-preserving approach is different from symmetry restoration based on projection. First-step space- (and space-spin-) symmetry preservation in atoms is shown to have a small effect on the atomization energies of molecules, quantifying earlier observations by Fertig and Kohn. Thus, the standard Kohn-Sham way of calculating atomization energies, with self-consistent symmetry breaking to minimize the energy, is justified at least for the common cases where the molecules cannot break symmetry. Unless symmetry breaking is allowed in the molecule, SCAN strongly underestimates the atomization energy of strongly correlated singlet C2.

Simple hydrogenic estimates for the exchange and correlation energies of atoms and atomic ions, with implications for density functional theory.

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-Z expansions for the corresponding energies of neutral atoms with atomic number Z and electron number N = Z, which are correct to the leading order (-0.221Z5/3 and -0.021Z ln Z, respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to Ex(N, Z) ≈ -0.354N2/3Z (as known before only for Z ≫ N ≫ 1) and Ec ≈ -0.02N ln N. These asymptotic estimates are most correct for atomic ions with large N and Z ≫ N, but we find that they are qualitatively and semi-quantitatively correct even for small N and N ≈ Z. The large-N asymptotic behavior of the energy is pre-figured in small-N atoms and atomic ions, supporting the argument that widely predictive approximate density functionals should be designed to recover the correct asymptotics. It is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to Z in the Z → ∞ limit for any fixed N.