Charge-state dynamics of barium ions in high-pressure xenon and its implications for Barium-Tagging in $0νββ$ searches

Technical Summary: Charge-State Dynamics of Barium Ions in High-Pressure Xenon

Problem Statement
The detection of neutrinoless double-beta decay ($0\nu\beta\beta$) in high-pressure xenon time projection chambers (TPCs), such as the NEXT experiment, relies heavily on Barium Tagging (BaTa) to achieve background-free sensitivity. This technique involves identifying the daughter barium atom ($^{136}\text{Ba}$) produced in the decay $^{136}\text{Xe} \to ^{136}\text{Ba}^{2+} + 2e^- + (2\bar{\nu}_e)$. However, the experimental feasibility of various BaTa strategies depends critically on the charge state of the daughter ion immediately following the decay. While some approaches assume the ion remains as $\text{Ba}^{2+}$ (allowing for fluorescent molecular capture), others assume rapid electron capture to $\text{Ba}^+$ or neutral Ba (enabling laser-induced fluorescence). Despite the importance of this parameter, quantitative estimates regarding the survival probability of the $\text{Ba}^{2+}$ ion in high-pressure xenon environments have been scarce. Most proposed techniques implicitly assume a specific charge state without rigorous physical validation of the recombination dynamics.

Methodology
The author develops a simplified physical model to estimate the survival probability of the $\text{Ba}^{2+}$ ion in the high-pressure xenon environment (specifically 15 bar at room temperature) characteristic of the NEXT experiment. Rather than attempting a complete microscopic simulation of thermalization and charge exchange, the work derives physically motivated bounds based on order-of-magnitude estimates.

The methodology proceeds through the following steps:

1. Ionization Track Characterization: The model estimates the spatial distribution of electrons generated by the decay. It calculates the mean free path ($\lambda_{ion} \approx 25\,\mu\text{m}$) for the primary $\beta$-electrons to cause the first ionization, establishing the characteristic distance ($r_0$) between the daughter ion and the initial electron cloud.
2. Electron Thermalization: The paper analyzes the energy loss of secondary $\delta$-electrons. It determines that these electrons thermalize extremely rapidly ($\tau_{th} \sim 10^{-13} - 10^{-12}\,\text{s}$) via collisions with xenon atoms, evolving from dense ionization clusters into a dilute cloud of thermal electrons.
3. Recombination Mechanism Evaluation: The study evaluates potential recombination channels:
   • Binary Mechanisms: Radiative recombination, transient molecular complex formation, and collisional charge transfer are analyzed. The paper concludes these are either inefficient due to low electron density or energetically unfavorable (e.g., charge transfer is endothermic by $\sim 2\,\text{eV}$).
   • Three-Body Recombination: This mechanism ($\text{Ba}^{2+} + e^- + \text{Xe} \to \text{Ba}^+ + \text{Xe}$) is identified as the dominant candidate. It relies on a neutral xenon atom to stabilize the recombined ion.
4. Kinetic Estimation: The author estimates the local electron density ($n_e \sim 10^7\,\text{cm}^{-3}$) surrounding the ion after diffusion and compares the Onsager radius ($r_O \approx 1.1\,\mu\text{m}$) and field-balance radius ($r^* \approx 2.4\,\mu\text{m}$) to the ionization distance. This confirms that thermalized electrons can diffuse into the Coulomb-dominated region of the ion before being swept away by the external drift field.
5. Timescale Calculation: Using a conservative order-of-magnitude estimate for the three-body recombination coefficient ($\alpha_3 \sim 10^{-25} - 10^{-24}\,\text{cm}^6\,\text{s}^{-1}$) derived from weakly ionized plasma data, the characteristic recombination time is calculated.

Key Results

• Dominance of Three-Body Recombination: Binary recombination channels are found to be strongly suppressed under NEXT operating conditions. Three-body recombination assisted by neutral xenon atoms constitutes the physically plausible mechanism for converting $\text{Ba}^{2+}$ to $\text{Ba}^+$.
• Recombination Timescale: The estimated characteristic time for the conversion $\text{Ba}^{2+} \to \text{Ba}^+$ is in the range of 0.25 to 2.5 milliseconds.
• Comparison to Detection Times: This millisecond timescale is comparable to the characteristic times required for charge transport and event reconstruction in the NEXT experiment.
• Electron Availability: Despite the initial separation between the ion and the ionization track, thermalized electrons can diffuse into the ion's Coulomb capture radius ($r \lesssim r_O$) within nanoseconds, ensuring a sufficient population of electrons is available for recombination before macroscopic charge extraction occurs.

Significance and Claims
The paper asserts that the charge state of the daughter barium ion should not be regarded as a fixed property but as a dynamical quantity sensitive to the microscopic conditions following the decay. The primary claim is that the assumption of a persistent $\text{Ba}^{2+}$ state is not supported by this analysis; instead, the conversion to $\text{Ba}^+$ is a physically plausible process occurring on timescales relevant to detector operation.

Consequently, the author argues that recombination effects must be explicitly considered in the design and experimental implementation of BaTa techniques. This work provides a physical framework suggesting that strategies relying on the exclusive detection of $\text{Ba}^{2+}$ (such as specific fluorescent molecular indicators) may face challenges if the ion recombines before detection. Conversely, the results support the viability of laser-induced fluorescence techniques targeting $\text{Ba}^+$. The paper concludes that future theoretical and experimental studies must account for these dynamics to optimize BaTa strategies in high-pressure xenon detectors.