Derive and Solve the Lane-Emden Equation for a Polytropic Gas Sphere

[MatinSAR]

Homework Statement

1. Assume a polytropic equation of state for the gas: $P=K \rho ^{\gamma}$
2. Using the mechanical part of the stellar equation (mass conservation and
hydrostatic equilibrium) , drive the below equation. $$ \frac{\gamma K}{r^2} \frac{d}{dr} \left[ r^2 \rho^{\gamma - 2} \frac{d\rho}{dr} \right] = -4 \pi G \rho $$ Using the specified change of variables, derive the Lane-Emden equation and solve it for the given boundary conditions when $𝑛 = 0$. Then, interpret the resulting solution.

Relevant Equations

Mass conservation and hydrostatic equilibrium equation: $$\frac{dM(r)}{dr} = 4\pi r^2 \rho(r)$$ $$\frac{dP(r)}{dr} = -\frac{GM(r)\rho(r)}{r^2} $$

We start with hydrostatic equilibrium: $$ \frac{dP(r)}{dr} = K \gamma \rho ^ {\gamma -1} \dfrac {d \rho}{dr} =-\frac{GM(r)\rho(r)}{r^2}$$ $$ K \gamma \rho ^ {\gamma -1} \dfrac {d \rho}{dr} = - \dfrac {G \rho (r) }{r^2} \int 4 \pi r^2 \rho (r) \, dr$$ $$ r^2 K \gamma \rho ^ {\gamma -2} \dfrac {d \rho}{dr} = -G \int 4 \pi r^2 \rho (r)$$ $$\dfrac { \gamma K}{r^2} \dfrac {d}{dr} \left[ r^2 \rho ^ {\gamma -2} \dfrac {d \rho}{dr} \right] = -4 \pi G \rho$$ I think my answer to part 2 is correct. I'm not sure about following parts of the question:
3. Set $\gamma \equiv \frac{(n+1)}{n}$ and change the variables as $$ \rho(r) \equiv \rho_c [D_n(r)]^n, \quad \text{where} \quad 0 \leq D_n \leq 1 \quad r \equiv \lambda_n \xi \quad \lambda_n \equiv \left[ (n+1) \left( \frac{K \rho_c^{(1-n)/n}}{4 \pi G} \right) \right]^{1/2}$$ 4. Now find the Lane-Emden equation $$\frac{1}{\xi^2} \frac{d}{d\xi} \left[ \xi^2 \frac{d D_n}{d\xi} \right] = -D_n^n$$ 5.Assuming proper boundary conditions solve the Lane-Emden equation analytically for n=0.
7.Can you interpret the results? What the results imply about the size of the star?

I thnik should use the change of variables in part 3 to obtain the equation in part 4, then set $n=0$ and solve it. Then, use the above boundary conditions to find $D(\xi)$. If I were right until now, I can use $r \equiv \lambda_n \xi$ to talk about the size of the star. Am I right?

 

[pasmith]

I thnik should use the change of variables in part 3 to obtain the equation in part 4, then set $n=0$ and solve it. Then, use the above boundary conditions to find $D(\xi)$.

That is exactly what the question is asking you to do. Note that the actual boundary conditions at the origin are $D_n(0) = 1$ and $D_n'(0) = 0$; $\xi_1$ is then the first positive zero of $D_n$.

Without finding $D_0$, what do the definitions $$\rho(r) = \rho_c [D_n(r)]^n$$ and $$\lambda_n \propto (n+1)^{1/2}\rho_c^{(1-n)/(2n)}$$ suggest about the case $n = 0$ for $\rho_c > 0$?