As discussed in detail in our previous work, the bias towards HS imposed by the Hubbard term (that penalizes LS more than HS) would yield negative values of all the adiabatic energy differences computed in this work, regardless of the molecule and framework choice (see comparison between the density-corrected and the standard DFT+U method in Table and Table ) and larger bias are computed for stronger ligand field molecules.