From 06b1bbf2bfc31e2204bc772a0f8717b3010dfca5 Mon Sep 17 00:00:00 2001
From: Nan Shi <shinan@fusion.gat.com>
Date: Thu, 4 Dec 2025 11:30:26 -0800
Subject: [PATCH] Implement new fast thermalization source

---
 src/physics/nuclear.jl | 7 +++++++
 1 file changed, 7 insertions(+)

diff --git a/src/physics/nuclear.jl b/src/physics/nuclear.jl
index 576c1f67..62f7c201 100644
--- a/src/physics/nuclear.jl
+++ b/src/physics/nuclear.jl
@@ -240,6 +240,13 @@ function fusion_reaction_source(
     ion.particles = reactivity
     ion.fast_particles_energy = eV
 
+    if out == :He4
+        k += 1
+        ion = resize!(s1d.ion, k)[k]
+        ion_element!(ion, out; fast=false)
+        ion.particles = reactivity
+    end
+
     return nothing
 end