Next: Status of this document Up: Precompound model. Previous: Reaction initial state.   Contents

Subsections

# Simulation of pre-compound reaction

The precompound stage of nuclear reaction is considered until nuclear system is not an equilibrium state. Further emission of nuclear fragments or photons from excited nucleus is simulated using an equilibrium model.

## Statistical equilibrium condition

In the state of statistical equilibrium, which is characterized by an eqilibrium number of excitons , all three type of transitions are equiprobable. Thus is fixed by . From this condition we can get

 (28.1)

## Level density of excited (n-exciton) states

To obtain Eq. () it was assumed an equidistant scheme of single-particle levels with the density , where is the level density parameter, when we have the level density of the -exciton state as

 (28.2)

## Transition probabilities

The partial transition probabilities changing the exciton number by is determined by the squared matrix element averaged over allowed transitions and the density of final states , which are really accessible in this transition. It can be defined as following:

 (28.3)

The density of final states were derived in paper [2] using the Eq. () for the level density of the -exciton state and later corrected for the Pauli principle and indistinguishability of identical excitons in paper [3]:
 (28.4)

 (28.5)

and
 (28.6)

where and it was taken to be equal zero. To avoid calculation of the averaged squared matrix element it was assumed [1] that transition probability is the same as the probability for quasi-free scattering of a nucleon above the Fermi level on a nucleon of the target nucleus, i. e.
 (28.7)

In Eq. () the interaction volume is estimated as , with the De Broglie wave length corresponding to the relative velocity , where is nucleon mass and fm.

The averaging in is further simplified by

 (28.8)

For we take approximation:
 (28.9)

where factor was introduced to take into account the Pauli principle. It is given by
 (28.10)

for and
 (28.11)

for .

The free-particle proton-proton and proton-neutron interaction cross sections are estimated using the equations [4]:

 (28.12)

and
 (28.13)

where cross sections are given in mbarn.

The mean relative kinetic energy is needed to calculate and the factor was computed as , where mean kinetic energies of projectile nucleons and target nucleons , respecively.

Combining Eqs. () - () and assuming that are the same for transitions with and we obtain for another transition probabilities:

 (28.14)

and
 (28.15)

## Emission probabilities for nucleons

Emission process probability has been choosen similar as in the classical equilibrium Weisskopf-Ewing model [5]. Probability to emit nucleon in the energy interval is given

 (28.16)

where is the inverse (absorption of nucleon ) reaction cross section, and are nucleon spin and reduced mass, the factor takes into account the condition for the exciton to be a proton or neutron, and are level densities of nucleus after and before nucleon emission are defined in the evaporation model, respectively and is the excitation energy of nucleus after fragment emission.

## Emission probabilities for complex fragments

It was assumed [1] that nucleons inside excited nucleus are able to "condense" forming complex fragment. The "condensation" probability to create fragment consisting from nucleons inside nucleus with nucleons is given by

 (28.17)

where and are fragment and nucleus volumes, respectively. The last equation was estimated [1] as the overlap integral of (constant inside a volume) wave function of independent nucleons with that of the fragment.

During the prequilibrium stage a "condense" fragment can be emitted. The probability to emit a fragment can be written as [1]

 (28.18)

where
 (28.19)

is the single-particle density for complex fragment , which is obtained by assuming that complex fragment moves inside volume in the uniform potential well whose depth is equal to be , and the factor garantees correct isotopic composition of a fragment .

## The total probability

This probability is defined as

 (28.20)

where total emission probabilities to emit fragment can be obtained from Eqs. () and () by integration over :
 (28.21)

## Calculation of kinetic energies for emitted particle

The equations () and () are used to sample kinetic energies of emitted fragment.

## Parameters of residual nucleus.

After fragment emission we update parameter of decaying nucleus:

 (28.22)

Here is the evaporated fragment four momentum.

Next: Status of this document Up: Precompound model. Previous: Reaction initial state.   Contents