Solar Structure & Models

1.1 Hydrostatic Equilibrium

The Sun is a self-gravitating ball of plasma in near-perfect mechanical equilibrium. At every radial shell, the inward pull of gravity is exactly balanced by the outward pressure gradient. This is the most fundamental equation of stellar structure, and its validity is confirmed by the fact that the Sun's radius changes by less than one part in 10⁵ over a solar cycle.

Derivation 1: Hydrostatic Equilibrium from Force Balance

Consider a thin spherical shell of thickness \(dr\) at radius \(r\)with density \(\rho(r)\). The shell has mass per unit area \(dm = \rho \, dr\).

Step 1. The gravitational force per unit area pulling the shell inward is:

where \(M(r)\) is the mass enclosed within radius \(r\).

Step 2. The net pressure force per unit area on the shell is the difference between the pressure on the inner and outer faces:

Step 3. In static equilibrium, \(f_{\text{grav}} + f_{\text{press}} = 0\), so:

Central pressure estimate. Replacing \(M(r)\) with \(M_\odot\),\(r\) with \(R_\odot\), and \(\rho\) with the mean density\(\bar{\rho} = 3M_\odot / (4\pi R_\odot^3) \approx 1410\) kg/m\(^3\):

The actual central pressure from the Standard Solar Model is \(P_c \approx 2.34 \times 10^{16}\) Pa, which is about 500 times larger because the density is strongly concentrated toward the center (\(\rho_c \approx 1.5 \times 10^5\) kg/m\(^3\)).

1.2 Mass Continuity

Derivation 2: Mass Conservation Equation

The mass enclosed within radius \(r\) increases as we move outward by the mass in each spherical shell. This is the simplest of the four structure equations, but it couples the mass profile to the density profile.

Step 1. The mass of a thin spherical shell of radius \(r\) and thickness \(dr\) is:

Step 2. Dividing both sides by \(dr\):

Step 3. The boundary conditions are \(M(0) = 0\) at the center and\(M(R_\odot) = M_\odot = 1.989 \times 10^{30}\) kg at the surface.

The mass continuity equation tells us that about 90% of the solar mass is contained within the inner 50% of the radius, reflecting the strong central concentration.

1.3 Energy Transport

Derivation 3: Radiative Temperature Gradient

Energy is transported from the core to the surface by either radiation or convection. In the radiative zone (from about 0.25 to 0.71 \(R_\odot\)), photons carry the energy outward through a random walk process.

Step 1. The radiative flux at radius \(r\) is given by Fick's law of diffusion applied to photons:

where \(\kappa\) is the Rosseland mean opacity, \(a\) is the radiation constant, and \(c\) is the speed of light.

Step 2. The luminosity at radius \(r\) is \(L(r) = 4\pi r^2 F_{\text{rad}}\). Solving for the temperature gradient:

Step 3. The Schwarzschild criterion determines whether convection operates: if the actual temperature gradient exceeds the adiabatic gradient,

then the layer is convectively unstable. In the Sun, the outer 29% by radius (the convection zone) satisfies this criterion due to the increased opacity from partially ionized hydrogen and helium.

In the convection zone, the temperature gradient is very nearly adiabatic:\(\nabla = \nabla_{\text{ad}} = 0.4\) for a monatomic ideal gas (\(\gamma = 5/3\)). The mixing-length theory (MLT) of Bohm-Vitense (1958) parameterizes the convective flux in terms of a single parameter \(\alpha_{\text{MLT}} \approx 1.5\text{--}2.0\).

1.4 Energy Generation

Derivation 4: Luminosity Equation

The luminosity increases outward as nuclear reactions generate energy in the core. The energy generation rate \(\varepsilon\) depends sensitively on temperature, making the luminosity strongly concentrated in the inner \(\sim 0.25 R_\odot\).

Step 1. In a thin shell of mass \(dM = 4\pi r^2 \rho \, dr\), the energy generated per unit time is:

Step 2. Dividing by \(dr\):

Step 3. For the pp chain, the energy generation rate scales as:

The \(T^4\) dependence means that a 10% increase in temperature roughly doubles the energy generation rate, confining most nuclear burning to the hottest central region.

The boundary conditions are \(L(0) = 0\) and \(L(R_\odot) = L_\odot = 3.828 \times 10^{26}\) W. The nuclear timescale \(\tau_{\text{nuc}} = E_{\text{nuc}} / L_\odot \approx 10^{10}\) years sets the main-sequence lifetime of the Sun.

1.5 The Standard Solar Model

Derivation 5: Equation of State and Closure

The four structure equations contain five unknowns: \(P, T, \rho, M, L\) as functions of\(r\). We need an equation of state (EOS) to close the system.

Step 1. For an ideal gas plus radiation pressure:

where \(\mu\) is the mean molecular weight. For a fully ionized gas of mass fractions\(X\) (hydrogen), \(Y\) (helium), \(Z\) (metals):

Step 2. For the solar center with \(X_c \approx 0.34\), \(Y_c \approx 0.64\):

Step 3. The complete Standard Solar Model (Bahcall et al.) is obtained by integrating the coupled ODEs from center to surface, with boundary conditions:

The modern SSM predictions (Vinyoles et al. 2017) are validated to better than 0.1% in the sound speed profile by helioseismology, making the Sun the best-tested stellar model in astrophysics.

Applications & Physical Scales

Numerical Simulation

This simulation integrates the four stellar structure equations for a simplified solar model, producing radial profiles of pressure, temperature, density, mass, and luminosity.

Summary

Full Derivation: Hydrostatic Equilibrium from First Principles

We derive the equation of hydrostatic equilibrium by considering the Newtonian force balance on a thin spherical shell element within the star. This is the most fundamental constraint governing the mechanical structure of all self-gravitating bodies.

Step 1: Define the Shell Element

Consider a thin spherical shell at radius \(r\) with infinitesimal thickness \(dr\). The shell subtends a solid angle \(d\Omega\) and has a cross-sectional area element\(dA = r^2 d\Omega\). The mass of this shell element is:

Step 2: Gravitational Force

By the shell theorem (Newton, Principia Book I, Prop. 71), only the mass interior to the shell exerts a net gravitational pull. The gravitational acceleration at radius \(r\) is:

The inward gravitational force on the shell element per unit area is:

Step 3: Pressure Gradient Force

The pressure on the inner face of the shell is \(P(r)\) (pushing outward) and on the outer face is \(P(r + dr) = P(r) + (dP/dr) \, dr\) (pushing inward). The net outward pressure force per unit area is:

Step 4: Newton's Second Law

For a static star, the acceleration is zero. Setting the net force to zero:

Dividing through by \(dr\) gives the equation of hydrostatic equilibrium:

Step 5: Virial Theorem Connection

Multiplying both sides by \(4\pi r^3\) and integrating from center to surface gives the virial theorem. Start with:

The left side integrates by parts (noting \(P(R) = 0\)):

Identifying the gravitational potential energy \(\Omega = -\int_0^M \frac{GM(r)}{r} dM\)and using \(P = (\gamma - 1)\rho u\) for the internal energy density \(u\):

Full Derivation: Energy Transport and Convection Onset

Radiative Temperature Gradient from Photon Diffusion

Step 1. In the stellar interior, photons undergo a random walk with mean free path \(\ell = 1/(\kappa \rho)\). The radiative energy flux follows from Fick's law applied to the photon energy density \(u_{\text{rad}} = a T^4\):

Step 2. Substitute \(\ell = 1/(\kappa\rho)\) and evaluate\(d(aT^4)/dr = 4aT^3 \, dT/dr\):

Step 3. The luminosity at radius \(r\) is carried entirely by radiation:\(L(r) = 4\pi r^2 F_{\text{rad}}\). Solving for \(dT/dr\):

Schwarzschild Criterion for Convection Onset

Step 4. A fluid element displaced upward by \(\delta r\) expands adiabatically. Its temperature after displacement is:

The surrounding temperature is:

Step 5. If the element is hotter than its surroundings, it is less dense and continues to rise (convectively unstable). This requires:

Step 6. In terms of the dimensionless temperature gradients\(\nabla = d\ln T / d\ln P\):

In the Sun, the radiative gradient exceeds \(\nabla_{\text{ad}}\) for\(r > 0.713 R_\odot\) due to the sharp increase in opacity from partial ionization of H and He. This defines the base of the convection zone, a critical boundary confirmed by helioseismology to be at \(r_{\text{cz}} = (0.713 \pm 0.001) R_\odot\).

Solar Interior Cross-Section

The Sun's interior is divided into three main zones with distinct energy transport mechanisms and physical conditions:

Standard Solar Model Profiles

This simulation plots T(r), \(\rho\)(r), P(r), and L(r) profiles from a calibrated Standard Solar Model, including the convection zone boundary and core-envelope transition.