Problem 22 - An Introduction to Modern Astrophysics

Question

(a) Assuming that Eq. ( \(\quad 51\) ) is valid for any arbitrary distance from the center of the Galaxy, show that the amount of dark matter interior to a radius \(r\) is given by the expression $$M_{r}=4 \pi \rho_{0} a^{2}\left[r-a \tan ^{-1}\left(\frac{r}{a}\right)\right].$$ $$\rho(r)=\frac{\rho_{0}}{1+(r / a)^{2}},$$ (b) If \(5.4 \times 10^{11} \mathrm{M}_{\odot}\) of dark matter is located within \(50 \mathrm{kpc}\) of the Galactic center, determine \(\rho_{0}\) in units of \(\mathrm{M}_{\odot} \mathrm{kpc}^{-1}\). Repeat your calculation if \(1.9 \times 10^{12} \mathrm{M}_{\odot}\) is located within \(230 \mathrm{kpc}\) of the Galactic center. Assume that \(a=2.8 \mathrm{kpc}\).

Step-by-Step Solution

Answer

\(\rho_{0}\) is calculated for each case based on the total mass and corresponding radius, using the derived mass formula.

Step 1: Understand the Given Density Function

We are provided with a density function: \( \rho(r) = \frac{\rho_{0}}{1+(r/a)^{2}} \). This describes how the density of dark matter changes with distance \(r\) from the center of the galaxy.

Step 2: Derive the Mass Formula

To find the total mass \(M_{r}\) within a radius \(r\), we integrate the density function over a spherical volume. The mass element \(dm\) is the product of the density and the differential volume element \(dV = 4\pi r^2 dr\). Thus, \[ dm = \rho(r) \cdot 4\pi r^2 dr. \] Integrate this from \(0\) to \(r\): \[ M_{r} = \int_{0}^{r} \rho(r') \cdot 4\pi r'^2 \, dr' = \int_{0}^{r} \frac{\rho_{0} \cdot 4\pi r'^2}{1 + (r'/a)^2} \, dr'. \]

Step 3: Solve the Integral

We solve the integral: \[ M_{r} = 4\pi \rho_{0} \int_{0}^{r} \frac{r'^2}{1 + (r'/a)^2} \, dr'. \] Use substitution \( u = \frac{r'}{a}\), \( du = \frac{1}{a} \, dr' \), so \( dr' = a \, du \). The limits change from \(0\) to \(\frac{r}{a}\). The integral becomes: \[ 4\pi \rho_{0} a^3 \int_{0}^{r/a} \frac{u^2}{1 + u^2} \, du. \] This integral can be solved using standard techniques: \( \int \frac{u^2}{1+u^2} \, du = u - \tan^{-1}(u) \). Evaluate from \(0\) to \(r/a\): \[ M_{r} = 4\pi \rho_{0} a^3 \left[ \frac{r}{a} - \tan^{-1}\left(\frac{r}{a}\right) \right]. \] Simplifying this gives: \[ M_{r} = 4 \pi \rho_{0} a^2 \left[ r - a \tan^{-1}\left( \frac{r}{a} \right) \right]. \]

Step 4: Verify the Result

This derived expression matches the given formula for \( M_{r} \), confirming that the derivation is correct: \[ M_{r}=4 \pi \rho_{0} a^{2}\left[r-a \tan ^{-1}\left(\frac{r}{a}\right)\right]. \]

Step 5: Calculate \(\rho_{0}\) for \(5.4 \times 10^{11} \mathrm{M}_{\odot}\) in 50kpc

Substitute \( M_{r} = 5.4 \times 10^{11} \mathrm{M}_{\odot} \), \( r = 50 \mathrm{kpc} \), and \( a = 2.8 \mathrm{kpc} \) into the mass formula: \[ 5.4 \times 10^{11} = 4 \pi \rho_{0} (2.8)^2 \left[ 50 - 2.8 \tan^{-1}\left( \frac{50}{2.8} \right) \right]. \] Solve for \( \rho_{0} \).

Step 6: Calculate \(\rho_{0}\) for \(1.9 \times 10^{12} \mathrm{M}_{\odot}\) in 230kpc

Similarly, substitute \( M_{r} = 1.9 \times 10^{12} \mathrm{M}_{\odot} \), \( r = 230 \mathrm{kpc} \), and \( a = 2.8 \mathrm{kpc} \) into the mass formula: \[ 1.9 \times 10^{12} = 4 \pi \rho_{0} (2.8)^2 \left[ 230 - 2.8 \tan^{-1}\left( \frac{230}{2.8} \right) \right]. \] Solve for \( \rho_{0} \).

Step 7: Numerically Solve for \(\rho_{0}\)

Plug the values into a calculator or computation software to solve the equations from steps 5 and 6 to find \(\rho_{0}\) for each case.

Key Concepts

Mass Distribution

In astrophysics, the concept of mass distribution plays a crucial role in understanding how objects like dark matter and stars are scattered across a galaxy. When we talk about mass distribution in the context of dark matter, it refers to how the total mass is spread out within a certain radius from the center of a galaxy. This helps astronomers determine not just how much dark matter is present, but also its influence on the movement of visible matter and the galaxy itself.

The density function is instrumental here, as it tells us how much mass is concentrated at different points in the galaxy. For dark matter, the density function is typically more concentrated around the center and decreases as you move further away. This is different from the distribution of visible matter, like stars, which tends to be more spread out. Understanding mass distribution helps in predicting gravitational effects and contributes to broader cosmological insights.

• Mass distribution shows us where mass is located within a galaxy.
• This helps predict gravitational effects and galaxy dynamics.
• Dark matter density decreases as you move away from the center.

Galaxy Density Function

The galaxy density function is a mathematical model that describes how the density of a particular component, like dark matter, changes with distance from the center of a galaxy. In our exercise, the function provided is:\[\rho(r) = \frac{\rho_{0}}{1 + (r/a)^2}\]This equation shows that the density at any distance \(r\) from the galaxy's center depends on two parameters: \(\rho_{0}\), which represents the core density, and \(a\), which is a scale factor related to the size of the central concentration.The density function is inversely related to the square of the distance. This means the further you go from the center, the density decreases significantly. Thus, it models how gravity pulls matter towards the center, compacting it more heavily at the core.

• Density decreases with distance, showing less mass in outer regions.
• It's important for describing how mass is distributed vis-à-vis distance.
• Helps in calculating gravitational pull exerted by dark matter.

Spherical Volume Integration

Spherical volume integration is a method used to calculate the total mass within a certain radius by integrating the density function over a spherical volume. When doing these calculations, especially when you deal with astronomical models, spherical coordinates are essential because celestial objects are frequently spherically symmetrical.To perform spherical integration, the mass element \(dm\) is expressed as the product of the density \(\rho(r)\) and the differential volume element \(dV\), which is given by \(4\pi r^2 dr\). Consequently, the integral to calculate mass \(M_{r}\) within a radius \(r\) becomes:\[M_{r} = \int_{0}^{r} \rho(r') \cdot 4\pi r'^2 \, dr'\]Substituting the provided density function into this equation allows us to solve for the total mass, by using integration techniques and occasionally variable substitutions, to manage the limits effectively. This results in understanding how much mass is contained and how it influences the galaxy's configuration.

• Essential for determining the total mass in a spherical region of space.
• Considered standard method for space-based calculations.
• Calculations consider the object’s symmetry, making it accurate for galaxies.

Astrophysics Calculations

Astrophysics calculations involve several mathematical methods to decipher the vast, complex universes' properties. Whether handling distances in kiloparsecs or masses in solar units, these calculations require converting physical quantities into mathematical expressions that explain cosmic phenomena.For example, calculating the core density \(\rho_{0}\) requires solving an equation like:\[M_{r} = 4 \pi \rho_{0} a^2 \left[r - a \tan^{-1}\left( \frac{r}{a} \right) \right]\]Used in conjunction with given conditions like specific mass values within certain radii, like dark matter mass within 50 kpc or 230 kpc, these calculations aid the determination of important parameters such as galaxy structure or mass distribution. Tools like computational software make solving these equations feasible, allowing astronomers to focus on interpreting what these numbers mean for understanding the universe.

• Converts physical insights into quantitative data.
• Predicts conditions and behaviors in galaxies and other celestial bodies.
• Requires transformation of complex equations into actionable insights.