We study the main astrophysical properties of differentially rotating neutron stars described as stationary and axisymmetric configurations of a moderately stiff $\Gamma=2$ polytropic fluid. The high level of accuracy and of stability of our relativistic multidomain pseudo-spectral code enables us to explore the whole solution space for broad ranges of the degree of differential rotation, but also of the stellar density and oblateness. Staying within an astrophysicaly motivated range of rotation profiles, we investigate the characteristics of neutron stars with maximal mass for all types of families of differentially rotating relativistic objects identified in a previous article Ansorg, Gondek-Rosinsla, Villain (2009). We find that the maximum mass depends on both the degree of differential rotation and on the type of solution. It turns out that the maximum allowed mass can be up to 4 times higher than what it is for non-rotating stars with the same equation of state. Such values are obtained for a modest degree of differential rotation but for one of the newly discovered type of solutions. Since such configurations of stars are not that extreme, this result may have important consequences for the gravitational wave signal to expect from coalescing neutron star binaries or from some supernovae events.

We present the results of a study of the low-mass X-ray binary 4U 1636-536. We have performed temporal analysis of all available RXTE/ASM, Swift/BAT and MAXI data. We have confirmed the previously discovered quasi-periodicity of ~45 d present during ~2004, however we found it continued to 2006. At other epochs, the quasi-periodicity is only transient, and the quasi-period, if present, drifts. We have then applied a time-dependent accretion disc model to the interval with the significant X-ray quasi-periodicity. For our best model, the period and the amplitude of the theoretical light curve agree well with that observed. The modelled quasi-periodicity is due to the hydrogen thermal-ionization instability occurring in outer regions of the accretion disc. The model parameters are the average mass accretion rate (estimated from the light curves), and the accretion disc viscosity parameters, for the hot and cold phases. Our best model gives relatively low values of viscosity parameter for cold phase 0.01 and for hot phase 0.03.

Kilohertz QPOs can be used as a probe of the inner regions of accretion disks in compact stars and hence also of the properties of the central object. Most models of kHz QPOs involve epicyclic frequencies to explain their origin. We compute the epicyclic frequencies of nearly circular orbits around rotating strange quark stars. The MIT bag model is used to model the equation of state of quark matter and the uniformly rotating stellar configurations are computed in full general relativity. The vertical epicyclic frequency and the related nodal precession rate of inclined orbits are very sensitive to the oblateness of the rotating star. For slowly rotating stellar models of moderate and high mass strange stars, the sense of the nodal precession changes at a certain rotation rate. At lower stellar rotation rates the orbital nodal precession is prograde, as it is in the Kerr metric, while at higher rotation rates the precession is retrograde, as it is for Maclaurin spheroids. Thus, qualitatively, the orbits around rapidly rotating strange quark stars are affected more strongly by the effects of stellar oblateness than by the effects of general relativity. We show that epicyclic and orbital frequencies calculated numerically for small mass strange stars are in very good agreement with analytical formulae for Maclaurin spheroids.

Using the StarTrack binary population synthesis code we model the population of double neutron stars in the Galaxy. We include a detailed treatment of the spin evolution of each pulsar due to processes such as spin-down and spin-up during accretion events as well as magnetic field decay. We also model the spatial distribution of double neutron stars by including their natal kicks and subsequent propagation in the Galactic gravitational potential. This synthetic pulsar population is compared to the observed sample of double neutron stars taking into account the selection effects of detection in the radio band, to determine the most likely evolutionary parameters. With these parameters we determine the properties of the double neutron star binaries detectable in gravitational waves by the high frequency interferometers LIGO and VIRGO. In particular, we discuss the distributions of chirp masses and mass ratios in samples selected by their radio or gravitational wave emission.

The current gravitational wave detectors have reached their operational sensitivity and are nearing detection of compact object binaries. In the coming years, we expect that the Advanced LIGO/VIRGO will start taking data. At the same time, there are plans for third generation ground-based detectors such as the Einstein Telescope, and space detectors such as DECIGO. We discuss the eccentricity distribution of inspiral compact object binaries during they inspiral phase. We analyze the expected distributions of eccentricities at three frequencies that are characteristic of three future detectors: Advanced LIGO/VIRGO (30 Hz), Einstein Telescope (3 Hz), and DECIGO (0.3 Hz). We use the StarTrack binary population code to investigate the properties of the population of compact binaries in formation. We evolve their orbits until the point that they enter a given detector sensitivity window and analyze the eccentricity distribution at that time. We find that the eccentricities of BH-BH and BH-NS binaries are quite small when entering the Advanced LIGO/VIRGO detector window for all considered models of binary evolution. Even in the case of the DECIGO detector, the typical eccentricities of BH-BH binaries are below 10^{-4}, and the BH-NS eccentricities are smaller than 10^{-3}. Some fraction of NS-NS binaries may have significant eccentricities. Within the range of considered models, we found that a fraction of between 0.2% and 2% NS-NS binaries will have an eccentricity above 0.01 for the Advanced LIGO/VIRGO detectors. For the ET detector, this fraction is between 0.4% and 4%, and for the DECIGO detector it lies between 2% and 27%.

We present the first relativistic calculations of the final phase of inspiral of a binary system consisting of two stars built predominantely of strange quark matter (strange quark stars). We study the precoalescing stage within the Isenberg-Wilson-Mathews approximation of general relativity using a multidomain spectral method. A hydrodynamical treatment is performed under the assumption that the flow is either rigidly rotating or irrotational, taking into account the finite density at the stellar surface -- a distinctive feature with respect to the neutron star case. The gravitational-radiation driven evolution of the binary system is approximated by a sequence of quasi-equilibrium configurations at fixed baryon number and decreasing separation. We find that the innermost stable circular orbit (ISCO) is given by an orbital instability both for synchronized and irrotational systems. This constrasts with neutron stars for which the ISCO is given by the mass-shedding limit in the irrotational case. The gravitational wave frequency at the ISCO, which marks the end of the inspiral phase, is found to be 1400 Hz for two irrotational 1.35 Msol strange stars and for the MIT bag model of strange matter with massless quarks and a bag constant B=60 MeV/fm^3. Detailed comparisons with binary neutrons star models, as well as with third order Post-Newtonian point-mass binaries are given.

We present calculations of quasiequilibrium sequences of irrotational binary neutron stars based on realistic equations of state (EOS) for the whole neutron star interior. Three realistic nuclear EOSs of various softness and based on different microscopic models have been joined with a recent realistic EOS of the crust, giving in this way three different EOSs of neutron-star interior. Computations of quasiequilibrium sequences are performed within the Isenberg-Wilson-Mathews approximation to general relativity. For all evolutionary sequences, the innermost stable circular orbit (ISCO) is found to be given by mass-shedding limit (Roche lobe overflow). The EOS dependence on the last orbits is found to be quite important: for two 1.35 M_sol neutron stars, the gravitational wave frequency at the ISCO (which marks the end of the inspiral phase) ranges from 800 Hz to 1200 Hz, depending upon the EOS. Detailed comparisons with 3rd order post-Newtonian results for point-mass binaries reveals a very good agreement until hydrodynamical effects (dominated by high-order functions of frequency) become important, which occurs at a frequency ranging from 500 Hz to 1050 Hz, depending upon the EOS.

We use the well tested StarTrack binary population synthesis code to examine the properties of the population of compact object binaries. We calculate the distribution of masses and mass ratios taking into account weights introduced by observability in gravitational waves during inspiral. We find that in the observability weighted distribution of double neutron star binaries there are two peaks: one for nearly equal mass systems, and one for the systems consisting of a low and a high mass neutron star, $q=0.6-0.7$. The observability weighted distribution of black hole neutron star binaries is concentrated on the systems with the mass ratio $q=0.3-0.5$, while for the double black hole binaries the observability weighted distribution is dominated by the massive, nearly equal mass binaries with $q>0.7$.

We study the viscosity driven (Jacobi-like) bar mode instability of rapidly rotating strange stars in general relativity. A triaxial, "bar shaped" compact star could be an efficient source of continuous wave gravitational radiation in the frequency range of the forthcoming interferometric detectors. We locate the secular instability point along several constant baryon mass sequences of uniformly rotating strange stars described by the MIT bag model. Contrary to neutron stars, strange stars with T/|W| (the ratio of the rotational kinetic energy to the absolute value of the gravitational potential energy) much lower than the corresponding value for the mass-shed limit can be secularly unstable to bar mode formation if shear viscosity is high enough to damp out any deviation from uniform rotation. The instability develops for a broad range of gravitational masses and rotational frequencies of strange quark stars. It imposes strong constraints on the lower limit of the frequency at the innermost stable circular orbit around rapidly rotating strange stars. The above results are robust for all linear self-bound equations of state assuming the growth time of the instability is faster than the damping timescale. We discuss astrophysical scenarios where triaxial instabilities (r-mode and viscosity driven instability) could be relevant in strange stars described by the standard MIT bag model of normal quark matter. Taking into account actual values of viscosities in strange quark matter and neglecting the magnetic field we show that Jacobi-like instability cannot develop in any astrophysicaly interesting temperature windows. The main result is that strange quark stars described by the MIT bag model can be accelerated to very high frequency in Low Mass X-ray binaries if the strange quark mass is ~ 200 MeV or higher.

The X-ray emission of RXJ1856.5-3754 has been found to coincide to unprecedented accuracy with that of a blackbody, of radius 5.8+-0.9 km for the measured parallax distance of 140 pc (Burwitz et al. 2001, Drake et al. 2002). If the emission is uniform over the whole surface of a non-rotating star, the mass of the star cannot exceed 0.75+-0.12 solar mass regardless of its composition. If the compact object is a quark star described by the MIT-bag equation of state (a ``strange star''), the mass is no more than 0.3 solar mass. Comparably small masses are also obtained for the X-ray bursters Aql X-1 and KS1731-260 for some fits to their spectra.

The spin-up of a neutron star crucially depends on the maximum orbital frequency around it, as do a host of other high energy accretion phenomena in low mass X-ray binaries, including quasi-periodic oscillations (QPOs) in the X-ray flux. We compare the maximum orbital frequencies for MIT-bag quark stars and for neutron stars modeled with the FPS equation of state. The results are based on relativistic calculations of constant baryon sequences of uniformly rotating strange star models, and are presented as a function of the stellar rotational frequency. The marginally stable orbit is present outside quark stars for a wide range of parameters, but outside the neutron stars it is present only for the highest values of mass. This allows a discrimination between quark stars and neutron stars in the resonance theory of kHz QPOs.

We perform some numerical study of the secular triaxial instability of rigidly rotating homogeneous fluid bodies in general relativity. In the Newtonian limit, this instability arises at the bifurcation point between the Maclaurin and Jacobi sequences. It can be driven in astrophysical systems by viscous dissipation. We locate the onset of instability along several constant baryon mass sequences of uniformly rotating axisymmetric bodies for compaction parameter $M/R = 0-0.275$. We find that general relativity weakens the Jacobi like bar mode instability, but the stabilizing effect is not very strong. According to our analysis the critical value of the ratio of the kinetic energy to the absolute value of the gravitational potential energy $(T/|W|)_{\rm crit}$ for compaction parameter as high as 0.275 is only 30% higher than the Newtonian value. The critical value of the eccentricity depends very weakly on the degree of relativity and for $M/R=0.275$ is only 2% larger than the Newtonian value at the onset for the secular bar mode instability. We compare our numerical results with recent analytical investigations based on the post-Newtonian expansion.

Observations of kHz quasi-periodic oscillations (QPOs) in the X-ray fluxes of low-mass X-ray binaries (LMXBs) have been used in attempts to constrain the external metric of the compact members of these binaries, as well as their masses and the equation of state of matter at supranuclear denisties. We compute the maximum orbital frequency of stable circular motion around uniformly rotating strange stars described by the MIT bag model. The calculations are performed for both normal and supramassive constant baryon mass sequences of strange stars rotating at all possible rates. We find the lower limits on the maximum orbital frequency and discuss them for a range of masses and for all rotational frequencies allowed in the model considered. We show that for slowly and moderately rotating strange stars the maximum value of orbital frequency can be a good indicator of the mass of the compact object. However, for rapidly rotating strange stars the same value of orbital frequency in the innermost stable circular orbit is obtained for stars with masses ranging from that of a planetoid to about three solar masses. At sufficiently high rotation rates of the strange star, the rotational period alone constrains the stellar mass to a surprisingly narrow range.

It has been suggested that observations of quasi periodic oscillations (QPOs) in the X-ray fluxes from low mass X-ray binaries can be used to constrain the mass of the compact object and the equation of state of its matter. A specific suggestion that the kHz QPO frequency saturates at the maximum orbital frequency has been widely considered. We examine rotating compact stars described by a new model of strange quark matter (Dey et al. 1998). We calculate the maximum orbital frequencies for both normal and supramassive strange stars described by the Dey model, and present these frequencies for sequences of equilibrium models with constant baryon mass. The maximum orbital frequencies for these compact objects are always higher than the kHz QPO frequencies observed to date.

We show that in the limit of low masses (M<0.1M solar), the orbital and rotational frequencies in exact relativistic numerical models of rotating, axially symmetric, quark stars coincide with those for Maclaurin spheroids. In particular, when the eccentricity of the spheroid exceeds a critical value (e>0.83458318), circular orbits in the equatorial plane are unstable for a range of orbital radii outside the stellar surface. The orbital frequency in the marginally stable orbit, and all other orbits, around Maclaurin spheroids goes to zero in the limit e=1. Maximum orbital frequencies in accreting Galactic X-ray sources may be a poor indication of the source mass.

We find constraints on minimum and maximum mass of ordinary neutron stars imposed by their early evolution (protoneutron star stage). We calculate models of protoneutron stars using a realistic standard equation of state of hot, dense matter valid for both supranuclear and subnuclear densities. Results for different values of the nuclear incompressibility are presented.

We study effects of the strange quark mass and of the QCD coupling constant in the MIT bag model of quark matter with lowest order quark-gluon interactions on the rapid rotation of strange stars. The influence of rotation on global parameters of strange stars is greater than in the case of the neutron star models. We show that independently of the strange quark mass and the coupling constans the ratio of the rotational kinetic energy to the absolute value of the gravitational potential energy for rotating strange stars is significantly higher than for ordinary neutron stars. This might indicate that rapidly rotating strange stars could be important sources of gravitational waves.

We compute numerical models of uniformly rotating strange stars (SS) in general relativity for the recently proposed QCD-based equation of state (EOS) of strange quark matter (Dey et al. 1998). Static models based on this EOS are characterised by a larger surface redshift than strange stars within the MIT bag model. The frequencies of the fastest rotating configurations described by Dey model are much higher than these for neutron stars (NS) and for the simplest SS MIT bag model. We determine a number of physical parameters for such stars and compare them with those obtained for NS. We construct constant baryon mass equilibrium sequences both normal and supramassive. Similarly to the NS a supramassive SS, prior to collapse to a black hole, spins up as it loses angular momentum. We find the upper limits on maximal masses and maximal frequencies of the rotating configurations. We show that the maximal rotating frequency for each of considered evolutionary sequences is never the Keplerian one. A normal and low mass supramassive strange stars gaining angular momentum always slows down just before reaching the Keplerian limit. For a high mass supramassive SS sequence the Keplerian configuration is the one with the lowest rotational frequency in the sequence. The value of $T/W$ for rapidly rotating SS of any mass is significantly higher than those for ordinary NS. For Keplerian configurations it increases as mass decreases. The results are robust for all linear self-bound equations of state.

We compute the maximum mass of moderately rotating strange stars as a function of the strange quark mass, of the QCD coupling constant, $\alpha_c$, and of the bag constant (vacuum energy density), $B$, in the MIT bag model of quark matter with lowest order quark-gluon interactions. For a fixed value of $B$, the maximum stellar mass depends only weakly on $\alpha_c$, and is independent of this coupling in the limit of massless quarks. However, if it is the value of the chemical potential of quark matter at zero pressure which is held constant, for example at the value corresponding to the stability limit of nucleons against conversion to quark matter, the maximum mass of the strange star is higher by up to 25% for $\alpha_c=0.6$, than for non-interacting quarks, and this may be relevant in the discussion of kHz QPO sources. The maximum mass of a non-rotating strange star could be sufficiently high to allow an orbital frequency as low as 1.0 kHz in the marginally stable orbit. However, for all $\alpha_c<0.6$, the stellar mass cannot exceed $2.6M_\odot$ at any rotational period $\ge1.6 $ms.

Exact calculations of innermost stable circular orbit (ISCO) around rotating strange stars are performed within the framework of general relativity. Equations of state (EOS) of strange quark matter based on the MIT Bag Model with massive strange quarks and lowest order QCD interactions, are used. The presence of a solid crust of normal matter on rotating, mass accreting strange stars in LMXBs is taken into account. It is found that, contrary to neutron stars, above some minimum mass a gap always separates the ISCO and stellar surface, independently of the strange star rotation rate. For a given baryon mass of strange star, we calculate the ISCO frequency as function of stellar rotation frequency, from static to Keplerian configuration. For masses close to the maximum mass of static configurations the ISCO frequencies for static and Keplerian configurations are similar. However, for masses significantly lower than the maximum mass of static configurations, the minimum value of the ISCO frequency is reached in the Keplerian limit. Presence of a solid crust increases the ISCO frequency for the Keplerian configuration by about ten percent compared to that for a bare strange star of the same mass.

The kilohertz quasi periodic oscillations (QPOs) discovered in several low mass X-ray binaries (LMXBs) by the Rossi X-ray Timing Explorer (XTE) are thought to occur at the orbital frequency in accretion discs whose inner edge corresponds to the innermost (marginally) stable orbit allowed by general relativity. These ideas have been applied to constrain the equation of state (e.o.s.) of the central neutron star. Here we discuss another possibility, that the central object is a strange star, and show how kHz QPOs constrain the e.o.s. of strange matter.