Detectors for Atomic Spectroscopy

A photomultiplier tube converts a single photon into a burst of roughly a million electrons. This happens in under ten nanoseconds. That single fact — enormous internal gain achieved before any external noise can be added — explains a lot. It explains why PMTs dominated atomic spectroscopy for sixty years. It also explains why replacing them with solid-state arrays took a genuine reinvention of detector electronics. It was not just a change of parts. The detector is the quiet partner in every atomic spectroscopy measurement. The plasma or flame does the atomizing and exciting. The monochromator or echelle does the wavelength sorting. But the detector decides how many wavelengths you can watch at once. It decides how low a concentration you can see above the noise floor. And it decides how fast you can finish a run of forty elements on a hundred samples. This matters because the choice between a scanning photomultiplier and a multichannel array is important. It is really a choice between two philosophies of measurement. Sequential detection with a single PMT measures one wavelength at a time. It moves a grating or a slit to hop from line to line. This approach is simple, cheap, and photometrically clean. But it is slow when the sample list and the element list both get long. Simultaneous detection with a charge-transfer array (CCD or CID) captures many wavelengths at once. Sometimes it captures the whole 165–800 nm range in one integration. This comes at the cost of a more complex readout chain. It also brings a differently shaped noise budget. Modern inductively coupled plasma optical emission spectrometers (ICP-OES) embody both philosophies. So do time-of-flight or scanning ICP mass spectrometers. A working analytical chemist needs to know which one a given problem demands. The physics underneath is a story of two very different approaches. Both turn photons into a countable, storable signal. A PMT is a vacuum-tube current amplifier. It consists of a photocathode, dynode chain, and anode, all riding on a high-voltage divider. It produces an analog current or a train of discrete pulses. External electronics then integrate or count this output. A CCD or CID is a semiconductor charge-integration device. Photons generate electron–hole pairs in a silicon lattice. The charge simply accumulates in a potential well — a tiny capacitor. It is read out once, at the end of the exposure. In the CID case, it can be read non-destructively many times during the exposure. One device amplifies as it goes; the other stores and waits. That distinction drives everything downstream. It shapes dynamic range, dark current, linearity, cost. It also shapes the entire architecture of a "simultaneous" versus "sequential" instrument. This article works through several key topics. It covers the physical structure of the photomultiplier and of charge-transfer arrays. It covers the governing equations for gain, quantum efficiency, dynode statistics, and shot-noise-limited signal-to-noise ratio. It includes three fully worked numerical examples. These span a textbook PMT gain calculation, a real ICP-OES detection-limit comparison, and a competition-level noise-budget problem. Finally, it surveys where each detector type is deployed today. This ranges from trace-metal environmental labs to semiconductor-fab process control.

Key takeaways
  • A PMT multiplies a single photoelectron through N dynodes, each with gain δ, giving overall gain G=δNG = \delta^N, typically 10510^5–10710^7.
  • CCDs and CIDs integrate photogenerated charge in silicon potential wells; the CID's non-destructive readout of each pixel is what enables true random-access, simultaneous multi-element operation.
  • Sequential (scanning PMT) instruments trade speed for photometric simplicity and superb dynamic range; simultaneous (array) instruments trade some per-line optimization for enormous throughput.
  • Dark current, read noise, and full-well capacity — not just quantum efficiency — set the practical detection limit of an array detector.
  • Photon shot noise is the fundamental floor for both detector families once amplifier and readout noise are minimized; SNR=np\mathrm{SNR} = \sqrt{n_p} for a Poisson-limited photon count npn_p.
  • Blooming, charge-transfer efficiency, and anti-blooming drains are CCD-specific artifacts that a CID's pixel architecture largely avoids.
  • Key idea: the detector's architecture — amplify-then-read (PMT) versus integrate-then-read (array) — dictates whether an instrument scans one line at a time or captures a whole spectrum simultaneously, and that single design choice cascades into every figure of merit that follows.

Where this lives today

Modern ICP-OES instruments split cleanly along this article's central axis. High-throughput environmental and industrial labs running EPA Method 200.7/6010 panels on hundreds of water or soil digests per day almost universally use simultaneous CID- or CCD-based (often now CMOS-hybrid) spectrometers capable of measuring 20–70 elements from a single 15–60 s plasma exposure, because the economics of sample throughput dominate. Ultratrace geochemistry and semiconductor-grade reagent qualification, where a handful of elements must be measured at parts-per-trillion levels, still sometimes favor dedicated sequential scanning instruments (or high-resolution sector-field ICP-MS, whose detector electronics — electron multipliers and Faraday cups — echo the same amplify-then-read versus integrate-then-read dichotomy) for their superior per-line dynamic range and lower stray-light contamination. Photomultiplier tubes persist essentially unchanged in flame atomic absorption spectroscopy (FAAS) and in graphite furnace AAS, in scintillation counting for radiochemistry, and in the reference channels of many double-beam UV-Vis spectrophotometers, precisely because a single, well-characterized wavelength with excellent dynamic range is exactly what those techniques need. CCD arrays dominate Raman spectroscopy, laser-induced breakdown spectroscopy (LIBS), and echelle-based simultaneous ICP-OES, where imaging an entire two-dimensional order-sorted spectrum onto one chip is the whole point of the optical design. CIDs, though a smaller commercial niche today, remain valued in process-monitoring and metallurgical ICP-OES applications where a few very bright matrix lines (Fe, major alloy constituents) must be measured alongside trace dopants in the same transient plasma burst — the non-destructive, variable-integration readout is tailor-made for exactly that combination. In biological and clinical mass spectrometry, related charge-integration and electron-multiplier hybrids (e.g., discrete-dynode electron multipliers in ICP-MS, and time-to-digital converters in time-of-flight instruments) inherit the same Poisson-statistics and gain-variance considerations developed here, linking detector physics from optical atomic spectroscopy directly to trace elemental mass spectrometry used in Inductively Coupled Plasma Spectrometry and beyond.

Ask the tutor: How would the noise budget and detection-limit formulas change for a hybrid CMOS active-pixel sensor used in newer ICP-OES instruments, compared with the classic CCD and CID architectures covered here?

The physical picture

A photomultiplier tube begins with a photocathode: a thin semiconductor film (bialkali Sb-Rb-Cs, or Ag-O-Cs for red-extended response) deposited on the inside of an evacuated glass envelope. An incoming photon of energy hνh\nu ejects a photoelectron via the photoelectric effect, provided hνh\nu exceeds the material's work function. That single electron is then accelerated by a potential difference of roughly 100 V toward the first dynode, a curved electrode coated with a secondary-emissive material (BeO, GaP, or Cs3_3Sb). On impact, the primary electron liberates several secondary electrons — the dynode gain δ\delta, typically 3–6 per stage. This process repeats through 9–12 dynodes, each held at successively higher potential by a resistive voltage-divider network, so that the single photoelectron becomes an avalanche of 10510^5–10710^7 electrons by the time it reaches the anode, all within about 5–20 ns. The anode current pulse is either integrated as an analog photocurrent (the traditional mode) or, in dim-light applications, counted as discrete pulses in photon-counting mode, where each detected photon produces a resolvable pulse above a discriminator threshold. A charge-coupled device (CCD) instead uses a two-dimensional array of metal-oxide-semiconductor (MOS) capacitors fabricated on a doped silicon substrate. A photon absorbed in the depletion region promotes an electron from the valence to the conduction band, creating an electron–hole pair; the electron is trapped in the local potential well beneath a gate electrode, and holes are swept away. Charge accumulates linearly with exposure time and photon flux until the well nears its full-well capacity (typically 10510^5–10610^6 electrons per pixel for scientific-grade CCDs). Readout is destructive and sequential: at the end of integration, charge packets are shifted row by row into a serial register and then clocked out pixel by pixel to an on-chip amplifier and analog-to-digital converter (ADC). Because readout empties the whole array, a CCD cannot be queried mid-exposure without losing the accumulated signal — a serious limitation for the very different dynamic ranges of major and trace emission lines in atomic spectra. A charge-injection device (CID) solves exactly that problem. Each CID pixel has two capacitively coupled electrodes so the stored charge can be sensed non-destructively — by measuring the change in voltage induced when charge is shifted between the two electrodes within the pixel — without removing it from the well. This allows random-access readout: any individual pixel, or any subset of pixels, can be interrogated at any time during the integration, as many times as desired, and the accumulated charge is only actually cleared ("injected" into the substrate) when the experimenter chooses. In atomic emission work this means a bright ionic line and a faint trace-element line on the same array can each be read out with an integration time individually matched to their intensity — a capability with no CCD equivalent, since a CCD's single global readout forces one compromise integration time for the whole chip (unless multiple non-destructive reads are averaged, a technique CCDs also use but less flexibly).

Figure
Photomultiplier tube: photon in, amplified electron pulse out
hν photocathode e⁻ (photoelectron) D1 D2 D3 D4 D5 D6 anode I = G·i₀ to ADC / discriminator δ ≈ 4 per dynode, N = 9–12 stages Gain G = δ^N, typically 10⁵–10⁷; rise time ~5–20 ns
Figure
CCD vs CID pixel architecture and readout timing
CCD: destructive, sequential readout bright well shift → serial register → single amp → ADC whole array cleared each readout — cannot query mid-exposure CID: non-destructive, random-access bright, read often faint, read once each pixel: two coupled electrodes, sense ΔV without emptying well bright line read short/often; faint line integrated long — same frame Simultaneous ICP-OES: echelle-dispersed spectrum imaged onto 2-D array, each (order, position) pixel maps to one wavelength — hundreds of lines per exposure.

The core ideas

Photomultiplier gain. For a dynode chain of NN stages, each with secondary-emission ratio δ\delta (the average number of secondary electrons ejected per incident electron), the overall current gain is

G=δNG = \delta^{N}

where:

  • GG = overall electron multiplication factor (dimensionless), the ratio of anode current pulse charge to photocathode photoelectron charge
  • δ\delta = secondary emission coefficient per dynode stage, itself a function of interstage voltage VdV_d roughly as δ∝Vd k\delta \propto V_d^{\,k} with k≈0.7k \approx 0.7–0.80.8
  • NN = number of dynode stages (commonly 9–12)

Reading: because δ\delta itself grows with the interstage voltage, the total PMT gain is an extremely steep function of the applied high voltage — doubling VdV_d can increase GG by an order of magnitude or more, which is why PMT power supplies must be stabilized to millivolt-level precision.

Anode signal and shot-noise-limited SNR. The mean anode current from a steady photon flux is

ia=G e Φ ηi_a = G\, e\, \Phi\, \eta

where:

  • iai_a = mean anode current (A)
  • ee = electron charge, 1.602×10−191.602\times10^{-19} C
  • Φ\Phi = photon arrival rate at the photocathode (photons s−1^{-1})
  • η\eta = photocathode quantum efficiency (photoelectrons per incident photon, dimensionless, typically 0.1–0.3 in the visible)

The number of photoelectrons collected in an integration time tt is npe=ηΦtn_{pe} = \eta \Phi t, a Poisson-distributed random variable, so the photon-shot-noise-limited signal-to-noise ratio is

SNR=npenpe=npe=ηΦt\mathrm{SNR} = \frac{n_{pe}}{\sqrt{n_{pe}}} = \sqrt{n_{pe}} = \sqrt{\eta \Phi t}

where:

  • SNR\mathrm{SNR} = ratio of mean signal to its standard deviation, assuming Poisson (shot-noise) statistics dominate over dark current and readout noise
  • the square-root dependence means doubling SNR requires quadrupling either the flux or the integration time

Reading: this equation is the ultimate floor for any photon-counting detector, PMT or array; everything else — dark counts, dynode statistics, read noise — only adds noise on top of this irreducible Poisson term, never removes it.

Array detector signal-to-noise ratio. For a CCD or CID pixel collecting SS photoelectrons of signal, with dark current contributing DD electrons and read noise RR electrons (rms) per readout, the total noise combines in quadrature (each independent noise source is uncorrelated, so variances add):

SNR=SS+D+R2\mathrm{SNR} = \frac{S}{\sqrt{S + D + R^{2}}}

where:

  • SS = signal electrons accumulated during integration (Poisson variance =S=S)
  • DD = dark-current electrons accumulated over the same integration time (Poisson variance =D=D)
  • RR = read noise per readout event, root-mean-square electrons (Gaussian, added as R2R^2 since it is already a standard deviation, not a Poisson mean)

Reading: at high signal levels the SS term dominates and SNR→S\mathrm{SNR}\to\sqrt{S} (photon-shot-noise limited, matching the PMT expression above); at very low signal the R2R^2 term dominates, and read noise — not photon statistics — sets the practical detection limit, which is exactly why cooled, low-read-noise CCDs matter for trace analysis.

Derivation ADynode-chain gain and its variance — why PMT pulse-height distributions are broad

Goal: show that the overall PMT gain is G=δNG=\delta^N and that the relative variance of the output pulse charge is dominated by fluctuations at the first dynode, not the later ones — the key fact behind PMT pulse-height spectra and the choice of high first-dynode voltage.

Step 1 — single-stage multiplication. A single photoelectron striking dynode 1 produces a random number k1k_1 of secondary electrons, with mean ⟨k1⟩=δ\langle k_1\rangle = \delta and variance σ12\sigma_1^2 (secondary emission is well modeled as a Poisson-like process, so to first approximation σ12≈δ\sigma_1^2 \approx \delta, though real dynodes show slightly super-Poissonian statistics; a common empirical factor is σ12=b δ\sigma_1^2 = b\,\delta with b≈1b\approx 1–1.51.5). Each of those k1k_1 electrons then independently produces a further generation at dynode 2, and so on.

Step 2 — mean output after N stages. By induction, the mean number of electrons after stage nn is ⟨kn⟩=δ n\langle k_n\rangle = \delta^{\,n}, because each stage multiplies the expected population by δ\delta regardless of the details of the distribution:

⟨kN⟩=δ⋅δ⋯δ  (N times)=δN=G\langle k_N \rangle = \delta \cdot \delta \cdots \delta \;(N \text{ times}) = \delta^{N} = G

This confirms the gain formula stated above purely from the multiplicative structure of cascaded amplification.

Step 3 — propagation of variance through the cascade. For a branching (cascade multiplication) process, the relative variance of the final output, σN2/⟨kN⟩2\sigma_N^2/\langle k_N\rangle^2, obeys the recursive relation (a standard result for cascaded random multiplication, sometimes called the Burgess variance theorem):

σN2⟨kN⟩2=σ12δ2+1δ(σ22δ2+1δ(⋯))\frac{\sigma_N^2}{\langle k_N\rangle^2} = \frac{\sigma_1^2}{\delta^2} + \frac{1}{\delta}\left(\frac{\sigma_2^2}{\delta^2} + \frac{1}{\delta}\Big(\cdots\Big)\right)

Expanding this geometric-like series with σn2≈δ\sigma_n^2 \approx \delta (Poisson single-stage statistics) for every stage gives

σN2⟨kN⟩2≈1δ(1+1δ+1δ2+⋯ )=1δ⋅11−1/δ=1δ−1\frac{\sigma_N^2}{\langle k_N\rangle^2} \approx \frac{1}{\delta}\left(1 + \frac{1}{\delta} + \frac{1}{\delta^2} + \cdots\right) = \frac{1}{\delta}\cdot\frac{1}{1-1/\delta} = \frac{1}{\delta - 1}

Step 4 — interpretation. The relative variance of the fully amplified pulse converges to 1/(δ−1)1/(\delta-1), a constant set almost entirely by the first dynode's gain δ\delta, because each subsequent stage's fractional contribution is suppressed by an extra factor of 1/δ1/\delta. For a typical δ≈4\delta \approx 4–5, the relative standard deviation of the single-photoelectron pulse-height distribution is 1/(δ−1)≈0.5\sqrt{1/(\delta-1)} \approx 0.5–0.6 — a broad, roughly 50% spread — which is exactly the characteristic single-photoelectron peak shape seen in PMT pulse-height spectra, and the reason instrument designers run the first dynode at unusually high voltage (to push δ1\delta_1 as high as possible and narrow this dominant term).

Go deeper · Derivation B · 1 minDetection limit from the Poisson-plus-read-noise budget
Derivation BDetection limit from the Poisson-plus-read-noise budget

Goal: derive an expression for the minimum detectable analyte concentration cLc_L in terms of blank signal, blank noise, and calibration sensitivity, then show how the array SNR expression feeds directly into it — connecting the microscopic noise sources to the macroscopic figure of merit reported on every spectrometer datasheet.

Step 1 — the IUPAC/analytical convention. The limit of detection is conventionally defined as the analyte concentration giving a net signal equal to kk (commonly k=3k=3) times the standard deviation of the blank signal, sbls_{bl}:

SL−Sˉbl=k sblS_L - \bar{S}_{bl} = k\, s_{bl}

where SLS_L is the signal at the detection limit and Sˉbl\bar S_{bl} the mean blank signal.

Step 2 — introduce calibration sensitivity. If the calibration curve is linear over the relevant range with slope mm (sensitivity, signal units per unit concentration), then

SL−Sˉbl=m cLS_L - \bar{S}_{bl} = m\,c_L

Combining with Step 1,

cL=k sblmc_L = \frac{k\, s_{bl}}{m}

Step 3 — express sbls_{bl} using the detector noise budget. The blank signal itself is a small photoelectron count SblS_{bl} sitting on dark current DD and read noise RR, so by the array SNR expression derived above, the noise (not signal) on the blank is

sbl=Sbl+D+R2s_{bl} = \sqrt{S_{bl} + D + R^{2}}

For a PMT operated in analog mode, replace DD with the dark-current-equivalent photoelectron count and RR with amplifier/electrometer noise; the same algebra applies.

Step 4 — assemble and interpret. Substituting,

cL=kSbl+D+R2mc_L = \frac{k\sqrt{S_{bl} + D + R^{2}}}{m}

Two limiting regimes fall directly out of this expression. When R2≫Sbl+DR^2 \gg S_{bl}+D (read-noise-limited, typical of a warm, fast-readout array), cL∝R/mc_L \propto R/m: the detection limit is fixed by electronics and can only be improved by better electronics or higher sensitivity mm, not by integrating longer. When Sbl+D≫R2S_{bl}+D \gg R^2 (shot-noise or dark-current limited, typical of cooled CCDs or PMTs at long integration), cL∝Sbl+D /mc_L \propto \sqrt{S_{bl}+D}\,/m, and because SblS_{bl} and DD both grow linearly with integration time tt while mm (proportional to accumulated analyte signal) also grows linearly with tt, cL∝1/tc_L \propto 1/\sqrt{t} — the familiar rule that detection limits improve with the square root of integration time, but only once you are past the read-noise floor.

Worked examples

Example 1 — Foundational: PMT anode current and photon-counting rate. A PMT with photocathode quantum efficiency η=0.22\eta = 0.22 at 589 nm (sodium D line) has N=10N = 10 dynodes, each with secondary emission coefficient δ=4.2\delta = 4.2. The photocathode receives a steady photon flux Φ=5.0×106\Phi = 5.0\times10^{6} photons s−1^{-1} from a flame emission source. Overall gain:

G=δN=4.210G = \delta^{N} = 4.2^{10}

Compute stepwise: 4.22=17.644.2^2 = 17.64; 4.24=17.642=311.24.2^4 = 17.64^2 = 311.2; 4.25=311.2×4.2=13074.2^5 = 311.2\times4.2 = 1307; 4.210=13072=1.708×1064.2^{10} = 1307^2 = 1.708\times10^{6}.

G≈1.71×106G \approx 1.71\times10^{6}

Mean anode current:

ia=G e η Φ=(1.71×106)(1.602×10−19 C)(0.22)(5.0×106 s−1)i_a = G\,e\,\eta\,\Phi = (1.71\times10^{6})(1.602\times10^{-19}\ \text{C})(0.22)(5.0\times10^{6}\ \text{s}^{-1})

Multiply progressively: 1.71×106×1.602×10−19=2.740×10−131.71\times10^6 \times 1.602\times10^{-19} = 2.740\times10^{-13} C. Then 2.740×10−13×0.22=6.027×10−142.740\times10^{-13}\times 0.22 = 6.027\times10^{-14}. Finally 6.027×10−14×5.0×106=3.01×10−76.027\times10^{-14}\times 5.0\times10^{6} = 3.01\times10^{-7}.

ia≈3.0×10−7 A=300 nAi_a \approx 3.0\times10^{-7}\ \text{A} = 300\ \text{nA}

This is a comfortably measurable current with a simple current-to-voltage (transimpedance) amplifier — no photon counting needed at this flux. Photon-shot-noise SNR for a 1.0 s integration: photoelectron count npe=ηΦt=0.22×5.0×106×1.0=1.10×106n_{pe} = \eta\Phi t = 0.22\times5.0\times10^{6}\times1.0 = 1.10\times10^{6} photoelectrons.

SNR=npe=1.10×106≈1.05×103\mathrm{SNR} = \sqrt{n_{pe}} = \sqrt{1.10\times10^{6}} \approx 1.05\times10^{3}

A shot-noise-limited relative standard deviation of about 1/1049=0.095%1/1049 = 0.095\% — excellent precision, consistent with why flame and ICP emission measurements with adequate photon flux routinely achieve sub-percent RSD. Limiting-case check: if the dynode voltage were reduced so that δ→1\delta \to 1 (no multiplication), G→110=1G\to 1^{10}=1 and iai_a collapses to the bare photocathode current of ηeΦ≈1.76×10−13\eta e \Phi \approx 1.76\times10^{-13} A — far below what any practical electrometer can resolve without excessive averaging, which is precisely why PMTs, not bare photodiodes, are required for weak-line atomic emission and absorption work.

Example 2 — Applied: comparing PMT-sequential and CID-simultaneous ICP-OES detection limits. An ICP-OES method for manganese at 257.61 nm is run two ways on the same plasma and optics train, differing only in the detector. Configuration A (scanning PMT): the monochromator dwells on the Mn line for t=2.0t=2.0 s. The blank (pure matrix, no Mn) gives a mean photoelectron count Sbl=4.0×104S_{bl} = 4.0\times10^{4} with negligible dark counts and negligible amplifier noise (PMT dark counts at this wavelength and cooled housing are D≈50D\approx 50 electron-equivalents, negligible here). Sensitivity from a calibration curve is m=8.0×104m = 8.0\times10^{4} counts per (mg L−1^{-1}). Blank noise (shot-noise dominated since D,R2≪SblD, R^2 \ll S_{bl}):

sbl=Sbl=4.0×104=2.00×102 countss_{bl} = \sqrt{S_{bl}} = \sqrt{4.0\times10^{4}} = 2.00\times10^{2}\ \text{counts}

Detection limit at k=3k=3:

cL=k sblm=3×2.00×1028.0×104 counts/(mg L−1)=7.5×10−3 mg L−1c_L = \frac{k\,s_{bl}}{m} = \frac{3\times2.00\times10^{2}}{8.0\times10^{4}\ \text{counts}/(\text{mg L}^{-1})} = 7.5\times10^{-3}\ \text{mg L}^{-1}

Configuration B (CID array, same integration time, but read noise now matters): the same optical signal falls on a CID pixel, but the shorter effective single-pixel light-gathering geometry (an echelle image is spread over more resolution elements) gives a smaller Sbl=1.6×104S_{bl} = 1.6\times10^{4} electrons, plus dark current D=2.0×102D = 2.0\times10^{2} electrons and read noise R=25R = 25 electrons (rms) per non-destructive read, with sensitivity accordingly scaled to m=3.2×104m = 3.2\times10^{4} counts per (mg L−1^{-1}) (same fraction of total light, smaller absolute pixel signal).

sbl=Sbl+D+R2=1.6×104+2.0×102+625s_{bl} = \sqrt{S_{bl} + D + R^{2}} = \sqrt{1.6\times10^{4} + 2.0\times10^{2} + 625}
sbl=1.6825×104=1.297×102 electronss_{bl} = \sqrt{1.6825\times10^{4}} = 1.297\times10^{2}\ \text{electrons}
cL=3×1.297×1023.2×104=1.216×10−2 mg L−1≈1.2×10−2 mg L−1c_L = \frac{3\times1.297\times10^{2}}{3.2\times10^{4}} = 1.216\times10^{-2}\ \text{mg L}^{-1} \approx 1.2\times10^{-2}\ \text{mg L}^{-1}

Interpretation: the PMT configuration gives a somewhat better single-line detection limit (7.5 µg L−1^{-1} vs. 12 µg L−1^{-1}) because it collects more total light onto one large-area photocathode and adds essentially no electronic noise, whereas the array pixel, though only modestly noisier, splits the same optical throughput across many resolution elements and adds a genuine read-noise term. In practice this per-line penalty of roughly a factor of 1.6 in detection limit is precisely why simultaneous CID/CCD instruments are called "compromise" detectors — but the array configuration measures all 40+ other elements in the same 2.0 s exposure, while the PMT would need 40 separate 2.0 s dwells (80 s total, plus grating-move overhead), so total analysis throughput favors the array by well over an order of magnitude for multi-element work. Limiting-case check: if read noise R→0R\to 0 and dark current D→0D\to 0 for the array (an idealized, perfectly cooled, zero-noise sensor), sbl→Sbl=1.6×104=126.5s_{bl}\to\sqrt{S_{bl}} = \sqrt{1.6\times10^4}=126.5, giving cL=3×126.5/3.2×104=1.19×10−2c_L = 3\times126.5/3.2\times10^4 = 1.19\times10^{-2} mg L−1^{-1} — barely different from the computed 1.2×10⁻² mg L⁻¹ value above, confirming that in this example the array is already close to shot-noise-limited and the modest detection-limit penalty relative to the PMT is mostly a light-collection-per-pixel effect, not a read-noise artifact.

Example 3 — Exam/competition level: choosing integration strategy under a fixed photon budget, and the crossover between read-noise-limited and shot-noise-limited regimes. A CCD pixel monitoring a weak emission line has dark current Dr=15D_r = 15 electrons s−1^{-1} (already dominated by thermal generation despite cooling) and a fixed read noise R=8R = 8 electrons rms per readout (independent of exposure time, since it is generated once per readout event, not once per unit time). The photon signal arrives at a rate Sr=40S_r = 40 electrons s−1^{-1} at the pixel (post quantum efficiency). Full-well capacity is W=1.2×105W = 1.2\times10^{5} electrons. (a) Find the exposure time t∗t^{*} at which shot noise from signal-plus-dark equals the read noise contribution — the crossover from read-noise-limited to shot-noise-limited operation. Signal and dark electrons accumulated after time tt: S(t)=SrtS(t) = S_r t, D(t)=DrtD(t)=D_r t. Shot-noise variance from these Poisson sources is S(t)+D(t)=(Sr+Dr)tS(t)+D(t) = (S_r+D_r)t. Set this equal to R2R^2:

(Sr+Dr) t∗=R2(S_r + D_r)\,t^{*} = R^{2}
t∗=R2Sr+Dr=8240+15=6455=1.164 st^{*} = \frac{R^{2}}{S_r+D_r} = \frac{8^2}{40+15} = \frac{64}{55} = 1.164\ \text{s}

So for exposures shorter than about 1.16 s, read noise dominates the noise budget; for longer exposures, Poisson (shot + dark) noise dominates. (b) Find the maximum single-exposure time before the pixel saturates (reaches full well), and hence the maximum SNR obtainable from a single, unsaturated exposure.

tsat=WSr+Dr=1.2×10555=2.182×103 st_{sat} = \frac{W}{S_r+D_r} = \frac{1.2\times10^{5}}{55} = 2.182\times10^{3}\ \text{s}

At this exposure, signal electrons S(tsat)=40×2182=8.727×104S(t_{sat}) = 40\times2182 = 8.727\times10^4 and dark electrons D(tsat)=15×2182=3.273×104D(t_{sat}) = 15\times2182 = 3.273\times10^4; check the sum 8.727×104+3.273×104=1.200×105=W8.727\times10^4+3.273\times10^4 = 1.200\times10^5 = W ✓ (consistent by construction).

SNRmax=S(tsat)S(tsat)+D(tsat)+R2=8.727×1041.200×105+64\mathrm{SNR}_{max} = \frac{S(t_{sat})}{\sqrt{S(t_{sat})+D(t_{sat})+R^2}} = \frac{8.727\times10^{4}}{\sqrt{1.200\times10^{5}+64}}
SNRmax=8.727×1041.20006×105=8.727×104346.4=251.9\mathrm{SNR}_{max} = \frac{8.727\times10^{4}}{\sqrt{1.20006\times10^{5}}} = \frac{8.727\times10^{4}}{346.4} = 251.9

(c) Now suppose total available observing time is fixed at T=6.0×103T = 6.0\times10^{3} s and the pixel must be read out (incurring a fresh R=8R=8 e⁻ hit each time) every time it is emptied. Compare one long exposure of duration tsat=2182t_{sat}=2182 s repeated n=T/tsat≈2.75n = T/t_{sat} \approx 2.75 times (round to n=3n=3 separate full exposures, each shorter than tsatt_{sat} to fit exactly, i.e. t=T/3=2000t = T/3 = 2000 s) against co-adding many short, read-noise-dominated exposures of t=1.0t=1.0 s each (n=6000n=6000 reads). *Three long exposures, t=2000t=2000 s each:* per-exposure signal S=40×2000=8.0×104S=40\times2000=8.0\times10^4, dark D=15×2000=3.0×104D=15\times2000=3.0\times10^4. Total after co-adding 3 frames: Stot=2.4×105S_{tot}=2.4\times10^5, Dtot=9.0×104D_{tot}=9.0\times10^4, and total read-noise variance =nR2=3×64=192= n R^2 = 3\times64=192 (read noise adds in quadrature once per readout, and there are 3 readouts).

SNR=2.4×1052.4×105+9.0×104+192=2.4×1053.3019×105=2.4×105574.6=417.7\mathrm{SNR} = \frac{2.4\times10^{5}}{\sqrt{2.4\times10^{5}+9.0\times10^{4}+192}} = \frac{2.4\times10^{5}}{\sqrt{3.3019\times10^{5}}} = \frac{2.4\times10^{5}}{574.6} = 417.7

Six thousand 1.0 s exposures co-added: per-exposure signal S=40S=40, dark D=15D=15; totals Stot=40×6000=2.4×105S_{tot}=40\times6000=2.4\times10^5 electrons, Dtot=15×6000=9.0×104D_{tot}=15\times6000=9.0\times10^4 electrons (identical totals, since total exposure time is the same 6.0×1036.0\times10^3 s either way), but now read-noise variance accumulates n=6000n=6000 times:

nR2=6000×64=3.84×105n R^{2} = 6000\times 64 = 3.84\times10^{5}
SNR=2.4×1052.4×105+9.0×104+3.84×105=2.4×1057.14×105=2.4×105845.0=284.0\mathrm{SNR} = \frac{2.4\times10^{5}}{\sqrt{2.4\times10^{5}+9.0\times10^{4}+3.84\times10^{5}}} = \frac{2.4\times10^{5}}{\sqrt{7.14\times10^{5}}} = \frac{2.4\times10^{5}}{845.0} = 284.0

Interpretation: for identical total signal and dark-current accumulation, coadding many short reads (284.0) gives markedly worse SNR than a few long exposures (417.7), purely because each readout re-injects a full R2=64R^2=64 e⁻² of variance. This is the rigorous, numerical justification for the standard practice in low-light array spectroscopy (and imaging) of using the longest single exposure the full well and background will allow, rather than many short ones — and it is exactly the argument used to justify a CID's or CCD's minimum practical read count in trace ICP-OES work. Limiting-case check: if read noise R→0R\to 0, both strategies give identical SNR — SNR=Stot/Stot+Dtot=2.4×105/3.3×105=417.7\mathrm{SNR}=S_{tot}/\sqrt{S_{tot}+D_{tot}} = 2.4\times10^5/\sqrt{3.3\times10^5}=417.7 regardless of how the total time is divided — confirming that the entire penalty for many short exposures traces to the nR2nR^2 term, exactly as the algebra shows, and vanishes when R=0R=0.

Limiting cases, scope & common traps

Regime / caseBehaviorWhy
δ→1\delta \to 1 (no dynode gain)G→1G\to1; PMT signal collapses to bare photocathode current, unusable at low fluxmultiplication is the entire reason PMTs beat photodiodes for dim atomic lines
Read-noise-limited array (R2≫S,DR^2 \gg S,D)cL∝R/mc_L \propto R/m; longer integration does not helpmust reduce RR (better ADC/preamp) or increase sensitivity mm, not integration time
Shot-noise-limited array or PMT (S,D≫R2S,D \gg R^2)cL∝t/t=1/tc_L \propto \sqrt{t}/t = 1/\sqrt tclassic "detection limit improves with square root of time" rule, valid only past the read-noise floor
CCD near full wellrisk of blooming — charge overflows into neighboring pixels/columns, corrupting adjacent wavelengthsanti-blooming drains sacrifice some full well and linearity to prevent this
CID at high signalpixel individually reset/read at shorter time before saturating, no blooming into neighborsnon-destructive per-pixel access is architecturally immune to blooming
Sequential PMT, many elementsthroughput poor (each element needs its own dwell + grating move) but each measurement gets full photometric attention and superb dynamic range per linebest when element list is short or when the very lowest detection limits on a handful of lines are the priority
Simultaneous array, many elementsone exposure captures the whole spectral window; huge throughput gain for multi-element workbest for routine multi-element screening (environmental, industrial QC) even at a modest per-line detection-limit cost
PMT in photon-counting modebest sensitivity at very low flux, since pulse discrimination rejects sub-threshold noise pulsesrequires flux low enough that pulses don't pile up (dead-time correction needed)

The chemist behind it

The photomultiplier tube's essential architecture — a photoemissive cathode followed by cascaded secondary-emission stages — was developed at RCA in the 1930s, principally through the work of Vladimir Zworykin along with colleagues including Joseph Rajchman and George Morton, building on Zworykin's earlier iconoscope and photocell research for television. Zworykin, a Russian-American engineer trained under Boris Rosing in St. Petersburg before emigrating to the United States, was primarily driven by the demands of practical television camera and low-light imaging technology, not analytical chemistry — yet the PMT his group perfected became, within two decades, the single most important detector in flame and atomic absorption spectroscopy, and it remains in routine use in modern instruments to this day, more than eighty years after its invention. The later shift toward multichannel solid-state detection owes its analytical-chemistry translation largely to workers such as James Winefordner and collaborators at the University of Florida, who through the 1970s–1990s systematically evaluated and championed CID and CCD arrays for atomic emission spectrometry, publishing the noise-budget framework (signal, dark current, read noise combining in quadrature) that underlies essentially every modern array detector specification sheet. Winefordner's group's work on detector-limited figures of merit — bridging solid-state physics and quantitative trace analysis — established the vocabulary ("read-noise-limited," "shot-noise-limited," "charge-transfer efficiency") that analytical spectroscopists still use to specify and compare instruments.

Origins · 1 minVladimir K. Zworykin and how this was discovered
VK
Vladimir K. Zworykin
1888–1982 · RCA Laboratories, Princeton, New Jersey, developer of the modern photomultiplier and iconoscope television camera tube

History
- 1887 — Heinrich Hertz observes the photoelectric effect, the physical basis of every photocathode.
- 1905 — Einstein's photon quantization explanation of the photoelectric effect, E=hν−ϕE=h\nu-\phi, underlying photocathode response.
- 1934–1936 — Zworykin, Rajchman, and RCA colleagues develop the first practical multi-dynode photomultiplier tube.
- 1969 — Boyle and Smith at Bell Labs invent the charge-coupled device, initially conceived as a memory element.
- 1970s–1990s — CID and CCD arrays adapted and characterized for simultaneous atomic emission spectrometry, notably by Winefordner and coworkers, establishing the modern noise-budget framework for array detectors.
Why does a PMT need such a high, stable voltage supply?

Because gain scales as G=δNG=\delta^N and the per-stage secondary emission coefficient δ\delta itself increases with interstage voltage (roughly δ∝Vd0.7–0.8\delta\propto V_d^{0.7\text{–}0.8}), the overall gain is an extremely steep, nonlinear function of the total applied high voltage, often scaling as Vtotal6V_{total}^{6}–Vtotal8V_{total}^{8} across a typical 9–12 stage tube. A mere 1% drift in supply voltage can shift gain by several percent, directly corrupting quantitative absorbance or emission-intensity measurements. This is why PMT high-voltage supplies for analytical instruments are regulated far more tightly (often to millivolt precision) than supplies for other lab equipment, and why manufacturers specify gain stability, not just absolute gain, as a key performance parameter.

What exactly is "blooming" and why don't CIDs suffer from it?

Blooming occurs when a CCD pixel's potential well fills past its full-well capacity during integration; excess photogenerated charge has nowhere to go and spills into neighboring pixels along the charge-transfer direction, corrupting the signal at nearby wavelengths — a serious problem when a strong matrix or major-element line sits near a faint trace line in an emission spectrum. CCDs mitigate this with anti-blooming drain structures that shunt excess charge to the substrate, at some cost to fill factor and linearity. CIDs largely avoid the problem because each pixel can be individually and non-destructively read and reset well before saturation, and because the two-electrode CID pixel architecture does not require charge to be physically shifted through neighboring pixels during readout, removing the principal blooming pathway.

Is photon-counting mode better than analog current mode for a PMT?

Photon-counting mode — discriminating and counting individual anode pulses above a threshold — gives superior sensitivity at very low light levels because it rejects sub-threshold electronic noise pulses that would otherwise be integrated into an analog current measurement, effectively achieving a cleaner Poisson-statistics signal. Its limitation is pulse pile-up: at high photon flux, individual pulses begin to overlap in time and the counter undercounts, requiring a dead-time correction, and eventually analog current mode becomes both simpler and more accurate. Most modern instruments switch modes (or use logarithmic/dual-mode electrometers) depending on measured light level, choosing counting for trace work and analog integration for major-constituent lines.

How do detection limits reported for simultaneous versus sequential ICP-OES actually compare in practice?

Manufacturer-published detection limits for well-optimized simultaneous CID/CCD-based instruments are typically within a factor of 1.5–3 of dedicated scanning-PMT instruments for the same element and line, consistent with the noise-budget comparison worked through in this article's examples. The gap has narrowed substantially since the 1990s as CCD/CMOS read noise and dark current have improved and as echelle optics have made better use of the full 2-D detector area. For routine multi-element panels, the modest per-line detection-limit penalty of simultaneous detection is almost always outweighed by the throughput gain, which is why simultaneous instruments now dominate commercial ICP-OES sales despite sequential instruments retaining a niche for ultratrace, few-element applications.

Why is the relative variance of a PMT's single-photoelectron pulse dominated by the first dynode?

Because each dynode stage contributes a fractional variance term that is suppressed by an additional factor of 1/δ1/\delta relative to the previous stage in the cascade, the sum of contributions from the second, third, and later dynodes forms a rapidly converging geometric series, while the first dynode's contribution enters unsuppressed. The closed-form result, relative variance ≈1/(δ−1)\approx 1/(\delta-1), is set almost entirely by the first stage's own gain and statistics. This is why photomultiplier designers deliberately run the first dynode at a higher-than-average interstage voltage — a larger δ1\delta_1 narrows the single-photoelectron pulse-height distribution, improving pulse-height discrimination in photon-counting applications.

Can a single instrument use both a PMT and an array detector?

Yes — many practical spectrometers are hybrids. Some ICP-OES systems dedicate an array to the bulk simultaneous spectral window while adding a scanning PMT channel (or a movable exit slit and secondary photomultiplier) for a small number of especially demanding trace lines needing the best possible detection limit or the widest linear dynamic range. Double-monochromator and dual-detector AAS/ICP-OES combination instruments similarly pair PMTs for absorption channels with arrays for simultaneous emission background monitoring. The underlying principle, amplify-then-read versus integrate-then-read, is chosen independently for each measurement channel based on which figure of merit — dynamic range and per-line sensitivity, or throughput and simultaneity — matters most for that specific measurement.

What role does quantum efficiency play separately from gain?

Quantum efficiency η(λ)\eta(\lambda) sets how many photoelectrons (PMT) or electron–hole pairs (array) are generated per incident photon; gain GG (PMT) or full-well capacity and readout amplification (array) then determines how faithfully and with how much added noise that initial photoelectric event is converted into a measurable electrical signal. A detector can have superb quantum efficiency but poor overall performance if gain is too low (signal lost in downstream amplifier noise) or too high (saturation, nonlinearity); conversely a modest quantum efficiency with excellent gain and readout electronics, as in a well-designed PMT, can outperform a higher-QE but noisier device at low light levels. Comparing detectors on quantum efficiency alone, without considering the full noise budget derived in this article, is a common and consequential error.

Before this