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.

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.

We carry out the study for finite nuclei, infinite nuclear matter and neutron star properties with the newly developed relativistic force named as the Institute Of Physics Bhubaneswar-I(IOPB-I). Using this force, we calculate the binding energies, charge radii and neutron skin thickness for some selected nuclei. From the ground state properties of superheavy element i.e. Z=120, it is noticed that considerable shell gaps appear at neutron numbers N=172, 184 and 198, showing the magicity of these numbers. The low density behavior of the equation of state for pure neutron matter is compatible with other microscopic models. Along with the nuclear symmetry energy, its slope and curvature parameters at the saturation density are consistent with those extracted from various experimental data. We calculate the neutron star properties with the equation of state composed of nucleons and leptons in $\it beta-equilibrium$ which are in good agreement with the X-ray observations by Steiner and N\"{a}ttil\"{a}. We find that the maximum mass of the neutron star to be 2.15$M_{\odot}$ and stellar radius 11.936 km . Moreover, the radius and tidal deformability of a {\it canonical} neutron star mass 1.4$M_\odot$ come out to be 13.242 km and 3.910$\times$10$^{36}$ g cm$^2$ s$^2$ respectively within this parameter set.

The spinodal instabilities in hot asymmetric nuclear matter and some important critical parameters derived thereof are studied using six different families of relativistic mean-field (RMF) models. The slopes of the symmetry energy coefficient vary over a wide range within each family. The critical densities and proton fractions are more sensitive to the symmetry energy slope parameter at temperatures much below its critical value ($T_c\sim$14-16 MeV). The spread in the critical proton fraction at a given symmetry energy slope parameter is noticeably larger near $T_c$, indicating that the warm equation of state of asymmetric nuclear matter at sub-saturation densities is not sufficiently constrained. The distillation effects are sensitive to the density dependence of the symmetry energy at low temperatures which tend to wash out with increasing temperature.

We examine the correlations of neutron star radii with the nuclear matter incompressibility, symmetry energy, and their slopes, which are the key parameters of the equation of state (EoS) of asymmetric nuclear matter. The neutron star radii and the EoS parameters are evaluated using a representative set of 24 Skyrme-type effective forces and 18 relativistic mean field models, and two microscopic calculations, all describing 2$M_\odot$ neutron stars. Unified EoSs for the inner-crust-core region have been built for all the phenomenological models, both relativistic and non-relativistic. Our investigation shows the existence of a strong correlation of the neutron star radii with the linear combination of the slopes of the nuclear matter incompressibility and the symmetry energy coefficients at the saturation density. Such correlations are found to be almost independent of the neutron star mass in the range $0.6\text{-}1.8M_{\odot}$. This correlation can be linked to the empirical relation existing between the star radius and the pressure at a nucleonic density between one and two times saturation density, and the dependence of the pressure on the nuclear matter incompressibility, its slope and the symmetry energy slope. The slopes of the nuclear matter incompressibility and the symmetry energy coefficients as estimated from the finite nuclei data yield the radius of a $1.4M_{\odot}$ neutron star in the range $11.09\text{-}12.86$ km.

We study the diversities in the properties of the neutron stars arising due to the different choices for the cross-coupling between various mesons which governs the density dependence of the nuclear symmetry energy in the extended relativistic mean-field(RMF) model. For this purpose, we obtain two different families of the extended RMF model corresponding to different non-linear cross-coupling term in the isovector part of the effective Lagrangian density. The lowest order contributions for the $\delta$ mesons are also included. The different models within the same family are so obtained that they yield wide variation in the value of neutron-skin thickness in the $^{208}$Pb nucleus. These models are employed to compute the neutron star properties such as, core-crust transition density, radius and red shift at canonical mass ($1.4M_{\odot}$), tidal polarizability parameter, and threshold mass required for the enhanced cooling through direct Urca process. Most of the neutron star properties considered are significantly different(10\%-40\%) for the different families of models at a smaller neutron-skin thickness ($\sim 0.15$ fm) in the $^{208}$Pb nucleus.

Empirically determined values of the nuclear volume and surface symmetry energy coefficients from nuclear masses are expressed in terms of density distributions of nucleons in heavy nuclei in the local density approximation. This is then used to extract the value of the symmetry energy slope parameter $L$. The density distributions in both spherical and well deformed nuclei calculated within microscopic framework with different energy density functionals give $L = 59.0 \pm 13.0$ MeV. Application of the method also helps in a precision determination of the neutron skin thickness of nuclei that are difficult to measure accurately.

The density dependence of nuclear symmetry energy remains poorly constrained. Starting from precise empirical values of the nuclear volume and surface symmetry energy coefficients and the nuclear saturation density, we show how in the ambit of microscopic calculations with different energy density functionals, the value of the symmetry energy slope parameter $L$ alongwith that for neutron skin can be put in tighter bounds. The value of $L$ is found to be $L$= 64$\pm $5 MeV. For $^{208}$Pb, the neutron skin thickness comes out to be 0.188 $\pm $0.014 fm. Knowing $L$, the method can be applied to predict neutron skins of other nuclei.

The possibility of existence of hyperons in the recently measured $2M_\odot$ pulsar PSRJ1614-2230 is explored using a diverse set of nuclear equations of state calculated within the relativistic mean-field models. Our results indicate that the nuclear equations of state compatible with heavy-ion data allow the hyperons to exist in the PSRJ1614-2230 only for significantly larger values for the meson-hyperon coupling strengths. The maximum mass configurations for these cases contain sizable hyperon fractions ($\sim 60%$) and yet masquared their counterpart composed of only nucleonic matter.

The recent high-resolution measurement of the electric dipole (E1) polarizability (alphad) in 208Pb [Phys. Rev. Lett. 107, 062502 (2011)] provides a unique constraint on the neutron-skin thickness of this nucleus. The neutron-skin thickness (rskin) of 208Pb is a quantity of critical importance for our understanding of a variety of nuclear and astrophysical phenomena. To assess the model dependence of the correlation between alphad and rskin, we carry out systematic calculations for 208Pb, 132Sn, and 48Ca based on the nuclear density functional theory (DFT) using both non-relativistic and relativistic energy density functionals (EDFs). Our analysis indicates that whereas individual models exhibit a linear dependence between alphad and rskin, this correlation is not universal when one combines predictions from a host of different models. By averaging over these model predictions, we provide estimates with associated systematic errors for rskin and alphad for the nuclei under consideration. We conclude that precise measurements of rskin in both 48Ca and 208Pb---combined with the recent measurement of alphad---should significantly constrain the isovector sector of the nuclear energy density functional.

The stable configurations of non-rotating and rotating hybrid stars composed of colour superconducting quark matter core are constructed using several equations of state (EOSs). We use a set of diverse EOSs for the nuclear matter which represents the low density phase. The EOSs at higher densities correspond to the quark matter in the colour superconducting phase and are computed within the NJL-like model for different values of the scalar diquark and vector current couplings strengths. The phase transition to the quark matter is computed by a Maxwell construction. We find that the stability of the hybrid stars are mainly governed by the behaviour of the EOSs for the colour superconducting quark matter. However the compositions of hybrid star are sensitive to the EOS of the nuclear matter. The value of the critical rotation frequency for the hybrid star depends strongly on the EOS of the nuclear matter as well as that for the colour superconducting quark matter. Our results indicate that the EOS for the colour superconducting quark matter can be obtained, by adjusting the parameters of the NJL model, to yield the stable configurations of the hybrid star having the maximum mass $\sim 1.5M_\odot$ in the non-rotating limit and the critical rotation frequency $\sim 1$ kHz.

We construct static and mass-shedding limit sequences of hybrid stars, composed of colour flavour locked (CFL) quark matter core, for a set of equations of state (EOSs). The EOS for the hadronic matter is obtained using appropriately calibrated extended field theoretical based relativistic mean-field model. The MIT bag model is employed to compute the EOSs of the CFL quark matter for different values of the CFL gap parameter in the range of $50 - 150\text{MeV}$ with the deconfinement phase transition density ranging from $4\rho_0 - 6\rho_0$ ($\rho_0 = 0.16\text{fm}^{-3}$). We find, depending on the values of the CFL gap parameter and the deconfinement phase transition density, the sequences of stable configurations of hybrid stars either form third families of the compact stars or bifurcate from the hadronic sequence. The hybrid stars have masses $1.0 - 2.1 M_\odot$ with radii $9 - 13.5$ km. The maximum values of mass shedding limit frequency for such hybrid stars are $1 -2$ kHz. For the smaller values of the CFL gap parameter and the deconfinement phase transition density, mass-radius relationships are in harmony with those deduced by applying improved hydrogen atmosphere model to fit the high quality spectra from compact star X7 in the globular cluster 47 Tucanae. We observed for some cases that the third family of compact stars exist in the static sequence, but, disappear from the mass-shedding limit sequence. Our investigation suggests that the third family of compact stars in the mass-shedding limit sequence is more likely to appear, provided they have maximum mass in the static limit higher than their second family counterpart composed of pure hadronic matter.

Correlations in the properties of the static compact stars (CSs) and the ones rotating with the highest observed frequency of 1122Hz are studied using a large set of equations of state (EOSs). These EOSs span various approaches and their chemical composition vary from the nucleons to hyperons and quarks in $\beta$-equilibrium. It is found that the properties of static CS, like, the maximum gravitational mass $M_{\rm max}^{\rm stat}$ and radius $R_{1.4}^{\rm stat}$ corresponding to t he canonical mass and supramassive or non-supramassive nature of the CS rotating at 1122 Hz are strongly correlated. In particular, only those EOSs yield the CS rotating at 1122Hz to be non-supramassive for which $\left (\frac{M_{\rm max}^{\rm stat}}{M_\odot}\right )^{1/2} \left (\frac{10{\rm km}}{R_{1.4}^{\rm stat}})^{3/2} $ is greater than unity. Suitable parametric form which can be used to split the $M_{\rm max}^{\rm stat}$ $-$ $R_{1.4}^{\rm stat}$ plane into the regions of different supramassive nature of the CS rotating at 1122Hz is presented. Currently measured maximum gravitational mass 1.76$M_\odot$ of PSR J0437-4715 suggests that the CS rotating at 1122Hz can be non-supramassive provided $R_{1.4}^{\rm stat} \leqslant 12.4$ km.

We study the properties of non-rotating and rotating neutron stars for a new set of equations of state (EOSs) with different high density behaviour obtained using the extended field theoretical model. The high density behaviour for these EOSs are varied by varying the $\omega-$meson self-coupling and hyperon-meson couplings in such a way that the quality of fit to the bulk nuclear observables, nuclear matter incompressibility coefficient and hyperon-nucleon potential depths remain practically unaffected. We find that the largest value for maximum mass for the non-rotating neutron star is $2.1M_\odot$. The radius for the neutron star with canonical mass is $12.8 - 14.1$ km provided only those EOSs are considered for which maximum mass is larger than $1.6M_\odot$ as it is the lower bound on the maximum mass measured so far. Our results for the very recently discovered fastest rotating neutron star indicate that this star is supra massive with mass $1.7 - 2.7M_\odot$ and circumferential equatorial radius $12 - 19$ km.