FAIMP3 strategy still underestimates the energy gaps, showing that also when the densities of the metal and the ligands are iteratively optimized to each others the pure electrostatic approach by FAIMP model is not good enough for accurate predictions. On the other hand, NOGF3 method, where the density of the ligands is adjusted to the metal density by the iterative procedure, shows crystal field gaps almost coincident with the CAHF/CASCI-SO results for the majority of the Tb(III), Dy(III), and Er(III) complexes in all the range of frequencies. Slightly bigger discrepancies are shown by no 3 -paah ligands, where an average overestimation of the energy gaps of about +40% and +17% is found for Tb(III) and Dy(III) ions, respectively, and an underestimation of about 26% for Er(III) metal. Differently, NOGF3 strategy do not lead to a significant improvement in the result for Ho(III) complexes.