We describe the software package FELIX that solves the equations of the time-dependent generator coordinate method (TDGCM) in N-dimensions (N $geq$ 1) under the Gaussian overlap approximation. The numerical resolution is based on the Galerkin finite
element discretization of the collective space and the Crank-Nicolson scheme for time integration. The TDGCM solver is implemented entirely in C++. Several additional tools written in C++, Python or bash scripting language are also included for convenience. In this paper, the solver is tested with a series of benchmarks calculations. We also demonstrate the ability of our code to handle a realistic calculation of fission dynamics.
The constrained Hartree-Fock-Bogoliubov method is used with the Gogny interaction D1S to calculate potential energy surfaces of fissioning nuclei ${}^{226}$Th and ${}^{256,258,260}$Fm up to very large deformations. The constraints employed are the ma
ss quadrupole and octupole moments. In this subspace of collective coordinates, many scission configurations are identified ranging from symmetric to highly asymmetric fragmentations. Corresponding fragment properties at scission are derived yielding fragment deformations, deformation energies, energy partitioning, neutron binding energies at scission, neutron multiplicities, charge polarization and total fragment kinetic energies.
