Downloads & Free Reading Options - Results

We present a general formalism to treat slowly rotating general relativistic superfluid neutron stars. As a first approximation, their matter content can be described in terms of a two-fluid model, where one fluid is the neutron superfluid, which is believed to exist in the core and inner crust of mature neutron stars, and the other fluid represents a conglomerate of all other constituents (crust nuclei, protons, electrons, etc.). We obtain a system of equations, good to second-order in the rotational velocities, that determines the metric and the matter variables, irrespective of the equation of state for the two fluids. In particular, allowance is made for the so-called entrainment effect, whereby the momentum of one constituent (e.g. the neutrons) carries along part of the mass of the other constituent. As an illustration of the developed framework, we consider a simplified equation of state for which the two fluids are described by different polytropes. We determine numerically the effects of the two fluids on the rotational frame-dragging, the induced changes in the neutron and proton densities and the inertial mass, as well as the change in shape of the star. We further discuss issues regarding conservation of the two baryon numbers, the mass-shedding (Kepler) limit and chemical equilibrium.

This article starts by providing an introductory overview of the theoretical mechanics of rotating neutron stars as developped to account for the frequency variations, and particularly the discontinuous glitches, observed in pulsars. The theory suggests, and the observations seem to confirm, that an essential role is played by the interaction between the solid crust and inner layers whose superfluid nature allows them to rotate independently. However many significant details remain to be clarified, even in much studied cases such as the Crab and Vela. The second part of this article is more technical, concentrating on just one of the many physical aspects that needs further development, namely the provision of a satisfactorily relativistic (local but not microscopic) treatment of the effects of the neutron superfluidity that is involved.

Causality alone suffices to set a lower bound on the period of rotation of relativistic stars as a function of their maximum observed mass. That is, by assuming a one-parameter equation of state (EOS) that satisfies v_sound < c and that allows stars with masses as large as the largest observed neutron-star mass, M_sph^max, we find P[ms] > 0.282 + 0.196 ( M_sph^max/M_odot-1.442). The limit does not assume that the EOS agrees with a known low-density form for ordinary matter, but if one adds that assumption, the minimum period is raised by a few percent. Thus the current minimum period of uniformly rotating stars, set by causality, is 0.28ms (0.29ms for stars with normal crust). The minimizing EOS yields models with a maximally soft exterior supported by a maximally stiff core. An analogous upper limit set by causality on the maximum mass of rotating neutron stars requires a low-density match and the limit depends on the matching density, epsilon_m. We recompute it, obtaining a slightly revised value, M_rot^max \sim 6.1( 2 * 10^14 g/cm^3 epsilon_m )^1/2 M_odot.

The influence of the rotation on the total masses and radii of the neutron stars are calculated by the Hartle's slow rotation formalism, while the equation of state is considered in a relativistic $\sigma-\omega$ model. Comparing with the observation, the calculating result shows that the double neutron star binaries are more like hyperon stars and the neutron stars of X-ray binaries are more like traditional neutron stars. As the changes of the mass and radius to a real neutron star caused by the rotation are very small comparing with the total mass and radius, one can see that Hartle's approximate method is rational to deal with the rotating neutron stars. If three property values: mass, radius and period are observed to the same neutron star, then the EOS of this neutron star could be decided entirely.

The influence of the rotation on the total masses and radii of the neutron stars are calculated by the Hartle's slow rotation formalism, while the equation of state is considered in a relativistic $\sigma-\omega$ model. Comparing with the observation, the calculating result shows that the double neutron star binaries are more like hyperon stars and the neutron stars of X-ray binaries are more like traditional neutron stars. As the changes of the mass and radius to a real neutron star caused by the rotation are very small comparing with the total mass and radius, one can see that Hartle's approximate method is rational to deal with the rotating neutron stars. If three property values: mass, radius and period are observed to the same neutron star, then the EOS of this neutron star could be decided entirely.

We present a new derivation of the equations governing the oscillations of slowly rotating relativistic stars. Previous investigations have been mostly carried out in the Regge-Wheeler gauge. However, in this gauge the process of linearizing the Einstein field equations leads to perturbation equations which as such cannot be used to perform numerical time evolutions. It is only through the tedious process of combining and rearranging the perturbation variables in a clever way that the system can be cast into a set of hyperbolic first order equations, which is then well suited for the numerical integration. The equations remain quite lengthy, and we therefore rederive the perturbation equations in a different gauge, which has been first proposed by Battiston et al. (1970). Using the ADM formalism, one is immediately lead to a first order hyperbolic evolution system, which is remarkably simple and can be numerically integrated without many further manipulations. Moreover, the symmetry between the polar and the axial equations becomes directly apparent.

The axial mode oscillations are examined for relativistic rotating stars with uniform angular velocity. Using the slow rotation formalism and the Cowling approximation, we have derived the equations governing the r-mode oscillations up to the second order with respect to the rotation. In the lowest order, the allowed range of the frequencies is determined, but corresponding spatial function is arbitrary. The spatial function can be decomposed in non-barotropic region by a set of functions associated with the differential equation of the second-order corrections. The equation however becomes singular in barotropic region, and a single function can be selected to describe the spatial perturbation of the lowest order. The frame dragging effect among the relativistic effects may be significant, as it results in rather broad spectrum of the r-mode frequency unlike in the Newtonian first-order calculation.

The axial mode oscillations are examined for relativistic rotating stars with uniform angular velocity. Using the slow rotation formalism and the Cowling approximation, we have derived the equations governing the r-mode oscillations up to the second order with respect to the rotation. In the lowest order, the allowed range of the frequencies is determined, but corresponding spatial function is arbitrary. The spatial function can be decomposed in non-barotropic region by a set of functions associated with the differential equation of the second-order corrections. The equation however becomes singular in barotropic region, and a single function can be selected to describe the spatial perturbation of the lowest order. The frame dragging effect among the relativistic effects may be significant, as it results in rather broad spectrum of the r-mode frequency unlike in the Newtonian first-order calculation.

Non-axisymmetric oscillations of differentially rotating stars are studied using both slow rotation and Cowling approximation. The equilibrium stellar models are relativistic polytropes where differential rotation is described by the relativistic j-constant rotation law. The oscillation spectrum is studied versus three main parameters: the stellar compactness $M/R$, the degree of differential rotation $A$ and the number of maximun couplings $\ell_{\rm max}$. It is shown that the rotational splitting of the non-axisymmetric modes is strongly enhached by increasing the compactness of the star and the degree of differential rotation. Finally, we investigate the relation between the fundamental quadrupole mode and the corotation band of differentially rotating stars.

We investigate the innermost stable circular orbit (ISCO) of a test particle moving on the equatorial plane around rotating relativistic stars such as neutron stars. First, we derive approximate analytic formulas for the angular velocity and circumferential radius at the ISCO making use of an approximate relativistic solution which is characterized by arbitrary mass, spin, mass quadrupole, current octapole and mass $2^4$-pole moments. Then, we show that the analytic formulas are accurate enough by comparing them with numerical results, which are obtained by analyzing the vacuum exterior around numerically computed geometries for rotating stars of polytropic equation of state. We demonstrate that contribution of mass quadrupole moment for determining the angular velocity and, in particular, the circumferential radius at the ISCO around a rapidly rotating star is as important as that of spin.

Whether analytic exact vacuum(electrovacuum) solutions of the Einstein(Einstein-Maxwell) field equations can accurately describe or not the exterior spacetime of compact stars remains still an interesting open question in Relativistic Astrophysics. As an attempt to establish their level of accuracy, the radii of the Innermost Stable Circular Orbits (ISCOs) of test particles given by analytic exterior spacetime geometries have been compared with the ones given by numerical solutions for neutron stars (NSs) obeying a realistic equation of state (EoS). It has been so shown that the six-parametric solution of Pach\'on, Rueda, and Sanabria (2006) (hereafter PRS) is more accurate to describe the NS ISCO radii than other analytic models. We propose here an additional test of accuracy for analytic exterior geometries based on the comparison of orbital frequencies of neutral test particles. We compute the Keplerian, frame-dragging, as well as the precession and oscillation frequencies of the radial and vertical motions of neutral test particles for the Kerr and PRS geometries; then we compare them with the numerical values obtained by Morsink and Stella (1999) for realistic NSs. We identify the role of high-order multipole moments such as the mass quadrupole and current octupole in the determination of the orbital frequencies especially in the rapid rotation regime. The results of this work are relevant to cast a separatrix between black hole (BH) and NS signatures as well as probe the nuclear matter EoS and NS parameters from the Quasi-Periodic Oscillations (QPOs) observed in Low Mass X-Ray Binaries.

Non-axisymmetric oscillations of rapidly rotating relativistic stars are studied using the Cowling approximation. The oscillation spectra have been estimated by Fourier transforming the evolution equations describing the perturbations. This is the first study of its kind and provides information on the effect of fast rotation on the oscillation spectra while it offers the possibility in studying the complete problem by including spacetime perturbations. Our study includes both axisymmetric and non-axisymmetric perturbations and provides limits for the onset of the secular bar mode rotational instability. We also present approximate formulae for the dependence of the oscillation spectrum from rotation. The results suggest that it is possible to extract the relativistic star's parameters from the observed gravitational wave spectrum.

We solve for the time-dependent dynamics of axisymmetric, general relativistic magnetohydrodynamic winds from rotating neutron stars. The mass loss rate is obtained self-consistently as a solution to the MHD equations, subject to a finite thermal pressure at the stellar surface. We consider both monopole and dipole magnetic field geometries and we explore the parameter regime extending from low magnetization (low-sigma_o), almost thermally-driven winds to high magnetization (high-sigma_o), relativistic Poynting-flux dominated outflows. We compute the angular momentum and rotational energy loss rates as a function of sigma_o and compare with analytic expectations from the classical theory of pulsars and magnetized stellar winds. In the case of the monopole, our high-sigma_o calculations asymptotically approach the analytic force-free limit. If we define the spindown rate in terms of the open magnetic flux, we similarly reproduce the spindown rate from recent force-free calculations of the aligned dipole. However, even for sigma_o as high as ~20, we find that the location of the Y-type point (r_Y), which specifies the radius of the last closed field line in the equatorial plane, is not the radius of the light cylinder R_L = c/omega (R = cylindrical radius), as has previously been assumed in most estimates and force-free calculations. Instead, although the Alfven radius at intermediate latitudes quickly approaches R_L as sigma_o exceeds unity, r_Y remains significantly less than R_L. Because r_Y < R_L, our calculated spindown rates thus exceed the classic ``vacuum dipole'' rate. We discussthe implications of our results for models of rotation-powered pulsars and magnetars, both in their observed states and in their hypothesized rapidly rotating initial state.

We study the magnetosphere of a slowly rotating magnetized neutron star subject to toroidal oscillations in the relativistic regime. Under the assumption of a zero inclination angle between the magnetic moment and the angular momentum of the star, we analyze the Goldreich-Julian charge density and derive a second-order differential equation for the electrostatic potential. The analytical solution of this equation in the polar cap region of the magnetosphere shows the modification induced by stellar toroidal oscillations on the accelerating electric field and on the charge density. We also find that, after decomposing the oscillation velocity in terms of spherical harmonics, the first few modes with $m=0,1$ are responsible for energy losses that are almost linearly dependent on the amplitude of the oscillation and that, for the mode $(l,m)=(2,1)$, can be a factor $\sim8$ larger than the rotational energy losses, even for a velocity oscillation amplitude at the star surface as small as $\eta=0.05 \ \Omega \ R$. The results obtained in this paper clarify the extent to which stellar oscillations are reflected in the time variation of the physical properties at the surface of the rotating neutron star, mainly by showing the existence of a relation between $P\dot{P}$ and the oscillation amplitude. Finally, we propose a qualitative model for the explanation of the phenomenology of intermittent pulsars in terms of stellar oscillations that are periodically excited by star glitches.

We study the final state of the gravitational collapse of uniformly rotating supramassive neutron stars by axisymmetric simulations in full general relativity. The rotating stars provided as the initial condition are marginally stable against quasiradial gravitational collapse and its equatorial radius rotates with the Kepler velocity (i.e., the star is at the mass-shedding limit). To model the neutron stars, we adopt the polytropic equations of state for a wide range of the polytropic index as $n=2/3$, 4/5, 1, 3/2 and 2. We follow the formation and evolution of the black holes, and show that irrespective of the value of $n (2/3\leq n \leq 2)$, the final state is a Kerr black hole and the disk mass is very small ($ < 10^{-3}$ of the initial stellar mass).

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.

The conversion of neutron star to strange star is argued to be a two step process. The first process involves the deconfinement of nuclear to two-flavour quark matter. The GR results shows that the propagating front breaks up into fragments which propagate with different velocities along different directions. The time taken for this conversion to happen is of the order of few $ms$. This calculation indicates the inadequacy of non-relativistic (NR) or even Special Relativistic (SR) treatments for these cases.

We investigate the viscosity driven instability in rotating relativistic stars by means of an iterative approach. We focus on polytropic rotating equilibrium stars and impose an m=2 perturbationin the lapse. We vary both the stiffness of the equation of state and the compactness of the star to study those effects on the value of the threshold. For a uniformly rotating star, the criterion T/W, where T is the rotational kinetic energy and W is the gravitational binding energy, mainly depends on the compactness of the star and takes values around 0.13 ~ 0.16, which differ slightly from that of Newtonian incompressible stars (~ 0.14). For differentially rotating stars, the critical value of T/W is found to span the range 0.17 - 0.25. This is significantly larger than the uniformly rotating case with the same compactness of the star. Finally we discuss a possibility of detecting gravitational waves from viscosity driven instability with ground-based interferometers.

Neutrino asymmetry in general relativistic radiative spacetime exterior to spinning stars is investigating by making use of Newmann-Penrose (NP) spin coefficient formalism. It is shown that neutrino current depends on the direction of rotation of the star. The solution is obtained in test field approximation where the neutrinos do not generate gravitational fields.

We study non-linear axisymmetric pulsations of rotating relativistic stars using a general relativistic hydrodynamics code under the assumption of a conformal flatness. We compare our results to previous simulations where the spacetime dynamics was neglected. The pulsations are studied along various sequences of both uniformly and differentially rotating relativistic polytropes with index N = 1. We identify several modes, including the lowest-order l = 0, 2, and 4 axisymmetric modes, as well as several axisymmetric inertial modes. Differential rotation significantly lowers mode frequencies, increasing prospects for detection by current gravitational wave interferometers. We observe an extended avoided crossing between the l = 0 and l = 4 first overtones, which is important for correctly identifying mode frequencies in case of detection. For uniformly rotating stars near the mass-shedding limit, we confirm the existence of the mass-shedding-induced damping of pulsations, though the effect is not as strong as in the Cowling approximation. We also investigate non-linear harmonics of the linear modes and notice that rotation changes the pulsation frequencies in a way that would allow for various parametric instabilities between two or three modes to take place. We assess the detectability of each obtained mode by current gravitational wave detectors and outline how the empirical relations that have been constructed for gravitational wave asteroseismology could be extended to include the effects of rotation.

We present the first direct comparison of codes based on two different numerical methods for constructing rapidly rotating relativistic stars. A code based on the Komatsu-Eriguchi-Hachisu (KEH) method (Komatsu et al. 1989), written by Stergioulas, is compared to the Butterworth-Ipser code (BI), as modified by Friedman, Ipser and Parker. We compare models obtained by each method and evaluate the accuracy and efficiency of the two codes. The agreement is surprisingly good. A relatively large discrepancy recently reported (Eriguchi et al. 1994) is found to arise from the use of two different versions of the equation of state. We find, for a given equation of state, that equilibrium models with maximum values of mass, baryon mass, and angular momentum are (generically) all distinct and either all unstable to collapse or are all stable. Our implementation of the KEH method will be available as a public domain program for interested users.

We present long-term (~10^4 M) axisymmetric simulations of differentially rotating, magnetized neutron stars in the slow-rotation, weak magnetic field limit using a perturbative metric evolution technique. Although this approach yields results comparable to those obtained via nonperturbative (BSSN) evolution techniques, simulations performed with the perturbative metric solver require about 1/4 the computational resources at a given resolution. This computational efficiency enables us to observe and analyze the effects of magnetic braking and the magnetorotational instability (MRI) at very high resolution. Our simulations demonstrate that (1) MRI is not observed unless the fastest-growing mode wavelength is resolved by more than about 10 gridpoints; (2) as resolution is improved, the MRI growth rate converges, but due to the small-scale turbulent nature of MRI, the maximum growth amplitude increases, but does not exhibit convergence, even at the highest resolution; and (3) independent of resolution, magnetic braking drives the star toward uniform rotation as energy is sapped from differential rotation by winding magnetic fields.

The r-mode instability in rotating relativistic stars has been shown recently to have important astrophysical implications (including the emission of detectable gravitational radiation, the explanation of the initial spins of young neutron stars and the spin-distribution of millisecond pulsars and the explanation of one type of gamma-ray bursts), provided that r-modes are not saturated at low amplitudes by nonlinear effects or by dissipative mechanisms. Here, we present the first study of nonlinear r-modes in isentropic, rapidly rotating relativistic stars, via 3-D general-relativistic hydrodynamical evolutions. Our numerical simulations show that (1) on dynamical timescales, there is no strong nonlinear coupling of r-modes to other modes at amplitudes of order one -- unless nonlinear saturation occurs on longer timescales, the maximum r-mode amplitude is of order unity (i.e., the velocity perturbation is of the same order as the rotational velocity at the equator). An absolute upper limit on the amplitude (relevant, perhaps, for the most rapidly rotating stars) is set by causality. (2) r-modes and inertial modes in isentropic stars are predominantly discrete modes and possible associated continuous parts were not identified in our simulations. (3) In addition, the kinematical drift associated with r-modes, recently found by Rezzolla, Lamb and Shapiro (2000), appears to be present in our simulations, but an unambiguous confirmation requires more precise initial data. We discuss the implications of our findings for the detectability of gravitational waves from the r-mode instability.

In this paper we compute general-relativistic polytropic models simulating rigidly rotating, pulsating neutron stars. These relativistic compact objects, with a radius of $\sim 10 \, \mathrm{km}$ and mass between $\sim 1.4$ and $3.2$ solar masses, are closely related to pulsars. We emphasize on computing the change in the pulsation eigenfrequencies owing to a rigid rotation, which, in turn, is a decisive issue for studying stability of such objects. In our computations, we keep rotational perturbation terms of up to second order in the angular velocity.

The paper considers the spectrum of axial perturbations of slowly uniformly rotating general relativistic stars in the framework of Y. Kojima. In a first step towards a full analysis only the evolution equations are treated but not the constraint. Then it is found that the system is unstable due to a continuum of non real eigenvalues. In addition the resolvent of the associated generator of time evolution is found to have a special structure which was discussed in a previous paper. From this structure it follows the occurrence of a continuous part in the spectrum of oscillations at least if the system is restricted to a finite space as is done in most numerical investigations. Finally, it can be seen that higher order corrections in the rotation frequency can qualitatively influence the spectrum of the oscillations. As a consequence different descriptions of the star which are equivalent to first order could lead to different results with respect to the stability of the star.

We study the effects of rotation on the torsional modes of oscillating relativistic stars with a solid crust. Earlier works in Newtonian theory provided estimates of the rotational corrections for the torsional modes and suggested that they should become CFS unstable, even for quite low rotation rates. In this work, we study the effect of rotation in the context of general relativity using elasticity theory and in the slow-rotation approximation. We find that the Newtonian picture does not change considerably. The inclusion of relativistic effects leads only to quantitative corrections. The degeneracy of modes for different values of $m$ is removed, and modes with $\ell=m$ are shifted towards zero frequencies and become secularly unstable at stellar rotational frequencies $\sim$ 20-30 Hz.

A numerical study of nonlinear r-modes in isentropic, rapidly rotating relativistic stars, via 3-D general-relativistic hydrodynamical evolutions, is presented. On dynamical timescales, we find no evidence for strong coupling of r-modes to other modes at amplitudes of order one or larger. Therefore, unless nonlinear saturation sets in on longer timescales, the maximum r-mode amplitude is of order one. An absolute upper limit on the amplitude is set by causality. Our simulations also show that r-modes and inertial modes in isentropic stars are discrete modes, with no evidence for the existence of a continuous part in the spectrum.

We present a new numerical code to compute non-axisymmetric eigenmodes of rapidly rotating relativistic stars by adopting spatially conformally flat approximation of general relativity. The approximation suppresses the radiative degree of freedom of relativistic gravity and the field equations are cast into a set of elliptic equations. The code is tested against the low-order f- and p-modes of slowly rotating stars for which a good agreement is observed in frequencies computed by our new code and those computed by the full theory. Entire sequences of the low order counter-rotating f-modes are computed, which are susceptible to an instability driven by gravitational radiation.

We study the magnetic braking and viscous damping of differential rotation in incompressible, uniform density stars in general relativity. Differentially rotating stars can support significantly more mass in equilibrium than nonrotating or uniformly rotating stars. The remnant of a binary neutron star merger or supernova core collapse may produce such a "hypermassive" neutron star. Although a hypermassive neutron star may be stable on a dynamical timescale, magnetic braking and viscous damping of differential rotation will ultimately alter the equilibrium structure, possibly leading to delayed catastrophic collapse. Here we consider the slow-rotation, weak-magnetic field limit in which E_rot < < E_mag < < W, where E_rot is the rotational kinetic energy, E_mag is the magnetic energy, and W is the gravitational binding energy of the star. We assume the system to be axisymmetric and solve the MHD equations in both Newtonian gravitation and general relativity. Toroidal magnetic fields are generated whenever the angular velocity varies along the initial poloidal field lines. We find that the toroidal fields and angular velocities oscillate independently along each poloidal field line, which enables us to transform the original 2+1 equations into 1+1 form and solve them along each field line independently. The incoherent oscillations on different field lines stir up turbulent-like motion in tens of Alfven timescales ("phase mixing"). In the presence of viscosity, the stars eventually are driven to uniform rotation, with the energy contained in the initial differential rotation going into heat. Our evolution calculations serve as qualitative guides and benchmarks for future, more realistic MHD simulations in full 3+1 general relativity.

We consider the slowly rotating relativistic stars with a uniform angular velocity in the scalar-tensor gravity, and examine the rotational effect around such compact objects. For this purpose, we derive a 2nd order differential equation describing the frame dragging in the scalar-tensor gravity and solve it numerically. As a result, we find that the total angular momentum is proportional to the angular velocity even in the scalar-tensor gravity. We also show that one can observe the spontaneous scalarization in rotational effects as well as the other stellar properties, if the cosmological value of scalar field is zero. On the other hand, if the cosmological value of scalar field is nonzero, the deviation from the general relativity can be seen in a wide range of the coupling constant. Additionally, we find that the deviation from the general relativity becomes larger with more massive stellar models, which is independent of the cosmological value of scalar field. Thus, via precise observations of astronomical phenomena associated with rotating relativistic stars, one may be possible to probe not only the gravitational theory in the strong-field regime, but also the existence of scalar field.

We study oscillations of slowly rotating relativistic barotropic as well as non-barotropic polytropic stars in the Cowling approximation, including first order rotational corrections. By taking into account the coupling between the polar and axial equations, we find that, in contrast to previous results, the $m=2$ $r$ modes are essentially unaffected by the continuous spectrum and exist even for very relativistic stellar models. We perform our calculations both in the time and frequency domain. In order to numerically solve the infinite system of coupled equations, we truncate it at some value $l_{\rm max}$. Although the time dependent equations can be numerically evolved without any problems, the eigenvalue equations possess a singular structure, which is related to the existence of a continuous spectrum. This prevents the numerical computation of an eigenmode if its eigenfrequency falls inside the continuous spectrum. The properties of the latter depend strongly on the cut-off value $l_{\rm max}$ and it can consist of several either disconnected or overlapping patches, which are the broader the more relativistic the stellar model is. By discussing the dependence of the continuous spectrum as a function of both the cut-off value $l_{\rm max}$ and the compactness $M/R$, we demonstrate how it affects the inertial modes. Through the time evolutions we are able to show that some of the inertial modes can actually exist inside the continuous spectrum, but some cannot. For more compact and therefore more relativistic stellar models, the width of the continuous spectrum strongly increases and as a consequence, some of the inertial modes, which exist in less relativistic stars, disappear.

For relativistic stars rotating slowly and differentially with a positive angular velocity, some properties in relation to the positiveness of the rate of rotational dragging and of the angular momentum density are derived. Also, a new proof for the bounds on the rotational mass-energy is given.

Neutron stars that are cold enough should have two or more superfluids/supercondutors in their inner crusts and cores. The implication of superfluidity/superconductivity for equilibrium and dynamical neutron star states is that each individual particle species that forms a condensate must have its own, independent number density current and equation of motion that determines that current. An important consequence of the quasiparticle nature of each condensate is the so-called entrainment effect, i.e. the momentum of a condensate is a linear combination of its own current and those of the other condensates. We present here the first fully relativistic modelling of slowly rotating superfluid neutron stars with entrainment that is accurate to the second-order in the rotation rates. The stars consist of superfluid neutrons, superconducting protons, and a highly degenerate, relativistic gas of electrons. We use a relativistic $\sigma$ - $\omega$ mean field model for the equation of state of the matter and the entrainment. We determine the effect of a relative rotation between the neutrons and protons on a star's total mass, shape, and Kepler, mass-shedding limit.

Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations and from first integrals and integrability conditions of the magnetohydrodynamic equilibrium equations. After a brief introduction of the formulation of the problem, we present the first results for highly deformed magnetized rotating compact stars.

We review here recent numerical results concerning the acceleration and collimation of relativistic outflows from neutron stars, and we discuss their implications for models of magnetars in their hypothesized initial high rotation states. New results with different injection conditions and extending to much larger distance from the central star are also presented and discussed, showing that, only in the case of weakly magnetized winds, collimated outflows are indeed possible. We finally comment on the possibility that newly born magnetars might trigger GRBs, and we show that numerical results do not support such hypothesis.

We study nonaxisymmetric perturbations of rotating relativistic stars. modeled as perfect-fluid equilibria. Instability to a mode with angular dependence $\exp(im\phi)$ sets in when the frequency of the mode vanishes. The locations of these zero-frequency modes along sequences of rotating stars are computed in the framework of general relativity. We consider models of uniformly rotating stars with polytropic equations of state, finding that the relativistic models are unstable to nonaxisymmetric modes at significantly smaller values of rotation than in the Newtonian limit. Most strikingly, the m=2 bar mode can become unstable even for soft polytropes of index $N \leq 1.3$, while in Newtonian theory it becomes unstable only for stiff polytropes of index $N \leq 0.808$. If rapidly rotating neutron stars are formed by the accretion-induced collapse of white dwarfs, instability associated with these nonaxisymmetric, gravitational-wave driven modes may set an upper limit on neutron-star rotation. Consideration is restricted to perturbations that correspond to polar perturbations of a spherical star. A study of axial perturbations is in progress.

During a Type-I burst, the turbulent deflagation front may excite waves in the neutron star ocean and upper atmosphere with frequencies, $\omega \sim 1$ Hz. These waves may be observed as highly coherent flux oscillations during the burst. The frequencies of these waves changes as the upper layers of the neutron star cool which accounts for the small variation in the observed QPO frequencies. In principle several modes could be excited but the fundamental buoyant $r-$mode exhibits significantly larger variability for a given excitation than all of the other modes. An analysis of modes in the burning layers themselves and the underlying ocean shows that it is unlikely these modes can account for the observed burst oscillations. On the other hand, photospheric modes which reside in a cooler portion of the neutron star atmosphere may provide an excellent explanation for the observed oscillations.

During a Type-I burst, the turbulent deflagation front may excite waves in the neutron star ocean and upper atmosphere with frequencies, $\omega \sim 1$ Hz. These waves may be observed as highly coherent flux oscillations during the burst. The frequencies of these waves changes as the upper layers of the neutron star cool which accounts for the small variation in the observed QPO frequencies. In principle several modes could be excited but the fundamental buoyant $r-$mode exhibits significantly larger variability for a given excitation than all of the other modes. An analysis of modes in the burning layers themselves and the underlying ocean shows that it is unlikely these modes can account for the observed burst oscillations. On the other hand, photospheric modes which reside in a cooler portion of the neutron star atmosphere may provide an excellent explanation for the observed oscillations.

Axial oscillations relevant to the r-mode instability are studied with slow rotation formalism in general relativity. The approximate equation governing the oscillations is derived with second-order rotational corrections. The equation contains an effective 'viscosity-like' term, which originates from coupling to the polar g-mode displacements. The term plays a crucial role on the resonance point, where the disturbance on the rotating stars satisfies a certain condition at the lowest order equation. The effect is significant for newly born hot neutron stars, which are expected to be subject to the gravitational radiation driven instability of the r-mode.

In order to examine the rotational effect around neutron star in tensor-vector-scalar (TeVeS) theory, we consider the slowly rotating relativistic stars with a uniform angular velocity. As a result, we find that similar to the case in general relativity (GR), the angular momentum is proportional to the angular velocity. Additionally, as the value of coupling constant $K$ becomes higher, the frame dragging in TeVeS becomes quite different distribution from that in GR, where we can also see the deviation even in the interior of star. While with smaller value of $K$, although the frame dragging approaches to that expected in GR, the induced vector field due to the rotation does not vanish and still exists. Thus, through the observations associated with relativistic object, one could be possible to distinguish the gravitational theory in strong field regime even in the case that the value of coupling constant $K$ is quite small.

We investigate properties of r-modes characterized by regular eigenvalue problem in slowly rotating relativistic polytropes. Our numerical results suggest that discrete r-mode solutions for the regular eigenvalue problem exist only for restricted polytropic models. In particular the r-mode associated with l=m=2, which is considered to be the most important for gravitational radiation driven instability, do not have a discrete mode as solutions of the regular eigenvalue problem for polytropes having the polytropic index N > 1.18 even in the post-Newtonian order. Furthermore for a N=1 polytrope, which is employed as a typical neutron star model, discrete r-mode solutions for regular eigenvalue problem do not exist for stars whose relativistic factor M/R is larger than about 0.1. Here M and R are stellar mass and stellar radius, respectively.

The frequencies of r-mode oscillations of rotating neutron stars can be useful for guiding and interpreting gravitational wave and electromagnetic observations. The frequencies of slowly rotating, barotropic, and non-magnetic Newtonian stars are well known, but subject to various corrections. After making simple estimates of the relative strengths of these corrections we conclude that relativistic corrections are likely to be the most important. For this reason we extend the formalism of K. H. Lockitch, J. L. Friedman, and N. Andersson [Phys. Rev. D 68, 124010 (2003)], who consider relativistic polytropes, to the case of realistic equations of state. This formulation results in perturbation equations which are solved using a spectral method. We find that for realistic equations of state the r-mode frequency ranges from 1.39 to 1.57 times the spin frequency of the star when the relativistic compactness parameter (M/R) is varied over the astrophysically motivated interval 0.11 to 0.31. The results presented here are relevant to the design of gravitational wave and electromagnetic r-mode searches, and following a successful r-mode detection could help constrain the high density equation of state.

R-modes of neutron stars could be a source of gravitational waves for ground based detectors. If the precise frequency $\sigma$ is known, guided gravitational wave searches with enhanced detectability are possible. Because of its physical importance many authors have calculated the r-mode frequency. For the dominant mode, the associated gravitational wave frequency is 4/3 times the angular velocity of the star $\Omega$, subject to various corrections of which relativistic and rotational corrections are the most important. This has led several authors to investigate the dependence of the r-mode frequency on factors such as the relativistic compactness parameter ($M/R$) and the angular velocity of stars with different equations of state. The results found so far, however, are almost independent of the equation of state. Here we investigate the effect of rapid rotation and differential rotation on $\sigma$. We evolve the perturbation equations using the Cowling approximation by applying finite differencing methods to compute the r-mode frequency for a series of rotating neutron stars with polytropic equations of state. We find that rotational effects in the r-mode frequency can be larger than relativistic effects for rapidly spinning stars with low compactness, reducing the observed frequency. Differential rotation also acts to decrease $\sigma$, but its effect inscreases with the compactness. The results presented here are relevant to the design of gravitational wave and electromagnetic r-mode searches.

Quasi-toroidal oscillations in slowly rotating stars are examined in the framework of general relativity. The oscillation frequency to first order of the rotation rate is not a single value even for uniform rotation unlike the Newtonian case. All the oscillation frequencies of the r-modes are purely neutral and form a continuous spectrum limited to a certain range. The allowed frequencies are determined by the resonance condition between the perturbation and background mean flow. The resonant frequency varies with the radius according to general relativistic dragging effect.

These lecture notes are intended to introduce the theory of rotating stars in general relativity. The focus is put on the theoretical foundations, with a detailed discussion of the spacetime symmetries, the choice of coordinates and the derivation of the equations of structure from the Einstein equation. The global properties of rotating stars (mass, angular momentum, redshifts, orbits, etc.) are also introduced.

Two pseudo-Newtonian potentials, which approximate the angular and epicyclic frequencies of the relativistic accretion disk around rotating (and counter rotating) compact objects, are presented. One of them, the Logarithmically Modified Potential, is a better approximation for the frequencies while the other, the Second-order Expanded potential, also reproduces the specific energy for circular orbits in close agreement with the General Relativistic values. These potentials may be included in time dependent hydrodynamical simulations to study the temporal behavior of such accretion disks.

We develop a numerical scheme for obtaining the $r$-mode oscillations of rapidly rotating relativistic stars. In the present scheme, we neglect all metric perturbations and only take account of the dynamics of the fluid in the background spacetime of the unperturbed star (the relativistic Cowling approximation). We also assume the star is barotropic, i.e. neutrally stable against convection under the assumption of adiabatic oscillations. Our numerical scheme is based on the Yoshida-Eriguchi formulation for the analysis of the general relativistic f-mode oscillations in the Cowling approximation and a general relativistic generalization of the Karino-Yoshida-Yoshida-Eriguchi's numerical scheme for obtaining oscillations of rapidly rotating Newtonian stars. By this new numerical scheme, frequencies of the $r$-mode oscillations are obtained as functions of the ratio of the rotational energy to the absolute value of the gravitational energy $T/|W|$ along sequences of polytropic equilibrium stars whose ratio of the pressure to the total energy density at the center of the star and the polytropic index are kept constant. It is found that the dimensionless oscillation frequency $\sigma/\Omega$ is a linearly decreasing function of $T/|W|$, where $\sigma$ and $\Omega$ are the oscillation frequency and the angular velocity of the star measured in an inertial frame at spatial infinity. We also find that oscillation frequencies of the $r$-modes are highly dependent on the relativistic factor $M/R$ of the star as already found in previous studies in which the slow rotation approximation has been used. Here $M$ and $R$ denote the mass and radius of the star, respectively.

We compute zero-frequency (neutral) quasi-normal f-modes of fully relativistic and rapidly rotating neutron stars, using several realistic equations of state (EOSs) for neutron star matter. The zero-frequency modes signal the onset of the gravitational radiation-driven instability. We find that the l=m=2 (bar) f-mode is unstable for stars with gravitational mass as low as 1.0 - 1.2 M_\odot, depending on the EOS. For 1.4 M_\odot neutron stars, the bar mode becomes unstable at 83 % - 93 % of the maximum allowed rotation rate. For a wide range of EOSs, the bar mode becomes unstable at a ratio of rotational to gravitational energies T/W \sim 0.07-0.09 for 1.4 M_\odot stars and T/W \sim 0.06 for maximum mass stars. This is to be contrasted with the Newtonian value of T/W \sim 0.14. We construct the following empirical formula for the critical value of T/W for the bar mode, (T/W)_2 = 0.115 - 0.048 M / M_{max}^{sph}, which is insensitive to the EOS to within 4 - 6 %. This formula yields an estimate for the neutral mode sequence of the bar mode as a function only of the star's mass, M, given the maximum allowed mass, M_{max}^{sph}, of a nonrotating neutron star. The recent discovery of the fast millisecond pulsar in the supernova remnant N157B, supports the suggestion that a fraction of proto-neutron stars are born in a supernova collapse with very large initial angular momentum. Thus, in a fraction of newly born neutron stars the instability is a promising source of continuous gravitational waves. It could also play a major role in the rotational evolution (through the emission of angular momentum) of merged binary neutron stars, if their post-merger angular momentum exceeds the maximum allowed to form a Kerr black hole.

We investigate the properties of relativistic $r$-modes of slowly rotating neutron stars by using a relativistic version of the Cowling approximation. In our formalism, we take into account the influence of the Coriolis like force on the stellar oscillations, but ignore the effects of the centrifugal like force. For three neutron star models, we calculated the fundamental $r$-modes with $l'=m=2$ and 3. We found that the oscillation frequency $\bar\sigma$ of the fundamental $r$-mode is in a good approximation given by $\bar\sigma\approx \kappa_0 \Omega$, where $\bar\sigma$ is defined in the corotating frame at the spatial infinity, and $\Omega$ is the angular frequency of rotation of the star. The proportional coefficient $\kappa_0$ is only weakly dependent on $\Omega$, but it strongly depends on the relativistic parameter $GM/c^2R$, where $M$ and $R$ are the mass and the radius of the star. All the fundamental $r$-modes with $l'=m$ computed in this study are discrete modes with distinct regular eigenfunctions, and they all fall in the continuous part of the frequency spectrum associated with Kojima's equation (Kojima 1998). These relativistic $r$-modes are obtained by including the effects of rotation higher than the first order of $\Omega$ so that the buoyant force plays a role, the situation of which is quite similar to that for the Newtonian $r$-modes.

We present a self-consistent, relativistic model of rapidly rotating neutron stars describing their exterior gravitational field. This is achieved by matching the new solution of Einstein's field equations found by Manko et al. (2000) and the numerical results for the interior of neutron stars with different equations of state calculated by Cook et al. (1994). This matching process gives constraints for the choice of the five parameters of the vacuum solution. Then we investigate some properties of the gravitational field of rapidly rotating neutron stars with these fixed parameters.