An improved numerical method to compute neutron/gamma deexcitation cascades starting from a high spin state
Abstract
Numerous nuclear processes involve the deexcitation of a compound nucleus through the emission of several neutrons, gammarays and/or conversion electrons. The characteristics of such a deexcitation are commonly derived from a total statistical framework often called “Hauser–Feshbach” method. In this work, we highlight a numerical limitation of this kind of method in the case of the deexcitation of a high spin initial state. To circumvent this issue, an improved technique called the Fluctuating Structure Properties (FSP) method is presented. Two FSP algorithms are derived and benchmarked on the calculation of the total radiative width for a thermal neutron capture on ^{238}U. We compare the standard method with these FSP algorithms for the prediction of particle multiplicities in the deexcitation of a high spin level of ^{143}Ba. The gamma multiplicity turns out to be very sensitive to the numerical method. The bias between the two techniques can reach 1.5 γγ/cascade. Lastly, the uncertainty of these calculations coming from the lack of knowledge on nuclear structure is estimated via the FSP method.
 Authors:
 CEA, DEN, DER, SPRC, Saint Paul lez Durance (France); Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
 CEA, DEN, DER, SPRC, Saint Paul lez Durance (France)
 Publication Date:
 Research Org.:
 Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
 Sponsoring Org.:
 USDOE
 OSTI Identifier:
 1342046
 Report Number(s):
 LLNLJRNL668707
Journal ID: ISSN 00104655
 Grant/Contract Number:
 AC5207NA27344
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Computer Physics Communications
 Additional Journal Information:
 Journal Volume: 201; Journal Issue: C; Journal ID: ISSN 00104655
 Publisher:
 Elsevier
 Country of Publication:
 United States
 Language:
 English
 Subject:
 73 NUCLEAR PHYSICS AND RADIATION PHYSICS; Hauser–Feshbach; statistical model; deexcitation; cascade; FIFRELIN
Citation Formats
Regnier, D., Litaize, O., and Serot, O. An improved numerical method to compute neutron/gamma deexcitation cascades starting from a high spin state. United States: N. p., 2015.
Web. doi:10.1016/j.cpc.2015.12.007.
Regnier, D., Litaize, O., & Serot, O. An improved numerical method to compute neutron/gamma deexcitation cascades starting from a high spin state. United States. doi:10.1016/j.cpc.2015.12.007.
Regnier, D., Litaize, O., and Serot, O. 2015.
"An improved numerical method to compute neutron/gamma deexcitation cascades starting from a high spin state". United States.
doi:10.1016/j.cpc.2015.12.007. https://www.osti.gov/servlets/purl/1342046.
@article{osti_1342046,
title = {An improved numerical method to compute neutron/gamma deexcitation cascades starting from a high spin state},
author = {Regnier, D. and Litaize, O. and Serot, O.},
abstractNote = {Numerous nuclear processes involve the deexcitation of a compound nucleus through the emission of several neutrons, gammarays and/or conversion electrons. The characteristics of such a deexcitation are commonly derived from a total statistical framework often called “Hauser–Feshbach” method. In this work, we highlight a numerical limitation of this kind of method in the case of the deexcitation of a high spin initial state. To circumvent this issue, an improved technique called the Fluctuating Structure Properties (FSP) method is presented. Two FSP algorithms are derived and benchmarked on the calculation of the total radiative width for a thermal neutron capture on 238U. We compare the standard method with these FSP algorithms for the prediction of particle multiplicities in the deexcitation of a high spin level of 143Ba. The gamma multiplicity turns out to be very sensitive to the numerical method. The bias between the two techniques can reach 1.5 γγ/cascade. Lastly, the uncertainty of these calculations coming from the lack of knowledge on nuclear structure is estimated via the FSP method.},
doi = {10.1016/j.cpc.2015.12.007},
journal = {Computer Physics Communications},
number = C,
volume = 201,
place = {United States},
year = 2015,
month =
}

CONTRIBUTION TO THE EXTRAGALACTIC GAMMARAY BACKGROUND FROM THE CASCADES OF VERY HIGH ENERGY GAMMA RAYS FROM BLAZARS
As very high energy (VHE) photons propagate through the extragalactic background light (EBL), they interact with the soft photons and initiate electromagnetic cascades of lower energy photons and electrons. The collective intensity of a cosmological population emitting at VHEs will be attenuated at the highest energies through interactions with the EBL and enhanced at lower energies by the resulting cascade. We calculate the cascade radiation created by VHE photons produced by blazars and investigate the effects of cascades on the collective intensity of blazars and the resulting effects on the extragalactic gammaray background. We find that cascade radiation greatly enhancesmore » 
A general, numerical method for solving twobody, bound state BetheSalpeter equations: finiteenergy solutions of a spin0 and spin1/2 constituent bound state by scalar electrodynamics
Finiteenergy, boundstate solutions are obtained to the BetheSalpeter equation describing a spin0 boson and a spin1/2 fermion with unequal masses that interact via scalar electrodynamics. The equation is solved in the ladder approximation using a numerical method that yields solutions to many, if not all such, twobody, boundstate BetheSalpeter equations.