Abekta

The Encyclopédie of CASSA

User Tools

Site Tools

You are here: root » Courses » AST 402: Stars and Interstellar Medium » Stellar Structure
courses:ast402:stellar-structure

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Next revision		Previous revision
courses:ast402:stellar-structure [2026/06/07 20:02] – created shuvocourses:ast402:stellar-structure [2026/06/07 22:08] (current) – [Stellar Models and Simulation] shuvo
Line 1: Line 1:
 ====== Stellar Structure ====== ====== Stellar Structure ======
  
-### **Hydrostatic Equilibrium**+===== Hydrostatic Equilibrium =====
  
 **Hydrostatic equilibrium** is the fundamental condition that describes the mechanical stability of a star, where the inward pull of gravity is precisely balanced by an outward pressure gradient. This balance ensures that the star remains static and does not undergo rapid expansion or collapse. Mathematically, for a spherically symmetric star, this is expressed by the **hydrostatic equilibrium equation**: **Hydrostatic equilibrium** is the fundamental condition that describes the mechanical stability of a star, where the inward pull of gravity is precisely balanced by an outward pressure gradient. This balance ensures that the star remains static and does not undergo rapid expansion or collapse. Mathematically, for a spherically symmetric star, this is expressed by the **hydrostatic equilibrium equation**:
Line 7: Line 7:
 where $P$ is pressure, $r$ is the radial distance from the center, $M_r$ is the interior mass within radius $r$, $\rho$ is the local density, and $g$ is the local acceleration of gravity. This equation dictates that **pressure must increase toward the center** of the star to support the weight of the overlying layers. where $P$ is pressure, $r$ is the radial distance from the center, $M_r$ is the interior mass within radius $r$, $\rho$ is the local density, and $g$ is the local acceleration of gravity. This equation dictates that **pressure must increase toward the center** of the star to support the weight of the overlying layers.
  
-### **Pressure Equation of State**+**Pressure Equation of State**
  
-A **pressure equation of state** relates the internal pressure of stellar material to its density, temperature, and chemical composition. The general **pressure integral** derived from the momentum transferred by particles (massive or massless) is:+A pressure equation of state relates the internal pressure of stellar material to its density, temperature, and chemical composition. The general **pressure integral** derived from the momentum transferred by particles (massive or massless) is:
 $$P = \frac{1}{3} \int_0^\infty p n_p v_p dp$$ $$P = \frac{1}{3} \int_0^\infty p n_p v_p dp$$
 where $p$ is momentum, $v$ is velocity, and $n_p$ is the number density of particles per momentum interval. In most stellar interiors, the total pressure is the sum of **gas pressure** ($P_g$) and **radiation pressure** ($P_{rad}$): where $p$ is momentum, $v$ is velocity, and $n_p$ is the number density of particles per momentum interval. In most stellar interiors, the total pressure is the sum of **gas pressure** ($P_g$) and **radiation pressure** ($P_{rad}$):
Line 15: Line 15:
 where $k$ is Boltzmann’s constant, $\mu$ is the mean molecular weight, $m_H$ is the mass of a hydrogen atom, and $a$ is the radiation constant. where $k$ is Boltzmann’s constant, $\mu$ is the mean molecular weight, $m_H$ is the mass of a hydrogen atom, and $a$ is the radiation constant.
  
-### **Ideal Gas Law**+**Ideal Gas Law**
  
-The **ideal gas law** applies to non-degenerate, non-relativistic stellar gases where particles interact only through perfectly elastic collisions. Derived using the **Maxwell-Boltzmann velocity distribution**, it is expressed as:+The ideal gas law applies to non-degenerate, non-relativistic stellar gases where particles interact only through perfectly elastic collisions. Derived using the **Maxwell-Boltzmann velocity distribution**, it is expressed as:
 $$P_g = n k T = \frac{\rho kT}{\mu m_H}$$ $$P_g = n k T = \frac{\rho kT}{\mu m_H}$$
 where $n$ is the particle number density. The **mean molecular weight** ($\mu$) depends on the ionization state and composition of the gas. For a **completely ionized gas**, it is approximated by: where $n$ is the particle number density. The **mean molecular weight** ($\mu$) depends on the ionization state and composition of the gas. For a **completely ionized gas**, it is approximated by:
Line 23: Line 23:
 where $X, Y, \text{ and } Z$ are the mass fractions of hydrogen, helium, and metals, respectively. where $X, Y, \text{ and } Z$ are the mass fractions of hydrogen, helium, and metals, respectively.
  
-### **Fermi-Dirac and Bose-Einstein Statistics**+**Fermi-Dirac and Bose-Einstein Statistics**
  
-When classical assumptions fail, **quantum statistics** are required to describe particle distributions: +When classical assumptions fail, **quantum statistics** are required to describe particle distributions:\\ 
-*   **Fermi-Dirac Statistics:** Apply to **fermions** (e.g., electrons, protons, neutrons), which obey the **Pauli exclusion principle**. In extremely dense environments like white dwarfs, **electron degeneracy pressure** becomes dominant. + 
-    *   **Non-relativistic degeneracy:** $P \approx \frac{(3\pi^2)^{2/3}}{5} \frac{\hbar^2}{m_e} \left( \frac{Z}{A} \frac{\rho}{m_H} \right)^{5/3}$. +Fermi-Dirac Statistics: Apply to **fermions** (e.g., electrons, protons, neutrons), which obey the **Pauli exclusion principle**. In extremely dense environments like white dwarfs, **electron degeneracy pressure** becomes dominant.\\ 
-    *   **Relativistic degeneracy:** $P \approx \frac{(3\pi^2)^{1/3}}{4} \hbar c \left( \frac{Z}{A} \frac{\rho}{m_H} \right)^{4/3}$. + 
-*   **Bose-Einstein Statistics:** Apply to **bosons** (e.g., photons), which do not obey the exclusion principle. This leads to the formulation of **radiation pressure**:+Non-relativistic degeneracy: $P \approx \frac{(3\pi^2)^{2/3}}{5} \frac{\hbar^2}{m_e} \left( \frac{Z}{A} \frac{\rho}{m_H} \right)^{5/3}$. 
 +Relativistic degeneracy: $P \approx \frac{(3\pi^2)^{1/3}}{4} \hbar c \left( \frac{Z}{A} \frac{\rho}{m_H} \right)^{4/3}$.\\ 
 + 
 +Bose-Einstein Statistics: Apply to **bosons** (e.g., photons), which do not obey the exclusion principle. This leads to the formulation of **radiation pressure**:
 $$P_{rad} = \frac{1}{3} a T^4$$ $$P_{rad} = \frac{1}{3} a T^4$$
 which can dominate in very massive, hot stars. which can dominate in very massive, hot stars.
  
-### **Stellar Energy Sources**+===== Stellar Energy Sources =====
  
-Stars are powered by two primary energy sources: **gravitation** and **nuclear fusion**. +Stars are powered by two primary energy sources: **gravitation** and **nuclear fusion**.\\ 
-*   **Gravitational Energy:** Released during contraction, with the total potential energy approximated as $U_g \approx -\frac{3}{5} \frac{GM^2}{R}$. According to the **virial theorem**, half is radiated and half heats the interior. + 
-*   **Nuclear Energy:** Principally the fusion of light nuclei into heavier ones, releasing energy based on $E=mc^2$. The **nuclear energy generation rate** ($\epsilon$) is often modeled as a power law:+**Gravitational Energy:** Released during contraction, with the total potential energy approximated as $U_g \approx -\frac{3}{5} \frac{GM^2}{R}$. According to the **virial theorem**, half is radiated and half heats the interior.\\ 
 +**Nuclear Energy:** Principally the fusion of light nuclei into heavier ones, releasing energy based on $E=mc^2$. The **nuclear energy generation rate** ($\epsilon$) is often modeled as a power law:
 $$\epsilon = \epsilon_0 X_i X_x \rho^\alpha T^\beta$$ $$\epsilon = \epsilon_0 X_i X_x \rho^\alpha T^\beta$$
 where $\beta$ varies significantly by reaction type (e.g., $\beta \approx 4$ for the pp chain and $\beta \approx 20$ for the CNO cycle). where $\beta$ varies significantly by reaction type (e.g., $\beta \approx 4$ for the pp chain and $\beta \approx 20$ for the CNO cycle).
  
-### **Timescales**+**Timescales:**
  
-Three primary timescales characterize stellar life: +Three primary timescales characterize stellar life:\\ 
-1.  **Dynamic (Free-fall) Timescale:** The time for a cloud to collapse under gravity if pressure is removed: $t_{ff} = \sqrt{\frac{3\pi}{32} \frac{1}{G \rho_0}}$. +1.  **Dynamic (Free-fall) Timescale:** The time for a cloud to collapse under gravity if pressure is removed: $t_{ff} = \sqrt{\frac{3\pi}{32} \frac{1}{G \rho_0}}$.\\ 
-2.  **Thermal (Kelvin-Helmholtz) Timescale:** The time to radiate away a star's total gravitational energy: $t_{KH} \approx \frac{3}{10} \frac{GM^2}{RL}$.+2.  **Thermal (Kelvin-Helmholtz) Timescale:** The time to radiate away a star's total gravitational energy: $t_{KH} \approx \frac{3}{10} \frac{GM^2}{RL}$.\\
 3.  **Nuclear Timescale:** The time to exhaust nuclear fuel: $t_{nuc} \approx \frac{E_{nuc}}{L}$. For the Sun, $t_{nuc} \approx 10^{10}$ years. 3.  **Nuclear Timescale:** The time to exhaust nuclear fuel: $t_{nuc} \approx \frac{E_{nuc}}{L}$. For the Sun, $t_{nuc} \approx 10^{10}$ years.
  
-### **Quantum Tunneling**+**Quantum Tunneling:**
  
 Classical physics suggests stellar cores are too cool for nuclei to overcome the **Coulomb barrier** through thermal motion alone. **Quantum mechanical tunneling** allows particles to penetrate this barrier even when their kinetic energy is insufficient. The probability of tunneling is exponential: Classical physics suggests stellar cores are too cool for nuclei to overcome the **Coulomb barrier** through thermal motion alone. **Quantum mechanical tunneling** allows particles to penetrate this barrier even when their kinetic energy is insufficient. The probability of tunneling is exponential:
Line 54: Line 58:
 where $b$ is a constant related to particle charges and masses. This probability combined with the Maxwell-Boltzmann distribution defines the **Gamow peak**, the narrow energy range where most nuclear reactions occur. where $b$ is a constant related to particle charges and masses. This probability combined with the Maxwell-Boltzmann distribution defines the **Gamow peak**, the narrow energy range where most nuclear reactions occur.
  
-### **Nucleosynthesis**+**Nucleosynthesis:**
  
-**Nucleosynthesis** is the sequence of nuclear reactions that transform elements: +Nucleosynthesis is the sequence of nuclear reactions that transform elements: 
-*   **Hydrogen Burning:** Occurs via the **proton-proton (pp) chains** (dominant in low-mass stars) or the **CNO cycle** (dominant in massive stars). +Hydrogen Burning: Occurs via the **proton-proton (pp) chains** (dominant in low-mass stars) or the **CNO cycle** (dominant in massive stars).\\ 
-*   **Helium Burning:** Occurs via the **triple-alpha process** ($3\alpha \to ^{12}C$) at temperatures $\sim 10^8$ K. +Helium Burning: Occurs via the **triple-alpha process** ($3\alpha \to ^{12}C$) at temperatures $\sim 10^8$ K.\\ 
-*   **Advanced Burning:** Successive stages (carbon, oxygen, neon, and silicon burning) produce elements up to the **iron peak**. Elements heavier than iron are produced via the **s-process** (slow neutron capture) or **r-process** (rapid neutron capture).+Advanced Burning: Successive stages (carbon, oxygen, neon, and silicon burning) produce elements up to the **iron peak**. Elements heavier than iron are produced via the **s-process** (slow neutron capture) or **r-process** (rapid neutron capture).
  
-### **Energy Transport and Thermodynamics**+===== Energy Transport and Thermodynamics =====
  
 Energy is transported from the core to the surface via **radiation, convection, or conduction**. Energy is transported from the core to the surface via **radiation, convection, or conduction**.
-*   **Radiative Transport:** Driven by the radiation pressure gradient:+**Radiative Transport:** Driven by the radiation pressure gradient:
 $$\frac{dT}{dr} = -\frac{3}{4ac} \frac{\kappa \rho}{T^3} \frac{L_r}{4\pi r^2}$$ $$\frac{dT}{dr} = -\frac{3}{4ac} \frac{\kappa \rho}{T^3} \frac{L_r}{4\pi r^2}$$
 where $\kappa$ is the opacity. where $\kappa$ is the opacity.
-*   **Convection:** Occurs when the temperature gradient is **superadiabatic**. The **Schwarzschild criterion** for convection (neglecting composition changes) is:+ 
 +**Convection:** Occurs when the temperature gradient is **superadiabatic**. The **Schwarzschild criterion** for convection (neglecting composition changes) is:
 $$\left| \frac{dT}{dr} \right|_{act} > \left| \frac{dT}{dr} \right|_{ad}$$ $$\left| \frac{dT}{dr} \right|_{act} > \left| \frac{dT}{dr} \right|_{ad}$$
 where the adiabatic gradient for an ideal gas is $\frac{dT}{dr}|_{ad} = -\left(1 - \frac{1}{\gamma}\right) \frac{\mu m_H}{k} \frac{GM_r}{r^2}$. where the adiabatic gradient for an ideal gas is $\frac{dT}{dr}|_{ad} = -\left(1 - \frac{1}{\gamma}\right) \frac{\mu m_H}{k} \frac{GM_r}{r^2}$.
-*   **Thermodynamics:** The **first law** ($dU = dQ - dW$) relates internal energy changes to heat and work. For an adiabatic process, $P V^\gamma = \text{constant}$. 
  
-### **Stellar Models and Simulation**+**Thermodynamics:** The **first law** ($dU = dQ - dW$) relates internal energy changes to heat and work. For an adiabatic process, $P V^\gamma = \text{constant}$. 
 + 
 +===== Stellar Models and Simulation ===== 
 + 
 +Stellar models are constructed by numerically solving the four fundamental differential equations (hydrostatic equilibrium, mass conservation, energy generation, and energy transport). \\ 
 + 
 +**Numerical Modeling:** The star is divided into **concentric shells (zones)**, and differential equations are replaced by **difference equations**. Solutions require matching **boundary conditions** at the center ($M_r \to 0, L_r \to 0$) and surface ($P \to 0, T \to 0$).\\
  
-**Stellar models** are constructed by numerically solving the four fundamental differential equations (hydrostatic equilibrium, mass conservation, energy generation, and energy transport).  +**Vogt-Russell Theorem:** States that a star's mass and composition uniquely determine its structure and evolution.\\ 
-*   **Numerical Modeling:** The star is divided into **concentric shells (zones)**, and differential equations are replaced by **difference equations**. Solutions require matching **boundary conditions** at the center ($M_r \to 0, L_r \to 0$) and surface ($P \to 0, T \to 0$). +**Polytropes:** Simplified analytical models where pressure is a power of density ($P = K\rho^{(n+1)/n}$); they are solved using the **Lane-Emden equation**:
-*   **Vogt-Russell Theorem:** States that a star's mass and composition uniquely determine its structure and evolution. +
-*   **Polytropes:** Simplified analytical models where pressure is a power of density ($P = K\rho^{(n+1)/n}$); they are solved using the **Lane-Emden equation**:+
 $$\frac{1}{\xi^2} \frac{d}{d\xi} \left( \xi^2 \frac{dD_n}{d\xi} \right) = -D_n^n$$ $$\frac{1}{\xi^2} \frac{d}{d\xi} \left( \xi^2 \frac{dD_n}{d\xi} \right) = -D_n^n$$
 where $\xi$ is a dimensionless radius and $D_n$ is a dimensionless density function. where $\xi$ is a dimensionless radius and $D_n$ is a dimensionless density function.
courses/ast402/stellar-structure.1780884161.txt.gz · Last modified: by shuvo
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki