A new software implementation of the Oslo method with rigorous statistical uncertainty propagation

The Oslo method comprises a set of analysis techniques designed to extract nuclear level density and average $\gamma$-decay strength from a set of excitation-energy tagged $\gamma$-ray spectra. Here we present a new software implementation of the entire Oslo method, called OMpy. We provide a summary of the theoretical basis and derive the essential equations used in the Oslo method. In addition to the functionality of the original analysis code, the new implementation includes novel components such as a rigorous method to propagate uncertainties throughout all steps of the Oslo method using a Monte Carlo approach. The resulting level density and $\gamma$-ray strength have to be normalized to auxiliary data. The normalization is performed simultaneously for both quantities, thus preserving all correlations. The software is verified by the analysis of a synthetic spectrum and compared to the results of the previous implementation, the oslo-method-software.

[1]  Axel,et al.  ELECTRIC DIPOLE GROUND STATE TRANSITION WIDTH STRENGTH FUNCTION AND 7 MEV PHOTON INTERACTIONS. Technical Report No. 30 , 1962 .

[2]  M. Strayer,et al.  The Nuclear Many-Body Problem , 2004 .

[3]  A. Görgen,et al.  Experimental level densities of atomic nuclei , 2015, 1503.03300.

[4]  P. Endt,et al.  Investigation of the 56Fe(n, γ)57Fe and 58Fe(n, γ)59Fe reactions , 1980 .

[5]  T. Tveter,et al.  On the relation between the statistical γ-decay and the level density in 162Dy , 1995 .

[6]  N. Bohr.  Neutron Capture and Nuclear Constitution , 1936, Nature.

[7]  C. Johnson.  Systematics of strength function sum rules , 2015, 1506.04700.

[8]  G. Choudalakis.  Fully Bayesian Unfolding , 2012, 1201.4612.

[9]  M. Guttormsen,et al.  Novel techniques for constraining neutron-capture rates relevant for r-process heavy-element nucleosynthesis , 2019, Progress in Particle and Nuclear Physics.

[10]  R. Chrien,et al.  The nuclear structure of 239U , 1984 .

[11]  P. Dirac.  The Quantum Theory of the Emission and Absorption of Radiation , 1927 .

[12]  R. Chrien,et al.  Observation of the M1 giant resonance by resonance averaging in 106Pd , 1987 .

[13]  V. Weisskopf,et al.  Theoretical Nuclear Physics , 1953 .

[14]  E. Somersalo,et al.  Statistical inverse problems: discretization, model reduction and inverse crimes , 2007 .

[15]  T. Massey,et al.  Test of nuclear level density inputs for Hauser-Feshbach model calculations , 2007, 0704.0916.

[16]  A. G. W. Cameron,et al.  A COMPOSITE NUCLEAR-LEVEL DENSITY FORMULA WITH SHELL CORRECTIONS , 1965 .

[17]  B. V. Kheswa,et al.  Low-energy enhancement and fluctuations of γ-ray strength functions in 56,57Fe: test of the Brink–Axel hypothesis , 2016, 1612.04231.

[18]  L. Bernstein,et al.  RAINIER: A simulation tool for distributions of excited nuclear states and cascade fluctuations , 2017, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment.

[19]  Uhl,et al.  Test of gamma-ray strength functions in nuclear reaction model calculations. , 1990, Physical review. C, Nuclear physics.

[20]  A. Merloni,et al.  X-ray spectral modelling of the AGN obscuring region in the CDFS: Bayesian model selection and catalogue , 2014, 1402.0004.

[21]  F. Feroz,et al.  MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics , 2008, 0809.3437.

[22]  T. Ericson.  A statistical analysis of excited nuclear states , 1959 .

[23]  Tveter,et al.  Observation of Fine Structure in Nuclear Level Densities and gamma -Ray Strength Functions. , 1996, Physical review letters.

[24]  B. A. Brown,et al.  Enhanced low-energy γ-decay strength of Ni70 and its robustness within the shell model , 2018, Physical Review C.

[25]  T. Ericson.  The statistical model and nuclear level densities , 1960 .

[26]  Rainer Storn,et al.  Differential Evolution – A Simple and Efficient Heuristic for global Optimization over Continuous Spaces , 1997, J. Glob. Optim..

[27]  M. P. Hobson,et al.  Importance Nested Sampling and the MultiNest Algorithm , 2013, The Open Journal of Astrophysics.

[28]  Gaël Varoquaux,et al.  The NumPy Array: A Structure for Efficient Numerical Computation , 2011, Computing in Science & Engineering.

[29]  J. E. Midtbø,et al.  Verification of detailed balance for γ absorption and emission in Dy isotopes , 2018, Physical Review C.

[30]  M. Guttormsen,et al.  The first generation of γ-rays from hot nuclei , 1987 .

[31]  F. Käppeler,et al.  Evidence for M1 scissors resonances built on the levels in the quasicontinuum of 163Dy. , 2004, Physical review letters.

[32]  T. Egidy,et al.  Nuclear level densities and level spacing distributions: Part II , 1988 .

[33]  Chrien,et al.  Test of photon strength functions by a method of two-step cascades. , 1992, Physical review. C, Nuclear physics.

[34]  France,et al.  First application of the Oslo method in inverse kinematics , 2018, The European Physical Journal A.

[35]  E. Melby,et al.  Extraction of level density and gamma strength function from primary gamma spectra , 1999, nucl-ex/9910009.

[36]  G. Bartholomew,et al.  Gamma-Ray Strength Functions , 1973 .

[37]  T. Tveter,et al.  The unfolding of continuum γ-ray spectra , 1996 .

[38]  M. Martini,et al.  E1 and M1 strength functions from average resonance capture data , 2017 .

[39]  S. Kailas,et al.  RIPL – Reference Input Parameter Library for Calculation of Nuclear Reactions and Nuclear Data Evaluations , 2009 .

[40]  T. Ross,et al.  First observation of low-energy γ-ray enhancement in the rare-earth region , 2016, 1602.05824.

[41]  M. J. D. Powell,et al.  An efficient method for finding the minimum of a function of several variables without calculating derivatives , 1964, Comput. J..

[42]  A. D. Ayangeakaa,et al.  Examination of the low-energy enhancement of the γ -ray strength function of Fe 56 , 2018, 1801.09668.

[43]  K. Takahashi,et al.  The r-process of stellar nucleosynthesis: Astrophysics and nuclear physics achievements and mysteries , 2007, 0705.4512.