The Sun, Stellar Interiors & Nuclear Energy Generation

The Sun is a middle-aged, G2V spectral dwarf star containing $99.86\%$ of the total mass of the Solar System. Understanding its physical parameters provides the canonical baseline ('solar units') against which all other stars are measured.

Global Solar Parameters

Solar Composition by Mass

Spectroscopic abundance analyses and helioseismic inversions establish the initial zero-age solar composition as:

Due to $4.57\text{ Gyr}$ of nuclear burning, the modern core has burned hydrogen into helium, shifting its central abundances to $X_c \approx 0.34$, $Y_c \approx 0.64$.

Johannes Kepler formulated three empirical laws of planetary motion based on Tycho Brahe's precision observational data. Isaac Newton later derived these laws from first principles using universal gravitation and Newtonian mechanics in a reduced two-body framework.

1. Kepler's First Law (The Law of Ellipses)

Every planet moves in an elliptical orbit with the Sun situated at one of the two foci.

In polar coordinates $(r, \theta)$ centered on the Sun, the equation of an ellipse is:

where $a$ is the semi-major axis, $e$ is the orbital eccentricity ($0 \le e < 1$), $\theta$ is the true anomaly (angle from perihelion), and $p = a(1-e^2)$ is the semi-latus rectum. The perihelion distance is $r_p = a(1-e)$ and the aphelion distance is $r_a = a(1+e)$.

2. Kepler's Second Law (The Law of Equal Areas)

A line joining a planet and the Sun sweeps out equal areas during equal intervals of time.

This is a direct mathematical consequence of the conservation of orbital angular momentum $\vec{L}$ in a central force field:

where $\mu = \frac{m_1 m_2}{m_1 + m_2}$ is the reduced mass. Thus, a planet travels fastest at perihelion ($v_p = \sqrt{\frac{GM}{a}\frac{1+e}{1-e}}$) and slowest at aphelion ($v_a = \sqrt{\frac{GM}{a}\frac{1-e}{1+e}}$).

3. Kepler's Third Law (The Harmonic Law)

The square of the orbital period $P$ of a planet is proportional to the cube of its semi-major axis $a$.

Integrating the areal velocity $\frac{dA}{dt} = \frac{L}{2\mu}$ over one full orbital period $P$ yields the area of an ellipse $A = \pi a b = \pi a^2 \sqrt{1-e^2}$:

Substituting $L/\mu = \sqrt{G(M_1 + M_2)a(1-e^2)}$ yields Newton's generalized form of Kepler's 3rd Law:

When mass is in solar masses ($M_\odot$), period in Earth years ($\text{yr}$), and distance in Astronomical Units ($\text{AU}$), this simplifies to $P^2 = a^3 / (M_1 + M_2)$.

The Vis-Viva Equation

Conservation of total specific orbital mechanical energy $\mathcal{E} = \frac{1}{2}v^2 - \frac{GM}{r} = -\frac{GM}{2a}$ yields the fundamental vis-viva equation relating instantaneous orbital speed $v$ to radius $r$:

Since the 1995 discovery of 51 Pegasi b, the detection of exoplanets has transformed astronomy. Two dominant quantitative techniques account for the vast majority of discoveries: the spectroscopic radial velocity method and photometric transit observations.

1. The Radial Velocity (Doppler Wobble) Method

A planet of mass $m_p$ orbiting a host star of mass $M_*$ ($m_p \ll M_*$) in an orbit inclined at angle $i$ relative to the sky plane causes the star to reflexively orbit the common center of mass. The observable line-of-sight radial velocity of the star varies sinusoidally:

where $\gamma$ is the systemic velocity and $K$ is the semi-amplitude of the radial velocity:

For a circular orbit ($e=0$) with $m_p \ll M_*$:

Radial velocity observations directly determine the minimum planetary mass $m_p \sin i$. For a Jupiter-mass planet at $1\text{ AU}$ around a solar-type star, $K \approx 28.4\text{ m/s}$; for an Earth-mass planet at $1\text{ AU}$, $K \approx 8.9\text{ cm/s}$, requiring extreme spectrometer precision (e.g., ESPRESSO, HARPS).

2. Photometric Transit Method

If the planetary orbital plane is oriented nearly edge-on ($i \approx 90^\circ$), the planet transits the stellar disk once per orbit. The geometric transit probability is:

During transit, the fractional decrease in stellar flux (the transit depth) is directly proportional to the ratio of geometric surface areas:

For a Jupiter-sized planet ($R_p \approx 0.1 R_\odot$), $\delta \approx 1\%$. For an Earth-sized planet ($R_p \approx 0.009 R_\odot$), $\delta \approx 8.4 \times 10^{-5} \approx 84\text{ ppm}$ (parts per million), which space telescopes like Kepler and TESS routinely measure.

Combining radial velocity ($m_p \sin i$ with $i$ determined from transit duration) and transit measurements ($R_p$) uniquely yields the planet's true mass, physical radius, and mean bulk density $\bar{\rho}_p = \frac{m_p}{\frac{4}{3}\pi R_p^3}$, distinguishing rocky terrestrial worlds from water worlds, gas giants, and puffy planets.

A star is a stable, self-gravitating plasma configuration. Its macroscopic structure is governed by an exact mechanical balance between inward gravitational force and outward thermal and radiation pressure gradients.

Equation of Hydrostatic Equilibrium

Consider an infinitesimal cylindrical fluid element of cross-sectional area $dA$, height $dr$, and density $\rho(r)$ situated at radius $r$ inside a spherically symmetric star. The mass of the element is $dm = \rho(r) dA dr$. The forces acting radially are:

1. Inward gravitational force: $dF_g = -\frac{G M(r) dm}{r^2} = -\frac{G M(r) \rho(r) dA dr}{r^2}$
2. Net outward pressure force: $dF_p = [P(r) - P(r + dr)] dA = -\frac{dP}{dr} dr dA$

Setting the sum of forces to zero ($dF_g + dF_p = 0$) and dividing by $dA dr$ gives the fundamental equation of hydrostatic equilibrium:

Paired with the mass conservation equation:

Derivation of the Virial Theorem for Stars

Multiply both sides of the hydrostatic equation by $4\pi r^3 dr$ and integrate from the center ($r=0$) to the stellar surface ($r=R$):

The right-hand side is precisely the total gravitational potential energy of the star: $\Omega = -\int_0^{M_*} \frac{G M(r)}{r} dM(r) < 0$.

Integrating the left-hand side by parts:

Equating both sides: $3 \int_V P\,dV = -\Omega$.

For a non-relativistic classical monoatomic ideal gas, the pressure is related to the thermal kinetic energy density $u_k$ by $P = \frac{2}{3} u_k$. Therefore:

Thermodynamic Implications of the Virial Theorem

The total mechanical energy of the star is:

This profound result has critical consequences for stellar physics:

• Negative Heat Capacity: As a star loses energy by radiating luminosity $L = -\frac{dE}{dt} > 0$, its total energy $E$ becomes more negative, which means $\Omega$ decreases (the star contracts), and $K$ increases ($dK = -dE > 0$). Contracting stars heat up!
• Kelvin-Helmholtz Timescale: The timescale over which a star can shine solely by releasing gravitational potential energy is:
  $$\tau_{\text{KH}} = \frac{|\Omega|}{2 L} \approx \frac{3 G M^2}{10 R L} \sim 3 \times 10^7\text{ years for the Sun}$$
  Because Earth's geological record proves ages $> 4\text{ Gyr}$, gravitational contraction cannot power the Sun, necessitating nuclear fusion.

Energy generated in the deep stellar core must be transported outward to the surface. Stars utilize three potential mechanisms: conduction (negligible except in degenerate white dwarfs), radiative diffusion, and fluid convection.

1. Radiative Transport & The Diffusion Approximation

Inside a star, the photon mean free path $\ell_{\text{mfp}} = \frac{1}{\kappa \rho} \sim 1\text{ mm} - 1\text{ cm}$ is microscopic compared to the stellar radius. Photons undergo trillions of random walk scatterings in near-perfect local thermodynamic equilibrium (LTE). The radiative flux $F_{\text{rad}}$ is governed by Fick's law of diffusion for radiation energy density $u = a_{\text{rad}} T^4$:

Since the total luminosity flowing through a sphere of radius $r$ is $L(r) = 4\pi r^2 F_{\text{rad}}$, we obtain the equation of radiative temperature gradient:

where $\kappa$ is the Rosseland mean opacity ($\text{cm}^2\text{ g}^{-1}$), arising from electron Thomson scattering ($\kappa_{\text{es}} = 0.40\text{ cm}^2\text{ g}^{-1}$ for pure hydrogen) and Kramers' free-free / bound-free opacity ($\kappa \propto \rho T^{-7/2}$).

2. The Schwarzschild Criterion for Convective Instability

Consider a fluid parcel displaced adiabatically upward by distance $\Delta r > 0$. The parcel expands to match the surrounding ambient pressure ($P_p = P_a$), but retains its own entropy. If the resulting parcel density is lower than the surrounding ambient density ($\rho_p < \rho_a$), the parcel experiences a positive Archimedean buoyant force and continues accelerating upward—the layer is convectively unstable.

Comparing temperature gradients gives the Schwarzschild criterion:

Expressed in terms of the logarithmic temperature gradient $\nabla \equiv \frac{d\ln T}{d\ln P}$:

Convection is triggered whenever:

• The opacity $\kappa$ is extremely high (e.g., in cool outer stellar envelopes where hydrogen/helium recombine, driving $\nabla_{\text{rad}} \propto \kappa$ up).
• The nuclear energy generation is intensely concentrated at the center (e.g., in massive stars where the CNO cycle produces huge core $L(r)/r^2$).

Consequently, the Sun possesses a radiative core ($0 \le r \le 0.71 R_\odot$) surrounded by an outer convective envelope ($0.71 R_\odot \le r \le R_\odot$). Massive stars ($M > 1.3 M_\odot$) feature convective cores and radiative envelopes.

At classical temperatures of $T_c \sim 1.5 \times 10^7\text{ K}$, the average thermal kinetic energy of protons is $k T \approx 1.3\text{ keV}$. However, the electrostatic Coulomb repulsive potential barrier between two approaching protons separated by nuclear radius $r_n \approx 1\text{ fm}$ is:

Classically, zero protons possess enough energy to fuse. Stellar fusion is possible exclusively through quantum mechanical wave tunneling.

The Gamow Peak

The fusion cross-section $\sigma(E)$ decomposes into the geometrical de Broglie area $\pi \lambda^2 \propto 1/E$, the quantum mechanical barrier penetration probability $P(E)$, and an intrinsic nuclear reaction structure factor $S(E)$:

where $E_G = 2 m_r c^2 (\pi \alpha Z_1 Z_2)^2$ is the Gamow energy ($m_r$ is the reduced mass). The total nuclear reaction rate per unit volume is the Maxwell-Boltzmann average:

The integrand is the product of a falling Maxwellian tail $\exp(-E/kT)$ and a rising quantum tunneling factor $\exp(-\sqrt{E_G/E})$. The product forms a sharp, localized Gaussian known as the Gamow Peak centered at energy:

1. The Proton-Proton ($p$-$p$) Chain

Dominates in stars with $M \le 1.3 M_\odot$ ($T_c < 1.8 \times 10^7\text{ K}$). Net reaction: $4\,{}^1\text{H} \to {}^4\text{He} + 2e^+ + 2\nu_e + 26.73\text{ MeV}$.

1. PP-I Initiation: $${}^1\text{H} + {}^1\text{H} \to {}^2\text{H} + e^+ + \nu_e \quad (Q = 1.442\text{ MeV}, \tau \sim 10^{10}\text{ yr - weak force bottleneck})$$ $${}^2\text{H} + {}^1\text{H} \to {}^3\text{He} + \gamma \quad (Q = 5.494\text{ MeV}, \tau \sim 1\text{ s})$$ $${}^3\text{He} + {}^3\text{He} \to {}^4\text{He} + 2\,{}^1\text{H} \quad (Q = 12.86\text{ MeV}, \tau \sim 10^6\text{ yr})$$
2. PP-II Branch ($14\%$ in Sun): Produces ${}^7\text{Be}$ and ${}^7\text{Li}$, terminating in ${}^4\text{He}$.
3. PP-III Branch ($0.02\%$ in Sun): Produces high-energy ${}^8\text{B}$ neutrinos ($E_\nu \le 14\text{ MeV}$) detected by Super-Kamiokande and SNO.

The energy generation rate scales as $\epsilon_{pp} \propto \rho X^2 T^4$.

2. The CNO (Carbon-Nitrogen-Oxygen) Cycle

Uses pre-existing carbon, nitrogen, and oxygen as nuclear catalysts:

Because the Coulomb barrier between protons and $Z=6, 7, 8$ nuclei is much higher ($E_G$ is large), the CNO cycle has a ferocious temperature sensitivity: $\epsilon_{\text{CNO}} \propto \rho X Z_{\text{CNO}} T^{17}$. Above $T \approx 1.8 \times 10^7\text{ K}$ ($M \gtrsim 1.3 M_\odot$), the CNO cycle completely surpasses the $p$-$p$ chain.

Above the dense convective envelope lies the stratified solar atmosphere: the photosphere, chromosphere, transition region, and corona.

1. The Photosphere & Limb Darkening

The photosphere is the optical surface layer (optical depth $\tau_\nu \approx 2/3$, thickness $\sim 400\text{ km}$) from which photons escape into space. Because we look along slant paths near the apparent edge (limb) of the solar disk, our line of sight penetrates to a shallower, cooler physical depth compared to disk center. The resulting limb darkening is described by the Eddington-Barbier approximation:

where $\theta$ is the angle between the emergent ray and the surface normal.

2. The Solar Corona & Coronal Heating Problem

Above the chromosphere ($T \sim 10^4\text{ K}$) and a razor-thin transition region lies the corona, where temperatures paradoxically surge to $T \sim 1-3 \times 10^6\text{ K}$, emitting energetic soft X-rays. This thermal inversion violates simple conductive/radiative thermodynamic transfer from the $5778\text{ K}$ photosphere. Modern solar magnetohydrodynamics (MHD) attributes coronal heating to:

• Magnetic Reconnection & Nanoflares: Tangling and sudden topological reconnection of complex magnetic flux tubes driving localized bursts of magnetic energy dissipation (Parker nanoflare model).
• MHD Alfvén Wave Dissipation: Upward-propagating Alfvén waves excited by convective granulation turbulent churning that dissipate energy nonlinearly in the low-density coronal plasma.

3. The Solar Wind: Parker's Hydrodynamic Solution

Eugene Parker (1958) demonstrated that a static, isothermal corona cannot maintain zero pressure at spatial infinity ($P(\infty) > 0$), requiring a continuous, supersonic hydrodynamic expansion known as the solar wind. The radial momentum equation for steady isothermal spherical expansion:

Using mass conservation $\dot{M} = 4\pi r^2 \rho v = \text{const}$, this transforms into the celebrated de Laval nozzle equation:

where $r_c = \frac{G M_\odot}{2 c_s^2}$ is the critical sonic radius. The physical solution starts subsonic ($v < c_s$) near the Sun, accelerates smoothly through the critical point $r = r_c$ where $v = c_s$, and becomes highly supersonic ($v \sim 400-800\text{ km s}^{-1}$) throughout the heliosphere.

4. Helioseismology

Convective turbulence excites millions of resonant acoustic standing waves ($p$-modes) that propagate throughout the solar interior and reflect at the surface. By Doppler-mapping surface velocity oscillations (periods $\sim 5\text{ minutes}$), helioseismologists invert the acoustic dispersion relation $\omega^2 = k_h^2 c_s^2$ to measure the internal sound speed profile $c_s(r) = \sqrt{\gamma P/\rho}$ to within $0.1\%$, mapping the base of the convective envelope ($0.713 R_\odot$) and validating the Standard Solar Model.