“…To simplify the numerical calculations and following [50,51], in the present work we use central correlation functions g(i,j ) independent on spin and isospin. Moreover, it has been shown [52,53] that this central correlation functions, in which are included the main contributions of the 1 S 0 and 3 S 1 channels, are weakly dependent on the density, and can be approximated [50,52,53] by a Heaviside step function θ (r ij − r c ), with r c = 0.6 fm in all the considered density range. Note that the average procedure has to be performed for each term involved in the cyclic permutation in Eq.…”