Comment on the paper Eur. Phys. J. A (2019) 55:150
The conclusions of the study published as Eur. Phys. J. A (2019) 55:150 questioning the adequacy of the recently proposed Gogny D1M* interaction for finite-nuclei calculations using harmonic oscillator (HO) basis are revised. The existence of an instability in finite nuclei for D1M* when coordinate-space methods are used to solve the HF equations (as shown in Eur. Phys. J. A(2019) 55:150) is independently confirmed using a computer code based on a quasi-local approximation (QLA) to the HF energy density with finite-range forces. We confirm that the most affected quantity in the coordinate-space calculation is the spatial density at the origin. Our study reveals that some findings concerning these instabilities are not easy to reconcile with the arguments used in Eur. Phys. J. A (2019) 55:150. For instance, some nuclei such as $^4$He and $^{40}$Ca, which diverge in HF mesh-point calculations performed with D1M*, become perfectly stable when Coulomb force is switched off. We have also found instabilities in some nuclei when the D1M interaction is used. Finally, a connection between the occupancy of $s$-orbitals near the Fermi level and the appearance of instabilities is observed. Several convergence and stability studies are performed with HO basis of different sizes and oscillator parameters to demonstrate the robustness of the D1M* results for finite nuclei when the HO basis is used.