“…( 35), we use m = 1.08 according to the HNC approximation [38], the soft MSA [22,39,40] signifies the mean spherical approximation (MSA) [1,41] for soft spheres, and the smeared interaction potential is identified with the DCF in the Bessel smoothed model [35] that supposes a smeared charge distribution expressed by the modified Bessel function K 1 (see Appendix A.2 for the details). [14,[20][21][22][23][24] In the real-space representation of the electrostatic DCF given by eq. ( 34), we have…”