Numerical Modeling of Scission Neutrons Emitted During Nuclear Fission

We have considered a problem related to nuclear fission, namely the emission of scission neutrons. To study this process we have used the sudden approximation and also a dynamical approach. Numerically, this implies the solution of the Schroedinger equation, both stationary and time-dependent. To describe deformed nuclear shapes, we have used a model based on cylindrical coordinates. The Hamiltonian is discretized by using finite difference approximations of the derivatives. The eigenstates of the two-dimensional stationary Schroedinger equation are obtained as solutions of an algebraic eigenvalue problem, using an appropriate procedure (Arnoldi). In the sudden approach, the physical characteristics of the phenomenon are obtained from the wavefunctions before and after scission. In the more realistic evolutional (diabatic) model, the initial wavefunctions are propagated in time, solving the temporal Schroedinger equation by a Crank-Nicolson scheme