Numerical approach for high precision 3-D relativistic star models

A multi-domain spectral method for computing very high precision 3-D stellar models is presented. The boundary of each domain is chosen in order to coincide with a physical discontinuity (e.g. the star's surface). In addition, a regularization procedure is introduced to deal with the infinite derivatives on the boundary that may appear in the density field when stiff equations of state are used. Consequently all the physical fields are smooth functions on each domain and the spectral method is absolutely free of any Gibbs phenomenon, which yields to a very high precision. The power of this method is demonstrated by direct comparison with analytical solutions such as MacLaurin spheroids and Roche ellipsoids. The relative numerical error reveals to be of the order of $10^{-10}$. This approach has been developed for the study of relativistic inspiralling binaries. It may be applied to a wider class of astrophysical problems such as the study of relativistic rotating stars too.Comment: Minor changes, Phys. Rev. D in pres

Inertial modes in slowly rotating stars : an evolutionary description

We present a new hydro code based on spectral methods using spherical coordinates. The first version of this code aims at studying time evolution of inertial modes in slowly rotating neutron stars. In this article, we introduce the anelastic approximation, developed in atmospheric physics, using the mass conservation equation to discard acoustic waves. We describe our algorithms and some tests of the linear version of the code, and also some preliminary linear results. We show, in the Newtonian framework with differentially rotating background, as in the relativistic case with the strong Cowling approximation, that the main part of the velocity quickly concentrates near the equator of the star. Thus, our time evolution approach gives results analogous to those obtained by Karino {\it et al.} \cite{karino01} within a calculation of eigenvectors. Furthermore, in agreement with the work of Lockitch {\it et al.} \cite{lockandf01}, we found that the velocity seems to always get a non-vanishing polar part.Comment: 36 pages, 27 figures, accepted for publication in Phys. Rev. D (discussion added in the introduction

Directionality effects in the transfer of X-rays from a magnetized atmosphere: Beam pulse shape

A formalism is presented for radiation transfer in two normal polarization modes in finite and semiinfinite plane parallel uniform atmospheres with a magnetic field perpendicular to the surface and arbitrary propagation angles. This method is based on the coupled integral equations of transfer, including emission, absorption, and scattering. Calculations are performed for atmosphere parameters typical of X-ray pulsars. The directionality of the escaping radiation is investigated for several cases, varying the input distributions. Theoretical pencil beam profiles and X-ray pulse shapes are obtained assuming the radiation is emitted from the polar caps of spinning neutron stars. Implications for realistic models of accreting magnetized X-ray sources are briefly discussed