pith. sign in

arxiv: 2510.06595 · v2 · submitted 2025-10-08 · ⚛️ physics.comp-ph

Integration of Silica in G4CMP for Phonon Simulations: Framework and Tools for Material Integration

Caitlyn Stone-Whitehead , Israel Hernandez , Connor Bray , Allison Davenport , Spencer Fretwell , Abigail Gillespie , Joren Husic , Mingyu Li
show 5 more authors
Andrew Marino Kyle Leach Bismah Rizwan Wouter Van De Pontseele Grace Wagner
This is my paper

Pith reviewed 2026-05-18 09:45 UTC · model grok-4.3

classification ⚛️ physics.comp-ph
keywords G4CMPphonon simulationssilicacryogenic detectorsBeESTsuperconducting detectorsmaterial integrationGeant4
0
0 comments X

The pith

A technical formalism and Python tools enable adding custom materials like silica to G4CMP for phonon simulations in cryogenic detectors.

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

This paper presents a technical formalism and Python tools that let users incorporate phonon simulations for their own materials into the G4CMP Geant4 extension. The method is shown through a case study that extracts and implements silica phonon transport properties for modeling substrate backgrounds in BeEST-style experiments with superconducting tunnel junctions. A sympathetic reader would care because the work lowers the barrier for accurate simulation of low-energy events in superconducting qubit and low-threshold detector systems. Making the tools available on the G4CMP repository supports broader community use for background studies in searches for physics beyond the standard model.

Core claim

The paper claims that defining a technical formalism for material integration into G4CMP, paired with Python-based tools, allows phonon and charge dynamics to be simulated in arbitrary custom materials, with silica phonon properties extracted and applied as a concrete demonstration for cryogenic detector background modeling.

What carries the argument

The technical formalism for implementing phonon simulations of custom materials in G4CMP, supported by Python tools for material integration.

If this is right

  • Enables detailed simulation of phonon transport in silica for BeEST-style experiments.
  • Supplies open tools so other researchers can add their own materials to G4CMP.
  • Improves modeling of substrate backgrounds in superconducting qubit and low-threshold detector setups.
  • Extends G4CMP utility for condensed matter simulations in cryogenic environments.

Where Pith is reading between the lines

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

  • The same integration approach could apply to other substrate materials used in quantum sensors or dark matter searches.
  • Validated silica models may help isolate and reduce specific background contributions in ongoing BeEST data runs.
  • Future tests could benchmark the implementation against measured phonon spectra from real silica samples under cryogenic conditions.

Load-bearing premise

The phonon transport properties extracted for silica are sufficiently accurate and complete to model substrate backgrounds in the target cryogenic detector geometry.

What would settle it

Direct comparison of G4CMP phonon propagation predictions using the new silica model against experimental measurements of energy loss or mean free paths in silica at millikelvin temperatures.

Figures

Figures reproduced from arXiv: 2510.06595 by Abigail Gillespie, Allison Davenport, Andrew Marino, Bismah Rizwan, Caitlyn Stone-Whitehead, Connor Bray, Grace Wagner, Israel Hernandez, Joren Husic, Kyle Leach, Mingyu Li, Spencer Fretwell, Wouter Van De Pontseele.

Figure 1
Figure 1. Figure 1: Total density of states (DOS) for α￾Quartz silica (blue line). (Calculation using Quan￾tum Espresso compared to reference (red line) cal￾culation [24]) DOS and band structure calculations with refer￾ence information ( [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The phonon density of states with contri [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Top: phonon caustic image for α-Quartz (SiO2) obtained from G4CMP. The crystal direction [1¯100] is at the center of the pattern, oriented out￾of-page. Bright regions indicate directions of high phonon flux. Bottom: phonon caustic image for α￾Quartz measured in Ref. [28], where the crystal direction [1¯100] is out of page and scan range is ±56◦ . The circular spot at the bottom of the image was reported as… view at source ↗
read the original abstract

Superconducting detectors with sub-eV energy resolution have demonstrated success setting limits on Beyond the Standard Model (BSM) physics due to their unique sensitivity to low-energy events. G4CMP, a Geant4-based extension for condensed matter physics, provides a comprehensive toolkit for modeling phonon and charge dynamics in cryogenic materials. This paper introduces a technical formalism to support the superconducting qubit and low-threshold detector community in implementing phonon simulations in custom materials into the G4CMP. As a case study, we present the results of a detailed analysis of silica phonon transport properties relevant for simulating substrate backgrounds in Beryllium Electron capture in Superconducting Tunnel junctions (BeEST)-style experiments using G4CMP. Additionally, Python-based tools were developed to aid users in implementing their own materials and are available on the G4CMP repository.

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 / 0 minor

Summary. The manuscript introduces a technical formalism and associated Python tooling for integrating custom materials with phonon transport parameters into the G4CMP Geant4 extension. It demonstrates the approach via a case study extracting and implementing silica phonon properties relevant to substrate background modeling in BeEST-style cryogenic detectors, with the tools released on the G4CMP repository.

Significance. If the formalism and tools function as described, the work could meaningfully lower the barrier for the superconducting-qubit and low-threshold-detector communities to perform material-specific phonon simulations. The explicit release of reusable Python code is a concrete strength that supports reproducibility and extension beyond the silica example.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript, the recognition of its potential impact on the superconducting-qubit and low-threshold-detector communities, and the recommendation for minor revision. We appreciate the emphasis placed on the release of reusable Python code for reproducibility.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper's core contribution is the introduction of a technical formalism, Python tooling, and integration framework for adding custom materials to G4CMP phonon simulations, with the silica analysis serving only as a concrete demonstration case study. No load-bearing derivations, equations, or predictions are presented that reduce by construction to parameters fitted inside the same work or to self-citations whose validity depends on the present results. The work is framed as software and implementation support rather than a closed physical derivation, making the contribution self-contained against external benchmarks such as the existing G4CMP codebase and literature phonon parameters.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central contribution rests on the assumption that Geant4's existing phonon transport machinery can be extended by user-supplied material parameters without introducing new physics; no free parameters, axioms, or invented entities are explicitly introduced in the abstract.

axioms (1)
  • domain assumption Geant4 phonon and charge transport models remain valid when new material tables are supplied.
    Implicit in any G4CMP material extension.

pith-pipeline@v0.9.0 · 5715 in / 1133 out tokens · 24596 ms · 2026-05-18T09:45:18.204439+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

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

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

31 extracted references · 31 canonical work pages

  1. [1]

    GEANT4—a simulation toolkit,

    S. Agostinelliet al., “GEANT4—a simulation toolkit,” Nucl. Instrum. Methods Phys. Res. A, vol. 506, no. 3, pp. 250–303, Jul. 2003

    work page 2003

  2. [2]

    G4CMP: A Geant4 extension for phonon and charge transport in cryogenic crystals,

    M. Kelseyet al., “G4CMP: A Geant4 extension for phonon and charge transport in cryogenic crystals,”Nucl. Instrum. Methods Phys. Res. A, vol. 1050, p. 168473,

    work page

  3. [3]

    Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment1027, 166251 (2022) https://doi.org/10.1016/j.nima

    [Online]. Available: https://doi.org/10.1016/j.nima. 2023.168473

    work page doi:10.1016/j.nima 2023

  4. [4]

    Projected sensitivity of the supercdms snolab experiment,

    R. Agneseet al., “Projected sensitivity of the supercdms snolab experiment,”Physical Review D, vol. 95, no. 8, Apr. 2017. [Online]. Available: http://dx.doi.org/10.1103/ PhysRevD.95.082002

    work page 2017

  5. [5]

    Modeling athermal phonons in novel materials using the G4CMP simulation toolkit,

    I. Hernandez, R. Linehan, R. Khatiwada, K. Anyang, D. Baxter, G. Bratrud, E. Figueroa-Feliciano, L. Hsu, M. Kelsey, and D. Temples, “Modeling athermal phonons in novel materials using the G4CMP simulation toolkit,” Nucl. Instrum. Methods Phys. Res. A, vol. 1073, p. 170172, Apr. 2025. [Online]. Available: https://doi.org/10. 1016/j.nima.2024.170172

    work page arXiv 2025

  6. [6]

    Signal processing and spectral modeling for the beest experiment,

    I. Kimet al., “Signal processing and spectral modeling for the beest experiment,”Physical Review D, vol. 111, no. 5, p. 052010, 2025. [Online]. Available: https://doi.org/10.1103/PhysRevD.111.052010

    work page doi:10.1103/physrevd.111.052010 2025

  7. [7]

    Spontaneous decay rates of LA phonons in quasi-isotropic solids,

    S. Tamura, “Spontaneous decay rates of LA phonons in quasi-isotropic solids,”Phys. Rev. B, vol. 31, no. 4, pp. 2574–2582, Feb. 1985

    work page 1985

  8. [8]

    Third-order elastic constants of ge, mgo, and fusedsio 2,

    E. H. Bogardus, “Third-order elastic constants of ge, mgo, and fusedsio 2,”Journal of Applied Physics, vol. 36, no. 8, pp. 2504–2513, 1965. [Online]. Avail- able: https://pubs.aip.org/aip/jap/article/36/8/2504/508649/ Third-Order-Elastic-Constants-of-Ge-MgO-and-Fused

    work page 1965

  9. [9]

    First-principles calculations of second- and third-order elastic constants for single crystals of arbitrary symmetry,

    J. Zhao, J. M. Winey, and Y . M. Gupta, “First-principles calculations of second- and third-order elastic constants for single crystals of arbitrary symmetry,”Phys. Rev. B, vol. 75, no. 9, p. 094105, Mar. 2007

    work page 2007

  10. [10]

    High-temperature phase transitions and elasticity of silica polymorphs,

    D. L. Lakshtanov, S. V . Sinogeikin, and J. D. Bass, “High-temperature phase transitions and elasticity of silica polymorphs,”Phys. Chem. Miner ., vol. 34, no. 1, pp. 11– 22, 2007

    work page 2007

  11. [11]

    Third-order elastic constants of fused silica,

    W. T. Yost and M. A. Breazeale, “Third-order elastic constants of fused silica,”J. Appl. Phys., vol. 44, no. 4, pp. 1909–1913, 1973

    work page 1909

  12. [12]

    ‘Third-order’ elastic coefficients,

    R. F. S. Hearmon, “‘Third-order’ elastic coefficients,” Acta Crystallographica, vol. 6, no. 4, pp. 331–340, Apr 1953. [Online]. Available: https://doi.org/10.1107/ S0365110X53000909

    work page 1953

  13. [13]

    F. I. Fedorov,General Theory of Elastic Waves in Crystals Based on Comparison with an Isotropic Medium. Boston, MA: Springer US, 1968, pp. 169–209. [Online]. Available: https://doi.org/10.1007/978-1-4757-1275-9 5

    work page doi:10.1007/978-1-4757-1275-9 1968

  14. [14]

    Multiphonon excitations from dark matter scattering in crystals,

    B. Campbell-Deem, P. Cox, S. Knapen, T. Lin, and T. Melia, “Multiphonon excitations from dark matter scattering in crystals,”Physical Review D, vol. 101, no. 3, Feb. 2020. [Online]. Available: http://dx.doi.org/10.1103/ PhysRevD.101.036006

    work page 2020

  15. [15]

    Srivastava,The Physics of Phonons

    G. Srivastava,The Physics of Phonons. Taylor & Francis, 1990. [Online]. Available: https://books.google. com/books?id=OE-bHd2gzVgC

    work page 1990

  16. [16]

    Isotope scattering of dispersive phonons in ge,

    S.-i. Tamura, “Isotope scattering of dispersive phonons in ge,”Phys. Rev. B, vol. 27, pp. 858–866, Jan

    work page

  17. [17]

    Available: https://link.aps.org/doi/10.1103/ PhysRevB.27.858

    [Online]. Available: https://link.aps.org/doi/10.1103/ PhysRevB.27.858

    work page

  18. [18]

    Isotope scattering of large-wave-vector phonons in gaas and insb: Deformation-dipole and overlap-shell models,

    ——, “Isotope scattering of large-wave-vector phonons in gaas and insb: Deformation-dipole and overlap-shell models,”Phys. Rev. B, vol. 30, pp. 849–854, Jul

    work page

  19. [19]

    Available: https://link.aps.org/doi/10.1103/ PhysRevB.30.849

    [Online]. Available: https://link.aps.org/doi/10.1103/ PhysRevB.30.849

    work page

  20. [20]

    Alloy scattering of phonons,

    R. Gurunathan, R. Hanus, and G. Snyder, “Alloy scattering of phonons,”Materials Horizons, vol. 7, no. 6, pp. 1452– 1456, Jun. 2020, publisher Copyright: © 2020 The Royal Society of Chemistry

    work page 2020

  21. [21]

    Estimation of the isotope effect on the lattice thermal conductivity of group iv and group iii-v semiconductors,

    D. T. Morelli, J. P. Heremans, and G. A. Slack, “Estimation of the isotope effect on the lattice thermal conductivity of group iv and group iii-v semiconductors,”Phys. Rev. B, vol. 66, p. 195304, Nov 2002. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevB.66.195304

    work page doi:10.1103/physrevb.66.195304 2002

  22. [22]

    Thermal transport in defective and disordered materials,

    R. Hanus, R. Gurunathan, L. Lindsay, M. T. Agne, J. Shi, S. Graham, and G. Jeffrey Snyder, “Thermal transport in defective and disordered materials,”Applied Physics Reviews, vol. 8, no. 3, p. 031311, 08 2021. [Online]. Available: https://doi.org/10.1063/5.0055593

    work page doi:10.1063/5.0055593 2021

  23. [23]

    Beyond the tamura model of phonon-isotope scattering,

    N. H. Protik and C. Draxl, “Beyond the tamura model of phonon-isotope scattering,”Phys. Rev. B, vol. 109, p. 165201, Apr 2024. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevB.109.165201

    work page doi:10.1103/physrevb.109.165201 2024

  24. [24]

    QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials,

    P. Giannozziet al., “QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials,”J. Phys. Condens. Matter, vol. 21, no. 39, p. 395502, 2009

    work page 2009

  25. [25]

    First principles phonon calcula- tions in materials science,

    A. Togo and I. Tanaka, “First principles phonon calcula- tions in materials science,”Scr . Mater ., vol. 108, pp. 1–5, 2015

    work page 2015

  26. [26]

    G4CMP PhononModeDOS: Phonon mode dos tutorial,

    C. StoneWhitehead, “G4CMP PhononModeDOS: Phonon mode dos tutorial,” https://github.com/cstonewhitehead/ G4CMP PhononModeDOS/tree/main/PhononModeDOS Tutorial, 2025, gitHub repository, accessed July 13, 2025

    work page 2025

  27. [27]

    Ab-initio phonon calculation for sio 2 / p6 222 (180) / materials id 6922,

    A. Togo, “Ab-initio phonon calculation for sio 2 / p6 222 (180) / materials id 6922,” 2023, mDR dataset, published May 14, 2023. [Online]. Available: https://mdr.nims.go.jp/ datasets/93d8ebae-5983-49de-b34c-61e055f72d0b

    work page 2023

  28. [28]

    G4cmp — caustics ex- ample,

    M. Kelsey and contributors, “G4cmp — caustics ex- ample,” https://github.com/kelseymh/G4CMP/tree/master/ examples/caustics, 2023, accessed: 2025-05-25

    work page 2023

  29. [29]

    G4cmp crystalmaps: Silica elas- tic constants and sound velocities,

    C. Stone-Whitehead, “G4cmp crystalmaps: Silica elas- tic constants and sound velocities,” https://github.com/ kelseymh/G4CMP/tree/master/CrystalMaps, 2025, includes silica-related material property data files for use with G4CMP. Accessed: 2025-05-25

    work page 2025

  30. [30]

    G4cmp consortium: Cryogenic phonon modeling toolkit,

    G4CMP Consortium, “G4cmp consortium: Cryogenic phonon modeling toolkit,” https://confluence.slac.stanford. edu/x/Qk6fDQ, 2024, sLAC Confluence portal for G4CMP documentation and collaboration. Accessed: 2025-05-25

    work page 2024

  31. [31]

    Phonon focusing in piezoelectric crystals: Quartz and lithium niobate,

    G. L. Koos and J. P. Wolfe, “Phonon focusing in piezoelectric crystals: Quartz and lithium niobate,”Phys. Rev. B, vol. 30, no. 6, pp. 3130–3140, Sep. 1984. [Online]. Available: https://doi.org/10.1103/PhysRevB.30.3130

    work page doi:10.1103/physrevb.30.3130 1984