☣️ Cyclotron Radioisotope Production Yield Simulator

The Cyclotron — Spiraling Charged Particles to Megaelectronvolt Energies

A cyclotron accelerates charged particles — typically protons (H⁺ or H⁻) or deuterons (²H⁺/⁻) for medical isotope production — using the elegant combination of a static magnetic field for confinement and an alternating radiofrequency electric field for repeated acceleration, spiraling the beam outward through hundreds of orbits until it reaches the energy needed to trigger a nuclear reaction in a target.

• 30–100 MHz: RF frequency range (alternating dee voltage)
• 1–2 T: Magnetic field strength (static, vertical confinement)
• 10–18 MeV: Typical medical proton energy (compact self-shielded units)
• ~100s: Orbits to full energy (spiral path, radius grows each pass)

The isochronous cyclotron principle

A charged particle moving in a uniform magnetic field B follows a circular path due to the Lorentz force, with cyclotron (orbital) frequency:

f = qB / (2πm)

Remarkably, this frequency is independent of the particle's orbital radius or speed (in the non-relativistic limit) — a proton at low energy near the center and the same proton after many acceleration cycles near the outer edge of the machine both complete an orbit in the same time, provided the field is uniform. This "isochronous" property is what allows a cyclotron to accelerate particles using a single, fixed RF frequency rather than needing to continuously ramp the frequency as in a synchrotron.

Two hollow D-shaped electrodes ("dees") are separated by a narrow gap and connected to an RF voltage source oscillating at the cyclotron frequency. Each time the particle crosses the gap between the dees, the electric field has reversed polarity from the previous crossing, so the particle is accelerated — never decelerated — gaining a fixed increment of kinetic energy on every half-orbit. As kinetic energy increases, orbital radius grows (r = mv/qB), producing the characteristic outward spiral trajectory.

At relativistic energies the isochronous condition breaks down slightly because relativistic mass increase reduces the cyclotron frequency; modern cyclotrons compensate with a magnetic field that increases slightly with radius (azimuthally varying field, AVF cyclotrons), maintaining isochronism up to several hundred MeV.

Negative-ion acceleration and stripping extraction

Most modern medical cyclotrons accelerate negative hydrogen ions (H⁻: a proton with two bound electrons) rather than bare protons (H⁺). This choice enables an elegant and highly efficient extraction method:

• H⁻ ions are accelerated through the same spiral process as positive ions, but curve in the opposite direction for a given magnetic field orientation • At the desired final energy and radius, the H⁻ beam is passed through an extremely thin carbon foil (a few micrograms/cm² thick) • The foil strips both electrons from each H⁻ ion in a fraction of a nanosecond, instantly converting it to a bare proton (H⁺) • Because the newly stripped proton has the same momentum but opposite charge sign, the magnetic field bends it in the opposite direction — curving it out of the machine and directly toward an external target or beamline with no additional electrostatic deflection hardware needed

This stripping-foil extraction is highly efficient (>95%) and allows simultaneous extraction of two independent beams by using two stripping foils at different azimuthal positions — enabling parallel production of two different isotopes from a single cyclotron.

Beam focusing and vacuum requirements

Maintaining a tightly collimated, stable beam over hundreds of orbits requires both vertical focusing (to counteract the tendency of particles to drift out of the median plane) and a high-vacuum environment (typically ~10⁻⁶ to 10⁻⁷ torr) to prevent beam loss from collisions with residual gas molecules.

Vertical focusing in AVF cyclotrons is achieved through the same azimuthally varying magnetic field geometry used for relativistic isochronism compensation — alternating higher- and lower-field sectors ("hills" and "valleys") produce a net restoring force that keeps the beam confined to the median plane throughout acceleration, a design pioneered in the 1950s–60s that enabled cyclotrons to reach much higher energies than the original uniform-field machines.

Target Bombardment — Turning Stable Isotopes into PET Radiotracers

The accelerated beam is directed onto an isotopically enriched target, where nuclear reactions transmute stable target nuclei into radioactive product nuclei. Each reaction has a characteristic energy threshold and an excitation function (cross-section vs. energy) that determines the optimal bombarding energy for maximum yield.

• 2.57 MeV: ¹⁸O(p,n)¹⁸F threshold (typical run at 11–18 MeV)
• 109.8 min: ¹⁸F half-life (most widely used PET isotope)
• 20.4 min: ¹¹C half-life (from ¹⁴N(p,α)¹¹C)
• 2.04 min: ¹⁵O half-life (from ¹⁴N(d,n)¹⁵O)

Fluorine-18 production — the workhorse reaction

The dominant clinical PET isotope, ¹⁸F, is produced almost universally via the ¹⁸O(p,n)¹⁸F reaction: a proton strikes an oxygen-18 nucleus, is captured, and a neutron is emitted, transmuting ¹⁸O (stable, ~0.2% natural abundance) into ¹⁸F (radioactive, t½ = 109.8 min, decays by positron emission).

The target is isotopically enriched ¹⁸O-water (typically >95–98% enriched, since natural water is almost entirely ¹⁶O), loaded into a small-volume (1–3 mL) target chamber. Enrichment is essential and expensive — ¹⁸O-water target liquid is recovered and recycled after each production run wherever possible to control costs.

The reaction threshold energy is 2.57 MeV, but the cross-section rises well above threshold and typical medical cyclotrons bombard at 11–18 MeV proton energy to maximize the integrated yield across the beam's energy-loss path through the target (protons continuously lose energy as they traverse the target volume, so the reaction proceeds across a range of energies from the incident value down toward threshold).

Carbon-11, oxygen-15, and other short-lived isotopes

Beyond fluorine-18, cyclotrons produce several other PET isotopes with much shorter half-lives, demanding on-site or near-site production immediately before use:

• Carbon-11 (t½ = 20.4 min): produced via ¹⁴N(p,α)¹¹C — protons bombard nitrogen gas (often with a small oxygen or hydrogen additive to control the chemical form of the product, e.g. ¹¹CO₂ or ¹¹CH₄), used to label tracers like ¹¹C-methionine or ¹¹C-Pittsburgh compound B (amyloid imaging) • Oxygen-15 (t½ = 2.04 min): produced via ¹⁴N(d,n)¹⁵O using a deuteron beam on nitrogen gas, used for ¹⁵O-water blood flow studies — but its extremely short half-life requires the cyclotron to be essentially adjacent to the PET scanner • Nitrogen-13 (t½ = 9.97 min): produced via ¹⁶O(p,α)¹³N, used for ¹³N-ammonia myocardial perfusion imaging

These short-lived isotopes cannot be transported any meaningful distance and fundamentally constrain hospital-based cyclotron siting — the production facility must be within the same building complex as the imaging suite.

Excitation functions and optimal bombarding energy

Every nuclear reaction has a cross-section σ(E) — the probability of reaction per unit target atom density per unit beam fluence — that varies strongly and non-monotonically with the incident particle energy, typically rising from zero at threshold, peaking at some characteristic energy (often tens of MeV above threshold for (p,n) and (p,α) reactions), and declining or plateauing thereafter as competing reaction channels open up.

Cyclotron operators select a bombarding energy that balances several factors: maximizing the energy-integrated cross-section (the excitation function convolved with the proton's energy-loss profile through the target thickness), avoiding energy regions where competing reactions produce unwanted radioactive or stable contaminant isotopes, and respecting the practical energy range available from the specific cyclotron model installed.

Saturation Yield Curves — Why Longer Irradiation Has Diminishing Returns

As a target is bombarded, radioactive product atoms accumulate — but they are simultaneously decaying. The activity produced follows a characteristic saturation curve that asymptotically approaches a maximum value, meaning that beyond roughly two to three half-lives of the product isotope, additional irradiation time yields comparatively little extra activity.

• 50%: Yield at 1 half-life (of saturation activity)
• 75%: Yield at 2 half-lives (of saturation activity)
• 87.5%: Yield at 3 half-lives (of saturation activity)
• 93.75%: Yield at 4 half-lives (of saturation activity)

Deriving the saturation yield equation

During bombardment, the rate of change of the number of radioactive product nuclei N(t) balances production against radioactive decay:

dN/dt = R − λN(t)

where R is the constant production rate (proportional to beam current, target thickness, and reaction cross-section) and λ = ln(2)/t½ is the decay constant of the product isotope. Solving this first-order differential equation with N(0)=0 gives:

A(t) = A_sat · (1 − e^(−λt))

where A(t) = λN(t) is the activity at time t, and A_sat = R is the saturation activity — the maximum activity achievable at infinite irradiation time, when production and decay reach equilibrium.

This exponential-approach-to-saturation behavior means the fraction of saturation activity reached depends only on irradiation time expressed as a multiple of the product's half-life, not on the absolute half-life value itself:

• t = 1×t½: 1−2⁻¹ = 50% • t = 2×t½: 1−2⁻² = 75% • t = 3×t½: 1−2⁻³ = 87.5% • t = 4×t½: 1−2⁻⁴ = 93.75% • t = 5×t½: 1−2⁻⁵ = 96.9%

Because ¹⁸F has a half-life of 109.8 minutes, routine clinical production runs are typically limited to roughly 60–120 minutes (about 0.5–1× half-life) — a pragmatic balance between achievable batch activity and total production-plus-processing time, since chemical synthesis, quality control, and delivery all consume additional time against the isotope's decay clock.

The full yield formula — linking physics to activity

The saturation activity itself depends on the beam and target parameters through:

A_sat = σ · (I/q) · N_target · (1 − e^(−μx))

where σ is the reaction cross-section (cm²), I/q is the beam particle flux (particles/second, from beam current I divided by particle charge q), N_target is the areal density of target nuclei (atoms/cm²) as seen by the beam, and the exponential term accounts for the fraction of the beam actually stopped in (and reacting with) the target of thickness x.

Combining this with the time-dependent saturation factor gives the full expression for activity produced after irradiation time t:

A(t) = σ · (I/q) · N_target · (1 − e^(−λt))

This single formula captures why cyclotron operators care about beam current (linear yield scaling), target design (thickness and enrichment determine N_target), and irradiation time relative to the product half-life (the saturation factor) as the three primary levers for maximizing batch activity.

Practical implications for production scheduling

The diminishing-returns shape of the saturation curve has direct operational consequences for radiopharmacy scheduling:

• For very short-lived isotopes like ¹⁵O (t½=2 min) or ¹³N (t½=10 min), irradiation runs are necessarily brief (a few minutes) since going much beyond 2–3 half-lives provides negligible additional yield while consuming scarce cyclotron beam time • For longer-lived isotopes like ¹⁸F, production facilities often run irradiations of roughly 0.5–2× half-life, balancing achievable end-of-bombardment activity against total cyclotron occupancy time (since the same cyclotron may need to produce multiple batches or multiple isotopes per day) • Doubling irradiation time from 1 to 2 half-lives only increases yield from 50% to 75% of saturation (a 1.5× gain, not 2×) — while doubling beam current increases yield linearly (2× gain) — making beam current a generally more efficient lever than extended irradiation time once past roughly 1–2 half-lives

Beam Current, Target Engineering, and Specific Activity

Yield scales linearly with beam current in the ideal case, but real targets face physical limits — heat deposition, pressure buildup from radiolysis, and window/foil degradation — that cap the practical beam current a given target design can sustain. Target and beam-line engineering is as important to final isotope yield as the underlying nuclear physics.

• 10–150 μA: Typical medical beam current (proton beam on target)
• several hundred mCi – multi-Ci: ¹⁸F batch yield (typical run) (depends on current, time, target)
• 1–3 mL: Target volume (¹⁸O-water) (liquid target, silver/niobium body)
• Havar foil: Target window material (cobalt-chromium alloy, high strength)

Beam current — linear yield scaling with a practical ceiling

From the saturation yield formula, activity produced scales directly and linearly with beam current I — doubling beam current doubles the production rate and, correspondingly, the saturation activity achievable for a given irradiation time. This makes higher beam current one of the most direct ways to increase batch yield without extending run time.

However, higher beam current deposits proportionally more power into the target (beam power = current × energy, so a 60 μA, 16 MeV beam deposits roughly 1 kW of power into a target volume of just a few milliliters). This intense, localized heating must be managed by:

• Efficient target cooling (helium gas cooling of the front window, water/liquid cooling of the target body) • Target geometry designed to maintain liquid-phase target material despite substantial heating (elevated operating pressure to raise the boiling point of the ¹⁸O-water target liquid) • Managing radiolytic gas generation — intense radiation exposure decomposes water into hydrogen and oxygen gas, which must be accommodated by target headspace and pressure-relief design without venting valuable ¹⁸O-enriched target material

Liquid target engineering for fluorine-18

The dominant target design for ¹⁸F production is a small-volume (typically 1–3 mL) liquid ¹⁸O-water target housed in a body machined from silver or niobium — chosen for their combination of high thermal conductivity, corrosion resistance to the acidic radiolysis products formed in irradiated water, and minimal activation by the proton beam (avoiding production of unwanted long-lived radioactive contaminants in the target body itself).

The target is sealed at the beam-entrance side by a thin metal foil, most commonly Havar (a cobalt-chromium-nickel-iron superalloy) chosen for its high tensile strength, allowing it to be made thin enough (tens of micrometers) to minimize beam energy loss and heating while still safely containing the pressurized target contents (routinely tens of atmospheres during irradiation).

After irradiation, the ¹⁸F-containing water is transferred (typically by helium gas pressure or vacuum) to an automated radiochemistry synthesis module for conversion into the final radiopharmaceutical (e.g., ¹⁸F-FDG via nucleophilic substitution chemistry), while the recovered, now-depleted ¹⁸O-water is collected for re-enrichment and reuse given its high cost.

Specific activity — why it matters beyond total yield

Specific activity (or molar activity) describes the radioactivity per unit mass (or per mole) of the compound of interest — a measure of how "diluted" the radioactive atoms are by chemically identical but non-radioactive ("cold") atoms of the same element.

Even a target enriched to >98% ¹⁸O still contains trace natural fluorine contamination (from target body corrosion, water purity, or target system materials) that competes with radioactive ¹⁸F during the subsequent radiochemical synthesis, diluting specific activity. For imaging applications like FDG-PET, where the injected tracer mass is negligible relative to physiological glucose pools, modest specific activity is generally adequate.

However, for receptor-binding radiotracers (e.g., neuroreceptor ligands), high specific activity is critical: if too much non-radioactive "cold" carrier compound is present alongside the radiolabeled tracer, it competitively occupies and saturates the limited number of receptor binding sites in vivo, distorting the measured binding signal and potentially causing unwanted pharmacological effects from the mass of unlabeled compound itself. Achieving no-carrier-added (NCA) or carrier-free production — minimizing all sources of stable isotopic or chemical contamination — is therefore a major target and radiochemistry engineering priority for such tracers.

Medical Cyclotron Facilities — Energy Ranges, Shielding, and Siting

Medical PET-isotope production is dominated by a class of compact, self-shielded cyclotrons operating in the 10–18 MeV range, purpose-built for hospital or radiopharmacy siting. Larger, higher-energy cyclotrons serve regional production centers for SPECT isotopes and bulk radioisotope manufacturing, each with distinct shielding, cost, and logistics requirements.

• 10–18 MeV: Medical PET cyclotron energy (compact, self-shielded units)
• 30+ MeV: Regional isotope centers (higher yield, SPECT isotopes)
• $2–5M: Compact system cost (installed, self-shielded)
• ~2 hours: ¹⁸F transport radius (limited by 109.8 min half-life)

Compact self-shielded medical cyclotrons

Commercial medical cyclotrons designed specifically for PET isotope production — such as the GE PETtrace, Siemens Eclipse, and IBA Cyclone series — typically accelerate protons to 10–18 MeV, sufficient to efficiently drive the (p,n) and (p,α) reactions used for ¹⁸F, ¹¹C, ¹³N, and ¹⁵O production without the added cost and shielding burden of higher-energy machines.

These systems are engineered as "self-shielded" units: the cyclotron vault incorporates its own integrated radiation shielding (thick concrete, sometimes combined with steel and borated polyethylene layers to attenuate both the intense fast-neutron flux generated by (p,n) reactions and secondary gamma radiation), allowing installation within a standard hospital basement or dedicated radiopharmacy building without requiring a separate, massive bunker structure — dramatically reducing facility construction cost and footprint compared to older, unshielded research cyclotrons.

Higher-energy regional production centers

Larger cyclotrons operating at 30 MeV and above serve a different role: bulk, high-volume production of a broader range of isotopes for regional or national distribution, including many SPECT isotopes (e.g., ²⁰¹Tl via ²⁰²Hg(p,2n)²⁰¹Tl, ⁶⁷Ga via ⁶⁸Zn(p,2n)⁶⁷Ga) that require higher bombarding energies to access the relevant (p,2n) or (p,3n) reaction channels efficiently.

These facilities benefit from higher achievable beam currents and energies to maximize batch yields for isotopes with half-lives measured in days rather than minutes (e.g., ⁶⁷Ga t½=3.3 days, ²⁰¹Tl t½=73 hours), where the longer half-life permits shipment over much greater distances — making centralized, high-throughput production economically favorable compared to distributed hospital-based cyclotrons.

Siting logistics and cumulative cost

The short half-lives of the most commonly used PET isotopes fundamentally dictate facility siting. ¹⁸F, with a 109.8-minute half-life, tolerates ground or short-flight transport within roughly a 2-hour radius before decay erodes too much of the delivered activity to be clinically useful — supporting a regional radiopharmacy model serving multiple hospitals within driving distance of a single cyclotron. Carbon-11 (20.4 min) and oxygen-15 (2.04 min) tolerate essentially no transport at all, requiring the cyclotron to be physically adjacent to (often directly beneath) the PET imaging suite that will use them.

A fully installed, self-shielded compact medical cyclotron facility — including the cyclotron itself, target systems, automated radiochemistry synthesis modules, quality control laboratory, and hot-cell infrastructure for safe handling — typically costs in the range of $2–5 million, a substantial capital investment that is generally justified only by sufficient local clinical PET imaging volume (or by serving as a regional isotope hub for several nearby imaging centers) to keep the cyclotron in near-continuous production use.

Because of self-shielding and compact bunker design, many modern medical cyclotron installations achieve full regulatory compliance for siting within occupied hospital buildings — a capability that would have been unthinkable with the room-sized, minimally shielded research cyclotrons of the mid-20th century.