FLUKA-Based Optimization of Muon Production Target Design for a Muon Collider Demonstrator

This study utilizes FLUKA simulations to optimize the geometry and material of an 8 GeV proton beam target within a 5 T solenoid for a muon collider demonstrator, balancing secondary particle yield and emittance against thermal limitations to guide future front-end system designs.

The Gist

Imagine you are trying to build a high-speed train station, but instead of regular trains, you are trying to catch muons. Muons are like tiny, super-fast, and very shy particles that are incredibly hard to catch. They are the "gold" of future particle physics, but they only exist for a split second before they vanish.

To get them, you have to smash a beam of protons (like a high-speed cannonball) into a solid block of material called a target. This collision creates a shower of new particles, including pions, which quickly decay into the muons we want.

The problem? It's a messy explosion. The muons fly off in all directions, and many of them get lost before we can grab them. This paper is essentially a blueprint for designing the perfect "catching net" (the target and the magnetic field) to grab as many muons as possible without melting the net.

Here is the breakdown of the research using simple analogies:

1. The Setup: The Proton Cannon and the Magnetic Net

Think of the experiment as a giant slingshot.

• The Cannon: An 8 GeV proton beam (a stream of super-fast protons) hits a target.
• The Net: Surrounding the target is a massive magnet (a solenoid) that acts like a funnel. It's 2 meters long and has a magnetic field as strong as 5 Tesla (about 100,000 times stronger than a fridge magnet). Its job is to grab the chaotic spray of particles and straighten them out into a neat beam.

2. The Challenge: The "Black Box" Problem

The researchers used a super-computer simulation called FLUKA to predict what would happen. However, the standard software was like a basic weather app: it could tell you "it's raining," but it couldn't tell you exactly how hard the drops were hitting the ground or where they were landing.

To get the detailed data they needed, the author had to write custom code (like writing your own weather sensor app). This allowed them to track exactly where every single particle went and how much energy it carried, which is crucial for designing the machine.

3. The Experiment: Tweaking the Target

The team asked two main questions: "What shape should the target be?" and "What should it be made of?"

A. The Shape (Geometry)

Imagine you are trying to catch water from a fire hose with a bucket.

• The Bucket Width (Radius): They tested buckets of different widths. They found that making the bucket wider didn't really help catch more water; it just made the splash zone a bit bigger. The "messiness" (emittance) of the beam stayed about the same.
• The Bucket Depth (Length): They tested buckets of different depths. A longer bucket (target) caught slightly more water and organized it better, but it also made the water take longer to get through.
• The Heat: The biggest worry was heat. When the proton cannon hits the target, it gets hot—like a car engine running too hard. They found that while the shape changed the heat distribution slightly, the area right where the beam hits gets hot no matter what.

B. The Material (The "Bucket" Substance)

They tested six different materials to see which one worked best. Think of this as choosing between a wooden bucket, a steel bucket, or a plastic bucket.

• Inconel (The Heavy Metal): This material was the star player. It produced the most muons and kept the beam fairly organized. It's like a bucket that catches the most water without breaking.
• Beryllium (The Light Feather): This material stayed the coolest. Because it is very light and not dense, the proton cannon didn't hit as many atoms, so less heat was generated. It's like a bucket that barely warms up, but it doesn't catch as many muons as Inconel.
• Tungsten (The Heavy Hitter): This produced a lot of neutrons (a type of radiation), which is like having a lot of dangerous debris flying around that you don't want.

4. The "Melting" Warning

The simulation showed how hot the target would get. However, the computer model (FLUKA) is a bit like a safety alarm that goes off too early. It predicts the heat will be higher than it actually will be because it doesn't account for the material cooling itself down (like a fan blowing on a hot engine). So, while the numbers look scary, the real-world target might actually survive better than the computer predicts.

The Bottom Line

This paper is the first step in building a Muon Collider. The researchers learned that:

1. Shape matters, but not as much as you think: Making the target slightly longer helps organize the beam, but making it wider doesn't help much.
2. Material is a trade-off: You have to choose between Inconel (great at making muons, but gets hot) and Beryllium (stays cool, but makes fewer muons).
3. The Future: Now that they have these "intuition" maps, the next step is to build a real target that can survive the heat and the radiation, likely using a mix of these ideas to create a machine that can finally collide muons at high energies.

In short: They figured out the best way to aim a proton cannon at a block of metal to harvest the most valuable particles without melting the block, paving the way for the next generation of physics experiments.

Technical Summary

1. Problem Statement

Muon colliders represent a promising frontier in high-energy physics, offering the potential for high-energy collisions with minimal synchrotron radiation. However, their feasibility is bottlenecked by the efficiency of muon production. Since muons are generated via the decay of pions, which are created when a high-energy proton beam strikes a fixed target, the design of this target is critical.

The core challenges addressed in this study are:

• Optimization Trade-off: Maximizing the yield of secondary particles (pions and muons) while minimizing beam emittance (to ensure beam quality).
• Thermal Survivability: Managing the intense heat deposition from an 8 GeV proton beam to prevent target failure and extend operational lifespan.
• Simulation Limitations: Standard simulation tools often lack the granularity to extract specific spatial and momentum distributions of particles exiting a target, requiring custom solutions.

2. Methodology

The study utilizes the FLUKA Monte Carlo particle transport code to simulate an 8 GeV proton beam striking a target positioned at the center of a 2-meter long, 5 Tesla solenoid. The solenoid is designed to capture secondary particles with high transverse momentum.

Simulation Setup

• Beam: 8 GeV proton beam.
• Statistics: 100,000 primary protons per simulation run.
• Target Parameters: Varied in radius (0.3–2.1 cm), length (39–78 cm), and material (Graphite, Beryllium, Inconel, Tungsten, etc.).
• Scoring: Yield of $\pi^+$ and $\mu^+$, beam emittance, and temperature rise ($\Delta T$) per bunch (assuming $10^{13}$ protons/bunch).

Technical Innovations & Custom Routines

A significant portion of the methodology involved overcoming limitations in the FLUKA graphical interface (Flair):

1. Custom User Routines: Standard scoring cards were insufficient for detailed spatial/momentum analysis. The author implemented custom Fortran routines:
   • mgdraw.f and fluscw.f: To extract detailed particle information (spatial distribution and momentum) exiting the target.
   • magfld.f and source.f: To define custom magnetic fields and particle sources.
2. Magnetic Field Modeling: Since Flair cannot natively define complex solenoidal fields, two approaches were tested:
   • Analytical Approximation: Using magfld.f to implement a Biot-Savart derived axial field approximation. This was accurate near the axis but degraded radially.
   • Field Map Import: Generating a high-resolution field map (>4,000 points) using G4beamline and importing it into FLUKA via the MGNDATA card. This method successfully accounted for fringe fields and provided a more accurate, flexible field representation.

3. Key Contributions

• Development of Custom FLUKA Workflows: Demonstrated the necessity and implementation of user routines (mgdraw, fluscw) to extract specific phase-space data beyond default capabilities.
• Hybrid Magnetic Field Integration: Validated a workflow combining G4beamline for field generation and FLUKA for particle transport, solving the issue of importing complex 3D magnetic maps into FLUKA.
• Systematic Parametric Study: Provided a comprehensive dataset on how target geometry (radius/length) and material properties (atomic number $Z$, density) simultaneously impact particle yield, beam emittance, and thermal load.

4. Key Results

A. Geometry Optimization (Graphite Target)

• Radius Variation (0.3 – 2.1 cm):
  • Yield/Emittance: Increasing the target radius resulted in only modest increases in beam spot size. The emittance remained relatively constant, with variations falling within statistical fluctuations.
  • Temperature: Temperature rise near the beam path was nearly identical across radii. Differences only appeared far from the beam trajectory.
• Length Variation (39 cm vs. 78 cm):
  • Yield/Emittance: Longer targets (78 cm) produced a more centrally concentrated phase-space distribution, resulting in slightly lower emittance compared to shorter targets. However, shorter targets produced smaller beam spot sizes.
  • Temperature: Temperature rise near the beam entry was consistent regardless of length.

B. Material Optimization

Six materials were tested with lengths corresponding to two interaction lengths for each material.

• Yield and Emittance:
  • Pion and muon spot sizes and phase-space distributions were nearly identical across all materials.
  • Inconel emerged as the superior performer, showing a noticeably higher pion and muon yield with relatively low emittance compared to other materials.
• Thermal Performance:
  • Beryllium exhibited the smallest temperature rise due to its low density, resulting in fewer proton-proton interactions and lower energy deposition.
  • Tungsten (High-Z): While not explicitly detailed in the yield summary, the paper notes that high-Z targets produce significantly more neutrons, a critical factor for shielding and activation.

C. Thermal Simulation Limitations

The study notes that FLUKA calculates energy deposition but does not simulate fluid dynamics or thermodynamic heat transfer. Therefore, the reported temperature rise ($\Delta T$) serves as a reasonable upper bound. Future work requires coupling FLUKA with dedicated thermal-fluid solvers (e.g., ANSYS) for accurate structural survivability analysis.

5. Significance and Conclusion

This work establishes a foundational framework for the front-end design of a muon collider demonstrator.

• Design Intuition: It provides initial intuition that while geometric changes (radius/length) have modest effects on emittance, they influence beam compactness and time structure.
• Material Selection: The study highlights a critical trade-off: Inconel offers the best particle production efficiency, while Beryllium offers the best thermal survivability.
• Future Roadmap: The paper concludes that while FLUKA is excellent for particle interactions and providing thermal upper bounds, the next phase of design must integrate dedicated thermal and structural simulations to define realistic operating limits and ensure the target can survive the extreme conditions of a muon collider.

The research successfully bridges the gap between theoretical particle production and practical engineering constraints, guiding the optimization of target systems for next-generation lepton colliders.