Effect of the Gradient of the Spin-Polarization in Density Functional Approximations

The construction of non-empirical density functional approximations is typically guided by the satisfaction of exact constraints. An important constraint is the recovery of the gradient expansion for slowly varying electron densities. In prior constructions of semilocal density functional approximations, the $\nabla \zeta$-dependent terms in the gradient expansion of the correlation have been dropped, where $\zeta$ is the relative spin polarization. We propose a scheme by which such terms can be reintroduced into already constructed functionals without significantly affecting other constraints and norms. We implement this scheme on the Strongly Constrained and Appropriately Normed (SCAN) functional to construct a $\nabla \zeta$-corrected version of SCAN. The resulting functional is shown to provide improvements in transition-metal atoms and molecules without significantly affecting SCAN's accurate description of $sp$-systems. For the binding energy curve of the chromium dimer Cr$_2$, the SCAN underbinding is fully corected at large bond lengths and reduced at short bond lengths.