Based on the general analysis of the grand canonical partition function in the 𝑆-matrix framework, a method is presented to calculate the equation of state of dilute warm nuclear matter. The result is a model-independent virial series for the pressure and density that systematically includes contributions from all the ground and excited states of all the stable nuclear species and their scattering channels. The multiplicity distribution of these species to keep the matter in statistical equilibrium is found out and then the pressure, incompressibility and the symmetry energy of the system are evaluated. The calculated symmetry energy coefficients are found to be in fair agreement with the recent experimental data.

The nature of equation of state for the neutron star matter is crucially governed by the density dependence of the nuclear symmetry energy. We attempt to probe the behaviour of the nuclear symmetry energy around the saturation density by exploiting the empirical values for volume and surface symmetry energy coefficients extracted from the precise data on the nuclear masses.

The nuclear symmetry energy at a given density measures the energy transferred in converting symmetric nuclear matter into the pure neutron matter. The density content of nuclear symmetry energy remains poorly constrained. Our recent results for the density content of the nuclear symmetry energy, around the saturation density, extracted using experimental data for accurately known nuclear masses, giant resonances and neutron-skin thickness in heavy nuclei are summarized.