Molecular effects in low-energy muon transfer from muonic hydrogen to oxygen

This paper presents an efficient computational model that incorporates molecular structure effects to accurately determine low-energy muon transfer cross-sections from muonic hydrogen to molecular oxygen, thereby significantly improving agreement with theoretical calculations and aiding the analysis of the FAMU experiment.

The Gist

The Big Picture: A Cosmic Pinball Game

Imagine a game of cosmic pinball. In this game, we have a tiny, heavy particle called a muon (it's like a super-heavy electron) that has grabbed onto a hydrogen atom (a proton). We call this a "muonic hydrogen" atom.

Now, imagine this muonic hydrogen is zooming around inside a cloud of gas that is mostly hydrogen but has a little bit of oxygen mixed in (like the air we breathe, but mostly hydrogen).

The goal of the FAMU experiment (the big project this paper is about) is to measure the internal structure of the proton (the heart of the hydrogen atom) with extreme precision. To do this, they need to know exactly how the muonic hydrogen behaves when it bumps into oxygen.

Here is the problem: When the muonic hydrogen hits an oxygen molecule, it doesn't just bounce off. It often steals the muon and gives it to the oxygen atom. This is called "muon transfer."

The scientists need to know the speed at which this theft happens. If they get the speed wrong, their measurements of the proton's structure will be wrong.

The Old Way vs. The New Way

The Old Way (The "Frozen" Model):
In previous studies, scientists treated the oxygen molecule like a frozen statue. They assumed the oxygen atoms inside the molecule were glued together and couldn't move or vibrate. They also assumed the muonic hydrogen was moving at a perfectly predictable, average speed (like a calm crowd of people walking).

Based on this "frozen statue" idea, they calculated how often the muon transfer happens. But when they compared their math to the actual data from the FAMU experiment, the numbers didn't quite line up. It was like trying to predict how a ball bounces off a trampoline by pretending the trampoline is a solid concrete floor.

The New Way (The "Dancing" Model):
This paper says: "Wait a minute! Oxygen molecules aren't statues. They are dancers."

1. The Oxygen is Alive: An oxygen molecule ($O_2$) is two atoms holding hands. They are constantly spinning, vibrating, and jiggling. When the muonic hydrogen approaches, it's not hitting a frozen target; it's hitting a wiggly, spinning target.
2. The Muon is Wobbly: The muonic hydrogen atoms themselves aren't moving in a perfect, calm line. Because they are bumping into other gas molecules, their speed distribution is a bit "wobbly" and uneven, not a perfect smooth curve.

The authors built a new, super-complex computer model that treats the oxygen molecule like a dancing partner and the muonic hydrogen like a wobbly skater.

The "Recipe" for the New Model

To fix the math, the team did three main things:

1. They looked at the "Dance Moves": They calculated how the oxygen molecule spins and vibrates at different temperatures. They realized that the internal motion of the oxygen atoms changes the speed at which the muon gets stolen.
2. They smoothed out the "Wobble": They used a sophisticated method (called a Boltzmann-Lorentz kinetic equation) to track exactly how the muonic hydrogen atoms lose energy and change direction as they bounce around the gas.
3. They Re-calculated the "Theft Rate": By feeding these realistic "dance" and "wobble" details into their equations, they derived a new cross-section.
   • Analogy: Think of the "cross-section" as the size of the target. If the target is a frozen statue, it's a certain size. But if the target is a spinning, wiggling dancer, the "effective" size changes depending on how fast they are spinning.

The Results: Why It Matters

When they compared their new "Dancing Model" to the old "Frozen Model":

• The Old Model predicted the muon transfer would happen most often at a specific low speed (around 63 "meV" of energy).
• The New Model showed that because the oxygen is wiggling, the peak of the transfer actually happens at a slightly higher speed (around 73 "meV").

Why is this a big deal?
The theoretical predictions (the math done by other scientists using quantum mechanics) had been predicting a peak around 73 "meV" all along. The old experimental data (based on the frozen model) disagreed with the theory.

By accounting for the "dance" of the oxygen molecule, the experimental data now matches the theory perfectly.

The Bottom Line

This paper is a crucial update for the FAMU experiment.

• Before: The scientists were trying to measure the proton's size, but their "ruler" (the muon transfer rate) was slightly bent because they ignored the fact that oxygen molecules wiggle.
• Now: They have straightened the ruler. They have a much more accurate map of how muons move between hydrogen and oxygen.

This ensures that when the FAMU experiment finally measures the Zemach radius (a specific measurement of the proton's "fuzziness" or size), the result will be incredibly precise. It's like finally realizing you were measuring a spinning top with a ruler meant for a still block, and now you can measure the top correctly.

In short: They stopped treating oxygen like a rock and started treating it like a living, breathing, spinning molecule, and that small change fixed a big problem in our understanding of the universe's smallest building blocks.

Technical Summary

1. Problem Statement

The paper addresses the need for precise modeling of the muon transfer reaction from muonic hydrogen ($p\mu$) to molecular oxygen ($O_2$):
$$p\mu^- + O_2 \rightarrow p + (O\mu^-) + O$$
This reaction is critical for the FAMU experiment, which aims to measure the hyperfine splitting in the ground state of muonic hydrogen to determine the proton's Zemach radius.

Key Challenges Identified:

• Energy Dependence: The rate of muon transfer varies significantly with the kinetic energy of the $p\mu$ atom.
• Molecular Structure Neglect: Previous experimental analyses (e.g., Ref. [13]) treated the oxygen molecule as a collection of "frozen" nuclei, ignoring the internal degrees of freedom (vibrational and rotational states) and the motion of the oxygen nucleus within the molecule.
• Discrepancy: Theoretical calculations of the transfer rate showed discrepancies with experimental data derived under the "frozen nucleus" assumption, particularly regarding the position of the resonance peak in the energy dependence.
• Non-Maxwellian Distribution: The kinetic energy distribution of thermalized $p\mu$ atoms in an $H_2/O_2$ mixture deviates from a standard Maxwell-Boltzmann distribution due to inelastic collisions with $H_2$ molecules.

2. Methodology

The authors developed a comprehensive computational model to extract the muon transfer cross-section from experimental data while accounting for molecular structure.

A. Kinetic Modeling of the FAMU Experiment

• Boltzmann-Lorentz Equation: The propagation of $p\mu$ atoms in the gas mixture is modeled using a kinetic equation for the energy and spin distribution ($n_F(E, T, t)$).
• Multigroup Method: The energy range is discretized into intervals to solve the differential equations numerically.
• Stationary Distribution: The model calculates the stationary energy distribution of $p\mu$ atoms after they thermalize. It confirms that this distribution deviates from Maxwell-Boltzmann due to energy loss via rotational excitation of $H_2$ molecules, though this deviation has a negligible effect on the final transfer rate extraction.

B. Derivation of the Transfer Rate Equation

• Three-Body Process: The authors distinguish between the translational motion of the $O_2$ molecule and the internal motion of the oxygen nucleus within the molecule.
• Convolution Integral: Instead of assuming a fixed target, they derive a 3D integral equation (Eq. 12) relating the observable transfer rate $\Lambda(T)$ to the unknown cross-section $\sigma(v)$. This equation convolves:
  1. The velocity distribution of $p\mu$ atoms ($f_A$).
  2. The velocity distribution of $O_2$ molecules ($f_M$).
  3. The velocity distribution of the oxygen nucleus relative to the molecule's center of mass ($f_N$), which accounts for ro-vibrational states.
• Regularization: The resulting integral equation is ill-posed. The authors solve it using Tikhonov regularization via Truncated Singular Value Decomposition (TSVD) applied to a discretized linear system derived from Gauss-Hermite quadrature.

C. Data Sources

• Experimental data from the FAMU collaboration (measurements of transfer rates at temperatures $70\text{ K} \leq T \leq 336\text{ K}$).
• Theoretical momentum distributions of $O_2$ nuclei from Ref. [26], accounting for temperature-dependent population of excited ro-vibrational states.

3. Key Contributions

1. Incorporation of Molecular Structure: The study is the first to rigorously include the internal degrees of freedom (nuclear motion within $O_2$) in the extraction of the muon transfer cross-section from experimental data.
2. Improved Computational Model: Development of an efficient kinetic model that accounts for the non-Maxwellian energy distribution of $p\mu$ atoms and the specific collision dynamics in $H_2/O_2$ mixtures.
3. Resolution of Theory-Experiment Discrepancy: By correcting the "frozen nucleus" approximation, the authors bridge the gap between experimental data and advanced theoretical calculations.
4. Validation of Approximations: The Appendix demonstrates that the interaction between $p\mu$ and $O_2$ does not significantly alter the thermal population of the $O_2$ rotational states, justifying the use of equilibrium distributions in the kernel calculation.

4. Key Results

• Cross-Section Extraction: The authors successfully extracted the muon transfer cross-section $\sigma(v)$ as a function of relative velocity.
• Shift in Resonance Peak:
  • Previous work (ignoring molecular effects) placed the peak transfer rate at 63 meV.
  • The new model, accounting for the motion of the oxygen nucleus, shifts the peak to 73 meV.
  • This shift significantly improves the agreement with theoretical predictions (specifically Ref. [9] and Ref. [12]), which predicted a peak in this higher energy range due to molecular effects "stretching" the curve toward higher collision energies.
• Rate Calculation: The energy-dependent transfer rate $\lambda(E)$ was calculated and normalized to Liquid Hydrogen Density (LHD). The results show a ~7% temperature dependence between 80 K and 300 K.
• Uncertainty Analysis: The study quantifies statistical uncertainties from experimental data and estimates systematic uncertainties arising from the lack of data at temperatures $T > 336\text{ K}$. The cross-section is considered reliable for relative velocities up to $5000\text{ m/s}$ (approx. 0.1 eV).

5. Significance

• FAMU Experiment Optimization: The derived cross-sections and transfer rates are essential for accurately modeling the time distribution of X-rays in the FAMU experiment. This allows for better optimization of the laser excitation parameters and more precise extraction of the hyperfine splitting frequency.
• Proton Structure: By improving the accuracy of the muon transfer modeling, the uncertainty in the determination of the proton Zemach radius is reduced, contributing to the resolution of the "proton radius puzzle."
• General Applicability: The computational techniques and the upgraded model are designed to be applicable to other muon transfer processes (e.g., to Carbon, as mentioned in Ref. [16]), making this a foundational work for muonic atom physics in gaseous mixtures.

In summary, this paper upgrades the understanding of low-energy muon transfer by moving beyond the "frozen nucleus" approximation, demonstrating that molecular structure effects are crucial for aligning experimental data with quantum mechanical theory.