Modeling the light response of an optically readout GEM based TPC for the CYGNO experiment

This paper presents a model developed by the CYGNO collaboration that accurately predicts the gain-reduction effect caused by space-charge buildup in optically readout GEM-based TPCs, validated by experimental data from a two-liter prototype with percent-level precision.

The Gist

Imagine you are trying to catch a very shy, tiny ghost (a dark matter particle or a neutrino) that only leaves a faint, fleeting footprint in a room filled with invisible gas. To see this ghost, you need a camera so powerful it can count individual photons of light, and a trap so sensitive it can amplify a single whisper into a shout.

This paper is about building and understanding that super-sensitive trap for the CYGNO experiment. Here is the story of how they did it, broken down into simple concepts.

1. The Goal: Catching Ghosts in a Gas Room

The scientists are building a giant, cubic-meter-sized box filled with a special gas mixture (Helium and Tetrafluoromethane). When a mysterious particle from space bumps into the gas, it knocks electrons loose, creating a tiny trail of ionization.

The problem? That trail is incredibly faint. To see it, they need to amplify it. They use a technology called GEMs (Gas Electron Multipliers). Think of a GEM as a microscopic honeycomb made of metal. When an electron enters a hole in the honeycomb, it gets zapped by electricity, creating a cascade of thousands of new electrons. It's like a snowball rolling down a hill, picking up more snow until it becomes an avalanche.

2. The Camera: Seeing the Light

Usually, detectors count the electrons directly. But the CYGNO team is doing something different: they are taking pictures of the light produced by the avalanche.

• The Analogy: Imagine the avalanche of electrons is a crowd of people running through a tunnel. As they run, they rub against the walls and create sparks (light). Instead of counting the people, the scientists are using a super-sensitive camera (a scientific CMOS sensor) to take a photo of the sparks.
• Why? This allows them to see the shape of the trail with incredible detail, helping them figure out exactly where the particle came from and what it was.

3. The Problem: The "Traffic Jam" Effect

The scientists wanted to make the avalanche as big as possible (a gain of 100,000 to 1,000,000) to see the faintest signals. But they noticed a weird glitch.

When the "snowball" of electrons gets too big and too crowded inside the tiny GEM holes, the detector stops working perfectly. The signal gets weaker than it should be.

• The Analogy: Imagine a narrow hallway (the GEM hole). If one person runs through, they move fast. But if 10,000 people try to run through at the exact same time, they bump into each other, creating a traffic jam.
• The Physics: The electrons are negatively charged, and the "traffic jam" creates a cloud of positive ions (the people left behind). This cloud acts like a shield, blocking the electric field that is supposed to push the electrons forward. The result? The avalanche gets smaller than expected. This is called Gain Saturation.

4. The Discovery: Diffusion is the Hero

The team realized that the distance the electrons travel before hitting the GEMs matters.

• The Scenario: If a particle hits the gas right next to the GEMs, the electrons arrive in a tight, dense clump. This causes a massive traffic jam, and the gain drops significantly.
• The Twist: If the particle hits the gas far away, the electrons drift for a long time. During this drift, they spread out (diffuse), like a drop of ink spreading in a glass of water.
• The Result: When the electrons finally arrive at the GEMs, they are spread out over a wider area. They don't crowd the holes as much. The traffic jam is less severe, and the detector works much better!

5. The Solution: A Mathematical Recipe

The scientists built a computer model (a mathematical recipe) to predict exactly how much the signal would drop based on how crowded the electrons were.

• They tested this with a special source (Iron-55) that shoots X-rays into the gas.
• They moved the source closer and farther away from the detector to change how spread out the electrons were.
• The Outcome: Their model predicted the detector's behavior with 96% accuracy. It successfully explained why the signal gets stronger when the electrons are spread out and weaker when they are bunched up.

Why Does This Matter?

This paper is a crucial step toward building the final, giant detector for the CYGNO experiment.

• For the Future: To find dark matter, the detector needs to be huge and incredibly sensitive. If the scientists didn't understand this "traffic jam" effect, their giant detector might give them wrong answers or miss the ghosts entirely.
• The Takeaway: They have proven that by understanding how electrons spread out and how they crowd together, they can build a detector that sees the invisible universe with crystal-clear precision. They turned a confusing glitch into a predictable tool.

In short: They built a super-camera for gas, realized the electrons get too crowded and block their own signal, figured out that spreading them out fixes the problem, and wrote a rulebook to predict exactly how it all works.

Technical Summary

1. Problem Statement

The CYGNO collaboration is developing a cubic-meter-scale gaseous Time Projection Chamber (TPC) for dark matter (WIMP) and neutrino detection. To achieve the necessary high spatial granularity for reconstructing low-energy (few keV) particle tracks, the project utilizes an optical readout system combined with a Gas Electron Multiplier (GEM) amplification stage.

The Core Issue:
To detect low-energy events, the detector must operate at very high gas gains ($10^5$–$10^6$). However, in this regime, a non-linear dependence of the detector's response on the spatial density of the collected charge has been observed. Specifically, the gain is reduced when the charge density in the GEM holes is high. This "gain saturation" effect, likely caused by space-charge buildup (positive ions shielding the electric field) within the multiplication channels, complicates the accurate reconstruction of energy and position, particularly for events where the ionization cloud is dense or localized.

2. Methodology

Experimental Setup (GIN Prototype):

• Detector: A 2.2-liter prototype named GIN (Gas Ionization Node), featuring a drift volume of 230 mm and a triple-GEM stack (10x10 cm² active area).
• Gas Mixture: Helium/CF$_4$ (60/40), chosen for its suitability in detecting low-mass WIMPs (due to Helium's light nucleus) and its strong electroluminescence in the visible spectrum (CF$_4$).
• Readout:
  • Optical: A scientific CMOS (sCMOS) camera (ORCA-Fusion) coupled via a lens system to capture photons emitted during the avalanche in the GEMs.
  • Supplementary: Two Photomultiplier Tubes (PMTs) for triggering.
• Source: A $^{55}$Fe source (5.9 keV X-rays) was used to generate mono-energetic electron recoils. The source was moved to various distances ($z$) from the GEM plane (4 cm to 22 cm) to vary the drift distance and, consequently, the transverse diffusion of the electron cloud.

Data Analysis Strategy:

1. Event Reconstruction: An algorithm based on DBSCAN identified clusters of pixels corresponding to photon interactions.
2. Selection: Cuts were applied on cluster integral (light intensity), length, and "slimness" (shape factor) to isolate $^{55}$Fe events from background (muons, alphas, noise).
3. Variable Extraction:
   • Spot Size ($\sigma$): Measured via Gaussian fits to the light profile to quantify electron diffusion.
   • Light Integral ($I_{SC}$): The total collected light, proportional to the gas gain.
4. Modeling: A phenomenological model was developed to describe the gain saturation effect based on space-charge screening.

3. Key Contributions

A. Characterization of Diffusion and Gain Dependence
The study quantified the relationship between drift distance and the transverse diffusion coefficient ($D_T$). It confirmed that as the drift distance increases, the electron cloud spreads (increasing $\sigma$), which reduces the charge density entering the GEM holes.

B. Observation of Gain Saturation
The data revealed a clear trend: Gas gain increases with drift distance.

• At short drift distances (small $\sigma$), the charge density in the GEM holes is high, leading to significant space-charge screening and lower effective gain.
• At long drift distances (large $\sigma$), the charge is more diffuse, reducing the screening effect and allowing the gain to approach its theoretical maximum.

C. Development of a Gain Saturation Model
The authors derived a mathematical model to explain this behavior:

• Mechanism: Positive ions produced in the avalanche migrate slowly, creating a local electric field ($E_p$) that opposes the applied field ($E_{GEM}$), effectively reducing the field accelerating electrons.
• Equation: The model modifies the standard Townsend equation to include a screening parameter ($\beta'$) dependent on the number of primary electrons ($n_{in}$) and the charge density.
  $$G = \frac{\tilde{G}}{1 + \beta' n_{in} (\tilde{G} - 1)}$$
  Where $\tilde{G}$ is the unsaturated gain.
• Integration: The model accounts for the 3D Gaussian distribution of the electron cloud, linking the charge density in a single GEM hole to the transverse diffusion width ($\sigma$).

4. Results

• Diffusion Coefficients: The transverse diffusion coefficient was measured as $D_T = (125 \pm 5) \, \mu\text{m}/\sqrt{\text{cm}}$ at 1.0 kV/cm, consistent with Garfield simulations and previous measurements.
• Model Fit: The proposed gain saturation model was fitted to the experimental data across three different GEM voltages (420 V, 430 V, 440 V) and a range of drift distances.
  • The fit yielded a reduced $\chi^2$ of 0.7 with a p-value of 0.81.
  • The model successfully reproduced the gain behavior over nearly one order of magnitude (from $\sim 2 \times 10^5$ to $7 \times 10^5$).
• Precision: The model predicts the gain with a relative precision of 4% (RMS of residuals), demonstrating that the gain dependence on charge density is well-understood and predictable.
• Parameters: The fit provided physical parameters for the gas mixture, including the ionization coefficient ($\tilde{\alpha}$) and the screening parameter ($p$).

5. Significance

• Validation of Optical Readout: The study confirms that optical readout of GEM-based TPCs is viable for high-granularity detection, provided that non-linear gain effects are modeled.
• Simulation Tool: The derived model serves as a critical tool for the CYGNO collaboration to simulate detector responses for different energy depositions and drift distances. It allows for the correction of gain saturation in data analysis, ensuring accurate energy reconstruction for rare event searches.
• Design Optimization: Understanding the space-charge effects allows for the optimization of future TPC designs (e.g., adjusting GEM geometry or gas mixtures) to minimize gain non-linearity, potentially enabling linear response even at very high ionization densities.
• Dark Matter Sensitivity: By accurately modeling the detector response, the experiment can better distinguish between nuclear recoils (WIMP signals) and electron recoils (background), enhancing the sensitivity of the CYGNO experiment to low-mass dark matter candidates.