jimmy-wong/langevin

Source file structure

src/cubic_spline.cpp

calculate the splined value in the gap of the grids.

src/inertia.cpp

calculate the inertia tensor of the compound nuclei.

src/dissipative.cpp

calculate the dissipative tensor of the compound nuclei.

src/random.cpp

generate the random force on the nuclei.

src/shape.cpp

class shape: calculate the surface normal speed and other shape related terms. A_Block is the function of $$\frac{\partial}{\partial q_l}\int_z^{z_max}\rho(\zeta',q)^2d\zete'$$

src/store.cpp

store the current PES.

src/Ronge_kutta.cpp

solve the Langevin equation using Ronge-Kutta method.

src/langevin.cpp

main program

Program illustration

This program is used to simulate the dynamics of the fission process on the potential energy surface (PES) in the framework of Langevin equation. Currently, we only consider the temperature independent PES because the excitation energy of the compound nuclei is very low. The temperature of low excitation energy is at most 1MeV, which will not bring too much correction to the PES. Also we perform the temperature independent inertia and dissipative tensor on this surface.