The calculations use the full potential, relativistic, linear augmented plane wave method with the LSDA+U approach for consideration of the local 4f electronic orbitals.
With full SOC the BR ranges from 1.5 to 6.0 in going from Gd to Er.