---
title: Fusion Neutron-Driven Transmutation
url: https://www.emergentmind.com/topics/fusion-neutron-driven-transmutation
type: topic
---

# Fusion Neutron-Driven Transmutation

Fusion neutron-driven transmutation is the modification of nuclide inventories by fusion-produced neutrons, most commonly the \(14.1\ \mathrm{MeV}\) neutrons from deuterium–tritium fusion, \(d+t\rightarrow \alpha(3.5\ \mathrm{MeV})+n(14.1\ \mathrm{MeV})\). In fusion systems these neutrons simultaneously drive unwanted composition changes in structural materials, such as tungsten plasma-facing armor, and intentional conversion of feedstocks into transuranic burn products, stable metals, medical radioisotopes, and \(\beta^{-}\) emitters. Because the birth spectrum is sufficiently hard to open threshold reactions such as \((n,2n)\), \((n,p)\), and \((n,\alpha)\), fusion transmutation is treated through coupled neutron transport, depletion, blanket design, and radiation-effects materials analysis [1612.03892, 2511.02814].

## 1. Reaction physics and transmutation kinetics

The basic physical principle is that a target nuclide in a neutron field undergoes reactions at a rate determined by the local spectral flux and the energy-dependent microscopic cross section. In the continuous-energy form used across transport and depletion studies, a reaction channel is written as
\[
R_i(r)=\int_0^\infty \phi(r,E)\,\sigma_i(E)\,dE,
\]
and the isotopic inventory evolves through Bateman-type balance equations with both decay and neutron-induced source and sink terms [2109.08741, 1612.03892]. In a monoenergetic or one-group approximation the same idea is often written as \(R=\Phi\,\sigma\,N\), or \(dC_i/d\Phi=\sigma_iN_W\) for tungsten transmutation products [1706.03252, 2512.09242].

The hard D–T spectrum opens reaction pathways that are weak or inaccessible in thermal or typical fission spectra. Representative examples include \(^{198}\mathrm{Hg}(n,2n){}^{197}\mathrm{Hg}\rightarrow{}^{197}\mathrm{Au}\), with a threshold \(E_{\mathrm{th}}\approx 9.1\ \mathrm{MeV}\); \(^{131}\mathrm{Xe}(n,p){}^{131}\mathrm{I}\); \(^{102}\mathrm{Ru}(n,\alpha){}^{99}\mathrm{Mo}\); and \(^{148}\mathrm{Nd}(n,2n){}^{147}\mathrm{Nd}\rightarrow{}^{147}\mathrm{Pm}\) [2507.13461, 2511.02814, 2605.20260]. In tungsten, neutron capture and \((n,2n)\) reactions produce tantalum and rhenium, while further capture on Re isotopes builds up osmium [1706.03252].

For isotope-production studies, the produced activity under constant flux is commonly written as
\[
A(t)=N_0\,\sigma(E_n)\,\Phi\,\bigl(1-e^{-\lambda t}\bigr),
\]
with saturation activity \(A_s=N_0\sigma\Phi\) [2511.02814]. For multistep chains, such as mercury-to-gold transmutation, population dynamics are described by chained depletion equations with neutron-driven transfers between adjacent isotopes and radioactive decay of the terminal precursor [2604.02590].

## 2. Transport, depletion, and spatially resolved modeling

Fusion transmutation calculations are typically performed by coupling neutron transport to inventory evolution. A detailed tungsten study computes \(\phi(r,E)\) with continuous-energy Monte Carlo using MCNP6 and then evolves local nuclide densities with FISPACT-II in each spatial cell [1612.03892]. In that implementation, FISPACT-II represents \(\phi(E)\) and \(\sigma(E)\) on a 660-group structure up to \(30\ \mathrm{MeV}\) (UKAEA-709) and applies probability-table self-shielding corrections through
\[
RR_i=\sum_{g=1}^G SSF_i(g)\,\sigma_i(g)\,\phi(g),
\]
where \(SSF_i(g)\le 1\) corrects resonance dilution within a group [1612.03892]. The same study resolves a \(2\ \mathrm{cm}\) tungsten armor shell into \(0.1\ \mathrm{mm}\) radial layers, with ADVANTG biasing used to achieve \(\sim 1\%\) statistical precision in most energy bins [1612.03892].

Blanket-scale transmutation studies use analogous workflows with different transport and depletion engines. An ARC-class mercury blanket study uses Paramak for geometry generation, DAGMC export, OpenMC v0.14 transport on a \(90^\circ\) toroidal sector with reflective boundary conditions, JENDL-5 nuclear data, and OpenMC’s CoupledOperator over \(365\) days in \(24\) timesteps while “harvesting” gold at each step [2507.13461]. Molten-salt transmutation studies solve fixed-source, \(k\)-eigenvalue, and subcritical problems with MCNP 6.2, SCALE/ORIGEN-S, and OpenMC, characterizing source amplification by
\[
M=\frac{1}{1-k_{\rm eff}},\qquad k_{\rm eff}<1
\]
[2109.08741].

A recurring result is that spatial resolution is not optional. In DEMO neutronics, the total neutron flux varies from \(\sim 8.3\times10^{14}\ \mathrm{n\,cm^{-2}\,s^{-1}}\) at the first-wall armor to \(\sim 3.9\times10^{13}\ \mathrm{n\,cm^{-2}\,s^{-1}}\) at the blanket backplate, and both dpa and gas-production rates follow that depth dependence [1311.5079]. In tungsten and blanket feedstocks alike, local moderation, backscatter, self-shielding, and geometrical interfaces set the actual transmutation field [1612.03892].

## 3. Tungsten as the benchmark fusion transmutation problem

Tungsten is widely treated as the reference case because it is the leading plasma-facing material candidate and because its transmutation products strongly affect service behavior. In spatially heterogeneous models of fusion armor, low-\(Z\) materials behind or within tungsten—especially water and steel—moderate and back-scatter neutrons, strongly increasing the eV–keV flux near interfaces [1612.03892]. In one representative configuration, the MeV flux falls only modestly, from \(\sim 3.2\times10^{14}\ \mathrm{n\,cm^{-2}\,s^{-1}}\) at the plasma side to \(\sim 2.5\times10^{14}\ \mathrm{n\,cm^{-2}\,s^{-1}}\) at the moderator side, yet the local \(^{186}\mathrm{W}(n,\gamma)\) rate can rise by up to a factor of \(7\) from front to back because moderated neutrons repopulate the giant capture resonances around \(1\)–\(30\ \mathrm{eV}\), with a dominant resonance at \(E_r\approx 18.8\ \mathrm{eV}\) [1612.03892].

After \(2\ \mathrm{fpy}\) of D–T irradiation in that scenario, the local rhenium concentration varies from \(C_{\rm Re}\approx 1200\ \mathrm{appm}\) \((0.12\ \mathrm{at.}\%)\) at the plasma-facing surface to \(\approx 3000\ \mathrm{appm}\) \((0.30\ \mathrm{at.}\%)\) at mid-thickness and \(\approx 5800\ \mathrm{appm}\) \((0.58\ \mathrm{at.}\%)\) at the back, a full-profile variation of \(\Delta C_{\rm Re}\approx 4600\ \mathrm{appm}\) [1612.03892]. A \(4\ \mathrm{mm}\) coolant layer directly behind tungsten boosts back-of-W Re by \(\approx 30\%\) over the reference case, while embedded \(4\ \mathrm{mm}\) internal water channels create sharp local Re peaks at the water–W interfaces [1612.03892]. By comparison, neglecting self-shielding factors over-predicts Re by only \(25\)–\(50\%\), which is smaller than the spatially induced variation [1612.03892].

Temperature broadening is present but weaker. Recomputed Re-vs-depth profiles at \(600\ \mathrm{K}\) and \(900\ \mathrm{K}\) show only a modest upward shift of \(\sim 2\)–\(5\%\) relative to room temperature, implying that Doppler broadening is subdominant to spatial heterogeneity in these geometries [1612.03892]. Experimentally, this matters because homogenized or blanket-averaged spectra can be badly misleading: comparison with HFR irradiation of tungsten at Petten shows that only a detailed geometry-plus-flux model reproduces the measured \(\sim 1.3\ \mathrm{at.}\%\) Re, whereas a blanket spectrum predicts \(\sim 5\ \mathrm{at.}\%\) [1612.03892].

The metallurgical consequences extend beyond composition. A first-principles matrix-Hamiltonian and thermodynamic-integration study of W–Re–Os–Vac and related systems shows that vacancies are decorated by Re and Os, but not by Ta, with the distinction traced to the sign of the Warren–Cowley short-range-order parameter \(\alpha_n^{AB}\) and to vacancy–solute binding energetics [2102.01614]. At \(T=800\ \mathrm{K}\) and \(1200\ \mathrm{K}\), Monte Carlo snapshots show Re and especially Os clustering around vacancy clusters, whereas Ta remains randomly mixed [2102.01614]. Under HFR-relevant conditions, the same framework predicts voids containing concentrated Re and Os, consistent with TEM and APT observations of Re “clouds,” Os “rods,” and core–shell clustering in irradiated tungsten [2102.01614].

## 4. Blanket-scale objectives: waste burning, isotope production, and elemental conversion

Fusion neutron-driven transmutation has been proposed for four distinct blanket-scale objectives: reducing minor actinides, producing radionuclides, manufacturing stable high-value products, and coupling these outputs to tritium-breeding blanket design [2109.08741, 2507.13461, 2511.02814, 2605.20260].

| Objective | Representative system or pathway | Reported result |
|---|---|---|
| Minor actinide transmutation | CAIL-driven molten salt transmutator, \(S_0\simeq 5\times10^{15}\ \mathrm{n/s}\), \(k_{\rm eff}\simeq 0.98\) | \(\simeq 7\ \mathrm{kg}\) MA/year; \(>90\%\) reduction of Am and Cm over \(2\)–\(3\) years |
| Gold synthesis | ARC-class tokamak with \(210\ \mathrm{mm}\) Hg/Li channel | \(^{197}\mathrm{Au}\) net yield \(\approx 2909\ \mathrm{kg/yr}\) at \(1500\ \mathrm{MW_{th}}\); scales to \(\approx 1.94\ \mathrm{t/GW_{th}/yr}\); TBR \(=1.19\) |
| Medical isotopes | Stable feedstocks irradiated by D–T neutrons | \(^{131}\mathrm{I}\): \(2.1/8.8\), \(^{177}\mathrm{Lu}\): \(2.4/10.0\), \(^{47}\mathrm{Sc}\): \(6.7/64.0\) g/(MW\(_{\rm th}\) yr) for natural/enriched feedstock |
| Betavoltaic isotopes | \(^{148}\mathrm{Nd}(n,2n)\rightarrow{}^{147}\mathrm{Pm}\) blanket | \(1424\ \mathrm{kg/yr}\) and \(1.3\times10^9\ \mathrm{Ci/yr}\) per GW\(_{\rm th}\) |

In molten-salt waste transmutation, the source is a distributed array of beam–target D–T neutron heads combined with a well-mixed FLiBe or FLiNaK core carrying Pu, Am, and Cm fluorides [2109.08741]. For \(S_0\simeq 5\times10^{15}\ \mathrm{n/s}\) and \(k_{\rm eff}\simeq 0.98\) \((M\simeq 50)\), the simulated thermal power is \(\simeq 1.42\ \mathrm{MW_{th}}\) and the transmutation rate is \(\simeq 7\ \mathrm{kg}\) MA/year in the static phase, while periodic in-line batch processing yields \(>90\%\) reduction of Am and Cm over \(2\)–\(3\) years [2109.08741].

For mercury-to-gold conversion, a \(1500\ \mathrm{MW_{th}}\) ARC-class tokamak with a circulating liquid alloy of \(85\ \mathrm{at}\%\) Hg \((90\%\ ^{198}\mathrm{Hg})\) and \(15\ \mathrm{at}\%\) Li \((90\%\ ^6\mathrm{Li})\) in a \(\sim 210\ \mathrm{mm}\) multiplier channel produces \(^{197}\mathrm{Au}\) at \(\approx 2909\ \mathrm{kg/yr}\), while preserving a simulated total TBR of \(1.19\), with \(0.694\) from the Hg/Li channel and \(0.490\) from FLiBe [2507.13461]. Because \((n,2n)\) on Hg contributes neutron multiplication, the study reports essentially zero penalty on net thermal power and continuous delivery of \(600\ \mathrm{MW_e}\) for a \(40\%\)-efficient \(1500\ \mathrm{MW_{th}}\) plant [2507.13461]. A related remediation study extends the concept to successive \((n,2n)\) steps on all stable mercury isotopes until \(^{197}\mathrm{Hg}\) or \(^{197m}\mathrm{Hg}\), emphasizing that at sufficiently high flux all stable Hg isotopes become eligible for conversion to gold [2604.02590].

For radionuclides, high-energy fusion neutrons enable both same-element and proton-changing routes. OpenMC slab studies at flux \(10^{14}\ \mathrm{n\,cm^{-2}\,s^{-1}}\) and \(20\ \mathrm{cm}\) thickness report production rates in g/(MW\(_{\rm th}\) yr) for multiple medical isotopes, including \(^{99}\mathrm{Mo}/^{99\mathrm{m}}\mathrm{Tc}\), \(^{131}\mathrm{I}\), \(^{177}\mathrm{Lu}\), \(^{161}\mathrm{Tb}\), \(^{47}\mathrm{Sc}\), and \(^{195m}\mathrm{Ir}\) [2511.02814]. A complementary economic study gives a specific low-gain example: a \(3\ \mathrm{MW}\) fusion system transmuting \(^{102}\mathrm{Ru}\rightarrow{}^{99}\mathrm{Mo}\) at \(Q_{\mathrm{plas}}=0.1\) and \(\Phi_0\approx 2\times10^{12}\ \mathrm{n/cm^2\,s}\) yields \(0.8\ \mathrm{kg/yr}\) and has \( \mathrm{NPV}>0\) for \(Q_{\mathrm{plas}}\gtrsim 0.02\) at \(\mathrm{CAPEX}\lesssim \$0.5\ \mathrm{B}\) [2512.09242].

For \(\beta^{-}\) emitters, blanket calculations identify \(^{147}\mathrm{Pm}\), \(^{63}\mathrm{Ni}\), \(^{39}\mathrm{Ar}\), and \(^{137}\mathrm{Cs}\) as prominent products [2605.20260]. In a tokamak blanket using enriched \(^{148}\mathrm{Nd}\) with \(^{6}\mathrm{Li}\), OpenMC simulations show \(^{147}\mathrm{Pm}\) production above \(1000\ \mathrm{kg/yr}\) per GW\(_{\rm th}\) with TBR \(\ge 1.1\), and the tabulated value for a fully enriched \(^{148}\mathrm{Nd}\) case is \(1424\ \mathrm{kg/yr}\), equivalent to \(\sim 1.3\times10^{9}\ \mathrm{Ci/yr}\) [2605.20260].

## 5. Source concepts, control architectures, and flux shaping

The transmutation problem depends strongly on the neutron-source concept. For materials qualification and fusion-prototypical irradiation, the IFMIF-Lite compact-cyclotron concept accelerates \(125\ \mathrm{mA}\) of deuterons to \(40\ \mathrm{MeV}\) onto a \(25\ \mathrm{mm}\) flowing lithium jet, producing a D–Li neutron spectrum peaked near \(14\ \mathrm{MeV}\) with a representative yield of \(\approx 1.3\times10^{16}\ \mathrm{n/s}\) and flux \(\approx 1.0\times10^{13}\ \mathrm{n/cm^2\,s}\) at \(r=10\ \mathrm{cm}\) [2302.09011]. In that framework, specimen positions \(5\)–\(15\ \mathrm{cm}\) downstream of the Li jet define a \(50\times 50\times 50\ \mathrm{mm^3}\) irradiation cavity, and multi-beam incidence is used to reduce the flux gradient across the test volume to below \(25\%\) [2302.09011].

Distributed neutron production has been proposed for subcritical molten-salt transmutators. In the CAIL concept, a high-repetition-rate array of coherently combined femtosecond fiber lasers irradiates nanometric deuterated foils, producing deuterons that strike a tritium-loaded catcher; the quoted scaling is \(Y_1\simeq 10^7\ \mathrm{n/shot}\) at \(f\simeq 10^5\ \mathrm{s^{-1}}\), or \(\simeq 10^{12}\ \mathrm{n/s}\) per laser head, scalable to \(\sim 10^3\) heads and \(S_0\simeq 10^{15}\)–\(10^{16}\ \mathrm{n/s}\) [2109.08741]. Source placement is optimized through a phyllotactic pattern, while real-time spectroscopy and an Evolutionary Algorithm coupled with a Neural Network tune individual source powers, processing schedules, and salt flow-rates between sub-volumes [2109.08741]. After optimization, Monte Carlo tallies show fast-spectrum flux uniformity within \(\pm 10\%\) [2109.08741].

Muon-catalyzed fusion and short-pulse laser sources extend the concept into high-flux or ultrashort-pulse regimes. In \(\mu\)CF, the absence of external plasma heating reduces wall loading to the \(\alpha\)-heating component alone, leading to an estimated allowable first-wall neutron flux of \(\Phi_0\simeq 2\times10^{15}\ \mathrm{n/cm^2\,s}\) for a \(10\ \mathrm{MW/m^2}\) heat-flux limit [2511.20951]. An example with \(10^8\) muons/s, \(N_{\mathrm{fus},\mu}=150\), and a \(10\ \mathrm{g}\) \(^{226}\mathrm{Ra}\) feedstock yields up to \(0.5\ \mathrm{mg}\) of \(^{225}\mathrm{Ac}\) per year [2511.20951]. Short-pulse DD fusion, although not a \(14\ \mathrm{MeV}\) D–T source, provides quasi-monoenergetic \(2.45\ \mathrm{MeV}\) neutrons with \(<20\ \mathrm{ps}\) pulse duration, \(\sim 100\ \mu\mathrm{m}\) source size, and local peak flux \(>10^{22}\ \mathrm{cm^{-2}\,s^{-1}}\) in \(0.1\ \mathrm{sr}\) bins, while LWFA-driven photonuclear sources trade lower peak brightness for higher repetition rates and a \(36\times\) advantage in cumulative capture events for short-lived isomers [2605.18968]. This suggests a division between blanket-scale continuous transmutation and pulsed fast-neutron transmutation experiments.

## 6. Materials consequences, qualification metrics, and unresolved constraints

Fusion transmutation cannot be separated from displacement damage and gas production. In DEMO first-wall conditions, integrated MCNP–FISPACT modeling gives for pure Fe at the outboard equatorial armor a helium production rate of \(\simeq 140\ \mathrm{appm\,fpy^{-1}}\), \(R_{\mathrm{He}}\approx 3.8\times10^{11}\ \mathrm{atoms\,cm^{-3}\,s^{-1}}\), and \(\simeq 150\ \mathrm{dpa\,fpy^{-1}}\), whereas tungsten in first-wall or divertor armor shows \(\simeq 4\)–\(9\ \mathrm{appm\,fpy^{-1}}\) He and \(\simeq 10\)–\(20\ \mathrm{dpa\,fpy^{-1}}\) [1311.5079]. Using DFT-based grain-boundary helium insertion energies, the same study estimates \(t_{\mathrm{emb}}\approx 5.5\ \mathrm{fpy}\) for Fe under first-wall conditions, but \(>700\ \mathrm{fpy}\) for W and Ta [1311.5079].

At the atomistic level, transmutation products feed back on tungsten behavior. First-principles calculations for W–Re, W–Os, and W–Ta alloys report that all transmutation elements lower the Pugh ratio \(G/B\), with pure W at \(G/B=0.56\) and values around \(0.50\)–\(0.53\) at \(x=0.3\), implying increased ductility under the assumptions of the calculation [1706.03252]. Re and Os reduce the formation energy of H and He in W and attract He at substitutional sites, whereas Ta repels He; transmutation elements also alter migration barriers and can increase the diffusion rate of H and He in W [1706.03252]. Coupled with the vacancy-decoration results for Re and Os, this links transmutation chemistry directly to swelling, embrittlement, and bubble evolution [2102.01614, 1706.03252].

Several engineering constraints recur across blanket proposals. One is nuclear-data uncertainty: a \(\pm 10\%\) neutron-flux uncertainty or \(\pm 20\%\) cross-section error produces linear variation in the predicted Au yield in the Hg blanket study, and a separate \(\beta^{-}\)-emitter analysis notes \(10\)–\(20\%\) uncertainties in \(14\ \mathrm{MeV}\) cross sections from ENDF/B-VIII.0 [2507.13461, 2605.20260]. Another is chemical and structural compatibility. Mercury has higher vapor pressure than Pb or Li, must satisfy strict containment because of toxicity and Minamata-Convention controls, and requires qualification of W, V-alloys, and steel in contact with Hg and Li [2507.13461]. The same study identifies \(^{204}\mathrm{Tl}\) (\(t_{1/2}=3.78\ \mathrm{y}\), \(\sim 100\ \mathrm{g/yr}\)) as the longest-lived radionuclide in the activation inventory, with the Au product cooling to \(<\)Class A waste in \(\sim 7\ \mathrm{y}\) and to background-equivalent levels in \(\sim 14\ \mathrm{y}\) [2507.13461].

Qualification sources introduce their own caveats. IFMIF-Lite’s D–Li neutron spectrum is slightly softer below \(14\ \mathrm{MeV}\) than pure D–T, so spectral weighting functions are required when equating damage and transmutation rates to reactor conditions [2302.09011]. For high-activity isotope blankets, chemical separation in in-blanket streams must accommodate dose rates on the order of \(10^7\)–\(10^9\ \mathrm{Ci}\), while tritium handling, permeation control, corrosion, swelling, and activation remain non-trivial [2605.20260].

A persistent misconception is that average or homogenized neutron spectra are adequate for design. The tungsten literature shows the opposite: fine spatial heterogeneity in moderators adjacent to W armor can produce local Re variations of \(\pm 0.2\)–\(0.3\ \mathrm{at.}\%\), and only detailed geometry-plus-flux models reproduce experimental transmutation levels [1612.03892]. A plausible implication is that fusion neutron-driven transmutation should be treated not as a single blanket-average quantity, but as a multiscale problem coupling local spectral shaping, resonance physics, self-shielding, isotope extraction, and defect-mediated microstructural evolution.

Source: https://www.emergentmind.com/topics/fusion-neutron-driven-transmutation