where the a, b, and c parameters are fitted to the QM calculations. Here a should be positive to obtain a stable system (minima as opposed to maxima at n = n0) and n0 is the number of explicit electrons for a neutral system (explicit because we utilize pseudopotentials). The quadratic form of F(n) is strictly verified in our calculations.