pith. sign in

arxiv: 2606.05385 · v2 · pith:K62IOTZ5new · submitted 2026-06-03 · ⚛️ physics.ins-det · hep-ex· nucl-ex

GPU optical photon Monte Carlo for noble liquid detectors: validation against Geant4 in a large liquid argon TPC benchmark

Pith reviewed 2026-06-28 03:13 UTC · model grok-4.3

classification ⚛️ physics.ins-det hep-exnucl-ex
keywords optical photon Monte CarloGPU accelerationliquid argon TPCwavelength shiftingGeant4 validationnoble liquid detectorsphoton transportCUDA OptiX
0
0 comments X

The pith

A GPU optical photon Monte Carlo tool matches Geant4 results to subpercent accuracy in a large liquid argon TPC benchmark while providing over 1000 times speedup.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces Simphony, a GPU-accelerated optical photon simulation built on Opticks with CUDA and OptiX, implementing wavelength shifting. It validates this against Geant4 in a 14.7 kt LAr TPC with two-stage wavelength-shifting shell, showing subpercent agreement in detected photon counts for electrons, muons, and protons, and matching spectra. The tool achieves 80 million photons per second on an RTX 4090, yielding 1053 times speedup in optical transport compared to single-thread Geant4, making full optical Monte Carlo practical for detector design and ML datasets.

Core claim

Simphony implements a GPU version of the Geant4 G4OpWLS wavelength-shifting model and returns Monte Carlo truth for detected hits with low per-photon overhead. We validate Simphony against Geant4 11.3.2 in a simplified 14.7 kt liquid argon Time Projection Chamber benchmark with a two-stage wavelength-shifting shell and idealized photon counting detector. For three paired 2.5 GeV electron simulations, each producing about 61 M optical photons, the integrated detected-photon ratio agrees with Geant at the subpercent level. The detected arrival time and wavelength spectra give χ²/ndf values of 0.98 and 1.08. Contained muon and near-Cerenkov-threshold proton samples give R_N=1.0017±0.0008 and R_

What carries the argument

GPU optical photon propagation using OptiX ray tracing combined with CUDA kernels for wavelength shifting and hit detection, allowing stacked event processing for high throughput.

If this is right

  • Detector development studies can incorporate full optical photon Monte Carlo without excessive computational cost.
  • Labeled optical response datasets for machine learning can be generated efficiently at scale.
  • Reconstruction algorithms and particle identification can benefit from detailed optical truth information.
  • The method supports repeated geometry-dependent simulations needed for optimization.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach could be extended to other noble liquid or scintillator detectors with similar optical processes.
  • Hybrid workflows combining GPU optics with CPU-based particle transport might reduce overall simulation times further.
  • Real-time optical response feedback during detector design iterations becomes feasible with this speed.

Load-bearing premise

The validation uses a simplified benchmark geometry with idealized photon counting and two-stage wavelength shifting, which may not capture complexities of real detector surfaces or material variations.

What would settle it

Running the same simulations in a geometry that includes measured surface reflectivities, material inhomogeneities, or additional optical processes and finding detected photon ratios deviating by more than a few percent would challenge the validation.

read the original abstract

Optical photon Monte Carlo simulation is a computational bottleneck for noble liquid Time Projection Chambers. Design studies require repeated, geometry dependent simulations of timing, wavelength shifting, and optical response, while reconstruction and particle identification workflows need labeled optical datasets. We present Simphony, a GPU optical simulation tool, formerly EIC-Opticks, built on Opticks with CUDA and NVIDIA OptiX. Simphony implements a GPU version of the Geant4 G4OpWLS wavelength-shifting model and returns Monte Carlo truth for detected hits with low per-photon overhead. We validate Simphony against Geant4 11.3.2 in a simplified 14.7 kt liquid argon Time Projection Chamber benchmark with a two-stage wavelength-shifting shell and idealized photon counting detector. For three paired 2.5 GeV electron simulations, each producing about 61 M optical photons, the integrated detected-photon ratio agrees with Geant at the subpercent level. The detected arrival time and wavelength spectra give $\chi^2/\mathrm{ndf}$ values of 0.98 and 1.08. Contained muon and near-Cerenkov-threshold proton samples give $R_N=1.0017\pm0.0008$ and $R_N=1.0005\pm0.0014$, confirming agreement for distinct source topologies. On an NVIDIA RTX 4090, a stacked launch of four 2.5 GeV electron events transports 243 M optical photons in $3.03\pm0.06$ s, giving $80.2\pm1.6$ M photons s$^{-1}$. Relative to a single-thread Geant reference and including GPU overheads and host-device transfers, the optical transport speedup is $1053\pm55$; the end-to-end wall time acceleration is $89\pm5$. These results show that Simphony can make explicit optical photon Monte Carlo practical for detector development studies and for generating machine learning optical response datasets.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The paper introduces Simphony (formerly EIC-Opticks), a GPU-based optical photon Monte Carlo tool built on Opticks with CUDA and NVIDIA OptiX that implements the Geant4 G4OpWLS wavelength-shifting model. It validates the implementation against Geant4 11.3.2 in a simplified 14.7 kt liquid argon TPC benchmark geometry with a two-stage wavelength-shifting shell and idealized photon-counting detector. For paired 2.5 GeV electron simulations (~61 M photons each), the detected-photon ratio agrees at the sub-percent level; detected arrival-time and wavelength spectra yield χ²/ndf of 0.98 and 1.08. Contained-muon and near-Cerenkov-threshold proton samples give R_N = 1.0017 ± 0.0008 and 1.0005 ± 0.0014. On an RTX 4090, stacked launches of four electron events transport 243 M photons in 3.03 ± 0.06 s (80.2 ± 1.6 M photons s⁻¹), yielding optical-transport speedup 1053 ± 55 and end-to-end wall-time acceleration 89 ± 5 relative to single-thread Geant4.

Significance. If the reported agreement holds, the work removes a major computational bottleneck for noble-liquid TPC design and reconstruction studies by making explicit optical-photon Monte Carlo practical at scale. The multi-metric, multi-topology validation (detected-photon ratios, χ² on spectra, R_N for electrons/muons/protons) with explicit uncertainties and direct comparison to an independent external code (Geant4 11.3.2) provides strong evidence that the GPU G4OpWLS implementation reproduces the reference for the stated benchmark. The performance numbers, including host-device overheads, are directly usable for planning large-scale simulations and ML training-set generation.

minor comments (2)
  1. §3 (benchmark geometry): the description of the two-stage wavelength-shifting shell and idealized photon-counting detector would benefit from an explicit statement of the optical surface properties and absorption lengths used, to facilitate exact reproduction.
  2. Figure 4 (time and wavelength spectra): the binning and normalization conventions for the χ² calculation are not stated; adding this detail would allow readers to verify the reported χ²/ndf values directly.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of the manuscript, detailed summary of the validation results, and recommendation to accept.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper implements a GPU port of Geant4's G4OpWLS model and validates it via direct numerical comparison of photon counts, arrival-time spectra, wavelength spectra, and topology-specific ratios against independent runs of Geant4 11.3.2. No parameters are fitted inside the paper, no predictions are derived from internal fits, and no load-bearing steps reduce to self-citations or self-definitions; all central claims rest on external code-to-code agreement in a fixed benchmark geometry.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

This is a software implementation and validation study against an established reference; no new physical parameters, axioms, or invented entities are introduced.

pith-pipeline@v0.9.1-grok · 5916 in / 1326 out tokens · 62034 ms · 2026-06-28T03:13:45.051611+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

28 extracted references · 25 canonical work pages · 3 internal anchors

  1. [1]

    Agostinelliet al.,Geant4–a simulation toolkit, Nucl

    S. Agostinelliet al.,Geant4–a simulation toolkit, Nucl. Instrum. Meth. A506(2003) 250

  2. [2]

    Allisonet al.,Recent developments in Geant4, Nucl

    J. Allisonet al.,Recent developments in Geant4, Nucl. Instrum. Meth. A835(2016) 186

  3. [3]

    Foreman, R

    W. Foreman, R. Acciarri, J. A. Asaadi, W. Badgett, F. d. M. Blaszczyk, R. Bouabid, C. Bromberg, R. Carey, F. Cavanna, J. I. Cevallos Aleman, A. Chatterjee, J. Evans, A. Falcone, W. Flanagan, B. T. Fleming, D. Garcia-Gomez, B. Gelli, T. Ghosh, R. A. Gomes, E. Gramellini, R. Gran, P. Hamilton, C. Hill, J. Ho, J. Hugon, E. Iwai, E. Kearns, E. Kemp, T. Kobila...

  4. [4]

    X. Ning, W. Shi, C. Zhang, C. Riccio and J. H. Jo,Self-compensating Light Calorimetry with Liquid Argon Time Projection Chamber for GeV Neutrino Physics, Phys. Rev. D111(2025) 032007, arXiv:2410.04603, doi:10.1103/PhysRevD.111.032007. – 22 –

  5. [5]

    Capozzi, S

    F. Capozzi, S. W. Li, G. Zhu and J. F. Beacom,DUNE as the Next-Generation Solar Neutrino Experiment, Phys. Rev. Lett.123(2019) 131803, arXiv:1808.08232, doi:10.1103/PhysRevLett.123.131803

  6. [6]

    W. Shi, X. Ning, D. Pershey, F. Marinho, C. Riccio, J. H. Jo, C. Zhang and F. Cavanna,Physics Prospects with MeV Neutrino-Argon Charged Current Interactions using Enhanced Photon Detection in Future LArTPCs, Phys. Rev. D112(2025) 012019, arXiv:2502.18498, doi:10.1103/f3t7-5vmd

  7. [7]

    G. Zhu, S. W. Li and J. F. Beacom,Developing the MeV potential of DUNE: Detailed considerations of muon-induced spallation and other backgrounds, Phys. Rev. C99(2019) 055810, arXiv:1811.07912, doi:10.1103/PhysRevC.99.055810

  8. [8]

    Abiet al.,Supernova neutrino burst detection with the Deep Underground Neutrino Experiment, Eur

    DUNE Collaboration, B. Abiet al.,Supernova neutrino burst detection with the Deep Underground Neutrino Experiment, Eur. Phys. J. C81(2021) 423, arXiv:2008.06647, doi:10.1140/epjc/s10052-021-09166-w

  9. [9]

    Abed Abudet al.,Supernova Pointing Capabilities of DUNE, Phys

    DUNE Collaboration, A. Abed Abudet al.,Supernova Pointing Capabilities of DUNE, Phys. Rev. D 111(2025) 092006, arXiv:2407.10339, doi:10.1103/PhysRevD.111.092006

  10. [10]

    Measuring the supernova unknowns at the next-generation neutrino telescopes through the diffuse neutrino background

    K. Møller, A. M. Suliga, I. Tamborra and P. B. Denton,Measuring the supernova unknowns at the next-generation neutrino telescopes through the diffuse neutrino background, JCAP05(2018) 066, arXiv:1804.03157, doi:10.1088/1475-7516/2018/05/066

  11. [11]

    Abed Abudet al.,DUNE Phase II: Scientific Opportunities, Detector Concepts, Technological Solutions, JINST19(2024) P12005, arXiv:2408.12725, doi:10.1088/1748-0221/19/12/P12005

    DUNE Collaboration, A. Abed Abudet al.,DUNE Phase II: Scientific Opportunities, Detector Concepts, Technological Solutions, JINST19(2024) P12005, arXiv:2408.12725, doi:10.1088/1748-0221/19/12/P12005

  12. [12]

    Marinho,APEX: Optimized vertical drift PDS for DUNE FD3, arXiv:2503.06032 (2025)

    F. Marinho,APEX: Optimized vertical drift PDS for DUNE FD3, arXiv:2503.06032 (2025)

  13. [13]

    SBND Collaboration,Scintillation light in SBND: simulation, reconstruction, and expected performance of the photon detection system, Eur. Phys. J. C84(2024) 1046, arXiv:2406.07514, doi:10.1140/epjc/s10052-024-13306-3

  14. [14]

    Garcia-Gamez, P

    D. Garcia-Gamez, P. Green and A. M. Szelc,Predicting transport effects of scintillation light signals in large-scale liquid argon detectors, Eur. Phys. J. C81(2021) 349, arXiv:2010.00324, doi:10.1140/epjc/s10052-021-09119-3

  15. [15]

    W. Mu, A. I. Himmel and B. Ramson,Photon detection probability prediction using one-dimensional generative neural network, Mach. Learn.: Sci. Technol.3(2022) 015033, arXiv:2109.07277, doi:10.1088/2632-2153/ac4a2c

  16. [16]

    E. L. Snider and G. Petrillo,LArSoft: toolkit for simulation, reconstruction and analysis of liquid argon TPC neutrino detectors, J. Phys.: Conf. Ser.898(2017) 042057, doi:10.1088/1742-6596/898/4/042057

  17. [17]

    LArSoft Collaboration,PDFastSimPAR_module.cc: semi-analytical fast optical simulation module, LArSoft source documentation,https: //code-doc.larsoft.org/docs/latest/html/PDFastSimPAR__module_8cc_source.html, accessed 2026

  18. [18]

    M. Lei, K. V. Tsang, S. Gasiorowski, C. Li, Y. Nashed, G. Petrillo, O. Piazza, D. Ratner and K. Terao, Implicit neural representation as a differentiable surrogate for photon propagation in a monolithic neutrino detector, arXiv:2211.01505

  19. [19]

    Blyth,Opticks: GPU optical photon simulation for particle physics with NVIDIA OptiX, EPJ Web Conf.214(2019) 02027, doi:10.1051/epjconf/201921402027

    S. Blyth,Opticks: GPU optical photon simulation for particle physics with NVIDIA OptiX, EPJ Web Conf.214(2019) 02027, doi:10.1051/epjconf/201921402027

  20. [20]

    S. C. Blyth,Opticks: GPU optical photon simulation via NVIDIA OptiX, EPJ Web Conf.295(2024) 11014, doi:10.1051/epjconf/202429511014. – 23 –

  21. [21]

    Wenzel, S.-Y

    H. Wenzel, S.-Y. Jun, K. Genser and F. De Figueiredo,CaTS: Integration of Geant4 and Opticks, EPJ Web Conf.295(2024) 11004, doi:10.1051/epjconf/202429511004

  22. [22]

    Y. Li, A. Davis, S. Easo, K. Evans, E. M. Gersabeck, M. Gersabeck, L. G. Girardey, R. Nandakumar and A. J. Pointer,GPU-based optical photon simulation for the LHCb RICH 1 detector, Eur. Phys. J. C 83(2023) 1036, arXiv:2307.10823, doi:10.1140/epjc/s10052-023-12158-7

  23. [23]

    Galgoczi, K

    G. Galgoczi, K. Kauder, M. Potekhin, S. Rahman, D. Smirnov and T. Wenaus,GPU acceleration of optical photon propagation in low photon yield applications: Opticks for the Electron Ion Collider, arXiv:2512.06061

  24. [24]

    Szydagis, N

    M. Szydagis, N. Barry, K. Kazkaz, J. Mock, D. Stolp, M. Sweany, M. Tripathi, S. Uvarov, N. Walsh and M. Woods,NEST: A comprehensive model for scintillation yield in liquid xenon, JINST6(2011) P10002, doi:10.1088/1748-0221/6/10/P10002

  25. [25]

    M. Szydagis et al.,A review of basic energy reconstruction techniques in liquid xenon and argon detectors for dark matter and neutrino physics using NEST, Instruments5(2021) 13, doi:10.3390/instruments5010013

  26. [26]

    NEXT Collaboration,Performance of an Optical TPC Geant4 Simulation with Opticks GPU-Accelerated Photon Propagation, Eur. Phys. J. C85(2025) 910, arXiv:2502.13215, doi:10.1140/epjc/s10052-025-14612-0

  27. [27]

    Babiczet al.,A measurement of the group velocity of scintillation light in liquid argon, JINST15 (2020) P09009, arXiv:2002.09346

    M. Babiczet al.,A measurement of the group velocity of scintillation light in liquid argon, JINST15 (2020) P09009, arXiv:2002.09346

  28. [28]

    The scintillation of liquid argon

    T. Heindl, T. Dandl, M. Hofmann, R. Krücken, L. Oberauer, W. Potzel, J. Wieser and A. Ulrich,The scintillation of liquid argon, EPL91(2010) 62002, arXiv:1511.07718. A One-dimensional projections for muon and proton primaries The per-bin Simphony/Geant4 ratio analysis of figures 6 for the electron is repeated here for the muon and proton primaries of table...