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
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.
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
- 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.
Referee Report
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)
- §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.
- 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
We thank the referee for their positive assessment of the manuscript, detailed summary of the validation results, and recommendation to accept.
Circularity Check
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
Reference graph
Works this paper leans on
-
[1]
Agostinelliet al.,Geant4–a simulation toolkit, Nucl
S. Agostinelliet al.,Geant4–a simulation toolkit, Nucl. Instrum. Meth. A506(2003) 250
2003
-
[2]
Allisonet al.,Recent developments in Geant4, Nucl
J. Allisonet al.,Recent developments in Geant4, Nucl. Instrum. Meth. A835(2016) 186
2016
-
[3]
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]
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]
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]
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]
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
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevc.99.055810 2019
-
[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]
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]
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
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/1475-7516/2018/05/066 2018
-
[11]
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]
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]
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]
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]
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]
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]
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
2026
- [18]
-
[19]
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]
S. C. Blyth,Opticks: GPU optical photon simulation via NVIDIA OptiX, EPJ Web Conf.295(2024) 11014, doi:10.1051/epjconf/202429511014. – 23 –
-
[21]
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]
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]
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]
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]
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]
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]
M. Babiczet al.,A measurement of the group velocity of scintillation light in liquid argon, JINST15 (2020) P09009, arXiv:2002.09346
-
[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...
work page internal anchor Pith review Pith/arXiv arXiv 2010
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.