---
title: The Standard Solar Model
draft: false
module: Stellar Structure
moduleNumber: 4
lessonNumber: 4
order: 404
summary: >
  The standard solar model integrates the structure equations for one solar mass
  and calibrates the composition and convection parameter to reproduce the Sun's
  observed luminosity, radius, and age. Helioseismology tests the model's sound
  speed through the Sun's acoustic p-mode oscillations, and the model predicts a
  neutrino flux by production channel. The measured deficit, the solar-neutrino
  problem, is resolved by matter-enhanced flavor oscillation, confirmed when SNO
  measured the total flux across all flavors.
topics: [Stellar Structure]
sources:
  - book: Carroll & Ostlie
    ref: "Ch. 11 — The Sun; §11.1 The Solar Interior, §11.2 The Solar Atmosphere"
  - book: Maoz
    ref: "Ch. 3 — Stellar Physics"
  - book: PDG
    ref: "Review of Particle Physics — Neutrino Masses, Mixing, and Oscillations"
---

The Sun is the one star whose interior can be tested against theory in detail: its
mass, radius, luminosity, age, and surface composition are all known to high
precision, its oscillation spectrum probes the sound speed at every depth, and its
core produces a neutrino flux that arrives at Earth in eight minutes. The **standard
solar model** is the result of integrating the four structure equations for one solar
mass with the best available microphysics, calibrated so that a model of the Sun's age
reproduces the observed luminosity and radius. This lesson describes how the model is
built and calibrated, how **helioseismology** tests it through the Sun's acoustic
oscillations, how it predicts the **solar-neutrino** flux channel by channel, and how
the long-standing deficit between prediction and measurement — the **solar-neutrino
problem** — was resolved by neutrino flavor oscillation.

## Building and calibrating the standard solar model

The standard solar model is not a fit but a constrained integration. The total mass
is fixed at $M_\odot = 1.989 \times 10^{30}\ \text{kg}$, and the structure equations
are integrated from an initial zero-age composition, then evolved forward to the solar
age $t_\odot = 4.57 \times 10^9\ \text{yr}$. Three inputs are adjusted so that the
evolved model matches three observations.

- **Initial helium abundance $Y_0$** is tuned so the model reaches the observed
  luminosity $L_\odot = 3.828 \times 10^{26}\ \text{W}$ at the solar age. Helium has no
  strong spectral lines at photospheric temperatures, so $Y_0$ cannot be measured
  directly and must be calibrated.
- **Mixing-length parameter $\alpha$** is tuned so the model reaches the observed radius
  $R_\odot = 6.957 \times 10^8\ \text{m}$; $\alpha$ controls the efficiency of
  convection in the outer envelope, which sets how far the model puffs out.
- **Initial metal abundance $Z_0$** is set from the photospheric metal-to-hydrogen
  ratio $Z/X$ measured spectroscopically, corrected for the gravitational settling of
  heavy elements over the Sun's lifetime.

The calibrated model is then a genuine prediction for every other quantity: the
central temperature $T_c \approx 1.57 \times 10^7\ \text{K}$, the central density
$\rho_c \approx 1.5 \times 10^5\ \text{kg}\,\text{m}^{-3}$, the central pressure
$P_c \approx 2.34 \times 10^{16}\ \text{Pa}$, the depth of the convection zone at
$0.713\,R_\odot$, and the entire interior run of pressure, temperature, density, and
luminosity.[^co-interior]

The energy generation is concentrated in the innermost core: the pp chain dominates,
and because $\epsilon \propto \rho X^2 T^{\sim 4}$ for the pp chain, the luminosity is
essentially complete by $0.25\,R_\odot$, inside which about half the solar mass and
almost all the nuclear burning reside. The temperature falls from $1.57 \times 10^7\
\text{K}$ at the center to about $5772\ \text{K}$ at the photosphere, and the density
falls by ten orders of magnitude over the same span.

$$
% caption: The calibrated solar model: energy generation is confined to the inner
% core, so the luminosity saturates within a quarter radius while T, rho, P fall outward.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,0) -- (8.4,0) node[right, black!70] {radius fraction};
\draw[->, black] (0,0) -- (0,4.4) node[above, black!70] {normalized value};
\foreach \x/\lab in {0/0,2/0.25,4/0.5,6/0.75,8/1.0}
  \draw[black] (\x,0.06) -- (\x,-0.06) node[below, black] {\lab};
% energy generation rate: sharp central peak
\draw[acc, very thick] (0.1,4.0) .. controls (0.7,3.6) and (1.4,0.9) .. (2.2,0.3)
  .. controls (3.0,0.08) and (4.0,0.03) .. (8.0,0.02);
\node[acc, anchor=north west] at (0.5,3.9) {energy generation};
% luminosity: rises and saturates by 0.25 R
\draw[acc, very thick, densely dotted] (0.1,0.1) .. controls (1.0,3.0) and (1.9,3.85) .. (3.2,3.95)
  -- (8.0,3.95);
\node[acc, anchor=south] at (5.0,3.6) {luminosity};
% temperature: falls monotonically
\draw[black, very thick, densely dashed] (0.1,3.6) .. controls (2.4,2.4) and (4.6,1.0) .. (8.0,0.2);
\node[black, anchor=west] at (2.5,1.4) {temperature};
% convection zone marker at 0.713 R -> x = 5.7
\draw[black, densely dashed] (5.7,0) -- (5.7,3.2);
\node[black, anchor=south, rotate=90] at (5.5,1.8) {convection zone base};
\end{tikzpicture}
$$

## Helioseismology and the p-mode oscillations

The Sun oscillates. Its surface shows a superposition of millions of standing acoustic
waves with periods clustered near five minutes, detected as Doppler shifts of
photospheric lines that oscillate at frequencies around $\nu \approx 3\ \text{mHz}$.
These are **p-modes**, pressure (sound) waves trapped in the solar interior, and their
frequencies encode the internal structure the way the tone of an organ pipe encodes its
length and the speed of sound inside it.

A sound wave launched downward from the surface refracts: the sound speed $c_s =
\sqrt{\Gamma_1 P/\rho}$ rises inward as the temperature rises, so the wave bends back
toward the surface at a **lower turning point** whose depth increases with the
horizontal wavelength of the mode. Above, the wave reflects near the surface where the
density scale height becomes comparable to the wavelength. Each mode is thus trapped in
a resonant cavity between an inner turning point and the surface, and only discrete
frequencies fit. Low-degree modes penetrate to the core; high-degree modes are confined
to the outer layers. Measuring thousands of mode frequencies and inverting the
resonance conditions reconstructs the sound-speed profile $c_s(r)$ throughout the
interior.[^co-helio]

The inversions confirm the standard solar model's sound speed to better than $0.5\%$
through most of the radius, pin the base of the convection zone at $r = 0.713\,R_\odot$,
and measure the near-uniform rotation of the radiative interior. The one persistent
discrepancy, the **solar abundance problem**, is a $\sim 1\%$ sound-speed mismatch just
below the convection zone that appeared when revised (lower) photospheric metal
abundances were adopted, and it remains unresolved. Helioseismology promoted the solar
model from a plausible construction to a quantitatively tested one.

$$
% caption: A p-mode is trapped between an inner turning point, where rising sound speed
% refracts the wave back up, and the surface; deeper modes probe the core.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
% solar disk
\draw[black] (0,0) circle (3.0);
\node[black, anchor=south] at (0,3.05) {surface};
% shallow ray cavity (high degree)
\draw[acc, very thick] (60:3.0) arc (140:220:1.5) ;
\node[acc, anchor=west] at (1.7,1.9) {shallow mode};
% deep ray cavity reaching near core (low degree)
\draw[acc!70, very thick, densely dashed] (200:3.0) .. controls (-1.5,-1.0) and (0.0,-0.4) .. (0.4,-0.2)
  .. controls (0.9,0.0) and (1.6,-0.4) .. (250:3.0);
\node[acc!80, anchor=north] at (0,-0.9) {deep mode probes core};
% inner turning point marker
\fill[black] (0.4,-0.2) circle (2pt);
\node[black, anchor=south west] at (0.45,-0.2) {turning point};
\end{tikzpicture}
$$

## The solar-neutrino flux prediction

Every completed pp-chain fuses four protons into one helium-4 nucleus, releasing two
electron neutrinos. The neutrinos escape the core immediately and stream to Earth
essentially unabsorbed, carrying a direct measurement of the nuclear reactions in the
core at the moment they occur — unlike photons, which take $\sim 10^5$ years to diffuse
out. The total flux is fixed by the luminosity: $L_\odot$ requires a definite
proton-fusion rate, hence a definite neutrino production rate, giving a predicted flux at
Earth of about $6 \times 10^{10}\ \text{neutrinos}\,\text{cm}^{-2}\,\text{s}^{-1}$.

The neutrinos arrive in distinct spectral components, one per production reaction, with
sharply different energies:

- **pp neutrinos** from $p + p \to d + e^+ + \nu_e$: a continuous spectrum up to
  $0.42\ \text{MeV}$, the dominant flux ($\sim 91\%$), fixed almost model-independently
  by the luminosity.
- **$^7$Be neutrinos** from electron capture on beryllium-7: two monoenergetic lines at
  $0.86$ and $0.38\ \text{MeV}$.
- **$^8$B neutrinos** from $^8\text{B} \to {}^8\text{Be} + e^+ + \nu_e$: a continuous
  spectrum to $\sim 15\ \text{MeV}$, a tiny fraction of the flux ($\sim 0.01\%$) but the
  most temperature-sensitive, scaling as roughly $T_c^{24}$, so the high-energy $^8$B
  flux is the sharpest probe of the central temperature.

The high energy of the $^8$B neutrinos made them the first to be detected, in the
Homestake chlorine experiment and later in the water-Cherenkov detectors Kamiokande and
Super-Kamiokande, while the gallium experiments SAGE and GALLEX reached down to the
dominant low-energy pp flux.[^pdg-nu]

$$
% caption: The predicted solar-neutrino spectrum by production channel: the dominant
% low-energy pp continuum, the Be-7 lines, and the rare high-energy B-8 tail.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,0) -- (8.6,0) node[right, black!70] {log energy in MeV};
\draw[->, black] (0,0) -- (0,4.2) node[above, black!70] {log rate};
\foreach \x/\lab in {0.6/0.1,3.0/1.0,5.4/3.0,7.8/15}
  \draw[black] (\x,0.06) -- (\x,-0.06) node[below, black] {\lab};
% pp continuum: high flux, low energy, dome shape ending at 0.42 MeV
\draw[acc, very thick] (0.4,2.4) .. controls (1.4,3.6) and (2.0,3.6) .. (2.6,2.0) -- (2.6,0.05);
\node[acc, anchor=south] at (1.4,3.4) {pp};
% Be-7 lines (two vertical lines)
\draw[black, very thick] (2.2,0.05) -- (2.2,2.3);
\draw[black, very thick] (3.0,0.05) -- (3.0,2.7);
\node[black, anchor=south] at (3.0,2.75) {Be-7};
% B-8 continuum: low flux, high energy, extends far right
\draw[black, very thick, densely dashed] (3.6,0.4) .. controls (5.0,1.0) and (6.4,1.0) .. (7.8,0.1);
\node[black, anchor=south] at (6.0,1.0) {B-8};
\end{tikzpicture}
$$

## The solar-neutrino problem and its resolution

Every experiment measured fewer neutrinos than the standard solar model predicted. The
chlorine experiment saw about one-third of the predicted rate; the water-Cherenkov
detectors saw about half of the $^8$B flux; the gallium experiments saw about $60\%$ of
the low-energy flux. The deficit was robust, energy-dependent, and persisted for three
decades. This was the **solar-neutrino problem**: either the solar model overpredicted
the core temperature and reaction rates, or something happened to the neutrinos in
transit.

Because the $^8$B flux scales as $T_c^{24}$, a solar-model explanation required lowering
$T_c$ by only $\sim 5\%$ to cut the $^8$B flux in half — but helioseismology measured the
core sound speed, hence the temperature, in agreement with the standard model to well
under a percent, closing off that escape. The alternative is that neutrinos change flavor
in flight. The detectors above were sensitive mainly or only to electron neutrinos
$\nu_e$; if the $\nu_e$ produced in the core converted partly into $\nu_\mu$ and
$\nu_\tau$ before reaching Earth, the electron-flavor count would fall below the total.

The resolution was confirmed by the **Sudbury Neutrino Observatory (SNO)**, a heavy-water
detector that measured the $^8$B flux through two channels simultaneously:

- the **charged-current** reaction $\nu_e + d \to p + p + e^-$, sensitive only to
  electron neutrinos;
- the **neutral-current** reaction $\nu_x + d \to p + n + \nu_x$, equally sensitive to all
  three flavors.

The charged-current rate recovered only about a third of the predicted flux, matching the
earlier deficit; but the neutral-current rate, counting all flavors, recovered the **full**
standard-solar-model flux within errors. The neutrinos were not missing — two-thirds of
them had changed flavor. The solar model was vindicated, and the missing neutrinos became
direct evidence that neutrinos have mass and mix between flavors.[^pdg-nu]

## Matter-enhanced flavor conversion

The conversion is stronger than simple vacuum oscillation because of the **MSW effect**:
electron neutrinos acquire an extra effective mass in matter through their coherent
forward scattering off electrons, which the other flavors do not experience. In the dense
solar core the matter term dominates the flavor evolution; as a neutrino travels outward
through the falling electron density, it passes adiabatically through a resonance that
converts a high-energy $\nu_e$ almost entirely into a heavier mass eigenstate. This leaves
a strongly **energy-dependent** survival probability:

- **Low-energy pp neutrinos** ($< 1\ \text{MeV}$) are below the MSW resonance and undergo
  ordinary vacuum-averaged oscillation, surviving with probability $P_{ee} \approx 1 -
  \tfrac{1}{2}\sin^2 2\theta_{12} \approx 0.55$.
- **High-energy $^8$B neutrinos** ($> 5\ \text{MeV}$) pass through the matter resonance and
  survive with the smaller probability $P_{ee} \approx \sin^2\theta_{12} \approx 0.3$.

The transition between the two regimes near a few MeV is the signature prediction of the
MSW mechanism, and its measurement across the pp, $^7$Be, and $^8$B energies matches the
mixing parameters $\sin^2\theta_{12} \approx 0.307$ and $\Delta m_{21}^2 \approx 7.5
\times 10^{-5}\ \text{eV}^2$ determined independently by reactor experiments.[^pdg-nu]

$$
% caption: The electron-neutrino survival probability drops from the vacuum value at low
% energy to the matter-dominated value above the MSW resonance near a few MeV.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,0) -- (8.6,0) node[right, black!70] {neutrino energy in MeV};
\draw[->, black] (0,0) -- (0,4.2) node[above, black!70] {survival probability};
\foreach \y/\lab in {0/0,2.2/0.5,4.0/0.9}
  \draw[black] (-0.08,\y) -- (0.08,\y) node[left, black, xshift=-3pt] {\lab};
\foreach \x/\lab in {0.6/0.1,3.0/1.0,5.4/5,7.8/15}
  \draw[black] (\x,0.06) -- (\x,-0.06) node[below, black] {\lab};
% survival probability curve: ~0.55 low, dropping to ~0.3 high
\draw[acc, very thick] (0.4,2.5) -- (2.6,2.45) .. controls (3.6,2.35) and (4.6,1.5) .. (5.6,1.35)
  -- (8.0,1.3);
\node[acc, anchor=south] at (1.6,2.55) {pp: vacuum};
\node[acc, anchor=south] at (6.6,1.4) {B-8: matter};
% resonance region marker
\draw[black, densely dashed] (4.6,0) -- (4.6,3.4);
\node[black, anchor=south, rotate=90] at (4.4,1.9) {MSW resonance};
\end{tikzpicture}
$$

## The solar atmosphere: granulation and the activity cycle

Above the convection zone the model connects to the observable **solar atmosphere**. The
top of the convection zone breaks the surface as **granulation**: a shifting pattern of
bright cells about $1000\ \text{km}$ across, the tops of convective upflows, ringed by
dark lanes of cooler descending gas, each cell lasting minutes. The granulation is the
direct surface signature of the convection that mixing-length theory parametrizes.

Overlaid on this is the **magnetic activity cycle**. The differential rotation of the
convection zone winds and amplifies the Sun's magnetic field, which erupts through the
surface as sunspots, cool ($\sim 3800\ \text{K}$) magnetically suppressed regions that
appear in an eleven-year cycle, migrating from mid-latitudes toward the equator as the
cycle proceeds. The field reverses polarity each cycle, giving a full 22-year magnetic
period. The activity cycle modulates the ultraviolet output, the solar wind, and the flare
and coronal-mass-ejection rate, coupling the calibrated interior model to the space-weather
environment of the Solar System.[^co-atmos]

## Summary

The standard solar model integrates the structure equations for $1\,M_\odot$, tuning the
initial helium and metal abundances and the mixing-length parameter to reproduce
$L_\odot$, $R_\odot$, and the surface composition at the solar age, and it then predicts a
central temperature of $1.57 \times 10^7\ \text{K}$ with the luminosity generated inside
$0.25\,R_\odot$. Helioseismology tests the model through the Sun's p-mode acoustic
oscillations, confirming the sound speed to under a percent and fixing the convection-zone
base at $0.713\,R_\odot$. The model's neutrino-flux prediction, channel by channel,
appeared deficient in every early experiment, but SNO's simultaneous measurement of the
electron-flavor and all-flavor $^8$B fluxes showed the total matched the model while
two-thirds had converted flavor, resolving the solar-neutrino problem through
matter-enhanced oscillation and confirming neutrino mass. This closes the stellar-structure
module; the [nuclear reactions](/astrophysics-cosmology/nuclear-astrophysics/thermonuclear-reaction-rates-and-the-gamow-peak)
that power the core and set $\epsilon$ are developed in the next module.

[^co-interior]: Carroll & Ostlie, §11.1 — The Solar Interior: the calibrated standard solar model, the central conditions, and the interior structure.
[^co-helio]: Carroll & Ostlie, §11.1 — helioseismology, the p-mode oscillations, and the inversion for the interior sound speed and convection-zone depth.
[^co-atmos]: Carroll & Ostlie, §11.2 — The Solar Atmosphere: granulation, sunspots, and the magnetic activity cycle.
[^pdg-nu]: Particle Data Group, Review of Particle Physics — Neutrino Masses, Mixing, and Oscillations: the solar-neutrino flux, the MSW matter effect, and the mixing parameters. https://pdg.lbl.gov