Performing Bayesian Analyses With AZURE2 Using BRICK: An Application to the 7Be System

Phenomenological R-matrix has been a standard framework for the evaluation of resolved resonance cross section data in nuclear physics for many years. It is a powerful method for comparing different types of experimental nuclear data and combining the results of many different experimental measurements in order to gain a better estimation of the true underlying cross sections. Yet a practical challenge has always been the estimation of the uncertainty on both the cross sections at the energies of interest and the fit parameters, which can take the form of standard level parameters. Frequentist (chi(2)-based) estimation has been the norm. In this work, a Markov Chain Monte Carlo sampler, emcee, has been implemented for the R-matrix code AZURE2, creating the Bayesian R-matrix Inference Code Kit (BRICK). Bayesian uncertainty estimation has then been carried out for a simultaneous R-matrix fit of the He-3 (alpha,gamma)Be-7 and He-3 (alpha,alpha)He-3 reactions in order to gain further insight into the fitting of capture and scattering data. Both data sets constrain the values of the bound state alpha-particle asymptotic normalization coefficients in Be-7. The analysis highlights the need for low-energy scattering data with well-documented uncertainty information and shows how misleading results can be obtained in its absence.

Automated summary

A Markov Chain Monte Carlo sampler emcee has been implemented, creating the Bayesian R-matrix Inference Code Kit (BRICK). Bayesian uncertainty estimation has been carried out for simultaneous R-matrix fit of the He-3 (alpha,gamma)Be-7 and He-3 (alpha,alpha)He-3 reactions to gain further insight.