Neutron Reconstruction via Blips in Liquid Argon Time Projection Chambers

This paper presents a simulation-based proof-of-concept demonstrating that isolated MeV-scale energy deposits ("blips") from neutron inelastic scattering in Liquid Argon Time Projection Chambers can be used to identify and reconstruct the direction and energy of final-state neutrons, thereby enhancing neutrino interaction studies such as neutrino-antineutrino separation.

The Gist

Imagine you are trying to solve a giant, invisible 3D puzzle inside a giant tank of super-cold liquid argon. This tank is a Time Projection Chamber (LArTPC), and its job is to catch neutrinos—ghostly particles that pass through the Earth almost without touching anything.

When a neutrino finally hits an argon atom, it explodes into a shower of other particles. Scientists usually track the "loud" and "bright" pieces of this explosion: charged particles like protons and electrons that leave clear, long trails.

But there's a problem. A huge chunk of the energy from the neutrino crash often goes into neutrons. Neutrons are the "ghosts of ghosts." They have no electric charge, so they don't leave a trail. In most current experiments, scientists just ignore them, assuming they aren't there. This is like trying to figure out the weight of a truck by only weighing the driver, ignoring the cargo in the back. You end up with a very wrong answer.

This paper is a "proof-of-concept" study that says: "We can actually see the ghosts, if we know what to look for."

Here is how they do it, explained through simple analogies:

1. The "Blip" (The Ghost's Footprint)

Since neutrons don't leave a trail, how do we see them?
When a fast neutron smashes into an argon atom, it doesn't just bounce off; it excites the atom, like hitting a bell. The atom then "rings" by releasing a tiny burst of energy (a gamma ray). This energy zips around, hits an electron, and creates a tiny, isolated spark of electricity in the detector.

The scientists call these tiny sparks "Blips."

• The Analogy: Imagine a dark room where a ghost walks by. You can't see the ghost, but every time it passes a candle, it flickers the flame for a split second. You can't see the ghost, but you can map its path by looking at where the candles flickered.
• The Reality: The "candles" are the argon atoms, and the "flickers" are the Blips.

2. Cleaning Up the Mess (Filtering the Noise)

The problem is that the detector is noisy. Other things besides neutrons cause flickers (blips).

• The Lepton Noise: Sometimes, the main charged particles (like electrons or muons) kick up dust (called delta-rays) that looks like a blip.
• The Background Noise: Natural radioactivity in the tank (like Argon-39) or cosmic rays from space can also cause flickers.

The paper proposes a set of "rules" to filter out the noise:

• The "Don't Stand Too Close" Rule: If a flicker is right next to the main particle trail, it's probably just dust kicked up by that particle. Ignore it.
• The "Too Small" Rule: If the flicker is too tiny (less than 0.6 MeV), it's probably just background radioactivity. Ignore it.
• The "Too Far" Rule: If the flicker is too far away from the crash site, it's probably a stray cosmic ray. Ignore it.

After applying these rules, the remaining "blips" are very likely to be the footprints of neutrons.

3. Reconstructing the Ghost (What We Learn)

Once they have the clean list of blips, they can do two amazing things:

A. Counting the Ghosts (Identification)
If a neutrino crash produces a lot of blips, it's a strong sign that neutrons were involved. If there are no blips, there were likely no neutrons.

• The Result: They can tell if a neutrino event had neutrons about 70% of the time. This is huge because previously, they were guessing 0% of the time.

B. Guessing the Direction and Energy (Reconstruction)
By looking at where the blips are located, they can draw a line from the crash site to the cluster of blips.

• The Direction: This tells them which way the neutrons were flying.
• The Energy: The total "brightness" (energy) of all the blips combined gives a rough estimate of how much energy the neutrons carried.
• The Analogy: If you see a cluster of broken glass on the floor, you can guess where the vase fell and how hard it hit the ground, even if you didn't see the vase fall.

4. Why Does This Matter? (The "Superpower")

Why go through all this trouble? Because it solves two major physics mysteries:

1. Telling Neutrinos from Anti-Neutrinos
Neutrinos and anti-neutrinos are twins that look almost identical, but they behave differently.

• Neutrinos tend to produce more protons.
• Anti-neutrinos tend to produce more neutrons.
• The Magic: By counting the blips (neutrons), scientists can now tell the difference between a neutrino and an anti-neutrino much better than before. This is crucial for understanding why the universe is made of matter and not antimatter.

2. Fixing the Energy Math
Because neutrons carry away so much energy, ignoring them makes the energy calculations wrong. By adding the "blip energy" back into the math, the scientists get a much more accurate picture of the original neutrino's energy.

The Bottom Line

This paper is like finding a new pair of glasses. For years, scientists looked at neutrino crashes with one eye closed (ignoring neutrons). This study shows that by looking for the tiny, scattered "blips" left behind by neutrons, we can open that second eye.

It's not perfect yet (the resolution is a bit fuzzy, like a low-resolution photo), but it proves the concept works. The authors suggest that with better software (like AI) and more data, we could eventually see these "ghosts" with crystal clarity, unlocking new secrets about the universe.

Technical Summary

1. Problem Statement

Neutrons are significant final-state particles in neutrino interactions, particularly in the sub-GeV energy regime where nuclear effects (Fermi motion, nucleon correlations, final-state interactions) dominate. However, current Liquid Argon Time Projection Chamber (LArTPC) physics analyses largely ignore neutrons due to detection difficulties.

• Detection Challenge: In LArTPCs, low-energy neutrons have a long interaction length (~30 m) due to destructive interference in the neutron-argon cross-section at ~57 keV, making capture on Argon (which produces detectable $\gamma$-rays) unlikely within the detector volume.
• Information Loss: The failure to reconstruct neutrons leads to significant information loss regarding the initial neutrino kinematics (energy and direction), hindering precision measurements for CP violation, mass ordering, and cross-section studies.
• Existing Limitations: Previous attempts to identify neutrons via secondary protons (e.g., in MicroBooNE) have shown low efficiency (~3% at 250 MeV).

2. Methodology

The authors present a simulation-based proof-of-concept study to reconstruct neutron systems using "blips"—isolated, MeV-scale energy deposits resulting from the de-excitation of Argon nuclei after neutron inelastic scattering.

Simulation Framework

• Event Generation: Used FLUKA with the PEANUT model for neutrino-argon interactions (sub-GeV atmospheric flux). Alternative models (GENIE-hA, GENIE-INCL) and transport codes (Geant4 with Bertini/NeutronHP) were used for systematic cross-checks.
• Particle Transport: A two-stage transport was implemented: FLUKA for nuclear interactions and LArSoft/Geant4 for electron transport to generate "truth-level" blips.
• Blip Response: Realistic detector response was applied using MicroBooNE blip response matrices (smearing true energy and applying efficiency losses) to simulate reconstructed blip energies and positions.

Signal Selection and Background Rejection

To isolate neutron-induced blips, the study applied a series of topological and energy cuts:

1. Lepton Rejection:
   • $\nu_e/\bar{\nu}_e$: Removed blips within a $45^\circ$ cone of the electromagnetic shower.
   • $\nu_\mu/\bar{\nu}_\mu$: Applied a Point of Closest Approach (PoCA) cut, requiring blips to be $>20$ cm from the muon track.
2. Vertex and Capture Rejection:
   • Removed blips $<20$ cm from the primary neutrino interaction vertex (to reject primary de-excitation $\gamma$-rays).
   • Removed blips $<40$ cm from the muon track endpoint (to reject $\mu^-$ capture-at-rest products).
3. Energy Cut: Applied a minimum reconstructed energy of 0.6 MeV to suppress the dominant background from intrinsic $^{39}$Ar $\beta$-decays (endpoint 0.56 MeV).
4. Fiducial Cut: Restricted blips to within 2.5 m of the interaction vertex to minimize cosmogenic and surface radioactivity backgrounds.

Reconstruction Strategy

Instead of reconstructing individual neutrons (which is difficult due to the spread of blips from neutron free paths), the study reconstructs the neutron system (the ensemble of all primary neutrons):

• Identification: Based on blip multiplicity. Events with more blips are classified as "neutron-rich" ($Xn$), while those with fewer are "neutron-poor" ($0n$).
• Direction Reconstruction: The system direction is estimated by summing vectors from the primary vertex to each selected blip ($\vec{V}_b = \sum \vec{V}_j$).
• Energy Reconstruction: The system energy is estimated by the naive sum of reconstructed blip energies ($E_{blip} = \sum E_j$).

3. Key Contributions

• First Quantification of Directional Reconstruction: This is the first study to demonstrate the feasibility of reconstructing the direction of a neutron system in a LArTPC using blip topology.
• System-Level Approach: Proposes treating neutrons as a collective system rather than individual particles, bypassing the difficulty of separating overlapping neutron signatures.
• Physics Application Framework: Demonstrates how blip-based neutron information can be integrated into standard neutrino analysis workflows to improve specific physics goals.

4. Key Results

Neutron Identification Performance

• Using a simple cut on blip multiplicity (optimized at >3 blips), the study achieves:
  • Efficiency: ~70% for identifying events with primary neutrons.
  • Purity: ~80% (meaning ~20% of identified "neutron-rich" events are actually neutron-free, contaminated by non-primary neutrons or backgrounds).
  • Discrimination Power: The Area Under the Curve (AUC) for ROC analysis ranges from 0.82 to 0.85 across different neutrino flavors.

Kinematic Reconstruction

• Direction: At 100 MeV neutron kinetic energy, the angular resolution is approximately $40^\circ$. This is sufficient for differential cross-section measurements.
• Energy: The energy resolution is approximately 50% at 100 MeV. The reconstructed energy accounts for only ~50% of the true neutron energy due to losses in elastic scattering, binding energy, and sub-threshold charge. However, a clear correlation exists between blip sum and true energy.

Physics Applications

• Neutrino-Antineutrino Separation:
  • Since antineutrinos produce more neutrons than neutrinos, blip information significantly enhances sample purity.
  • For atmospheric $\bar{\nu}_e$, purity improved from 28% (using proton cuts alone) to 43% (adding blip cuts), with 49% efficiency.
  • For BNB Reverse Horn Current (RHC) beams, $\bar{\nu}$ purity reached 68% in the $0pXn$ category.
• Kinematic Reconstruction: Adding neutron system information to standard lepton/hadron reconstruction reduces the bias in neutrino energy estimation and improves angular resolution (e.g., improving direction resolution from $\sim100^\circ$ to $\sim90^\circ$ for sub-300 MeV neutrinos).

5. Significance and Future Outlook

• Impact on Precision Era: This study provides a pathway to recover missing energy and momentum in sub-GeV neutrino interactions, which is critical for next-generation oscillation experiments (like DUNE and Hyper-K) aiming for precision CP violation measurements.
• Scalability: The current results are based on simple cuts and truth-level analogs. The authors argue that performance will improve significantly with:
  • Advanced AI/ML techniques to exploit multi-dimensional correlations between blips and event topology.
  • Better modeling of neutron-argon interactions and nuclear de-excitation.
  • Utilization of scintillation light for timing and background rejection.
• Systematic Uncertainties: The study identifies modeling uncertainties (FLUKA vs. Geant4, GENIE variants) as a major challenge, noting ~5% differences in efficiency and ~25% in background rejection. Future work requires dedicated experimental data to tune these models.

Conclusion: The paper successfully demonstrates that "blips" in LArTPCs are a viable and powerful tool for neutron reconstruction. Even with a simplistic approach, it offers substantial improvements in neutrino-antineutrino separation and kinematic reconstruction, laying the groundwork for more sophisticated, AI-driven neutron analysis in future LArTPC experiments.