# Thermonuclear fusion rates for tritium + deuterium using Bayesian methods

## TL;DR

This paper employs Bayesian methods to accurately determine thermonuclear fusion rates for the tritium-deuterium reaction, incorporating comprehensive parameter uncertainties and challenging previous claims of electron screening effects.

## Contribution

It introduces a Bayesian R-matrix model that includes all relevant parameters and uncertainties, providing more reliable fusion rate estimates for astrophysics and fusion energy applications.

## Key findings

- Reaction rate uncertainties are between 0.2% and 0.6%.
- Reaction rates differ by up to 2.9% from previous results.
- No evidence of electron screening effects was found.

## Abstract

The $^3$H(d,n)$^4$He reaction has a large low-energy cross section and will likely be utilized in future commercial fusion reactors. This reaction also takes place during big bang nucleosynthesis. Studies of both scenarios require accurate and precise fusion rates. To this end, we implement a one-level, two-channel R-matrix approximation into a Bayesian model. Our main goals are to predict reliable astrophysical S-factors and to estimate R-matrix parameters using the Bayesian approach. All relevant parameters are sampled in our study, including the channel radii, boundary condition parameters, and data set normalization factors. In addition, we take uncertainties in both measured bombarding energies and S-factors rigorously into account. Thermonuclear rates and reactivities of the $^3$H(d,n)$^4$He reaction are derived by numerically integrating the Bayesian S-factor samples. The present reaction rate uncertainties at temperatures between $1.0$ MK and $1.0$ GK are in the range of 0.2% to 0.6%. Our reaction rates differ from previous results by 2.9% near 1.0 GK. Our reactivities are smaller than previous results, with a maximum deviation of 2.9% near a thermal energy of $4$ keV. The present rate or reactivity uncertainties are more reliable compared to previous studies that did not include the channel radii, boundary condition parameters, and data set normalization factors in the fitting. Finally, we investigate previous claims of electron screening effects in the published $^3$H(d,n)$^4$He data. No such effects are evident and only an upper limit for the electron screening potential can be obtained.