Hartree-Fock-Bogoliubov study of quantum shell effects on the path to fission in $^{180}$Hg, $^{236}$U and $^{256}$Fm

Quantum shell effects stabilising fission fragments with various shapes have been invoked as a factor determining the distribution of nucleons between the fragments at scission. Shell effects also induce asymmetric shapes in the nucleus on its way to fission well before the final fragments are (pre)formed. These shell effects are studied in fission of $^{180}$Hg, $^{236}$U and $^{256}$Fm with constrained Hartree-Fock-Bogoliubov calculations using the D1S parametrisation of the Gogny interaction. Strutinsky shell energy correction and single-particle energy level density near the Fermi surface are computed. Several neutron and proton shell effects are identified as drivers towards asymmetric fission. Shell effects are also used to identify the preformation of the fragments in the later stage of fission.


Introduction
Quantum shell effects are usually associated with a large energy gap of a few MeV in the single-particle energy spectrum near the Fermi level.Such gaps provide additional binding, and thus more stability, to the system which can be spherical (as in the so-called "magic" nuclei) or deformed.In the fission process, the nucleus experiences a competition between Coulomb repulsion and nuclear attraction which normally deforms the system up to the formation of two fragments with similar masses.On its way to fission, however, the system may encounter shell effects that drive it out of the symmetric fission path, i.e., favouring pear shapes to the nucleus, leading to mass asymmetric fragments [1,2].The system thus follows one or several fission valleys in a "potential energy surface" (PES) representing the potential energy as a function of collective coordinates used to characterise the deformation of the system.Its dynamics is then guided by the topology of the PES along with the associated inertia.
Quantum shells stabilising fission fragments with spherical [3,4], cigar [5], and pear [6] shapes have been invoked as a factor determining the distribution of nucleons between the fragments at scission.These shells, however, are found only in the later stage of the fission process when pre-fragments are already formed, while other shell effects are expected in the fisa Present address: CEA, DES, Cadarache, F-13108 Saint Paul les Durance, France b cedric.simenel@anu.edu.ausioning nucleus.For instance, the decrease in energy of neutron orbitals has been invoked as a mechanism inducing asymmetric shapes near the second saddle point of actinides [7] as well as in neutron deficient mercury region [8].In addition, a third hyperdeformed minimum, produced by shell effects, is usually found in the fission path of actinides [9].
Shell energy corrections favouring specific elongations and asymmetries on the way to fission can be identified with the Strutinsky method [10,11].Naturally, the shell-energy correction and single-particle level density are intrinsically related as such corrections are the result of energy gaps appearing in the single-particle energy spectra [12].Deformed shell effects on the way to fission can then be investigated through their effect on shell energy correction, or directly from single-particle energy level density near the Fermi surface.
In this work, shell effects are studied in fission of 180 Hg, 236 U and 256 Fm.These nuclei are located in different mass regions of interest for nuclear fission studies [1,2] as they are expected to exhibit different shell effects in the final fragments [6,13,14,15,16].PES are obtained with the constrained Hartree-Fock-Bogoliubov (HFB) method.In addition to the total potential energy of the deformed system used to construct the PES, the HFB calculations give access to the single-particle energy levels that are used to compute the shell-energy corrections from the Strutinsky method as well as the single-particle energy level density near the Fermi surface.The method and numerical details are provided in section 2. Resulting PES and shell effects are discussed in section 3. Conclusions are drawn in section 4.
2 Theoretical method and numerical details 2.1 Constrained Hartree-Fock-Bogoliubov calculations Quantities of interest are studied within the HFB approximation under constraints using the D1S parametrisation [17] of the Gogny nuclear interaction [18].Details of the procedure can be found in Ref. [19].Axial (z−axis), time reversal and simplex symmetries are preserved in this study.As a result of the time reversal symmetry, only systems with an even number of protons and neutrons are computed.
A good representation of PES is usually obtained by imposing quadrupole Q 20 (cigar shape) and octupole Q 30 (pear shape) moments, while leaving free the higher multipole moments in an energy minimisation procedure [20,21].PES are generated by constraining the quadrupole moment and the octupole moment where ρ(r) is the density of nucleons.Higher multipole moments are not constrained in the PES calculations.The PES are built by minimising the HFB energy for each (Q 20 , Q 30 ) deformation with a 2 b quadrupole moment step and a 2 b 3/2 octupole moment step (1 barn ≡ 1 b = 10 −28 m 2 ).The resulting constrained HFB states have no internal excitation.Thus, the present calculations are only valid at low energy.Indeed, pairing energy and shell effects are affected by excitation energy and thus must be accounted for in fission studies of warm nuclei [22,23].The Fock space is spanned in the harmonic oscillator basis whose number of major shells is 14 for fissioning nuclei and fragments (see Appendix A).The HFB calculations are performed with Lagrange parameters to constrain the number of protons and neutrons in the system.These Lagrange parameters are chemical potentials, or equivalently Fermi energies marking the transition between occupied and unoccupied states.The HFB generalised Hamiltonian is represented by a matrix containing the usual Hartree-Fock single-particle Hamiltonian h in its diagonal and the pairing field in the off-diagonal terms.In this work, the single-particle energies are defined as the eigenvalues of h.The stability of the results is studied in Appendix A in which the size of the harmonic oscillator basis has been raised to N = 15 for 256 Fm.Its barrier height has been calculated up to N = 24.
Discontinuities [20,24,25,26] are often associated with a sudden change of shape in the system.They can be identified through a measure of the integrated difference in the density distributions between two neighbouring points in the PES [24], through an evaluation of overlaps between neighbouring states of the PES [27], or through the analysis of unconstrained multipole moments, among which the triaxial quadrupole moment Q 22 and the hexadecapole moment Q 40 are expected to be the most important [26].We have used the latter method with Q 40 only as for an axial code, Q 22 is constrained to be zero.Smoothing these discontinuities is beyond the scope of this work and could be achieved using methods proposed in Refs.[24,28].
It is also possible to find more than one local minimum for a given set of constraints in the construction of the PES.To avoid this issue, one can perform a retropropagation scheme and keep the lowest minimum [25].Another approach, which is the one we adopted, is to keep the system in the fission valley as long as possible (see discussion in [26]).
The resulting PES of 180 Hg, 236 U and 256 Fm are shown in Fig. 1(a,b,c) and are discussed in Sec. 3.

Definitions of the fission path, of scission and of the pre-fragments
In principle the path to fission is the one that leaves the action of the system stationary.In this work, what is called "fission path" is obtained by assuming that the system follows a minimum energy path on its way to fission.Although this approximation is expected to break down when non-adiabatic effects become significant (e.g., near scission) as well as near discontinuities, it is sufficient for the purpose of this work that is to investigate the role of shell effects.Note that, although collective inertia is calculated (see Appendix C), it is not used in the path calculation.
Several definitions of scission have been used by different groups (see discussion in [21]).Here, it is assumed that scission is achieved when the spatial density in the neck becomes lower than 0.08 fm −3 (approximatively half the nuclear saturation density).Calculations along the fission paths and within the PES are stopped at the last point before reaching scission when increasing the elongation.This defines the scission point and the pre-scission line, respectively.Although this simple definition of scission could lead to ambiguities in characterising the fragments (e.g., due to their entanglement), it is sufficient in the present work that is devoted to shell effects within the PES up to the scission region, as close as possible to the nascent fragments.
For a given HFB configuration of the fissioning nucleus, the spatial density minimum in the neck along the deformation axis is used to identify left and right pre-fragments.Integrating the spatial density on the left and on the right enables us to get pre-fragment mass and charge numbers and their individual multipole moment deformations for λ ∈ [2,6].Y λ 0 are the spherical harmonics and θ and φ are the polar and azimuthal angles.Particle numbers and these multipole moments are then used as constraints to calculate pre-fragments alone.

Shell energy corrections and single-particle level density near the Fermi surface
Details of the shell energy corrections method of Strutinsky [10,11] can be found in [29].The width over which the smoothing of the single particle level scheme is performed is chosen as γ = 41A −1/3 , A being the number of neutrons or protons.
The density of energy levels in the system determines important quantities that affect fission mechanisms, including energy dissipation [23], fragment excitation [30], and shell effects [12].In the case of shell effects, however, only the level density near the Fermi level is relevant.To investigate the evolution of shell effects along the fission path, we introduce an energy level density that counts the number of neutron or proton single-particle states in an energy window ∆ E centred in ε n,p 0 that are defined as halfway between the two levels surrounding the neutron and proton Fermi energy, respectively.∆ E should be of the order of typical shell energy gaps so that η n,p produces small values in the presence of such a gap.As the nucleon energy spectra are less compressed in lighter nuclei, we use ∆ E = 3 MeV in 180 Hg and 2.5 MeV in the other nuclei.To avoid rapid changes in η n,p when a nucleon energy level enters or leaves the energy window, we also weight the contribution of the levels by a function f of the nucleon energy e i that is equal to 1 at ε n,p 0 and linearly decreases to 0 at the edges of the energy window.The effect of ∆ E and of the smoothing function are studied in Appendix B.
Figure 2 shows how η n and η p are computed for a state of 180 Hg on its fission path [represented by a star in Figs 1(a,d,g)].
In particular, we see that a large energy gap near the neutron Fermi level leads to a small value of η n .The level densities η n and η p can then be used as an indication of the presence of neutron and proton shell effects, respectively.In turn, these shell effects are expected to produce a shell correction energy (as computed with the Strutinsky method) and also affect the inertia (and thus the dynamics) of the system as it evolves on the PES [31].
Figures 1(a-c) show the PES for 180 Hg, 236 U and 256 Fm.In each case the fission is mass asymmetric as indicated by the non-zero values of Q 30 along the one-dimensional fission path obtained by following the path of local minimum energy in the potential energy surface until scission.The predicted fragment mass asymmetries at the end of the fission path are in good agreement with the centroids of the experimental fragment mass distributions [32,33].No discontinuity is observed along the fission path except a small one at Q 20 100 b in 180 Hg and a larger one at Q 20 90 b in 256 Fm.In both cases, they correspond to "jumps" in the octupole moment Q 30 as seen from the PES in Figs.1(a) and 1(c).We checked that, for the main discontinuity in 256 Fm, Q 40 evolves smoothly with Q 30 [See Fig. 7(c) in Appendix C].As inertia plays an important role in cold fission, the perturbative masses have been calculated using the adiabatic time-dependent Hartree-Fock-Bogoliubov (ATDHFB) and generator coordinate method (GCM) approximations.They are presented in Appendix C.
The existence of a third minimum in fission path of actinides is currently debated [34,35,36,37].Here, no third minimum has been found on the 236 U asymmetric path.However, 256 Fm experiences a strong first barrier (about 10.5 MeV) followed by a shallow second minimum at ∼ 52 b and a shallow third minimum at ∼ 72 b.The second and third minima are separated by a ∼ 1 MeV barrier.
Figure 1 also shows the level densities η n,p for neutrons (d-f) and protons (g-i) as well as Strutinsky shell energy correction from the neutron (j-l) and proton (m-o) HFB single particle energy spectra for the same three nuclei and ranges of deformations as in the PES.Both methods indicate that, in each system, when mass symmetry is broken the fission paths follow successive regions associated with shell energy correction and, similarly, small neutron and/or proton level densities.We see that several shell effects are at play along the fission path and that both proton and neutron shells may take turn in driving the system towards mass asymmetric fission.
Note also that in the Hartree-Fock-Bogoliubov theory used in this work, static pairing correlations are fully accounted for and as the finite range Gogny interaction is used, no cut off is required for the pairing energy.Thus the results on the level density function and the Strutinski shell energy correction based on the PES analysis both include a full treatment of the static pairing correlations.Ambiguities may arise when comparing pre-scission configurations with individual studies of fragments as different Lagrange parameters have to be used to constrain the number of particles in both cases.However, our study focuses essentially on the role of shell effects in pre-scission configurations and thus we do not expect our conclusions to be affected by such issues.
The pairing interaction, which occurs essentially near the Fermi level, is attractive and increases with single-particle energy level density.It is then in competition with shell effects that favour low level densities.The fact that the fission paths follow regions of low level density shows that the dominant effect comes from shell effects rather than pairing energy.
Comparing figures 1(d-i) with figures 1(j-o), we see that the variations of the energy correction display structures over the 2D surfaces that are sometimes similar to the ones in the near Fermi level density surfaces but also differ in some cases.Focussing on the 1D fission paths, both methods predict a dominance of proton shell effects at Q 20 > 70 b in 236 U [Figs.1(h,n)] and at Q 20 > 110 b in 256 Fm [Figs.1(i,o)].The level density method, however, predicts a region of low level density from neutron shell effects at Q 20 ∼ 100 b and Q 30 ∼ 20 b 3/2 [Fig.1(f)] that could be responsible for the early asymmetry in 256 Fm, though the latter is not associated with significant shell correction energy [Fig.1(l)].Similarly, the level density method predicts a neutron gap in 236 U at 100 < Q 20 < 130 b [Fig.1(e)] that is not observed in the shell correction energy [Fig.1(k)].As for 180 Hg, while both methods agree on the dominance of neutron shell effects in the later stage of fission [Fig.1(d,j)], only the level density exhibits structures that are clear enough so that we can interpret the discontinuity at Q 20 ∼ 100 b as due to a change from proton to neutron shell effects.
In the following, we base our analysis of shell effects on the level density plots in Figs.1(d-i).In 180 Hg, both proton and neutron shells are present when the system first acquires a non-zero octupole moment Q 30 [Figs.1(b,g In the case of 236 U, this confirms that, on the one hand, neutron shell effects allow the nucleus to leave the symmetric path to find a lower (asymmetric) second barrier [7,12], and, on the other hand, proton shell effects are present near scission [6].We also see that these two mechanisms do not provide a complete picture.Indeed, in 236 U, we observe several additional proton and neutron shell effects on the asymmetric fission path [Figs.1(e,h)], showing that the connection between the first and last shell effects is not trivial.This observation is in contradiction with an earlier interpretation (in the case of 232 Th) that a single shell effect due to a partially formed spherical Sn-like fragment is responsible for lowering the level density from the asymmetric saddle to scission [38].
Note that shell effects are not limited to the fission path (which is only the most probable one) and may affect other regions of the PES which could influence fragment mass and charge distributions.For instance, by comparing Figs.1(b), 1(e) and 1(h) we see that in 236 U, proton shell effects produce another fission valley which is symmetric (Q 30 0), while the extension of the PES around Q 20 ∼ 200 − 220 b and Q 30 ∼ 30 − 50 b 3/2 is clearly due to neutron shell effects.In addition, a shallow symmetric valley in 256 Fm [Fig.1(c)], produced by spherical shell effects in the fragments leading to low proton level density [Fig.1(i)], is observed (see also Fig. 3).
In principle, both the level density and the shell energy correction methods are able to predict at which stage the final fragments are pre-formed in the system, while still interacting via the strong force through the neck.Fragments can be considered as being pre-formed when the nucleons occupy the states that are necessary to build their final microscopic structure, even though the shapes of the pre-formed fragments may not be the final ones.If a shell energy gap, characterised by a low level density at the Fermi level, is present and extends up to scission, we can then consider the fragments to be pre-formed in this region of the PES.Indeed, the presence of the gap means, in principle, that no state is crossing the Fermi level, thus all the "required" states for the final fragments are already occupied [39] (at least for the specie of nucleons -neutrons or protons -that exhibits the gap).For example, the neutron shell gap that constrains the fission path of 180 Hg near scission appears at Q 20 100 b, at which point the fragments can be considered to be pre-formed.
The methods of identifying shell effects from shell energy corrections and from the level density provide information on the entire system, such as which of proton and/or neutron shells are the main drivers towards fission.Once the fragments are pre-formed, however, they have their own microscopic structure which we can use to investigate which of the fragments exhibits shell effects fixing the final mass asymmetry [6,15,16].A similar analysis is presented in figure 3. The fragments are identified by integrating the proton and neutron densities on each side of the neck, at the most probable scission configuration (i.e., at the end of the fission path).This only gives an average number of protons and neutrons.However, the shell structure of a fragment is expected to show only small variations by adding or removing one nucleon.This allows us to determine the proton and neutron levels of the fragments from microscopic calculations of nuclei with even numbers of protons and neutrons.Nevertheless, knowing the numbers of protons and neutrons is not sufficient as the shell structure is affected by the deformation of the nuclei.It is therefore important to impose deformations in the calculations of the individual nuclei which are close to those of the fragments observed in the total system at scission.The fragments are then computed individually (without interaction between them), with constraints on their multipole moments which are optimised to reproduce the scission configuration.This allows us to compute the singleparticle spectra for each fragment.In figure 3, this is done by adjusting the multipole moments Q λ 0 of the simulated fragments with λ = 2 − 6.
Figure 3 shows that the low level densities in compound nuclei identified at scission are associated with shell gaps at the Fermi level in the fragments.A neutron shell gap at N = 46 is present in the light fragment ( 82 Kr) of 180 Hg, confirming that neutron shells drive this system near scission.A similar situation occurs with proton shells in both fragments of 236 U and 256 Fm with, in addition, the presence of neutron shell gaps in the heavy fragments which could also impact the fission process.The fact that shell effects can be found in both fragments does not necessarily mean, however, that both fragments have an equal influence in the final asymmetric configuration.Indeed, shell effects in the fission fragments could either be a cause or a consequence of the most likely state at scission.As mentioned earlier, shell effects provide an additional stability to the system, thus favouring a specific shape in which such shell effects can be found.At scission, each fragment can have a different deformation, e.g., one fragment can be compact and the other one elongated, as in the systems we have studied [see Fig. 3(a)].A leading shell effect could stabilise the shape of one fragment whereas the complementary fragment minimises the energy while preserving the compound nucleus elongation and asymmetry.Under this overall deformation constraint it is then possible that one fragment is responsible for the configuration at scission, thus driving the final asymmetry, while the other fragment "just" finds a deformed shell effect.
In actinide nuclei such as 236 U and 256 Fm, experimental [33,40,41,42,43] systematic studies indicate an influence of proton shell effects on the heavy fission fragment properties.One possible origin of these quantum shell effects is octupole correlations favouring pear shapes in the fragments [6].In lighter nuclei in the lead region (such as 180 Hg), the origin of the shell effects driving the final asymmetry at scission is still debated [14,15,16,44,45].Future studies of the role of shell effects in the fission paths should provide a deeper insight into this region.More generally, signatures indicating which shell effect (if any) is the main driver could be searched for in the evolution of the system in the fragment pre-formation region.Indeed, the same technique presented in Fig. 3 to analyse the fragments at scission can be used to investigate less elongated configurations, prior to scission.

Conclusions
The role of shell effects in fission has been investigated with the shell correction method of Strutinsky together with the singleparticle energy level density near the Fermi surface.Both methods show that the shape evolution of the fissioning nucleus is determined by several quantum shells that are driving the nucleus towards mass asymmetric fission well before the final fragments are (pre)formed.In the three nuclides we investigated ( 180 Hg, 236 U and 256 Fm), at least 1 to 2 neutron and 2 to 3 proton shell effects can be identified on the way to fission.Several (but not all) of these shell effects lead to significant shell energy corrections.The appearance of shell effects attributed to the fragments can also be used to sign their pre-formation.These pre-fragments at the end of the fission path are identified as 180 Hg→ 82 Kr+ 98 Ru, 236 U→ 96 Sr+ 140 Xe, and 256 Fm→ 112 Ru+ 144 Ba.The analysis of their single-particle structures indicates that several deformed shell gaps are present in these pre-fragments: N = 46, Z = 38, 44 with elongated shapes, and Z = 52, 56, N = 84, 88 with compact octupole shapes.Although our focus were on asymmetric fission, similar methods could also be used to investigate symmetric modes, e.g., 256 Fm→ 128 Sn+ 128 Sn that exhibits spherical shell gaps in the fragments.Furthermore, it would be interesting to study the microscopic origin of the new compact symmetric fission mode in light thorium isotopes recently observed experimentally [46].
Recent studies with time-dependent extensions to the meanfield approach that was used here (see [47] for a review) have also shown the importance of shell effects in quasi-fission reaction [48,49,50].Combined with shell effect identification methods as used in this work, such time-dependent approaches could be used to investigate the interplay between energy dissipation and shell effects that would impact the energy released in nuclear fission.A complete theoretical description of the fission process remains an important challenge which requires further development of fully microscopic approaches to tackle the complexity of the fission mechanism [51], as illustrated by the present observation that many proton and neutron shells are at play in each nuclides.
APPENDIX A: Convergence with the number of harmonic oscillator shells.
To test the convergence of our calculations with the number of harmonic oscillator shells (N = 14 in this work), the PES and level densities have been computed for N = 15 shells in the case of 256 Fm.Results are presented in Fig. 4(a) for the PES and Figs.4(c,f) for level densities.The results shown in Fig. 1(c,f,i) with N = 14 are reported in Figs.4(b,d,g) for comparison.The differences between level densities computed with N = 14 and N = 15 are small as shown in Figs.4(e,h), which is a good indicator of the convergence of our calculations with the number of oscillator shells.In particular, they do not exceed 0.5 and usually remain smaller than 0.2 except in small areas of the PES.

APPENDIX B: Comparison between several definitions of the density level near the Fermi energy
The energy level density near the Fermi level is defined as η n,p = ∑ i f (e i − ε n,p 0 ).The effect of the smoothing function f (defined in Fig. 2) and of the energy window ∆ E in which it is non-zero are illustrated in Figs. 5 and 6 for neutrons and protons, respectively.The smoothing is removed in Figs.5(ac) and Figs.6(a-c), while the energy window is doubled in Figs.5(g-i) and Figs.6(g-i).Results with the standard smoothing defined in the manuscript are reproduced in Figs.5(d-f) and Figs.6(d-f) for comparison.The same shell effects are visible with these three definitions, and thus our conclusions are not sensitive to (reasonable) variations of the energy window and the smoothing function.

APPENDIX C: Properties of the 1D fission paths
Characteristics of the 1D asymmetric paths are presented in Fig. 7.The evolution of the hexadecapole moment Q 40 along the fission path is plotted for (a) 180 Hg, (b) 236 U, and (c) 256 Fm as a function of the quadrupole moment Q 20 .No discontinuity is observed for this degree of freedom except a small one at Q 20 100 b in 180 Hg and a larger one at Q 20 90 b in 256 Fm.In both cases, they correspond to "jumps" in the octupole moment Q 30 as seen from the PES in Figs.1(a,c).The inset in Fig. 7(c) shows that, for the main discontinuity in 256 Fm, Q 40 evolves smoothly with Q 30 , indicating that no higher multipole moments are necessary to explain this discontinuity.The evolution of the neutron and proton level densities η n,p along the 1D paths extracted from Figs. 1(d-i) are plotted in Figs.7(d,e,f) for the same systems.
The perturbative mass parameters accounting for the inertia of the systems are shown in Figs.7(g,h,i).They have been computed via two standard methods: the adiabatic time-dependent Hartree-Fock-Bogoliubov (ATDHFB) approach and the generator coordinate method (GCM) (See Ref. [31] for details).Level crossing at the Fermi surface are known to produce peaks in the mass parameters, which is why we see more structure in the inertia evolution at small Q 20 .Naturally, the level density near the Fermi level is lowered by such crossing, while in regions of the PES affected by strong shell effects, no such crossing is expected.As a result, we expect correlations between mass parameters and level density at the Fermi surface.When shell effects are strong, such as in the 236 U neutrons at Q 20 ∼ 110 − 130 b in the fission path, we indeed see a peak in the mass parameter on each side of this region.

Fig. 1 .Fig. 2 .
Fig. 1. (a-c) The potential energy is plotted as a function of the quadrupole Q 20 and octupole Q 30 moments in 180 Hg, 236 U and 256 Fm.Fission paths are represented by thick green lines.Density contours along the fission path are also shown.The most probable fragments produced at scission are indicated by an arrow.The dotted segment in (c,f,i,l,o) indicates a discontinuity between the symmetric and asymmetric valleys (see, e.g.,[24,25,26]) when both valleys have the same energy (in this case, the ridge between both valleys is only ∼ 0.5 MeV).The (d-i) panels show the neutron (η n ; (d-f)) and proton (η p ; (g-i)) level densities near the Fermi energy, respectively.Shell effects are stronger for darker colours.The star (a,d,g) indicates the configuration used to produce the single particle level scheme in figure2.Finally, neutron (j-l) and proton (m-o) Strutinsky shell energy corrections are plotted from the HFB single particle energy spectra.
)].In 236 U, the first asymmetry is induced by neutrons [Fig.1(e)], while in 256 Fm it is due to protons [Fig.1(i)].Indeed, in the latter, the low η p at Q 20 60 − 90 b lowers the energy ridge between the symmetric path and a nascent asymmetric one.Near scission, the final asymmetry between the fragments of 180 Hg is determined by neutron shells [Fig.1(d)], while in 236 U and 256 Fm, it is dominated by proton ones [Figs.1(h,i)].Strong shell effects are also observed at intermediate configurations along the fission path, such as neutron shells at Q 20 ∼ 120 b in 236 U [Fig. 1(e)] and 100 b in 256 Fm [Fig.1(f)].

Fig. 3 .
Fig. 3. (a) The density distributions of the nuclei on their fission path just before scission are shown in the top panels.Below each panel are shown the density distributions for the "closest" fragments with even numbers of protons and neutrons.(b) The resulting proton and neutron energy levels e i and occupation numbers g n,p i (non-transparent part of horizontal lines) near the Fermi energy are shown for each fragment.The numbers in the energy gaps correspond to the number of particles that can occupy the single-particle states below the gap.
Further convergence tests up to N = 16 have been made at specific configurations for 256 Fm, namely the ground state and the barrier height, giving E GS = −1891.998MeV (N = 14), −1892.412MeV (N = 15), and -1892.710MeV (N = 16) for ground-state energy and B = 11.220MeV (N = 14), 11.248 MeV (N = 15) and 11.265 MeV (N = 16) for the barrier B = E saddle − E GS .When increasing further the number of major shells these numbers are found to be stable with: B = 11.236MeV (N = 20) and B = 11.175MeV (N = 24).
The large peak observed at ∼ 90 b in 256 Fm [Fig.7(i)] and the smaller one at ∼ 100 b in 180 Hg [Fig.7(a)] are associated with the "jumps" in the PES [Figs.1(c,a)] and in Q 40 [Figs.7(c,a)].

Fig. 5 .
Fig. 5.The neutron level density near the Fermi energy is computed without smoothing function within the energy windows ∆ E = 3 MeV in (a) 180 Hg and ∆ E = 2.5 MeV in (b) 236 U and (c) 256 Fm.Panels (d-f) are the same as in Fig. 1.Level densities in panels (g-i) are obtained with the smoothing function applied on energy windows increased by a factor 2.

Fig. 7 .
Fig. 7. Hexadecapole moments (a-c), smoothed level densities (d-f) and perturbative masses (g-i) as a function of quadrupole moment along the 1D asymmetric fission paths of the three nuclei.The inset of panel (c) shows the evolution of the hexadecapole moment as a function of the octupole moment in the region of the discontinuity.See Appendix C for details.