Simulation of semiconductor detectors in 3D with SolidStateDetectors.jl

This paper introduces SolidStateDetectors.jl, an open-source Julia package designed for efficient, parallel 3D simulation of electric fields, charge carrier drifts, and signal pulses in semiconductor detectors, with a specific focus on germanium devices.

The Gist

Imagine you are trying to listen to a single, tiny whisper in a massive, echoing cathedral. That's essentially what scientists do when they use Germanium detectors to hunt for rare cosmic events, like dark matter or the mysterious "neutrinoless double-beta decay." These detectors are incredibly sensitive, but to understand what they are hearing, you need to know exactly how the sound (or in this case, electric signals) travels through the building.

This paper introduces a new, open-source software tool called SolidStateDetectors.jl (let's call it "SSD" for short) that acts like a super-powered virtual simulator for these detectors.

Here is a breakdown of what the paper is about, using everyday analogies:

1. The Problem: The "Black Box" Detector

High-purity Germanium detectors are like giant, frozen blocks of pure crystal. When a particle hits them, it creates a tiny spark of electricity (electrons and holes) that drifts through the crystal to be read by sensors.

• The Challenge: The path these sparks take isn't a straight line. It's like trying to predict how a leaf will drift down a river with hidden currents, rocks, and varying water depths.
• The Old Way: Previous software was either closed off (like a locked black box you couldn't change), couldn't handle complex 3D shapes, or ignored the environment surrounding the detector (like the metal case or liquid argon it sits in).

2. The Solution: The "Digital Twin"

The authors built SSD to create a perfect digital twin of a real detector.

• The Blueprint: You can tell the software exactly what the detector looks like, what it's made of, and what's around it (like a copper shell or liquid argon). It uses a method called "Constructive Solid Geometry," which is like building with digital LEGO blocks to create complex shapes.
• The Map: Before simulating any events, the software calculates the Electric Field. Think of this as mapping the invisible "wind" inside the detector that pushes the electrons. It also calculates "Weighting Potentials," which are like a map showing how much a specific sensor will "hear" a spark if it passes by.

3. How It Works: The Race

Once the map is ready, the software simulates what happens when a particle hits the detector:

1. The Start: A particle hits the crystal, creating a cloud of electrons and "holes" (missing electrons).
2. The Drift: The software tracks these clouds as they race toward the electrodes. It calculates their speed step-by-step, just like a GPS tracking a car's speed based on traffic and road conditions.
3. The Signal: As the clouds move, they induce a signal (a pulse) on the sensors. The software sums up these signals to create a "pulse shape"—the unique waveform that tells scientists where the hit happened and what kind of particle it was.

4. Why It's Special: The "Julia" Engine

The software is written in a programming language called Julia.

• The Analogy: Imagine a race car engine. Most simulation tools are like a reliable sedan—slow but steady. SSD is built with a Formula 1 engine. It is designed to be incredibly fast and can run many simulations at the same time (parallel processing). This allows scientists to test thousands of scenarios quickly, which is crucial for designing the next generation of massive detectors.

5. Testing the Engine: The "Sound Check"

The authors didn't just build the software; they tested it against real data from a real detector (a segmented Germanium crystal).

• The Test: They compared the "virtual pulses" generated by the software with the "real pulses" recorded in the lab.
• The Result: The match was excellent. The software could predict the shape of the electrical signal so accurately that it could distinguish between:
  • Surface events: Particles hitting the outside skin of the detector (often background noise).
  • Bulk events: Particles hitting the deep interior (the signal they want).
• The "A/E" Trick: They used a specific metric called the A/E parameter (Amplitude divided by Energy). Think of this as a "voice pitch test." Single-site events (good signals) have a different pitch than multi-site events (background noise). The simulation predicted this pitch perfectly, proving it can help scientists filter out the noise in future experiments.

6. Why This Matters

This tool is a game-changer for the LEGEND experiment and others searching for the secrets of the universe.

• Optimization: Before cutting a massive, expensive crystal into a detector, scientists can use SSD to simulate different cuts. They can ask, "If I cut it this way, will the electric field get stuck in a corner?" This saves money and ensures the final detector works perfectly.
• Background Rejection: By accurately simulating how background noise looks, scientists can teach their AI algorithms to ignore it, making their search for rare events much more sensitive.

Summary

In short, SolidStateDetectors.jl is a fast, flexible, and open-source "flight simulator" for particle detectors. It allows scientists to build, test, and perfect their detectors in a virtual world before they ever touch a real crystal, ensuring that when they finally listen for the whispers of the universe, they know exactly what they are hearing.

Technical Summary

1. Problem Statement

High-purity germanium (HPGe) detectors are critical for nuclear and particle physics experiments, particularly in neutrinoless double-beta decay ($0\nu\beta\beta$) searches (e.g., LEGEND, Gerda, Majorana) and gamma-ray tracking (e.g., AGATA, GRETINA). These experiments rely heavily on pulse-shape analysis (PSA) to reconstruct interaction positions and reject background events (such as multi-site interactions).

Existing simulation tools for HPGe detectors suffer from several limitations:

• They often neglect the detector environment (surrounding materials, cryostats, liquid argon immersion), which significantly alters electric fields near passivated surfaces.
• Many are not open-source or are difficult to extend for custom geometries and impurity profiles.
• They lack the flexibility to handle complex 3D geometries efficiently while utilizing detector symmetries.
• There is a need for a tool that can optimize crystal cutting and detector design by simulating active impurity profiles and undepleted regions.

2. Methodology

The authors introduce SolidStateDetectors.jl, an open-source software package written in the Julia programming language, designed to simulate electric fields, charge carrier drift, and induced signals in 3D.

Key Technical Components:

• Geometry Definition: Uses Constructive Solid Geometry (CSG) defined in structured text configuration files. It supports both Cartesian and cylindrical coordinate systems and allows for the inclusion of the detector's surroundings (e.g., copper shields, liquid argon).
• Electric Field Calculation:
  • Solves Gauss's law ($\nabla(\epsilon_0 \epsilon_r(r) \nabla \Phi(r)) = -\rho(r)$) numerically.
  • Uses a Successive Over-Relaxation (SOR) algorithm on an adaptive grid with red-black division for parallelization.
  • Supports adaptive grid refinement: The grid starts coarse and refines iteratively based on potential gradients, ensuring high precision where needed without excessive computational cost.
  • Handles undepleted regions by setting impurity density to zero in specific volumes, allowing the simulation of depletion zone development.
• Weighting Potentials: Calculated for each electrode using the same SOR algorithm to determine induced signals via the Shockley-Ramo theorem.
• Charge Drift & Signal Induction:
  • Simulates the drift of electron and hole clouds using user-defined drift velocity models (defaulting to AGATA Detector Library models).
  • Calculates drift paths step-by-step (default 1 ns) based on interpolated electric fields.
  • Computes induced signals on electrodes using weighting potentials.
  • Includes models for charge trapping (in bulk or near surfaces) and allows for virtual volumes to simulate surface effects.
• Integration: Designed to work with Geant4 for radiation interaction simulation. Geant4 provides energy deposition "hits," which are then clustered and drifted within SolidStateDetectors.jl.

3. Key Contributions

• Full 3D Simulation with Symmetry: Capable of full 3D calculations for segmented detectors while efficiently utilizing symmetries (e.g., periodic boundary conditions in $\phi$) for simpler geometries.
• Environment Modeling: Uniquely capable of simulating the impact of surrounding materials (vacuum, copper shells, liquid argon) on electric fields, which is crucial for surface event background rejection.
• Julia Implementation: Leverages Julia's high-level syntax and high-performance capabilities (just-in-time compilation, parallelism) to execute complex 3D field calculations and drift simulations efficiently.
• Modularity and Extensibility: The package is open-source and modular, allowing users to easily adapt it for different semiconductor types (e.g., silicon) or custom impurity profiles.
• Validation Framework: Includes automated tests comparing numerical results against analytical solutions for standard geometries (coaxial, parallel plate).

4. Results and Verification

The package was validated using data from an n-type segmented BEGe (Broad Energy Germanium) detector with a point-contact core and four-fold segmentation.

• Field Calculation Verification:
  • Compared numerical results against analytical solutions for an infinite true-coaxial detector.
  • Achieved extremely low RMS differences ($2.1 \times 10^{-6}$ V for cylindrical coordinates), confirming the accuracy of the SOR solver and adaptive grid.
  • Capacitance calculations matched analytical values within 0.03% to 1%.
• Environment Impact: Simulations showed that surrounding materials (e.g., liquid argon vs. vacuum) significantly alter the electric potential near passivated surfaces, a critical factor for background modeling.
• Pulse Shape Comparison:
  • Surface Events: Simulated "super-pulses" (averaged waveforms) for 81 keV photons matched experimental data well, validating electron drift models.
  • Bulk Events: Simulated pulses for bulk interactions (using Compton-scattered 662 keV photons) agreed with data. The inclusion of differential cross-talk modeling significantly improved the agreement for mirror pulses in non-collecting segments.
• Pulse-Shape Analysis (PSA):
  • The simulation successfully reproduced the $A/E$ parameter distributions (ratio of current amplitude to energy) for Double Escape (DEP), Single Escape (SEP), and Full Energy Peaks (FEP).
  • Background Rejection: When applying an $A/E$ cut to achieve 90% survival for single-site DEP events, the simulation predicted a multi-site rejection rate of 10.5%, closely matching the experimental data of 12.3%.

5. Significance

• Experimental Optimization: The tool is essential for the LEGEND experiment, allowing for the optimization of crystal cutting to maximize detector mass while maintaining reasonable bias voltages and ensuring full depletion.
• Background Reduction: By accurately simulating surface effects and environmental interactions, the package aids in developing better background rejection algorithms for rare-event searches.
• Community Resource: As an open-source, extensible Julia package, it lowers the barrier to entry for high-fidelity detector simulation, replacing proprietary or rigid legacy codes.
• Future Roadmap: The paper outlines upcoming features in version 1.0, including charge cloud self-repulsion/diffusion, stochastic surface trapping models, and GPU acceleration, positioning the software as a long-term standard for semiconductor detector simulation.