Ultra-cold neutron simulation framework for the free neutron lifetime experiment $τ$SPECT

To address the neutron lifetime puzzle and characterize systematic uncertainties in the fully magnetic trap of the $\tau$SPECT experiment, the authors developed a dedicated simulation framework based on the PENTrack software, enhanced with tools for flexible configuration and data visualization, which successfully aligns with experimental data from the Paul Scherrer Institute.

The Gist

The Big Mystery: The "Ghost" Neutron

Imagine you have a bag of marbles that are slowly disappearing. You want to know exactly how long it takes for half of them to vanish. This is the free neutron lifetime.

Neutrons are the building blocks of atoms. When they are stuck inside an atom, they are stable. But when they are floating alone in space (a "free" neutron), they eventually decay into other particles. Scientists have been trying to time this decay for decades.

The Problem: There are two main ways to time these marbles, and they give different answers.

1. The "Beam" Method: Counting how many marbles fall out of a moving stream.
2. The "Bottle" Method: Putting marbles in a jar and counting how many are left after an hour.

Currently, these two methods disagree by about 10 seconds. This is a huge mystery in physics, often called the "Neutron Lifetime Puzzle." If the answer is wrong, it might mean our understanding of the universe's beginning (the Big Bang) or the fundamental rules of particle physics is broken.

The New Experiment: τSPECT

To solve this, a team of scientists built a new experiment called τSPECT. Instead of a glass jar, they use a magnetic bottle.

Think of neutrons like tiny, invisible magnets. The τSPECT experiment uses powerful magnets to create a "force field" trap. Because the neutrons are so cold (moving very slowly, like a sleepy snail), they bounce off this magnetic force field instead of hitting the walls. This avoids the problem of neutrons getting stuck or absorbed by the physical walls of a jar, which was causing errors in previous experiments.

The Challenge: Simulating the Invisible

You can't just watch a single neutron. It's too small, too fast, and too quiet. To understand how the experiment works, the scientists needed a virtual reality simulator.

They built a digital twin of the entire experiment. Imagine a video game where you can spawn thousands of invisible, super-slow marbles and watch how they bounce around inside a magnetic force field.

The "Video Game" Engine (PENTrack):
The core of their simulation is a piece of software called PENTrack. It's like a physics engine in a video game, but instead of calculating how a car crashes, it calculates how a neutron bounces off a magnetic field, spins, and eventually decays.

The Two "Helper" Tools

The scientists realized that running this simulation was like trying to drive a race car without a steering wheel or a dashboard. So, they built two custom tools to make it easier:

1. The "Remote Control" (penconf):
   Imagine you want to test what happens if you raise the floor of the experiment by 10 centimeters, or if you turn the magnets up to 11. Instead of rewriting the entire computer code every time, this tool lets you just type in the new settings (like changing the speed on a thermostat). It automatically sets up the simulation for you.

2. The "Movie Maker" (penplot):
   Once the simulation runs, you get millions of lines of numbers. That's boring. This tool takes those numbers and turns them into 3D movies, graphs, and animations. It lets the scientists see the neutrons bouncing around, so they can spot where things go wrong.

How the Simulation Works (The Story of a Neutron)

Here is the journey of a simulated neutron in their framework:

1. Birth: The neutron is born in a giant nuclear reactor (the PSI source in Switzerland). It's hot and fast.
2. Cooling Down: It travels through a "mirror tunnel" (a beamline) where it cools down until it's moving slower than a walking human (under 8 meters per second).
3. The Spin Flip: This is the tricky part. Neutrons have a "spin" (like a top spinning). To get trapped in the magnetic bottle, they need to be spinning in the right direction. The simulation includes special "spin-flipper" units (like a DJ spinning a record) that flip the neutron's spin so it can be caught by the magnetic trap.
4. The Trap: The neutron enters the magnetic bottle. It bounces around inside the invisible magnetic walls.
5. The Cleanup: Sometimes, a few "rogue" neutrons get trapped with too much energy. They might eventually escape, making the scientists think the neutron decayed when it actually just ran away. The simulation helps design a "cleaning" phase where a detector scoops up these energetic rogues before the real timing starts.
6. The Count: After a set time, a detector goes into the trap and counts how many neutrons are left.

The "Aha!" Moment

The scientists compared their computer simulation with the real experiment.

• The Good News: The simulation matched the real data almost perfectly for most things. It proved the magnetic trap works as designed.
• The Bad News: The simulation predicted more neutrons would survive than actually did.
• The Detective Work: By tweaking the simulation, they realized the "copper shield" inside the real experiment was acting like a sticky trap, absorbing neutrons it was supposed to bounce off. This discovery tells them they need to polish that shield to make the experiment more accurate.

Why Does This Matter?

This paper isn't just about writing code; it's about building a time machine for data.

By perfecting this simulation, the scientists can:

• Find the errors: They can spot tiny systematic mistakes (like the sticky copper) before they ruin the experiment.
• Solve the puzzle: Once they fix these errors, they hope to get a precise measurement of the neutron's life.
• Understand the Universe: If they solve the "Neutron Lifetime Puzzle," it could tell us why the universe is made of matter instead of just energy, and how the first stars were born.

In short: They built a super-accurate digital twin of a neutron trap to help them catch the "ghost" of the universe's most elusive particle, ensuring that when they finally time its life, the answer is the truth.

Technical Summary

1. Problem Statement

The precise determination of the free neutron lifetime ($\tau_n$) is a critical observable in particle physics and cosmology. It is essential for:

• Determining the first element of the Cabibbo–Kobayashi–Maskawa (CKM) quark mixing matrix ($V_{ud}$).
• Calculating primordial element abundances from Big-Bang Nucleosynthesis.

Currently, two primary measurement techniques—the beam method and the bottle method—yield incompatible results, creating the "neutron lifetime puzzle." The bottle method, which uses Ultra-Cold Neutrons (UCNs) trapped in a volume, is susceptible to systematic uncertainties caused by neutron interactions with material walls (absorption and up-scattering).

The $\tau$SPECT experiment aims to resolve this by using a fully magnetic trap to confine neutrons, thereby eliminating material wall interactions. To achieve a target sensitivity of $\sigma_{\tau_n} \leq 0.3$ s, the experiment requires a rigorous understanding of UCN dynamics, storage efficiency, and systematic error sources.

2. Methodology

The authors developed a dedicated end-to-end simulation framework to model the entire lifecycle of UCNs in the $\tau$SPECT experiment, from production at the PSI source to detection.

Core Software Architecture

The framework is built upon the existing C++ Monte Carlo software PENTrack, which simulates UCN dynamics in complex geometries and electromagnetic fields. It was significantly enhanced with two custom Python companion tools:

1. penconf: A configuration tool that allows for flexible, parametrizable upstream generation of PENTrack input files. It handles complex experimental settings (e.g., shutter timing, spin-flipper positions, detector speeds) and manages moving geometries by activating/deactivating pre-defined CAD volumes.
2. penplot: A downstream tool for analyzing, visualizing, and animating simulation data, including 3D plots and trajectory animations.

Key Simulation Components

• UCN Source Modeling: The framework integrates the Paul Scherrer Institute (PSI) UCN source. To optimize computational efficiency, a "Virtual Source" technique was implemented. Instead of simulating neutrons from the solid deuterium target every time, a global simulation records neutron properties at a virtual detector plane (after the beamport). Subsequent $\tau$SPECT-specific simulations sample from this distribution, drastically reducing computation time while maintaining statistical accuracy.
• Magnetic Trap Simulation:
  • Fields: Magnetic fields are generated externally using magpylib and imported as tricubic interpolation meshes into PENTrack.
  • Components: The trap consists of superconducting coils (longitudinal field) and a permanent Halbach octupole magnet (transverse field).
  • Potential: The simulation calculates the combined magnetic and gravitational potential ($V_{trap} \approx 48$ neV) to determine which neutrons are storable.
• Spin-Flipping Units (SFUs): The framework models two adiabatic spin-flipping units located at specific longitudinal positions.
  • Single Spin-Flip (sSF): Converts High Field Seekers (HFS) to Low Field Seekers (LFS), increasing total energy by $\approx +25$ neV.
  • Double Spin-Flip (dSF): Converts LFS $\to$ HFS $\to$ LFS, reducing total energy by $\approx -96$ neV.
  • Efficiency: Dedicated sub-simulations calculate spin-flip efficiency based on coil geometry and RF current, which are then applied as probabilities in the main simulation to avoid computationally expensive real-time spin tracking.
• Cleaning and Detection: The simulation models the "cleaning" phase where a detector absorbs marginally trapped neutrons (those with energy slightly above $V_{trap}$) before the storage phase begins.

3. Key Contributions

• Integrated Framework: Creation of a unified simulation environment linking the PSI UCN source, the beamline, the magnetic trap, and the detection system.
• Virtual Source Technique: Development of a sampling method that decouples the source simulation from the experiment simulation, enabling high-statistics runs without the computational cost of tracing neutrons from the source crystal.
• Dynamic Geometry Handling: Implementation of a method to simulate moving components (spin-flippers, shutters, detectors) by toggling the state of static CAD volumes, allowing for the modeling of complex experimental cycles (filling, cleaning, storing, counting).
• Systematic Uncertainty Analysis: The framework allows for the isolation and study of specific loss mechanisms, such as surface roughness, material absorption, and the population of marginally trapped neutrons.

4. Results

• Validation against Data: The simulation results show excellent agreement with experimental data from $\tau$SPECT at PSI.
  • Beamline Monitoring: Simulated UCN arrival times and storage time constants ($\tau_s \approx 37.8$ s) match measured values ($\approx 37.2$ s) within uncertainties.
  • Filling Time Scan: The simulation correctly predicts an optimal filling time of $\approx 25$ s, matching experimental scans.
• Discrepancy Identification: When using nominal material properties, the simulation predicted nearly three times more stored neutrons than observed experimentally.
  • Diagnosis: By adjusting the loss parameters of the copper shield (increasing the imaginary part of the potential, $W$, to simulate high loss), the simulation amplitude matched the experimental data.
  • Implication: This suggests the copper shield has higher UCN loss rates than expected, likely due to surface roughness or contamination, highlighting a critical systematic uncertainty.
• Energy Spectrum Analysis: The simulation successfully models the energy distribution of trapped neutrons, distinguishing between fully storable neutrons and marginally trapped ones that leak out over time.

5. Significance

• Resolving the Neutron Lifetime Puzzle: By providing a high-fidelity tool to identify and quantify systematic uncertainties (specifically material losses and marginally trapped neutrons), this framework is essential for $\tau$SPECT to achieve its goal of an independent, high-precision measurement of $\tau_n$.
• Experimental Optimization: The framework enables "virtual commissioning," allowing researchers to test new hardware configurations, spin-flipping strategies, and cleaning protocols before physical implementation.
• Community Resource: The authors plan to release a clean version of the framework on GitLab, providing a valuable resource for other UCN experiments to model their own setups and reduce systematic errors.
• Future Directions: The framework will be used to study depolarization rates, field inhomogeneities, and source fluctuations, further refining the precision of the neutron lifetime measurement.