We explore the density variation of the correlation coefficient of the key parameters of the nuclear equation of state (EoS) with the bulk and crustal properties of neutron stars. The analysis was performed using two diverse sets of nuclear effective interaction theories based on nonrelativistic Skyrme-Hartree Fock model and relativistic mean field model. We find that the commonly studied EoS parameters, namely the isoscalar incompressibility of symmetric nuclear matter K (ρ) and the isovector slope of symmetry energy L (ρ), reveal consistently maximum correlation with the radius, tidal deformability, and moment of inertia all around twice the saturation density. We find even more tighter and robust correlations beyond the saturation density for constructed parameter η= KL²^1/3 allowing the possibility to impose stringent constraints on high-density K (ρ) and L (ρ). Extensive correlation analysis of the EoS parameters with the radius and tidal deformability bounds from the gravitational wave events and recent pulsar observations allow us to provide reliable constraints on the central values of K (ρ₀) 240 MeV and L (ρ₀) 48 MeV at saturation density and K (1. 6ρ₀) 332^+88-₅₀ MeV and L (1. 6ρ₀) 122^+26-₁₈ MeV at 1. 6 times the saturation density. The crust-core transition density and the crustal fraction of moment of inertia are shown to correlate moderately with L (ρ) and η (ρ) near the subsaturation density.
Alam et al. (Mon,) studied this question.