arxiv: 2605.13987 · v1 · submitted 2026-05-13 · ⚛️ physics.ins-det

Recognition: 2 theorem links

· Lean Theorem

GPU-Accelerated Simulation of 3D Diamond Sensors Using the TeRABIT Infrastructure

Lucio Anderlini , Clarissa Buti , Elia Eredi , Giovanni Passaleva

Authors on Pith no claims yet

Pith reviewed 2026-05-15 02:33 UTC · model grok-4.3

classification ⚛️ physics.ins-det

keywords 3D diamond sensorstime-dependent weighting potentialsGPU accelerationRamo-Shockley theoremsignal simulationlaser graphitizationradiation detectorsTeRABIT infrastructure

0 comments

The pith

GPU porting of a third-order PDE solver for time-dependent weighting potentials speeds up 3D diamond sensor simulations.

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

This paper develops a simulation method for diamond detectors with laser-graphitized 3D electrodes by extending the Ramo-Shockley theorem to include signal propagation via time-dependent weighting potentials. These potentials come from solving a third-order partial differential equation obtained as a quasi-static approximation to Maxwell's laws. The numerical solution of that equation, previously the main bottleneck, is handled with fundamental solutions for boundary conditions and spectral methods for the detector bulk. The solver has been moved to GPUs and distributed over multiple sites through the TeRABIT infrastructure, cutting computation time enough to support what-if explorations of electrode geometry. The work targets applications in high-energy physics, nuclear medicine, and dosimetry where time resolutions better than 100 ps are now achievable in fabrication.

Core claim

The central claim is that an innovative numerical solver for time-dependent weighting potentials—obtained by solving a third-order PDE derived from the quasi-static approximation of Maxwell's laws using fundamental solutions to set boundary conditions and spectral methods across the detector volume—has been ported to GPUs and distributed across computing sites with the TeRABIT HPC infrastructure. This change removes the previous computational barrier and makes practical the full modeling of intertwined effects from energy deposition fluctuations, carrier transport, signal propagation, and readout electronics inside full-carbon 3D diamond sensors.

What carries the argument

The GPU-accelerated solver for the third-order PDE of time-dependent weighting potentials, using fundamental solutions to impose boundary conditions and spectral methods to extend the solution through the detector bulk.

If this is right

Signal formation simulations can now treat propagation effects inside the sensor in a theoretically consistent way.
Computation time reductions make systematic studies of how electrode geometry influences timing performance feasible.
The framework supports disentangling contributions from energy deposition, transport, propagation, and electronics in a single model.
Applications in high-energy physics, nuclear medicine, and dosimetry gain access to faster iteration on sensor designs.

Where Pith is reading between the lines

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

The same GPU solver could be applied to other wide-bandgap materials with 3D electrode structures to test generality beyond diamond.
Integration with Monte Carlo codes for primary energy deposition would create end-to-end sensor design loops.
If the model matches data, it could guide fabrication toward electrode layouts that push time resolution below current 100 ps records.

Load-bearing premise

The quasi-static approximation of Maxwell's laws remains accurate enough to model signal propagation and readout effects inside the diamond volume without a full time-dependent electromagnetic treatment.

What would settle it

A direct comparison of the model's predicted time resolution and pulse shapes against experimental measurements on fabricated diamond sensors with known electrode geometries and varying graphitization parameters.

read the original abstract

Diamond detectors with electrodes orthogonal to the surface, engraved via laser-induced graphitization, are full-carbon sensors of interest for a wide range of applications, spanning from High Energy Physics to Nuclear Medicine and dosimetry. In recent years, significant progress has been made in graphitization techniques, enabling the fabrication of lower-resistance electrodes. This has resulted in faster sensors, achieving time resolutions better than 100 ps. However, simulating signal formation in these devices remains a challenge. The effects of fluctuations in energy deposition, carrier transport, signal propagation, and readout electronics intertwine in a way that is non-trivial to disentangle. We have developed an innovative simulation approach based on an extension of the Ramo-Shockley theorem, modeling propagation effects in a theoretically sound manner by introducing time-dependent weighting potentials. These are obtained by solving a third-order partial differential equation derived as a quasi-static approximation of Maxwell's laws. The numerical solution of this equation emerged as the main challenge of the new approach. In this contribution, we discuss an innovative solver that uses fundamental solutions to impose boundary conditions and spectral methods to extend the solution to the bulk of the diamond detector. We report on how the solver has recently been ported to GPUs and distributed across multiple computing sites, leveraging the TeRABIT HPC Bubbles and the InterLink protocol. This drastically reduces time-to-insight and effectively enables what-if studies on sensor geometry.

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.

Desk Editor's Note private letter to a colleague

The paper's real advance is the GPU-distributed spectral solver for time-dependent weighting potentials in 3D diamond sensors, but it lacks any validation data or checks on the quasi-static approximation.

read the letter

The key takeaway is that this work ports a spectral solver for time-dependent weighting potentials in 3D diamond sensors to GPUs using the TeRABIT infrastructure, making geometry optimization studies practical. They extend the Ramo-Shockley theorem with a third-order PDE derived from quasi-static Maxwell equations and solve it by combining fundamental solutions at boundaries with spectral methods in the volume. This is useful because diamond sensors with laser-graphitized electrodes are getting fast enough that standard simulations struggle with the intertwined effects of charge transport and signal propagation. The GPU distribution across sites reduces time-to-insight, which is a concrete engineering win for the field. The approach looks technically coherent. It avoids fitting parameters and builds directly on established electrostatics with a controlled extension for time dependence. The citation pattern seems fine since it references the base theorem without circularity. Where it falls short is validation. The description stops at the method and the speed-up; there are no results shown against real data or against full electromagnetic simulations. For signals under 100 ps, the quasi-static approximation might miss important effects like finite propagation speed, and without an error estimate or benchmark the reliability for design work remains open. That part needs more evidence before the claims about enabling reliable what-if studies can be taken at face value. This paper is for instrumentation physicists who model diamond detectors. A reader interested in numerical methods for sensor simulation would get value from the solver details and the HPC implementation. It deserves peer review because the technical contribution is specific and the problem it addresses is current, even though revisions will likely be needed to add validation data.

Referee Report

2 major / 0 minor

Summary. The paper describes the development of a numerical solver for time-dependent weighting potentials in 3D diamond sensors, obtained by solving a third-order PDE that extends the Ramo-Shockley theorem under a quasi-static approximation to Maxwell's equations. The solver employs fundamental solutions for boundary conditions and spectral methods for the bulk, and has been ported to GPUs and distributed via the TeRABIT infrastructure to reduce computation time and enable what-if studies on sensor geometry for devices targeting sub-100 ps time resolution.

Significance. If the numerical approach and quasi-static approximation prove accurate, the work would provide a practical tool for optimizing fast diamond detectors used in high-energy physics and dosimetry, where disentangling energy deposition, transport, and propagation effects is otherwise intractable. The GPU acceleration and multi-site distribution represent a concrete engineering advance that could lower barriers to iterative sensor design.

major comments (2)

[Abstract] Abstract: The central claim that the solver 'enables reliable what-if studies' lacks any supporting validation data, error analysis, comparison to measurements, or benchmarks against full-wave electromagnetic solutions. Without these, the reliability for <100 ps signals cannot be assessed.
[Abstract] Abstract and method description: The quasi-static approximation yielding the third-order PDE neglects retardation and radiation terms; no validity criterion (e.g., sensor size versus wavelength at relevant frequencies) or error estimate is supplied to confirm its accuracy inside the diamond volume for the targeted time resolutions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback. Our manuscript focuses on the numerical implementation of the spectral solver for time-dependent weighting potentials and its GPU acceleration via the TeRABIT infrastructure. We address the two major comments below and will revise the manuscript to strengthen the presentation of the method's scope and limitations.

read point-by-point responses

Referee: [Abstract] Abstract: The central claim that the solver 'enables reliable what-if studies' lacks any supporting validation data, error analysis, comparison to measurements, or benchmarks against full-wave electromagnetic solutions. Without these, the reliability for <100 ps signals cannot be assessed.

Authors: We agree that the abstract's phrasing overstates the current evidence for reliability. The manuscript demonstrates substantial reductions in computation time through GPU porting and multi-site distribution, which in practice allow iterative geometry studies that were previously intractable. However, direct validation data, error analysis, or comparisons to measurements and full-wave solutions are not included, as the contribution centers on the solver architecture and performance. In the revised manuscript we will moderate the abstract claim to 'facilitates what-if studies' and add a new subsection on numerical convergence of the spectral method together with error estimates for representative 3D geometries. revision: yes
Referee: [Abstract] Abstract and method description: The quasi-static approximation yielding the third-order PDE neglects retardation and radiation terms; no validity criterion (e.g., sensor size versus wavelength at relevant frequencies) or error estimate is supplied to confirm its accuracy inside the diamond volume for the targeted time resolutions.

Authors: The quasi-static approximation is adopted because typical 3D diamond sensor dimensions (a few mm) are much smaller than the electromagnetic wavelength at frequencies corresponding to 100 ps signals (wavelengths of order 10 cm in diamond). We acknowledge that the manuscript provides neither an explicit validity criterion nor quantitative error estimates for the neglected terms. In the revision we will insert a short derivation of the validity condition (sensor size / wavelength ≪ 1) and supply order-of-magnitude error bounds obtained by comparing the quasi-static solution against the leading retardation correction for the relevant frequency range. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation chain is self-contained numerical extension of established theorem

full rationale

The paper starts from the established Ramo-Shockley theorem and introduces time-dependent weighting potentials via a third-order PDE derived as a quasi-static approximation to Maxwell's equations. The core contribution is the numerical solver (fundamental solutions for boundary conditions plus spectral methods for the bulk), which is then ported to GPUs and distributed via TeRABIT. No step reduces the final result to a fitted parameter, self-referential definition, or self-citation chain; the GPU acceleration and what-if studies are direct computational outcomes, not predictions forced by construction. The approach remains independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the quasi-static Maxwell approximation to derive the third-order PDE and on the numerical stability of the fundamental-solution plus spectral-method solver; no free parameters or invented entities are introduced in the abstract.

axioms (1)

domain assumption Quasi-static approximation of Maxwell's laws accurately describes signal propagation in diamond detectors
Invoked to obtain the third-order PDE for time-dependent weighting potentials

pith-pipeline@v0.9.0 · 5562 in / 1268 out tokens · 28419 ms · 2026-05-15T02:33:01.241201+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

time-dependent weighting potentials... obtained by solving a third-order partial differential equation derived as a quasi-static approximation of Maxwell’s laws
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

extension of the Ramo-Shockley theorem... correlation function H(x,t)

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

15 extracted references · 15 canonical work pages · 1 internal anchor

[1]

Michal Pomorski , title = "

work page
[2]

High Mobility Diamonds and Particle Detectors

Pernegger, Heinz , publisher =. High Mobility Diamonds and Particle Detectors. Physics and Applications of CVD Diamond , chapter =. doi:https://doi.org/10.1002/9783527623174.ch12 , url =. https://onlinelibrary.wiley.com/doi/pdf/10.1002/9783527623174.ch12 , year =

work page doi:10.1002/9783527623174.ch12
[3]

Operation and performance of the active target of PADME

Oliva, F. Operation and performance of the active target of PADME. Nucl. Instrum. Meth. A. 2020. doi:10.1016/j.nima.2019.162354

work page doi:10.1016/j.nima.2019.162354 2020
[4]

Performance of 3D diamond detectors in small field dosimetry: the impact of pixel size

Keida Kanxheri and Cinzia Talamonti and Silvio Sciortino and Stefano Lagomarsino and Maria Ionica and Mirco Caprai and Francesco Moscatelli and Leonello Servoli , keywords =. Performance of 3D diamond detectors in small field dosimetry: the impact of pixel size. Phys. Med. , volume =. 2022 , issn =. doi:https://doi.org/10.1016/j.ejmp.2022.09.006 , url =

work page doi:10.1016/j.ejmp.2022.09.006 2022
[5]

Fabrication and First Full Characterisation of Timing Properties of 3D Diamond Detectors

Anderlini, Lucio and Bellini, Marco and Corsi, Chiara and Lagomarsino, Stefano and Lucarelli, Chiara and Passaleva, Giovanni and Sciortino, Silvio and Veltri, Michele. Fabrication and First Full Characterisation of Timing Properties of 3D Diamond Detectors. Instruments. 2021. doi:10.3390/instruments5040039

work page doi:10.3390/instruments5040039 2021
[6]

and others

Bachmair, F. and others. A 3D diamond detector for particle tracking. Nucl. Instrum. Meth. A. 2015. doi:10.1016/j.nima.2015.03.033

work page doi:10.1016/j.nima.2015.03.033 2015
[7]

Extended theorems for signal induction in particle detectors VCI 2004

Riegler, W. Extended theorems for signal induction in particle detectors VCI 2004. Nucl. Instrum. Meth. A. 2004. doi:10.1016/j.nima.2004.07.129

work page doi:10.1016/j.nima.2004.07.129 2004
[8]

Currents induced by electron motion

Ramo, Simon. Currents induced by electron motion. Proc. Ire. 1939. doi:10.1109/JRPROC.1939.228757

work page doi:10.1109/jrproc.1939.228757 1939
[9]

Currents to conductors induced by a moving point charge

Shockley, W. Currents to conductors induced by a moving point charge. J. Appl. Phys. 1938. doi:10.1063/1.1710367

work page doi:10.1063/1.1710367 1938
[10]

Induced signals in resistive plate chambers

Riegler, Werner , doi =. Induced signals in resistive plate chambers. Nucl. Instrum. Meth. A , keywords =

work page
[11]

Gatti and G

E. Gatti and G. Padovini and V. Radeka , doi =. Signal evaluation in multielectrode radiation detectors by means of a time dependent weighting vector. Nuclear Instruments and Methods in Physics Research , month =. 1982

work page 1982
[12]

Riegler, W. , doi =. An application of extensions of the Ramo–Shockley theorem to signals in silicon sensors. Nucl. Instrum. and Meth A , month =. arXiv , arxivId =:1812.07570 , issn =

work page internal anchor Pith review Pith/arXiv arXiv
[13]

Resistive electrodes and particle detectors: Modelling and measurements of novel detector structures

Janssens, Djunes. Resistive electrodes and particle detectors: Modelling and measurements of novel detector structures. 2024

work page 2024
[14]

Optimization of 3D diamond detectors with graphitized electrodes based on an innovative numerical simulation

Anderlini, Lucio and Bombini, Alessandro and Buti, Clarissa and Janssens, Djunes and Lagomarsino, Stefano and Passaleva, Giovanni and Veltri, Michele. Optimization of 3D diamond detectors with graphitized electrodes based on an innovative numerical simulation. JINST. 2026. doi:10.1088/1748-0221/21/01/P01008. arXiv:2505.12874

work page doi:10.1088/1748-0221/21/01/p01008 2026
[15]

Anderlini and M

L. Anderlini and M. Bellini and C. Corsi and S. Lagomarsino and C. Lucarelli and G. Passaleva and S. Sciortino and M. Veltri , keywords =. A 4D diamond detector for HL-LHC and beyond , journal =. 2022 , issn =. doi:https://doi.org/10.1016/j.nima.2022.167230 , url =

work page doi:10.1016/j.nima.2022.167230 2022