Experiment: Artificial Transmutation and the Instruments of Nuclear Physics

Contents
  1. Conventions and constants
  2. Rutherford: the first artificial transmutation (1919)
  3. Cockcroft and Walton: lithium split by accelerated protons (1932)
  4. Lawrence and Livingston: the cyclotron (1932)
  5. Wilson: the expansion cloud chamber (1912)
  6. Glaser: the bubble chamber (1952)
  7. Charpak: the multiwire proportional chamber (1968)
  8. Summary of the evidence

Two questions had to be answered before particle physics could exist as a quantitative science. Can one nucleus be turned into another at will, and can the products be counted one at a time? This chapter reports the experiments that answered them. The first half follows the reaction: the transmutation of nitrogen by natural alpha particles, and then the disintegration of lithium by protons from a machine — the first nuclear reaction produced entirely by human means and, more importantly for this book, the first quantitative test of the mass–energy equivalence of Relativistic Dynamics in a nuclear reaction. The second half follows the instruments, from Wilson's cloud chamber through Glaser's bubble chamber to Charpak's multiwire proportional chamber, which replaced the photograph with an electronic address and made the modern event rate possible.

The chapter is the experimental floor of Part XI — Quantum Field Theory and the Standard Model. Its nuclear reactions belong to Nuclear Forces and Nuclear Structure, where the binding-energy systematics they sample are developed and where Chadwick's discovery of the neutron — made in 1932 with the same scintillation technique used here, in the same laboratory, in the same year as the lithium disintegration — is treated in full [Chadwick:1932b]. Its instruments are the ancestors of every detector in Cosmic Rays and Astroparticle Physics, Experiment: Deep Inelastic Scattering and Experiment: The Higgs Boson Discovery. Five of the six experiments below were recognised by a Nobel Prize in Physics: Wilson in 1927, Lawrence in 1939, Cockcroft and Walton in 1951, Glaser in 1960 and Charpak in 1992.

Conventions and constants

Three conventions are used throughout, and they are worth stating before any number appears, because every one of the primary papers below uses different ones.

Nuclear masses are the evaluated atomic masses of the AME 2020 tables [Wang:2021], and are converted to energy with the CODATA 2022 value of the unified atomic mass constant,

\begin{equation}\tag{117.1} u\,c^{2}=931.494\,\mathrm{MeV}\ec\qquad u=1.66053907\times 10^{-27}\,\mathrm{kg}\ec \end{equation}

to the precision needed here [Mohr:2025]. Because atomic rather than nuclear masses are used, the electrons balance on the two sides of every reaction written below, and no electron-mass correction is required.

Energies are quoted in electronvolts and in joules. The conversion is exact, since the elementary charge is a defining constant of the SI [Mohr:2025]:

\begin{equation}\tag{117.2} 1\,\mathrm{eV}=1.602176634\times 10^{-19}\,\mathrm{J}\ec\qquad 1\,\mathrm{MeV}=1.602176634\times 10^{-13}\,\mathrm{J}\ep \end{equation}

Two further combinations recur, and are taken from Equations (107.1) and (107.2) of Nuclear Forces and Nuclear Structure, where they are fixed once from the same adjustment: \(\hbar c=197.32698\,\mathrm{MeV}\,\mathrm{fm}\) and \(e^{2}/(4\pi\varepsilon_{0})=1.4399645\,\mathrm{MeV}\,\mathrm{fm}\). Both are ordinary SI quantities written in convenient units — the first is \(3.1615\times 10^{-26}\,\mathrm{J}\,\mathrm{m}\), the second \(2.3071\times 10^{-28}\,\mathrm{J}\,\mathrm{m}\) — and \(\hbar\) and \(c\) appear explicitly in every equation below.

Ranges are the third convention and the one most easily misread. The experimenters of 1912–1932 quoted the penetrating power of a particle as the thickness of air at atmospheric pressure that would stop it, and inserted solid absorbers whose stopping power had been calibrated against air. A length is not the physical quantity: what stops a particle is the number of electrons it must pass, so the invariant statement is an areal mass density. At \(288.15\,\mathrm{K}\) and \(101.325\,\mathrm{kPa}\) the ideal-gas law with the mean molar mass of dry air, \(M_{\text{air}}=28.9647\times 10^{-3}\,\mathrm{kg}/\mathrm{mol}\) [Navas:2024], gives \(\rho_{\text{air}}=pM_{\text{air}}/(RT) =1.2250\,\mathrm{kg}/\mathrm{m}^{3}\), so

\begin{equation}\tag{117.3} 10^{-2}\,\mathrm{m}\ \text{of air} =1.2250\times 10^{-2}\,\mathrm{kg}/\mathrm{m}^{2}\ec \end{equation}

and every range quoted below is given both ways. Imperial units, in which two of the primary papers are partly written, are converted on sight and not reproduced.

Remark 117.1 (Reading this chapter against the natural-units literature).

The nuclear and detector literature almost universally sets \(\hbar=c=1\), so that an energy, a momentum and a mass are all quoted in \(\mathrm{MeV}\) and a length in \(/\mathrm{MeV}\). This book does not: \(\hbar\) and \(c\) are carried explicitly, a momentum is written \(2.5\times 10^{8}\,\mathrm{eV}/c\) and a mass \(7\,\mathrm{GeV}/c^{2}\). To map a formula from that literature onto one here, restore the unique powers of \(\hbar\) and \(c\) that make the two sides dimensionally equal in SI; the combinations that do the work are the two quoted above, and the statement “a length of \(1\,\mathrm{fm}\) corresponds to an energy of \(197.3\,\mathrm{MeV}\)” means \(\hbar c/1\,\mathrm{fm}\) and nothing more. No derivation in this chapter is carried out in those units.

Rutherford: the first artificial transmutation (1919)

Tests Phenomenon 117.2.

Apparatus

The apparatus is a scintillation counter with a gas cell in front of it, and nothing else. A rectangular brass box about \(0.18\,\mathrm{m}\) long, with an internal cross-section of roughly \(0.06\,\mathrm{m}\times0.02\,\mathrm{m}\), was closed at one end by a silver foil whose stopping power had been calibrated as equivalent to about \(6\times 10^{-2}\,\mathrm{m}\) of air, that is \(7.4\times 10^{-2}\,\mathrm{kg}/\mathrm{m}^{2}\) by Equation (117.3). Immediately beyond the foil sat a screen of zinc sulfide, viewed through a microscope of low power by an observer who had dark-adapted for half an hour in a light-tight room.

Inside the box, on a slide that could be moved along its length without breaking the seal, a metal disc carried an active deposit collected from radon gas. The alpha emitter in that deposit is radium C\('\), known today as \(^{214}\mathrm{Po}\): a nuclide of half-life \(1.64\times 10^{-4}\,\mathrm{s}\) which is fed continuously by the longer-lived members of the chain above it and which therefore behaves as a steady source. Its alpha particle is among the most energetic available from any natural source in useful quantity — thorium C\('\) (\(^{212}\mathrm{Po}\)), reached through the thorium active deposit and worked with at the Cavendish in the same period, is more energetic still, at \(8.785\,\mathrm{MeV}\) by the same mass balance [Wang:2021] — and its energy is not an experimental input but a consequence of a mass difference. The alpha decay energy from the mass evaluation is \(Q_{\alpha}=7.8335\,\mathrm{MeV}\) [Wang:2021], and the recoiling \(^{210}\mathrm{Pb}\) takes a share fixed by momentum conservation, so the alpha particle carries

\begin{equation}\tag{117.4} E_{\alpha}=Q_{\alpha}\, \frac{m\!\left({}^{210}\mathrm{Pb}\right)} {m\!\left({}^{210}\mathrm{Pb}\right)+m\!\left({}^{4}\mathrm{He}\right)} =7.687\,\mathrm{MeV}=1.2316\times 10^{-12}\,\mathrm{J}\ec \end{equation}

which is the number quoted in every table of the period and the number used throughout this section. Its range in air at \(288.15\,\mathrm{K}\) and \(101.325\,\mathrm{kPa}\) is \(7.0\times 10^{-2}\,\mathrm{m}\), that is \(8.6\times 10^{-2}\,\mathrm{kg}/\mathrm{m}^{2}\).

The box could be filled with a dry gas of choice at a controlled pressure — oxygen, carbon dioxide, dry air, or nitrogen — and evacuated and refilled between runs. Gas purity is not a refinement here but the entire experiment: the effect looked for is the appearance of hydrogen nuclei, and hydrogen is the commonest contaminant in any vacuum system. The gas was therefore dried over phosphorus pentoxide, and the comparison between gases, rather than the absolute rate in any one of them, carries the argument.

The detector was the dark-adapted eye of the observer. A zinc sulfide screen converts the energy a stopping ion deposits into a flash of visible photons, and the flash from a single particle is at the edge of what the eye can see; the technique works only at rates of a few counts per minute, only in complete darkness, and only for observers who alternate in short shifts. Its virtue is that it counts single particles, and in 1919 nothing else did.

Procedure

Absorbing foils were interposed between source and screen and their stopping power expressed as an equivalent thickness of air, so that a count could be plotted against the distance the emerging particles were able to travel. Enough absorber was always present to stop every primary alpha particle: with the silver foil alone accounting for \(6\times 10^{-2}\,\mathrm{m}\) of air equivalent and the alpha range being \(7.0\times 10^{-2}\,\mathrm{m}\), a short additional path in the gas suffices, and whatever produced a scintillation on the screen therefore had a longer range than the alpha particles that made it.

The experiment is then a comparison between gases, made under four controls, each of which excludes a specific alternative.

First, the choice of filling gas. Oxygen and carbon dioxide contain no nitrogen and no hydrogen, and are chemically comparable stopping media; they are the null runs. Dry air is one fifth oxygen and four fifths nitrogen; pure nitrogen is the signal run. If the effect is proportional to the amount of nitrogen present it cannot be a property of the source, of the box, or of the screen.

Second, the source-to-screen distance and the gas pressure were varied independently. Both change the amount of gas traversed, and a genuine product of the gas must respond to their product and not to either separately.

Third, the absorber thickness was scanned. A particle with a definite energy has a definite range, so the count as a function of absorber falls to zero at a characteristic thickness; reading that end point is how an energy was measured before magnetic analysis was available in this energy range.

Fourth, the residual effect with the box evacuated was measured and subtracted. Hydrogen occluded in the source mount itself produces long-range particles by exactly the elastic mechanism the experiment must exclude, and that background is not small.

Observations and data

With oxygen or carbon dioxide in the box, no scintillations survived beyond the alpha range. With dry air, and more strongly with pure nitrogen, scintillations persisted at absorber thicknesses far beyond it — out to an air equivalent of about \(0.40\,\mathrm{m}\), that is \(0.49\,\mathrm{kg}/\mathrm{m}^{2}\) — and beyond the greatest range attainable by a hydrogen nucleus knocked on elastically by an alpha particle. The effect was small: a handful of scintillations per minute against an enormous number of incident alphas, which is why it took the full apparatus of scintillation counting to see it at all, and why the comparison between fillings rather than any absolute rate is the published result.

Table 117.1 collects the quantities the argument turns on. The historical rows are the observations; the computed rows are evaluated here from the mass evaluation [Wang:2021] and the kinematics derived below, and are what the observations have to be tested against. The counting uncertainties on the original scintillation tallies are Poissonian on tallies of order ten per run, that is of order \(30\,\mathrm{\%}\) per point; the ranges quoted from the paper are end points read off an absorption curve and are good to perhaps \(10\,\mathrm{\%}\). Neither is quoted to more precision than that here, and the original run-by-run tallies are not reproduced: what the paper establishes, and what is reproduced, is the ordering of the end points.

QuantityValue (eV)Value (SI)Kind, source
Alpha energy of radium~C$'$ ($^{214}\mathrm{Po}$)\(7.687\,\mathrm{MeV}\)\(1.2316\times 10^{-12}\,\mathrm{J}\)calc., [Wang:2021]
Range of those alphas in air\(8.6(9)\times 10^{-2}\,\mathrm{kg}/\mathrm{m}^{2}\)obs. 1919, [Rutherford:1919]
Longest range seen in nitrogen\(0.49(5)\,\mathrm{kg}/\mathrm{m}^{2}\)obs. 1919, [Rutherford:1919]
Elastic ceiling for a struck proton\(4.941\,\mathrm{MeV}\)\(7.916\times 10^{-13}\,\mathrm{J}\)calc., Equation (117.6)
Reaction energy $Q$ of ${}^{14}\mathrm{N}(\alpha,p){}^{17}\mathrm{O}$\(-1.1919\,\mathrm{MeV}\)\(-1.9096\times 10^{-13}\,\mathrm{J}\)calc., [Wang:2021]
Laboratory threshold\(1.5326\,\mathrm{MeV}\)\(2.4554\times 10^{-13}\,\mathrm{J}\)calc., Equation (117.11)
Proton energy at \(0\,^\circ\) to the beam\(5.929\,\mathrm{MeV}\)\(9.499\times 10^{-13}\,\mathrm{J}\)calc., Equation (117.12)
Proton energy at \(90\,^\circ\)\(4.423\,\mathrm{MeV}\)\(7.086\times 10^{-13}\,\mathrm{J}\)calc., Equation (117.12)
Proton energy at \(180\,^\circ\)\(3.299\,\mathrm{MeV}\)\(5.286\times 10^{-13}\,\mathrm{J}\)calc., Equation (117.12)
Two-branch forks among alpha tracks\(8(3)\) in \(415000\)obs. 1925, [Blackett:1925]
The range and energy comparison that carries Rutherford's argument. Rows marked obs. are the experimental values, quoted from the primary paper; rows marked calc. are computed in this chapter from the evaluated atomic masses and the elastic and two-body kinematics derived below, and involve no fitted quantity. Air equivalents are converted to areal density by Equation (117.3). The elastic ceiling is the largest energy an alpha particle can give a free proton in any collision whatever; the forward reaction proton exceeds it, and that is the whole observation. Bracketed uncertainties on the observed rows are the \(10\,\mathrm{\%}\) range-end-point reading error and the Poisson error on the fork count discussed in the text; the primary papers quote none, so these are estimates from the stated technique and are labelled as such. The computed rows inherit only the mass-evaluation uncertainty, which is of order \(10\,\mathrm{eV}\) and therefore falls six orders of magnitude below the last digit shown.

The uncertainty on the calculated rows deserves a word, because it is negligible and it is worth saying why. The AME 2020 evaluation quotes the masses of \(^{1}\mathrm{H}\), \(^{4}\mathrm{He}\), \(^{14}\mathrm{N}\) and \(^{17}\mathrm{O}\) to eleven significant figures, so the propagated uncertainty on the \(Q\) value of Equation (117.10) is at the level of electronvolts — six orders of magnitude below the \(0.1\,\mathrm{MeV}\) with which a range end point can be read. In this comparison the theory is exact and the measurement is not, which is the reverse of the usual situation and is what makes a nineteen-nineteen range measurement still worth reporting.

Interpretation

The long-range particles are hydrogen nuclei, and they cannot have come from hydrogen in the gas, for two independent reasons. Their range exceeds what an elastic collision can give — the ceiling in Table 117.1 is \(4.941\,\mathrm{MeV}\) and particles of longer range than that were counted — and the yield tracks the nitrogen content of the filling rather than the total amount of gas. Rutherford drew the conclusion in the plainest possible terms: the nitrogen atom is disintegrated.

What he could not establish was the fate of the alpha particle. Two possibilities remained open in 1919. The alpha might knock a proton out and continue, in which case four objects leave the collision point (alpha, proton, residual \(^{13}\mathrm{C}\), and nothing else) and a photograph of the event would show three tracks emerging. Or the alpha might be captured, in which case only two objects leave and a photograph shows two. That was settled six years later by Blackett, who photographed the collisions in a Wilson chamber (Section 117.5) and found, among some \(415000\) alpha tracks recorded on \(23000\) stereoscopic photographs, eight forks with exactly two branches and no third. A two-branch fork means the alpha particle did not survive the encounter: it was captured, and the reaction is

\begin{equation}\tag{117.5} {}^{14}\mathrm{N}+{}^{4}\mathrm{He}\longrightarrow {}^{17}\mathrm{O}+{}^{1}\mathrm{H}\ep \end{equation}

The counting is itself the measurement: eight events in \(415000\) tracks, each track traversing a few centimetres of air, is a yield of order \(10^{-5}\) per alpha, and a Poisson uncertainty of \(\sqrt{8}\approx3\) on eight events is why Blackett's number is a statement about the mechanism and not about the cross-section.

That the two branches also close the momentum triangle is what excludes the remaining alternative, that the alpha was captured and a neutral object emitted alongside the proton. In a stereoscopic pair the directions of both branches are measured, and the vector sum of the two final momenta must equal the incident momentum. Blackett's forks close; there is no missing momentum to assign to an unseen particle. The same argument run in reverse, on a reaction where the triangle does not close, is what Chadwick used in 1932 to establish the neutron [Chadwick:1932b], and it is the ancestor of every missing-momentum measurement in Experiment: The Higgs Boson Discovery.

This is the first transmutation of one element into another by deliberate human action, and it was done with a naturally occurring projectile whose energy the experimenter could not change. Removing that restriction is the content of Section 117.3.

Phenomenon 117.2 (Nitrogen is transmuted by alpha bombardment).

Bombarding nitrogen with alpha particles from a natural source ejects hydrogen nuclei whose range exceeds that of any proton an alpha particle can set in motion by an elastic collision [Rutherford:1919]. The alpha particle is absorbed in the process, so that the reaction is the capture Equation (117.5) and not a scattering [Blackett:1925]. Rests on Equation (117.5).

Derivation. Derives Phenomenon 117.2. Three independent statements have to be checked: that an elastic collision cannot produce the observed particles, that the capture reaction is energetically open to alphas of \(7.687\,\mathrm{MeV}\), and that when it is open the protons it produces do exceed the elastic ceiling.

The elastic ceiling. Take a projectile of mass \(m_{1}\) and kinetic energy \(E\) striking a free target of mass \(m_{2}\) at rest, non-relativistically — and here \(E/m_{1}c^{2}=7.687\,\mathrm{MeV}/3727.4\,\mathrm{MeV} =2.06\times 10^{-3}\), with the alpha-particle mass energy from [Mohr:2025], so the non-relativistic treatment is ample. Conservation of momentum and of kinetic energy in the head-on case give the target the energy

\begin{equation}\tag{117.6} E_{2}^{\max}=\frac{4m_{1}m_{2}}{\left(m_{1}+m_{2}\right)^{2}}\,E\ec \end{equation}

which is the largest transfer any elastic collision can achieve, at any angle, for any impact parameter, and independently of the force law: only the two masses enter. For an alpha particle on a proton the mass numbers alone would give \(16/25=0.64\), but the ceiling is a number the argument leans on, so it is evaluated from the same atomic masses as everything else: \(4m_{1}m_{2}/(m_{1}+m_{2})^{2} =0.64274\) and

\begin{equation}\tag{117.7} E_{2}^{\max}=0.64274\times7.687\,\mathrm{MeV} =4.941\,\mathrm{MeV}=7.916\times 10^{-13}\,\mathrm{J}\ep \end{equation}

Protons of longer range than this were seen, and no elastic process can make them. The same formula applied to an alpha on a nitrogen nucleus gives a factor \(0.6915\), so elastic recoils of nitrogen carry up to \(5.32\,\mathrm{MeV}\) — energetic, but a nitrogen nucleus of that energy has a range in air of a few millimetres and cannot reach the screen.

The reaction energy. Form the mass balance of Equation (117.5) from the evaluated atomic masses [Wang:2021]. Summing the two sides,

\begin{align} \tag{117.8} m\!\left({}^{14}\mathrm{N}\right)+m\!\left({}^{4}\mathrm{He}\right) &=18.005677\,u\ec\\ \tag{117.9} m\!\left({}^{17}\mathrm{O}\right)+m\!\left({}^{1}\mathrm{H}\right) &=18.006957\,u\ec \end{align}

so the products are heavier than the reactants by \(0.001280\,u=2.1247\times 10^{-30}\,\mathrm{kg}\) and the reaction energy is negative,

\begin{equation}\tag{117.10} Q=\left[m\!\left({}^{14}\mathrm{N}\right) +m\!\left({}^{4}\mathrm{He}\right) -m\!\left({}^{17}\mathrm{O}\right) -m\!\left({}^{1}\mathrm{H}\right)\right]c^{2} =-1.1919\,\mathrm{MeV}=-1.9096\times 10^{-13}\,\mathrm{J}\ep \end{equation}

An endothermic reaction needs more than \(\abs{Q}\) in the laboratory, because the centre of mass must keep moving and its kinetic energy is unavailable; demanding that the invariant available energy suffice to create the products at rest in their own centre of mass gives the threshold

\begin{equation}\tag{117.11} E_{\text{th}}=\abs{Q} \left(1+\frac{m\!\left({}^{4}\mathrm{He}\right)} {m\!\left({}^{14}\mathrm{N}\right)}\right) =1.5326\,\mathrm{MeV}=2.4554\times 10^{-13}\,\mathrm{J}\ep \end{equation}

Radium C\('\) supplies \(7.687\,\mathrm{MeV}\), five times the threshold, so the reaction is open.

The proton energy. Being open is not enough; the surplus has to end up in the proton in the observable amount. For a two-body reaction \(1+2\to3+4\) with target \(2\) at rest, non-relativistic throughout, conservation of energy and of momentum give the energy of product \(3\) at laboratory angle \(\theta\) as

\begin{equation}\tag{117.12} \sqrt{E_{3}}=v+\sqrt{v^{2}+w}\ec\qquad v=\frac{\sqrt{m_{1}m_{3}E_{1}}}{m_{3}+m_{4}}\cos\theta\ec\qquad w=\frac{m_{4}Q+\left(m_{4}-m_{1}\right)E_{1}}{m_{3}+m_{4}}\ec \end{equation}

where \(E_{1}\) is the projectile energy and the masses may be taken in any consistent unit, the ratios being all that enter. Both \(v^{2}\) and \(w\) carry the SI dimension of energy, so \(\sqrt{E_{3}}\) is an energy to the half power on both sides and the relation is dimensionally homogeneous.

Equation (117.12) is used six times in this chapter, so it is derived rather than quoted. Let \(p_{1},p_{3},p_{4}\) be the laboratory momenta and \(\theta\) the angle of product \(3\) to the beam. Momentum conservation gives \(\vect{p}_{4}=\vect{p}_{1}-\vect{p}_{3}\), hence

\begin{equation}\tag{117.13} p_{4}^{2}=p_{1}^{2}+p_{3}^{2}-2p_{1}p_{3}\cos\theta\ec \end{equation}

and energy conservation, with the reaction energy \(Q\) released into kinetic energy, gives \(E_{4}=E_{1}+Q-E_{3}\). Non-relativistically \(p^{2}=2mE\) throughout — legitimate here because \(E/mc^{2}\le2.1\times 10^{-3}\) for every particle involved — so Equation (117.13) becomes

\begin{equation}\tag{117.14} 2m_{4}\left(E_{1}+Q-E_{3}\right) =2m_{1}E_{1}+2m_{3}E_{3} -2\sqrt{4m_{1}m_{3}E_{1}E_{3}}\,\cos\theta\ep \end{equation}

Dividing by \(2\left(m_{3}+m_{4}\right)\) and collecting the terms in \(E_{3}\) and in \(\sqrt{E_{3}}\) turns this into \(E_{3}-2v\sqrt{E_{3}}-w=0\) with exactly the \(v\) and \(w\) of Equation (117.12); solving the quadratic in \(\sqrt{E_{3}}\) and discarding the negative root, which would make \(E_{3}\) negative whenever \(w>0\), gives \(\sqrt{E_{3}}=v+\sqrt{v^{2}+w}\). With \(m_{1}=m({}^{4}\mathrm{He})\), \(m_{2}=m({}^{14}\mathrm{N})\), \(m_{3}=m({}^{1}\mathrm{H})\), \(m_{4}=m({}^{17}\mathrm{O})\), \(E_{1}=7.687\,\mathrm{MeV}\) and \(Q\) from Equation (117.10), the bracket \(w\) is \(4.423\,\mathrm{MeV}\) — which is \(E_{p}(90\,^\circ)\) itself, since \(v\) vanishes there — and

\begin{equation}\tag{117.15} E_{p}(0\,^\circ)=5.929\,\mathrm{MeV}\ec\qquad E_{p}(90\,^\circ)=4.423\,\mathrm{MeV}\ec\qquad E_{p}(180\,^\circ)=3.299\,\mathrm{MeV}\ep \end{equation}

The forward protons exceed the elastic ceiling of \(4.941\,\mathrm{MeV}\) by \(0.988\,\mathrm{MeV}\), a margin of \(20\,\mathrm{\%}\) in energy and considerably more in range, since a proton's range grows faster than its energy. Backward protons do not exceed it, which is consistent: the effect appears in the forward hemisphere, is a small fraction of the total yield, and is exactly what a scintillation screen placed beyond an absorber selects.

All three statements come from the same mass table and from elementary kinematics, and none of them is available to a hydrogen-contamination explanation, which is bounded by Equation (117.6) whatever the amount of hydrogen present.

Remark 117.3 (What this experiment does not measure).

It is worth being exact about the limits, because the experiment is often described as though it had measured a cross-section. It did not. The observable is a count of scintillations above a range threshold, with an efficiency set by a human eye, a solid angle known to perhaps a factor of two, and a source strength known from an ionisation measurement. What is established at high confidence is a qualitative statement of a kind that no amount of precision could improve on — that particles of a range no elastic collision can produce are emitted from nitrogen and from no other gas tried — and the identification of the reaction had to wait six years for a different instrument. The quantitative rows of Table 117.1 are the modern calculation of what should have been seen, not a modern measurement of what was.

Primary references

The observation is Rutherford's [Rutherford:1919], the fourth and last paper of a series in the Philosophical Magazine on the collision of alpha particles with light atoms; the first three establish the scintillation technique and the elastic behaviour of hydrogen and of the other light gases, which is what makes the fourth readable as a controlled comparison rather than as an isolated claim. It is worth reading for its caution: Rutherford states the anomaly, the controls that exclude contamination, and explicitly declines to say what became of the alpha particle.

Blackett's paper [Blackett:1925] is the identification, and it is a different kind of document — an exhaustive scan of a photographic record, with the eight forks reproduced and measured. The two together are the pattern this chapter follows throughout: a counting experiment establishes that something happens, and an imaging experiment establishes what. Wilson's instrument [Wilson:1912], reported in Section 117.5, is what made the second half possible, and it is seven years older than the anomaly it resolved.

Cockcroft and Walton: lithium split by accelerated protons (1932)

Tests Phenomena 40.6, 117.4 and 117.5. Assuming Theorem 9.76.

Apparatus

The apparatus divides into a high-voltage generator, an accelerating tube, and a scintillation counter of the kind already described in Section 117.2.

The generator is the part that carries the names. A transformer supplied an alternating potential of peak value \(V_{0}\) of order \(100\,\mathrm{kV}\), and a cascade of capacitors and rectifiers converted it to a steady direct potential several times larger. The principle is worth stating because it is the reason the machine existed. In a cascade of \(n\) stages, the alternating column of capacitors is charged through the rectifiers on alternate half cycles, each stage of the direct column ends up at \(2V_{0}\) above the one below it, and the output is

\begin{equation}\tag{117.16} V_{\text{out}}=2nV_{0}\ep \end{equation}

With four stages and \(V_{0}=100\,\mathrm{kV}\) this is \(800\,\mathrm{kV}\), which is what the Cavendish machine reached; in practice \(100\,\mathrm{kV}\) to \(700\,\mathrm{kV}\) was used. No component ever holds more than \(2V_{0}\) across it, which is what makes the arrangement buildable at all, and the same circuit — now called a Cockcroft–Walton generator — is still the first stage of many accelerator chains.

The direct potential was applied across a sealed accelerating tube built of glass cylinders about a metre long and \(0.3\,\mathrm{m}\) in diameter, with intermediate electrodes so that the drop was shared and no single gap had to hold the whole voltage. That constraint is physical, not merely engineering: the dielectric strength of air at atmospheric pressure is about \(3\times 10^{6}\,\mathrm{V}/\mathrm{m}\), so \(800\,\mathrm{kV}\) requires a clearance of at least \(0.27\,\mathrm{m}\) in air even between ideal electrodes, and several times that between real ones with sharp edges and dust. It is the reason the electrostatic route to high energies has a ceiling of a few megavolts, and hence the reason for Section 117.4.

Protons were drawn from a hydrogen discharge at the top of the tube, a canal-ray source running at a few kilovolts, delivering a beam current of order \(10^{-5}\,\mathrm{A}\) at the target with the tube continuously pumped to a residual pressure of order \(10^{-3}\,\mathrm{Pa}\). The target was lithium metal, deposited on a backing and set at an angle to the beam so that the products could leave sideways. They left through a mica window whose stopping power had been calibrated against air, and were counted as scintillations on a zinc-sulfide screen under a microscope, exactly as in Section 117.2; further mica sheets placed in front of the screen measured their range.

Procedure

The measurement has three parts, and only the third identifies the reaction.

The accelerating potential was set and measured — by a sphere gap and a resistance chain — and the scintillation rate recorded as a function of it, so that the yield curve could be followed downwards towards the point where the effect disappeared into the background. The importance of pushing the curve downwards rather than upwards is that the interesting statement is a lower limit on where the reaction is still detectable, which is a statement about the barrier and not about the apparatus.

The range of the emitted particles was then measured by absorption in mica, converted to an equivalent range in air by the calibration, and hence to an energy through the range–energy relation for helium nuclei established with natural alpha sources. The particles' behaviour in absorbers — the shape of the absorption curve near its end point, and the specific ionisation, which follows Equation (115.15) — identified them as helium nuclei rather than protons or electrons.

Finally, two screens were placed on opposite sides of the target and their scintillations examined for coincidence, with an observer at each. This is the step that turns a detection into an identification. If the particles are emitted singly, the two observers see uncorrelated flashes at a rate set by chance; if they are emitted in pairs travelling in opposite directions, the coincidence rate exceeds the accidental rate by a large factor. The kinematics of the proposed reaction makes that test sharp, as Equation (117.24) below shows.

Observations and data

Scintillations appeared with accelerating potentials as low as about \(125\,\mathrm{kV}\), and the yield rose rapidly as the potential was raised. The particles had a range in air of about \(8.4\times 10^{-2}\,\mathrm{m}\), that is \(0.103\,\mathrm{kg}/\mathrm{m}^{2}\), which for a helium nucleus corresponds to a kinetic energy of about \(8.6\,\mathrm{MeV}\), that is \(1.38\times 10^{-12}\,\mathrm{J}\); their behaviour in absorbers was that of ordinary alpha particles. The coincidence test showed them emitted in pairs, in opposite directions.

Table 117.2 sets the measured quantities beside the predictions computed here. Their accuracy is limited on the experimental side by the reading of a range end point, good to a few per cent, and by the calibration of the range–energy relation for helium nuclei; on the calculated side it is limited by nothing that matters, the mass-evaluation uncertainty on the \(Q\) value being of order \(10\,\mathrm{eV}\).

QuantityValue (eV)Value (SI)Kind, source
Lowest potential showing the effect\(125(5)\,\mathrm{kV}\)obs. 1932, [Cockcroft:1932a]
Range of the products in air\(0.103(4)\,\mathrm{kg}/\mathrm{m}^{2}\)obs. 1932, [Cockcroft:1932c]
Energy inferred from that range\(8.6(3)\,\mathrm{MeV}\)\(1.38(5)\times 10^{-12}\,\mathrm{J}\)obs. 1932, [Cockcroft:1932c]
Mass lost, ${}^{7}\mathrm{Li}+{}^{1}\mathrm{H}\to2\,{}^{4}\mathrm{He}$\(3.0922\times 10^{-29}\,\mathrm{kg}\)calc., [Wang:2021]
Energy released, $Q=\Delta m\,c^{2}$\(17.3462\,\mathrm{MeV}\)\(2.7792\times 10^{-12}\,\mathrm{J}\)calc., Equation (117.21)
Alpha energy at \(90\,^\circ\) to the beam\(8.785\,\mathrm{MeV}\)\(1.4076\times 10^{-12}\,\mathrm{J}\)calc., Equation (117.12)
Alpha energy at \(0\,^\circ\)\(9.639\,\mathrm{MeV}\)\(1.5443\times 10^{-12}\,\mathrm{J}\)calc., Equation (117.12)
Alpha energy at \(180\,^\circ\)\(8.008\,\mathrm{MeV}\)\(1.2830\times 10^{-12}\,\mathrm{J}\)calc., Equation (117.12)
Departure of the pair from collinearity\(5.3\,^\circ\)calc., Equation (117.24)
Cockcroft and Walton's lithium disintegration: measured against calculated. The calculated alpha energies come from the two-body kinematics of Equation (117.12) applied to Equation (117.17) at a proton energy of \(300\,\mathrm{keV}\), which is representative of the potentials at which the range measurement was made; the mass difference is from the evaluated atomic masses. The last row is the reason the coincidence test works: at these energies the two alpha particles are back to back in the laboratory to within a few degrees. Bracketed uncertainties on the observed rows are estimates from the stated technique — a few per cent on a range end point and on its conversion to energy, and the sphere-gap calibration on the potential — since the primary papers quote none; the onset potential is in addition a limit of detection set by the scintillation background, not a physical threshold. The calculated rows carry the mass-evaluation uncertainty of order \(10\,\mathrm{eV}\), far below the last digit shown.

The second table is the sub-barrier behaviour, and it is computed rather than observed. Table 117.3 evaluates the semiclassical barrier transmission for a proton incident on a lithium nucleus at each of the accelerating potentials used, from the Gamow exponent of Equation (9.93). The point of tabulating it is that the numbers are not small in the way one expects of a tunnelling problem: the barrier here is low and thin, the exponent is of order unity to five, and a transmission of \(6\times 10^{-3}\) at the lowest potential is entirely consistent with a countable rate. It is the classical prediction — exactly zero at every potential in the table — that the observation contradicts.

$V$$E$ (centre of mass)$b$$2\gamma$$T$$Y/Y_{0}$
\(125\,\mathrm{kV}\)\(109.3\,\mathrm{keV}\)\(39.5\,\mathrm{fm}\)\(5.04\)\(6.5\times 10^{-3}\)\(1.00\)
\(150\,\mathrm{kV}\)\(131.2\,\mathrm{keV}\)\(32.9\,\mathrm{fm}\)\(4.32\)\(1.33\times 10^{-2}\)\(2.24\)
\(200\,\mathrm{kV}\)\(174.9\,\mathrm{keV}\)\(24.7\,\mathrm{fm}\)\(3.31\)\(3.64\times 10^{-2}\)\(7.02\)
\(250\,\mathrm{kV}\)\(218.6\,\mathrm{keV}\)\(19.8\,\mathrm{fm}\)\(2.63\)\(7.18\times 10^{-2}\)\(15.4\)
\(300\,\mathrm{kV}\)\(262.3\,\mathrm{keV}\)\(16.5\,\mathrm{fm}\)\(2.14\)\(1.18\times 10^{-1}\)\(27.5\)
\(400\,\mathrm{kV}\)\(349.8\,\mathrm{keV}\)\(12.4\,\mathrm{fm}\)\(1.46\)\(2.31\times 10^{-1}\)\(62.4\)
\(500\,\mathrm{kV}\)\(437.2\,\mathrm{keV}\)\(9.9\,\mathrm{fm}\)\(1.02\)\(3.60\times 10^{-1}\)\(109\)
\(700\,\mathrm{kV}\)\(612.1\,\mathrm{keV}\)\(7.1\,\mathrm{fm}\)\(0.48\)\(6.17\times 10^{-1}\)\(225\)
Barrier transmission for a proton on ${}^{7}\mathrm{Li}$, computed here in SI from Equation (9.93) with a contact radius $R=4.08\,\mathrm{fm}$ and the reduced mass $\mu=1.4626\times 10^{-27}\,\mathrm{kg}$. $V$ is the accelerating potential, $E$ the energy in the centre of mass, $b$ the classical turning point, $2\gamma$ the Gamow exponent and $T=\ee^{-2\gamma}$ the transmission. The last column is the relative thick-target yield of Equation (117.31), normalised to the lowest potential; it is the quantity a yield curve actually measures, and it rises far more steeply than $T$ alone. The classical value of $T$ is zero in every row, since the barrier height is \(1.06\,\mathrm{MeV}\). Nothing here is a measurement, so no row carries an experimental uncertainty; what the numbers inherit is the \(4\,\mathrm{\%}\) uncertainty of the contact radius $R$, which moves $2\gamma$ in the first row by about \(0.05\).

Interpretation

Two alpha particles emitted back to back, with nothing else, identifies the reaction uniquely as

\begin{equation}\tag{117.17} {}^{7}\mathrm{Li}+{}^{1}\mathrm{H}\longrightarrow {}^{4}\mathrm{He}+{}^{4}\mathrm{He}\ep \end{equation}

This is the first nuclear reaction driven by artificially accelerated particles: unlike Section 117.2, the projectile energy is now a knob on the apparatus, and the yield curve is a measurement rather than an accident of which radioactive source was to hand.

Its importance for this book, however, is quantitative. The energy released was measured directly, from the range of the products, and compared with the energy predicted by the difference of the masses — a prediction that has nothing to do with nuclear physics and everything to do with Phenomenon 40.6. The two agree. Cockcroft and Walton made the comparison in 1932 with the mass-spectrograph values then available, which gave a mass difference of about \(0.018\,u\) and hence an energy release close to \(17\,\mathrm{MeV}\), against the \(2\times8.6\,\mathrm{MeV}=17.2\,\mathrm{MeV}\) they measured; with the modern evaluation the prediction is \(17.3462\,\mathrm{MeV}\) (Equation (117.21)) and the kinematically correct comparison at \(90\,^\circ\) is \(8.785\,\mathrm{MeV}\) per alpha against a measured \(8.6\,\mathrm{MeV}\), a discrepancy of \(2\,\mathrm{\%}\) and entirely within the accuracy of a range measurement. This is the first quantitative confirmation of mass–energy equivalence in a nuclear reaction, and it was obtained from a table of atomic weights and a ruler laid against an absorption curve.

The comparison is worth appreciating for how little it assumes. The left-hand side of Equation (117.21) is measured by a mass spectrograph, an instrument that deflects ions in electric and magnetic fields and knows nothing of nuclear forces. The right-hand side is measured by the distance an alpha particle travels in mica. The only bridge between them is \(E=mc^{2}\). Had the two disagreed, no nuclear model would have been at fault; the relativistic energy–momentum relation of Relativistic Dynamics would have been.

A second reading of the same data was equally consequential. The reaction was observed at proton energies an order of magnitude below the height of the Coulomb barrier between a proton and a lithium nucleus, which classical mechanics forbids outright. The experiment is therefore also a direct confirmation of barrier penetration [Gamow:1928] [Gurney:1928] [Gurney:1929], whose other early evidence — alpha decay — is developed in Theorem 107.63 of Nuclear Forces and Nuclear Structure. It is a cleaner confirmation than alpha decay in one respect and a dirtier one in another. Cleaner, because the energy of the incident particle is set by the experimenter and can be varied, so the exponential dependence predicted by Equation (9.87) can be traced out rather than inferred from a correlation across nuclides. Dirtier, because the absolute yield involves the nuclear matrix element and the thick-target energy loss as well as the barrier, so only the shape of the curve is a clean test.

Phenomenon 117.4 (Lithium disintegrates into two alpha particles).

Lithium bombarded with protons accelerated through a few hundred kilovolts emits helium nuclei in pairs, in nearly opposite directions, each of kinetic energy near \(8.6\,\mathrm{MeV}\) [Cockcroft:1932a] [Cockcroft:1932c]. The total energy released is equal, within the accuracy of the measurement, to the mass lost by the system, converted at the rate \(E=mc^{2}\). Rests on Equations (40.6), (117.12) and (117.17).

Derivation. Derives Phenomenon 117.4. Write the mass balance of Equation (117.17) with the evaluated atomic masses [Wang:2021]:

\begin{align} \tag{117.18} m\!\left({}^{7}\mathrm{Li}\right)+m\!\left({}^{1}\mathrm{H}\right) &=7.0160034\,u+1.0078250\,u =8.0238284\,u\ec\\ \tag{117.19} 2\,m\!\left({}^{4}\mathrm{He}\right) &=2\times4.0026033\,u=8.0052065\,u\ep \end{align}

The system is lighter after the reaction than before by

\begin{equation}\tag{117.20} \Delta m=0.0186220\,u=3.0922\times 10^{-29}\,\mathrm{kg}\ec \end{equation}

and by Equation (117.1) the energy released is

\begin{equation}\tag{117.21} Q=\Delta m\,c^{2} =3.0922\times 10^{-29}\,\mathrm{kg} \times\left(2.99792458\times 10^{8}\,\mathrm{m}/\mathrm{s}\right)^{2} =2.7792\times 10^{-12}\,\mathrm{J}=17.3462\,\mathrm{MeV}\ep \end{equation}

Note that the electron count balances — three electrons on each side, since lithium is \(Z=3\) and each helium is \(Z=2\) — so atomic masses may be used without correction, and that \(\Delta m\) is \(0.23\,\mathrm{\%}\) of the total mass, which is why a mass spectrograph of the period could see it.

Now share it. Because the two products are identical, the reaction is symmetric in the centre of mass and each alpha particle takes half the available energy there; in the laboratory the pair is boosted along the beam. Applying Equation (117.12) with \(m_{1}=m({}^{1}\mathrm{H})\), \(m_{2}=m({}^{7}\mathrm{Li})\), \(m_{3}=m_{4}=m({}^{4}\mathrm{He})\), \(Q\) from Equation (117.21) and \(E_{1}=300\,\mathrm{keV}\):

\begin{equation}\tag{117.22} E_{\alpha}(0\,^\circ)=9.639\,\mathrm{MeV}\ec\quad E_{\alpha}(90\,^\circ)=8.785\,\mathrm{MeV}\ec\quad E_{\alpha}(180\,^\circ)=8.008\,\mathrm{MeV}\ec \end{equation}

the three summing correctly in pairs to \(Q+E_{1}=17.646\,\mathrm{MeV}\) for the forward–backward pair. The screens in the coincidence experiment sat transverse to the beam, where the prediction is \(8.785\,\mathrm{MeV}\) per particle, and the measured \(8.6\,\mathrm{MeV}\) agrees to \(2\,\mathrm{\%}\).

Two checks complete the argument. First, that the crude statement “each alpha takes half” is a good approximation and why: with \(E_{1}\ll Q\) the incident momentum is small beside the momenta of the products. An alpha particle of \(8.785\,\mathrm{MeV}\) has

\begin{equation}\tag{117.23} p_{\alpha}c=\sqrt{E_{\alpha}^{2}+2E_{\alpha}m_{\alpha}c^{2}} =256.1\,\mathrm{MeV}\ec\qquad p_{\alpha}=1.368\times 10^{-19}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}\ec \end{equation}

using \(m_{\alpha}c^{2}=3727.4\,\mathrm{MeV}\) [Mohr:2025], while a proton of \(300\,\mathrm{keV}\) has \(p_{p}c=23.7\,\mathrm{MeV}\), that is \(1.27\times 10^{-20}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}\). Second, the angular consequence: if one alpha leaves at \(90\,^\circ\), the other must carry the whole incident longitudinal momentum, so the pair departs from collinearity by

\begin{equation}\tag{117.24} \delta=\arctan\frac{p_{p}}{p_{\alpha}} =\arctan\frac{23.7\,\mathrm{MeV}/c}{256.1\,\mathrm{MeV}/c} =5.3\,^\circ\ep \end{equation}

Five degrees is small enough that two screens set at \(180\,^\circ\) both intercept the pair, which is why the coincidence test worked, and large enough to have been detectable had the apparatus been built to look for it.

Nothing in this argument is nuclear. The only physical input beyond bookkeeping is \(E=mc^{2}\) (Equation (40.6)), and the experiment measures the two sides of Equation (117.21) independently — the left from a mass spectrograph, the right from an absorption curve.

Phenomenon 117.5 (The reaction proceeds far below the Coulomb barrier).

Disintegration is observed with protons accelerated through as little as about \(125\,\mathrm{kV}\) [Cockcroft:1932a], an order of magnitude below the electrostatic barrier that a proton must surmount to reach a lithium nucleus. A classical particle of that energy never arrives; the observed yield is the transmission of the barrier by tunnelling [Gamow:1928] [Gurney:1928]. Rests on Theorem 9.76, Equation (9.93) and Equation (107.2).

Derivation. Derives Phenomenon 117.5. The barrier. Two nuclei of charges \(Z_{1}e\) and \(Z_{2}e\) touch when their centres are separated by \(R\approx r_{0}\left(A_{1}^{1/3}+A_{2}^{1/3}\right)\) with \(r_{0}\approx1.4\,\mathrm{fm}\), and up to that separation the only force acting is the Coulomb repulsion, so the barrier height is

\begin{equation}\tag{117.25} B=\frac{Z_{1}Z_{2}e^{2}}{4\pi\varepsilon_{0}R}\ep \end{equation}

With \(e^{2}/(4\pi\varepsilon_{0})=1.4400\,\mathrm{MeV}\,\mathrm{fm}\) from Equation (107.2), and for a proton on lithium, \(Z_{1}Z_{2}=3\) and \(R=1.4\,\mathrm{fm}\times(1+7^{1/3})=4.08\,\mathrm{fm}\), giving

\begin{equation}\tag{117.26} B=1.059\,\mathrm{MeV}=1.697\times 10^{-13}\,\mathrm{J}\ep \end{equation}

The classical statement. The relevant energy is the one available in the centre of mass, not the laboratory energy: with a proton of \(125\,\mathrm{keV}\) on a lithium nucleus,

\begin{equation}\tag{117.27} E=E_{\text{lab}}\,\frac{m_{\text{Li}}}{m_{\text{Li}}+m_{p}} =109.3\,\mathrm{keV}=1.751\times 10^{-14}\,\mathrm{J}\ec \end{equation}

a factor \(9.7\) below Equation (117.26). A classical proton of that energy turns around at \(b=Z_{1}Z_{2}e^{2}/(4\pi\varepsilon_{0}E)=39.5\,\mathrm{fm}\), nearly ten times the radius at which the nuclear force acts, and the yield would be exactly zero rather than merely small. There is no classical parameter to adjust: the observation is not of a rate lower than expected but of a rate where none was possible.

The quantum statement. The transmission of a Coulomb barrier between a contact radius \(R\) and a turning point \(b\) is a solved problem in this book and is not rebuilt here. Theorem 9.76 gives the semiclassical transmission as \(T\approx\ee^{-2\gamma}\) with \(\gamma=\hbar^{-1}\int\sqrt{2\mu(V-E)} \,\dd r\) taken between the turning points (Equation (9.87)), and Example 9.81 evaluates that integral for exactly this potential, in SI, obtaining

\begin{equation}\tag{117.28} \gamma=\frac{\sqrt{2\mu E}\;b}{\hbar} \left[\arccos\sqrt{\xi}-\sqrt{\xi\left(1-\xi\right)}\right]\ec\qquad \xi=\frac{R}{b}\ec \end{equation}

which is Equation (9.93). Only three inputs are needed to use it. The reduced mass of the proton–lithium system is

\begin{equation}\tag{117.29} \mu=\frac{m_{p}m_{\text{Li}}}{m_{p}+m_{\text{Li}}} =1.4626\times 10^{-27}\,\mathrm{kg}\ec\qquad \mu c^{2}=820.5\,\mathrm{MeV}\ec \end{equation}

built from the nuclear masses of [Wang:2021] with the electrons removed, since it is the nuclei that move. The turning point \(b\) and the ratio \(\xi\) follow from \(E\) as above. And \(\hbar\) is \(1.054572\times 10^{-34}\,\mathrm{J}\,\mathrm{s}\), which is what makes \(\gamma\) dimensionless: the numerator of Equation (117.28) is a momentum times a length, in \(\mathrm{kg}\,\mathrm{m}/\mathrm{s}\times\mathrm{m}\), and dividing by \(\mathrm{J}\,\mathrm{s}\) leaves a pure number.

At \(E=109.3\,\mathrm{keV}\) the momentum scale is \(\sqrt{2\mu E}=7.157\times 10^{-21}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}\), \(b=39.5\,\mathrm{fm}\), \(\xi=0.1032\), the bracket is \(0.9396\), and

\begin{equation}\tag{117.30} 2\gamma=5.04\ec\qquad T=\ee^{-5.04}=6.5\times 10^{-3}\ep \end{equation}

The whole of Table 117.3 is this calculation repeated.

The shape of the yield curve. \(T\) itself is not what a yield curve measures. A proton entering a thick lithium target loses energy continuously to atomic electrons and has a chance of reacting at every energy from its incident value down to zero, so the yield per incident proton is

\begin{equation}\tag{117.31} Y(E_{0})=n\int_{0}^{E_{0}} \frac{\sigma(E)}{\abs{\dd E/\dd x}}\,\dd E\ec\qquad \sigma(E)=\frac{S(E)}{E}\,\ee^{-2\gamma(E)}\ec \end{equation}

with \(n\) the number density of lithium nuclei and \(S(E)\) the astrophysical \(S\) factor of Equation (107.100), which is slowly varying because it is defined to absorb both the barrier and the geometrical \(1/E\). Evaluating Equation (117.31) with \(S\) constant and \(\abs{\dd E/\dd x}\) taken as constant over this narrow range gives the relative yields \(1.00\), \(2.24\), \(7.02\), \(15.4\), \(27.5\), \(62.4\), \(109\) and \(225\) at the eight accelerating potentials of Table 117.3, which run from \(125\,\mathrm{kV}\) to \(700\,\mathrm{kV}\); they are its last column. The measured curve does rise by more than two orders of magnitude across that range, and the reason it rises faster than \(T\) does is the integral, not the barrier.

Remark 117.6 (What the barrier factor does and does not explain).

The transmission in Table 117.3 is of order \(10^{-2}\), while the observed yield per incident proton is of order \(10^{-8}\). The six orders of magnitude between them are not a defect of the calculation and are not tunnelling. A proton of \(125\,\mathrm{keV}\) in lithium metal has a range of order \(10^{-6}\,\mathrm{m}\) and makes of order \(10^{5}\) collisions with atomic electrons before stopping, in the course of which it passes within nuclear distances of essentially no nucleus at all; the geometrical probability of an encounter, and the nuclear matrix element once an encounter occurs, supply the rest of the suppression. What the barrier factor explains is the shape of the curve and the fact that there is any yield whatever at \(125\,\mathrm{kV}\). Reading it as an absolute prediction would be exactly the error Remark 9.77 warns against for the semiclassical method in general: the exponent is the reliable output, the prefactor is not.

Remark 117.7 (The exponent is small here, and that is unusual).

Barrier problems in this book are usually deep in the semiclassical regime: the alpha decay of Example 9.81 has \(2\gamma=67\) and a transmission of \(7\times 10^{-30}\). Here \(2\gamma\) runs from \(5.0\) down to \(0.5\), so the thick-barrier form Equation (9.87) is being used near the edge of its validity, and the two turning points are not far apart on the Airy scale at the upper end of the table. Neither refinement of the prefactor is exact — Remark 9.77 is explicit that Equation (9.86) is not, and offers the uniform form \(T=\left(1+\ee^{2\gamma}\right)^{-1}\) as the alternative — but both move in the same direction, downwards, and by how much is the point. At the lowest row, \(2\gamma=5.04\): \(\ee^{-2\gamma}=6.47\times 10^{-3}\), Equation (9.86) gives \(6.45\times 10^{-3}\) and the uniform form \(6.43\times 10^{-3}\), so the correction is below \(1\,\mathrm{\%}\) where the physical argument is made. At the last row, \(2\gamma=0.48\), the same three read \(0.617\), \(0.463\) and \(0.381\): a correction of \(25\,\mathrm{\%}\) and \(38\,\mathrm{\%}\) respectively, in a row where the proton is nearly over the barrier and the notion of tunnelling has stopped being the right description. The rows that carry the physical argument are the first three, and there the approximation is sound.

Primary references

Three papers report this work and they should be read in order. The Nature letter [Cockcroft:1932a] is the announcement, four weeks after the observation, and it is where the \(125\,\mathrm{kV}\) figure and the identification of the products as alpha particles first appear. The first Proceedings paper [Cockcroft:1932b] is about the machine and not about the physics: it is the reference for the cascade generator, the accelerating tube and the beam, and it is the document that made the technique reproducible elsewhere. The second [Cockcroft:1932c] is the full report of the disintegration, with the yield curve, the absorption measurements and the coincidence test.

The tunnelling interpretation rests on Gamow [Gamow:1928] and, independently, Gurney and Condon [Gurney:1928] [Gurney:1929], all published four years earlier for alpha decay; the Cockcroft–Walton experiment was in fact undertaken because Gamow's calculation, applied in reverse, showed that a few hundred kilovolts might suffice. That is worth noticing: this is one of the few experiments of the period designed from a quantum-mechanical prediction rather than interpreted by one afterwards.

Lawrence and Livingston: the cyclotron (1932)

Tests Phenomenon 117.8.

Apparatus

A flat evacuated brass chamber, about \(0.28\,\mathrm{m}\) across — eleven inches, in the units of the paper — placed between the circular pole faces of an electromagnet whose field is perpendicular to its plane and uniform to a fraction of a per cent over the working region. Inside it, two hollow D-shaped electrodes, the “dees”, face each other across a narrow diametral gap and are connected to a radiofrequency oscillator. The interior of each dee is a field-free region, being a closed conductor: the whole accelerating action takes place in the gap, and the dees exist to shield the ion from the alternating field during the half revolution in which that field would decelerate it.

An ion source sits at the centre, a filament and a gas feed producing hydrogen ions in the middle of the gap. A deflecting electrode at the periphery, held at a few kilovolts, pulls the beam out of its last orbit and onto a collector. The chamber was pumped to a residual pressure of order \(10^{-4}\,\mathrm{Pa}\), which is set not by any technical convenience but by the physics: the ion travels a total path length of order

\begin{equation}\tag{117.32} L=\sum_{k=1}^{N}2\pi r_{k} =2\pi R\sum_{k=1}^{N}\sqrt{\frac{k}{N}} \approx\frac{4\pi}{3}NR \approx73\,\mathrm{m}\ec \end{equation}

inside the machine — using \(r_{k}=R\sqrt{k/N}\), since the energy grows linearly with the crossing count and the radius as its square root — and must not scatter out of resonance on the way. Put the requirement in areal density, which is what stops a particle. A proton of \(1.22\,\mathrm{MeV}\) has a range in air of about \(3.3\times 10^{-2}\,\mathrm{kg}/\mathrm{m}^{2}\), while \(73\,\mathrm{m}\) of air at \(1\,\mathrm{Pa}\) and \(288.15\,\mathrm{K}\) — a density of \(1.21\times 10^{-5}\,\mathrm{kg}/\mathrm{m}^{3}\) by Equation (117.3) scaled with pressure — presents only \(8.8\times 10^{-4}\,\mathrm{kg}/\mathrm{m}^{2}\). So \(1\,\mathrm{Pa}\) does not stop the beam outright; the pressure that would is about \(40\,\mathrm{Pa}\), and the residual pressure has to be far below even that, because a single elastic scatter through a small angle already takes the ion out of resonance. At \(10^{-4}\,\mathrm{Pa}\) the gas-kinetic mean free path \(k_{B}T/(\sqrt{2}\pi d^{2}p)\), with a molecular diameter \(d\approx3\times 10^{-10}\,\mathrm{m}\) and \(T=300\,\mathrm{K}\), is \(10^{2}\,\mathrm{m}\) — longer than the \(73\,\mathrm{m}\) orbit, and for a \(1\,\mathrm{MeV}\) proton, whose scattering cross-section is far below gas-kinetic, longer still. That is the margin the pump had to buy.

The magnet is the expensive component and the one that sets the energy. Everything else in the machine handles a few kilovolts.

Procedure

Ions released at the centre are bent by the magnetic field into a circle and cross the gap twice per revolution. If the oscillator frequency is matched to the revolution frequency, each crossing finds the gap field reversed and therefore accelerating, so the ion gains \(qV_{\text{gap}}\) at every crossing while its orbit grows as the square root of its energy.

The experiment consists of three measurements. First, tuning: the magnetic field is swept at fixed oscillator frequency and the collected current recorded, which produces a resonance curve whose peak locates the condition \(\omega_{c}=\omega_{\text{osc}}\) and whose width measures how far off resonance the machine will still accelerate. Second, the current at resonance, which measures how much beam survives the \(73\,\mathrm{m}\) of orbit of Equation (117.32). Third, the energy of the extracted beam, measured by absorption in foils against the known range–energy relation for protons — the same technique as Section 117.3, and the same one that provides the only independent check that the ions really did make a hundred and fifty turns.

The resonance curve is the crucial one, because it is what distinguishes this machine from a peculiar high-voltage device. A current that appears only in a narrow band of magnetic field, centred on the value that makes Equation (117.34) equal the oscillator frequency, cannot be explained by any static potential in the apparatus.

Observations and data

Protons of kinetic energy above \(1.2\,\mathrm{MeV}\) were obtained from a peak voltage on the dees of a few kilovolts, with beam currents of order \(10^{-9}\,\mathrm{A}\) at the collector. The ratio of the two voltages is the whole point: at about \(4\,\mathrm{kV}\) per gap crossing and two crossings per turn, the final energy corresponds to something like a hundred and fifty revolutions, and the accelerating voltage that has to be insulated anywhere in the machine is smaller than the beam energy in electronvolts by more than two orders of magnitude.

Table 117.4 collects the measured quantities and the machine parameters that follow from them. The field and frequency are not independent measurements but consequences of the resonance condition and the extraction radius, and are computed here; they are included because they are what the resonance curve actually determines.

QuantityValueKind, source
Chamber diameter\(0.28(1)\,\mathrm{m}\)obs. 1932, [Lawrence:1932]
Extraction radius\(0.114(5)\,\mathrm{m}\)obs. 1932, [Lawrence:1932]
Proton kinetic energy\(1.22(6)\,\mathrm{MeV}\), i.e. \(1.955(96)\times 10^{-13}\,\mathrm{J}\)obs. 1932, [Lawrence:1932]
Beam current at the collectororder \(10^{-9}\,\mathrm{A}\)obs. 1932, [Lawrence:1932]
Peak dee voltage\(4.0(5)\,\mathrm{kV}\)obs. 1932, [Lawrence:1932]
Proton momentum at extraction\(2.557(63)\times 10^{-20}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}\)calc.
Magnetic field required\(1.40(7)\,\mathrm{T}\)calc., Equation (117.33)
Revolution frequency\(21.3(11)\,\mathrm{MHz}\)calc., Equation (117.34)
Number of revolutions\(153(20)\)calc.
Total orbit length\(73(10)\,\mathrm{m}\)calc., Equation (117.32)
Relativistic frequency shift at extraction$1.30(6)\times 10^{-3}$calc.
Accumulated phase slip over the acceleration\(0.62(9)\,\mathrm{rad}\), i.e. \(36(5)\,^\circ\)calc., Equation (117.38)
The 1932 cyclotron: what was measured and what follows. The proton energy and the beam current are the reported observations; the magnetic field, the revolution frequency and the turn count are computed here from Equation (117.36) and Equation (117.34) with the quoted chamber size, and are the parameters the resonance condition fixes. The last two rows are the relativistic limit of the design, and are the reason no later machine was built this way. Bracketed uncertainties on the observed rows are estimates from the stated technique, the paper quoting none: \(5\,\mathrm{\%}\) on an energy read from a foil absorption curve, \(4\,\mathrm{\%}\) on a radius set by the position of the deflector, and an order of magnitude on the collected current. The computed rows inherit them — the momentum to \(2.5\,\mathrm{\%}\), the field and the frequency to \(5\,\mathrm{\%}\), the turn count, the orbit length and the phase slip to about \(14\,\mathrm{\%}\) — and are quoted with the propagated figure in brackets.

Interpretation

The cyclotron removed the ceiling that high-voltage technique imposes. Cockcroft and Walton's method needs the full energy to exist as a potential difference somewhere in the laboratory, and the dielectric strength of air sets a limit of a few million volts; the cyclotron needs only the gap voltage, and buys energy with turns instead of volts. The exchange is not free — it is paid for in orbit length, in vacuum, and in magnetic-field uniformity — but all three are technical problems with technical solutions, whereas insulation is a material property.

Everything downstream in this treatise depends on that trade: the beams of Experiment: Deep Inelastic Scattering, the colliding beams of Experiment: The Higgs Boson Discovery, and the reactor and accelerator neutrino sources of Experiment: Neutrino Oscillations are all descendants of it. So is a large part of nuclear astrophysics, since the reaction cross-sections that Stellar Structure and Nucleosynthesis needs — the cross-sections whose Gamow-peak energies are computed in Theorem 107.85 — are measured on machines of this lineage, at energies where the direct astrophysical rate is unmeasurably small.

The classical cyclotron has one limit, and it is relativistic. As the ion's energy grows its revolution frequency falls, the resonance detunes, and acceleration stops. The cure is to modulate the frequency (the synchrocyclotron) or to raise the field in step with the energy at fixed radius (the synchrotron); both are corrections demanded by Relativistic Dynamics and by nothing else, and they are the reason no accelerator built after the 1940s is a cyclotron in the original sense. The 1932 machine sits exactly at the boundary, which is what makes it worth computing the phase slip below rather than merely asserting the limitation.

Phenomenon 117.8 (Cyclotron resonance).

A non-relativistic ion of charge \(q\) and mass \(m\) moving in a uniform magnetic field of magnitude \(B\) circulates with an angular frequency that does not depend on its speed or on the radius of its orbit. A fixed-frequency alternating voltage across a single gap therefore stays in step with the ion for hundreds of revolutions and accelerates it at every crossing, so that the final energy is set by the magnet and not by the applied voltage [Lawrence:1932].

Derivation. Derives Phenomenon 117.8. For circular motion in a uniform field perpendicular to the orbit, the magnetic force supplies the centripetal acceleration,

\begin{equation}\tag{117.33} qvB=\frac{mv^{2}}{r} \quad\Longrightarrow\quad r=\frac{mv}{qB}=\frac{p}{qB}\ec \end{equation}

in SI a radius in metres when \(p\) is in \(\mathrm{kg}\,\mathrm{m}/\mathrm{s}\), \(q\) in coulombs and \(B\) in tesla. The angular frequency is therefore

\begin{equation}\tag{117.34} \omega_{c}=\frac{v}{r}=\frac{qB}{m}\ec \end{equation}

independent of \(v\) and of \(r\): the ion takes the same time to complete a small slow circle as a large fast one, which is exactly the property that lets one fixed frequency serve the whole acceleration. It is a peculiarity of the Lorentz force and of no other force law that this is so; a linear restoring force would give the same property, and an inverse-square attraction would not.

Numbers. At the extraction radius \(R=0.114\,\mathrm{m}\) a proton of \(1.22\,\mathrm{MeV}\) has momentum \(p=\sqrt{2m_{p}T}=2.557\times 10^{-20}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}\), using \(m_{p}=1.6726\times 10^{-27}\,\mathrm{kg}\) [Mohr:2025], so Equation (117.33) requires

\begin{equation}\tag{117.35} B=\frac{p}{eR} =\frac{2.557\times 10^{-20}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}} {1.602177\times 10^{-19}\,\mathrm{C}\times0.114\,\mathrm{m}} =1.400\,\mathrm{T}\ec \end{equation}

an ordinary iron-cored electromagnet, and Equation (117.34) then gives \(\omega_{c}/2\pi=eB/2\pi m_{p}=21.3\,\mathrm{MHz}\), a free-space wavelength of \(14\,\mathrm{m}\) and an entirely standard radiofrequency in 1932.

Energy from the magnet, not the voltage. Eliminating \(v\) between Equation (117.33) and the kinetic energy gives the extraction energy at the outer radius \(R\),

\begin{equation}\tag{117.36} T=\tfrac{1}{2}mv^{2}=\frac{q^{2}B^{2}R^{2}}{2m}\ec \end{equation}

which contains the gap voltage nowhere. Checking: \(e^{2}B^{2}R^{2}/2m_{p}=1.955\times 10^{-13}\,\mathrm{J}=1.22\,\mathrm{MeV}\), as observed. The number of crossings needed is \(N_{\text{cross}}=T/qV_{\text{gap}}=1.22\,\mathrm{MeV} /4\,\mathrm{keV}=305\), that is \(153\) revolutions, and the gap voltage merely fixes how long the ion must stay in the machine. Note the scaling: \(T\propto B^{2}R^{2}\), so the energy is quadratic in the size of the magnet, which is why the next machines were bigger rather than cleverer.

The limit. Relativistically the momentum is \(\gamma mv\) with \(\gamma=\left(1-v^{2}/c^{2}\right)^{-1/2}\), so Equation (117.34) becomes \(\omega=qB/\gamma m\): the frequency falls as the ion becomes energetic, and the resonance that Equation (117.34) guarantees is lost. For protons at \(1.22\,\mathrm{MeV}\) the shift is

\begin{equation}\tag{117.37} \gamma-1=\frac{T}{m_{p}c^{2}} =\frac{1.22\,\mathrm{MeV}}{938.27\,\mathrm{MeV}} =1.30\times 10^{-3}\ec \end{equation}

using the proton mass energy of [Navas:2024]. The consequence is cumulative and that is what has to be estimated. Suppose the energy grows linearly with turn number \(n\), so that \(\gamma_{n}-1\approx(n/N)(\gamma_{N}-1)\); each turn the ion arrives at the gap late by a fraction \(\gamma_{n}-1\) of a period, and after \(N\) turns the accumulated phase error is

\begin{equation}\tag{117.38} \Delta\phi=2\pi\sum_{n=1}^{N}\left(\gamma_{n}-1\right) \approx\pi N\left(\gamma_{N}-1\right) =\pi\times153\times1.30\times 10^{-3} =0.62\,\mathrm{rad}\ec \end{equation}

that is \(36\,^\circ\). Acceleration continues as long as the ion arrives within a quarter period of the correct phase, so \(36\,^\circ\) is tolerable and the machine works. At ten times the energy the same estimate gives \(\Delta\phi\approx62\,\mathrm{rad}\), since both \(N\) and \(\gamma_{N}-1\) grow with \(T\) and the slip therefore grows as \(T^{2}\); the ion would be decelerated for most of its life and never reach the rim. The classical cyclotron works precisely in the regime where Relativistic Dynamics is not yet needed, and it stops working exactly where it is.

Remark 117.9 (Why the resonance is not fragile).

It is easy to read Equation (117.38) as saying that the machine requires the field and the frequency to be matched to a part in \(10^{3}\), and that is not quite right. What matters is the accumulated phase, so a static mismatch \(\Delta\omega/\omega\) produces a slip \(2\pi N\Delta\omega/\omega\), and tolerating a quarter period over \(153\) turns means \(\abs{\Delta\omega/\omega}<1.6\times 10^{-3}\) — which is a demanding but attainable tolerance on an iron magnet, and is precisely the width of the resonance curve that the experiment measured. The measured width is therefore not an instrumental artefact but a direct reading of the turn count, and it is the strongest evidence in the paper that the ions really were accelerated many times rather than once.

Primary references

The single reference is Lawrence and Livingston's Physical Review paper [Lawrence:1932], whose title — the production of high-speed light ions without the use of high voltages — states the whole result. It should be read alongside [Cockcroft:1932b], published the same year, since the two describe opposite solutions to the same problem and each is explicit about the limitation of the other's approach. Neither paper reports a nuclear reaction; the cyclotron's first disintegrations followed within months, using the lithium target of Section 117.3 as the obvious first check that the beam was real.

Wilson: the expansion cloud chamber (1912)

Tests Phenomenon 117.10.

Apparatus

A shallow cylindrical chamber with glass walls and roof, of order \(0.16\,\mathrm{m}\) in diameter and a few centimetres deep, containing dust-free air saturated with water vapour. Its floor is a piston whose sudden withdrawal expands the gas; the expansion is completed in a time of order \(10^{-2}\,\mathrm{s}\), short compared with the time in which heat can reach the gas from the walls, so it is adiabatic and therefore cools the gas, leaving the vapour supersaturated.

Three further elements make it an instrument rather than a demonstration. The expansion is stopped mechanically at a preset ratio, because the ratio is the control variable and must be reproducible to better than a per cent. The chamber is illuminated from the side by a bright flash — the illumination must be transverse, since the droplets are seen by scattered light against a dark background — and photographed, in the developed form stereoscopically, so that the track can be reconstructed in three dimensions. And a clearing field of order \(10^{3}\,\mathrm{V}/\mathrm{m}\) between electrodes above and below the gas sweeps out the ions left by previous events between expansions, without which the chamber accumulates a fog of old tracks.

Dust is the enemy and its removal is the reason the instrument took fifteen years to develop. Every dust grain is a condensation nucleus, and a chamber with dust in it fogs at any expansion. Wilson removed it by repeated expansions with the condensate allowed to settle, each cycle sweeping the remaining particles down with the droplets they nucleated.

Procedure

The expansion ratio is the control. Below a threshold no condensation occurs at all; above a first threshold droplets form on ions; above a second, higher one the vapour condenses on everything and the chamber fogs uniformly. Working between the two thresholds makes ions, and only ions, into condensation nuclei — and the existence of a window between the thresholds, rather than a single onset, is the physical fact on which the instrument rests. It is derived below.

The operating cycle is then fixed by the physics of droplet growth. The chamber is expanded in synchronism with the passage of the radiation under study — in the developed form, triggered by a counter telescope that has already registered the particle. The illumination is flashed a few hundredths of a second later: long enough for the droplets to grow by condensation to a size that scatters visible light efficiently, which requires a radius of order \(10^{-6}\,\mathrm{m}\), and short enough that diffusion has not yet blurred the trail. A water molecule in air at \(270\,\mathrm{K}\) has a diffusion coefficient of order \(2\times 10^{-5}\,\mathrm{m}^{2}/\mathrm{s}\), so in \(2\times 10^{-2}\,\mathrm{s}\) an ion wanders \(\sqrt{2Dt}\approx9\times 10^{-4}\,\mathrm{m}\) — under a millimetre, which is why the track is still sharp at that delay and would not be at one second.

After the photograph the chamber is recompressed, the clearing field applied, and the cycle repeated, at a rate of a few per minute. That duty cycle — seconds of dead time for each millisecond of sensitivity — is the instrument's defining weakness and is what Section 117.7 eventually removed.

Observations and data

The photographs show the paths of individual ionising particles as sharp lines of droplets. Alpha particles give dense, straight, sharply terminated trails, occasionally with a sudden deflection near the end where a nuclear encounter has occurred; fast beta particles give thin, tortuous trails; X-rays give no trail at all but a scatter of short electron tracks starting at points along the beam.

That contrast is quantitative and is the first thing the instrument measures. The specific ionisation of a charged particle behaves as Equation (115.15), so a slow or multiply charged particle ionises far more densely than a fast singly charged one. For an alpha particle of \(7.687\,\mathrm{MeV}\) depositing its whole energy in \(7.0\times 10^{-2}\,\mathrm{m}\) of air, the mean deposit is \(1.10\times 10^{8}\,\mathrm{eV}/\mathrm{m}\), that is \(1.76\times 10^{-11}\,\mathrm{J}/\mathrm{m}\); at a mean expenditure of about \(34\,\mathrm{eV}\) — \(5.4\times 10^{-18}\,\mathrm{J}\) — per ion pair in air [Navas:2024], that is roughly \(3.2\times 10^{6}\,/\mathrm{m}\) of track. A minimum-ionising electron in the same gas loses about \(2.2\times 10^{5}\,\mathrm{eV}/\mathrm{m}\), that is \(3.5\times 10^{-14}\,\mathrm{J}/\mathrm{m}\), and makes about \(6.5\times 10^{3}\,/\mathrm{m}\). The ratio, some five hundred, is exactly the difference between a dense unbroken trail and a sparse dotted one, and it is why the two species are distinguishable by eye in a photograph.

The thresholds themselves are the chamber's calibration, and are collected in Table 117.5. The expansion ratios are Wilson's; the temperatures and supersaturations are computed here from them by Equations (117.39) and (117.40), with no free parameter — the initial temperature is the only input and the saturation vapour pressure cancels.

$r=V_{2}/V_{1}$$T_{2}$$S$What is observed, and what is predicted
\(1.250(5)\)\(268.1(4)\,\mathrm{K}\)\(4.11(11)\)droplets on negative ions; charged-nucleus threshold predicted at $S=4.47$
\(1.310(5)\)\(263.1(4)\,\mathrm{K}\)\(5.64(15)\)droplets on positive ions as well; same prediction, and the sign asymmetry is outside the model
\(1.380(5)\)\(257.7(4)\,\mathrm{K}\)\(8.08(23)\)dense fog on no nucleus at all; neutral barrier here is $34\,k_{B}T$, the onset of macroscopic homogeneous nucleation
The condensation window of a Wilson chamber. The expansion ratio $r=V_{2}/V_{1}$ is the experimental control [Wilson:1912]; $T_{2}$ and the supersaturation $S$ are computed here from Equations (117.39) and (117.40) for air initially saturated at $T_{1}=293.15\,\mathrm{K}$, using $\gamma=1.40$ and the Clausius–Clapeyron relation with a latent heat of \(2.50\times 10^{6}\,\mathrm{J}/\mathrm{kg}\). The last column is the critical supersaturation predicted by Equation (117.46) for a singly charged nucleus, or by the neutral nucleation barrier Equation (117.43) where marked; the agreement between the observed thresholds and the predictions is the content of Phenomenon 117.10. The expansion ratio is set mechanically and is reproducible to better than \(0.5\,\mathrm{\%}\), which is the bracketed uncertainty; $T_{2}$ and $S$ carry what that propagates to through Equations (117.39) and (117.40). The barrier in the last row is quoted with the constant $\sigma=0.0765\,\mathrm{N}/\mathrm{m}$ used throughout the derivation, and would fall to $36\,k_{B}T$ with the surface tension of supercooled water at \(257.7\,\mathrm{K}\); the difference is the honest size of the model's grip on the third figure.

The tabulated droplet densities of the original paper are not transcribed here. What is reproduced is the calibration that makes the instrument an instrument: the two thresholds, their ratio, and the prediction of both from the thermodynamics of a charged drop.

Interpretation

Before this apparatus, radiation was studied through the aggregate effects of very many particles — an ionisation current, a photographic blackening, a scintillation rate. The cloud chamber makes a single particle individually visible and individually measurable, and with a magnetic field applied it measures momentum and sign of charge from the curvature, through \(p_{\perp}=ZeB\rho\) (Equation (115.14)).

Everything that follows in this book from the inspection of individual events begins here: the two-branch forks that identified Equation (117.5) [Blackett:1925], the electron–positron pairs photographed by Blackett and Occhialini [Blackett:1933] and treated in Experiment: The Positron, and the cosmic-ray showers and strange particles of Cosmic Rays and Astroparticle Physics. In each case the decisive information is not a rate but a picture: how many branches, which way they curve, whether the momenta balance.

Its weaknesses defined the next two instruments. The gas is thin — a cubic metre of air holds \(1.2\,\mathrm{kg}\) of matter and a cubic metre of liquid holds a tonne — so interactions are rare and the chamber is a poor target; and the vapour must be re-expanded and allowed to settle, so the chamber is dead for most of every minute. Section 117.6 fixes the first, Section 117.7 the second.

Phenomenon 117.10 (Ion tracks are made visible by condensation).

In a gas supersaturated with vapour, droplets condense preferentially on the ions left by the passage of a charged particle, and on nothing else, provided the supersaturation lies between the threshold for condensation on ions and the threshold for condensation on neutral molecules. The trail of droplets marks the particle's path and can be photographed [Wilson:1912]. Rests on Equation (115.15).

Derivation. Derives Phenomenon 117.10. The claim has two halves — that an expansion produces a calculable supersaturation, and that a charge lowers the supersaturation at which a droplet can grow — and the instrument exists because the two thresholds so produced are different.

The supersaturation from the expansion. Air is expanded from volume \(V_{1}\) to \(V_{2}\) fast enough to be adiabatic and slowly enough to be reversible, so \(pV^{\gamma}\) is constant with \(\gamma=c_{p}/c_{V}=1.40\) for a diatomic gas, and \(TV^{\gamma-1}\) is constant. Writing \(r=V_{2}/V_{1}\),

\begin{equation}\tag{117.39} T_{2}=T_{1}\,r^{1-\gamma}\ep \end{equation}

The water vapour is a minor component and is carried along: its partial pressure falls both because the volume grew and because the temperature fell, by the ideal-gas law at fixed particle number, \(p_{v2}=p_{v1}\left(V_{1}/V_{2}\right)\left(T_{2}/T_{1}\right) =p_{v1}r^{-\gamma}\). The saturation vapour pressure falls far faster, by the Clausius–Clapeyron relation integrated with a constant latent heat \(L\) of vaporisation per unit mass and a specific gas constant \(R_{v}=R/M_{\mathrm{H_{2}O}} =461.5\,\mathrm{J}/\mathrm{kg}/\mathrm{K}\): \(p_{\text{sat}}(T_{2})=p_{\text{sat}}(T_{1}) \exp\left[-\left(L/R_{v}\right)\left(1/T_{2}-1/T_{1}\right)\right]\). Since the chamber starts saturated, \(p_{v1}=p_{\text{sat}}(T_{1})\) and that quantity cancels entirely: the supersaturation \(S=p_{v2}/p_{\text{sat}}(T_{2})\) is

\begin{equation}\tag{117.40} S=r^{-\gamma} \exp\left[\frac{L}{R_{v}} \left(\frac{1}{T_{2}}-\frac{1}{T_{1}}\right)\right]\ec \end{equation}

a function of the expansion ratio and the initial temperature alone. It requires no vapour-pressure table, which is what makes it a prediction rather than a fit; the one measured property it does need is the latent heat, and that enters only through the exponential. With \(L=2.50\times 10^{6}\,\mathrm{J}/\mathrm{kg}\) and \(T_{1}=293.15\,\mathrm{K}\), Equation (117.39) gives \(T_{2}=268.1\,\mathrm{K}\) at \(r=1.25\) and Equation (117.40) gives \(S=4.11\); the three rows of Table 117.5 are this calculation repeated.

Why a charge helps. Consider forming a spherical droplet of radius \(\rho\) from vapour at supersaturation \(S\), the droplet carrying a charge \(q\) at its centre. Three terms contribute to the free energy of formation. The surface costs \(4\pi\rho^{2}\sigma\) with \(\sigma\) the surface tension. The bulk gains, because the vapour is supersaturated: each molecule transferred from vapour to liquid lowers the free energy by \(k_{B}T\ln S\), and a droplet of radius \(\rho\) holds \(\tfrac{4}{3}\pi\rho^{3}n_{l}\) of them, with the number density of liquid water — mass density over molar mass, times the Avogadro constant, so that it is not an independent input either —

\begin{equation}\tag{117.41} n_{l}=\frac{1.000\times 10^{3}\,\mathrm{kg}/\mathrm{m}^{3} \times6.02214076\times 10^{23}\,/\mathrm{mol}} {18.0153\times 10^{-3}\,\mathrm{kg}/\mathrm{mol}} =3.343\times 10^{28}\,/\mathrm{m}^{3}\ep \end{equation}

And the electrostatic self-energy of the charge falls as the droplet grows around it, because the charge's field is being moved out of vacuum and into a medium of relative permittivity \(\varepsilon_{r}\). The geometry matters and is easy to get backwards: the droplet is the dielectric, so it is the region \(a<r<\rho\) between the bare ion of radius \(a\) and the droplet surface whose permittivity changes, while the field beyond \(\rho\) is unaltered. The energy stored there falls from \(\left(q^{2}/8\pi\varepsilon_{0}\right)\left(1/a-1/\rho\right)\) to \(\left(q^{2}/8\pi\varepsilon_{0}\varepsilon_{r}\right) \left(1/a-1/\rho\right)\), so relative to the bare ion the change is \(-\left(q^{2}/8\pi\varepsilon_{0}\right) \left(1-1/\varepsilon_{r}\right)\left(1/a-1/\rho\right)\), whose \(\rho\)-independent piece is a constant of the ion and drops out. Collecting the terms that depend on \(\rho\),

\begin{equation}\tag{117.42} \Delta G(\rho)=4\pi\rho^{2}\sigma -\tfrac{4}{3}\pi\rho^{3}n_{l}k_{B}T\ln S +\frac{q^{2}}{8\pi\varepsilon_{0}\rho} \left(1-\frac{1}{\varepsilon_{r}}\right)\ec \end{equation}

which is Thomson's expression [Thomson:1906], in SI, with every term in joules. Setting \(q=0\) recovers the classical nucleation problem: the surface term wins at small \(\rho\) and the bulk term at large, so there is a barrier at \(\rho_{c}=2\sigma/\left(n_{l}k_{B}T\ln S\right)\) of height

\begin{equation}\tag{117.43} W_{c}=\frac{16\pi\sigma^{3}} {3\left(n_{l}k_{B}T\ln S\right)^{2}}\ec \end{equation}

and nothing condenses unless a thermal fluctuation supplies \(W_{c}\).

The two thresholds. With \(q\neq0\) the third term dominates at small \(\rho\), so \(\Delta G\) rises to \(+\infty\) as the droplet shrinks back onto the bare ion — but it falls as \(\rho\) grows, with a downhill slope \(q^{2}\left(1-1/\varepsilon_{r}\right)/8\pi\varepsilon_{0}\rho^{2}\) that diverges as \(\rho\to0\), while the surface term's uphill slope \(8\pi\rho\sigma\) vanishes there. Near the ion the free energy is therefore downhill throughout, and the barrier, if one survives at all, is pushed out to larger \(\rho\) and lowered. Above some supersaturation it disappears altogether, and a droplet grows on the ion without any fluctuation at all. Setting \(\dd\Delta G/\dd\rho=0\) and asking when the resulting equation has no root gives the condition directly. Dividing \(\dd\Delta G/\dd\rho=0\) through by \(4\pi\rho^{2}\), the stationarity condition reads

\begin{equation}\tag{117.44} n_{l}k_{B}T\ln S=\frac{2\sigma}{\rho} -\frac{q^{2}\left(1-1/\varepsilon_{r}\right)} {32\pi^{2}\varepsilon_{0}\rho^{4}}\ep \end{equation}

The right-hand side rises from \(-\infty\), peaks, and falls as \(2\sigma/\rho\); its maximum is at

\begin{equation}\tag{117.45} \rho_{*}=\left[\frac{q^{2}\left(1-1/\varepsilon_{r}\right)} {16\pi^{2}\varepsilon_{0}\sigma}\right]^{1/3}\ec \end{equation}

where the second term is one quarter of the first, so the peak value is \(\tfrac{3}{4}\times2\sigma/\rho_{*}\). If the left-hand side exceeds that peak there is no stationary point at all, \(\Delta G\) decreases monotonically, and the droplet grows freely. The threshold is therefore

\begin{equation}\tag{117.46} \ln S_{\text{crit}}=\frac{3\sigma} {2\,n_{l}k_{B}T\,\rho_{*}}\ep \end{equation}

Numbers. At \(T=268.1\,\mathrm{K}\), with \(\sigma=0.0765\,\mathrm{N}/\mathrm{m}\), \(\varepsilon_{r}=88\) for water at that temperature, and \(q=e=1.602177\times 10^{-19}\,\mathrm{C}\), Equation (117.45) gives \(\rho_{*}=6.19\times 10^{-10}\,\mathrm{m}\) — about \(33\) water molecules — and Equation (117.46) gives \(S_{\text{crit}}=4.47\). The observed threshold for condensation on ions is \(S=4.11\), from the first row of Table 117.5. Meanwhile at that same supersaturation the neutral barrier Equation (117.43) is \(2.45\times 10^{-19}\,\mathrm{J}=66\,k_{B}T\), so the probability of a spontaneous fluctuation forming a droplet is of order \(\ee^{-66}\) and nothing whatever condenses on neutral molecules. At the third row, \(S=8.08\) and \(T=257.7\,\mathrm{K}\), the same barrier has fallen to \(1.21\times 10^{-19}\,\mathrm{J}=34\,k_{B}T\). The prefactor here is the number density of vapour molecules and not the liquid density \(n_{l}\): the partial pressure of the water vapour after the expansion is of order \(10^{3}\,\mathrm{Pa}\), so \(p/k_{B}T\) is of order \(10^{23}\,/\mathrm{m}^{3}\). With molecular attempt frequencies of order \(10^{12}\,/\mathrm{s}\) and a Boltzmann factor \(\ee^{-34}=1.7\times 10^{-15}\), the nucleation rate is of order \(10^{20}\,/\mathrm{m}^{3}/\mathrm{s}\) — that is \(10^{14}\) droplets per cubic centimetre per second, a dense fog immediately. That is the second threshold. The conclusion is insensitive to the prefactor by a wide margin: getting it wrong by five orders of magnitude, which is the size of the gap between the vapour and liquid densities, is worth a change of \(12\) in \(W_{c}/k_{B}T\) — and \(W_{c}/k_{B}T\) moves from \(66\) to \(34\) across the two rows.

The window between \(S=4.5\) and \(S=8\) is therefore not an empirical convenience but a consequence of Equation (117.42): the charge removes the barrier while the neutral vapour still has one. Operating inside it makes ions, and only ions, into condensation nuclei, which is the statement of Phenomenon 117.10.

Remark 117.11 (Where this derivation is honest and where it is not).

Two things in the calculation above deserve flagging rather than burying. First, \(\rho_{*}=0.62\,\mathrm{nm}\) contains about thirty molecules, and applying a macroscopic surface tension and a bulk permittivity to an object that size is not defensible in detail; that the answer comes out within \(10\,\mathrm{\%}\) of the measurement is better than the method deserves, and the right claim is that the mechanism and the order of magnitude are established, not the third figure. Second, Equation (117.42) depends on the charge only through \(q^{2}\) and therefore predicts the same threshold for positive and negative ions, whereas Wilson measured \(r=1.25\) for negative and \(r=1.31\) for positive, a difference of \(37\,\mathrm{\%}\) in \(S\). The asymmetry is real and is outside this model entirely: it involves the orientation of the permanent dipole of the water molecule at the droplet surface, which no theory built from a scalar \(\varepsilon_{r}\) can see. Reporting the model as though it accounted for both thresholds would misstate what it does.

Third, the thermophysical inputs are not cited, because this bibliography carries no source for them: the latent heat of vaporisation \(2.50\times 10^{6}\,\mathrm{J}/\mathrm{kg}\), the surface tension \(0.0765\,\mathrm{N}/\mathrm{m}\) and the static permittivity \(\varepsilon_{r}=88\) are standard tabulated properties of water taken at the stated temperature, and the reader should treat them as handbook values rather than as claims this book supports. They matter unequally: \(S\) depends on \(L\) through an exponential, \(\rho_{*}\) on \(\sigma^{-1/3}\) and only logarithmically on \(\varepsilon_{r}\), and \(W_{c}\) on \(\sigma^{3}\). The surface tension is the one that carries the third figure, and it is also the one that is least defensible at this droplet size, as the first point above says. Holding \(\sigma\) fixed at its \(268.1\,\mathrm{K}\) value across the table, rather than letting it run with temperature, is what separates the \(34\,k_{B}T\) quoted in the last row of Table 117.5 from the \(36\,k_{B}T\) the temperature-dependent value would give.

Primary references

Wilson's 1912 paper [Wilson:1912] is the instrument paper and reads like one: most of it is about dust, valves, timing and the photography, and the tracks appear at the end as the demonstration that the apparatus works. The physics of the condensation window is older — it is Wilson's own work of the eighteen-nineties on the two expansion thresholds — and the theory of the charged drop is Thomson's [Thomson:1906], given there in the chapter on the condensation of vapour on ions.

The papers that made the instrument famous are not Wilson's. Blackett's identification of the nitrogen reaction [Blackett:1925], Blackett and Occhialini's photographs of electron–positron pairs [Blackett:1933], and the cosmic-ray work of Cosmic Rays and Astroparticle Physics are all cloud-chamber results, and each is a case where the answer was visible in a photograph and inaccessible to any counting experiment then available.

Glaser: the bubble chamber (1952)

Tests Phenomenon 117.12.

Apparatus

A small, scrupulously clean glass bulb of a few cubic centimetres, filled with diethyl ether and held under pressure at a temperature above its normal boiling point of \(307.7\,\mathrm{K}\). In the working condition the liquid sat near \(413\,\mathrm{K}\), where its saturated vapour pressure is about \(1.48\,\mathrm{MPa}\), roughly fifteen atmospheres; the bulb was held above that pressure so that the liquid stayed liquid. Releasing the pressure to about \(0.10\,\mathrm{MPa}\) leaves the liquid superheated by \(\Delta p\approx1.38\,\mathrm{MPa}\): it is far above its boiling point at the prevailing pressure and yet, in the absence of nucleation sites, does not boil.

A radioactive source could be brought up to the bulb or removed, or interposed behind lead, and the onset of boiling was registered — in the first experiments by timing, later photographically.

Cleanliness is the apparatus. A superheated liquid is metastable against boiling everywhere, and the barrier that protects it, Equation (117.50) below, is enormous only for a homogeneous liquid. A dissolved gas bubble, a dust grain, or a scratch on the glass provides a surface of large radius of curvature at which the barrier is absent, so the vessel boils at once and no measurement is possible. Glass was chosen for its smoothness, the bulbs were sealed after prolonged degassing, and the ether was purified; the experimental difficulty of the whole subject is in that sentence.

Procedure

The bulb is taken through a cycle: compress, heat, expand into the superheated state, wait. The quantity measured is the delay between the expansion and the onset of boiling, recorded many times over with the source present and with it absent, so that a shift in the distribution of delays is a statistical statement rather than an anecdote.

Two features of that protocol matter. The delay is a random variable, because nucleation is a stochastic process — a barrier crossing has an exponentially distributed waiting time — so the observable is a distribution and not a number, and the comparison must be between distributions. And the null run is not “no source” but “source shielded”, which controls for everything about the thermal cycle while changing only the radiation reaching the liquid.

The degree of superheat was varied, since the barrier Equation (117.50) depends on it as \(\left(\Delta p\right)^{-2}\) and therefore very steeply: too little superheat and even a radiation track cannot nucleate, too much and the liquid boils spontaneously and the chamber is useless. Finding the window is the same experimental problem as finding Wilson's condensation window in Section 117.5, and has the same two-sided structure.

Observations and data

Superheated ether that stayed liquid for a long time when shielded boiled promptly when an ionising source was brought near, and the distribution of delays shifted with the source strength. The nucleation is therefore caused by the radiation. Photographs of bubbles strung out along the paths of the ionising particles — the bubble chamber proper — followed shortly after.

Table 117.6 gives the thermodynamic state of the liquid in the working condition, with the critical bubble it implies. The temperature and the pressures are the operating point of the experiment; the surface tension, critical radius and barrier are computed here, and the surface tension in particular is an estimate rather than a measurement, obtained by scaling the room-temperature value with \(\left(1-T/T_{c}\right)^{1.26}\) to the working temperature, which is the corresponding-states form [Guggenheim:1945]. It is flagged as such because the barrier depends on its cube.

QuantityValueKind
Normal boiling point of diethyl ether\(307.7\,\mathrm{K}\)tabulated, uncited
Critical point\(466.7\,\mathrm{K}\), \(3.64\,\mathrm{MPa}\)tabulated, uncited
Working temperature\(413(5)\,\mathrm{K}\)obs. 1952, [Glaser:1952]
Saturated vapour pressure there\(1.48\,\mathrm{MPa}\)est.
Pressure after release\(0.10(1)\,\mathrm{MPa}\)obs. 1952, [Glaser:1952]
Superheat $\Delta p$\(1.38\,\mathrm{MPa}\)est.
Surface tension at \(413\,\mathrm{K}\)\(3.9\times 10^{-3}\,\mathrm{N}/\mathrm{m}\)est., [Guggenheim:1945]
Critical bubble radius $r_{c}$\(5.7\times 10^{-9}\,\mathrm{m}\)calc., Equation (117.48)
Nucleation barrier $W_{c}$\(5.2\times 10^{-19}\,\mathrm{J}\), i.e. \(3.3\,\mathrm{eV}\)calc., Equation (117.50)
$W_{c}/k_{B}T$\(91.5\)calc.
Energy a minimum-ionising particle deposits over $2r_{c}$\(1.6\,\mathrm{eV}\), i.e. \(2.6\times 10^{-19}\,\mathrm{J}\)est.
Energy a \(5\,\mathrm{MeV}\) alpha deposits over $2r_{c}$\(1.6\times 10^{3}\,\mathrm{eV}\), i.e. \(2.6\times 10^{-16}\,\mathrm{J}\)est.
The superheated state in Glaser's bulbs, and the critical bubble it implies. Rows marked est. are computed here from corresponding-states scaling and the Clausius–Clapeyron relation and are order-of-magnitude statements, not measurements; the barrier depends on $\sigma^{3}$ and inherits three times the relative uncertainty of $\sigma$. The last two rows are the comparison that decides the mechanism: the barrier is far too high for a thermal fluctuation and far too low to stop a charged particle's energy deposit. The two temperatures and the released pressure are the operating point and are quoted with the uncertainty an operating point of this kind carries; the rows marked est. have no meaningful bracket, since their uncertainty is the validity of the scaling and not a measurement error, and the caveat above is the honest statement of it. The two rows marked tabulated are handbook properties of diethyl ether for which this bibliography carries no source.

The measured delay distributions are not transcribed. What the table establishes is the quantitative frame in which they are read, and the frame is one-sided. A barrier of \(91.5\,k_{B}T\) is essentially never crossed by a thermal fluctuation — the derivation below puts homogeneous nucleation at about one event per clean cubic centimetre per week — while it is crossed at once where a heavily ionising particle has passed. The delays of seconds to minutes actually observed in the shielded runs are therefore not homogeneous nucleation: they are set by residual heterogeneous sites and by the ambient background radiation, which is why the cleanliness of the bulb is the apparatus and why the comparison is made between distributions rather than against a computed spontaneous rate.

Interpretation

The bubble chamber is the cloud chamber with the phase transition run backwards, and the exchange is decisive: the sensitive medium is now a liquid, roughly a thousand times denser than a gas, so it is a target as well as a detector. Air at \(288.15\,\mathrm{K}\) has a density of \(1.225\,\mathrm{kg}/\mathrm{m}^{3}\) and diethyl ether about \(714\,\mathrm{kg}/\mathrm{m}^{3}\), a factor \(580\); liquid hydrogen, at \(70.8\,\mathrm{kg}/\mathrm{m}^{3}\) [Navas:2024], is still a factor \(58\). Interactions that a cloud chamber would see once a month happen continuously, and the whole event — production point, decay point, and every charged secondary — is recorded in one picture.

Filled with liquid hydrogen the chamber becomes a pure proton target, which is what made it the workhorse of hadron physics; the Gargamelle heavy-liquid chamber's photographs of neutrino interactions with no outgoing muon or electron [Hasert:1973b] are the discovery of the weak neutral current, and hence part of the evidential base of Experiment: Deep Inelastic Scattering and of the electroweak theory. That result is a good illustration of the instrument's character: the signature is the absence of a track, which is a statement no counting experiment could make and only a complete picture of the final state can support.

Its limitation is the one Charpak removed. A bubble chamber cannot be triggered: the liquid is sensitive only during the few milliseconds after the expansion, so it must be expanded on a schedule synchronised with the beam and photographed whether or not anything interesting happened, and the film must then be scanned and measured by human beings. The rate is set by mechanics and by scanning tables, not by physics.

Phenomenon 117.12 (Ionising radiation nucleates bubbles in a superheated liquid).

A liquid held above its boiling point without boiling begins to boil promptly when ionising radiation passes through it, and the bubbles form along the path of the ionising particle [Glaser:1952]. Rests on Equation (115.15).

Derivation. Derives Phenomenon 117.12. Mechanical equilibrium. A spherical bubble of radius \(r\) in a liquid holds its vapour at a pressure exceeding the liquid's by the Laplace term: balancing the force on a hemisphere, the pressure difference across the surface times the cross-section \(\pi r^{2}\) equals the surface tension \(\sigma\) times the circumference \(2\pi r\), so

\begin{equation}\tag{117.47} p_{\text{vap}}-p_{\text{liq}}=\frac{2\sigma}{r}\ec \end{equation}

in SI a pressure in pascals when \(\sigma\) is in \(\mathrm{N}/\mathrm{m}\) and \(r\) in metres. In a liquid superheated by \(\Delta p=p_{\text{vap}}-p_{\text{liq}}>0\), Equation (117.47) is satisfied at exactly one radius,

\begin{equation}\tag{117.48} r_{c}=\frac{2\sigma}{\Delta p} =\frac{2\times3.9\times 10^{-3}\,\mathrm{N}/\mathrm{m}} {1.38\times 10^{6}\,\mathrm{Pa}}=5.7\times 10^{-9}\,\mathrm{m}\ec \end{equation}

and that equilibrium is unstable: a bubble smaller than \(r_{c}\) has too much curvature, is crushed and collapses, while one larger than \(r_{c}\) finds the vapour pressure winning and grows without limit.

The barrier. Nothing happens spontaneously because reaching \(r_{c}\) costs work. The free energy of a bubble of radius \(r\) is the surface term \(4\pi r^{2}\sigma\) less the work \(\tfrac{4}{3}\pi r^{3}\Delta p\) done against the pressure difference,

\begin{equation}\tag{117.49} W(r)=4\pi r^{2}\sigma-\tfrac{4}{3}\pi r^{3}\Delta p\ec \end{equation}

whose derivative \(8\pi r\sigma-4\pi r^{2}\Delta p\) vanishes at \(r=r_{c}\), confirming Equation (117.48), and whose value there is the barrier

\begin{equation}\tag{117.50} W_{c}=\frac{16\pi\sigma^{3}}{3\left(\Delta p\right)^{2}} =5.2\times 10^{-19}\,\mathrm{J}=3.3\,\mathrm{eV}\ep \end{equation}

Compared with \(k_{B}T=5.70\times 10^{-21}\,\mathrm{J}\) at \(413\,\mathrm{K}\) this is \(91.5\,k_{B}T\), so the Boltzmann factor is \(\ee^{-91.5}=1.8\times 10^{-40}\): with of order \(10^{28}\) molecular sites per cubic metre and attempt frequencies of order \(10^{12}\,/\mathrm{s}\), the homogeneous nucleation rate is of order \(1.8\,/\mathrm{m}^{3}/\mathrm{s}\), that is \(1.8\times 10^{-6}\,/\mathrm{s}\) in a clean cubic centimetre — one event a week, against a cycle time of seconds. That is why the liquid can be kept superheated at all. It is also why the observed delays cannot be homogeneous nucleation: the calculation is not a prediction of the shielded-run delay distribution but a demonstration that homogeneous nucleation contributes nothing to it.

Why a track nucleates. A charged particle passing through deposits energy in a narrow cylinder along its path, in localised amounts far larger than thermal. Take a minimum-ionising particle in ether: the mass stopping power is about \(0.2\,\mathrm{MeV}\,\mathrm{m}^{2}/\mathrm{kg}\), so at a liquid density of \(714\,\mathrm{kg}/\mathrm{m}^{3}\) — a handbook value for which this bibliography carries no source — it loses about \(1.43\times 10^{8}\,\mathrm{eV}/\mathrm{m}\), that is \(2.29\times 10^{-11}\,\mathrm{J}/\mathrm{m}\), and over the critical diameter \(2r_{c}=11.3\,\mathrm{nm}\) it deposits about \(1.6\,\mathrm{eV}\), that is \(2.6\times 10^{-19}\,\mathrm{J}\) — which is below \(W_{c}\). A \(5\,\mathrm{MeV}\) alpha particle, with a range in ether of some tens of micrometres, deposits about \(10^{3}\) times more per unit length and delivers of order \(1.6\,\mathrm{keV}\), that is \(2.6\times 10^{-16}\,\mathrm{J}\), over the same distance, far above \(W_{c}\).

Two consequences follow and both are observed: the boiling delay collapses when a sufficiently ionising source is present, and the bubbles appear along the track, because that is where the energy was deposited. Note also from Equation (117.50) that the barrier falls as \(\left(\Delta p\right)^{-2}\), steeply with increasing superheat, which is why the chamber must be expanded into a narrow window and photographed at once — the operational counterpart of the two thresholds in Phenomenon 117.10.

Remark 117.13 (The energy budget is not as simple as the comparison above).

The comparison at the end of the derivation is the standard one and it is too generous. Not all the energy a particle deposits within a critical diameter is available to the nucleus: most of it goes into electronic excitation and into heating a region much larger than \(r_{c}\), and only the fraction thermalised inside the critical volume, within the time the bubble needs to reach \(r_{c}\), can pay for the barrier. The efficiency is of order a few per cent, which is why real chambers require a deposit one to two orders of magnitude above \(W_{c}\), and why the bubble density along a track is a usable measure of \(\dd E/\dd x\) — Equation (115.15) — rather than a threshold effect. The correct treatment follows the local temperature field produced by the track and asks when a supercritical bubble appears in it; it is a heat-conduction problem with a source, and it is what fixes the operating superheat of a real chamber.

Derivation pending.

The thermal-spike calculation of bubble nucleation: the temperature field around an ionising track from the heat equation with a line source, the fraction of the deposited energy thermalised within the critical radius before the bubble can grow, and the resulting threshold in \({[}\dd E/\dd x{]}\) as a function of the superheat. It belongs in Appendix A; the derivation above establishes the barrier and the mechanism but treats the energy delivery as instantaneous and perfectly local.

Primary references

Glaser's paper [Glaser:1952] is a two-page letter, and it is about the delay statistics of boiling in a superheated liquid; the track photographs that give the instrument its name came afterwards. It is worth reading as an example of an experiment whose published result is a shift in a distribution — there is no picture of a bubble in it — and whose importance is entirely in what the shift implies.

The instrument's scientific yield is documented elsewhere in this book rather than in Glaser's paper. The Gargamelle neutral-current photographs [Hasert:1973b] are the case discussed above; the hydrogen-chamber results on hadron spectroscopy underlie Experiment: Deep Inelastic Scattering, and the technique's decline is recorded in Section 117.7.

Charpak: the multiwire proportional chamber (1968)

Tests Phenomenon 117.14.

Apparatus

A plane of fine parallel anode wires, of the order of \(20\,\mu\mathrm{m}\) in diameter and spaced \(2\,\mathrm{mm}\) apart, stretched midway between two parallel cathode planes about \(8\,\mathrm{mm}\) away on either side, the whole in a gas mixture of a noble gas with a polyatomic quenching additive. The wires are held at a positive potential of a few kilovolts with respect to the cathodes.

Each of those dimensions is doing work. The wire must be thin because the field at its surface, and hence the avalanche gain, depends on its radius — Equation (117.54) — and a thick wire gives no multiplication at usable voltages. The spacing must be small because it sets the position resolution directly. The cathode gap must be several times the spacing, or the cells stop being independent and the field near each wire is distorted by its neighbours. And the quenching additive is not optional: a pure noble gas de-excites by emitting ultraviolet photons that liberate photoelectrons from the cathode metres away, so the avalanche propagates and the counter goes into permanent discharge. A polyatomic molecule absorbs those photons into vibrational modes and dissipates them, confining the avalanche to the wire that started it.

The decisive feature, however, is not the mechanics but the readout: every wire has its own amplifier and its own discriminator, so the chamber emits a pattern of logic pulses rather than an image. In 1968 that was newly affordable, and it is the reason the paper was written when it was.

Procedure

A charged particle crossing the gas leaves a track of electron–ion pairs — of order \(100\) of them in \(8\,\mathrm{mm}\) of argon at atmospheric pressure, at a mean expenditure of about \(26\,\mathrm{eV}\), that is \(4.2\times 10^{-18}\,\mathrm{J}\), per pair [Navas:2024]. The electrons drift towards the nearest anode wire through the low-field region between the cathodes, and in the last few tens of micrometres, where the field rises as \(1/r\) and exceeds the threshold for ionisation by collision, each initiates an avalanche whose gain is set by the wire voltage and the gas.

The measurements reported are those any new counter owes. The gas gain was measured against anode voltage, by comparing the collected charge with the known primary ionisation. The detection efficiency was measured against voltage with an external trigger — a scintillator telescope defining a particle's passage — which produces the characteristic plateau, a range of voltage over which efficiency is near unity and insensitive to the setting, and finding that plateau is what makes a counter usable. The time resolution was measured against the same trigger. And, the result that was not obvious in advance, the pulses on wires adjacent to the struck one were recorded, to establish that a chamber with wires two millimetres apart really does behave as a row of independent counters and not as one large one.

Observations and data

Each wire behaved as an independent proportional counter, with adjacent wires showing negligible cross-talk despite being only millimetres apart. Detection efficiency was high over a plateau of several hundred volts, the position was given to within the wire spacing, and both the time resolution — of the order of the maximum electron drift time computed in Table 117.7, which is what bounds it — and the rate a single wire could sustain were better than any photographic technique by orders of magnitude.

Uncertainties are owed on all of this and the paper quotes none, so what Table 117.7 carries are estimates from the mechanical tolerance of a wound wire plane and from the stability of a high-voltage supply, each flagged in the table. The efficiency and the time resolution are reported here as the paper reports them, in words and without a figure, because assigning them a numerical uncertainty would mean inventing the measurement they belong to; what the argument below needs from them is only that the plateau exists and that the resolution is drift-limited, and both are qualitative statements.

Table 117.7 gives the geometry, the reported performance, and the electrostatic quantities those imply, computed here from Equations (117.53) and (117.54) and from Equations (117.56) and (117.57). Table 117.8 is the avalanche gain against anode voltage, computed from Equation (117.55) with representative Townsend coefficients for argon; it is included because the steepness of that dependence is the whole reason a proportional chamber has an operating plateau at all, and because it is a closed-form consequence of the \(1/r\) field.

QuantityValueKind, source
Anode wire diameter\(20(1)\,\mu\mathrm{m}\)obs. 1968, [Charpak:1968]
Wire spacing $s$\(2.00(5)\,\mathrm{mm}\)obs. 1968, [Charpak:1968]
Anode–cathode gap $L$\(8.0(1)\,\mathrm{mm}\)obs. 1968, [Charpak:1968]
Anode potential\(3\,\mathrm{kV}\) to \(4\,\mathrm{kV}\), stable to \(10\,\mathrm{V}\)obs. 1968, [Charpak:1968]
Detection efficiency, over a plateau several hundred volts widenear unityobs. 1968, [Charpak:1968]
Capacitance per unit length of one wire\(3.47\,\mathrm{pF}/\mathrm{m}\)calc., Equation (117.53)
Line charge at \(3\,\mathrm{kV}\)\(1.04\times 10^{-8}\,\mathrm{C}/\mathrm{m}\)calc.
Field at the wire surface\(1.87\times 10^{7}\,\mathrm{V}/\mathrm{m}\)calc., Equation (117.54)
Radius within which multiplication occurs\(62\,\mu\mathrm{m}\)calc.
Primary ion pairs in the gapabout \(100\)est.
Signal charge at gain \(10^{4}\)\(1.6\times 10^{-13}\,\mathrm{C}\)calc.
Position resolution\(0.58\,\mathrm{mm}\)calc., Equation (117.56)
Maximum electron drift time\(160\,\mathrm{ns}\)calc.
Ion signal time constant $t_{0}$\(1.8\,\mathrm{ns}\)calc., Equation (117.57)
Total ion collection time\(1.14\,\mathrm{ms}\)calc., Equation (117.58)
The multiwire proportional chamber: geometry, performance and the electrostatics that follow. Geometry and performance are the reported values [Charpak:1968]; the capacitance, line charge, surface field, multiplication radius, position resolution and drift and ion times are computed here in the derivation below. The signal charge assumes a gain of \(10^{4}\). Bracketed uncertainties on the observed geometry are estimates from the mechanical tolerance of a wound wire plane, the paper quoting none; the computed rows are exact consequences of that geometry and inherit it, the capacitance and everything built on it to about \(1.5\,\mathrm{\%}\). The efficiency row carries no bracket for the reason given in the text, and the primary ion-pair count is an order-of-magnitude estimate while the gain is a chosen reference value.
Anode potentialSurface fieldMultiplication radiusGain $M$
\(2.0\,\mathrm{kV}\)\(1.25\times 10^{7}\,\mathrm{V}/\mathrm{m}\)\(42\,\mu\mathrm{m}\)\(1.5\times 10^{1}\)
\(2.5\,\mathrm{kV}\)\(1.56\times 10^{7}\,\mathrm{V}/\mathrm{m}\)\(52\,\mu\mathrm{m}\)\(6.8\times 10^{1}\)
\(3.0\,\mathrm{kV}\)\(1.87\times 10^{7}\,\mathrm{V}/\mathrm{m}\)\(62\,\mu\mathrm{m}\)\(3.6\times 10^{2}\)
\(3.5\,\mathrm{kV}\)\(2.18\times 10^{7}\,\mathrm{V}/\mathrm{m}\)\(73\,\mu\mathrm{m}\)\(2.1\times 10^{3}\)
\(4.0\,\mathrm{kV}\)\(2.50\times 10^{7}\,\mathrm{V}/\mathrm{m}\)\(83\,\mu\mathrm{m}\)\(1.3\times 10^{4}\)
\(4.5\,\mathrm{kV}\)\(2.81\times 10^{7}\,\mathrm{V}/\mathrm{m}\)\(94\,\mu\mathrm{m}\)\(8.1\times 10^{4}\)
Avalanche gain against anode potential, computed here from Equation (117.55) for the geometry of Table 117.7 in argon at \(101.3\,\mathrm{kPa}\), with Townsend coefficients $A=9.0\,/\mathrm{m}/\mathrm{Pa}$ and $B=135\,\mathrm{V}/\mathrm{m}/\mathrm{Pa}$ representative of argon, and a multiplication threshold $E^{*}=3\times 10^{6}\,\mathrm{V}/\mathrm{m}$. The gain rises by nearly four orders of magnitude for a factor \(2\) in voltage, which is why the operating point is chosen on an efficiency plateau and not by tuning the gain.

Interpretation

This is the change of era. A cloud or bubble chamber produces a photograph, and a photograph must be developed, scanned and measured by a person; the experiment's rate is therefore set by human labour. A bubble chamber cycles at a few hertz at best, and a scanner works through pictures at a rate of order one per minute, so a large experiment of the nineteen-sixties accumulated of order \(10^{6}\) pictures over years and employed rooms full of people to measure them. A multiwire chamber produces a number, immediately, on a wire, and a computer can read it; a single wire sustains of order \(10^{5}\,/\mathrm{s}\) and the limit is the electronics.

Triggering becomes possible with it — the apparatus can be told to record only events satisfying a condition, evaluated in the time between one particle and the next — and with it the whole strategy of searching for rare processes in an overwhelming background, which is how every result in Experiment: The Higgs Boson Discovery was obtained. Modern detectors record of order \(10^{9}\) interactions per second and select a handful for storage, and the reason they can is that the sensitive element gives an electronic address rather than an image.

The device also founded a family. Measuring the drift time of the electrons, rather than merely which wire fired, turns the wire spacing into a far finer coordinate: an electron drifting at \(5\times 10^{4}\,\mathrm{m}/\mathrm{s}\) covers \(50\,\mu\mathrm{m}\) in a nanosecond, so timing to a nanosecond gives a resolution an order of magnitude better than Equation (117.56) — that is the drift chamber. Extending the drift volume to metres and reading the arrival pattern on an endplate gives the time projection chamber; replacing wires with etched electrodes on a substrate gives the micropattern detectors of current practice. All of them are read out as numbers, and all of them descend from this experiment.

Phenomenon 117.14 (Each anode wire is an independent counter).

In a plane of parallel anode wires only millimetres apart between two cathode planes, the avalanche produced by an ionising particle develops on a single wire and induces a signal there, with negligible response on its neighbours. The chamber therefore reports position as the address of the wire that fired, at a rate limited by electronics rather than by photography [Charpak:1968]. Rests on Theorem 58.10.

Derivation. Derives Phenomenon 117.14. The field. Model one wire as a long straight conductor of radius \(a\) carrying charge \(\lambda\) per unit length. Gauss's theorem over a coaxial cylinder gives

\begin{equation}\tag{117.51} E(r)=\frac{\lambda}{2\pi\varepsilon_{0}r}\ec \end{equation}

in \(\mathrm{V}/\mathrm{m}\) when \(\lambda\) is in \(\mathrm{C}/\mathrm{m}\), and hence a potential varying only logarithmically, \(V\propto\ln(r/a)\). The periodic geometry has to be summed, and the sum is worth doing because it fixes \(\lambda\) and hence every electrostatic number in this section. Put the wires along \(y\) at \(y=ns\), \(n\in\Z\), each carrying \(\lambda\), and let \(x\) measure distance from the wire plane. Superposing the logarithmic potentials of the whole row and using the product formula \(\sin\left(\pi w\right)/\left(\pi w\right) =\prod_{n\ge1}\left(1-w^{2}/n^{2}\right)\), applied to the complex combination \(w=\left(y+\ii x\right)/s\), gives the closed form

\begin{equation}\tag{117.52} V(x,y)=-\frac{\lambda}{4\pi\varepsilon_{0}} \ln\left[2\left(\cosh\frac{2\pi x}{s} -\cos\frac{2\pi y}{s}\right)\right] +\text{const}\ep \end{equation}

Two limits of Equation (117.52) are all that is needed. Close to one wire, \(r=\sqrt{x^{2}+(y-ns)^{2}}\ll s\), expanding both hyperbolic and circular functions to second order gives \(2(\cosh-\cos)\to\left(2\pi r/s\right)^{2}\) and hence \(V\to-\left(\lambda/2\pi\varepsilon_{0}\right)\ln\left(2\pi r/s\right)\), which recovers Equation (117.51) and shows that the row looks like a single wire from close up. Far from the plane, \(\abs{x}\gg s\), the cosine is negligible against \(\cosh(2\pi x/s)\to \tfrac{1}{2}\ee^{2\pi\abs{x}/s}\), and the potential becomes linear, \(V\to-\left(\lambda/2\varepsilon_{0}\right)\abs{x}/s\): the wire plane looks like a charged sheet, as it must. That second limit is the reason the cathode gap has to be several times the wire spacing: at the \(L/s=4\) of this chamber the terms dropped are smaller than the one kept by a factor \(2\times 10^{-11}\), while at \(L/s=1\) they would not be negligible and the cells would stop being independent. Subtracting the second limit at the cathode, \(\abs{x}=L\), from the first at the wire surface, \(r=a\), gives the working potential \(V=\left(\lambda/2\pi\varepsilon_{0}\right) \left[\pi L/s-\ln\left(2\pi a/s\right)\right]\), so that the capacitance per unit length of one wire is

\begin{equation}\tag{117.53} C=\frac{\lambda}{V}=\frac{2\pi\varepsilon_{0}} {\dfrac{\pi L}{s}-\ln\dfrac{2\pi a}{s}} =3.47\times 10^{-12}\,\mathrm{F}/\mathrm{m} \end{equation}

for \(a=10\,\mu\mathrm{m}\), \(s=2\,\mathrm{mm}\) and \(L=8\,\mathrm{mm}\), so at \(V=3\,\mathrm{kV}\) the line charge is \(\lambda=CV=1.04\times 10^{-8}\,\mathrm{C}/\mathrm{m}\) and the field at the wire surface is

\begin{equation}\tag{117.54} E(a)=\frac{\lambda}{2\pi\varepsilon_{0}a} =1.87\times 10^{7}\,\mathrm{V}/\mathrm{m}\ep \end{equation}

Two facts follow from the \(1/r\) dependence, and between them they are the whole device.

First: the multiplication is local. Avalanche multiplication requires an electron to gain the ionisation energy of the gas between collisions, that is a field above some threshold \(E^{*}\) of order \(3\times 10^{6}\,\mathrm{V}/\mathrm{m}\) in argon at atmospheric pressure; by Equation (117.51) the region where \(E>E^{*}\) is a cylinder of radius \(\lambda/2\pi\varepsilon_{0}E^{*}=62\,\mu\mathrm{m}\), a few wire radii. The multiplication is therefore confined to the immediate neighbourhood of one wire, thirty times closer to it than to its neighbours, and the gain does not depend on where in the gas volume the original ionisation occurred — which is what makes the signal proportional to the deposited charge, and what keeps the wires independent.

The gain itself follows in closed form. The first Townsend coefficient — the number of ionising collisions an electron makes per metre — is well represented by \(\alpha=Ap\exp\left(-Bp/E\right)\) with \(p\) the pressure, and the multiplication is \(M=\exp\int\alpha\,\dd r\) taken inward from \(E=E^{*}\) to \(E=E(a)\). Changing variable from \(r\) to \(E\) using Equation (117.51), so that \(r=k/E\) with \(k=\lambda/2\pi\varepsilon_{0}\) and \(\dd r=-k\,\dd E/E^{2}\), the integral is elementary:

\begin{equation}\tag{117.55} \ln M=Apk\int_{E^{*}}^{E(a)}\frac{\ee^{-Bp/E}}{E^{2}}\,\dd E =\frac{Ak}{B} \left(\ee^{-Bp/E(a)}-\ee^{-Bp/E^{*}}\right)\ep \end{equation}

With \(A=9.0\,/\mathrm{m}/\mathrm{Pa}\) and \(B=135\,\mathrm{V}/\mathrm{m}/\mathrm{Pa}\) for argon, \(p=1.013\times 10^{5}\,\mathrm{Pa}\) and \(k=\lambda/2\pi\varepsilon_{0}\) proportional to \(V\) through Equation (117.53), this generates Table 117.8: a gain of \(15\) at \(2\,\mathrm{kV}\) and \(1.3\times 10^{4}\) at \(4\,\mathrm{kV}\). The extreme steepness is structural — \(\ln M\) is roughly linear in \(V\), so \(M\) is exponential in it — and it is why the operating point of a proportional chamber is chosen by finding the efficiency plateau rather than by setting a gain.

Second: the cells are defined by the field, not by walls. The drift region between the cathodes is a low-field region in which electrons simply move to the nearest wire along the field lines, and the division of the plane into cells belongs to the field geometry: no walls are needed between adjacent counters, which is why the wires can be packed at millimetre pitch. A particle's coordinate is then the address of the wire that fired, uniformly distributed over the cell of width \(s\), so the resolution is the standard deviation of a uniform distribution,

\begin{equation}\tag{117.56} \sigma_{x}=\frac{s}{\sqrt{12}} =\frac{2\,\mathrm{mm}}{3.464} =0.58\,\mathrm{mm}\ep \end{equation}

Third: the signal is fast because it is induced by the ions. The electrons produced in the avalanche are collected in a fraction of a nanosecond and induce almost nothing, having been created within a few micrometres of the wire; the signal is the induced charge from the positive ions drifting away from the wire through the whole logarithmic field. An ion of mobility \(\mu\) obeys \(\dd r/\dd t=\mu E(r)=\mu k/r\), so \(r^{2}(t)=a^{2}+2\mu kt\) and the induced charge grows as \(\ln\left(1+t/t_{0}\right)\) with

\begin{equation}\tag{117.57} t_{0}=\frac{a^{2}}{2\mu k} =\frac{\pi\varepsilon_{0}a^{2}}{\mu\lambda} =1.8\,\mathrm{ns}\ec \end{equation}

for \(\mu=1.5\times 10^{-4}\,\mathrm{m}^{2}/\mathrm{V}/\mathrm{s}\) [Navas:2024]. The signal therefore rises on a nanosecond scale even though full ion collection is slow. The same \(r^{2}(t)\) run out to the cathode, at \(r=L\) and not \(L/2\) — \(L\) is the whole anode-to-cathode gap of Table 117.7, which is also the distance the electron drift time of that table is computed over — gives

\begin{equation}\tag{117.58} t_{\max}=\frac{L^{2}-a^{2}}{2\mu k} \approx\frac{\pi\varepsilon_{0}L^{2}}{\mu\lambda} =1.14\,\mathrm{ms}\ec \end{equation}

six hundred thousand times \(t_{0}\), since \(t_{\max}/t_{0}=(L/a)^{2}\). Differentiating the signal with a short time constant keeps the fast leading edge and discards the long tail, which is what makes the coordinate available promptly enough to be used in a trigger. The residual limit on rate is the space charge of those slow ions, which is what fixes the \(10^{5}\,/\mathrm{s}\) per wire quoted above.

Remark 117.15 (What is representative and what is measured here).

The Townsend coefficients used in Equation (117.55) and the ion mobility used in Equation (117.57) are representative values for argon at atmospheric pressure, not values measured in Charpak's chamber, which ran on a gas mixture whose composition changes both by factors of order unity. They are used to show that the closed form reproduces the observed orders of magnitude — a gain of \(10^{4}\) at a few kilovolts, a multiplication region of tens of micrometres, a nanosecond leading edge — and the argument does not depend on their precise values. The measurement of electron transport in argon on which such coefficients ultimately rest is Townsend and Bailey's [Townsend:1922].

Primary references

The paper is Charpak and collaborators' [Charpak:1968], in Nuclear Instruments and Methods, and it is short. Its significance is easy to miss on a first reading because it claims very little: that a plane of closely spaced wires can be operated with each wire as an independent proportional counter, with the gains, efficiencies and time resolutions tabulated. The claim that was not obvious — and that the paper's referees would have been right to doubt — is the independence of the wires at millimetre spacing, which Equation (117.54) makes plausible and the measurements established.

What followed is not in this paper and is not any single paper. The technique displaced the bubble chamber over the following decade, and the detectors of Experiment: Deep Inelastic Scattering, Experiment: Neutrino Oscillations and Experiment: The Higgs Boson Discovery are all, in their tracking elements, its descendants.

Summary of the evidence

ExperimentYearResult
Rutherford1919first artificial transmutation; nitrogen yields protons of up to \(5.93\,\mathrm{MeV}\) against an elastic ceiling of \(4.94\,\mathrm{MeV}\)
Cockcroft–Walton1932first accelerator-driven nuclear reaction, and the first quantitative nuclear test of $E=mc^{2}$: alpha particles of about \(8.6\,\mathrm{MeV}\) measured against \(8.79\,\mathrm{MeV}\) predicted from a mass difference of \(3.09\times 10^{-29}\,\mathrm{kg}\)
Lawrence–Livingston1932protons of \(1.22\,\mathrm{MeV}\) from \(4\,\mathrm{kV}\) per gap, in \(153\) turns at \(1.40\,\mathrm{T}\)
Wilson1912individual charged particles made visible as droplet trails, in a supersaturation window between $S=4.1$ and $S=8.1$
Glaser1952the sensitive medium becomes a liquid \(580\) times denser than air, hence target and detector at once
Charpak1968position read as an electronic address to \(0.58\,\mathrm{mm}\), at \(10^{5}\,/\mathrm{s}\) per wire instead of a photograph
The reactions and the instruments reported in this chapter. The first two rows are nuclear reactions, the third is the machine that freed the projectile energy from the insulation limit, and the last three are the detectors that made reactions countable one at a time. The three detectors are best read by what each removed: the invisibility of the single particle, the thinness of the target, and the photograph.

Table 117.9 contains no fitted parameter and no theoretical input beyond mass–energy equivalence, elementary electrostatics, the semiclassical barrier transmission of Section 9.6, and classical nucleation theory. The first two rows establish that nuclei can be rearranged at will and that the energy bookkeeping of Relativistic Dynamics holds when they are; the third makes the projectile energy an adjustable quantity rather than a gift of nature; the last three establish that the rearrangements can be counted, one at a time, at rates that rose by roughly six orders of magnitude across the half-century the table spans. Neither half is of any use without the other, which is why they are reported together here rather than split between a nuclear chapter and a detector appendix.

Two threads run through all six and are worth naming, because they recur in every later experiment in this part. The first is that the decisive measurement is almost never the one the experiment appears to be making. Rutherford's result is a comparison between gases, not a count; Cockcroft and Walton's is a coincidence rate, not a yield; Lawrence's is the width of a resonance curve, which is what proves the turn count; Glaser's is a shift in a distribution of delays. In each case the number that carries the physics is a ratio or a correlation, chosen because it is insensitive to the quantities the experimenter cannot control.

The second is that every one of the six instruments is a threshold device operated inside a deliberately chosen window. The cloud chamber must sit between two supersaturations, the bubble chamber between two superheats, the proportional chamber between the onset of multiplication and the onset of discharge, the scintillation counter between the eye's sensitivity and its fatigue, the cyclotron between resonance and relativistic detuning, the cascade generator between usable voltage and breakdown. Locating those windows is what the Procedure of each experiment above actually describes, and it is the part of experimental physics that no theory supplies.