Design Guide › Physics theory
Cyclotron physics theory design rules
77 of the guide’s 2170 rules carry the physics-theory tag.
Rules that are theoretical results rather than practice: stripping lifetimes, measurement accuracy floors, kinematic broadening, and detector efficiency models.
Each rule keeps its formula where the source gives one, a verbatim quote, a page-level
citation, and a stable identifier (dg-NNNN) that resolves here and on the
all-in-one guide. Where an editorial note says
“the reference machine”, its parameters are on the
guide’s front page.
By applicability level: level 1 (3) · level 2 (21) · level 3 (30) · level 4 (23) — levels rank breadth, never license to skip (method). Related domains, by shared rules: Beam dynamics (22), Ion source (19), Vacuum (18), RF (14), Magnet (8). To combine tags or levels, open this domain in the filterable view.
Verify before use. Every rule here is a source extract in the vocabulary of the editorial methodology — faithful to its cited page, not an independently validated engineering requirement. Re-read any rule that drives a real design decision at the cited page before committing metal, money, or high voltage to it. The editorial note under each quote is this site’s extrapolation to a tabletop machine, not something the source said: an editor’s judgement, audited for overreach, never a citation.
-
Lorentz (magnetic) stripping of H- has rest-frame lifetime tau = (A1/E)*exp(A2/E) with A1=2.714e-6 s*V/m, A2=4.474e9 V/m, E=gamma*beta*c*B - negligible below a few MeV even at 4 T.
tau = (A1/E)*exp(A2/E), E = gamma*beta*c*BSource quote & editorial note
a 4 T magnetic field for the maximum achievable energy of 8.5 MeV in the AMIT cyclotron corresponds to a beam-rest-frame electric field of E = 160 MV/m. This entails a marginal beam fraction loss per unit length of 1.42e-6 m-1
Calvo et al., Beam Stripping Interactions in Compact Cyclotrons — PRAB 24, 090101 (2021) — p. 7-8, 14
Editorial note, tabletop extrapolation: At 0.889 T and 500 keV the rest-frame field is ~9 MV/m, where the exponential makes the lifetime effectively infinite - Lorentz stripping can be ignored entirely for a next machine.
-
Know the accuracy floor of the cited absorber-based measurement: range-energy data and straggling limited the most-probable-energy determination to a few hundred keV, with the high-energy portion nearly as good, and the low-energy portion involving considerably greater uncertainty.
Source quote & editorial note
These factors limit the accuracy of determination of the most probable energy to a few hundred kilovolts. The high energy portion of the energy distribution can be determined with almost equivalent accuracy
Editorial note, tabletop extrapolation: The asymmetry - high-energy side of an absorber spectrum better determined than the low-energy tail - is the shape to remember, but the historical few-hundred-keV floor belongs to that apparatus: for a PIPS-plus-degrader setup, build the detector-and-degrader response matrix and quote separate uncertainties for mode, upper edge and tail rather than scaling ORNL's numbers.
-
Opening a spectrograph's in-plane angular acceptance costs kinematic broadening of peaks when scattering off light nuclei; the cited instrument accepted that trade for large-angle reach.
Source quote & editorial note
In scattering from light nuclei, this introduces appreciable kinematic broadening of peaks as the entrance aperture is opened.
Editorial note, tabletop extrapolation: For a next machine's Rutherford-scattering station the same knob exists: closing the entrance aperture trades count rate for resolution. State the broadening honestly - dE ~ |dE/dtheta|*dtheta depends on beam energy, angle and kinematics in general; only the FRACTIONAL elastic broadening at fixed masses and angle drops out energy-independent - so compute it for the actual geometry rather than quoting a mass-ratio shortcut.
-
A simple single-scattering model predicts organic-scintillator neutron efficiency usefully: eff = (1 - E0/En)*(1 - exp(-nH*sigma_np*l)). The source calls the first factor exact and response-shape-independent - which holds given its implicit assumptions: isotropic center-of-mass n-p scattering (uniform recoil spectrum), a single hydrogen scatter, and a sharp threshold; find E0 per threshold setting from a calibration reaction.
eff = (1 - E0/En)(1 - exp(-n_H*sigma_np(En)*l)); valid as a first-order model for thin hydrogenous scintillators; degrades with n-p anisotropy at higher energy, carbon interactions, proton escape, multiple scatteringSource quote & editorial note
The first factor in (1) is exact. It does not depend on the exact shape of the response curve (pulse-height vs proton recoil energy) of the scintillator.
Fulbright et al., A Fast Neutron Time of Flight System for Use with Cyclotrons — NYO-9360 (1962) — p. 15
Editorial note, tabletop extrapolation: A lovely teaching derivation - a two-factor closed form students can test against a calibration reaction - and the transferable habit is identifying which factor rests on which assumption. Validate against calibration or Monte Carlo for the actual detector and energy range before leaning on it.
-
Choose design beam energy from the reaction excitation curve: the thesis justifies its 150 keV deuteron design energy by placing the d(d,n)3He cross-section 'near the maximum' of its plotted curve, calling the reaction exothermic with - as printed - '2.227 MeV released for each deuteron pair'. [Source-internal error, surfaced: 2.227 MeV is approximately the DEUTERON BINDING energy; the D(d,n)3He Q-value is 3.27 MeV and D(d,p)3H is 4.03 MeV. And the D-D cross-section keeps rising well beyond 150 keV - 'near the maximum' holds only within the thesis's plotted range.]
Source quote & editorial note
The maximum energy of the current cyclotron, using (6), is 150 keV for deuterons, which puts the cross section near the maximum. ... [the d(d,n)3He reaction is] exothermic with 2.227 MeV released for each deuteron pair.
Cressman, The Design and Construction of a Small Cyclotron — Houghton College thesis (2006) — p. 49
Editorial note, tabletop extrapolation: Working back from the excitation curve converts machine energy from a bragging number into a requirement - the transferable design move. D-D is the standout low-energy neutron reaction because its cross-section is already usable near 100 keV; take Q-values and cross-sections from live evaluated data (per site policy), not from the thesis's figures.
-
Practical relativistic ceiling as the thesis cites it: beyond raising electrode voltage, relativity can be countered by shaping the magnetic field, and 'the actual relativistic limit for magnetic resonators accelerating protons is closer to 25 MeV' - a historical (Rose-era) estimate for that machine class, not a universal constant.
Source quote & editorial note
The actual relativistic limit for magnetic resonators accelerating protons is closer to 25 MeV.
Cressman, The Design and Construction of a Small Cyclotron — Houghton College thesis (2006) — p. 28
Editorial note, tabletop extrapolation: Contextualizes the phase-slip table (dg-1468): field shaping buys real headroom beyond the uniform-field estimate - how much depends on the field design, so don't carry a fixed multiplier. For any machine, calculate cumulative phase slip from the actual energy gain per turn and B(r) rather than trusting an energy-class exemption; slip can bite below 1 MeV when the gain per turn is small.
-
The source's stated approximation for cyclotron-mass-spectrometer resolution is R ≈ 3 n H, with n the number of in-phase turns before extraction and H the RF harmonic number — so resolution is bought with more turns or a higher harmonic, each carrying its cost elsewhere in the design (the source's center-region compromise, dg-1556).
R ~ 3 * n * H (n = turns before extraction, H = RF harmonic)Source quote & editorial note
a mass resolution of about 1800 is needed to separate 14 C from 13 CH. The resolution of a CMS is approximately: R ≈ 3 x n x H, where n is the number of turns that in-phase particles make in a synchronous field before extraction and H is the rf harmonic number
Editorial note, tabletop extrapolation: Editorial note, tabletop extrapolation: an estimating relation, not a law — the coefficient depends on the phase-slip criterion that defines an in-phase turn. Useful for order-of-magnitude estimates of how sharply a small machine discriminates species or off-resonance drive; derive the real number from a phase-history calculation for the actual field and RF program.
-
Rose-type phase integral (as applied to the ISU cyclotron, 1961): with u = sin(theta) the phase lag, tune-down profile deltaB(r) = B1 - B(r), field index n = -(r/B)(dB/dr), and V0 the peak dee-to-ground voltage, du/dr = pi*e*r*B*deltaB*(1-n)/(2*m*V0). Integrating from the measured B(r) gives the phase-lag curve for any initial phase on the accelerating branch (-pi/2 < theta < pi/2); the modeled solution remains admissible while -1 < u < 1, u = +/-1 being the model's phase-loss boundary.
u = (pi*e/(2*m*V0)) * integral_0_to_r [ r*B*(B1-B)*(1-n) ] dr + u0, with u = sin(theta)Source quote & editorial note
du/dr = πerB∆B(1−n)/(2mV0). This equation gives the rate of change of the sine of the phase lag, θ, as a function of r and the magnetic field, B. Integration gives u = (πe/2mV0) ∫ rB∆B(1−n) dr + u0. (1) From this result the phase of the proton can be obtained at any radius if the initial phase lag and the magnetic field are known
Editorial note, tabletop extrapolation: Editorial note, tabletop extrapolation: this single quadrature over the measured field map predicts phase history without tracking orbits, and it runs in a spreadsheet — under the model's assumptions (nonrelativistic centered orbits, continuous acceleration, initial phase restricted to the accelerating branch). The right first tool for choosing frequency and trim before any trajectory code is written; it bounds phase admissibility only — vertical loss, radial loss and scattering are separate budgets.
-
Median-plane field expansion used for axial motion (ISU, 1963): from curl B = 0 and div B = 0 in the dee box, with pole symmetry giving Br = 0 on the median plane, the small-z approximation is Br = -z*dB/dr with B taken independent of z - reducing the axial equation of motion to z'' - (e*v/m)*(dB/dr)*z = 0 driven entirely by the median-plane gradient.
Br = -z*dB/dr (small z); axial equation m*z'' = -r*phidot*e*BrSource quote & editorial note
Since Br is equal to zero on the median plane, the following approximation is valid for small axial displacements
Editorial note, tabletop extrapolation: This is why a median-plane-only field survey suffices for a first-cut axial-focusing model - the off-plane field follows from Maxwell to first order in z. The 1963 authors also flagged its limit - the approximation degrades for large axial amplitudes and unknown off-plane field shape.
-
The Rutgers group point out that an E x B channel embedded in a cyclotron cannot resolve the q/m ambiguity between fully ionized deuterium (2H+) and helium (4He++), because the cyclotron itself acts as a velocity filter at each radius: at fixed magnetic field both species have the same resonant frequency and the same angular velocity, hence the same velocity at the channel entrance, though not the same energy.
f_cyc = (B/2*pi)*(q/m)Source quote & editorial note
It might be expected that the combined effect of the cyclotron's resonant acceleration and our embedded Wien Filter's velocity selection might separate the mass ambiguity. However, this is not the case, as the cyclotron itself acts a velocity filter at each of its radii. […] It also holds that at any given radius, both the deuterium and helium cover the same angular distance and thus must have the same angular velocity to keep in step with the oscillating RF voltage. Now it is easily seen that the velocity of either deuterium or helium will be the same at the entrance to the deflection channel. Note, while the two ions have the same velocity, they obviously do not have the same energy.
Editorial note, tabletop extrapolation: Directly relevant to a hydrogen-fed tabletop machine, where the contaminant species of interest are H2+ and H3+ rather than deuterium and helium: an internal E x B channel will separate q/m = 1 from q/m = 1/2, but will not distinguish two species sharing a q/m. Species identification has to come from elsewhere (gas fill, source chemistry, or a downstream detector), a limit worth knowing before building the channel as a diagnostic.
-
Off-harmonic operation observed and rationalized on the Rutgers 12-inch: because a cyclotron only resonantly accelerates at odd integer harmonics, operating near but not on an odd harmonic can still give a successfully accelerated beam provided the integrated phase slippage over all revolutions is less than 180 degrees before the target or extraction point — and since higher DEE voltage means fewer revolutions to reach a given energy, the tolerable phase slippage per turn increases with DEE voltage.
Source quote & editorial note
This result is not understood, as only integer odd harmonic numbers support magnetic resonance acceleration. At an even harmonic, when acceleration occurs at a gap crossing, deceleration must occur at the subsequent crossing, yielding zero net accelerator per revolution. In the region between an even and odd harmonic, there is a balance of acceleration and phase slippage which the ions encounter. Operating a cyclotron near, but not on, an odd harmonic, can still lead to a successful resonantly accelerated beam, provided that the integrated phase slippage over all revolutions is less than 180 before hitting the target or extraction point. The greater the DEE voltage, the fewer the number of ion revolutions are needed to achieve the desired energy, thus the tolerance of phase slippage per turn increases with DEE voltage.
Editorial note, tabletop extrapolation: Directly relevant to low-dee-voltage machines, in mirror image: many hundreds of turns means very little tolerable slip per turn, which is an operational argument for dee voltage beyond simple turn-count. State the physics as the source's gap phasing gives it: odd harmonics are the resonant condition for this conventional geometry, and near-harmonic operation can survive if the bunch stays inside the accelerating phase window — the 180-degree integrated-slip figure is an approximate span, conditional on where in phase the ions start and which way they slip. The reported 2.22 and 4.25 harmonic numbers are stated by the authors as "not understood" — an open anomaly, not a result. ("accelerator per revolution" is the source's typo for "acceleration".)
-
Complete transverse stability in a constant-gradient (weak-focusing) cyclotron requires 0 < n < 1, where the field index n = -(r/B)(dB/dr); the axial tune is nu_z = sqrt(n) and the radial tune is nu_x = sqrt(1-n), both from the Kerst-Serber equation.
n = -(r/B)(dB/dr); d2z/dt2 + n w^2 z = 0; d2x/dt2 + w^2 (1-n) x = 0; nu_z = sqrt(n); nu_x = sqrt(1-n)Source quote & editorial note
Complete transverse stability. It has thus been shown for axial stability, n must be greater than 0, and for radial stability n must be less than 1. Total transverse stability exists in the region of: 0 < n <1
Editorial note, tabletop extrapolation: The design inequality for a weak-focusing machine, derived in this report from scratch: away from the central region, 0 < n < 1 buys simultaneous linear axial and radial stability (at r = 0 itself n = 0, as the source's own next passage states — the center is handled by other means, dg-1729/dg-1835). It is a LOCAL linear-stability window: resonances (dg-1682), acceleration and field errors still get their say. Sign convention: this document's leading minus makes n > 0 a falling field; the companion AVF paper uses k = d ln⟨B⟩/d ln R with opposite sign, so reconcile n = −k before mixing formulas.
-
For the Rutgers 12-inch AVF work the axial tune is written nu_z^2 = -k + F(1+tan^2 xi) and the radial tune nu_r^2 = 1+k, where k = d ln<B> / d ln R is the average field index, F is the rms flutter (the rms azimuthal variation of the vertical field) and xi is the instantaneous angle the sector edge makes with the orbit.
nu_z^2 = -k + F(1+tan^2 xi); nu_r^2 = 1+k; k = d ln<B> / d ln RSource quote & editorial note
AVF focusing can be used to supplement weak focusing. In this context, the weak focusing comes from the average radial gradient’s field index, denoted as k, where: k = d ln〈B〉/ d ln R . The tune is proportional to the relative focusing strength. Following the treatment of J.J. Livingood,[6] one can write the axial tune in terms of the average field index, flutter, and the instantaneous edge angle: νz² = -k + F(1+tan²ξ) The radial tune is written as νr² = 1+k … The rms variation of the vertical field is called flutter and is denoted as F. The azimuthal magnetic field component, Bθ, is also proportional to the flutter.
Editorial note, tabletop extrapolation: The design equation for combining a weak-focusing taper with AVF sectors, showing the two contributions add. Two convention traps, both resolved here: (1) this paper calls F 'the rms variation' — for the linear-in-F tune formula to be the standard Livingood form, F must be the MEAN-SQUARE fractional variation ⟨((B−⟨B⟩)/⟨B⟩)²⟩, i.e. the square of the rms fraction, exactly as the same program's later paper defines it (F² there = ⟨…²⟩, tune quadratic in its F; dg-1746) — reconcile against Livingood before numeric use; (2) k = d ln⟨B⟩/d ln R is NEGATIVE for a falling field, opposite in sign to the magnet study's n, so n = −k. The tan²ξ factor is why edge angle is a powerful and dangerous knob — it grows without bound (dg-1697).
-
The Rutgers 12-inch neutron work used the d-d reaction, described by the author as having a broadly peaked cross section at a mere 180 keV, with the d(d,n)He3 reaction producing 2.45 MeV neutrons quasi-isotropically for an incident beam in the 180 keV regime. This is a DEUTERON beam on a deuterated target — not a proton reaction.
Source quote & editorial note
Many nuclear reactions produce neutrons, but perhaps the simplest is d-d reaction, with a broadly peaked cross section at a mere 180 keV. With an incident energy beam, in the regime of 180 keV, the reaction d(d,n)He3 reaction produces 2.45 MeV neutrons quasi-isotropically.
Koeth, Neutron Production with a 12-Inch Cyclotron (2017) — p. 1
Editorial note, tabletop extrapolation: The species distinction is the load-bearing radiological fact: deuterons on a deuterated target make fast neutrons at only ~180 keV, while common stable targets have (p,n) thresholds above 1 MeV (7Li(p,n) at ~1.88 MeV is among the lowest) — so a sub-MeV proton machine's neutron picture hinges on verifying the beam really is protons (deuterium contamination opens the D–D channel), what the beam actually strikes, and the truthful maximum energy; no blanket neutron-free claim follows. The cross-section characterization and the quasi-isotropic 2.45 MeV figure are the source's own: at finite beam energy the neutron energy is angle-dependent, and quantitative cross-section shapes should be taken from evaluated data at design time, not from this description. The source states its laboratory move was what provided the radiological controls to permit fast-neutron generation (abstract, p.1).
-
For spiral-edged AVF sectors the vertical tune obeys nu_z^2 = -k + F(1 + tan^2 xi), where F is the flutter (mean field variation at fixed radius), k the average negative field index, and xi the edge angle; the form is convenient for Archimedean spirals r = a*theta^(1/n), for which the Rutgers paper states tan xi = d(theta)/dr.
nu_z^2 = -k + F(1 + tan^2 xi); Archimedean spiral r = a*theta^(1/n); edge angle (from radial): tan xi = r*d(theta)/dr = n*theta [source prints tan xi = d(theta)/dr, which is not dimensionless — corrected 2026-09-05, site wave-18 audit; verify conventions against Livingood, the paper's ref 10, before numerical use]Source quote & editorial note
This form is convenient for sectors defined by an Archimedean spiral, r = aθ^(1/n), for which tan ξ = dθ/dr.
Editorial note, tabletop extrapolation: PDF p.3 = printed p.293 (the exponent 1/n is printed as a superscript; the quote transcribes it inline). The design equation for trading spiral tightness against vertical tune before cutting steel — with one correction applied: as printed, tan ξ = dθ/dr is not dimensionless; the standard edge-angle relation is tan ξ = r·dθ/dr, which for the stated Archimedean spiral evaluates to nθ. The flutter and approximation conventions are the paper's; verify against Livingood (its own ref [10]) before using the tune expression numerically.
-
A weak-focusing cyclotron only meets the cyclotron condition at one point in the ion's flight from source to target; the accumulated error is tolerable as long as the overall integrated phase slippage stays under 90 degrees, and raising the accelerating dee voltage reduces the number of turns and hence the accumulated slippage. Alternatively, starting the ions in a field that is too high lets the slippage run one way, meet the condition midway, then reverse to net zero.
Source quote & editorial note
This error is acceptable, as long as the overall integrated phase slippage is less than 90°.
Editorial note, tabletop extrapolation: PDF p.1 = printed p.299. The governing constraint for any non-isochronous tabletop machine, under the source's convention: keep the integrated phase slippage inside the source's 90-degree budget, remembering the whole phase TRAJECTORY matters — a net-zero final slip does not save a beam that left the accelerating window mid-flight. Low dee voltage hurts twice (more turns against the same budget), and the deliberate start-above-nominal-field trick is best read as centering the phase excursion, not as a free correction.
-
In the Rutgers 12-inch cyclotron's weak-focusing field the axial tune is nu_z = sqrt(n) and the radial tune nu_x = sqrt(1-n), with total transverse stability for 0 < n < 1; coupling resonances further exclude n = 0.2, 0.36 and 0.5 (and higher values).
n = -(r/B)(dB/dr); nu_z = sqrt(n); nu_x = sqrt(1-n)Source quote & editorial note
Values of n=0.2, 0.36, 0.5 (and others yet higher) need to be avoided.
Editorial note, tabletop extrapolation: PDF p.2 = printed p.370. The explicit forbidden-n list a weak-focusing pole-tip designer rarely sees written down: inside 0 < n < 1, the taper must also avoid 0.2 (Qx = 2Qz), 0.36 and 0.5. A well-chosen profile keeps n below 0.2 for nearly the whole acceleration (the source's own following prescription); the design questions are where n(r) crosses what, and how fast — compute or map n(r) rather than assuming which resonances are in play.
-
The axial betatron period Tz relates to the ion revolution period T0 by Tz = T0/sqrt(n), so it takes 1/sqrt(n) revolutions to complete one vertical betatron oscillation and the betatron phase advances by sqrt(n) of a period per revolution.
T_z = T_0 / sqrt(n)Source quote & editorial note
Thus for a given n it takes 1/√n ion revolutions to complete one vertical betatron oscillation
Editorial note, tabletop extrapolation: PDF p.2 = printed p.370. What converts a photograph into a number: count N revolutions between same-phase vertical maxima and νz ≈ 1/N — and, under the smooth azimuthally-symmetric weak-focusing approximation, n ≈ 1/N². It is an average over the interval, not a point measurement. On this machine's gently tapered poles νz ≈ 0.09 mid-radius (Fig. 3a), i.e. about 11 turns per oscillation, comfortably resolvable on its radial-draw images; treat that as Rutgers calibration context, not a class-typical value.
-
The author places the limit for H- acceleration at about 70 MeV: beyond it H- ions are stripped by the magnetic field (motional electric field), forcing a lower field and therefore a larger, more expensive cyclotron. [source prints 'magnetic dissipation'; the standard term is Lorentz/magnetic stripping (dissociation)]
Source quote & editorial note
The limitation for H¯ acceleration is the energy of about 70 MeV, as further acceleration would result in magnetic dissipation of H¯ ions, and require reduction of magnetic field
Karamyshev, A New Design of PET Cyclotron — IPAC2022, paper THPOMS016, JACoW (2022) — p. 3
Editorial note, tabletop extrapolation: The 70 MeV figure is the author's limit for this magnet regime. Lorentz stripping depends on energy and field together through the rest-frame electric field E′ = γβcB, so a tabletop H- design compares a cited stripping lifetime at its own E′ with the transit time rather than assuming the loss negligible; residual-gas stripping is a separate loss mechanism.
-
Treat the Red Book (ORNL-6086, Barnett 1990) cross sections as evaluated recommendations based almost exclusively on experiment, each carrying its own estimated uncertainty derived from the original sources' quoted errors and their mutual consistency; the compilation states these accuracies were estimated conservatively as upper limits, and for some reactions no uncertainty could be estimated at all.
Source quote & editorial note
In all cases, the estimated accuracies of the recommended data are believed to have been conservatively estimated, and are intended to represent upper limits.
Editorial note, tabletop extrapolation: A beam-survival estimate inherits the per-reaction accuracy line printed under each table (10 % to 45 % for the reactions relevant here, or "unknown"); carry that band through to the predicted survival fraction rather than quoting a single number. That the stated bands are upper limits is the compilation's own judgment.
-
For protons in residual H2 the charge-state loss process tabulated is electron capture, H+ + H2 -> H + H2+. Its cross section (Red Book A-28) peaks near 9.6e-16 cm2 at 4-7 keV/amu, falls to 2.5e-16 at 40 keV, 2.9e-17 at 100 keV and 1.8e-18 at 200 keV; the compilation assigns 10 % accuracy above 100 eV/amu and unknown accuracy below.
sigma_c(H+ + H2): 4 keV 9.33e-16, 7 keV 9.57e-16, 10 keV 8.66e-16, 20 keV 5.79e-16, 40 keV 2.50e-16, 70 keV 7.78e-17, 100 keV 2.91e-17, 200 keV 1.76e-18 cm2 (eV/amu = lab eV for H+)Source quote & editorial note
Accuracy: E < 100 eV/amu - Unknown; E > 100 eV/amu - 10%
Editorial note, tabletop extrapolation: Because the cross section drops about thirtyfold between 10 and 100 keV, the proton loss in a 100 keV-1 MeV machine is concentrated in the early, low-energy turns near the centre; a survival integral must therefore weight each turn by sigma at that turn's energy rather than use one cross section for the whole spiral. Raising the energy gain per turn shortens exactly the part of the path where sigma is largest.
-
Use the Chebyshev fits only inside the stated Emin-Emax range and check each reaction's printed fit-quality statement before replacing the table with the fit; for H+ + H2 electron capture the fit deviates by 15.5 % rms and up to 41.3 % (at 10 eV/amu), larger than the 10 % data accuracy, so interpolate the table in log-log where precision matters.
A-28 fit (PDF p. 44): Emin 2.6 eV/amu, Emax 4.0e6 eV/amu, rms 15.5 %, max 41.3 % at 10 eV/amu; A-58 H2+ + H2 fit (p. 74): rms 3.8 %, max 7.1 % at 8.1e4 eV/amu; G-46 fit (p. 629): rms 0.9 %, max 1.3 % at 1e4 eV/amuSource quote & editorial note
For convenience, least-squares Chebyshev polynomial fits to the recommended cross-sections and rate coefficients are presented for each reaction. These polynomial fits are indicated graphically, and a statement is made concerning the accuracy of each fit to the recommended data.
Editorial note, tabletop extrapolation: A loss calculator coded from the fitting coefficients is convenient, but the fit discrepancy is an additional approximation uncertainty on top of the data accuracy (the two are not simply additive); for the proton capture fit the largest deviations sit at eV energies irrelevant to the accelerated beam, so check the fit against a few table points over the machine's own 1–100 keV range before preferring either the fit or log-log interpolation.
-
Look up molecular ions at energy per amu, not laboratory energy: the Red Book tabulates every reaction in eV/amu (equal projectile velocity), so an H2+ at lab energy E is read at E/2 and an H3+ at E/3. For H2+ + H2 electron capture (A-58) the cross section is 5.5e-16 cm2 at 2 keV/amu, 3.4e-16 at 10 keV/amu and 4.6e-17 at 50 keV/amu, with stated accuracy 30 % below 1 keV/amu and 20 % above.
E_amu = E_lab / A (A = 1, 2, 3 for H+, H2+, H3+); A-58 sigma_c(H2+ + H2): 1e3 5.39e-16, 2e3 5.49e-16, 5e3 4.81e-16, 1e4 3.39e-16, 2e4 2.04e-16, 5e4 4.58e-17, 1e5 6.43e-18 cm2Source quote & editorial note
Accuracy: H2+ + H2 - 30% E < 1x10^3 eV/amu; 20% E > 1x10^3 eV/amu
Editorial note, tabletop extrapolation: In a fixed-field cyclotron an H2+ reaches half the lab energy of a proton at the same radius, i.e. a quarter of the energy per amu, so over the whole spiral it sits much nearer the 1-10 keV/amu region where loss cross sections are largest. Reading an H2+ table at lab energy instead of per-amu energy misplaces the ion by a factor of two in velocity-equivalent energy.
-
For H2+ in H2 use the total destruction cross section (Red Book G-46), not electron capture alone: total destruction, the sum of all reactions that destroy the molecular ion, stays near 8.1-8.5e-16 cm2 from 1.5 to 15 keV/amu and falls to 3.8e-16 at 50 keV/amu (20 % accuracy), exceeding capture alone by 1.5x at 2 keV/amu, 2.4x at 10 keV/amu and about 8x at 50 keV/amu.
G-46 Q_d(H2+ + H2): 1.5e3 8.12e-16, 2e3 8.18e-16, 4e3 8.30e-16, 7e3 8.15e-16, 1e4 8.18e-16, 1.5e4 8.46e-16, 2e4 7.75e-16, 4e4 4.69e-16, 5e4 3.81e-16 cm2 (eV/amu); ratio Q_d / sigma_c(A-58) = 1.5 (2e3), 2.4 (1e4), 3.8 (2e4), 8.3 (5e4)Source quote & editorial note
This cross section is the sum of cross sections for all reactions that destroy the H2+ molecular ion in passage through H2.
Editorial note, tabletop extrapolation: An H2+ beam disappears by dissociation as well as by capture, and dissociation grows in relative importance with energy; a survival estimate built on the capture table alone can understate the H2+ loss several-fold in the tens-of-keV/amu range. The same page notes that large variations in ion-source operating conditions changed the measured cross sections by up to 10 %.
-
H3+ in H2 is destroyed with a total cross section (Red Book G-44) that rises from 3.4e-16 cm2 at 1.6 keV/amu to 7.1e-16 at 15 keV/amu and is 5.6e-16 at 35 keV/amu (20 % accuracy); across 2-20 keV/amu this is smaller than the H2+ total destruction (G-46) at equal energy per amu (about half at 2 keV/amu, 80-90 % of it at 10-20 keV/amu), and the two are comparable near 30 keV/amu.
G-44 Q_d(H3+ + H2): 1.6e3 3.39e-16, 2e3 3.78e-16, 4e3 5.07e-16, 7e3 6.10e-16, 1e4 6.65e-16, 1.5e4 7.06e-16, 2e4 6.81e-16, 3e4 5.93e-16, 3.5e4 5.58e-16 cm2 (eV/amu); ratio G-44/G-46 = 0.46 (2e3), 0.61 (4e3), 0.81 (1e4), 0.88 (2e4)Source quote & editorial note
This cross section is the sum of cross sections for all reactions that destroy the H3+ molecular ion in passage through H2.
Editorial note, tabletop extrapolation: Species comparison for a survival calculation: at the same energy per amu the three hydrogen ions have loss cross sections of the same order (several 1e-16 cm2 at a few keV/amu). The species difference in survival therefore comes from each ion's energy-per-amu history along its spiral together with its path length and turn count — integrate n·Q_d(E/A) over each species' actual orbit rather than inferring survival from cross-section magnitude alone. The table covers only 1.6–35 keV/amu (single reference); outside that range the compilation gives no H3+ total.
-
The Red Book states that the H2+ electron-capture cross section depends on the vibrational level of the incident ion, with the effect growing as energy decreases; its v = 0-3 data (A-60, from Campbell et al., accuracy unknown) show v >= 1 about 1.4-1.6x v = 0 at 4 eV/amu, with the four levels still spread by about 25 % (5.97-7.43e-16 cm2) at 0.5 keV/amu, the top of the tabulated range.
A-60 (PDF p. 76) sigma_c(H2+(v) + H2) at 4 eV/amu: v=0 7.65e-16, v=1 1.20e-15, v=2 1.22e-15, v=3 1.10e-15 cm2; at 500 eV/amu: 7.20e-16, 7.43e-16, 5.97e-16, 6.59e-16 cm2Source quote & editorial note
It is well known that the H2+ electron capture cross sections are dependent on the vibrational levels of the incident H2+ ion. See the following set of data for the magnitude of this effect which increases as the energy decreases.
Editorial note, tabletop extrapolation: The vibrational-state data stop at 0.5 keV/amu, so for the accelerated beam the effect enters as an additional uncertainty on the H2+ loss rather than a correction; it matters most in the source and the first gap crossings, where H2+ energies per amu are lowest.
-
At thermal to few-eV energies H2+ in H2 converts to H3+ by H2+ + H2 -> H3+ + H with a very large cross section (Red Book H-10: 1.2e-14 cm2 at 0.01 eV/amu, 1.9e-15 at 1 eV/amu, 2.2e-17 at 7 eV/amu) and a thermal rate coefficient of (2.0 +/- 0.15)e-9 cm3/s; the compilation's accuracy for the cross section is +/- 45 %.
k_thermal = 2.0e-9 cm3/s ; conversion rate per H2+ = k * n_H2 ; H-10 sigma: 1e-2 1.16e-14, 1e-1 6.23e-15, 1 1.87e-15, 2 8.51e-16, 4 2.09e-16, 7 2.22e-17 cm2 (eV/amu)Source quote & editorial note
(6) Rate coefficient k at thermal energy for production of H3+ is: k = [2.0 ± .15]x10^-9 cm3/s
Editorial note, tabletop extrapolation: This reaction acts predominantly before acceleration — in the source plasma and extraction region, where H2 density is highest and H2+ is slow — shifting the extracted species mix toward H3+ as source pressure rises. The cross section falls steeply but the table ends at 7 eV/amu rather than at zero, so it is usually much less important after the first substantial acceleration; include it when estimating H2+ transmission through any dense, low-energy path.
-
Phelps (1990) states that it did not attempt to assign accuracy estimates to its recommended hydrogen-ion cross sections, that it indicated areas of uncertainty and where interpolations and extrapolations were made, and that every tabulated entry should be considered uncertain by at least +/-5 %, whatever the number of figures printed (all three statements in the introduction on the cited page).
Source quote & editorial note
Although some entries in the tables are given to several significant figures, all entries should be considered uncertain to at least ± 5%.
Editorial note, tabletop extrapolation: Treat +/-5 % as a floor, not an error bar. Where Phelps and the Red Book overlap (H+, H2+, H3+ in H2 up to 10 keV laboratory energy) the spread between the two compilations is a useful cross-check; Phelps states on the same page that it made extensive use of the Barnett compilations, so agreement between them is not fully independent confirmation.
-
Convert units before mixing compilations: Phelps tabulates laboratory energies of the projectile (target H2 at rest) and cross sections in units of 1e-20 m2 (= 1e-16 cm2), whereas the Red Book uses eV/amu and cm2; a 10 keV H2+ in Phelps is a 5 keV/amu entry in the Red Book.
1e-20 m2 = 1e-16 cm2 ; E_amu = E_lab/A ; target H2 at restSource quote & editorial note
Unless otherwise specified, all energies are laboratory energies EL rather than relative, center-of-mass, or "collision" energies.
Editorial note, tabletop extrapolation: With the target taken at rest, the cross section applies directly with the beam speed: a keV ion moves hundreds of times faster than room-temperature H2, so a gas-kinetic sqrt(2) relative-speed factor (appropriate when both partners are thermal) does not belong in the beam mean-free-path formula.
-
For H+ + H2 charge transfer (H2+ product, i.e. fast H formation), Phelps gives 4.5e-20 m2 at 1 keV, 6.8e-20 at 1.78 keV and 8.1e-20 m2 at 10 keV laboratory energy, and states that between 2 and 10 keV its values approach the fast-H production tabulated by Barnett et al. (the 1977 ORNL compilation); against the 1990 Red Book A-28 the 10 keV values agree within ~7 % (8.1 vs 8.66e-16 cm2), inside the Red Book's 10 % band.
Phelps Table 1 (PDF p. 5), H2+ product column: 562 eV 2.3, 1 keV 4.5, 1.78 keV 6.8, 3.16 keV 8.2, 5.62 keV 8.4, 10 keV 8.1 (1e-20 m2); the 7 % comparison with ORNL-6086 A-28 is the reviewer's, not Phelps'sSource quote & editorial note
For energies from 2 to 10 keV, our recommended values approach the tabulated values for fast H production from Barnett et al.
Editorial note, tabletop extrapolation: Proton capture near the cross-section peak is a well-determined input to a survival estimate, with the two compilations agreeing inside the Red Book's 10 % band (and not fully independently, see dg-1916). The hydrogen density along the path — gauge gas correction, pressure gradients — is an additional uncertainty to be quantified for the machine alongside the cross-section band, not ranked above or below it without an uncertainty budget.
-
For H2+ + H2 the two compilations split capture and dissociation differently and agree less closely on the total: at 10 keV lab (5 keV/amu) Phelps gives charge transfer (slow H2+) 7.0e-20 m2 plus dissociation to H+ + H 2.44e-20 m2 (9.4e-16 cm2 combined), against Red Book capture 4.8e-16 cm2 and total destruction ~8.2e-16 cm2 (log-log interpolation, not a printed value); Phelps states, citing Latimer et al., that the fast-H2+ destruction cross section is only slightly larger than charge transfer over 1-20 keV.
Phelps Table 4 (PDF p. 11) at 10 keV lab, 1e-20 m2: slow H2+ 7.0, v=0-1 0.215, H+ + H 2.44, Ly-a 0.72, H-a 0.157, e 1.27; Red Book 5e3 eV/amu: A-58 4.81e-16 (printed), G-46 ~8.2e-16 cm2 (reviewer's log-log interpolation between printed 4e3 8.30e-16 and 7e3 8.15e-16; 8.24e-16)Source quote & editorial note
The cross section for the destruction of fast H2+ is slightly larger than that shown for charge transfer at energies from 1 to 20 keV.
Editorial note, tabletop extrapolation: For a survival estimate use a total destruction cross section rather than a single channel: the capture-only numbers differ by ~45 % between the two sources (7.0e-20 m2 vs 4.8e-16 cm2), outside the Red Book's 20 % band, while the totals (9.4 vs ~8.2e-16 cm2) differ by ~15 %, inside it. Phelps's recommended table ends at 10 keV laboratory energy (the destruction measurements it cites run to 20 keV); the Red Book tabulates its total through 50 keV/amu.
-
For slow H2+ in H2 the dominant process is H3+ formation, not charge transfer: Phelps lists H3+ + H at 78e-20 m2 at 0.1 eV against 0.28e-20 for charge transfer (slow H2+), with the two crossing between 3.16 and 4.22 eV laboratory energy; above that charge transfer dominates, peaking at 11.0e-20 m2 near 7.5 eV, while H3+ formation falls to 0.42e-20 at 10 eV.
Phelps Table 4 (PDF p. 11, 1e-20 m2) H3+ / slow H2+: 0.1 eV 78 / 0.28; 1 eV 22.7 / 1.53; 3.16 eV 8.4 / 5.9; 4.22 eV 5.4 / 8.5; 7.5 eV 1.3 / 11.0; 10 eV 0.42 / 10.8; 17.78 eV 0.015 / 10.2Source quote & editorial note
The dominant cross section for low-energy H2+ in H2 is that for the formation of H3+ + H.
Editorial note, tabletop extrapolation: H2+ ions below roughly 3–4 eV in the source or extraction gap convert mainly to H3+; above the crossover charge transfer dominates and H3+ formation falls to about 0.15 % of the charge-transfer value by 18 eV in the table (0.015e-20 vs 10.2e-20 m2). Its absolute probability after acceleration follows from that cross section times the gas column the ions traverse, and the loss physics of the preceding rules takes over.
-
Do not simply add product-channel cross sections to obtain a molecular-ion loss cross section: Phelps tabulates H3+ + H2 by product (fast H+, fast H2+, fast H, fast H2, ...), and reports that the total destruction cross section measured by Williams, Geddes and Gilbody above 2.5 keV is approximately an appropriately weighted sum of those channels; use a total destruction table for survival.
Phelps Table 6 at 10 keV lab (3.3 keV/amu), 1e-20 m2: fast H+ 1.27, fast H2+ 1.17, fast H 7.8, fast H2 3.75; Red Book G-44 total at 3.3 keV/amu ~4.7e-16 cm2 (reviewer's log-log interpolation between printed 2e3 3.78e-16 and 4e3 5.07e-16; 4.67e-16)Source quote & editorial note
for EL > 2.5 keV and is approximately equal to an appropriately weighted sum of the reaction cross sections shown.
Editorial note, tabletop extrapolation: A single dissociative event yields several fast fragments, so unweighted product cross sections overcount destruction; the fast-H column alone (7.8e-20 m2 at 10 keV) already exceeds the Red Book total (~4.7e-16 cm2 at the same velocity). For H3+ survival use the Red Book G-44 total.
-
Collisional dissociation of H2+ depends strongly on the ion's vibrational excitation (data of Guyon et al. at 8-32 eV and of Lindsay, Yousif and Latimer at 1 keV, as reviewed by Phelps); Guyon et al. suggest this explains why low-energy H+ formation cross sections measured with H2+ made by 80-eV electron impact are much larger than those from H2+ formed in known vibrational levels.
Source quote & editorial note
Data showing a strong dependence of the cross section for collisional dissociation of H2+ on the degree of vibrational excitation have been obtained
Editorial note, tabletop extrapolation: A filament or discharge source produces H2+ in a spread of vibrational states, so the dissociation part of the H2+ loss in a working machine may differ from values measured with state-selected ions; carry this as an unquantified state-distribution uncertainty on H2+ survival — its direction depends on which states the source populates relative to the reference data. The Red Book's separate note that ion-source operating conditions shifted its measured H2+ totals by up to 10 % (G-46, p. 629, carried on dg-1912) is the only scale indication in these sources.
-
Exposure in mR overweights soft photons: Fewell and Shuping note that comparing tubes by mR/mAs is not fair to the W anode because its spectrum has fewer low-energy photons, which are more efficient at ionizing air, and they supply photon-fluence conversions for that reason.
Source quote & editorial note
Furthermore, an intensity comparison in terms of mR/mAs is really not fair to the W anode, because its spectrum has fewer of the low energy photons which are more efficient at ionizing air.
Editorial note, tabletop extrapolation: Exposure in roentgens weights soft photons heavily — the point of the authors' remark. For shielding, the quantity that matters behind a wall is the hard tail of the actual spectrum: use a measured or conservatively modelled source spectrum and transport it through the real wall, including hard-tail, fluorescence and scatter contributions. Unfiltered tube output is not a conservative bound for a chamber, and filtered-beam data are not automatically realistic for it.
-
Gas focusing works by ion space charge: positive ions made by beam collisions drift out only at thermal speed and accumulate in the beam. Johnson's estimate for a 1 degree initial divergence at 2 x 10^-5 A is a radial field of about 1 V/cm, supplied if there are about four positive ions per beam electron; a crossed-beam test in his tube measured a deflection equivalent to about 3 V/cm.
E_r(required) ~ 1 V/cm for 1 deg divergence, I = 2e-5 A; ~4 ions per beam electron; measured ~3 V/cm (crossed beams)Source quote & editorial note
This field strength is produced, with beams of the ordinary intensity, if there are four positive ions for each electron in the stream, a condition which seems not unreasonable.
Johnson, A Low Voltage Cathode Ray Oscillograph — Bell System Technical Journal 1(2), 142–153 (1922) — p. 7
Editorial note, tabletop extrapolation: Johnson's four-ions-per-electron figure belongs to his assumed beam (1 degree divergence, 2 x 10^-5 A) and stays there. The transferable point is that positive-ion compensation can oppose electron space-charge spreading; how much compensation a given beamline needs is found experimentally or from a geometry- and species-specific ionization and transport model, since the ion yield per electron depends on energy, gas, pressure and loss mechanisms rather than falling monotonically with beam speed.
-
A 26 keV electron beam crossing a low-density helium flow in the NASA Langley 22-inch helium tunnel produced visible fluorescence whose brightness the authors report as dependent on the local gas density; the free-stream number densities in that study were 6.4 x 10^23 to 14.3 x 10^23 per m^3 at 2.2-2.5 K.
n = 6.4e23 - 14.3e23 m^-3 (p.6, image-read); p = n k T -> 2.65-5.92 kPa at 300 K, and roughly 19-49 Pa at the tunnel's own 2.2-2.5 K (referee arithmetic, corrected 2026-10-25: the first edition of this card mis-stated the room-temperature equivalent by a factor of 1,000)Source quote & editorial note
It has been found that the electron beam efficiently excites the test gas to produce visible fluorescence throughout the flow field with a brightness which is dependent upon the local gas density.
Editorial note, tabletop extrapolation: Beam-induced emission scales with the number of target atoms per unit path, so beam brightness in a fill gas is set by number density (pressure and temperature together), not by pressure alone. The tunnel's number densities correspond to 2.65–5.92 kPa at room temperature — three decades above a few-pascal visible-beam fill — so the Langley brightness regime is not that of a room-temperature bench fill, and the brightness–density relation reported is the source's for its own regime. The measurement is in helium; the source has no neon data.
-
Trace nitrogen can dominate the visible emission from an electron beam in helium: with 40-400 ppm N2 in the Langley helium tunnel, N2+ first-negative emission appeared strongly at every observation point and was the main radiator far from the beam (upstream and downstream); its intensity was too high for direct beam excitation (the He 5016 A to N2+ (0,0) excitation cross-section ratio is about 2 at 26 keV), and the authors cite two suggested mechanisms, transfer from helium metastables or from He2+.
N2 impurity 40-400 ppm; sigma(He 5016 A)/sigma(N2+ (0,0)) ~ 2 at 26 keVSource quote & editorial note
the intensity of the molecular nitrogen radiation observed is too high to be attributed to direct excitation by the electron beam.
Editorial note, tabletop extrapolation: Under the Langley conditions (26 keV beam, cryogenic flowing helium), 40–400 ppm N2 produced strong N2+ emission and was the main radiator away from the beam, apparently through an indirect helium-mediated process the authors suggest but do not establish; atomic helium stayed the primary radiator within the beam. A spectrometer check of the emission is a sensitive impurity monitor whose sensitivity and concentration scale have to be calibrated for the actual fill and beam; the suggested transfer reactions are helium-specific, so any impurity effect in neon must be measured separately.
-
Where the light comes from depends on distance from the beam: in the Langley study atomic helium was the primary radiator inside the electron beam, molecular helium just downstream, and N2+ far from the beam; atomic-helium lines extended several beam-core diameters because of electron scattering and diffusion of resonance radiation, and He2 emission persisted for more than 500 microseconds downstream.
He2 afterglow > 500 microseconds (p.9); halo ~ several core diameters (p.9)Source quote & editorial note
(3)Atomic helium was identified as the primary radiator in the electron beam, whereas molecular helium was the primary radiator near the beam.
Editorial note, tabletop extrapolation: The apparent width of a luminous beam can exceed the electron core: in the Langley helium flow, scattered electrons and trapped resonance radiation spread atomic-line emission over several core diameters, and afterglow species drifted with the gas flow. Size a visible beam from the bright core, not the halo; the halo extent in neon is not measured in this source.
-
Brode's absorption coefficient a is the total effective electron-loss area of all molecules in 1 cm^3 of gas at 1 mm Hg, defined by I = I0 exp(-a x p) with x in cm and p in mm Hg; its reciprocal at a given pressure is the mean free path for removal from the beam, and dividing a by 3.65 x 10^16 (molecules per cm^3 at 1 mm Hg) gives the cross-section per molecule.
I = I0 exp(-a x p); lambda = 1/(a p) [cm, p in mm Hg]; sigma = a / 3.65e16 cm^2. Image-read: digits 3.65 are sharp, exponent blurred but 16 is the only physical reading (ideal gas at 1 mm Hg: 3.54e16 cm^-3 at 0 C, 3.30e16 at 20 C); Brode's printed 3.65e16 is ~3% above the 0 C value (his constant as printed, possibly a 3.56 transposition - unconfirmed).Source quote & editorial note
It is the total effective absorbing area in square centimeters of all the molecules in a cubic centimeter of gas at a pressure of 1 mm of Hg
Brode, The Absorption Coefficient for Slow Electrons in Gases — Physical Review 25, 636 (1925) — p. 5
Editorial note, tabletop extrapolation: Converts Brode's curves into beam path lengths at the reference temperature of his constant: at 1 Pa (7.5 x 10^-3 mm Hg), lambda (cm) = 1/(a x 0.0075). Brode's 'absorption' counts any electron deflected out of the narrow beam past a finite collector, so it is an effective beam-removal (attenuation) length — not a total-scattering cross-section, since small-angle forward scattering within the collector's acceptance is not counted, and not an energy-loss length. For another gas temperature convert a to a cross-section per molecule and use lambda = 1/(n sigma).
-
Argon is strongly energy-selective for slow electrons: Brode's coefficient peaks sharply at about 12 V and the 2 V value is one-fifth of the peak; his Fig. 5 peak is about a = 55 cm^-1 mm Hg^-1 (Ramsauer's dotted curve peaks higher), falling to about 7 at 360 V (figure readings).
Ar: a_max ~ 55 at ~12 V; a(360 V) ~ 7 (Fig. 5, p.7, figure reading); lambda(1 Pa) = 1/(a x 0.0075) cmSource quote & editorial note
In argon, the coefficient reaches a sharp maximum at about 12 volts, the value for 2 volts being only one-fifth of the maximum value.
Brode, The Absorption Coefficient for Slow Electrons in Gases — Physical Review 25, 636 (1925) — p. 2
Editorial note, tabletop extrapolation: Referee arithmetic at 1 Pa argon: lambda ~ 2.4 cm near 12 eV and ~ 19 cm at 360 eV. At fixed pressure, operating where Brode's removal coefficient is lower lengthens the unscattered path; whether the beam is visible over that path is a separate question of excitation and emission cross sections, beam current and viewing geometry.
-
Helium has a broad maximum at about 4 V and is the most transparent of Brode's gases at tens to hundreds of volts: his Fig. 5 peak is about a = 22, falling to roughly 2-3 at 360 V (figure readings); methane peaks near 7.5 V (about 70).
He: a_max ~ 22 at 4 V; a(360 V) ~ 2-3 (Fig. 5, p.7, figure reading)Source quote & editorial note
In helium there is a maximum at 4 volts. Methane gave a curve similar to these, with a maximum at about 7.5 volts
Brode, The Absorption Coefficient for Slow Electrons in Gases — Physical Review 25, 636 (1925) — p. 2
Editorial note, tabletop extrapolation: Referee arithmetic at 1 Pa helium: lambda ~ 6 cm near 4 eV and roughly half a metre at a few hundred eV. Brode measured no neon; his introduction only cites Akesson (two maxima in Ne) and Ramsauer and Mayer (Ne coefficient nearly constant, decreasing slightly with velocity), so a neon mean free path must come from another source.
-
Molecular gases behave differently: Brode found hydrogen has no maximum (a rises steadily as energy falls to 2 V, about 50 there, about 10 near 36 V and about 6 near 170 V by figure reading), while nitrogen and carbon monoxide are nearly identical with a maximum near 18 V and a minimum near 9 V; he cautions that those extremum positions were not determined accurately.
H2: a ~ 50 (2 V), ~10 (36 V), ~6 (170 V); N2/CO: max ~18 V, min ~9 V, a ~ 30 (Fig. 6, p.8, figure reading)Source quote & editorial note
The curves for nitrogen and carbon monoxide are nearly identical, both having a maximum at about 18 volts and a minimum at about 9 volts. The curve for hydrogen shows no maximum, the coefficient increasing steadily down to 2 volts.
Brode, The Absorption Coefficient for Slow Electrons in Gases — Physical Review 25, 636 (1925) — p. 2
Editorial note, tabletop extrapolation: Residual air (N2) or hydrogen in a noble-gas fill adds beam loss below a few hundred volts in proportion to its partial pressure; in pure nitrogen at 1 Pa (a ~ 30) the removal length is about 4-5 cm at 10-40 eV (referee arithmetic). Keep air leaks and outgassing hydrogen low if the beam must survive tens of centimetres.
-
Berger and Seltzer's electron tables run from 10 keV upward and are least certain at the low end: when the electron velocity falls to about the atomic-electron velocity their stopping-power formula ceases to be valid, so the range integral is closed by linear interpolation to zero below 1 keV; their 10 keV ranges exceed Nelms' (5 keV cut) by 12% for Pb, 6% Ag, 2% Al and 1% H2, with differences of at most 0.5% at 50 keV.
range r = integral_0^T [ (1/rho) dE/dx ]_tot^-1 dT (c.s.d.a.)Source quote & editorial note
The ranges of Nelms are on the low side, ours on the high side, and the discrepancies give an indication of the error likely to result from lack of knowledge of low-energy stopping behavior.
Editorial note, tabletop extrapolation: Use these tables for the keV end of a low-energy electron beamline (ranges, stopping, bremsstrahlung yield), not for the tens-to-hundreds-of-eV regime they do not cover. The 10 keV range discrepancy the authors quantify is for H2, Al, Ag and Pb (1–12 %); for neon, argon and air it is unquantified here — compare against modern ESTAR-class data before assigning a percentage.
-
SP-3012 accuracy, as the authors estimate it: measured stopping powers carry 4 to 10% uncertainty and are presumably less accurate than the calculation; the dominant error is the mean excitation energy (up to 10% for very small Z), and the collision stopping power is expected to be good to 1-2%; bremsstrahlung cross-sections are quoted at 20% up to 2 MeV.
error(collision dE/dx) <= 1-2%; I_adj up to 10% for very small Z; bremsstrahlung +/-20% below 2 MeV (p.26)Source quote & editorial note
It is difficult to estimate very closely the error of the data in Table I.
Editorial note, tabletop extrapolation: The authors estimate roughly 1–2 % for their calculated collision stopping powers; that does not establish a few-percent range accuracy for neon, argon or air, because ranges also carry the explicitly uncertain low-energy closure near the 10 keV lower limit — check ranges there against modern data. The radiation-yield column at 10–20 keV inherits the ±20 % bremsstrahlung uncertainty and is an estimate.
-
SP-3012 Table I, neon (page image read; table OCR unreliable): at 10 keV the collision stopping power is 17.91 MeV cm^2/g, total 17.92, c.s.d.a. range 3.226 x 10^-4 g/cm^2 and radiation yield 2.19 x 10^-4; at 20 keV the range is 1.078 x 10^-3 g/cm^2 and stopping power 10.62.
rho_Ne(1 Pa, 293 K) = pM/RT ~ 8.3e-9 g/cm^3; L = r/rho; dE/dx = S rhoSource quote & editorial note
0.010 1.791E 01 6.784E-03 1.792E 01 3.226E-04 2.190E-04
Editorial note, tabletop extrapolation: Referee arithmetic for a neon fill at 1 Pa and 293 K (density ~ 8.3 x 10^-9 g/cm^3): a 10 keV electron loses only ~0.15 eV per cm and its path length to stop is ~390 m, so a 10 keV beam crosses a bench-scale vessel essentially without energy loss and its light is a thin-target signal; at 100 Pa the range is still ~4 m.
-
SP-3012 Table I, argon (page image read): at 10 keV the collision stopping power is 14.57 MeV cm^2/g, c.s.d.a. range 4.050 x 10^-4 g/cm^2 and radiation yield 4.425 x 10^-4 (twice neon's); at 20 keV the range is 1.324 x 10^-3 g/cm^2.
rho_Ar(1 Pa, 293 K) ~ 1.64e-8 g/cm^3Source quote & editorial note
0.010 1.457E 01 1.112E-02 1.458E 01 4.050E-04 4.425E-04
Editorial note, tabletop extrapolation: Referee arithmetic for argon at 1 Pa and 293 K (~1.64 x 10^-8 g/cm^3): 10 keV path length ~250 m, energy loss ~0.24 eV per cm. The higher radiation yield than neon matters only for X-ray estimates; about 4 x 10^-4 of a stopped 10 keV beam's energy goes to bremsstrahlung in argon itself.
-
SP-3012 Table I, air (page image read): at 10 keV the collision stopping power is 19.70 MeV cm^2/g, c.s.d.a. range 2.892 x 10^-4 g/cm^2 and radiation yield 1.463 x 10^-4; at 20 keV the range is 9.805 x 10^-4 g/cm^2.
rho_air(1 atm, 20 C) = 1.205e-3 g/cm^3Source quote & editorial note
0.010 1.970E 01 5.012E-03 1.971E 01 2.892E-04 1.463E-04
Editorial note, tabletop extrapolation: Referee arithmetic at atmospheric density (1.205 x 10^-3 g/cm^3): 10 keV electrons travel about 2.4 mm of path in air and 20 keV about 8.1 mm (c.s.d.a. path, so projected penetration is less), so a beam vented or escaping through a thin window stops within millimetres of air. The 1.46 x 10^-4 radiation yield is for electrons stopped in air only; bremsstrahlung from the exit window, the target or any other higher-Z material intercepting the beam is assessed separately, and the 20 keV range is the table's own row (page image).
-
In Woo's account, multipacting (two-surface secondary-electron resonance) relies only on electrode secondary emission, so it does not depend on the gas species, and it is independent of pressure once the electron mean free path exceeds the gap; his scaling analysis makes the breakdown voltage a function of fd alone for a fixed line impedance (double the size and halve the frequency, same voltage).
V_multipactor = F(fd) for fixed Z0 (or parallel plates); valid when mean free path > dSource quote & editorial note
multipacting relies only on secondary electron emission from the electrodes. Thus, it is not dependent on the type of gas in which the breakdown takes place.
Editorial note, tabletop extrapolation: Within Woo's collisionless, field-free two-surface model the gas species drops out and fd is the scaling parameter, so results from other machines scale by fd for a fixed impedance rather than by voltage alone. Whether a given gap actually multipacts still depends on the RF voltage band and phase, geometry, the electrodes' secondary-emission yield and surface condition, the pressure-dependent collision probability near the mean-free-path limit, and any magnetic or dc field — a gap inside a cyclotron magnet needs its own check (dg-1986, dg-1988).
-
Woo's fd-p/f plane separates the regimes: the mean-free-path limit (electron mean free path = d) divides multipacting (lower p/f) from ionization breakdown (higher p/f); multipacting is cut off below fd = 70 MHz-cm; the collision-frequency transition sits at p/f = 10^-3 torr/MHz; and below roughly 100 MHz-cm the electron-oscillation-amplitude limit raises breakdown. The lines indicate transitions, not abrupt changes, and the underlying data are for air.
mean-free-path limit: (p/f)(fd) = p d ~ 0.025 torr-cm (air, Fig. 7 reading); cutoff fd = 70 MHz-cm; collision-frequency transition p/f = 1e-3 torr/MHzSource quote & editorial note
The mean free path limit represents the condition that the electron mean free path is equal to d. The mean free path limit, therefore, separates the ionization and multipacting breakdown regions.
Editorial note, tabletop extrapolation: Referee arithmetic from Fig. 7 (p.19): the air mean-free-path line follows p x d ~ 0.025 torr-cm, so at 1 Pa (7.5 x 10^-3 torr) a gap below ~3 cm is on the collisionless side and at 0.1 Pa almost any bench gap is — a collision-regime estimate only. fd below 70 MHz-cm excludes Woo's field-free two-surface mode (a 3 cm gap at 9 MHz is fd ~ 27); fd above it merely permits multipacting, which also needs a resonant voltage band and sufficient secondary yield, and a cyclotron's magnetic field changes the trajectories and allows other modes. A different fill gas shifts the line with its electron mean free path; the report places multipacting generally below 10^-2 torr (p.9).
-
For diffusion-controlled RF ionization breakdown (fd > 100 MHz-cm) the minimum breakdown power falls at p/f of about 10^-3 torr/MHz, i.e. at a pressure in microns of mercury equal to the frequency in MHz; below that minimum, reduced collisional energy absorption and fewer ionizing collisions (with wall losses where the gap is finite) raise the breakdown level.
p_min [micron Hg] ~ f [MHz] (p/f ~ 1e-3 torr/MHz)Source quote & editorial note
This is equivalent to saying that each minimum occurs at a pressure given in microns (1 micron = 10^-3 torr) equal to the frequency given in MHz.
Editorial note, tabletop extrapolation: Referee arithmetic: the diffusion-controlled minimum would fall near 9 microns (~1.2 Pa) at 9 MHz and near 27 microns (~3.6 Pa) at 27 MHz, so a 0.1-1 Pa fill sits at or below it, on the side where ignition gets harder as pressure falls. The rule holds for fd > 100 MHz-cm (a gap above ~11 cm at 9 MHz); bench gaps at 9 MHz are fd < 100, where Woo's oscillation-amplitude limit raises breakdown further. Air data in coaxial lines.
-
Korolov & Donkó (2015) measured the DC Paschen curve of hydrogen between 7.5 cm stainless-steel electrodes 1 cm apart: minimum breakdown at (pL)* = 1.2 Torr cm with V_BR,min = 290 V; deuterium matched at high pressure but its minimum was about 25 V higher and its low-pL breakdown voltage significantly higher.
(pL)* = 1.2 Torr cm, V_min = 290 V (H2, stainless, L = 1 cm)Source quote & editorial note
The present measurement gives a location of the Paschen minimum at (pL)* = 1.2 Torr cm with VBR,min = 290 V
Editorial note, tabletop extrapolation: Referee arithmetic: a 0.1–1 Pa hydrogen fill across a few-cm gap is pL ~ 0.002–0.02 Torr cm, two to three decades left of this minimum, where the source gives no breakdown data. The Paschen/Townsend curve describes gas breakdown only — field emission, surface flashover and particulates can break down an electrode assembly well below any extrapolated gas-breakdown voltage — so it is not a hold-off or electrical-safety guarantee; test the actual assembly conservatively. Electrode material and cleanliness move the minimum and especially the low-pressure branch.
-
At 13.56 MHz the hydrogen RF breakdown curve of Korolov & Donkó (2015) lies below their DC curve and is multi-valued: its minimum (point A, about 160 V by figure reading, V_BR being the peak of the electrode potential difference) is near pL = 3 Torr cm; past turning point B the breakdown voltage increases suddenly up to turning point C (about 3 Torr cm, 350 V), beyond which it increases smoothly as pL falls; where two breakdown voltages exist at fixed pressure they measured by slowly raising pressure at constant RF voltage.
V_BR = 2 U0 (peak difference of push-pull drive); f = 13.56 MHz; L = 1 cmSource quote & editorial note
we find a turning point (point "B" in figure 4(a)) beyond which we experience a sudden increase of the breakdown voltage
Editorial note, tabletop extrapolation: In this 13.56 MHz, 1 cm hydrogen experiment the low-pL branch is multi-valued: ignition at fixed pressure under rising voltage and at fixed voltage under rising pressure can occur at different points. In comparable systems characterise ignition with both sweeps rather than assuming a unique threshold.
-
The RF breakdown shape is set by electron oscillation amplitude versus gap (Korolov & Donkó 2015): when the amplitude is smaller than the gap, charge reproduction is purely gas-phase; when it becomes comparable (lower pressure or lower frequency) electrons reach the electrodes and surface emission takes over; Githens' 11 MHz data shift turning point B to higher pL.
Source quote & editorial note
at lower frequency the electron oscillation amplitude is higher at otherwise same conditions, therefore the bending point "B" moves to higher pL at lower f
Editorial note, tabletop extrapolation: Referee inference, limited: Korolov & Donkó tie the shape of the RF curve to the electron oscillation amplitude relative to the gap, which depends on the RF field, frequency, collision rate and geometry — and, in a cyclotron's 0.5–1 T field, on the strong magnetization of electrons at MHz frequencies, which the unmagnetized 1 cm bench case does not include. The surface-governed regime therefore cannot be assumed for a magnetized dee gap: calculate the electron motion for the actual field, frequency, gas and gap, or test ignition in the magnetized gap. The paper's data are for hydrogen at L = 1 cm and pL of roughly 1–30 Torr cm.
-
RF breakdown data disagree strongly between laboratories (Korolov & Donkó 2015): for Lisovskiy & Yegorenkov's 13.56 MHz, 2 cm stainless-steel data set the authors report a factor of 1.8 (71 V vs 125 V) between the minimum breakdown voltages, which they attribute to systematic error or impurities in those measurements; their own model reproduced their hydrogen curve with an effective secondary yield of 0.008 and electron reflection 0.2, after cleaning electrodes with ~8 mA DC discharges of both polarities for about 10 minutes.
gamma* = 0.008, rho = 0.2 (text p.5; Fig. 4(a) legend 0.008, Fig. 4(b) legend prints 0.08, p.6); DC ramp 0.5-1.5 V/s (p.2)Source quote & editorial note
a factor of 1.8 is found (71 V vs. 125 V) between the minimum breakdown voltage values.
Editorial note, tabletop extrapolation: Two careful measurements of the same nominal system differed by a factor of about 1.8, which the authors attribute to systematic error or impurities in the other data set — the reason to characterise ignition in the actual cell, not a general uncertainty cap. The authors conditioned their stainless electrodes with ~8 mA dc discharges of both polarities for about 10 minutes (p. 2) — their preparation, with its own electrical and surface-modification hazards, not a universal prescription. The yield is quoted as 0.008 in the text and the Fig. 4(a) legend; the Fig. 4(b) legend prints 0.08.
-
Size test pulses to the resonator's fill time: the Siemens team used 50 µs pulses because the fill time of their Q = 2600, 72 MHz test resonator was nearly 10 µs, and the authors judged longer pulses to make no sense for this transient test.
Referee arithmetic: T_f = Q_L/(pi f) = 2600/(pi x 72.5e6 Hz) = 1.14e-5 s = 11.4 us (Q dimensionless, f in Hz -> s); consistent with the paper's 'nearly 10 us' if Q = 2600 is the loaded Q (the paper does not say loaded or unloaded). Example: Q_L = 62 at 9.1 MHz -> 62/(pi x 9.1e6) = 2.2 us; Q_L = 1000 at 9 MHz -> 35 us.Source quote & editorial note
50 us pulses used at these measurements. Use of longer pulses makes no sense because the fill time of the cavity is nearly 10 us.
Editorial note, tabletop extrapolation: The relation transfers and favours low frequency: at fixed loaded Q the fill time grows as 1/f, so a 6-30 MHz dee has a longer settling time than a 72 MHz cavity of the same Q. A soft-start or pulsed test on a tabletop dee should use pulses several fill times long, computed from that dee's measured loaded Q and frequency.
-
Choose the transistor operating class knowing its efficiency ceiling: class A conducts continuously and at best approaches 50 % efficiency at full output; class B conducts only when driven and ideally approaches 78.5 %, so each device conducts for about half of each cycle — a single-ended stage relies on a tuned load to extract the fundamental, while push-pull combines complementary half-cycles; push-pull class B keeps the efficiency and strongly reduces even harmonics, and the accelerator SSA modules Jacob describes run push-pull class AB with a slight gate bias to further linearise and minimise harmonic distortion.
Ideal class B efficiency = pi/4 = 0.785 (78.5 %); class A <= 50 %.Source quote & editorial note
In class B, the amplifier only consumes power when it is driven by RF. This provides a much better efficiency, approaching 78.5% for an ideal transistor.
Editorial note, tabletop extrapolation: Transfers to a tabletop dee drive: push-pull class AB LDMOS is the standard compromise. A high-Q dee filters near its fundamental but cannot be assumed to suppress harmonics that land on other cavity or feed-line modes, so the output network and the modes deserve a survey. Class A wideband amplifiers (common as bench drivers) draw at least 2 W of DC per RF watt at full output (efficiency <= 50 %), against ~1.4 W for a 70 %-efficient class-AB stage.
-
Push-pull operation suppresses even harmonics because the pair's transfer characteristic is odd, Ids(-Vgs) = -Ids(Vgs), which minimises H2 and the higher even harmonics; the odd harmonics are not removed by this symmetry (a consequence of the odd characteristic, not stated separately by the source).
Odd transfer function -> output contains only odd harmonics (ideal symmetry); residual H2 from device mismatch.Source quote & editorial note
the odd characteristic Ids(−Vgs) = −Ids(Vgs) minimizes the second harmonic H2 and higher even harmonics.
Editorial note, tabletop extrapolation: Explains why a push-pull LDMOS stage at 6-30 MHz can be expected to show a suppressed second harmonic relative to the odd ones; the third harmonic of a 9 MHz drive (27 MHz) may fall where a dee resonator or its feed line has other modes, so it deserves attention in output-filter design. The ESRF specification (same source, p16) nonetheless puts H3 lower than H2 for that system, so measured levels must come from the actual amplifier.
-
A non-directional n-way combiner presents each module, on its own, a reflection coefficient (1-n)/n, close to full reflection for large n, and couples the inputs to each other by 1/n; power combines fully only by constructive interference, so all inputs must have equal amplitude and phase, and any differential-mode power returns to the module outputs.
b_j = [(1-n)/n + (n-1)/n] a0 = 0 for equal inputs (Eq. 2); |b_(n+1)|^2 = n |a0|^2 (Eq. 3). Wave amplitudes in sqrt(W).Source quote & editorial note
So, the power adds up correctly at the output of a combiner only by constructive interference, i.e. if and only if the signals at their inputs have the same amplitude and phase.
Editorial note, tabletop extrapolation: Relevant only when a tabletop drive sums two or more modules: phase-matched cables and equal module gain decide combining efficiency, and a non-directional combiner of the kind the source analyses has input-port coupling and no intrinsic isolation, so one failing module disturbs the others. Inspect the actual combiner's scattering matrix and isolation/load arrangement (hybrids obtain isolation through a terminated port) rather than judging by the presence of resistors. A single-module drive avoids combiner effects.
-
Operate a solid-state amplifier near its nominal power: the ESRF/ELTA SSA reaches 57 % DC-to-RF efficiency at its 150 kW nominal but only about 47 % at two-thirds power, so over-dimensioning the amplifier costs efficiency; drain-voltage modulation can recover efficiency at reduced power, but module stability must then be checked over the whole drain-bias range.
Fig. 15 (image-verified): eta = 57 % at 150 kW, ~47 % at 100 kW; plotted P_dc ~263 kW at 150 kW RF and ~212 kW at 100 kW RF. Referee cross-check: P_dc = 150/0.57 = 263 kW and 100/0.47 = 213 kW -> endpoint secant dP_rf/dP_dc ~1.0, whereas the text and the plotted eta_diff curve give 'about 90 %' — a rounding-sensitive check (e.g. 48 % at 100 kW gives 0.90) of a local derivative against a secant, not a demonstrated inconsistency.Source quote & editorial note
It is worth noting that the efficiency drops fast when operating the amplifier at reduced power, and is only about 47% at 100 kW, i.e. at ⅔ of the nominal power.
Editorial note, tabletop extrapolation: Transfers as a question to measure for a tabletop class-AB LDMOS stage: back-off efficiency depends on bias, matching, supply and device, so take the selected pallet's efficiency map — a dee needing ~150 W from a 1 kW-capable pallet (an assumed example) may run far below peak efficiency and dissipate correspondingly more heat per RF watt; sizing the device or its supply voltage to the actual requirement helps. The percentages are this 352 MHz system's.
-
A high-Q resonator driven by a pulse reflects as follows: the leading edge sees a short (r = -1); after the fill time T_f = Q0/(pi f (1+beta)) the reflection approaches r = (beta-1)/(beta+1); at the falling edge a further transient of amplitude 1 adds to the near-steady reflection, so the reflected wave can peak above 1.5 times the incident amplitude, more than twice the incident power.
ESRF example (Fig. 18, image-verified): Q0 = 32500, beta = 3.8, f = 352.2 MHz -> T_f = 32500/(pi x 352.2e6 Hz x 4.8) = 6.1e-6 s (6.1 us; Q dimensionless, f in Hz). r_ss = 2.8/4.8 = 0.58 (dimensionless). Fig. 18(a) pulse runs from -20 to +20 us (40 us); r(t) = r_ss - (1 + r_ss) e^(-t/T_f): at 40 us r = 0.583 - 1.583 x e^(-6.5) = 0.58, at 20 us r = 0.583 - 1.583 x e^(-3.3) = 0.52. The text's '0.53 amplitude of the almost steady reflection' fits a 20 us pulse, not the 40 us plotted pulse; a small source-internal discrepancy (0.53 vs 0.58), not resolved here. Peak |r| = 1.5 (figure label) -> power 1.5^2 = 2.25x incident.Source quote & editorial note
the total transient reflection therefore peaks at more than 1.5 times the incoming pulse amplitude, corresponding to more than twice the input power;
Editorial note, tabletop extrapolation: Transfers to any tabletop dee driven by a soft-start, pulsed or interlock-tripped drive, at resonance and after the cavity has filled: switching the drive off is itself a reflected-power event — the emitted amplitude at switch-off is 2β/(1+β)·(1−e^(−t_p/T_f)), about 1 at β = 1 and above 1 for an over-coupled resonator once filled, smaller for a short pulse — so trips and foldback must not rely on reflected power falling when the drive is removed, and the device must tolerate the switch-off transient. Evaluate the actual coupling, detuning, pulse length and interlock response time; at β = 1 on resonance the steady reflection is zero, off resonance it is not.
-
Yonts's Child's-law estimate for the calutron: with no gas the space-charge-limited emission at 150 V across 0.5 cm is about 19 mA/cm^2 (as printed; the printed inputs give ~17), i.e. of the order of 10 mA from the ~0.5 cm^2 of a 140-mil filament facing the slot, whereas the struck arc easily draws 5 A (~10 A/cm^2); the ion sheath at the filament is what makes arc currents of amperes possible, so ions must reach the filament. The report then solves for the sheath thickness and prints x = 6.7 x 10^-2 cm, "67 mm" and "about 7 mm" (decimal points lost in the typescript; evidently 0.67 mm and 0.7 mm) [referee: its x^2 step is ×10 high; the same inputs give x ≈ 2.1 x 10^-2 cm ≈ 0.2 mm].
j = 2.334e-6·V^1.5/d^2 [j in A/cm^2, V in volts, d in cm] (printed "2 334 x 10-6"; decimals lost in the typescript). Units pass: V = 150 V → V^1.5 = 1837 V^1.5; d^2 = 0.25 cm^2 → j = 2.334e-6 × 1837 / 0.25 = 0.0172 A/cm^2 = 17 mA/cm^2 (with the report's own V^1.5 ≈ 1.9e3 it is 17.7 mA/cm^2; the printed 19 mA/cm^2 is not reproduced from the printed inputs). × 0.5 cm^2 → ~9 mA, "of the order of 10 ma" as stated. Sheath: x^2 = 2.334e-6 × 1.9e3 / 10 = 4.4e-4 cm^2 → x = 2.1e-2 cm = 0.21 mm. The report prints x^2 = "2.334 x 10^-7 x 1.9 x 10^3 = 44 x 10^-4" (the product is 4.4 x 10^-4), then x = 6.7 x 10^-2 cm (= sqrt of 44e-4), "67 mm" (p.11) and "about 7 mm" (p.12): one ×10 slip in x^2 plus dropped decimal points. Corrected sheath ≈ 0.2 mm.Source quote & editorial note
The calculations below show that it is necessary for these ions to pass to the filament if arc currents of several amperes are to be obtained […] For V = 150 volts and d = 0.5 cm (usual calutron conditions), the emission current density is approximately 19 ma per cm² of the filament. Since the total area of a 140-mil filament facing the collimating surface slot is just over 0.5 cm², space-charge-limited current is of the order of 10 ma. When charge gas is present and the arc is struck and an ion sheath is formed, arc currents of 5 amps are easily obtainable. This is a current of 10 amp/cm² of effective filament surface. [Decimal points and the cm² superscripts are lost in the typescript scan; restored here.]
Editorial note, tabletop extrapolation: The same physics governs a small filament source: at the calutron's 150 V across 0.5 cm, space-charge-limited emission into vacuum is of the order of milliamperes, and a different gap and voltage must be put through Child's law on their own; drawing much more current than that needs ion neutralization near the filament. The ion flux that permits the high current is the flux that sputters the filament, so emission current and filament erosion are coupled.
-
Yonts's empirical erosion fit for the calutron filament (TiCl4 charge): weight loss per ampere-hour of arc current scales as (arc power)^(3/5), J/I = A·P^(3/5) with A = 0.36, a form suggested by Townes's cold-cathode sputtering law; in conclusion, full diode operation minimized wear but reduced output 'to the point of uselessness', while modified diode gave about normal output with useful life increased by a factor of at least 2 (ORNL-1758 p.38 of the report) [ORNL-1525: a factor of 2 to 3].
J/I = A·P^(3/5) with A = 0.36 (printed "0 36"), J in mg/h, I in A, P in W → J/I in mg/(A·h), A in mg/(A·h·W^0.6). Units check vs Table IV run 1: P = 798.8 W → P^0.6 = 55.1 (table 55.2) → 0.36 × 55.1 = 19.9 mg/(A·h); measured 95.5 mg/h / 5.03 A = 19.0 mg/(A·h) (table "Actual J/I" 18.98). Consistent to ~5 %. [Source-internal: Table IV's computed J/I column (20.9, 19.6, 16.1 …) implies a coefficient ≈ 0.37, not 0.36; its run-9 P^3/5 entry 38.7 does not match 635.1^0.6 = 48.0.]Source quote & editorial note
It has also been shown that operation of the full diode type reduces filament wear to a minimum but also reduces ion output to the point of uselessness
Editorial note, tabletop extrapolation: The fitted coefficient and exponent belong to that calutron, charge material and operating range, and the fit is for the normalized loss J/I, so total erosion J has no simple sub-linear law in P unless arc current is held fixed. The transferable content is the method — plot weight loss per ampere-hour against arc power and fit — and the finding that the operating mode that minimizes wear is not the one that maximizes useful output, so the choice is a compromise set by measurement.
-
Yonts contrasts arc-source practice with the gas-filled diode: there the tube runs on the nearly vertical part of the emission curve with an arc drop near the gas ionization potential, and filaments last a very long time provided the positive-ion energy stays below a critical value that depends on the gas (he cites 22 V for mercury-vapour diodes), which is why oxide-coated or thoriated-tungsten filaments are usable in diodes.
Source quote & editorial note
The filaments give extremely long life, provided the velocity of the positive ions does not exceed some critical value (22 volts in mercury vapor diodes) which depends on the type of gas in the tube.
Editorial note, tabletop extrapolation: The physics behind low-voltage operation as Yonts presents it: below some ion energy, cathode sputtering is small enough that low-work-function cathodes survive. The 22 V is his reported operating limit for mercury-vapour diodes, not a universal threshold: sputtering depends on projectile and cathode material, incidence angle, charge state and surface condition, so any hydrogen-source limit must be established for the particular ion species, cathode material and sheath conditions; low ion energy only tends to reduce sputtering, and an oxide cathode can also fail by poisoning or chemical attack.
-
Smyth (Princeton, 1925) concluded from magnetic-deflection ion analysis that the only primary product of electron-impact ionization in hydrogen is H2+, with H+ arising by secondary dissociation of H2+ even when the electron energy is just above 16 V and H3+ appearing in large numbers as a secondary product at any voltage above 16 V — a 1925 interpretation under his apparatus conditions (see the note).
Source quote & editorial note
However, atomic ions H+ are readily produced by secondary dissociation of molecular ions H2+ even when the maximum electron energy is just above 16 volts.
Smyth, Primary and Secondary Products of Ionization in Hydrogen — Physical Review 25, 452 (1925) — p. 1
Editorial note, tabletop extrapolation: Smyth's conclusion is a 1925 reading of one apparatus. Modern data separate the channels: H2+ formation opens at the H2 ionization energy, 15.426 eV, and direct dissociative ionization H2 + e → H+ + H + 2e opens at an appearance energy of 18.078 eV (NIST Chemistry WebBook, hydrogen ion-energetics page, retrieved 2026-10-25), so H+ is not only a secondary product; H3+ is chiefly made by the ion–molecule reaction H2+ + H2 → H3+ + H. For a filament source in hydrogen the species mix is therefore set by electron energy (which primary channels are open) together with gas density, the field and the ions' path and residence time (how much secondary chemistry occurs).
-
In Smyth's apparatus (electron energy V1 = 25 V, ion-drift voltage V2 = 9 V) H2+ dominated at low pressure — at 0.0007 mm Hg the spectrum was essentially H2+ — while at 0.053 mm the H2+ peak was almost suppressed and H3+ and H+ were strong. Figure reading of Fig. 4: H2+ peaks near 0.002 mm and falls steeply; H3+ overtakes it near 0.012–0.015 mm and peaks near 0.03 mm; H+ (plotted ×10) keeps rising to 0.07 mm.
Source quote & editorial note
Curve B in Fig. 3 taken at p=.053 mm shows the H2+ peak almost suppressed but H3+ and H+ strong.
Smyth, Primary and Secondary Products of Ionization in Hydrogen — Physical Review 25, 452 (1925) — p. 6
Editorial note, tabletop extrapolation: The pressures are those of Smyth's ionization region with fixed fields, electron energy and path length: there, raising pressure suppressed H2+ and strengthened H3+ and H+. What transfers to another source is the collision-depth comparison (gas density times path length, or residence time), not the direction of the H+ change, which depends on the electron-energy distribution, ion collision energy and competing reactions and must be measured or modelled for the source in hand.
-
From fitting H2+ intensity against pressure, Smyth concluded that the H2+ ion is dissociated about 15 times more readily than a kinetic-theory collision rate would give — it travels on average only about 1/15 of the ordinary molecular free path before being destroyed.
λ_destroy = 1/(n2·p) with n2 ≈ 1100 per cm of path at 1 mm Hg (Smyth's fit in his apparatus, drift field V2/d2 = 9 V/0.8 cm). Units: n in collisions per cm at 1 mm Hg. Smyth takes 1/n0 ≈ 4√2 × molecular free path with n0 = 12.7 /cm, so his molecular rate ≈ 4√2 × 12.7 ≈ 72 /cm (free path ≈ 0.014 cm at 1 mm); 1100/72 ≈ 15. Example: p = 1e-3 mm Hg (0.13 Pa) → λ ≈ 1/1.1 cm ≈ 0.9 cm; p = 1e-4 mm Hg (0.013 Pa) → λ ≈ 9 cm.Source quote & editorial note
the distance an H2+ ion travels on an average before being dissociated by the proximity of another molecule is about 1/15 the ordinary molecular free path.
Smyth, Primary and Secondary Products of Ionization in Hydrogen — Physical Review 25, 452 (1925) — p. 10
Editorial note, tabletop extrapolation: The short destruction length means H2+ survives only where gas density times path length is small: on Smyth's fitted value, about 0.9 cm of path at 1e-3 mm Hg and about 9 cm at 1e-4 mm Hg (formula). Where a source's pressure and ion path put the product near or above one destruction length, the molecular-ion fraction will track pressure strongly; the fitted rate belongs to Smyth's fields and ion energies and must be re-measured for another geometry.
-
Smyth found H+ always relatively much stronger than H3+ at higher values of V2, the field through which ions drift after formation; he explains that faster ions associate with neutral molecules less readily while dissociation continues. Figure reading of Fig. 5 (V2 = 25 V, p ≈ 0.045 mm): H+ peak tallest, H3+ next, H2+ smallest. [Notation: Smyth's V1 is the electron energy, V2 the field across the ionization region and V3 the analyser accelerating voltage; Hogness and Lunn use (V1+V2) for electron energy, V3 for the ion draw-out (electron-retarding) field and V4 for acceleration, so Smyth's V2 corresponds to their V3 (they call it V3 when citing Smyth's Fig. 6, p.9 of their paper). The subscript on p.12 is hard to read in the scan; the "only condition varied" sentence fixes it as V2.]
Source quote & editorial note
It was found that H+ was always relatively much stronger at higher values of V2
Smyth, Primary and Secondary Products of Ionization in Hydrogen — Physical Review 25, 452 (1925) — p. 12
Editorial note, tabletop extrapolation: The transfer is a testable hypothesis, not a prediction: in a cyclotron source the local extraction or bias field at the source exit — distinct from the dee voltage itself, which is RF, phase-dependent and may be screened inside the source — changes ion residence time and so may alter the species mix. Faster removal can also preserve H2+ rather than shift H3+ toward H+. This is a 1925 tube result; the direction and size of any shift in a particular source must be measured.
-
Smyth reports that the apparent appearance potential of H+ or H3+ can be observed anywhere between 16 and 21 V depending on pressure and fields, and concludes that H3+ and H+ may appear as secondary products as soon as H2+ is formed; on his reading, measured 'critical potentials' for the secondary species reflect the tube conditions rather than separate primary processes.
Source quote & editorial note
an apparent critical potential for H+ or H3+ may be observed anywhere between 16 and 21 volts if conditions are properly chosen and not varied.
Smyth, Primary and Secondary Products of Ionization in Hydrogen — Physical Review 25, 452 (1925) — p. 15
Editorial note, tabletop extrapolation: A species threshold measured on one source by sweeping electron energy says as much about that source's pressure and geometry as about hydrogen; the 16–21 V range is Smyth's observed apparent-onset range. Modern thresholds separate the channels: H2+ opens at 15.426 eV and direct dissociative production of H+ at 18.078 eV (NIST Chemistry WebBook, hydrogen ion-energetics page, retrieved 2026-10-25), so electron energy does change which primary channels are open, while pressure and residence time set the secondary chemistry. Just above threshold the H2 ionization cross section rises with electron energy; it peaks and then falls, so more energy is not simply more yield.
-
Hogness and Lunn (California, 1925) measured the H+/H2+/H3+ mix from below 1e-5 mm to 0.006 mm Hg and electron energies to 60 V: at low pressure only H2+ forms, H3+ rises in proportion to pressure up to 60 %, and in their apparatus H+ did not exceed 4 % (abstract). [Source-internal: Table I lists H+ up to 4.6 % (at 39 × 1e-4 mm) and H3+ up to 57.8 %, against the abstract's "did not exceed 4 percent" and "reached 60".]
Source quote & editorial note
At low pressures only H2+ ions are formed; as the pressure is increased the percentage of H3+ increases in proportion to the pressure. In the apparatus used, the percentage of H+ increased with pressure but did not exceed 4 percent
Editorial note, tabletop extrapolation: This is the second independent 1925 measurement agreeing with Smyth that H2+ dominated the low-pressure spectrum and that H3+ was secondary. The low H+ fraction is specific to their tube, drawing-out field and sensitivity: H+ also has a direct dissociative-ionization channel from H2 (appearance energy 18.078 eV, NIST Chemistry WebBook, hydrogen ion-energetics page, retrieved 2026-10-25), open at their electron energies, and both papers show more H+ with stronger fields along the ion path.
-
Hogness and Lunn's Table I (pressure in units of 1e-4 mm Hg, electron energy 48 V): 2.0 → ~97 % H2+, 2 % H3+, ~1 % H+; 11-12 → ~80-87 % H2+, 12-18 % H3+; 28 → ~67 % H2+, ~30 % H3+; 56 (0.0056 mm) → 39 % H2+, ~58 % H3+, ~3 % H+. Figure reading of Fig. 2: H2+ and H3+ cross near 0.0045 mm (interpolating Table I gives ~0.0047 mm).
Units: 56 × 1e-4 mm Hg = 5.6e-3 Torr = 0.75 Pa; 2.0 × 1e-4 mm Hg = 2e-4 Torr = 0.027 Pa (1 Torr = 133.3 Pa)Source quote & editorial note
Electrometer readings for H+, H2+ and H3+ peaks, at various pressures.
Editorial note, tabletop extrapolation: The table gives the measured shape of the species-versus-pressure curve for one positive-ray tube. The crossover sits at a lower pressure than the figure reading of Smyth's Fig. 4 recorded under dg-2085 (H3+ overtaking H2+ near 0.012–0.015 mm); different ion path lengths are one plausible apparatus-dependent cause, not an established explanation, and a source's own crossover depends on its geometry and other parameters and must be measured.
-
When Hogness and Lunn raised the field drawing ions out (their V3, which corresponds to Smyth's V2) of the ionization chamber, so that ionization occurred only in its upper part and ions had a longer path, (H+ + H3+) rose relative to H2+ and the H+/H3+ ratio rose considerably; their 1925 interpretation is that H2+ ions gaining energy before collision break up more often, the H+ formed escaping instead of being captured into H3+.
Source quote & editorial note
When the field V3 for drawing the positive ions from the ionization chamber was increased to such an extent that only those electrons in the upper portion of that chamber had ionizing energy, then the intensity of (H+ + H3+) increased with respect to that of H2+ as might have been expected since the H2+ ions formed had greater chance of collision having had a longer average path to traverse. Moreover the ratio of the intensity of H+ to that of H3+ increased considerably. This can be explained by assuming that under these conditions H2+ breaks up more often in accordance with reaction (4). This would happen if when the H2+ ion acquires sufficient energy before colliding, the H+ formed by disruption were not trapped by the H2 molecule to form H3+ but passed on into the resolving chamber.
Editorial note, tabletop extrapolation: Ion path length through gas and the energy ions gain along it both move the species mix; in a cyclotron source these are set by where the filament's electrons ionize relative to the exit and by the field at the exit. The mechanism is the authors' historical reading (the fast H3+ channel in modern terms is H2+ + H2 → H3+ + H, not capture of H+); the effect is described qualitatively, and its size in a given source must be measured.
-
Hogness and Lunn's summary recipe for which hydrogen ion predominates in a discharge tube: H+ with a large potential difference plus some inert gas such as helium; H2+ at low pressure; H3+ at higher pressure and small potential difference. With helium and a little hydrogen the last line of their Table I reads 68.6 % H+, 31.4 % H2+ and no H3+ (table reading).
Source quote & editorial note
It was found that the percentage of H+ was much greater than in pure hydrogen at any pressure employed (see the last line of Table 1). […] For H+, large potential difference and the addition of some inert gas such as helium; for H2+, low pressures; for H3+, higher pressure and small potential difference.
Editorial note, tabletop extrapolation: This gives the direction of each knob for a hydrogen filament source: pressure, and the applied potential in their tube (the draw-out field V3), which is not interchangeable with a generic source field without identifying which electrodes and ion energies are affected. The helium-admixture result shows that collision partners matter as well as density; using it on a cyclotron would change gas load and discharge behaviour, and the paper does not test that.
-
Hogness and Lunn found that lowering electron energy from 48 V toward 18 V at fixed pressure raised the H3+ percentage, and suggest ("may be explained") that in their tube the stronger retarding field moved ionization farther from the exit gauze, lengthening the ion path, rather than the impact energy being responsible.
Source quote & editorial note
The increase in the percentage of H3+ as (V1+V2) was decreased may be explained as due to the characteristics of the discharge tube. The now greater retarding potential V3 decreased the velocity of those electrons that had penetrated into the lower part of the ionization chamber to such an extent that they could not ionize the gas. The H2+ ions were then formed at a greater average distance from the gauze J, and having had a longer average path to traverse had greater chance to collide.
Editorial note, tabletop extrapolation: In a filament source, changing electron energy can also change where ionization happens, so an apparent energy dependence of the species mix may be a geometry effect; whether an arc or bias scan moves the place of ion birth depends on the electrode geometry and must be checked rather than assumed.
-
Halbach's scaling argument for permanent magnets: scaling an electromagnet down at fixed field requires coil current density rising as 1/L, which cooling (or superconductor limits) eventually forbids, whereas permanent-magnet systems scale to any size with no loss of field; so, under similar geometry and the stated steady-state current-density limits, there is a critical size below which a permanent magnet gives the higher field.
Referee sketch of Halbach's verbal argument (not an equation in the paper): electromagnet at fixed B scaled by factor L: NI ~ B x gap ~ L, coil area ~ L^2 -> j = NI/area ~ 1/L. PM: B unchanged under uniform scaling (fields depend on geometry ratios and Br only).Source quote & editorial note
current density in the coils scales inversely proportional to the linear dimensions L of the magnet. Since superconductors have an upper limit for the current density j that can be carried, and dissipative coils have an upper limit for j due to the need to remove the dissipated power, j needs to be reduced below that prescribed by simple scaling when L reaches a certain small value that depends, of course, on many details of the magnet design. […] PM, on the other hand, can be scaled to any size without any loss in field strength. From this follows that when it is necessary that a magnetically significant dimension of a magnet is very small, a permanent magnet will always produce higher fields than an electromagnet.
Editorial note, tabletop extrapolation: For a small cyclotron the magnetically significant dimensions are the gap (centimetres) and the pole radius; the argument says where permanent magnets gain most. The critical size depends, per Halbach, on field, configuration and available material, so it must be evaluated for the specific design rather than assumed.
Cited in: Cyclotron Magnet Design
-
In a magnet with no free currents Ampere's law makes the line integral of H along a full infinite line (closed at infinity) zero: for a hybrid dipole with a hole drilled along its vertical centre line, used as a solenoid-like field source, the line integral of H along the whole centre line is identically zero, so the field on axis must reverse somewhere outside the gap.
Integral from -inf to +inf of H_z dz along the axis = 0 (no conduction current enclosed). Halbach 1980 (LBL-11393, Eq. 8) gives the same condition for solenoid-symmetry REC devices.Source quote & editorial note
for this kind of solenoidal field magnet, the line integral of H along the whole vertical centerline is, because of Ampere's Law, always identically zero.
Editorial note, tabletop extrapolation: A permanent-magnet cyclotron with an axial-injection bore through pole, magnet layer and yoke is this case: the gap field along the axis must be balanced by reversed axial field elsewhere on the axis, so the injection line passes through field reversals that a coil-driven magnet (whose axial H integral equals its ampere-turns) is not forced to impose. Injection optics through the bore should be computed with the full on-axis field.
Cited in: Cyclotron Magnet Design
-
Pure (iron-free) rare-earth magnets have a hard limit: for given geometric parameters the obtainable field is limited by the properties of the magnet material, so Halbach advises that high-field undulators with large gaps and period lengths be built with conventional technology.
Source spelling "charaterizes" (sic). Context: undulator Eqn. (5), stated as general for pure REC magnets.Source quote & editorial note
Eqn. (5) shows the same hard limit on the achievable field strength that charaterizes all pure REC magnets: for given geometric parameters, the obtainable field strength is limited by the properties of the REC material.
Editorial note, tabletop extrapolation: Halbach's hard limit is stated for iron-free REC magnets. A PM-energized cyclotron is a hybrid, in which high-permeability poles may permit flux concentration beyond the corresponding open-magnet field; whether they do, and what limits first — pole or return-yoke saturation, the magnet recoil line and demagnetization margin, leakage — has to come from a magnetic-circuit or FEM analysis early in the design.
Cited in: Cyclotron Magnet Design