The Standard Solar Model

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.713R_{\odot }$, and the entire interior run of pressure, temperature, density, and luminosity.¹

The energy generation is concentrated in the innermost core: the pp chain dominates, and because $ϵ\propto \rho X^2{T}^{\sim 4}$ for the pp chain, the luminosity is essentially complete by $0.25R_{\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.

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.²

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.713R_{\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.

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.³

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.³

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-\frac{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.³

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.⁴

Summary

The standard solar model integrates the structure equations for $1M_{\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.25R_{\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.713R_{\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 that power the core and set $ϵ$ are developed in the next module.