Problem 17 - An Introduction to Modern Astrophysics

Question

Starting with the Lane-Emden equation and imposing the necessary boundary conditions, prove that the \(n=0\) polytrope has a solution given by $$D_{0}(\xi)=1-\frac{\xi^{2}}{6}, \quad \text { with } \xi_{1}=\sqrt{6}$$

Step-by-Step Solution

Answer

The solution for the \(n=0\) polytrope is \(D_{0}(\xi)=1-\frac{\xi^{2}}{6}\) with first zero at \(\xi_1=\sqrt{6}\).

Step 1: Understand the Lane-Emden Equation

The Lane-Emden equation is a dimensionless form of the Poisson's equation for the gravitational potential of a spherically symmetric polytropic self-gravitating fluid. It is given by: \[ \frac{1}{\xi^2}\frac{d}{d\xi}\left(\xi^2 \frac{d \theta}{d \xi}\right) + \theta^n = 0 \] where \(\theta\) is the dimensionless density, \(\xi\) is the dimensionless radial coordinate, and \(n\) is the polytropic index. For \(n=0\), the equation simplifies.

Step 2: Simplify the Lane-Emden Equation for n=0

Substituting \(n=0\) into the equation, \(\theta^0=1\). The Lane-Emden equation simplifies to: \[ \frac{1}{\xi^2}\frac{d}{d\xi}\left(\xi^2 \frac{d \theta}{d \xi}\right) + 1 = 0 \] The equation transforms into a second-order ordinary differential equation (ODE).

Step 3: Apply the Boundary Conditions

For a polytropic sphere, the boundary conditions typically are: \(\theta(0) = 1\) for the central density, and \( \frac{d\theta}{d\xi}(0) = 0 \) for symmetry, ensuring a smooth curve at the center. These conditions help solve the ODE.

Step 4: Solve the ODE for n=0

Integrate the ODE by first integrating \( \frac{d}{d\xi}\left(\xi^2 \frac{d \theta}{d \xi}\right) = -\xi^2 \) with regards to \(\xi\), giving: \[ \xi^2 \frac{d \theta}{d \xi} = -\frac{\xi^3}{3} + C \] Using \( \frac{d\theta}{d\xi}(0) = 0 \), find that \(C=0\). Integrate again, \( \frac{d \theta}{d \xi} = -\frac{\xi}{3} \).

Step 5: Find \(\theta\)

Integrate \( \frac{d \theta}{d \xi} = -\frac{\xi}{3} \) to obtain \( \theta(\xi) = -\frac{\xi^2}{6} + C \). Apply \( \theta(0) = 1 \), so \(C = 1\). Thus, \( \theta(\xi) = 1 - \frac{\xi^2}{6} \).

Step 6: Determine the First Zero of \(\theta(\xi)\)

To find \(\xi_1\), the point where \( \theta(\xi) = 0 \), solve the equation: \[ 1 - \frac{\xi^2}{6} = 0 \] Simplifying gives \(\xi^2 = 6\), hence \(\xi = \sqrt{6}\).

Step 7: Verify the Solution

Check that the boundary conditions and the obtained form \(\theta(\xi)=1 - \frac{\xi^2}{6}\) satisfy the simplified Lane-Emden equation and the specific condition \(\xi_1 = \sqrt{6}\) is correct.

Key Concepts

Gravitational Potential

Gravitational potential is a measure of the potential energy per unit mass at a specific point in a gravitational field. Imagine it as the work needed to bring a unit mass from infinity to that point, without altering the mass distribution creating the field. In the context of spherically symmetric bodies, like stars, the Lane-Emden equation describes how this potential varies with radius when the body has a consistent density and pressure gradient.

For stars and similar objects, gravitational potential is crucial in understanding the balance between the inward gravitational forces and the pressure pushing outwards. If the potential is too high, the star might contract or collapse under its own gravity, while too low a potential can cause expansion. Solving the Lane-Emden equation helps model these conditions by providing an equation that accounts for how matter is distributed within the star. The equation links potential directly to the pressure and density, which are also linked through state equations, such as polytropic models. These models are simplified representations but are essential for understanding stellar structure.

Polytropic Index

The polytropic index, denoted as \( n \), is a parameter that characterizes a specific relationship between a star's pressure and density. In polytropic models, the index \( n \) defines how pressure \( P \), density \( \rho \), and specific energy \( \epsilon \) relate to each other through the equation:\[P = K \rho^{1 + \frac{1}{n}}\]where \( K \) is a constant depending on the specific characteristics of the star or fluid.

• When \( n = 0 \), pressure does not depend on density, meaning the star is in an incompressible state. This gives rise to a uniform density model, which simplifies solving the Lane-Emden equation.
• Higher values of \( n \) indicate more complex relationships where pressure and density interact strongly, altering the compactness and internal structure of a star.

For example, the Sun can be approximately described by a polytropic index of 3, implying a more intricate balance between its pressure and density, unlike the simple, uniform case of \( n=0 \) used in the stepping solution to the Lane-Emden equation.

Ordinary Differential Equation (ODE)

An Ordinary Differential Equation (ODE) is an equation involving functions and their derivatives. In the context of the Lane-Emden equation, an ODE is used to describe how density varies with radius in a star or other spherically symmetric body. The challenge lies in solving these ODEs to predict or describe physical phenomena accurately.The Lane-Emden equation is a key example of an ODE derived from the principles of hydrostatic equilibrium and mass conservation, which describe how forces balance in a star. It's a second-order, nonlinear ODE that can be tricky to solve.

• The solution involves simplifying assumptions and enforcing boundary conditions, like central density and symmetry at the core, to find a usable form for \( \theta(\xi) \).
• Solving such equations often requires integration, step-by-step application of conditions for continuity and smoothness, and stability across the star's radius.

Overall, understanding ODEs is crucial for anyone attempting to model or study astrophysical objects, providing a foundation for interpreting stellar structures and behaviors.