The $β$-decay spectrum of Tritiated graphene: combining nuclear quantum mechanics with Density Functional Theory

This paper presents a multi-methodological study combining Density Functional Theory with nuclear quantum mechanics to analyze how graphene substrates influence the $\beta$-decay spectrum of tritium, offering critical insights for optimizing hosting materials in future neutrino experiments.

The Gist

The Big Picture: Catching a Ghost with a Trampoline

Imagine you are trying to weigh a ghost. In the world of physics, neutrinos are these ghosts. They are tiny, invisible particles that barely interact with anything else. Scientists believe they have mass, but they don't know exactly how heavy they are.

To figure this out, scientists look at Tritium (a heavy version of hydrogen) as it decays. When Tritium decays, it turns into Helium, shoots out an electron, and releases a neutrino. By measuring the speed of that electron very precisely, scientists can calculate the weight of the missing neutrino.

The paper you asked about is about a specific experiment called PTOLEMY. Instead of using gas, this experiment plans to stick Tritium atoms onto a sheet of graphene (a material made of carbon atoms arranged in a honeycomb pattern, like a microscopic chicken wire).

The authors of this paper asked a crucial question: "What happens to the electron's speed when the Tritium is stuck to this carbon sheet, rather than floating freely in a vacuum?"

The Problem: The "Sudden" Change

To understand their answer, imagine a game of musical chairs, but with a twist.

1. The Setup (Before the Decay): A Tritium atom is sitting comfortably on the graphene sheet. It's holding hands with the carbon atoms. The electrons in the system are dancing around in a specific, happy pattern. This is the "ground state."
2. The Event (The Decay): Suddenly, the Tritium nucleus changes into a Helium nucleus. This happens incredibly fast—faster than a blink of an eye. It's like a person in a chair suddenly turning into a different person with a different weight and shape.
3. The Confusion (The Aftermath): Because the change happened so fast, the electrons don't have time to react. They are still dancing to the "Tritium music," even though the nucleus is now "Helium." This creates a chaotic, excited state.

The paper tries to figure out exactly how this chaos affects the electron that gets shot out.

The Three Scenarios (The "What Ifs")

The researchers used powerful computer simulations (called Density Functional Theory) to model three different ways this situation could play out:

• Scenario A: The "Freeze-Frame" (Sudden Approximation)
  Imagine taking a photo of the electrons right at the moment of the switch. The electrons are frozen in their old positions. In this scenario, the new Helium atom feels a very strong pull from the carbon sheet because the electrons haven't moved to shield it yet. It's like the Helium is a magnet suddenly appearing on a metal plate before the metal has time to adjust.
• Scenario B: The "Slow Follow" (Semi-Sudden Approximation)
  Imagine the electrons are a little faster. As the Helium moves, one electron decides to tag along with it immediately. Now, the Helium is a bit less "naked" and feels a slightly different pull from the sheet.
• Scenario C: The "Relaxed" (Adiabatic Approximation)
  Imagine the electrons have enough time to calm down and rearrange themselves perfectly around the new Helium atom. In this case, the Helium becomes a neutral, happy atom that doesn't want to stick to the sheet at all. It's like a guest who has settled in and decides to leave the party.

What They Found

The authors discovered that it matters which scenario is true.

• The Shape of the Signal: The speed of the outgoing electron creates a "spectrum" (a graph of energy). If the Helium stays stuck to the sheet (Scenarios A and B), the graph looks like a staircase with distinct steps. If the Helium flies away immediately (Scenario C), the graph looks like a smooth slide.
• The "End-Point": The most important part of the graph is the very top edge (the end-point), where the neutrino mass is hidden. The paper shows that the presence of the graphene sheet shifts this edge significantly compared to a vacuum.
• The "Kick": After the decay, the Helium atom gets a "kick" from the reaction. The authors simulated what happens next: the Helium bounces off the graphene sheet and flies away, transferring some energy to the carbon atoms (making them vibrate). They found that while this creates a lot of heat in their tiny computer model, in a real experiment, the sheet has time to cool down between decays.

Why This Matters

The paper concludes that you cannot ignore the graphene sheet.

If scientists build the PTOLEMY experiment and assume the Tritium behaves exactly like it does in empty space, they will get the wrong answer for the neutrino's mass. The graphene changes the rules of the game.

The authors have built a new "theoretical toolbox" that combines nuclear physics (the decay) with solid-state physics (the graphene sheet). They are essentially saying: "To catch the neutrino ghost, we first need to understand exactly how the carbon trampoline changes the dance of the electron."

Summary in a Nutshell

• Goal: Measure the weight of a neutrino using Tritium on a graphene sheet.
• Challenge: The graphene sheet changes how the Tritium decays and how the electron flies out.
• Method: The authors used supercomputers to simulate the decay under three different "time-speed" assumptions (frozen electrons, tag-along electrons, and relaxed electrons).
• Result: The graphene sheet creates a unique "signature" in the electron's energy that is very different from empty space. Ignoring this would ruin the experiment.
• Next Step: Future experiments need to use these new calculations to ensure they are measuring the neutrino correctly, not just the effect of the carbon sheet.

Technical Summary

Technical Summary: The $\beta$-Decay Spectrum of Tritiated Graphene

Problem Statement
Future neutrino mass experiments, such as PTOLEMY, aim to utilize Tritium ($^3$H) adsorbed on a solid graphene substrate to achieve high event rates and minimize energy losses. However, the presence of a solid-state substrate introduces complex many-body effects that alter the $\beta$-decay spectrum compared to the vacuum case. Specifically, the localization of Tritium near Carbon sites modifies its wave function, and the substrate's degrees of freedom actively participate in the reaction. Current theoretical frameworks often assume the interactions of Tritium and Helium ($^3$He) with graphene are identical or rely on standard adiabatic approximations that fail to capture the non-adiabatic nature of the nuclear transition ($T \to He$). This paper addresses the need for a rigorous theoretical description of the $\beta$-decay rate in the presence of a solid substrate to optimize hosting materials and accurately interpret spectral data for neutrino mass measurements.

Methodology
The authors employ a multi-methodological approach combining Density Functional Theory (DFT) with full quantum treatments of nuclear states to bridge high-energy particle physics and low-energy condensed matter physics. The framework follows the temporal sequence of the decay process (illustrated in Fig. 1):

1. Initial State (Pre-Decay): The interaction potential of Tritium on graphene is calculated using standard DFT within the Born-Oppenheimer (BO) approximation. The study considers various loading levels (3.1% to 100%), spatial distributions (symmetric, dimers), and magnetization states. Three prescriptions for handling nuclear coordinates are tested: full relaxation, fixed Carbon scaffold with relaxed Tritium, and fully frozen nuclei except for the moving Tritium.
2. The Decay Event (Non-Adiabatic Transition): Recognizing that the decay time scale ($t_\beta \approx 10^{-18}$ s) is much shorter than the electronic relaxation time ($t_r \approx 10^{-16} - 10^{-14}$ s), the BO approximation is invalid immediately post-decay. The authors implement a "sudden approximation" and a "semi-sudden approximation." In these schemes, the electronic density is "frozen" at the configuration of the initial Tritium state, while the nuclear charge changes from $Z=1$ (Tritium) to $Z=2$ (Helium). This requires non-standard DFT post-processing to calculate non-self-consistent potentials where the electronic density does not relax to the new nuclear configuration.
3. Final State (Post-Decay): The potential felt by the Helium nucleus is evaluated under three schemes:
   • Sudden: Electronic density frozen at the equilibrium Tritium position (resulting in a $He^{++}$ extraction).
   • Semi-sudden: Electronic density frozen but evaluated at the instantaneous Helium position (resulting in a $He^+$ extraction).
   • Adiabatic: Full electronic relaxation to the new ground state (resulting in neutral $He$ extraction).
4. Decay Rate Calculation: The $\beta$-decay spectrum is computed by solving the Schrödinger equation for the nuclear wave functions in the derived potentials. The matrix element is calculated using the initial Tritium ground state and the final Helium states (both bound and continuous).
5. Long-Term Dynamics: Born-Oppenheimer Molecular Dynamics (BO-MD) are performed to simulate the relaxation of the system after electronic relaxation, tracking the release of Helium and energy dissipation into the substrate.

Key Contributions

• Novel Theoretical Scheme: The paper outlines a computational framework that combines DFT with non-conventional "sudden" and "semi-sudden" approximations to treat the extreme non-adiabaticity of the $T \to He$ transition.
• Potential Characterization: The authors provide detailed interaction potentials for both Tritium and Helium on graphene, demonstrating that the potential landscape changes drastically depending on the nuclear charge and the treatment of electronic relaxation.
• Spectral Analysis: The study calculates the resulting $\beta$-electron spectra, explicitly accounting for the discrete bound states of the final Helium nucleus and the continuous spectrum, which differ significantly from vacuum predictions.

Results

• Potential Differences: The binding potential for Tritium is deep (binding energy $E_b \approx -0.5$ to $-4.3$ eV depending on loading and magnetization). Immediately after decay, the Helium potential varies significantly by scheme:
  • Sudden/Semi-sudden: The potentials are attractive, allowing for bound states of the Helium ion ($He^+$ or $He^{++}$).
  • Adiabatic: The potential is purely repulsive, as the neutral Helium atom is stable and unbound from the substrate.
• Spectral Features: The calculated electron spectra exhibit distinct features near the end-point region.
  • The "sudden" and "semi-sudden" schemes predict multiple discrete peaks in the near-end-point region, corresponding to transitions into different Helium bound states.
  • The separation between the end-point and the onset of the continuous spectrum differs between schemes. The semi-sudden onset is closer to the end-point than the sudden one.
  • The adiabatic scheme yields a spectrum with no bound state contributions, only a continuous spectrum.
• Energy Shifts: The difference between the predicted end-point energy in the graphene substrate and the vacuum expectation is extremely large (on the order of eV), far exceeding the projected energy resolution of the PTOLEMY experiment (~100 meV).
• Dynamics: BO-MD simulations indicate that Helium detaches from the graphene lattice within a few femtoseconds. The released kinetic energy (1.2–3.5 eV) dissipates into the substrate over picoseconds. The simulations suggest that under typical conditions, the release of Helium does not cause structural damage (vacancy creation) to the graphene, provided the kinetic energy does not exceed ~3 eV and the impact is not frontal on weak bonds.

Significance and Claims
The authors claim this work provides a crucial, conceptually new foundation for understanding $\beta$-decay in solid-state environments. By quantifying the substantial modifications to the decay spectrum caused by the substrate, the paper argues that ignoring these effects would lead to incorrect interpretations of neutrino mass experiments. The large predicted shift in the end-point energy serves as a clear, experimentally measurable signature of the substrate's influence.

The paper positions itself as a first step toward determining the physics reach of upcoming neutrino mass experiments. It acknowledges limitations, such as the simplification of the problem to a single-particle quantum mechanics approach and the neglect of long-range electrostatic corrections and full vibrational/electronic excitations in the initial state. The authors emphasize that their multi-methodological approach, while innovative, requires further systematic investigation, particularly regarding the decoupling of electronic and nuclear degrees of freedom and the statistical distribution of local environments in real, disordered samples. They conclude that while this study frames the problem, future work must address these theoretical uncertainties to fully constrain the neutrino mass.