It is usually thought that a single equation of state (EoS) model "correctly" represents cores of all compact stars. Here we emphasize that two families of compact stars, viz., neutron stars and strange stars, can coexist in nature, and that neutron stars can get converted to strange stars through the nucleation process of quark matter in the stellar center. From our fully general relativistic numerical computations of the structures of fast-spinning compact stars, known as millisecond pulsars, we find that such a stellar conversion causes a simultaneous spin-up and decrease in gravitational mass of these stars. This is a new type of millisecond pulsar evolution through a new mechanism, which gives rise to relatively lower mass compact stars with higher spin rates. This could have implication for the observed mass and spin distributions of millisecond pulsars. Such a stellar conversion can also rescue some massive, spin-supported millisecond pulsars from collapsing into black holes. Besides, we extend the concept of critical mass $M_{\rm cr}$ for the neutron star sequence (Berezhiani et al. 2003; Bombaci et al. 2004) to the case of fast-spinning neutron stars, and point out that neutron star EoS models cannot be ruled out by the stellar mass measurement alone. Finally, we emphasize the additional complexity for constraining EoS models, for example, by stellar radius measurements using X-ray observations, if two families of compact stars coexist.

The recent fast growth of a population of millisecond pulsars with precisely measured mass provides an excellent opportunity to characterize these compact stars at an unprecedented level. This is because the stellar parameter values can be accurately computed for known mass and spin rate and an assumed equation of state (EoS) model. For each of the 16 such pulsars and for a set of EoS models from nucleonic, hyperonic, strange quark matter and hybrid classes, we numerically compute fast spinning stable stellar parameter values considering the full effect of general relativity. This first detailed catalogue of the computed parameter values of observed millisecond pulsars provides a testbed to probe the physics of compact stars, including their formation, evolution and EoS. We estimate uncertainties on these computed values from the uncertainty of the measured mass, which could be useful to quantitatively constrain EoS models. We note that the largest value of the central density $\rho_{\rm c}$ in our catalogue is $\sim 5.8$ times the nuclear saturation density $\rho_{\rm sat}$, which is much less than the expected maximum value $13 \rho_{\rm sat}$. We argue that the $\rho_{\rm c}$-values of at most a small fraction of compact stars could be much larger than $5.8 \rho_{\rm sat}$. Besides, we find that the constraints on EoS models from accurate radius measurements could be significantly biased for some of our pulsars, if stellar $spinning$ configurations are not used to compute the theoretical radius values.

For a set of equation of state (EoS) models involving interacting strange quark matter, characterized by an effective bag constant (B_eff) and a perturbative QCD corrections term (a_4), we construct fully general relativistic equilibrium sequences of rapidly spinning strange stars for the first time. Computation of such sequences is important to study millisecond pulsars and other fast spinning compact stars. Our EoS models can support a gravitational mass (M_G) and a spin frequency at least up to approximately 3.0 solar mass and approximately 1250 Hz respectively, and hence are fully consistent with measured M_G and spin frequency values. This paper reports the effects of B_eff and a_4 on measurable compact star properties, which could be useful to find possible ways to constrain these fundamental quark matter parameters, within the ambit of our EoS models. We confirm that a lower B_eff allows a higher mass. Besides, for known M_G and spin frequency, measurable parameters, such as stellar radius, radius-to-mass ratio and moment of inertia, increase with the decrease of B_eff. Our calculations also show that a_4 significantly affects the stellar rest mass and the total stellar binding energy. As a result, a_4 can have signatures in evolutions of both accreting and non-accreting compact stars, and the observed distribution of stellar mass and spin and other source parameters. Finally, we compute the parameter values of two important pulsars, PSR J1614-2230 and PSR J1748-2446ad, which may have implications to probe their evolutionary histories, and for constraining EoS models.

We initiate numerical studies of differentially rotating magnetised (proto) neutron stars by studying - through construction from first principles - the coupling between an assumed differential rotation and an impressed magnetic field. For a perfect incompressible, homogeneous, non-dissipative fluid sphere immersed in an ambient plasma, we solve the (coupled) azimuthal components of the Navier-Stokes equation and the Maxwell induction equation. The assumed time--independent poloidal field lines get dragged by the rotating fluid and produce toroidal magnetic fields. Surface magnetic fields take away energy redistributing the angular momentum to produce rigid rotation along poloidal field lines. Due to absence of viscous dissipation, sustained torsional oscillations are set up within the star. However, the perpetual oscillations of neighbouring `closed' field lines get increasingly out of phase with time, leading to structure build up as in Liu & Shapiro (2004) implying the importance of taking into account diffusion (Spruit 1999) for realistic modeling.

Quantized orbital structures are typical for many aspects of classical gravity (Newton's as well as Einstein's). The astronomical phenomenon of orbital resonances is a well-known example. Recently, Rothman, Ellis and Murugan (2001) discussed quantized orbital structures in the novel context of a holonomy invariance of parallel transport in Schwarzschild geometry. We present here yet another example of quantization of orbits, reflecting both orbital resonances and holonomy invariance. This strong-gravity effect may already have been directly observed as the puzzling kilohertz quasi-periodic oscillations (QPOs) in the X-ray emission from a few accreting galactic black holes and several neutron stars.

We compute spectra from accretion disks around rapidly rotating neutron stars. The full effect of general relativity is considered for the structure calculation of the stars. We take into account the Doppler shift, gravitational redshift and light-bending effects in order to compute the observed spectra. To facilitate direct comparison with observations, a simple empirical function is presented which describes the numerically computed spectra well. This function can in principle be used to distinguish between the Newtonian spectra and the relativistic spectra. We also discuss the possibility of constraining neutron star's equation of state using our spectral models.

Structure parameters of rapidly rotating strange stars are compared with those for neutron stars. The salient differences are pointed out.

Disk accretion onto a neutron star possessing a weak surface magnetic field ($B \le 10^8$ G) provides interesting X-ray emission scenarios, and is relevant for understanding X-ray bursters and low-mass X-ray binaries (LMXBs). The standard (Newtonian) theory of disk-accretion predicts that the matter spiralling in from infinity loses one-half of its total gravitational energy in the extended disk, and the remainder in a narrow boundary layer girdling the neutron star. The ratio of the boundary layer luminosity to that from the disk ($E_{\rm BL}/E_{\rm D}$) is, therefore, unity. On incorporation of general relativity without rotation (Schwarzschild solution), $E_{\rm BL}/E_{\rm D}$ is seen to be as high as 6. We construct rotating sequences of neutron stars for three representative equations of state. We show here that for a neutron star rotating at a limit where centrifugal force balances the inward gravitational force, $E_{\rm BL}/E_{\rm D}\sim 0$.

We compute the temperature profiles of accretion discs around rapidly rotating strange stars, using constant gravitational mass equilibrium sequences of these objects, considering the full effect of general relativity. Beyond a certain critical value of stellar angular momentum ($J$), we observe the radius ($r_{\rm orb}$) of the innermost stable circular orbit (ISCO) to increase with J (a property seen neither in rotating black holes nor in rotating neutron stars). The reason for this is traced to the crucial dependence of $dr_{\rm orb}/dJ$ on the rate of change of the radial gradient of the Keplerian angular velocity at $r_{\rm orb}$ with respect to $J$. The structure parameters and temperature profiles obtained are compared with those of neutron stars, as an attempt to provide signatures for distinguishing between the two. We show that when the full gamut of strange star equation of state models, with varying degrees of stiffness are considered, there exists a substantial overlap in properties of both neutron stars and strange stars. However, applying accretion disc model constraints to rule out stiff strange star equation of state models, we notice that neutron stars and strange stars exclusively occupy certain parameter spaces. This result implies the possibility of distinguishing these objects from each other by sensitive observations through future X-ray detectors.

We present computed spectra, as seen by a distant observer, from the accretion disc around a rapidly rotating neutron star. Our calculations are carried out in a fully general relativistic framework, with exact treatment of rotation. We take into account the Doppler shift, gravitational redshift and light bending effects in order to compute the observed spectrum. We find that light bending significantly modifies the high-energy part of the spectrum. Computed spectra for slowly rotating neutron stars are also presented. These results would be important for modelling the observed X-ray spectra of Low Mass X-ray Binaries containing fast spinning neutron stars.

We calculate the temperature profiles of (thin) accretion disks around rapidly rotating neutron stars (with low surface magnetic fields), taking into account the full effects of general relativity. We then consider a model for the spectrum of the X-ray emission from the disk, parameterized by the mass accretion rate, the color temperature and the rotation rate of the neutron star. We derive constraints on these parameters for the X-ray source Cygnus X-2 using the estimates of the maximum temperature in the disk along with the disk and boundary layer luminosities, using the spectrum inferred from the EXOSAT data. Our calculations suggest that the neutron star in Cygnus X-2 rotates close to the centrifugal mass-shed limit. Possible constraints on the neutron star equation of state are also discussed.

General relativistic spectra from accretion disks around rotating neutron stars in the appropriate space-time geometry for several different equation of state, spin rates and mass of the compact object have been computed. The analysis involves the computation of the relativistically corrected radial temperature profiles and the effect of Doppler and gravitational red-shifts on the spectra. Light bending effects have been omitted for simplicity. The relativistic spectrum is compared with the Newtonian one and it is shown that the difference between the two is primarily due to the different radial temperature profile for the relativistic and Newtonian disk solutions. To facilitate direct comparison with observations, a simple empirical function has been presented which describes the numerically computed relativistic spectra well. This empirical function (which has three parameters including normalization) also describes the Newtonian spectrum adequately. Thus the function can in principle be used to distinguish between the two. In particular, the best fit value of one of the parameters ($\beta$-parameter) $\approx 0.4$ for the Newtonian case, while it ranges from 0.1 to 0.35 for relativistic case depending upon the inclination angle, EOS, spin rate and mass of the neutron star. Constraining this parameter by fits to future observational data of X-ray binaries will indicate the effect of strong gravity in the observed spectrum.

For a new equation of state of strange quark matter, we construct equilibrium sequences of rapidly rotating strange stars in general relativity. The sequences are the normal and supramassive evolutionary sequences of constant rest mass. We also calculate equilibrium sequences for a constant value of $\Omega$ corresponding to the most rapidly rotating pulsar PSR 1937 + 21. In addition to this, we calculate the radius of the marginally stable orbit and its dependence on $\Omega$, relevant for modeling of kilo-Hertz quasi-periodic oscillations in X-ray binaries.

Kilohertz quasi-periodic brightness oscillations (kHz QPOs) observed in certain X-ray burst sources may represent Keplerian frequencies in the inner regions of the accretion disk in such systems. If this assumption is strictly adhered to, we show here that if the central accretor in stellar X-ray burst sources is a strange star (made up of u, d and s quarks in beta equilibrium, referred to as strange matter) then the calculated QPO frequencies are reconcilable with the observed QPO frequencies (corresponding to the highest frequency of 1.22 kHz, observed so far from the source 4U 1636-53) only for particular values of the QCD-related parameters which describe the equation of state of strange matter. We demonstrate that QPO frequencies in the very high range (1.9-3.1) kHz can be understood in terms of a (non- magnetized) strange star X-ray binary (SSXB) rather than a neutron star X-ray binary (NSXB). Future discovery of such high frequency QPOs from X-ray burst sources will constitute a new astrophysical di- agnostic for identifying solar mass range stable strange stars in our galaxy.

Kilohertz quasi-periodic oscillations (QPOs) in low-mass X-ray binaries (LMXBs) may represent the orbital frequencies of the innermost Keplerian orbits around accreting neutron stars. Attempts have recently been made to derive constraints on the mass and the equation of state of the neutron star, by identifying the highest observed QPO frequency with the Keplerian frequency at the marginally stable orbit given by the general theory of relativity. These estimates have either neglected the effect of neutron star rotation or used an approximate treatment of rotation in general relativity. We rederive these constraints using a fully general relativistic formalism including the effect of rapid rotation. We also present constraints corresponding to the case where the innermost stable orbit touches the stellar surface.

For four newly suggested microscopic equations of state of neutron star matter, we construct equilibrium sequences of rapidly rotating neutron stars in general relativity. The sequences are the normal and supramassive evolutionary sequences of constant rest mass. We find that for these equations of state the maximum (gravitational) mass rotating models occur (in central density and rotation rate $\Omega$) before the maximum--$\Omega$ models. We calculate equilibrium sequences for a constant value of $\Omega$ corresponding to the most rapidly rotating pulsar PSR 1937+21. Also calculated is the radius of the marginally stable orbit and its dependence on $\Omega$.

We calculate the disk and boundary layer luminosities for accreting rapidly rotating neutron stars with low magnetic fields in a fully general relativistic manner. Rotation increases the disk luminosity and decreases the boundary layer luminosity. A rapid rotation of the neutron star substantially modifies these quantities as compared to the static limit. For a neutron star rotating close to the centrifugal mass shed limit, the total luminosity has contribution only from the extended disk. For such maximal rotation rates, we find that much before the maximum stable gravitational mass configuration is reached, there exists a limiting central density, for which particles in the innermost stable orbit will be more tightly bound than those at the surface of the neutron star. We also calculate the angular velocity profiles of particles in Keplerian orbits around the rapidly rotating neutron star. The results are illustrated for a representative set of equation of state models of neutron star matter.