“…The latter are very difficult to constrain using only ground state properties, thus the use of Landau response functions could be of help to define the correct interval of variation of such terms, thus avoiding instability regions [27,11]. However, the latter conditions are necessary but not sufficient to have a stable Skyrme pseudo-potential: as recently shown in [47,48], poorly-constrained coupling constants related to gradient terms can produce finite-size instabilities. These instabilities can be detected using the linear response formalism [49,50,51], with no extra computational cost and directly included into the optimization procedure [52].…”