“…An important observation from the results of Table 2 is that, contrarily to B3LYP, the BP86 functional leads to very high energy differences between the HS and BS states. As stressed by previous authors [32-36, 41-46, 48-55], the resulting values of the exchange coupling constant J, in the GGA case, are systematically larger in absolute values, than those expected for such systems [62,63,77,78]. We noticed the same overestimation using another GGA functional.…”