where \(R\)  is the neutron star radius and \(\lambda =2R^5k_2/3G\)  is the other tidal parameter that we use in our study, the so-called tidal deformability. The tidal Love number \(k_2\)  is given by [1]}, [2]}
\(k_2&=&\frac{8\beta ^5}{5}\left(1-2\beta \right)^2\left[2-y_R+(y_R-1)2\beta \right] \times \left[\frac{}{} 2\beta \left(6 -3y_R+3\beta (5y_R-8)\right) \right. \nonumber \\&+& 4\beta ^3 \left. \left(13-11y_R+\beta (3y_R-2)+2\beta ^2(1+y_R)\right)\frac{}{} \right.\nonumber \\&+& \left. 3\left(1-2\beta \right)^2\left[2-y_R+2\beta (y_R-1)\right] {\rm ln}\left(1-2\beta \right)\right]^{-1},\) 
