In this case, the density and kinetic energy densities read, for a given (spin,isospin) value, in agreement with [1]}, as
\(\rho _{q\sigma }(r)=\sum _{nl} \frac{2l+1}{4\pi }\varphi _{nl}^2(r)\) 
\(\tau _{q\sigma }(r)=\sum _{nl} \frac{2l+1}{4\pi }\left[\left(\partial _r\varphi _{nl}(r)\right)^2+\frac{l(l+1)}{r^2}\varphi _{nl}(r)^2\right]\) 
