The nuclear matter parameters define the nuclear equation of state (EoS), they appear as coefficients of expansion around the saturation density of symmetric and asymmetric nuclear matter. We review their correlations with several properties of finite nuclei and of neutron stars within mean-field frameworks. The lower order nuclear matter parameters such as the binding energy per nucleon, incompressibility and the symmetry energy coefficients are found to be constrained in narrow limits through their strong ties with selective properties of finite nuclei. From the correlations of nuclear matter parameters with neutron star observables, we further review how precision knowledge of the radii and tidal deformability of neutron stars in the mass range $1 - 2 M_\odot$ may help cast them in narrower bounds. The higher order parameters such as the density slope and the curvature of the symmetry energy or the skewness of the symmetric nuclear matter EoS are, however, plagued with larger uncertainty. From inter-correlation of these higher order nuclear matter parameters with lower order ones, we explore how they can be brought to more harmonious bounds.

We search for possible correlations between neutron star observables and thermodynamic quantities that characterize high density nuclear matter. We generate a set of model-independent equations of state describing stellar matter from a Taylor expansion around saturation density. Each equation of state which is a functional of the nuclear matter parameters is thermodynamically consistent, causal and compatible with astrophysical observations. We find that the neutron star tidal deformability and radius are strongly correlated with the pressure, the energy density and the sound velocity at different densities. Similar correlations are also exhibited by a large set of mean-field models based on non-relativistic and relativistic nuclear energy density functionals. These model independent correlations can be employed to constrain the equation of state at different densities above saturation from measurements of NS properties with multi-messenger observations. In particular, precise constraints on the radius of PSR J0030+0451 thanks to NICER observations would allow to better infer the properties of matter around two times the nuclear saturation density.

We confront admixture of dark matter inside neutron star using gravitational wave constraints coming from binary neutron star merger. We consider a relativistic mean field model including $\sigma-\omega-\rho$ meson interaction with NL3 parameterization. We study fermionic dark matter interacting with nucleonic matter via Higgs portal mechanism. We show that admixture of dark matter inside the neutron star soften the equation state and lower the value of tidal deformability. Gravitational wave GW170817 observation puts an upper bound on tidal deformability of a binary neutron star with low spin prior at 90\% confidence level, which disfavors stiff equation of state such as Walecka model with NL3 parameterization. However, we show that Walecka model with NL3 parameterization with a fermionic dark matter component satisfy the tidal deformability bound coming from the GW170817 observation.

Constraints set on key parameters of the nuclear matter equation of state (EoS) by the values of the tidal deformability, inferred from GW170817, are examined by using a diverse set of relativistic and non-relativistic mean field models. These models are consistent with bulk properties of finite nuclei as well as with the observed lower bound on the maximum mass of neutron star $\sim 2 ~ {\rm M}_\odot$. The tidal deformability shows a strong correlation with specific linear combinations of the isoscalar and isovector nuclear matter parameters associated with the EoS. Such correlations suggest that a precise value of the tidal deformability can put tight bounds on several EoS parameters, in particular, on the slope of the incompressibility and the curvature of the symmetry energy. The tidal deformability obtained from the GW170817 and its UV/optical/infrared counterpart sets the radius of a canonical $1.4~ {\rm M}_{\odot}$ neutron star to be $11.82\leqslant R_{1.4}\leqslant13.72$ km.