Data Analysis

Chapter 9 of 11 · 13 min read

Photometry toolkit

Magnitudes and fluxes, distance moduli, extinction and reddening, bolometric corrections, AB magnitudes and reddening-free indices.

Magnitudes and flux

Pogson's relation and its inverse:

m1−m2=−2.5log⁡10F1F2,F1F2=10−0.4(m1−m2).m_1 - m_2 = -2.5\log_{10}\frac{F_1}{F_2}, \qquad \frac{F_1}{F_2} = 10^{-0.4(m_1-m_2)}.

Magnitudes of unresolved stars do not add; their fluxes do:

m1+2=−2.5log⁡(10−0.4m1+10−0.4m2).m_{1+2} = -2.5\log\left(10^{-0.4m_1} + 10^{-0.4m_2}\right).

Going backwards, with a luminosity ratio L2/L1L_2/L_1 as in the 2023 eclipsing binaries,

m1=m1+2+2.5log⁡(1+L2/L1),m2=m1+2+2.5log⁡(1+L1/L2).m_1 = m_{1+2} + 2.5\log(1 + L_2/L_1), \qquad m_2 = m_{1+2} + 2.5\log(1 + L_1/L_2).

Distance modulus and extinction

m−M=5log⁡10dpc−5+A,d=10(m−M−A+5)/5 pc.m - M = 5\log_{10}d_{\rm pc} - 5 + A, \qquad d = 10^{(m - M - A + 5)/5}\ {\rm pc}.

Dust dims and reddens. The colour excess is E(B−V)=(B−V)−(B−V)0E(B-V) = (B-V) - (B-V)_0, with the intrinsic colour from the spectral type and luminosity class, and the extinction is

AV=RVE(B−V),RV≈3.1A_V = R_VE(B-V), \qquad R_V\approx3.1

in the diffuse Galactic ISM. Since EX−V=AX−AVE_{X-V} = A_X - A_V and AX→0A_X\to0 at long wavelengths, plotting EX−V/EB−VE_{X-V}/E_{B-V} against 1/λ1/\lambda and extrapolating to 1/λ=01/\lambda = 0 gives −RV-R_V (2015 P2).

If a source is seen in two bands with known absolute magnitudes, the difference of the two apparent distance moduli is a colour excess: μR−μI=AR−AI=E(R−I)\mu_R - \mu_I = A_R - A_I = E(R-I).

Bolometric quantities

The bolometric correction converts a band magnitude to the total: mbol=mV+BCm_{\rm bol} = m_V + BC, with BC from a table against temperature (most negative for very hot and very cool stars). Then

LL⊙=10−0.4(Mbol−Mbol,⊙),Mbol,⊙≈4.74,\frac{L}{L_\odot} = 10^{-0.4(M_{\rm bol} - M_{\rm bol,\odot})}, \qquad M_{\rm bol,\odot}\approx4.74,

and L=4πR2σT4L = 4\pi R^2\sigma T^4 links luminosity, radius and temperature. Temperatures come from colours through a calibration curve, which is how the 2015 Baade–Wesselink problem gets TT from V−RV-R.

AB magnitudes and spectral slopes

The AB system is defined by a flat fνf_\nu:

mAB=−2.5log⁡fν3631 Jy.m_{\rm AB} = -2.5\log\frac{f_\nu}{3631\ {\rm Jy}}.

Converting between frequency and wavelength densities uses fν ∣dν∣=fλ ∣dλ∣f_\nu\,|d\nu| = f_\lambda\,|d\lambda|, so fλ=fνc/λ2f_\lambda = f_\nu c/\lambda^2. A source with constant mABm_{\rm AB} therefore has fλ∝λ−2f_\lambda\propto\lambda^{-2} (β=−2\beta = -2).

For distant galaxies, observed wavelengths are (1+z)(1+z) times the rest-frame ones. Track every (1+z)(1+z) factor explicitly: in 2018 D1.4 one missing factor of 7.6 changes the conclusion.

Reddening-free indices

When the extinction of each star is unknown, combine two bands so that it cancels. The Wesenheit magnitude

WVI=V−RVI(V−I),RVI=AVE(V−I),W_{VI} = V - R_{VI}(V - I), \qquad R_{VI} = \frac{A_V}{E(V-I)},

is unchanged by reddening that follows the assumed law: dimming moves VV and V−IV-I together exactly so that WW stays fixed. In 2017 D1 the Wesenheit PL relation gives LMC distances with a scatter of 0.9 kpc against 4.1 kpc in V.

The surface brightness SV=mV+5log⁡θS_V = m_V + 5\log\theta is a similar trick. It depends almost only on temperature, and therefore on colour, so a colour predicts the angular diameter.

Atmospheric extinction

From the ground, m=m0+kXm = m_0 + kX, with airmass X≈sec⁡zX\approx\sec z for zenith distance z<60∘z<60^\circ. Plot mm against XX for one star over a night: the slope is the extinction coefficient kk and the intercept at X=0X = 0 is the magnitude above the atmosphere.

Key formulas

Pogson
m1−m2=−2.5log⁡10F1F2m_1 - m_2 = -2.5\log_{10}\frac{F_1}{F_2}
Distance modulus
m−M=5log⁡10dpc−5+Aλm - M = 5\log_{10} d_{\rm pc} - 5 + A_\lambda
Reddening
AV=RV E(B−V),RV≈3.1A_V = R_V\,E(B-V), \qquad R_V \approx 3.1
Luminosity from bolometric magnitude
LL⊙=10−0.4 (Mbol−Mbol,⊙)\frac{L}{L_\odot} = 10^{-0.4\,(M_{\rm bol} - M_{\rm bol,\odot})}
Wesenheit index
WVI=V−RVI (V−I)W_{VI} = V - R_{VI}\,(V - I)
Surface brightness
SV=mV+5log⁡10θmas  ⇒  θ=100.2(SV−mV) masS_V = m_V + 5\log_{10}\theta_{\rm mas} \;\Rightarrow\; \theta = 10^{0.2(S_V - m_V)}\ {\rm mas}

Common mistakes

  • Adding magnitudes instead of fluxes.
  • Forgetting extinction in the distance modulus, or subtracting it with the wrong sign.
  • Using the main-sequence intrinsic colour for a supergiant.
  • Dropping a (1 + z) factor when moving between observed and rest frames.

Practise this chapter

Past IOAA problems that use these skills.

Saved in this browser only.