“…As the Landau levels are completely destroyed, we see that even at this order in the approximation the magnetization is completely different from the expression (35). Actually, even Formulas (31) and (32) do not hold anymore, for such formulas were derived by assuming the linear dependence (32) of mi2 n,s on the integer n. As expression (24) of E n, ,s in terms of n displays no such linearity, and since that expression (24) is itself only an approximation for the highly nonlinear full expression (23), we conclude that the Landau-levels-based formalism for computing the magnetization M of neutron stars does not generally apply to the charged plasma in the core of highly dense neutron stars. The condition for the applicability of the formalism is to have the quartic anharmonic oscillator potential term in Equation ( 22) negligible, i.e., to have λ |e|B/h or to apply the formalism only to fermions close to the vicinity of the star's center.…”