
2
populated final states are parameterized in the exciton
model [9, 12, 13]. This technique enables a more realis-
tic spin transfer to the residual nucleus, while the whole
pre-equilibrium strength can be determined by the more
established exciton model framework.
Although this combined approach compensates defi-
cient information of angular momentum transfer in the
exciton model, it is insufficient to provide individual
contributions from different particle-hole configurations
to the total pre-equilibrium energy spectrum. It is
known that deformed nuclei at relatively low excitation
energies show collective behavior, which can be evalu-
ated by the Quasi-particle Random Phase Approxima-
tion (QRPA) [14, 15], as shown by Kerveno et al. [9]. This
collective excitation can be interpreted as an effective
enhancement in the partial state density for 1-particle-
1-hole configurations. Ergo, incorporating a collective
enhancement for the 1-particle-1-hole state density into
the exciton model may offer better modeling of the en-
tire nuclear reaction occurring in highly deformed nuclei
such as the actinides. Crucially, this procedure can be in-
tegrated into the Hauser-Feshbach theory which follows
the statistical decay of the residual nucleus.
In this paper we study this combined practical tech-
nique. We propose an increase to the 1-particle-1-hole
state density used in the exciton model and include it
in the Los Alamos statistical model framework, CoH3
[16, 17]. We study the impact of this enhancement in the
context of neutron-induced reactions on 239Pu. We use
feedback from Lawrence Livermore National Laboratory
(LLNL) pulsed-sphere neutron-leakage spectra to set the
magnitude of the enhancement factor and find that this
scale factor is significantly above unity. We present the
changes to the cross sections in the results section and
summarize our findings in the final section.
II. THEORY
A. Exciton model
We employ the two-component exciton model [18, 19],
which distinguishes neutron and proton in the particle-
hole configurations. This is denoted schematically in
Fig. 1. Since this model has been well established and ex-
tensively applied to particle emission data analysis, only
a brief description of some of the relevant parts of the
model is given below.
We denote the particle-hole configuration by c, which
abbreviates the number of particles and holes in the
neutron and proton shells as c≡(pν, hν, pπ, hπ). We
also define the total number of excitons, nν=pν+hν,
nπ=pπ+hπ, and nt=nν+nπ. For a particle hav-
ing z-protons and n-neutrons emitted in output channel
b, the residual configuration will be designated by cb,
that stands for pπ−zand pν−n. In the case of an
incident neutron on a target system with Z-protons and
N-neutrons, the composite system would be the nucleus
… To equilibrium …
3 excitons 5 excitons
Fermi Energy
FIG. 1. (Color Online) A schematic depiction of the first
few stages of the 2-component exciton model from an initial
excitation with a neutron. The particles, in this case nucleons
(neutrons and protons), are shown as filled circles with holes
indicated by open circles. The solid lines represent equally
spaced single-particle states.
— before compound nucleus formation — (Z,N+1), and
the residual system might be (Z,N) after emission of the
neutron, e.g. in the case of inelastic scattering.
For the pre-equilibrium nuclear reaction, (a, b), with
input channel, a, and output channel, b, the emission
rate of the outgoing particle bis written as
Wb(c, E, b) = 2sb+ 1
π2~3µbσCN
b(b)b
ω(cb, U)
ω(c, E)fFW ,(1)
where Eis the total energy of the composite system,
Uis the excitation energy in the residual nucleus, and
ω(c, E) is the composite state density at the excitation
energy E. A commonly used step function, fFW , is em-
ployed to limit the hole state configuration within the
potential depth [20]. The values of ,sb, and µb, denote
the emission energy, the intrinsic spin of particle b, and
the reduced mass respectively. The compound formation
cross section for the inverse reaction calculated by the
particle transmission coefficient is σCN
b(b).
The pre-equilibrium emission takes place at different
particle-hole configurations, which is characterized by the
occupation probability P(c) and its lifetime τ(c). The
observed energy-differential cross section is a convolution
of all the configurations
dσ
db
=σCN
a(a)X
c
P(c)τ(c)Wb(c, E, b),(2)
where σCN
ais the compound nucleus formation cross sec-
tion for channel a.
We employ the τ(c) calculation proposed by
Kalbach [21] and adopt the closed-form expression for
P(c). The most important ingredients of this model are
the single-particle state densities, g, and the effective av-
erage squared matrix element M2for the two-body inter-
action. The effective average squared matrix element is
considered as an adjustable model parameter in the exci-
ton model, and often phenomenologically parameterized