Journal of Applied Mathematics and Physics, 2014, 2, 27-31
Published Online April 2014 in SciRes. http://www.scirp.org/journal/jamp
How to cite this paper: Gontchar, I.I., et al. (2014) Dynamical Modeling of the Nuclear Fission Process at Low Excitation En-
ergies. Journal of Applied Mathematics and Physics, 2, 27-31. http://dx.doi.org/10.4236/jamp.2014.25004
Dynamical Modeling of the Nuclear Fission
Process at Low Excitation Energies
I. I. Gontchar1, M. V. Chushnyakova1,2, E. P. Oskin1, E. G. Demina1
1Physics and Chemistry Department, Omsk State Transport University, Omsk, Russia
2Physics Department, Omsk State Technical University, Omsk, Russia
Email: email@example.com om
Received December 2013
Two recipes for modeling the dynamics of the nuclear fission process are known in literature. The
underlying equations contain the driving, dissipative, and random forces. The two recipes are
mostly different in the prescriptions for the driving force. In this work we carefully compare these
driving forces and the resulting fission rates. It turns out that the rates may be very close or
strongly different depending on the value the shell correction to the nuclear deformation energy.
We give arguments in favor of one of the recipes.
Nuclear Fissio n, Dynamical Modeling, Stochastic Differential Equations, Fermi Gas Model
Numerous experiments related to the nuclear fission process [1,2], in particular those aiming for discovering
new superheavy elements [3-5], require practical and reliable numerical modeling of fission dynamics at low
excitation energies. The most important physical quantity characterizing the fission process, although not ob-
servable, is the fission rate which is by definition
Rt N Ntt
is the number of trajectories reached the scission point during the time interval
is the number of trajectories reached the scission point until the time moment
Presently the dynamical modeling of fission is mostly performed by solving numerically the stochastic differ-
rential equations with the white noise. In the simplest one dimensional case these equations read
2dpm pdtfdtT dW
I. I. Gontchar et al.
is the time interval during which the collective coordinate (deformation parameter)
and its con-
is the Wiener process whose increment
obeys the normal distribution with the variance
are the friction and inertia
parameters, respectively. For the sake of simplicity we consider
to be deformation independent
The driving force can be calculated by two methods. The first one is a generalization  of the thermody-
namical approach, which was developed for the one dimensional problem in [7,8] and later was extended to the
multidimensional case in . Note that this method is applied, except , for the high excitation energies when
the shell correction to the nuclear Potential Energy (
, PE) and to the single particle Level Density Parameter
, LDP) is believed does not play a role. The driving force calculated using this approach,
, is defined as a
proper derivative of a thermodynamical potential:
is the Helmholtz free energy,
is the nucleus entropy, and
is its total excitation energy. Within the
framework of the Fermi gas model, the entropy reads
The intrinsic excitation energy
is related to
and PE as follows:
()( )( )()
UqEVqEVq Vq=−=− −
are the smooth (liquid drop) part of PE and the shell correction, respectively.
According to  the LDP can be written as
()( )()( )
reaches unity when
becomes large enough:
. The deformation dependence of the smooth part of the LDP is defined as follows :
( )( )
aq aAaABq= +
is the nucleus mass number,
is the dimensionless nuclear surface area.
Using Equations (5), (7) and (8) we obtain for the driving force the following expression
The prime denotes the derivative with respect to the coordinate, the subscript
indicates that the deriva-
tive must be taken keeping
Within the framework of the second approach [11,12] the driving force
is calculated according the for-
mula, which (using our notations) reads
Equations (10) and (11) look very much different. One hardly can expect to find the fission rates similar if
are used in Equation (4). In order to check whether the rates are different indeed we construct a sche-
matic PE which reproduces the main distinct features of the PE of 236U-nucleus. This schematic
in Figure 1 along with
One sees clearly the double humped structure of the fission barrier which is well known for the actinide nuclei
We modeled the fission process applying Equations (3), (4) and the driving forces
. The PE of Fig-
ure 1 was used for the modeling. The initial conditions corresponded to the Brownian particles at rest at the left
minimum of the PE. Typical fission rates resulted from this modeling according to Equation (1) are shown in
Figure 2. Although the quasistationary rate
by some 20%, one can expect much larger dif-
ference considering Equations (10) and (11).
In order to figure out the reason for this surprisingly small difference let us transform Equation (10). The de-
rivative of the LDP reads:
I. I. Gontchar et al.
Figure 1. The schematic PE for 236U (curve
with dots) as a function of collective coordinate
along with its components: the smooth (solid
curve) and shell correction (dashed curve) parts.
Figure 2. The dynamical fission rates calcu-
lated for 236U according to Equation (1) versus
modeling time. The horizontal lines indicate
the quasistationary values of the rates.
∂⋅ ⋅ ⋅⋅⋅
Keeping in mind that
we convert Equation (10) to the form
turns out to be a small parameter not exceeding 0.1 in the region of the barrier at
MeV in particular because in the considered case
is in order of 2 MeV. Omitting in (13) the
terms obviously proportional to
we arrive at the formula
IIIS LS SL
is about 20 MeV, in Equation (14) the terms containing
are rather small. The last term is
is rather smooth (see Figure 1). We believe that this derivation explains relatively small
resulting in small difference between the corresponding fission rates in Figure 2.
We also have considered the case of a nucleus with the large shell correction at the ground state, namely the
double magic lead-208. The PE of Figure 3 was used for the modeling in this case. The absolute value of the
shell correction is now about 12 MeV. The way of modeling Figure 4. One sees that in this case
by factor of 2 whereas the conventional value of the intrinsic excitation energy at which the shell correction
ceases to manifest itself (and consequently
must converge) is about 60 MeV.
It is worthwhile to discuss and illustrate what does the smearing of the shell correction exactly mean. Within
the framework of the first recipe, it turns out that one can calculate the fission rate ignoring completely the shell
I. I. Gontchar et al.
correction in the PE (i.e. using
in Equation (7)) and in the LDP (i. e. using
in Equation (6)). The corresponding fission rate
must converge with
increases. This is seen in
Figure 5 where the quasistationary rates resulting from dynamical modeling are shown. In order to quantify the
we have calculated the fractional difference
. It has turned out
that the difference goes to zero reaching circa −10% at
MeV as it is known for the statistical rates.
Figure 3. The schematic PE for 208Pb (curve
with dots) as a function of collective coordi-
nate along with its components: the smooth
(solid curve) and shell correction (dashed
Figure 4. The dynamical fission rates calcu-
lated for 208Pb according to Equation (1) ver-
sus modeling time. The horizontal lines indi-
cate the quasistationary values of the rates.
Figure 5. The dynamical quasistationary fis-
sion rates calculated for 208Pb nucleus ac-
counting for the shell correction (squares) and
ignoring those (line with circles).
I. I. Gontchar et al.
Within the framework of the second recipe, the convergence is not reached because the reference point for the
excitation energy always include the potential energy with the shell correction (see Equation (6) of ).
To finalize, we have compared two approaches for calculating the driving force for the nuclear fission process
at low excitation energies when the shell effects are expected to be significant. We have found that in the case of
uranium-236 nucleus the quasistationary decay rates
resulting from these approaches are rather
close (the difference is about 20%). This is however just because for this nucleus the shell correction is small in
comparison with the typical energy
MeV controlling the smearing out the shell effects. For the lead-
208 nucleus with larger value of the shell correction, the difference between
reaches factor of 2.
This is significantly larger than the difference between the rates calculated within the frame work of the first ap-
proach with and without the shell correction. Since the first approach is based on the thermodynamical argu-
ments, we are inclined to make favor to it in comparison with the second one.
M. V. C. and E. G. D. are grateful to the Dmitry Zimin Foundation ‘Dynasty’ for financial support.
 Wagemans, C. (1991) The Nuclear Fission Process. CRC Press Inc.
 Hilscher, D. and Rossner, H. (1992) Dynamics of Nuclear Fission. Annals of Physics (France), 17, 471-552.
 Oganessian, Yu. (2013) Heaviest Nuclei. Nuclear Physics News, 23, 15-21.
 Itkis, I.M., et al. (2011) Fission and Quasifission Modes in Heavy-Ion-Induced Reactions Leading to the Formation of
Hs. Physical Review C, 83, 064613. http://dx.doi.org/10.1103/PhysRevC.83.064613
 Oganessian, Yu.Ts., et al. (2010) Synthesis of a New Element with Atomic Number Z = 117. Physical Review Letters,
104, 142502. http://dx.doi.org/10.1103/PhysRevLett.104.142502
 Pavlova, E.G. and Gontchar, I.I. (2012) Dissipative Statistical and Dynamical Fission Rates: Case of the Microcanoni-
cal Ensemble. Proceedings of the 4th International Conference of the Current Problems in Nuclear Physics and
Atomic Energy, Kyiv, 3-7 September 2012, 315-319.
 Gontchar, I.I., Fröbrich, P. and Pisc hasov, N.I. (1993) Consistent Dynamical and Statistical Description of Fission of
Hot Nuclei. Physical Review C, 47, 2228-2235. http://dx.doi.org/10.1103/PhysRevC.47.2228
 Gontchar, I.I. and Fröbrich, P. (1993) Nuclear Fission: Combining the Dynamical Langevin Equation with the Statisti-
cal Model. Nuclear Physics A, 551, 495-507. http://dx.doi.org/10.1016/0375-9474(93)90459-B
 Adeev, G.D., Karpov, A.V., Nadtochy, P.N. and Vanin, D.V. (2005) Multidimensional Stochastic Description of the
Excited Nuclei Fission Dynamics. Physics of Particles and Nuclei, 36, 733-820.
 Ignatyuk, A. V., Itkis, M.G., Okolovich, V.N., Smirenkin, G.N. and Tishin, A.S. (1975) Fission of Pre-Actinide Nuclei.
Excitation Functions for the (α, f) Reaction. Physics of Atomic Nuclei (Yadernaya Fizika), 21, 1185-1205.
 Aritomo, Y. and Oh ta, M. (2004) Dynamical Calculation for Fusion-Fission Probability in Superheavy Mass Region,
Where Mass Symmetric Fission Events Originate. Nuclear Physics A, 744, 3-14.
 Aritomo, Y. and Chiba, S. (2013) Fission Process of Nuclei at Low Excitation Energies with a Langevin Approach.
Physical Review C, 88, 044614. http://dx.doi.org/10.1103/PhysRevC.88.044614
 Brack, M., Damgaard, J., Jensen, A.S., Pauli , H.C., Strutinsky, V.M. and Wong, C.Y. (1972) Funny Hills: The Shell-
Correction Approach to Nuclear Sell Effects and Its Applications to the Fission Process. Reviews of Modern Physics,