Geant4 and FLUKA Simulations of a Cyclotron Based 30 MeV Proton-Beryllium Reaction: Benchmarking and Optimization of Neutron Fields

This paper presents a comparative simulation study using Geant4 and FLUKA to benchmark and optimize a 30 MeV proton-induced beryllium neutron source, evaluating neutron fluence, gamma dose, and moderation effects to guide the design of a modular irradiation station for thermal neutron fields.

The Gist

Imagine you are a chef trying to bake a very specific kind of cake. You need a steady stream of "neutron sprinkles" to make it work. Usually, big nuclear reactors are the industrial ovens that provide these sprinkles in massive quantities. But sometimes, you don't need a whole factory; you just need a small, precise pinch for a delicate experiment. You need a "kitchen-scale" neutron source.

This paper is about building and testing a new, custom "neutron sprinkler" using a particle accelerator (a machine that shoots tiny particles at high speeds) instead of a giant reactor.

Here is the story of how they built it, explained simply:

1. The Recipe: Shooting Protons at Beryllium

The team wanted to create neutrons by shooting a stream of protons (tiny, positively charged particles) at a block of Beryllium (a light, silvery metal).

• The Collision: When the protons hit the beryllium, it's like throwing a fast-moving billiard ball at a stack of other balls. The collision knocks neutrons loose, creating a spray of them.
• The Angle: They found that tilting the beryllium block at a 45-degree angle was the sweet spot. Think of it like angling a garden hose; if you point it straight at the ground, the water splashes everywhere. If you angle it just right, you get a better, more focused spray.
• The Safety Net: Since beryllium is toxic and gets very hot, they sandwiched it between a heat sink (to cool it down) and a thick lead shield (to block dangerous gamma rays, which are like invisible, high-energy X-rays).

2. The Simulation: Two Different GPS Systems

Before building the real machine, they had to predict exactly what would happen. To do this, they used two powerful computer programs: Geant4 and FLUKA.

• The Analogy: Imagine you are planning a road trip. You ask two different GPS apps (like Google Maps and Waze) for directions. Both use real data, but they have different algorithms and map styles. Sometimes they agree perfectly; other times, one says "turn left" and the other says "turn right."
• The Goal: The authors wanted to see if these two "GPS apps" would give them the same map for their neutron sprinkler. They ran thousands of virtual simulations to see how many neutrons would be produced, where they would go, and how much radiation dose they would create.
• The Result: The two programs agreed very well on the "low-energy" neutrons (the slow, gentle ones). However, for the "high-energy" neutrons (the fast, aggressive ones), FLUKA predicted slightly more than Geant4. This is important because it tells scientists to be careful when interpreting data from high-speed particles.

3. The Filter: Turning Fast Neutrons into Slow Ones

The neutrons coming out of the collision are like a chaotic crowd of runners sprinting in all directions. For many experiments (like testing how materials react to radiation), you need "thermal neutrons"—these are the slow, lazy walkers that move gently.

• The Analogy: To slow the runners down, the team built a giant moderator out of High-Density Polyethylene (HDPE). Think of this material as a giant, thick sponge or a dense forest.
• The Process: When the fast neutrons run into the plastic "forest," they bounce off the hydrogen atoms inside the plastic, losing speed with every bump.
• The Optimization: They tested different thicknesses of this plastic "forest." They found that wrapping the target completely in a thick block of plastic (12 cm) was much better than just putting a slab in front of it. It acted like a reflector, bouncing neutrons back into the stream and slowing them down efficiently. This boosted the number of useful "slow neutrons" by nearly double.

4. The Final Beam: A Gaussian Cloud

After all the filtering and slowing down, they looked at the final beam of neutrons.

• The Shape: Instead of a laser-like straight line, the beam spread out in a wide, soft cloud. It followed a Gaussian distribution, which is a fancy way of saying it looked like a bell curve.
• The Size: The beam was quite wide (about 30 cm across), which is great for experiments that need to cover a large area, like testing electronic boards for space missions (to see if cosmic rays will flip their switches).

Why Does This Matter?

This paper is a "user manual" for scientists who want to build their own neutron sources without needing a massive nuclear reactor.

• Safety & Cost: It shows how to use a standard hospital-style cyclotron (usually used for making medical isotopes) to create a safe, controllable neutron source.
• Reliability: By comparing the two computer programs, they gave future researchers a guide on how to trust their data.
• Versatility: They proved that by tweaking the angle of the target and the thickness of the plastic moderator, you can tune the machine to produce exactly the type of neutron field you need for your specific experiment.

In short, the team built a virtual prototype of a "neutron sprinkler," tested it with two different simulation engines to make sure the math was right, and figured out the best way to turn a chaotic spray of fast particles into a gentle, useful stream of slow neutrons.

Technical Summary

1. Problem Statement

The research addresses the need for reliable, alternative neutron sources for applications such as material hardness studies, neutron radiography, activation analysis, and Single Event Upset (SEU) testing for electronics in space environments. While nuclear reactors provide high flux, they are often unsuitable for low-flux research or facilities lacking reactor infrastructure.

The specific challenge is to optimize a compact, accelerator-driven neutron source using the $^{9}\text{Be}(p,n)^{9}\text{B}$ reaction driven by a 30 MeV proton cyclotron (IBA Cyclone 30 XP). A critical gap exists in understanding the discrepancies between major Monte Carlo simulation toolkits (Geant4 and FLUKA) regarding neutron production, transport, and moderation. Without benchmarking these differences, experimental design and safety assessments (shielding, dose) may be inaccurate.

2. Methodology

The authors employed a comparative simulation approach using Geant4 and FLUKA to model a proton-beryllium irradiation setup.

• Target Configuration:

  • Material: 100% pure Beryllium ($^{9}\text{Be}$).
  • Geometry: A 15 mm $\times$ 15 mm $\times$ 3 mm target.
  • Orientation: Rotated at a 45° angle relative to the beam axis to minimize backscattering and maximize forward yield.
  • Cooling: The target rests on an aluminum beam dump to facilitate water cooling and prevent beryllium blistering (hydrogen buildup). The effective proton path length was calculated to ensure the target thickness (3 mm) remains below the blistering threshold (5.7 mm range) even with the 45° tilt.
  • Shielding: The setup includes a 20 mm aluminum beam dump, a 5 cm lead shell for gamma attenuation, and a 1.5 m concrete wall for the irradiation chamber.

• Simulation Parameters:

  • Beam: 30 MeV protons with a 1 cm diameter and Gaussian energy spread ($\sigma = 500$ keV/c).
  • Physics Lists:
    • Geant4: Utilized the QGSP_BIC_HP physics list for inelastic interactions and HadronElasticPhysicsHP for elastic scattering, coupled with the G4NDL 3.16 library (ENDF/B-VII) for high-precision thermal neutron transport (< 19.5 MeV).
    • FLUKA: Used version 4-3.3 with the PRECISION preset, COALESCE (for light cluster formation), DECAYS, and EVAPORAT cards to model evaporation and decay processes.
  • Moderation: High-Density Polyethylene (HDPE) was used to thermalize neutrons. Various configurations were tested, including a 12 cm slab versus a 12 cm encapsulating box.
  • Analysis: Neutron fluence, angular distribution, energy spectra, and dose equivalents were scored using detectors (e.g., USRBDX, USRTRACK, USRBIN) at a probe point 70 cm downstream.

3. Key Contributions

• Toolkit Benchmarking: A rigorous comparison of Geant4 and FLUKA outputs for the $^{9}\text{Be}(p,n)$ reaction, establishing that while results agree at low energies, deviations occur at higher energies.
• Optimization of Target Geometry: Determination that a 45° tilt yields the highest neutron fluence compared to perpendicular or other angles.
• Moderation Strategy: Demonstration that encapsulating the target with HDPE (rather than placing a slab in the beam path) significantly increases the thermal neutron yield due to neutron reflection properties.
• Safety and Design Guidelines: Calculation of optimal target thickness to avoid hydrogen blistering and the design of a modular irradiation station capable of generating thermal neutron fields.

4. Key Results

• Target Orientation: The 45° tilt resulted in the highest total response ($3.3 \times 10^{-5}$ neutrons/cm$^2$/proton), outperforming the perpendicular (0°) orientation ($2.2 \times 10^{-5}$).
• Angular Distribution: Neutrons exhibited a quasi-isotropic distribution with a peak around $\theta \approx 70^\circ$ relative to the beam axis. High-energy neutrons showed more collimation, while low-energy neutrons traveled parallel to coordinate axes.
• Moderation Efficiency:
  • A 12 cm HDPE encapsulation achieved a thermal-to-total neutron ratio of 37.6%, compared to only 21.2% for a 12 cm slab.
  • The thermal neutron yield (encapsulation) was $4.46 \times 10^{-7}$ per primary proton.
• Toolkit Comparison:
  • Low Energy (< 1 MeV): Geant4 and FLUKA showed strong agreement in spectral fluence.
  • High Energy: FLUKA predicted higher neutron fluence than Geant4.
  • Thermal Neutrons: Both toolkits agreed on the thermalization efficiency when using the HP physics lists in Geant4 and the PRECISION settings in FLUKA.
• Beam Profile: The resulting neutron beam followed a Gaussian profile with a Full Width at Half Maximum (FWHM) significantly larger than the incident proton beam, spanning approximately 100 cm in the YZ plane.
• Dose Assessment: Prompt gamma and neutron dose equivalents were mapped, confirming the necessity of the lead and concrete shielding for safe operation.

5. Significance

This work provides a validated framework for designing accelerator-based neutron sources, specifically for facilities utilizing the IBA Cyclone 30 XP. By benchmarking Geant4 and FLUKA, the authors provide researchers with a guide to interpret simulation data accurately, acknowledging that slight variations in physics models can lead to different outcomes.

The study confirms that Beryllium is a viable, high-conductivity target material for 30 MeV protons, provided safety protocols for its toxicity and blistering risks are managed. The optimized configuration (45° tilt + HDPE encapsulation) offers a robust method for generating thermal neutron fields, enabling diverse research applications from material science to radiation hardening of electronics without the need for a nuclear reactor.