use Wien’s displacement law and the Stefan–Boltzmann law to estimate the radius of a star

Estimating the Radius of a Star using Wien’s Displacement Law and the Stefan–Boltzmann Law

Why this topic matters for Cambridge International AS & A Level Physics (9702)

  • AO1 – Knowledge & Understanding: recognise and state Wien’s law, the Stefan–Boltzmann law, the luminosity‑radius‑temperature relation and the definition of effective temperature.
  • AO2 – Application: combine observed peak wavelength, apparent magnitude and distance to calculate a star’s radius, including a bolometric correction.
  • AO3 – Analysis & Evaluation: discuss the black‑body assumption, interstellar extinction, limb darkening, propagate uncertainties and evaluate the reliability of the result.

Link to other syllabus sections

Syllabus topicConnection to this lesson
7 Waves – electromagnetic radiationWien’s law describes the wavelength at which a black‑body spectrum peaks; the relation c = λ f explains how spectrographs record λ.
5 Work, Energy & PowerStefan–Boltzmann law gives the radiant power per unit area (W m⁻²); multiplying by the star’s surface area yields the total power (luminosity).
14 TemperatureEffective temperature is obtained from the peak wavelength; it is the temperature a perfect black‑body would need to emit the observed total power.
12 Gravitation (A‑level)Once R is known, surface gravity follows from
g = GM / R², allowing a link to the mass‑radius relationship.
25 Astronomy & CosmologyThe trio (L, R, T) places a star on the Hertzsprung–Russell diagram (main‑sequence, giant, super‑giant).
22 Quantum Physics (A‑level)Black‑body radiation originates from Planck’s law; Wien’s law is the high‑temperature approximation used here.

Key constants (SI units)

ConstantSymbolValueUnits
Wien’s constantb2.898 × 10⁻³m·K
Stefan–Boltzmann constantσ5.670374419 × 10⁻⁸W·m⁻²·K⁻⁴
Solar luminosityL3.828 × 10²⁶W
Solar radiusR6.96 × 10⁸m
Solar absolute bolometric magnitudeMbol,⊙4.74mag

1. Theoretical background

1.1 Wave nature of light and spectral measurement

Electromagnetic radiation obeys c = λ f, where c is the speed of light, λ the wavelength and f the frequency. A spectrograph disperses the incoming light so that a detector records flux as a function of λ, producing the spectral‑energy distribution (SED) used in this unit.

1.2 Wien’s Displacement Law

For a perfect black‑body the wavelength at which the spectral radiance per unit wavelength is maximum is

\[

\lambda_{\max}= \frac{b}{T}

\]

where b = 2.898 × 10⁻³ m·K. A brief derivation (high‑temperature limit of Planck’s law) can be shown in a few steps:

  1. Planck’s law: \(B_{\lambda}(T)=\frac{2hc^{2}}{\lambda^{5}}\frac{1}{e^{hc/(\lambda kT)}-1}\).
  2. Set \(\frac{dB{\lambda}}{d\lambda}=0\) and solve for λ; the transcendental equation reduces to \(\lambda{\max}T \approx 2.898\times10^{-3}\) m·K for \(hc/(\lambda kT) \gg 1\).

Thus hotter stars have a shorter λmax.

1.3 Stefan–Boltzmann Law

The total radiant emittance (power per unit area) from a black‑body surface is

\[

j^{\star}=σT^{4}\qquad\bigl[\text{W·m}^{-2}\bigr]

\]

Multiplying by the surface area of a sphere gives the total luminosity:

\[

L = 4πR^{2}σT^{4}

\]

1.4 Effective temperature versus surface temperature

The temperature obtained from Wien’s law is the effective temperature – the temperature a perfect black‑body would need to emit the same total power as the star. It is not necessarily the local temperature at any point on the stellar surface, but the assumption of a uniform temperature is required for the simple \(L=4πR^{2}σT^{4}\) relation.

1.5 Surface gravity (link to gravitation)

When the stellar mass M is known, the surface gravity follows from Newton’s law:

\[

g = \frac{GM}{R^{2}}

\]

where G = 6.674 × 10⁻¹¹ N m² kg⁻². This provides a bridge to the A‑level gravitation topic.

1.6 Quantum background

Planck’s law (1900) solved the ultraviolet catastrophe and underpins the black‑body spectrum. Wien’s law is the high‑temperature (short‑wavelength) approximation of Planck’s law, which is why it is appropriate for most hot stars.

1.7 Assumptions & limitations (AO2/AO3)

  • Black‑body approximation: real stars have absorption lines, varying opacity and limb darkening.
  • Interstellar extinction: dust preferentially scatters blue light, shifting λmax to longer values and causing an underestimate of T.
  • Bolometric correction (BC): observations are usually made in a limited band (e.g., V). A BC converts a band‑limited absolute magnitude M to a bolometric magnitude Mbol. Typical values are tabulated by spectral type.
  • Geometric simplicity: the star is treated as a perfect sphere with a uniform temperature.

2. Propagating uncertainties (AO3)

If the measured peak wavelength and luminosity have uncertainties \(\Delta\lambda\) and \(\Delta L\), the fractional errors are

\[

\frac{\Delta T}{T}= \frac{\Delta\lambda}{\lambda_{\max}},\qquad

\frac{\Delta R}{R}= \frac12\left(\frac{\Delta L}{L}+4\frac{\Delta T}{T}\right)

\]

Consequences:

  • A 5 % error in λmax (hence T) contributes a 10 % error to R.
  • A 10 % error in L contributes only a 5 % error to R.

3. Step‑by‑step procedure (AO2)

  1. Obtain the spectrum (digital SED from a catalogue or your own observation). Note the instrument’s wavelength resolution.
  2. Determine the peak wavelength \(\lambda_{\max}\):

    • Fit a smooth curve (e.g., a Gaussian or a Planck curve) to the SED.
    • Read the wavelength of maximum flux; adopt \(\Delta\lambda =\) half the wavelength bin width or the fitting uncertainty.

  3. Calculate the effective temperature using Wien’s law:

    \[

    T = \frac{b}{\lambda_{\max}},\qquad

    \Delta T = T\frac{\Delta\lambda}{\lambda_{\max}}

    \]

  4. Find the bolometric luminosity \(L\):

    1. From the catalogue obtain the apparent magnitude \(m\) (usually V‑band) and the parallax \(\pi\) (mas). Convert to distance:

      \[

      d\;(\text{pc}) = \frac{1000}{\pi}

      \]

    2. Absolute magnitude:

      \[

      M = m - 5\log_{10}(d) + 5

      \]

    3. Apply the bolometric correction appropriate for the star’s spectral type (e.g., BC = –0.07 mag for a G2 V star). Then

      \[

      M_{\text{bol}} = M + \text{BC}

      \]

    4. Convert to luminosity relative to the Sun:

      \[

      \frac{L}{L{\odot}} = 10^{\;(M{\text{bol},\odot}-M_{\text{bol}})/2.5}

      \]

      where \(M_{\text{bol},\odot}=4.74\).

    5. Calculate the absolute luminosity:

      \[

      L = \left(\frac{L}{L{\odot}}\right)L{\odot}

      \]

      Propagate uncertainties from distance (parallax error) and photometric error to obtain \(\Delta L\).

  5. Compute the radius:

    \[

    R = \sqrt{\frac{L}{4\pi\sigma T^{4}}},\qquad

    \frac{\Delta R}{R}= \frac12\left(\frac{\Delta L}{L}+4\frac{\Delta T}{T}\right)

    \]

    Express the final answer in solar radii:

    \[

    \frac{R}{R_{\odot}} = \frac{R\;(\text{m})}{6.96\times10^{8}\ \text{m}}

    \]

  6. Optional – surface gravity (if the mass is known):

    \[

    g = \frac{GM}{R^{2}}

    \]

4. Practical activity (Paper 3/5 skill)

Goal: Using a real stellar spectrum, determine the star’s radius and write a concise evaluation of the result.

  1. Download a calibrated spectrum (flux vs. λ) for a bright star from an online database (e.g., SDSS, Gaia DR3).
  2. Using a spreadsheet or Python:

    • Plot the SED.
    • Locate \(\lambda_{\max}\) and record \(\Delta\lambda\) (half the wavelength bin width or fitting error).

  3. Calculate \(T\) and \(\Delta T\) with Wien’s law.
  4. From the same catalogue obtain:

    • Parallax \(\pi\) (mas) and its error.
    • Apparent V‑band magnitude \(m\) and its error.
    • Spectral type → read a suitable bolometric correction (BC) from a table.

    Follow the steps in §3‑4 to obtain \(L\), \(\Delta L\) and finally \(R\) and \(\Delta R\).

  5. Complete the Evaluation Checklist:

AspectComments / Evidence
Measurement uncertainties (λ, m, π)
Assumption of black‑body behaviour
Interstellar extinction (estimate E(B–V) if possible)
Bolometric correction accuracy
Possible systematic errors (limb darkening, instrumental response)
Suggested improvement (e.g., interferometry, eclipsing binary analysis)

Write a short (≈150 words) paragraph addressing the points above – this satisfies AO3.

5. Worked examples

5.1 The Sun (validation)

  • Peak wavelength (green light): \(\lambda_{\max}=5.00\times10^{-7}\ \text{m}\), \(\Delta\lambda=5.0\times10^{-9}\ \text{m}\).
  • Effective temperature:

    \[

    T_{\odot}= \frac{2.898\times10^{-3}}{5.00\times10^{-7}} = 5.80\times10^{3}\ \text{K},

    \qquad

    \frac{\Delta T}{T}=1\% \Rightarrow \Delta T\approx58\ \text{K}

    \]

  • Luminosity: \(L_{\odot}=3.828\times10^{26}\ \text{W}\) (negligible relative error for this illustration).
  • Radius:

    \[

    R_{\odot}= \sqrt{\frac{3.828\times10^{26}}{4\pi(5.670374419\times10^{-8})(5.80\times10^{3})^{4}}}

    = 6.96\times10^{8}\ \text{m}

    \]

    \[

    \frac{\Delta R}{R}= \frac12\bigl(0+4\times0.01\bigr)=2\%

    \Rightarrow \Delta R\approx1.4\times10^{7}\ \text{m}

    \]

  • Result: \(\displaystyle R_{\odot}= (6.96\pm0.14)\times10^{8}\ \text{m}\) – matches the accepted solar radius.

5.2 Betelgeuse (A‑level challenge)

Given:

  • \(\lambda_{\max}=1.0\times10^{-6}\ \text{m}\), \(\Delta\lambda=5.0\times10^{-8}\ \text{m}\) (5 % error).
  • Luminosity \(L = 1.2\times10^{5}\,L_{\odot}\) with \(\Delta L/L = 5\%\).

Calculations:

\[

T = \frac{2.898\times10^{-3}}{1.0\times10^{-6}} = 2.90\times10^{3}\ \text{K},

\qquad

\frac{\Delta T}{T}=5\%

\]

\[

R = \sqrt{\frac{1.2\times10^{5}\,L_{\odot}}{4\pi\sigma(2.90\times10^{3})^{4}}}

= 8.6\times10^{11}\ \text{m}

\]

\[

\frac{\Delta R}{R}= \frac12\bigl(0.05+4\times0.05\bigr)=0.125\;(12.5\%)

\]

\[

R = (8.6\pm1.1)\times10^{11}\ \text{m}

= (1\,240\pm160)\ R_{\odot}

\]

The huge radius places Betelgeuse firmly in the red‑supergiant region of the HR diagram.

6. Extension box – Connecting to the wider curriculum

Stellar classification & the HR diagram

The three observables (L, R, T) locate a star on the Hertzsprung–Russell diagram. Main‑sequence stars follow roughly \(L\propto M^{3.5}\) and have relatively small radii for their luminosities. Giants and super‑giants lie above the main sequence because their radii are vastly larger at comparable temperatures.

Black‑body radiation in the early universe

The Cosmic Microwave Background is a near‑perfect black‑body at 2.73 K. Using the Stefan–Boltzmann law one can calculate its energy density, linking this lesson to the quantum‑physics requirement (Planck’s law).

Other methods of measuring stellar radii

Interferometry (direct angular diameter → physical radius with distance).

Eclipsing binaries (light‑curve modelling gives relative radii).

Asteroseismology (oscillation frequencies constrain mean density → radius).

Discussing these methods reinforces AO3 – evaluation of experimental techniques.

7. Practice questions

  1. Betelgeuse – peak wavelength \(1.0\times10^{-6}\) m, luminosity \(1.2\times10^{5}\ L_{\odot}\).
    Estimate the radius in solar radii and comment on its position on the HR diagram.
  2. Exoplanet‑host star – \(\lambda{\max}=4.5\times10^{-7}\) m, apparent magnitude corresponding to \(0.85\ L{\odot}\).
    Calculate its radius and discuss whether it is likely to be on the main sequence.
  3. Explain qualitatively why, for a fixed luminosity, a hotter star must have a smaller radius than a cooler star. (Hint: use \(L\propto R^{2}T^{4}\).)
  4. Given uncertainties \(\Delta\lambda_{\max}=2\times10^{-9}\) m and \(\Delta L/L=0.07\), propagate the errors to find \(\Delta R/R\) for the star in question 2.

Suggested diagram

Sketch of a star showing radius R, effective temperature T, and a representative black‑body spectral‑energy distribution. The diagram should include arrows indicating the logical flow: observed spectrum → determine \(\lambda_{\max}\) → compute T (Wien) → obtain absolute (bolometric) magnitude → convert to luminosity L → calculate radius R (Stefan–Boltzmann).