Puzzling problem (using virial theorem)

[Findg]

Consider a star of Radius R and mass M, with a pressure gradient given by

\frac{dP}{dr} = \frac{4\pi}{3}G\rho²r exp(-\frac{rr}{\lambda\lambda})

where \rho is the central density. calculate the gravitational energy, using the Virial theorem. Show that in the limit \lambda « R this energy is given by

E = \frac{RGMM}{3R\lambda}for tecnical reasons:
MM = M²
rr = r2
\lambda\lambda = \lambda²

 

[member 553083]

Using the Virial theorem, the total gravitational energy of a star is given by:E = -\frac{1}{2}\int \rho \Phi \ dVwhere \Phi is the gravitational potential. Since we are dealing with an isotropic pressure gradient, the gravitational potential is given by:\Phi = -4\pi G \int_0^R \frac{\rho r'^2 exp(-\frac{r'^2}{\lambda^2})}{r} dr'Integrating this expression and substituting in the expression for the gravitational energy, we obtain:E = - \frac{2\pi G R^5 \rho}{15\lambda^2}In the limit \lambda << R, we can approximate the integral as:E \approx - \frac{8\pi G R^5 \rho}{15\lambda^2}Finally, substituting the expression for the central density \rho = \frac{3M}{4\pi R^3}, we obtain:E \approx \frac{RGMM}{3R\lambda}