We construct for the first time, the sequences of stable neutron star (NS) models capable of explaining simultaneously, the glitch healing parameters, $Q$, of both the pulsars, the Crab ($Q \geq 0.7$) and the Vela ($Q \leq 0.2$), on the basis of starquake mechanism of glitch generation, whereas the conventional NS models cannot give such consistent explanation. Furthermore, our models also yield an upper bound on NS masses similar to those obtained in the literature for a variety of modern equations of state (EOSs) compatible with causality and dynamical stability. The results obtained in this paper set an upper and lower bound on the energy of a gravitationally redshifted electron-positron annihilation line in the range of about 0.309 - 0.315 MeV from the Crab and in the range of about 0.425 - 0.439 MeV from the Vela pulsar.

We have obtained a criterion for spherically symmetric and static structures under hydrostatic equilibrium in general relativity (GR), which states that for a given value of $\sigma \equiv (P_0/E_0) \equiv $ the ratio of central pressure to central energy-density], the compaction parameter $u \equiv (M/a)$, where $M$ is the total mass and $a$ is the radius of the configuration] or the surface redshift of any regular configuration cannot exceed that of the corresponding homogeneous density sphere, that is, $u \leq u_h$, where $u_h$ is the compaction parameter of the homogeneous density sphere. By examining various exact solutions and surface). On the other hand, configurations having a finite density on the surface (that is, the self-bound structures) do not fulfill this criterion. This criterion puts a severe restriction on the static structures based upon the general relativistic field equations and consequently on the upper limit of mass, surface and central redshift and other physical parameters of spherically symmetric and static configurations.

We show that the mass-radius $(M-R)$ relation corresponding to the MIT bag models of strange quark matter (SQM) and the models obtained by Day et al (1998) do not provide the necessary and sufficient condition for dynamical stability for the equilibrium configurations, since such configurations can not even fulfill the necessary condition of hydrostatic equilibrium provided by the exterior Schwarzschild solution. These findings will remain unaltered and can be extended to any other sequence of pure SQM. This study explicitly show that although the strange quark matter might exist in the state of zero pressure and temperature, but the models of pure strange quark `stars' can not exist in the state of hydrostatic equilibrium on the basis of General Relativity Theory. This study can affect the results which are claiming that various objects like - RX J1856.5-3754, SAX J1808.4-3658, 4U 1728-34, PSR 0943+10 etc. might be strange stars.

The various changes in the `abstract' are mentioned in point (1) of the present erratum.

We show that the mass-radius $(M-R)$ relation corresponding to the stiffest equation of state (EOS) does not provide the necessary and sufficient condition of dynamical stability for the equilibrium configurations, since such configurations can not satisfy the `compatibility criterion'. In this connection, we construct sequences composed of core-envelope models such that, like the stiffest EOS, each member of these sequences satisfy the extreme case of causality condition, $v = c = 1$, at the centre. We, thereafter, show that the $M-R$ relation corresponding to the said core-envelope model sequences can provide the necessary and sufficient condition of dynamical stability only when the `compatibility criterion' for these sequences is `appropriately' satisfied. However, the fulfillment of `compatibility criterion' can remain satisfied even when the $M-R$ relation does not provide the necessary and sufficient condition of dynamical stability for the equilibrium configurations. In continuation to the results of previous study, these results explicitly show that the `compatibility criterion' {\em independently} provides, in general, the {\em necessary} and {\em sufficient} condition of hydrostatic equilibrium for any regular sequence. Beside its fundamental feature, this study can also explain simultaneously, both (the higher as well as lower) values of the glitch healing parameter observed for the Crab and the Vela-like pulsars respectively, on the basis of starquake model of glitch generation.

We present explicit examples to show that the `compatibility criterion' is capable of providing a {\em necessary} and {\em sufficient} condition for any regular configuration to be compatible with the state of hydrostatic equilibrium. This conclusion is drawn on the basis of the finding that the $M-R$ relation gives the necessary and sufficient condition for dynamical stability of equilibrium configurations only when the compatibility criterion for these configurations is appropriately satisfied. In this regard, we construct an appropriate sequence composed of core-envelope models on the basis of compatibility criterion, such that each member of this sequence satisfies the extreme case of causality condition $v = c = 1$ at the centre. The maximum stable value of $u \simeq 0.3389$ (which occurs for the model corresponding to the maximum value of mass in the mass-radius relation) and the corresponding central value of the local adiabatic index, $(\Gamma_1)_0 \simeq 2.5911$, of this model are found fully consistent with those of the corresponding {\em absolute} values, $u_{\rm max} \leq 0.3406$, and $(\Gamma_1)_0 \leq 2.5946$, which impose strong constraints on these parameters of such models. In addition to this example, we also study dynamical stability of pure adiabatic polytropic configurations on the basis of variational method for the choice of the `trial function', $\xi =re^{\nu/4}$, as well as the mass-central density relation, since the compatibility criterion is appropriately satisfied for these models. The results of this example provide additional proof in favour of the statement regarding compatibility criterion mentioned above.

Various neutron star (NS) models correspond to an unrealistic mass range $\leq 0.5M_\odot$ for the Vela pulsar, if the observational constraints of the glitch healing parameter in the starquake model ($Q = I_{\rm core}/I_{\rm total} \leq 0.2$), are imposed on these models. However, we show that these observational constraints yield a realistic mass range for NS models, corresponding to a core given by the stiffest equations of state (EOS) and the envelope is characterized by the well known EOS of adiabatic polytrope, if the continuity of the adiabatic speed of sound together with pressure, energy-density, and the two metric parameters is assured at the core-envelope boundary of the models and this boundary is worked out on the basis of the `compatibility criterion' for hydrostatic equilibrium. The models yield the surface redshift $z_R \simeq 0.6913$ and mass $M \simeq 2.153 M_\odot$ for the ``central'' weighted mean value, $Q = 0.12 \pm 0.07$, of the glitch healing parameter of the Vela pulsar. These values of mass and surface redshift can increase upto $M \simeq 2.196 M_\odot$ and $z_R \simeq 0.7568$ (which represents an ultra-compact object (UCO; $z_R \geq 0.73$)) respectively for the upper weighted mean value of $Q \simeq 0.19$. However, for the lower weighted mean value of $Q \simeq 0.05$, the mass and surface redshift can decrease upto the values of $M \simeq 2.052 M_\odot$ and $z_R \simeq 0.6066$ respectively. The observation of the lower bound on the energy of a $\gamma$-ray pulse at about 0.30 MeV from the Vela pulsar in 1984 is in excellent agreement with the results of this study, provided this energy could be interpreted as the energy of a gravitationally redshifted electron-positron annihilation radiation from the star's surface.

An analysis of insular solutions of Einstein's field equations for static, spherically symmetric, source mass, on the basis of exterior Schwarzschild solution is presented. Following the analysis, we demonstrate that the {\em regular} solutions governed by a self-bound (that is, the surface density does not vanish together with pressure) equation of state (EOS) or density variation can not exist in the state of hydrostatic equilibrium, because the source mass which belongs to them, does not represent the `actual mass' appears in the exterior Schwarzschild solution. The only configuration which could exist in this regard is governed by the homogeneous density distribution (that is, the interior Schwarzschild solution). Other structures which naturally fulfill the requirement of the source mass, set up by exterior Schwarzschild solution (and, therefore, can exist in hydrostatic equilibrium) are either governed by gravitationally-bound regular solutions (that is, the surface density also vanishes together with pressure), or self-bound singular solutions (that is, the pressure and density both become infinity at the centre).

Spiral galaxies are considered as static, spherically symmetric Dark Matter Configurations (DMC)in which non-zero rest-mass particles (NZRPs) move along appropriate trajectories. Using general relativity (GR), we show that a mass of dark matter about 127 - 212 $\times 10^{11} M_{\odot}$ is required inside the sphere of size $a \sim 10 $ kpc for agreement with the observed typical orbital velocity ranging from 150 to 250 km sec$^{-1}$. In our model, it is possible to obtain flat, slightly rising, and even declining rotation curves. The model self-evidently explains the reason why most of the spirals show almost flat or slightly rising rotation curves and very few show declining curves. The NZRPs follow trajectories, and are finally trapped in circles of minimum radii $\sim 0.03a - 0.005a$ with rotation velocities reaching $\sim 6544-8445$ km sec $^{-1}$. The velocity values are consistent with observations of central regions of various galaxies. This scenario may also lead to conditions suitable for evolution of relativistic supermassive stars and supermassive black holes at the centres of various galaxies. Based on this study, we obtain the density parameter $\Omega \sim 0.127 - 0.212 $, leading to an open model of the Universe with an age of $\sim 16.8 - 17.6 $ Gyr, which is significantly higher than the presently estimated age of globular clusters $\sim 13 - 15 $ Gyr.

If the causality condition [the speed of sound always remains less than that of light in vacuum, i. e., $v \leq c = 1$] is imposed on the spheres of homogeneous energy density, the `ratio of the specific heats', $\gamma \leq 2.59457$, constraints the compaction parameter, $u [\equiv (M/a)$, mass to size ratio in geometrized units] of the dynamically stable configurations $ \leq 0.34056 $ [corresponding to a surface redshift $(z_a) \leq 0.771$]. Apparently, The maximum value of $u$ obtained in this manner belongs to an absolute upper bound, and gives: (i) The maximum value for static neutron star masses as $5.4 M_\odot$, if we substitute the density at the surface of the configuration equal to the average nuclear density, $E = 2 \times 10^{14}$ g cm$^{-3}$ [e.g. {\it Nature}, {\bf 259}, 377 (1976)]. (ii) However, if the density of the static configuration is constrained to the value $1.072 \times 10^{14}$ g cm$^{-3}$, by imposing the empirical result that the minimum rotation period of the fastest rotating pulsar known to date, PSR 1937 + 21, is 1.558 ms, the maximum mass value for static neutron stars exceed upto $7.4 M_\odot$. These masses have important implications for the massive compact objects like Cyg X-1, Cyg XR-1, and LMC-X3 etc., which may not, necessarily, represent black holes. (iii) The minimum rotation periods for a static $1.442 M_\odot$ neutron star to be 0.3041 ms. (iv) A suitable stable model of ultra-compact objects [$u > (1/3)$] which has important astrophysical significance.