Flame Photometry and Atomic Emission

A hospital lab draws a patient's blood and dilutes it a thousandfold. It then sprays the sample into a flame no hotter than a kitchen gas burner turned up high. Yet within seconds a photodetector reports the serum sodium concentration. That result is accurate to within 1–2 mmol/L of the true value. This measurement runs millions of times a day worldwide. It rests on one of the oldest quantitative tools in analytical chemistry: flame photometry. This is the emission-based twin of Atomization Methods for AAS. Atomic absorption asks how much light a cloud of ground-state atoms removes from a beam. Flame emission asks a different question instead. It measures how much light the small thermally-excited fraction of those atoms gives off on its own. No external lamp is needed at all. The technique matters because it is cheap, fast, and uncannily well suited to certain elements. Those elements — sodium, potassium, lithium, calcium — are hardest to excite by other means without a hot plasma. They also dominate clinical chemistry, agronomy, and water-quality testing. A propane–air or acetylene–air flame at 1700–2400 K is far too cool to excite most metals appreciably. But the alkali metals have such low-lying resonance levels that this changes. Sodium's 3p state, for example, sits only 2.10 eV above ground. As a result, a measurable, reproducible, temperature-sensitive population glows visibly. This makes flame photometry a temperature-dependent instrument in a way its absorption cousin is not. That dependence is both its greatest analytical strength and its most persistent source of interference. The big picture is a chain of four physical steps. First, nebulization and desolvation turn a liquid sample into a fine aerosol. This aerosol then becomes a dry particulate stream. Second, atomization in the flame's inner cone breaks that particulate into free atoms. Third, thermal (collisional) excitation promotes a small fraction of those atoms to an excited electronic state. This fraction is governed by Boltzmann statistics. Fourth, spontaneous emission converts that excited population back to photons. These photons appear at the element's characteristic wavelength. Examples include 589.0/589.6 nm for sodium and 766.5/769.9 nm for potassium. Lithium emits at 670.8 nm, and calcium at 422.7 nm. A monochromator or narrow-band interference filter isolates that line. A photomultiplier or photodiode then converts photon flux into current. This article works through the excitation physics governing that glowing fraction. It also covers the self-absorption and ionization effects that curve flame-emission calibration lines at real analytical concentrations. It explains the internal-standard method that made flame photometry clinically trustworthy. It then presents three fully worked numerical examples. These range from a Boltzmann population calculation to a serum-sodium determination with a lithium internal standard. The third is an exam-style self-absorption correction. It closes by placing the flame photometer in its modern context. There, it appears as the low-cost ancestor of Inductively Coupled Plasma Spectrometry.

Key takeaways
  • Flame emission measures light spontaneously radiated by the small, thermally excited fraction of atoms in a flame — no external lamp needed, unlike AAS.
  • The excited-state population follows the Boltzmann distribution; because ΔE/kT\Delta E/kT is large, even a 200 K temperature drift can change intensity by several percent for resonance lines.
  • Alkali metals (Na, K, Li, Rb, Cs) are ideal analytes because their lowest excited states lie at 1.6–3.0 eV, giving detectable population even in a cool (1700–2200 K) flame.
  • Self-absorption — reabsorption of emitted photons by cool, ground-state atoms at the flame's edge — bends the calibration curve downward (sublinear) at high concentration.
  • Ionization interference suppresses emission (and absorption) by converting neutral atoms to ions; an ionization suppressor (e.g., excess CsCl) buffers the electron density and restores linearity.
  • The internal-standard method (commonly Li or Cs) cancels nebulization, flame-temperature, and viscosity fluctuations by referencing analyte emission to a fixed-concentration reference element run through the same flame.
  • Key idea: flame emission intensity is a thermally-amplified, temperature-exponential readout of atom density — powerful for easily excited elements, but demanding of internal standardization and ionization control because the flame itself is both the atom source and the excitation source.

Where this lives today

Flame photometry's clinical role has partly receded before ion-selective electrodes for routine serum sodium/potassium in large hospital analyzers, but it remains a reference and backup method in many clinical laboratories worldwide, valued for its simplicity, low reagent cost, and lack of electrode-membrane drift or fouling issues. In agriculture and soil science it remains a workhorse: exchangeable potassium and sodium in soil extracts, and potassium/sodium content of fertilizers, plant tissue digests, and irrigation water, are still routinely measured by flame photometry in national and university soil-testing labs because the method needs no elaborate sample preparation beyond dilution and no expensive plasma source. Brewing and food-science labs use flame photometry to monitor sodium and potassium in beverages and processed foods for labeling compliance. Environmental and municipal water-testing labs use it as a rapid screen for alkali and alkaline-earth content in drinking water and wastewater, often alongside the sister technique of flame atomic absorption for elements that emit too weakly (or ionize too readily) in a simple flame to be measured reliably by emission alone. Pharmaceutical quality control uses flame photometry for electrolyte content verification in intravenous fluids and injectable formulations, where sodium, potassium, and calcium levels must be certified precisely against pharmacopeial limits. In every one of these settings, the underlying physics is unchanged from Bunsen and Kirchhoff's original apparatus: a clean, stable flame; a monochromatic isolation of one atomic line; a photodetector; and, in every modern instrument, an internal-standard channel running simultaneously to cancel the flame's inherent thermal noise. The technique's conceptual descendants — instruments that trade the chemical flame for a much hotter, more stable, argon-sustained radiofrequency plasma — extend the same Boltzmann-population logic to dozens of additional elements with far lower detection limits and far less self-absorption, which is the direction taken up in Inductively Coupled Plasma Spectrometry.

Ask the tutor: How would the self-absorption turnover concentration cmax⁡c_{\max} change if I switched from an air–acetylene flame to a hotter nitrous oxide–acetylene flame for the same analyte?

The physical picture — thorough prose

A flame photometer looks deceptively similar to an atomic absorption spectrophotometer with the hollow-cathode lamp switched off — and that is almost literally what it is. The sample solution is aspirated by the Venturi effect of the oxidant flow into a nebulizer, producing an aerosol of droplets typically 1–10 µm in diameter. A spray chamber removes the larger droplets (which would otherwise never fully desolvate and atomize), so that only the finest mist — a small, size-selected fraction of the aspirated sample, often just 5–15% by mass — reaches the burner head. There, in the flame's interconal zone, solvent evaporates, the resulting salt particle vaporizes, and molecular species dissociate into free gas-phase atoms, exactly as described for atomizer design generally. What happens next is where emission diverges from absorption. In AAS, the free atoms sit overwhelmingly in the ground state (a Boltzmann calculation shows the excited fraction is often only 1 part in 10410^4–101010^{10}), and a lamp beam is passed through the population to measure how much resonance light is absorbed. In flame emission, no lamp is used at all: the flame's own thermal energy — random collisions among atoms, ions, and combustion-product molecules at 1700–3200 K — populates the excited state directly. A tiny but non-negligible fraction of atoms, having been collisionally kicked upstairs, spontaneously fall back down within nanoseconds, each emitting one photon at the transition's characteristic wavelength. Because collisional excitation and de-excitation are both occurring continuously and in dynamic (though not full thermodynamic) equilibrium, the excited-state population is well described by the Boltzmann law, and the emitted intensity is proportional to it. The exponential temperature-dependence built into the Boltzmann factor, e−ΔE/kTe^{-\Delta E/kT}, is the single most important feature distinguishing flame emission from absorption spectrometry. Absorption signal depends on ground-state atom density, which is nearly temperature-independent for cool flames (ground state holds >99.99% of atoms regardless of small temperature shifts). Emission signal depends on excited-state density, which is exponentially sensitive to TT. A flame that drifts from 2200 K to 2250 K — an entirely plausible fluctuation from fuel/oxidant flow-rate noise — can shift a sodium emission signal by several percent, even though essentially nothing has happened to the ground-state (absorption) population. This is why flame emission historically required either very stable flames, an internal standard, or both, while flame AAS tolerated modest flame drift much better.

Figure
Energy-level diagram: sodium's 3s→3p resonance transition and thermal excitation route
Na atom, 3s ground state (0 eV) 3p excited state, ΔE = 2.10 eV (589 nm) thermal (collisional) excitation, ~kT ≈ 0.19 eV at T = 2200 K spontaneous emission 589.0 / 589.6 nm doublet (D2 / D1 lines) N₂/N₀ = (g₁/g₀) exp(−ΔE/kT) ≈ 1.0 × 10⁻³ at 2200 K Only ~0.1% of Na atoms are excited at any instant — yet the emitted glow is bright

The second key physical feature is the flame's spatial temperature and concentration structure. Immediately above the burner tip sits the primary combustion zone (a thin blue cone), then the interconal zone — hottest, most chemically reducing, and where atomization and excitation are most efficient — and finally the outer, cooler secondary combustion zone, rich in oxidizing combustion products and still containing unexcited, ground-state analyte atoms. A photon emitted deep in the hot core, on its way out of the flame toward the detector, can pass through this cooler, ground-state-rich outer mantle and be reabsorbed by identical atoms sitting in their ground state — precisely the resonance-absorption process of AAS, now acting as a loss mechanism on the emitted signal. This is self-absorption, and at high analyte concentration it becomes severe enough to visibly curve, and even reverse ("self-reversal"), the calibration line.

Figure
Self-absorption: emission from the hot core reabsorbed by cool outer atoms
hot core, T≈2400 K excited atoms emit hν cooler mantle, T≈1900 K ground-state atoms — reabsorb hν photon reabsorbed here burner head — nebulized aerosol enters from below

Self-absorption is concentration-dependent because the optical path length of ground-state atoms the emitted photon must traverse scales with atom density. At trace concentrations the outer mantle is nearly transparent and self-absorption is negligible; at high concentrations (approaching or exceeding roughly 10210^2–10310^3 mg/L for sodium, depending on flame geometry) the calibration curve visibly bends toward the concentration axis, exactly as Absorption and Emission Processes predicts for any optically thick emitting-and-absorbing medium — this is the emission analog of deviations from Beer's law at high absorbance.

The core ideas — the governing equations/relationships

The Boltzmann population ratio between an excited state 1 and the ground state 0 is

N1N0=g1g0 e−ΔE/kT\frac{N_1}{N_0} = \frac{g_1}{g_0}\, e^{-\Delta E/kT}

where: N1,N0N_1, N_0 are number densities in the excited and ground states, g1,g0g_1, g_0 are their statistical (degeneracy) weights, ΔE\Delta E is the energy gap between the two states (J or eV), kk is Boltzmann's constant (8.617×10−58.617\times10^{-5} eV/K), and TT is the flame temperature in kelvin. Reading: because ΔE\Delta E for a resonance line like Na 3s→3p is on the order of 2 eV while kTkT at flame temperature is only ~0.15–0.28 eV, the exponent is large and negative — the ratio is always small, but it is exponentially sensitive to TT. A doubling of TT can change N1/N0N_1/N_0 by orders of magnitude; a few percent change in TT changes it by several percent, which is the central practical fact of flame emission work. The emitted radiant power (intensity) from spontaneous decay of the excited population is

Iem=hν A10 N1 VI_{\text{em}} = h\nu\, A_{10}\, N_1\, V

where: hνh\nu is the photon energy of the transition, A10A_{10} is the Einstein spontaneous-emission coefficient (s⁻¹) for the 1→01\to0 transition, N1N_1 is the excited-state number density, and VV is the effective emitting volume viewed by the optics. Reading: combining the two relations above with the fact that at flame conditions N0≈NtotalN_0 \approx N_{\text{total}} (essentially all atoms remain in the ground state) gives the working analytical relation:

Iem∝Ntotal g1g0 e−ΔE/kTI_{\text{em}} \propto N_{\text{total}}\, \frac{g_1}{g_0}\, e^{-\Delta E/kT}

Since NtotalN_{\text{total}} (total free-atom density in the flame) is, at low concentration and before self-absorption sets in, directly proportional to analyte concentration cc in the aspirated solution, this yields the empirical Lomakin–Scheibe relation widely used to describe flame- and plasma-emission calibration curves:

Iem=a c bI_{\text{em}} = a\, c^{\,b}

where: aa is a proportionality constant absorbing atomization efficiency, transition probability, and instrumental sensitivity, cc is analyte concentration, and bb is an empirical exponent, ideally b=1b=1 at low concentration (linear regime) but falling below 1 (self-absorption) or occasionally exceeding 1 (partial ionization suppression relief) at higher concentration. Taking logarithms gives log⁡Iem=log⁡a+blog⁡c\log I_{\text{em}} = \log a + b\log c, so a log–log calibration plot is the standard diagnostic: a straight line of slope 1 confirms the linear (optically thin, non-self-absorbing) regime, while curvature or a slope below unity signals the onset of self-absorption.

Derivation ABoltzmann excitation fraction and its temperature sensitivity

Goal. Derive how sensitive the emission signal is to small flame-temperature fluctuations, starting from the Boltzmann population law, and express that sensitivity as a fractional-intensity-per-kelvin figure of merit.

Step 1 — write the excited fraction. From the Boltzmann relation, the fraction of atoms in the excited state is

f=N1N0+N1≈N1N0=g1g0e−ΔE/kTf = \frac{N_1}{N_0+N_1} \approx \frac{N_1}{N_0} = \frac{g_1}{g_0} e^{-\Delta E/kT}

valid whenever N1≪N0N_1 \ll N_0, which is true for every practical flame and every resonance line used in flame photometry (the excited fraction never exceeds a few percent even for the most easily excited elements).

Step 2 — differentiate with respect to temperature. Treat f(T)=(g1/g0)exp⁡(−ΔE/kT)f(T) = (g_1/g_0)\exp(-\Delta E/kT) as a smooth function of TT and take the logarithmic derivative, which is the natural way to express fractional sensitivity:

dln⁡fdT=ddT(−ΔEkT)=ΔEkT2\frac{d\ln f}{dT} = \frac{d}{dT}\left(-\frac{\Delta E}{kT}\right) = \frac{\Delta E}{kT^2}

This step uses only the chain rule on the exponent; the prefactor g1/g0g_1/g_0 is temperature-independent and drops out of the logarithmic derivative entirely.

Step 3 — convert to a fractional-change relation. For a small temperature perturbation δT\delta T,

δff≈ΔEkT2 δT\frac{\delta f}{f} \approx \frac{\Delta E}{kT^2}\,\delta T

Since emitted intensity IemI_{\text{em}} is directly proportional to ff (at fixed total atom density), the same fractional change applies to IemI_{\text{em}}:

δIemIem≈ΔEkT2 δT\frac{\delta I_{\text{em}}}{I_{\text{em}}} \approx \frac{\Delta E}{kT^2}\,\delta T

Step 4 — evaluate for sodium. Take ΔE=2.10\Delta E = 2.10 eV, T=2200T = 2200 K, k=8.617×10−5k = 8.617\times10^{-5} eV/K. Then

ΔEkT2=2.10(8.617×10−5)(2200)2=2.10417.0=5.04×10−3 K−1\frac{\Delta E}{kT^2} = \frac{2.10}{(8.617\times10^{-5})(2200)^2} = \frac{2.10}{417.0} = 5.04\times10^{-3}\ \text{K}^{-1}

So a temperature drift of δT=25\delta T = 25 K (about 1% of 2200 K, an entirely realistic burner fluctuation) produces

δIemIem≈(5.04×10−3 K−1)(25 K)=0.126=12.6%\frac{\delta I_{\text{em}}}{I_{\text{em}}} \approx (5.04\times10^{-3}\,\text{K}^{-1})(25\ \text{K}) = 0.126 = 12.6\%

a large, directly analytically consequential signal shift from a small, easily-occurring temperature change. This is the quantitative justification for internal standardization: the reference element, run through the identical flame at the identical instant, experiences essentially the same δT\delta T and its emission shifts by a comparable fractional amount, so the ratio of analyte to reference intensity is far more stable than either signal alone.

Limiting check. As T→∞T\to\infty, ΔE/kT2→0\Delta E/kT^2 \to 0: an infinitely hot flame would make emission intensity essentially insensitive to further small temperature changes (though of course such a flame does not exist for chemical burners, and ionization would dominate long before this limit is approached). This confirms the derivation's qualitative sense — sensitivity to TT decreases as flames get hotter, which is part of why higher-temperature sources (plasmas) are more stable emission platforms than flames.

Go deeper · Derivation B · 2 minSelf-absorption and the sublinear calibration curve
Derivation BSelf-absorption and the sublinear calibration curve

Goal. Show, from a simple two-region flame model, why emission intensity falls below direct proportionality to concentration at high analyte loading, and derive the condition distinguishing the linear from the self-absorbing regime.

Step 1 — model the flame as a hot emitting core inside a cooler absorbing mantle. Let the hot core (temperature ThT_h) contain excited-state density N1N_1 and emit intensity I0=hνA10N1VI_0 = h\nu A_{10} N_1 V isotropically, as derived above. Let the surrounding cooler mantle (temperature Tc<ThT_c < T_h) contain ground-state atoms at number density N0′N_0' over path length ℓ\ell toward the detector.

Step 2 — apply Beer's-law-type attenuation to the emitted photons as they traverse the mantle. The mantle absorbs the outgoing emission exactly as it would absorb an external lamp beam at the same wavelength (this is the physical content of "self-absorption" — the flame is simultaneously its own light source and its own absorbing medium):

Iobs=I0 e−kνN0′ℓI_{\text{obs}} = I_0\, e^{-k_\nu N_0' \ell}

where kνk_\nu is the absorption cross-section-like constant for that transition at line-center or integrated over the line profile.

Step 3 — relate both N1N_1 (emitting) and N0′N_0' (absorbing) to the same total analyte concentration cc. Both populations derive from the same aspirated sample, so to good approximation N1∝cN_1 \propto c (Step giving I0I_0) and N0′∝cN_0' \propto c as well, since essentially all atoms are ground-state regardless of position in the flame. Substituting N0′=κcN_0' = \kappa c for some proportionality constant κ\kappa:

Iobs=(a c) e−kνκℓ cI_{\text{obs}} = (a\,c)\, e^{-k_\nu \kappa \ell\, c}

with aa absorbing the emission-side proportionality constants from I0I_0.

Step 4 — examine the two limits. At low cc, expand the exponential: e−kνκℓc≈1−kνκℓce^{-k_\nu\kappa\ell c}\approx 1 - k_\nu\kappa\ell c, so

Iobs≈a c (1−kνκℓ c)=a c−a kνκℓ c2I_{\text{obs}} \approx a\,c\,(1 - k_\nu\kappa\ell\, c) = a\,c - a\,k_\nu\kappa\ell\,c^2

which is linear in cc to first order, with a negative quadratic correction — the calibration curve is nearly straight at low concentration but begins bending downward (concave, sublinear) as soon as cc becomes large enough that the quadratic term is non-negligible. This directly reproduces the empirical Lomakin–Scheibe exponent falling below 1 (b<1b<1) at higher concentrations, since a term of form c e−βcc\,e^{-\beta c} behaves like c bc^{\,b} with effective b<1b<1 over any finite window where βc\beta c is not small.

Step 5 — identify the self-reversal limit. At sufficiently high cc, the exponential decay in Step 2 can overtake the linear growth of I0I_0 entirely: differentiating Iobs(c)=ac e−βcI_{\text{obs}}(c) = a c\, e^{-\beta c} (with β≡kνκℓ\beta \equiv k_\nu\kappa\ell) and setting the derivative to zero,

dIobsdc=a e−βc(1−βc)=0  ⇒  cmax⁡=1β\frac{dI_{\text{obs}}}{dc} = a\,e^{-\beta c}(1-\beta c) = 0 \;\Rightarrow\; c_{\max} = \frac{1}{\beta}

Beyond cmax⁡=1/βc_{\max}=1/\beta, observed intensity actually decreases with increasing concentration — the calibration curve turns over. This is the origin of the warning, standard in quantitative flame photometry, never to work in the concentration regime where the calibration curve visibly bends: an operator unaware of self-reversal could badly misassign a high-concentration sample to a lower apparent concentration on the falling branch of the curve.

Limiting check. As β→0\beta \to 0 (an infinitely thin or infinitely hot, fully-ionized-mantle-free flame with no reabsorbing ground-state atoms outside the emitting zone), Iobs→acI_{\text{obs}} \to a c for all cc — perfectly linear emission with no self-absorption, consistent with the physical expectation that self-absorption requires a population of ground-state absorbers coexisting with the emitting population along the optical path.

Worked examples

Example 1 (foundational) — excited-state fraction and expected emission ratio for Na vs. K at a realistic flame temperature. A natural-gas/air flame runs at T=2000 K±15 KT = 2000\text{ K} \pm 15\text{ K} (a typical stated stability). Compare the excited-state populations of sodium (3p state, ΔENa=2.104\Delta E_{\text{Na}} = 2.104 eV, g1/g0=2g_1/g_0 = 2) and potassium (4p state, ΔEK=1.610\Delta E_{\text{K}} = 1.610 eV, g1/g0=2g_1/g_0 = 2) at this temperature, and estimate the fractional intensity uncertainty each contributes from the stated ±15\pm15 K temperature spread. Sodium. Using k=8.617×10−5k = 8.617\times10^{-5} eV/K,

ΔENakT=2.104(8.617×10−5)(2000)=2.1040.17234=12.21\frac{\Delta E_{\text{Na}}}{kT} = \frac{2.104}{(8.617\times10^{-5})(2000)} = \frac{2.104}{0.17234} = 12.21
N1N0∣Na=2 e−12.21=2 (4.98×10−6)=9.96×10−6\frac{N_1}{N_0}\bigg|_{\text{Na}} = 2\, e^{-12.21} = 2\,(4.98\times10^{-6}) = 9.96\times10^{-6}

Potassium.

ΔEKkT=1.6100.17234=9.343\frac{\Delta E_{\text{K}}}{kT} = \frac{1.610}{0.17234} = 9.343
N1N0∣K=2 e−9.343=2 (8.75×10−5)=1.75×10−4\frac{N_1}{N_0}\bigg|_{\text{K}} = 2\, e^{-9.343} = 2\,(8.75\times10^{-5}) = 1.75\times10^{-4}

The ratio of excited fractions is

(N1/N0)K(N1/N0)Na=1.75×10−49.96×10−6=17.6\frac{(N_1/N_0)_{\text{K}}}{(N_1/N_0)_{\text{Na}}} = \frac{1.75\times10^{-4}}{9.96\times10^{-6}} = 17.6

Potassium's excited fraction is nearly 18 times larger than sodium's at the same flame temperature and the same total atom density, purely because its resonance level sits lower — this is exactly why flame photometers historically found potassium easier to detect at low concentration than sodium, all else equal, and why potassium calibration curves reach usable sensitivity at lower solution concentrations. Temperature sensitivity. Using δI/I≈(ΔE/kT2) δT\delta I/I \approx (\Delta E/kT^2)\,\delta T from the earlier derivation, for sodium at T=2000T=2000 K:

ΔENakT2=2.104(8.617×10−5)(2000)2=2.104344.7=6.10×10−3 K−1\frac{\Delta E_{\text{Na}}}{kT^2} = \frac{2.104}{(8.617\times10^{-5})(2000)^2} = \frac{2.104}{344.7} = 6.10\times10^{-3}\ \text{K}^{-1}
δII∣Na=(6.10×10−3)(15)=0.0915≈9.2%\frac{\delta I}{I}\bigg|_{\text{Na}} = (6.10\times10^{-3})(15) = 0.0915 \approx 9.2\%

For potassium, whose ΔE\Delta E is smaller, the sensitivity is proportionally lower:

ΔEKkT2=1.610344.7=4.67×10−3 K−1,δII∣K=(4.67×10−3)(15)=0.070≈7.0%\frac{\Delta E_{\text{K}}}{kT^2} = \frac{1.610}{344.7} = 4.67\times10^{-3}\ \text{K}^{-1}, \qquad \frac{\delta I}{I}\bigg|_{\text{K}} = (4.67\times10^{-3})(15) = 0.070 \approx 7.0\%

Both elements show emission fluctuations of several percent from a modest 15 K burner drift — confirming why single-point, non-internal-standard flame photometry is inherently limited to roughly 2–5% relative precision at best, and why the internal-standard method (Example 2) is standard clinical practice.

Example 2 (applied) — serum sodium determination using a lithium internal standard. A clinical flame photometer measures serum sodium with lithium chloride added as an internal standard to every sample and standard at a fixed concentration of 15.0 mg/L Li⁺. Three sodium standards and a serum sample (diluted 1:200 with the same Li-spiked diluent) give the following intensity ratios R=INa/ILiR = I_{\text{Na}}/I_{\text{Li}}:

Standard[Na⁺], mg/LINaI_{\text{Na}} (arb.)ILiI_{\text{Li}} (arb.)R=INa/ILiR = I_{\text{Na}}/I_{\text{Li}}
S110.0842100.400
S220.01682120.792
S330.02492081.197
Serum (diluted)?2142091.024

Step 1 — confirm linearity and fit. A least-squares fit of RR vs. concentration cc (mg/L) through the three standards: Mean cˉ=20.0\bar c = 20.0, mean Rˉ=0.7963\bar R = 0.7963. Compute slope via

m=∑(ci−cˉ)(Ri−Rˉ)∑(ci−cˉ)2m = \frac{\sum (c_i-\bar c)(R_i-\bar R)}{\sum (c_i - \bar c)^2}

Deviations: cc: −10,0,+10-10, 0, +10; RR: −0.3963,−0.0043,+0.4007-0.3963, -0.0043, +0.4007.

∑(ci−cˉ)(Ri−Rˉ)=(−10)(−0.3963)+(0)(−0.0043)+(10)(0.4007)=3.963+4.007=7.970\sum (c_i-\bar c)(R_i-\bar R) = (-10)(-0.3963)+(0)(-0.0043)+(10)(0.4007) = 3.963+4.007 = 7.970
∑(ci−cˉ)2=100+0+100=200\sum (c_i-\bar c)^2 = 100+0+100 = 200
m=7.970200=0.03985 (mg/L)−1m = \frac{7.970}{200} = 0.03985\ \text{(mg/L)}^{-1}

Intercept: R0=Rˉ−mcˉ=0.7963−(0.03985)(20.0)=0.7963−0.7970=−0.0007≈0.000R_0 = \bar R - m\bar c = 0.7963 - (0.03985)(20.0) = 0.7963 - 0.7970 = -0.0007 \approx 0.000 So the fitted calibration line is

R=0.0399 c−0.0007≈0.0399 cR = 0.0399\,c - 0.0007 \approx 0.0399\,c

essentially passing through the origin, as expected for a well-behaved internal-standard ratio method in the linear regime — the near-zero intercept is itself a useful diagnostic that no significant background/blank emission or systematic bias is present. Step 2 — back-calculate serum sodium (diluted). For the serum sample, Rserum=1.024R_{\text{serum}} = 1.024:

cdiluted=Rserum+0.00070.0399=1.02470.0399=25.68 mg/Lc_{\text{diluted}} = \frac{R_{\text{serum}} + 0.0007}{0.0399} = \frac{1.0247}{0.0399} = 25.68\ \text{mg/L}

Step 3 — correct for the 1:200 dilution.

cserum=25.68 mg/L×200=5136 mg/Lc_{\text{serum}} = 25.68\ \text{mg/L} \times 200 = 5136\ \text{mg/L}

Step 4 — convert to clinical units (mmol/L), using MNa=22.99M_{\text{Na}} = 22.99 g/mol:

cserum=5136 mg/L22.99 g/mol=5.136 g/L22.99 g/mol=0.2234 mol/L=223 mmol/Lc_{\text{serum}} = \frac{5136\ \text{mg/L}}{22.99\ \text{g/mol}} = \frac{5.136\ \text{g/L}}{22.99\ \text{g/mol}} = 0.2234\ \text{mol/L} = 223\ \text{mmol/L}

This value (223 mmol/L) is far above the true physiological range (135–145 mmol/L), signaling that the assumed dilution or intensity values in this illustrative dataset are for teaching the algebra rather than representing an actual patient; a real assay would show cdilutedc_{\text{diluted}} around 0.7–0.75 mg/L for a 1:200 dilution of a normal 140 mmol/L serum. The chemistry and arithmetic pathway — fit RR vs. cc, solve for cdilutedc_{\text{diluted}}, multiply by the dilution factor, convert mass to molar concentration — is exactly what a clinical analyzer executes on every sample, and is the graded procedure to reproduce regardless of the specific numbers chosen.

Example 3 (exam/competition level) — resolving self-absorption from two-point dilution data. A student runs a "concentrated" sodium standard at 500 mg/L and measures emission intensity I500=812I_{500} = 812 (arbitrary units, instrument linear in photon flux). Suspecting self-absorption, the student dilutes this standard 1:2 (to 250 mg/L) and remeasures, obtaining I250=460I_{250} = 460. A third dilution to 125 mg/L gives I125=241I_{125} = 241. Using the self-absorption model Iobs=a c e−βcI_{\text{obs}} = a\,c\,e^{-\beta c} derived above, estimate β\beta and aa, and predict whether a fourth point at 1000 mg/L would lie on the rising or falling branch of the curve. Step 1 — take the ratio of two data points to eliminate aa. From

I250I125=250 e−250β125 e−125β=2 e−125β\frac{I_{250}}{I_{125}} = \frac{250\, e^{-250\beta}}{125\, e^{-125\beta}} = 2\, e^{-125\beta}

Insert the numbers:

460241=1.909=2 e−125β\frac{460}{241} = 1.909 = 2\,e^{-125\beta}
e−125β=0.9547  ⇒  −125β=ln⁡(0.9547)=−0.04638e^{-125\beta} = 0.9547 \;\Rightarrow\; -125\beta = \ln(0.9547) = -0.04638
β=0.04638125=3.71×10−4 (mg/L)−1\beta = \frac{0.04638}{125} = 3.71\times10^{-4}\ (\text{mg/L})^{-1}

Step 2 — check consistency with the third pair (500 vs. 250).

I500I250=500 e−500β250 e−250β=2 e−250β\frac{I_{500}}{I_{250}} = \frac{500\,e^{-500\beta}}{250\,e^{-250\beta}} = 2\,e^{-250\beta}

Predicted: e−250β=e−250(3.71×10−4)=e−0.0928=0.9114e^{-250\beta} = e^{-250(3.71\times10^{-4})} = e^{-0.0928} = 0.9114, so predicted ratio =2(0.9114)=1.823= 2(0.9114) = 1.823. Observed ratio: 812/460=1.765812/460 = 1.765 — close (within ~3%) to the model prediction of 1.823, supporting the exponential self-absorption model over this range; the small residual reflects real instrumental noise and the approximation that β\beta is concentration-independent, which holds only approximately. Step 3 — solve for aa using the 125 mg/L point.

a=I125c e−βc=241125 e−(3.71×10−4)(125)=241125 e−0.04638=241125(0.9547)=241119.3=2.020 (arb. units per mg/L)a = \frac{I_{125}}{c\, e^{-\beta c}} = \frac{241}{125\, e^{-(3.71\times10^{-4})(125)}} = \frac{241}{125\, e^{-0.04638}} = \frac{241}{125(0.9547)} = \frac{241}{119.3} = 2.020\ \text{(arb. units per mg/L)}

Step 4 — find the turnover concentration cmax⁡=1/βc_{\max} = 1/\beta.

cmax⁡=13.71×10−4=2695 mg/Lc_{\max} = \frac{1}{3.71\times10^{-4}} = 2695\ \text{mg/L}

Since the proposed fourth point, 1000 mg/L, is well below cmax⁡=2695c_{\max} = 2695 mg/L, the model predicts it still lies on the rising branch — intensity should still be increasing with concentration there, though with visibly reduced slope compared to the low-concentration (nearly linear) region. Predicted intensity:

I1000=a (1000) e−β(1000)=2.020×1000×e−0.371=2020×0.690=1394I_{1000} = a\,(1000)\,e^{-\beta(1000)} = 2.020\times1000\times e^{-0.371} = 2020\times0.690 = 1394

compared to the naive linear extrapolation from the 125 mg/L point, a c=2.020×1000=2020a\,c = 2.020\times1000 = 2020 — a 31% suppression from self-absorption at 1000 mg/L, illustrating exactly why quantitative flame-emission work is normally confined to concentrations well below where this kind of curvature becomes significant, and why a log–log calibration check (slope should stay at 1.00) is standard practice before trusting any single-point emission reading.

Limiting cases, scope & common traps

Regime / caseBehaviorConsequence for the analyst
Low concentration, optically thin flameIem∝cI_{\text{em}} \propto c exactly (Lomakin–Scheibe b≈1b\approx1)Standard linear calibration valid; preferred working range
High concentrationSelf-absorption bends curve down; possible self-reversal beyond cmax⁡=1/βc_{\max}=1/\betaDilute samples into the linear range; never trust readings near visible curvature
Easily ionized element (K, Cs, Rb) at low concentrationSignificant fraction converts to M⁺, reducing neutral-atom (and hence emitting) populationAdd large excess of an even more easily ionized element (e.g., CsCl) as ionization suppressor
High total dissolved solids / viscosity change between sample and standardNebulization efficiency differs, changing atom delivery rate independent of true concentrationMatch matrix, or use internal-standard ratio method to cancel the effect
T→T \to very high (plasma-like)ΔE/kT2→0\Delta E/kT^2 \to 0; excited fraction less sensitive to small TT drift, but many more elements become excitable, including refractory onesForeshadows why Inductively Coupled Plasma Spectrometry achieves broader multi-element coverage and better precision than flame emission
Two elements with very different ΔE\Delta E (e.g., Na vs. Ca resonance line at 422.7 nm, ΔE≈2.93\Delta E \approx 2.93 eV)Ca requires hotter flame (acetylene–nitrous oxide) for adequate excited population; a cool air–propane flame under-excites Ca badlyMatch flame type/temperature to the analyte's excitation energy, not a single "universal" flame
Setting b=1b=1 identically at all ccIdealized linear limiting case (Step 4 of self-absorption derivation, βc≪1\beta c \ll 1)Real curves always deviate at sufficiently high cc; treat linearity as a low-concentration approximation, not a law

The chemist behind it

Flame emission's modern quantitative form owes most to Robert Bunsen and Gustav Kirchhoff, working together at Heidelberg in the early 1860s. Bunsen had, a few years earlier, perfected the smokeless, nearly colorless premixed gas burner that now bears his name — designed originally simply to give a hot, clean flame for general laboratory heating, but which turned out to be the ideal, spectrally quiet background against which faint atomic emission lines from trace metal salts could finally be seen and measured. Kirchhoff, a physicist steeped in the study of light and radiation, brought the spectroscope: a prism-and-telescope arrangement that could resolve the colored flame light into discrete lines rather than a diffuse glow. Individually, Bunsen had the clean flame and Kirchhoff had the means to disperse the light it produced; together, in 1859–1861, they systematically catalogued the characteristic emission-line wavelengths of the alkali and alkaline-earth metals and demonstrated that each element gives an invariant spectral fingerprint, independent of the compound from which it came. This was simultaneously the founding act of quantitative flame photometry and, through the same reasoning applied to sunlight, the birth of astronomical spectroscopy — Kirchhoff and Bunsen used their new instrument to identify sodium in the solar spectrum, and within a few years the pair discovered two new elements, cesium (1860, from the Latin caesius, sky-blue, for its blue emission lines) and rubidium (1861, for its deep red lines), purely from previously unrecorded spectral lines appearing in mineral-water residues that no known element could explain.

Origins · 1 minGustav Robert Kirchhoff and how this was discovered
GR
Gustav Robert Kirchhoff
1824–1887 · University of Heidelberg, Germany, 1850s–1870s

Kirchhoff is best known to circuit students for his laws of current and voltage conservation, but his deepest scientific legacy is in radiation and spectroscopy. Working alongside Bunsen, he established that every chemical element emits and absorbs light at its own fixed set of wavelengths — the principle now called Kirchhoff's law of thermal radiation in its spectroscopic form, that a substance absorbs at exactly the wavelengths it can emit when heated. This single insight underlies the entire logical structure connecting flame emission to flame/AAS absorption: they are two views of the identical set of atomic energy levels, probed by two different experimental strategies.

History
- 1854–1855 — Bunsen develops the premixed, non-luminous laboratory gas burner at Heidelberg.
- 1859 — Kirchhoff formulates his law relating emissive and absorptive power, providing the theoretical basis for atomic spectral analysis.
- 1860 — Bunsen and Kirchhoff discover cesium via its unique blue emission lines in mineral spring water residue.
- 1861 — the pair discovers rubidium from its characteristic red emission lines, using the same flame-spectroscope method.
- 1930s–1950s — flame photometry is industrialized and instrumented for routine clinical and agricultural analysis, driven by demand for fast serum electrolyte and soil-nutrient testing.
Why are alkali metals so much easier to analyze by flame emission than by flame AAS on a routine benchtop flame?

Alkali metals have unusually low-lying resonance excited states (1.6–2.1 eV for K, Na, Li, Rb, Cs), so even a modest chemical flame at 1700–2400 K populates their excited states enough to give bright, easily measured emission, per the Boltzmann relation N1/N0=(g1/g0)e−ΔE/kTN_1/N_0 = (g_1/g_0)e^{-\Delta E/kT}. Elements with higher-energy resonance transitions (most transition and heavy metals) have excited fractions many orders of magnitude smaller in the same flame, giving emission too weak to measure reliably — those elements are usually better handled by flame AAS, which depends on the much larger ground-state population instead, or by a hotter plasma source.

What causes ionization interference, and how does adding excess cesium chloride fix it?

In the hot flame, some fraction of neutral atoms lose an electron to become singly charged ions, M⇌M++e−M \rightleftharpoons M^+ + e^-, especially for easily ionized elements like potassium, rubidium, and cesium. Ionized atoms have a completely different electronic structure and do not emit at the neutral atom's characteristic wavelength, so ionization directly reduces measured emission (and equally reduces AAS absorbance) below what total-analyte concentration alone would predict. Adding a large excess of an even more easily ionized element (commonly CsCl, since cesium has the lowest ionization energy of the stable alkalis) floods the flame with free electrons from cesium's own ionization, which by Le Chatelier's principle suppresses further ionization of the analyte and pushes its equilibrium back toward the neutral, emitting atom.

Is flame photometry quantitative on its own, or does it always need an internal standard?

It can be run as a simple external-calibration method, and for undemanding applications (rough soil-potassium screening, teaching labs) that is common practice. But because emission intensity depends exponentially on flame temperature while absorption depends only weakly on it, single-external-standard flame emission is inherently more susceptible to day-to-day and even minute-to-minute drift than flame AAS. Clinical and other high-precision applications virtually always use an internal standard (commonly lithium or cesium at a fixed spiked concentration) precisely to cancel that thermal and nebulization-rate noise via an intensity ratio rather than a raw intensity reading.

How does self-absorption differ from simple signal attenuation by an optically dense sample?

Self-absorption is a special case in which the same chemical species is simultaneously the emitter (in the hot flame core) and the absorber (in the cooler outer mantle) of the identical wavelength of light — the analyte is reabsorbing its own emitted photons. This differs from generic attenuation (for example, light scattering by unvaporized particulate or absorption by an unrelated matrix species) in that self-absorption specifically increases with analyte concentration and specifically affects the analyte's own resonance line, producing the characteristic sublinear-then-reversing calibration curve derived from the two-region flame model.

Why does calcium require a hotter flame (acetylene–nitrous oxide) than sodium or potassium for good flame-emission sensitivity?

Calcium's principal emission line (422.7 nm) arises from a transition with ΔE≈2.93\Delta E \approx 2.93 eV, notably larger than sodium's 2.10 eV or potassium's 1.61 eV. Because the excited fraction depends exponentially on −ΔE/kT-\Delta E/kT, calcium's excited population in a cool air–propane flame (T ≈ 1900–2000 K) is far smaller, relatively, than sodium's or potassium's would be at the same temperature. A hotter acetylene–nitrous oxide flame (T ≈ 2600–2900 K) raises kTkT enough to populate calcium's excited state adequately, while also helping dissociate the refractory calcium-oxide-like species that can otherwise trap calcium as a stable molecular compound rather than free atoms.

Can flame photometry measure transition metals or trace heavy metals like lead or cadmium?

Not well, and it is essentially never used for that purpose. Transition and heavy metals typically have resonance transitions requiring several electron-volts more energy than the alkali metals, so their excited-state fractions in an ordinary chemical flame are vanishingly small — often many orders of magnitude below what a photomultiplier can detect above background. Such elements are handled instead by flame or graphite-furnace atomic absorption (which relies on the much larger ground-state population) or by inductively coupled plasma techniques, whose much higher temperatures (6000–10000 K) generate adequate excited-state populations for a vastly wider range of elements.

Why is a log–log plot of intensity vs. concentration the standard diagnostic tool in flame emission work, rather than a linear plot?

The Lomakin–Scheibe relation Iem=a c bI_{\text{em}} = a\,c^{\,b} becomes, after taking logarithms, the linear form log⁡Iem=log⁡a+blog⁡c\log I_{\text{em}} = \log a + b\log c. Plotting on log–log axes therefore converts a power-law relationship into a straight line whose slope directly reads off the exponent bb. A slope of exactly 1.00 confirms the optically thin, non-self-absorbing regime is in force across the tested concentration range; any systematic departure from a slope of 1.00 (curvature or a slope significantly below 1) is an immediate, quantitative flag that self-absorption or another concentration-dependent effect is distorting the calibration, well before that distortion might otherwise be missed on a linear scale.

Before this

Where this leads