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}