Benchmarking total energies with Hund's J terms in Hubbard-corrected spin-crossover chemistry

The effect of the Hund's J terms in various DFT+U+J corrections to semi-local spin-density functional theory is assessed for a series of four octahedrally coordinated Fe(II) spin-crossover molecules spanning the covalent end of the ligand field spectrum. We report values and analyze trends for the Hubbard U and Hund's J parameters determined via minimum-tracking linear response for all valence atomic subspaces and relevant spin states in these molecules. We then methodically apply them via simplified rotationally-invariant Hubbard functionals in search of the simplest combination to yield reliable adiabatic energy differences with respect to those obtained using CASPT2/CC. The observed failure of canonical, positively-signed Hund's J terms in furthering the already robust capacity of DFT+U to obtain accurate energetics prompts an evaluation of their limitations when seeking to account for the static correlation phenomena in such strongly covalent systems and suggests directions for their improvement.