|Program name||Package id||Status||Status date|
|Package ID||Orig. computer||Test computer|
|NEA-0648/01||IBM 3033||IBM 3033|
The CERBERO-3 code calculates a complete consistent set of binary cross sections, on the basis of the nuclear optical model and of the nuclear statistical model, for incident particles with spin 0 and 1/2. Neutron, proton, alpha particle and gamma-ray competitions can be considered. The present program is designed for use in the incident particle energy region where it may be assumed that shape elastic scattering and compound nucleus absorption are predominant, so that compound nucleus trans- mission coefficients are those of the optical model and the absorp- tion cross section is equivalent to the compound nucleus cross section. The program outputs are the total and compound nucleus cross sections, the shape elastic differential and integrated cross sections, the differential and integrated fluctuation cross sections for the excitation of discrete levels, the differential and inte- grated elastic cross sections, the continuum level excitation func- tions, the total inelastic scattering cross section and the capture cross section.
Only spherical local potentials are allowed.he scattering amplitudes and the other optical model quantities as well as the shape elastic and compound nucleus cross sections are calculated in accordance with the above described optical model potentials, by the usual methods.
The total cross section is given as the sum of the shape elastic and compound nucleus cross sections.
The integrated and differential cross sections for compound nucleus processes may be given according to the three following approximations:
The Hauser-Feshbach theory.
The Hauser-Feshbach theory corrected by the width fluctuation effect.
The Moldauer theory.
Integrated capture cross sections are given according to the various approximations 1, 2, 3 as a difference between the total compound nucleus cross section and the sum of the cross sections for compound elastic and all possible inelastic processes. Therefore, for the purpose of capture cross section calculations, the program should be used only as long as the ratio of the total radiative width to the capture width is roughly unity.
Whenever direct processes may be considered unimportant, the CERBERO program can work up to the threshold for the second particle emission.
All cross sections are calculated in the c.m.s. and are given in barns.
a) Optical model particle transmission coefficients are calculated according to a spherical optical model.
Optical model parameter sets can be given in input or selected by means of acronyms out of a library available internally to the code.
b) The correction for the width fluctuations has been introduced into channels leading both to discrete and to continuum level excitation.
c) Composite level density like Gilbert-Cameron is used but with constant in spin distribution K = .146, when not given in input. For low lying levels between Ecut and U(x), sigma(E)**2 is automatically interpolated between sigma(LEVELS)**2 and sigma(U(x))**2 = K.SQRT(aU(x)).exp(A,2/3). The value
sigma(LEVELS)**2 is obtained by maximum likelihood method to fit known discrete levels distribution.
Sigma(LEVELS)**2 is calculated by the code on the basis of adopted level for each nucleus involved. Alternatively sigma(LEVELS)**2 can also be given in input, in those particular nuclei where additional information is known above Ecut, as far as spin attribution is concerned.
d) Optionally a parity distribution p(pi)=exp(AE+B) according to ref. 1 can also be assumed provided A and B are given in input.
e) Gamma-ray transmission coefficients are calculated according to one or two Lorentzian curves for the E1 photoabsorption cross sections. Peak energy, half maximum width, peak cross section must be given in input for the E1 giant resonance.
The resulting total radiative width is spin and parity dependent. In principle it should not be normalized because the model proved to work satisfactorily. For the purpose of evaluation a normalization constant N (J and pi independent) can be given in input.
f) Q values are calculated from recent mass excess tables (internal to the code) provided by Wapstra in 1978 as a private communication.
The output are average resonance parameters like strength functions (from adopted optical model), radiative width and mean observed level spacing.
Angular distributions are given for compound, shape and total elastic. Total cross section and primary spectra are given for all involved particle and gamma-ray emissions.
Compound nucleus and total cross section from optical model are given at the end together with the percentual difference between compound nucleus cross section and the sum of the contribution of all channels via compound nucleus reaction mechanism.
|Package ID||Status date||Status|
|NEA-0648/01||22-SEP-1981||Tested at NEADB|
G. Reffo: "Theory and Application of Moment Methods in Many Fermion Systems" p. 167, edited by B.J. Dalton, S.M. Grimes, J.P. Vary, S.A. Williams - Plenum 1980.
G. Reffo, F. Fabbri: N.S.E. 66, 251 (1978).
For Optical Model Potential Calculations see for example D.T. Goldmann and C. L. Lubitz, KAPL-2163 (1961).
|Package ID||Computer language|
|File name||File description||Records|
|NEA0648_01.001||CERBERO SOURCE F4||3375|
|NEA0648_01.002||TODAY DUMMY ROUTINE||9|
|NEA0648_01.003||CERBERO SAMPLE PROBLEM INPUT||53|
|NEA0648_01.004||CERBERO SAMPLE PROBLEM OUTPUT||795|
|NEA0648_01.005||JCL AND INFORMATION||56|
Keywords: Hauser-Feshbach theory, absorption, capture, cross sections, elastic scattering, fluctuations, optical models, statistical models.