LAUSR

The results of the investigation of the core-envelope model presented in Negi et al. (Gen. Relativ. Gravit. 22:735, 1990 ) have been discussed in view of the reference (Negi et al. in Gen. Relativ. Gravit. 51:131, 2019 ). It is seen that there are significant changes in the results to be addressed. In addition, I have also calculated the gravitational binding energy, causality and pulsational stability of the structures which were not considered in Negi et al. (Gen. Relativ. Gravit. 22:735, 1990 ). The modified results have important consequences to model neutron stars and pulsars. The maximum neutron star mass obtained in this study corresponds to the mean value of the classical results obtained by Rhoades and Ruffini (Phys. Rev. Lett. 32:324, 1974 ) and the upper bound on neutron star mass obtained by Kalogera and Baym (Astrophys. J. 470:L61, 1996 ) and is much closer to the most recent theoretical estimate made by Sotani (Phys. Rev. C 95:025802, 2017 ). On one hand, when there are only few equations of state (EOSs) available in the literature which can fulfil the recent observational constraint imposed by the largest neutron star masses around 2 M ⊙ $M_{\odot }$ (Demorest et al. in Nature 467:1081, 2010 ; Antoniadis et al. in Science 340:6131, 2013 ; Cromartie et al. in Nat. Astron. 4:72, 2020 ), the present analytic models, on the other hand, can comfortably satisfy this constraint. Furthermore, the maximum allowed value of compactness parameter u $u$ ( ≡ M / a $\equiv M/a$ ; mass to size ratio in geometrized units) ≤ 0.30 $\leq 0.30$ obtained in this study is also consistent with an absolute maximum value of u max = 0.333 − 0.005 + 0.001 $u_{\mathrm{max}} = 0.333^{+0.001}_{-0.005}$ resulting from the observation of binary neutron stars merger GW170817 (see, e.g. Koliogiannis and Moustakidis in Astrophys. Space Sci. 364:52, 2019 ).