We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030$+$0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer (NICER) X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply-connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass $M$ and equatorial radius $R_\mathrm{eq}$ are, respectively, $1.34_{-0.16}^{+0.15}$ M$_{\odot}$ and $12.71_{-1.19}^{+1.14}$ km, whilst the compactness $GM/R_\mathrm{eq}c^2 = 0.156_{-0.010}^{+0.008}$ is more tightly constrained; the credible interval bounds reported here are approximately the $16\%$ and $84\%$ quantiles in marginal posterior mass.

Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mass $M$. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14}\ M_\odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density.

We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes.

Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EOS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EOS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EOS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass-radius likelihood function derived by Riley et al. (2019) using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451.

The global climate crisis poses new risks to humanity, and with them, new challenges to the practices of professional astronomy. Avoiding the more catastrophic consequences of global warming by more than 1.5 degrees requires an immediate reduction of greenhouse gas emissions. According to the 2018 United Nations Intergovernmental Panel report, this will necessitate a 45% reduction of emissions by 2030 and net-zero emissions by 2050. Efforts are required at all levels, from the individual to the governmental, and every discipline must find ways to achieve these goals. This will be especially difficult for astronomy with its significant reliance on conference and research travel, among other impacts. However, our long-range planning exercises provide the means to coordinate our response on a variety of levels. We have the opportunity to lead by example, rising to the challenge rather than reacting to external constraints. We explore how astronomy can meet the challenge of a changing climate in clear and responsible ways, such as how we set expectations (for ourselves, our institutions, and our granting agencies) around scientific travel, the organization of conferences, and the design of our infrastructure. We also emphasize our role as reliable communicators of scientific information on a problem that is both human and planetary in scale.

The advent of moderately high-resolution X-ray spectroscopy with Chandra and XMM promised to usher in a new age in the study of neutron stars: we thought we would study neutron stars like stars, with resolved absorption spectra revealing their surface chemical composition and physical conditions (e.g. surface gravity, pressure, temperature). Nature, however, did not cooperate in this endeavor, as observations of neutron stars have not revealed verified atomic absorption lines yet. In the near future, advancements in transition-edge sensors (TES) technology will allow for electron-volt-resolution spectroscopy combined with nanoseconds-precision timing. Combining these detectors with collector optics will also us to study neutron stars in much greater detail by achieving high-energy resolution with much larger collecting areas to uncover even weak spectral features over a wide range of the photon energies. Perhaps we will finally be able to study neutron stars like stars.

The X-ray emission of neutron stars and black holes presents a rich phenomenology that can lead us to a better understanding of their nature and to address more general physics questions: Does general relativity apply in the strong gravity regime? Is spacetime around black holes described by the Kerr metric? This white paper considers how we can investigate these questions by studying reverberation mapping and quasi-periodic oscillations in accreting systems with a combination of high-spectral and high-timing resolution. In the near future, we will be able to study compact objects in the X-rays in a new way: advancements in transition-edge sensors (TES) technology will allow for electron-volt-resolution spectroscopy combined with nanoseconds-precision timing.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

We investigate the highest-frequency kHz QPOs previously detected with RXTE in six neutron star (NS) low mass X-ray binaries. We find that the highest-frequency kHz QPO detected in 4U 0614+09 has a 1267 Hz 3$\sigma$ confidence lower limit on its centroid frequency. This is the highest such limit reported to date, and of direct physical interest as it can be used to constrain QPO models and the supranuclear density equation of state (EoS). We compare our measured frequencies to maximum orbital frequencies predicted in full GR using models of rotating neutron stars with a number of different modern EoS and show that these can accommodate the observed QPO frequencies. Orbital motion constrained by NS and ISCO radii is therefore a viable explanation of these QPOs. In the most constraining case of 4U 0614+09 we find the NS mass must be M$<$2.1 M$_{\odot}$. From our measured QPO frequencies we can constrain the NS radii for five of the six sources we studied to narrow ranges ($\pm$0.1-0.7 km) different for each source and each EoS.

The equation of state of cold supra-nuclear-density matter, such as in neutron stars, is an open question in astrophysics. A promising method for constraining the neutron star equation of state is modelling pulse profiles of thermonuclear X-ray burst oscillations from hotspots on accreting neutron stars. The pulse profiles, constructed using spherical and oblate neutron star models, are comparable to what would be observed by a next-generation X-ray timing instrument like ASTROSAT, NICER, or LOFT. In this paper, we showcase the use of an evolutionary optimization algorithm to fit pulse profiles to determine the best-fit masses and radii. By fitting synthetic data, we assess how well the optimization algorithm can recover the input parameters. Multiple Poisson realizations of the synthetic pulse profiles, constructed with 1.6 million counts and no background, were fitted with the Ferret algorithm to analyze both statistical and degeneracy-related uncertainty, and to explore how the goodness-of-fit depends on the input parameters. For the regions of parameter space sampled by our tests, the best-determined parameter is the projected velocity of the spot along the observer's line-of-sight, with an accuracy of $\le3$% compared to the true value and with $\le5$\% statistical uncertainty. The next best-determined are the mass and radius; for a neutron star with a spin frequency of 600\,Hz, the best-fit mass and radius are accurate to $\le5$%, with respective uncertainties of $\le7$% and $\le10$%. The accuracy and precision depend on the observer inclination and spot co-latitude, with values of $\sim1$ % achievable in mass and radius if both the inclination and co-latitude are > 60 degrees.

Fitting the thermal X-ray spectra of neutron stars (NSs) in quiescent X-ray binaries can constrain the masses and radii of NSs. The effect of undetected hot spots on the spectrum, and thus on the inferred NS mass and radius, has not yet been explored for appropriate atmospheres and spectra. A hot spot would harden the observed spectrum, so that spectral modeling tends to infer radii that are too small. However, a hot spot may also produce detectable pulsations. We simulated the effects of a hot spot on the pulsed fraction and spectrum of the quiescent NSs X5 and X7 in the globular cluster 47 Tucanae, using appropriate spectra and beaming for hydrogen atmosphere models, incorporating special and general relativistic effects, and sampling a range of system angles. We searched for pulsations in archival Chandra HRC-S observations of X5 and X7, placing 90\% confidence upper limits on their pulsed fractions below 16\%. We use these pulsation limits to constrain the temperature differential of any hot spots, and to then constrain the effects of possible hot spots on the X-ray spectrum and the inferred radius from spectral fitting. We find that hot spots below our pulsation limit could bias the spectroscopically inferred radius downward by up to 28\%. For Cen X-4 (which has deeper published pulsation searches), an undetected hot spot could bias its inferred radius downward by up to 10\%. Improving constraints on pulsations from quiescent LMXBs may be essential for progress in constraining their radii.

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

The Neutron-star Interior Composition Explorer (NICER) is an X-ray astrophysics payload that will be placed on the International Space Station. Its primary science goal is to measure with high accuracy the pulse profiles that arise from the non-uniform thermal surface emission of rotation-powered pulsars. Modeling general relativistic effects on the profiles will lead to measuring the radii of these neutron stars and to constraining their equation of state. Achieving this goal will depend, among other things, on accurate knowledge of the source, sky, and instrument backgrounds. We use here simple analytic estimates to quantify the level at which these backgrounds need to be known in order for the upcoming measurements to provide significant constraints on the properties of neutron stars. We show that, even in the minimal-information scenario, knowledge of the background at a few percent level for a background-to-source countrate ratio of 0.2 allows for a measurement of the neutron star compactness to better than 10% uncertainty for most of the parameter space. These constraints improve further when more realistic assumptions are made about the neutron star emission and spin, and when additional information about the source itself, such as its mass or distance, are incorporated.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

We calculate the rotational broadening in the observed thermal spectra of neutron stars spinning at moderate rates in the Hartle-Thorne approximation. These calculations accurately account for the effects of the second-order Doppler boosts as well as for the oblate shapes and the quadrupole moments of the neutron stars. We find that fitting the spectra and inferring the bolometric fluxes under the assumption that a star is not rotating causes an underestimate of the inferred fluxes and, thus, radii. The correction depends on the stellar spin, radius, and observer's inclination. For a 10 km neutron star spinning at 600 Hz, the rotational correction to the flux is ~1-4%, while for a 15 km neutron star with the same spin period, the correction ranges from 2% for pole-on sources to 12% for edge-on sources. We calculate the inclination-averaged corrections to inferred radii as a function of the neutron-star radius and mass and provide an empirical formula for the corrections. For realistic neutron star parameters (1.4 M$_\odot$, 12 km, 600 Hz), the stellar radius is on the order of 4% larger than the radius inferred under the assumption that the star is not spinning.

On the surface of a rapidly rotating neutron star, the effective centrifugal force decreases the effective acceleration due to gravity (as measured in the rotating frame) at the equator while increasing the acceleration at the poles due to the centrifugal flattening of the star into an oblate spheroid. We compute the effective gravitational acceleration for relativistic rapidly rotating neutron stars and show that for a star with mass $M$, equatorial radius $R_e$, and angular velocity $\Omega$, the deviations of the effective acceleration due to gravity from the nonrotating case take on a universal form that depends only on the compactness ratio $M/R_e$, the dimensionless square of the angular velocity $\Omega^2R_e^3/GM$, and the latitude on the star's surface. This dependence is universal, in that it has very little dependence on the neutron star's equation of state. The effective gravity is expanded in the slow rotation limit to show the dependence on the effective centrifugal force, oblate shape of the star and the quadrupole moment of the gravitational field. In addition, an empirical fit and simple formula for the effective gravity is found. We find that the increase in the acceleration due to gravity at the poles is of the same order of magnitude as the decrease in the effective acceleration due to gravity at the equator for all realistic value of mass, radius and spin. For neutron stars that spin with frequencies near 600 Hz the difference between the effective gravity at the poles and the equator is about 20%.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

The accreting millisecond pulsar XTE J1807-294 is studied through a pulse shape modeling analysis. The model includes blackbody and Comptonized emission from the one visible hot spot and makes use of the Oblate Schwarzschild approximation for ray-tracing. We include a scattered light contribution, which accounts for flux scattered off an equatorial accretion disk to the observer including time delays in the scattered light. We give limits to mass and radius for XTE J1807-294 and compare to limits determined for SAX J1808-3658 and XTE J1814-334 previously determined using similar methods. The resulting allowed region for mass-radius curves is small but is consistent with a mass-radius relation with nearly constant radius (~12 km) for masses between 1 and 2.5 solar masses.

The pulse shapes detected during multiple outbursts of SAX J1808 are analyzed in order to constrain the neutron star's mass and radius. We use a hot-spot model with a small scattered-light component to jointly fit data from two different epochs, under the restriction that the star's mass and radius and the binary's inclination do not change from epoch to epoch. All other parameters describing the spot location, emissivity, and relative fractions of blackbody to Comptonized radiation are allowed to vary with time. The joint fit of data from the 1998 "slow decay" and the 2002 "end of outburst maximum" epochs using the constraint i<90 degrees leads to the 3 sigma confidence constraint on the neutron star mass 0.8 M_sun < M < 1.7 M_sun and equatorial radius 5 km < R < 13 km. Inclinations as low as 41 degrees are allowed. The best-fit models with M > 1.0 M_sun from joint fits of the 1998 data with data from other epochs of the 2002 and 2005 outbursts also fall within the same 3 sigma confidence region. This 3 sigma confidence region allows a wide variety of hadronic equations of state, in contrast with an earlier analysis (Leahy et al 2008) of only the 1998 outburst data that only allowed for extremely small stars.

The accretion-powered (non-X-ray burst) pulsations of XTE J1814-338 are modeled to determine neutron star parameters and their uncertainties. The model is a rotating circular hot spot and includes: (1) an isotropic blackbody spectral component; (2) an anisotropic Comptonized spectral component; (3) relativistic time-delays and light-bending; and (4) the oblate shape of the star due to rotation. This model is the simplest possible model that is consistent with the data. The resulting best-fit parameters of the model favor stiff equations of state, as can be seen from the 3-sigma allowed regions in the mass-radius diagram. We analyzed all data combined from a 23 day period of the 2003 outburst, and separately analyzed data from 2 days of the outburst. The allowed mass-radius regions for both cases only allow equations of state (EOS) that are stiffer than EOS APR (Akmal et al. 1998), consistent with the large mass that has been inferred for the pulsar NGC 6440B (Freire et al. 2008). The stiff EOS inferred by this analysis is not compatible with the soft EOS inferred from a similar analysis of SAX J1808.

SAX J1808.4-3658 has a 2.5 millisecond neutron star rotation period and exhibits X-ray pulsations due to its rotating hot spot. Here we present an analysis of the pulse shapes of SAX J1808 during its 1998 outburst. The modeling of the pulse shape includes several effects, including gravitational light-bending, doppler effects and two spectral components with different emissivity. In addition we include the new effects of light-travel time-delays and the neutron star's oblate shape. We also consider two different data sets, with different selection in time period (1 day versus 19 days of data combined) and different energy binning and time resolution. We find that including time-delays and oblateness results in stronger restriction on allowed masses and radii. A second result is that the choice of data selection strongly affects the allowed masses and radii. Overall, the derived constraints on mass and radius favor compact stars and a soft equation of state.

We present a simple method for including the oblateness of a rapidly rotating neutron star when fitting X-ray light curves. In previous work we showed that the oblateness induced by rotation at frequencies above 300 Hz produces a geometric effect which needs to be accounted for when modelling light curves to extract constraints on the neutron star's mass and radius. In our model X-rays are emitted from the surface of an oblate neutron star and propagate to the observer along geodesics of the Schwarzschild metric for a spherical neutron star. Doppler effects due to rotation are added in the same manner as in the case of a spherical neutron star. We show that this model captures the most important effects due to the neutron star's rotation. We also explain how the geometric oblateness effect can rival the Doppler effect for some emission geometries.

Raytracing computations for light emitted from the surface of a rapidly rotating neutron star are carried out in order to construct light curves for accreting millisecond pulsars. These calculations are for realistic models of rapidly rotating neutron stars which take into account both the correct exterior metric and the oblate shape of the star. We find that the most important effect, comparing the full raytracing computations with simpler approximations currently in use, arises from the oblate shape of the rotating star. Approximating a rotating neutron star as a sphere introduces serious errors in fitted values of the star's radius and mass if the rotation rate is very large. However, for lower rotation rates acceptable mass and radius values can be obtained using the spherical approximation.

We present raytracing computations for light emitted from the surface of a rapidly-rotating neutron star in order to construct light curves for X-ray pulsars and bursters. These calculations are for realistic models of rapidly-rotating neutron stars which take into account both the correct exterior metric and the oblate shape of the star. We find that the most important effect arising from rotation comes from the oblate shape of the rotating star. We find that approximating a rotating neutron star as a sphere introduces serious errors in fitted values of the star's radius and mass if the rotation rate is very large. However, in most cases acceptable fits to the ratio M/R can be obtained with the spherical approximation.

The discovery of the first gravitationally redshifted spectral line from a neutron star (NS) by Cottam, Paerels and Mendez has triggered theoretical studies of the physics of atomic line formation in NS atmospheres. Chang, Bildsten and Wasserman showed that the hydrogenic Fe H$\alpha$ line formed above the photosphere of a bursting NS is intrinsically broad. We now include rotational broadening within general relativity and compare the resulting profile to that observed during Type I bursts from EXO 0748-676. We show that the fine structure splitting of the line precludes a meaningful constraint on the radius. Our fitting of the data show that the line forming Fe column is ${\rm log}_{10} (N_{\rm Fe, n=2}/{\rm cm^{-2}})=17.9_{-0.42}^{+0.27}$ and gravitational redshift $1+z =1.345_{-0.008}^{+0.005}$ with 95% confidence. We calculate the detectability of this spectral feature for a large range of spins and inclinations assuming that the emission comes from the entire surface. We find that at 300 (600) Hz only 10-20% (5-10%) of NSs would have spectral features as deep as that seen in EXO 0748-676.

We demonstrate that fitted values of stellar radius obtained by fitting theoretical light curves to observations of millisecond period X-ray pulsars can significantly depend on the method used to calculate the light curves. The worst-case errors in the fitted radius are evaluated by restricting ourselves to the case of light emitted and received in the equatorial plane of a rapidly-rotating neutron star. First, using an approximate flux which is adapted to the one-dimensional nature of such an emission region, we show how pulse shapes can be constructed using an exact spacetime metric and fully accounting for time-delay effects. We compare this to a method which approximates the exterior spacetime of the star by the Schwarzschild metric, inserts special relativistic effects by hand, and neglects time-delay effects. By comparing these methods, we show that there are significant differences in these methods for some applications, for example pulse timing and constraining the stellar radius. In the case of constraining the stellar radius, we show that fitting the approximate pulse shapes to the full calculation yields errors in the fitted radius of as much as about +/- 10%, depending on the rotation rate and size of the star as well as the details describing the emitting region. However, not all applications of pulse shape calculations suffer from significant errors: we also show that the calculation of the soft-hard phase lag for a 1 keV blackbody does not strongly depend on the method used for calculating the pulse shapes.

In a relativistic model of a neutron star, the star's mass is less than the mass of the individual component baryons. This is due to the fact that the star's negative binding energy makes a contribution to the star's total energy and its mass. A consequence of this relativistic mass deficit is that a neutron star that is accreting matter increases its mass at a rate which is slower than the mass of a baryon times the rate that baryons are accreted. This difference in the rate of change of the masses has a simple relation with the star's gravitational redshift. We show that this effect has the potential to be observed in binaries where the mass transfer is driven by angular momentum losses from the gravitational radiation emitted by the binary motion.

Rossby waves (r-modes) in rapidly rotating neutron stars are unstable because of the emission of gravitational radiation. We study saturation of this instability by nonlinear transfer of energy to stellar "inertial" oscillation modes. We present detailed calculations of stellar inertial modes in the WKB limit, their linear damping by bulk and shear viscosity, and the nonlinear coupling forces among these modes. The saturation amplitude is derived in the extreme limits of strong or weak driving by radiation reaction, as compared to the damping rate of low order inertial modes. We find the saturation energy is {\it extremely small}, at least four orders of magnitude smaller than that found by previous investigators. We discuss the consequences of this result for spin evolution of young neutron stars, and neutron stars being spun up by accretion in Low Mass X-ray Binaries.We also discuss the detection of these gravitational waves by LIGO.

We investigate the effect of a magnetic field on the global oscillation modes of a rotating fluid star in the magnetohydrodynamic approximation. We present general equations for the modification of any type of fluid mode due to a general magnetic field which is not aligned with the star's spin axis. In the case of any internal dipole magnetic field we derive the equations for the frequency corrections to the r-modes. We solve for the frequency correction explicitly for the case when the internal dipole field is force-free, including the uniform density case. In the weak-field limit, the spatial form of the r-mode velocity perturbation is unchanged, but the magnitude of the frequency in the rotating frame increases.

The r-modes of neutron stars can be driven unstable by gravitational radiation. While linear perturbation theory predicts the existence of this instability, linear theory can't provide any information about the nonlinear development of the instability. The subject of this paper is the weakly nonlinear regime of fluid dynamics. In the weakly nonlinear regime, the nonlinear fluid equations are approximated by an infinite set of oscillators which are coupled together so that terms quadratic in the mode amplitudes are kept in the equations of motion. In this paper, the coupling coefficients between the r-modes are computed. The stellar model assumed is a polytropic model where a source of buoyancy is included so that the Schwarzschild discriminant is nonzero. The properties of these coupling coefficients and the types of resonances possible are discussed in this paper. It is shown that no exact resonance involving the unstable $l=m=2$ r-mode occur and that only a small number of modes have a dimensionless coupling constant larger than unity. However, an infinite number of resonant mode triplets exist which couple indirectly to the unstable r-mode. All couplings in this paper involve the l>|m| r-modes which only exist if the star is slowly rotating. This work is complementary to that of Schenk et al (2002) who consider rapidly rotating stars which are neutral to convection.

We present calculations of the spin-down of a neutron star atmosphere due to hydrostatic expansion during a Type I X-ray burst. We show that (i) Cumming and Bildsten overestimated the spin-down of rigidly-rotating atmospheres by a factor of two, and (ii) general relativity has a small (5-10%) effect on the angular momentum conservation law. We rescale our results to different neutron star masses, rotation rates and equations of state, and present some detailed rotational profiles. Comparing with recent observations of large frequency shifts in MXB 1658-298 and 4U 1916-053, we find that the spin-down expected if the atmosphere rotates rigidly is a factor of two to three less than the observed values. If differential rotation is allowed to persist, we find that the upper layers of the atmosphere spin down by an amount comparable to the observed values; however, there is no compelling reason to expect the observed spin frequency to be that of only the outermost layers. We conclude that hydrostatic expansion and angular momentum conservation alone cannot account for the largest frequency shifts observed during Type I bursts.

A remarkable correlation between the centroid frequencies of quasi periodic oscillations, QPOs, (or peaked noise components) from low mass X-ray binaries, has been recently discovered by Psaltis, Belloni and van der Klis (1999). This correlation extends over nearly 3 decades in frequency and encompasses both neutron star and black hole candidate systems. We discuss this result in the light of the relativistic precession model, which has been proposed to interpret the kHz QPOs as well as some of the lower frequency QPOs of neutron star low mass X-ray binaries of the Atoll and Z classes. Unlike other models the relativistic precession model does not require the compact object to be a neutron star and can be applied to black hole candidates as well. We show that the predictions of the relativistic precession model match both the value and dependence of the correlation to a very good accuracy without resorting to additional assumptions.

General relativity predicts that a rotating body produces a frame-dragging (or Lense-Thirring) effect: the orbital plane of a test particle in a non-equatorial orbit precesses about the body's symmetry axis. In this paper we compute the precession frequencies of circular orbits around rapidly rotating neutron stars for a variety of masses and equations of state. The precession frequencies computed are expressed as numerical functions of the orbital frequency observed at infinity. The post-Newtonian expansion of the exact precession formula is examined to identify the relative magnitudes of the precession caused by the Lense-Thirring effect, the usual Newtonian quadrupole effect and relativistic corrections. The first post-Newtonian correction to the Newtonian quadrupole precession is derived in the limit of slow rotation. We show that the post-Newtonian precession formula is a good approximation to the exact precession close to the neutron star in the slow rotation limit (up to \sim 400 Hz in the present context). The results are applied to recent RXTE observations of neutron star low-mass X-ray binaries, which display kHz quasi-periodic oscillations and, within the framework of beat frequency models, allow the measurement of both the neutron star spin frequency and the Keplerian frequency of the innermost ring of matter in the accretion disk around it. For a wide range of realistic equations of state, we find that the predicted precession frequency of this ring is close to one half of the low-frequency (\sim 20 - 35 Hz) quasi-periodic oscillations seen in several Atoll sources.

Perturbations of rotating relativistic stars can be classified by their behavior under parity. For axial perturbations (r-modes), initial data with negative canonical energy is found with angular dependence $e^{im\phi}$ for all values of $m\geq 2$ and for arbitrarily slow rotation. This implies instability (or marginal stability) of such perturbations for rotating perfect fluids. This low $m$-instability is strikingly different from the instability to polar perturbations, which sets in first for large values of $m$. The timescale for the axial instability appears, for small angular velocity $\Omega$, to be proportional to a high power of $\Omega$. As in the case of polar modes, viscosity will again presumably enforce stability except for hot, rapidly rotating neutron stars. This work complements Andersson's numerical investigation of axial modes in slowly rotating stars.

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.