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

This paper presents a new formalism and Python-based tools within the G4CMP framework to enable phonon simulations in custom materials, demonstrated through a detailed analysis of silica phonon transport properties for BeEST-style superconducting detector experiments.

The Gist

Imagine you are trying to listen to a very faint whisper in a noisy room. In the world of physics, scientists use special "super-sensitive ears" called superconducting detectors to hear the tiniest whispers of energy from particles. These detectors are so good that they can spot events that are far weaker than what standard physics predicts (what the paper calls "Beyond the Standard Model" physics).

However, to trust what they hear, they need to know exactly how sound travels through the materials inside their detectors. If they don't understand how sound moves, they might mistake background noise for a real discovery.

Here is a simple breakdown of what this paper does:

1. The Problem: A Missing Map

The scientists use a giant digital simulation tool called Geant4 (think of it as a super-complex video game engine for particles). They added a special "mod" to this engine called G4CMP, which helps them simulate how phonons (tiny packets of sound/vibration) move through cold, solid materials.

But, there was a gap. The simulation didn't know how to handle silica (glass/sand), which is a common material used in these experiments. It's like having a map of a city that shows every street except the one you actually live on. Without the right rules for silica, the simulation couldn't accurately predict how vibrations travel through the glass layers in their detectors.

2. The Solution: Building a Rulebook for Glass

This paper is essentially a "user manual" or a "rulebook" for adding silica to the simulation. The authors didn't just guess; they did the heavy math to figure out exactly how silica behaves when it gets cold.

They broke the job down into four main steps, using some creative physics analogies:

• The Elastic Stiffness (The Springs): Imagine the atoms in silica are connected by invisible springs. The paper calculates exactly how stiff those springs are. They figured out how to translate real-world measurements of glass into the specific numbers the computer needs to know how "bouncy" or "stiff" the material is.
• The Sound Speed (The Highway): Different types of sound waves travel at different speeds. The authors mapped out how fast these "vibrational cars" drive through the glass, depending on which direction they are going.
• The Energy Breakdown (The Domino Effect): Sometimes, a high-energy vibration hits a wall and breaks into two smaller vibrations (like a big domino knocking over two smaller ones). The paper provides the math to predict how often this happens in silica.
• The Impurity Scattering (The Potholes): Real glass isn't perfect; it has tiny atomic "potholes" (isotopes) that scatter sound waves. The authors calculated how much these potholes slow down or scatter the vibrations.

3. The Test: The "Shadow" Experiment

How do you know your new rulebook is correct? You test it.

The authors simulated a scenario where they "shook" the bottom of a crystal and watched the "shadows" (called phonon caustics) appear on the top.

• The Analogy: Imagine shining a flashlight through a complex, faceted crystal onto a wall. You get a specific pattern of light and dark spots.
• The Result: They ran their new silica simulation and compared the resulting "light patterns" to real photos taken in a lab. The computer-generated patterns matched the real photos perfectly. This proved their new rules for silica were accurate.

4. The Gift to the Community

The most important part of this paper is that it didn't just solve the problem for themselves. They created Python tools (like a set of Lego instructions) that anyone else can use.

If another scientist wants to simulate a new material that isn't in the database yet, they can use these tools to calculate the necessary numbers and add that material to the simulation themselves. They also provided a tutorial on how to calculate the "vibrational fingerprint" (Density of States) of any material.

Summary

In short, this paper is a technical guide that taught a super-computer how to understand glass (silica). By figuring out exactly how sound travels through glass at freezing temperatures, they removed a major source of confusion for scientists looking for new physics. They validated their work by showing the computer's "shadows" matched real-life photos, and then they shared their "instruction manual" with the rest of the scientific community so others can do the same.

Technical Summary

Technical Summary: Integration of Silica in G4CMP for Phonon Simulations

Problem Statement
Superconducting detectors with sub-eV energy resolution are critical for setting limits on Beyond the Standard Model (BSM) physics, particularly for low-energy events. However, accurately modeling background events in these detectors requires precise simulation of phonon transport through various materials, including thermal oxide layers. The G4CMP framework, a Geant4-based extension for condensed matter physics, provides tools for modeling phonon and charge dynamics in cryogenic materials. A specific gap existed in the ability to simulate silica (silicon dioxide) within G4CMP, which is necessary for reducing systematic uncertainties in reconstructing nuclear recoil spectra for experiments such as the BeEST (Beryllium Electron capture in Superconducting Tunnel junctions) experiment. The challenge lay in deriving the necessary material parameters—specifically elastic constants, sound speeds, and scattering rates—for both amorphous silica and $\alpha$-Quartz to enable accurate phonon propagation simulations.

Methodology
The authors developed a technical formalism to calculate and implement the required material parameters for silica in G4CMP. The methodology proceeded through several distinct steps:

1. Elastic Constants and Sound Speeds: The study utilized existing experimental data for second-order (SOECs) and third-order (TOECs) elastic constants for silicon dioxide. Mode-dependent sound speeds were implemented based on experimental results for bulk materials.
2. Derivation of Lamé Parameters: A general formalism was established to express Lamé parameters ($\mu, \lambda$) and higher-order coefficients ($\alpha, \beta, \gamma$) as functions of SOECs and TOECs. By applying the symmetry groups specific to the crystal structures (trigonal space group $P3_221$ for $\alpha$-Quartz and isotropic approximation for amorphous silica), the authors derived specific expressions for these parameters.
3. Anharmonic Downconversion: Using the derived Lamé parameters and the Tamura formalism, the authors calculated the anharmonic downconversion rates. This process models the redistribution of a longitudinal phonon's energy into lower-energy phonons (transverse or mixed modes). A Python program was developed to automate these calculations, requiring only Lamé parameters, sound speeds, and density as inputs.
4. Isotopic Scattering: The authors calculated isotopic scattering rates, which account for phonon scattering due to natural isotopic mass variations. This was modeled using time-dependent perturbation theory, resulting in a scattering rate with a quartic dependence on phonon energy ($\Gamma_{isotopic} \propto \nu^4$). The mass defect coefficient ($\Gamma_{md}$) and volume per atom were calculated for both amorphous silica and $\alpha$-Quartz.
5. Phonon Density of States (DOS): To predict the relative probability of exciting acoustic phonon modes, the phonon DOS was calculated for $\alpha$-Quartz using Quantum Espresso (QE) and Phonopy. The contributions of individual acoustic modes (transverse slow, transverse fast, and longitudinal) were extracted. For amorphous silica, the DOS was approximated as equivalent to $\alpha$-Quartz due to the computational cost of calculating it for thermal oxide and the projected small impact on simulation results.
6. Validation: The implementation was validated by simulating phonon caustics (intensity patterns generated by an isotropic point source) for $\alpha$-Quartz. These simulations were compared against experimental bolometer images from previous literature for two specific crystal orientations ($[1\bar{1}00]$ and $[0001]$).

Key Contributions

• Material Implementation: The successful integration of silica (both amorphous and $\alpha$-Quartz) into the G4CMP framework, providing the necessary parameters for simulating substrate backgrounds in BeEST-style experiments.
• Technical Formalism: The derivation of a general method to express Lamé parameters and TOEC coefficients as functions of elastic constants for specific crystal symmetries, specifically applied to $\alpha$-Quartz.
• Tool Development: The creation of Python-based tools to aid users in calculating phonon parameters (specifically anharmonic downconversion rates and DOS) for custom materials. These tools and the silica implementation are made available on the G4CMP GitHub repository.
• Validation Data: The reproduction of experimental phonon caustic patterns for $\alpha$-Quartz, confirming the accuracy of the simulated phonon propagation.

Results

• Parameter Tables: The paper provides calculated values for Lamé parameters ($\mu, \lambda$) and TOEC coefficients ($\alpha, \beta, \gamma$) for both amorphous silica and two sets of elastic constants for $\alpha$-Quartz.
• Scattering Rates: The isotopic scattering rates were calculated as $1.24 \times 10^{-43} \text{ s}^3$ for amorphous silica and $9.60 \times 10^{-44} \text{ s}^3$ for $\alpha$-Quartz.
• DOS Parameters: For implementation in G4CMP, the relative DOS contributions at 1 THz were determined (transverse slow = 0.543, transverse fast = 0.397, longitudinal = 0.0595), with a maximum acoustic phonon energy ($\omega_a$) of 4.2 THz.
• Caustic Validation: The simulated caustic patterns for $\alpha$-Quartz successfully reproduced the experimental bolometer images for the $[1\bar{1}00]$ and $[0001]$ orientations, demonstrating that the model accurately captures the directional dependence of phonon flux. The amorphous case correctly demonstrated intensity dependent only on path length.

Significance
The paper claims that this work supports the superconducting detector community by enabling the implementation of phonon simulations in custom materials within G4CMP. Specifically, the integration of silica allows for the modeling of background events in thermal oxide layers and phonon propagation across these layers from Si substrate events. By providing a tutorial for calculating phonon DOS and general tools for material integration, the authors aim to aid the repeatability of this work and facilitate the simulation of new materials. The validation against experimental caustic images confirms the reliability of the framework for modeling phonon transport in crystalline silica. The authors note that the G4CMP consortium is available to support community members wishing to repeat this analysis or implement new materials.