Aims. Many recent observations of pulsars and magnetars can be interpreted in terms of neutron stars (NS) with multipole electromagnetic fields. As a first approximation, we investigate the multipole magnetic and electric fields in the environment of a rotating star when this environment is deprived of plasma. Methods. We compute a multipole expansion of the electromagnetic field in vacuum for a given magnetic field on the conducting surface of the rotating star. Then, we consider a few consequences of multipole fields of pulsars. Results. We provide an explicit form of the solution. For each spherical harmonic of the magnetic field, the expansion contains a finite number of terms. A multipole magnetic field can provide an explanation for the stable sub-structures of pulses, and they offer a solution to the problem of current closure in pulsar magnetospheres. Conclusions. This computation generalises the widely used model of a rotating star in vacuum with a dipole field. It can be especially useful as a first approximation to the electromagnetic environment of a compact star, for instance a neutron star, with an arbitrarily magnetic field.

Electron-positron pair production by collision of photons is investigated in view of application to pulsar physics. We compute the absorption rate of individual gamma-ray photons by an arbitrary anisotropic distribution of softer photons, and the energy and angular spectrum of the outgoing leptons. We work analytically within the approximation that 1 mc 2 /E > /E, with E and the gamma-ray and soft-photon maximum energy and mc 2 the electron mass energy. We give results at leading order in these small parameters. For practical purposes, we provide expressions in the form of Laurent series which give correct reaction rates in the isotropic case within an average error of $\sim$ 7%. We apply this formalism to gamma rays flying downward or upward from a hot neutron star thermally radiating at a uniform temperature of 10 6 K. Other temperatures can be easily deduced using the relevant scaling laws. We find differences in absorption between these two extreme directions of almost two orders of magnitude, much larger than our error estimate. The magnetosphere appears completely opaque to downward gamma rays while there are up to $\sim$ 10% chances of absorbing an upward gamma ray. We provide energy and angular spectra for both upward and downward gamma rays. Energy spectra show a typical double peak, with larger separation at larger gamma-ray energies. Angular spectra are very narrow, with an opening angle ranging from 10 --3 to 10 --7 radians with increasing gamma-ray energies.

In a previous paper, we derived the quantum states of a Dirac particle in a circular, intense magnetic field in the limit of low momentum perpendicular to the field with the purpose of giving a quantum description of the trajectory of an electron, or a positron, in a typical pulsar or magnetar magnetosphere. Here we continue this work by computing the radiation resulting from transitions between these states. This leads to derive from first principles a quantum theory of the so-called curvature and synchro-curvature radiations relevant for rotating neutron-star magnetospheres. We find that, within the approximation of an infinitely confined wave-function around the magnetic field and in the continuous energy-level limit, classical curvature radiation can be recovered in a fully consistent way. Further we introduce discrete transitions to account for the change of momentum perpendicular to the field and derive expressions for what we call quantum synchro-curvature radiation. Additionally, we express deconfinement and quantum recoil corrections.

Neutron-star magnetospheres are structured by very intense magnetic fields extending from 100 to 10 5 km traveled by very energetic electrons and positrons with Lorentz factors up to $\sim$ 10 7. In this context, particles are forced to travel almost along the magnetic field with very small gyro-motion, potentially reaching the quantified regime. We describe the state of Dirac particles in a locally uniform, constant and curved magnetic field in the approximation that the Larmor radius is very small compared to the radius of curvature of the magnetic field lines. We obtain a result that admits the usual relativistic Landau states as a limit of null curvature. We will describe the radiation of these states, that we call quantum curvature or synchro-curvature radiation, in an upcoming paper.

We present a new code, Lorene's Ghost (for Lorene's gravitational handling of spectral transport) developed to treat the problem of neutrino transport in supernovae with the use of spectral methods. First, we derive the expression for the nonrelativistic Liouville operator in doubly spherical coordinates (r, theta, phi, epsilon, Theta, Phi)$, and further its general relativistic counterpart. We use the 3 + 1 formalism with the conformally flat approximation for the spatial metric, to express the Liouville operator in the Eulerian frame. Our formulation does not use any approximations when dealing with the angular arguments (theta, phi, Theta, Phi), and is fully energy-dependent. This approach is implemented in a spherical shell, using either Chebyshev polynomials or Fourier series as decomposition bases. It is here restricted to simplified collision terms (isoenergetic scattering) and to the case of a static fluid. We finish this paper by presenting test results using basic configurations, including general relativistic ones in the Schwarzschild metric, in order to demonstrate the convergence properties, the conservation of particle number and correct treatment of some general-relativistic effects of our code. The use of spectral methods enables to run our test cases in a six-dimensional setting on a single processor.

We investigate whether there may be one or many companions orbiting at close distance to the light cylinder around the extremely intermittent pulsars PSR B1931+24 and PSR J1841-0500. These pulsars, behaving in a standard way when they are active, also "switch off" for durations of several days, during which their magnetospheric activity is interrupted or reduced. We constrained our analysis on eight fundamental properties of PSR B1931+24 that summarise the observations. We considered that the disruption/activation of the magnetospheric activity would be caused by the direct interaction of the star with the Alfv\'en wings emanating from the companions. We also considered the recurrence period of 70 days to be the period of precession of the periastron of the companions orbit. We analysed in which way the time scale of the "on/off" pseudo-cycle would be conditioned by the precession of the periastron and not by the orbital time scale, and we derived a set of orbital constraints that we solved. We then compared the model, based on PSR 1931+24, with the known properties of PSR 1841+0500. We conclude that PSR B1931+24 may be surrounded at a close distance to the star by a stream of small bodies of kilometric or sub-kilometric sizes that could originate from the tidal disruption of a body of moderate size that fell at a close distance to the neutron star on an initially very eccentric orbit. This scenario is also compatible with the properties of PSR J1841-0500, although the properties of PSR J1841-0500 are, by now, less constrained. These results raise new questions. Why are the asteroids not yet evaporated ? What kind of interaction can explain the disruption of the magnetospheric activity ? These questions are the object of two papers in preparation that will complete the present analysis.

We investigate whether one or many companions are orbiting the extremely intermittent pulsar PSR B1931+24. We constrained our analysis on previous observations of eight fundamental properties of PSR B1931+24. The most puzzling properties are the intermittent nature of the pulsar's activity, with active and quiet phases that alternate quasi-periodically; the variation of the slowing-down rate of its period between active and quiet phases; and because there are no timing residuals, it is highly unlikely that the pulsar has a massive companion. Here, we examine the effects that one putative companion immersed in the magnetospheric plasma or the wind of the pulsar might have, as well as the associated electric current distribution. We analysed several possibilities for the distance and orbit of this hypothetical companion and the nature of its interaction with the neutron star. We show that if the quasi-periodic behaviour of PSR B1931+24 was caused by a companion orbiting the star with a period of 35 or 70 days, the radio emissions, usually considered to be those of the pulsar would in that specific case be emitted in the companion's environment. We analysed four possible configurations and conclude that none of them would explain the whole set of peculiar properties of PSR 1931+24. We furthermore considered a period 70 days for the precession of the periastron associated to an orbit very close to the neutron star. This hypothesis is analysed in a companion paper.

We present a numerical method for handling the resolution of a general transport equation for radiative particles, aimed at physical problems with a general spherical geometry. Having in mind the computational time difficulties encountered in problems such as neutrino transport in astrophysical supernovae, we present a scheme based on full spectral methods in 6d spherical coordinates. This approach, known to be suited when the characteristic length of the dynamics is much smaller than the domain size, has the potential advantage of a global speedup with respect to usual finite difference schemes. An analysis of the properties of the Liouville operator expressed in our coordinates is necessary in order to handle correctly the numerical behaviour of the solution. This reflects on a specific (spherical) geometry of the computational domain. The numerical tests, performed under several different regimes for the equation, prove the robustness of the scheme: their performances also point out to the suitability of such an approach to large scale computations involving transport physics for mass less radiative particles.

We discuss whether the winding-up of the magnetic field by differential rotation in a new-born quark star can produce a sufficiently-high, energy, emission rate of sufficiently long duration to explain long gamma-ray bursts. In the context of magnetohydrodynamics, we study the torsional oscillations and energy extraction from a new-born, hot, differentially rotating quark star. The new-born compact star is a rapid rotator that produces a relativistic, leptonic wind. The star's torsional oscillation modulates this wind emission considerably when it is odd and of sufficient amplitude, which is relatively easy to reach. Odd oscillations may occur just after the formation of a quark star. Other asymmetries can cause similar effects. The buoyancy of wound-up magnetic fields is inhibited, or its effects are limited, by a variety of different mechanisms. Direct electromagnetic emission by the torsional oscillation in either an outside vacuum or the leptonic wind surrounding the compact object is found to be insignificant. In contrast, the twist given to the outer magnetic field by an odd torsional oscillation is generally sufficient to open the star's magnetosphere. The Poynting emission of the star in its leptonic environment is then radiated from all of its surface and is enhanced considerably during these open episodes, tapping at the bulk rotational energy of the star. This results in intense energy shedding in the first tens of minutes after the collapse of magnetized quark stars with an initial poloidal field of order of 10**14 Gauss, sufficient to explain long gamma-ray bursts.

A flexible spectral code for the study of general relativistic magnetohydrodynamics is presented. Aiming at investigating the physics of slowly rotating magnetized compact stars, this new code makes use of various physically motivated approximations. Among them, the relativistic anelastic approximation is a key ingredient of the current version of the code. In this article, we mainly outline the method, putting emphasis on algorithmic techniques that enable to benefit as much as possible of the non-dissipative character of spectral methods, showing also a potential astrophysical application and providing a few illustrative tests.

With (non-barotropic) equations of state valid even when the neutron, proton and electron content of neutron star cores is not in beta equilibrium, we study inertial and composition gravity modes of relativistic rotating neutron stars. We solve the relativistic Euler equations in the time domain with a three dimensional numerical code based on spectral methods, in the slow rotation, relativistic Cowling and anelastic approximations. Principally, after a short description of the gravity modes due to smooth composition gradients, we focus our analysis on the question of how the inertial modes are affected by non-barotropicity of the nuclear matter. In our study, the deviation with respect to barotropicity results from the frozen composition of non-superfluid matter composed of neutrons, protons and electrons, when beta equilibrium is broken by millisecond oscillations. We show that already for moderatly fast rotating stars the increasing coupling between polar and axial modes makes those two cases less different than for very slowly rotating stars. In addition, as we directly solve the Euler equations, without coupling only a few number of spherical harmonics, we always found, for the models that we use, a discrete spectrum for the $l = m = 2$ inertial mode. Finally, we find that, for non-barotropic stars, the frequency of this mode, which is our main focus, decreases in a non-negligible way, whereas the time dependence of the energy transfer between polar and axial modes is substantially different due to the existence of low-frequencies gravity modes.

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.

Binary black hole systems in the pre-coalescence stage are numerically constructed by demanding that the associated spacetime admits a helical Killing vector. Comparison with third order post-Newtonian calculations indicates a rather good agreement until the innermost stable circular orbit.

We present a new approach to the problem of binary black holes in the pre-coalescence stage, i.e. when the notion of orbit has still some meaning. Contrary to previous numerical treatments which are based on the initial value formulation of general relativity on a (3-dimensional) spacelike hypersurface, our approach deals with the full (4-dimensional) spacetime. This permits a rigorous definition of the orbital angular velocity. Neglecting the gravitational radiation reaction, we assume that the black holes move on closed circular orbits, which amounts to endowing the spacetime with a helical Killing vector. We discuss the choice of the spacetime manifold, the desired properties of the spacetime metric, as well as the choice of the rotation state for the black holes. As a simplifying assumption, the space 3-metric is approximated by a conformally flat one. The problem is then reduced to solving five of the ten Einstein equations, which are derived here, as well as the boundary conditions on the black hole surfaces and at spatial infinity. We exhibit the remaining five Einstein equations and propose to use them to evaluate the error induced by the conformal flatness approximation. The orbital angular velocity of the system is computed through a requirement which reduces to the classical virial theorem at the Newtonian limit.

We present the first results from a new method for computing spacetimes representing corotating binary black holes in circular orbits. The method is based on the assumption of exact equilibrium. It uses the standard 3+1 decomposition of Einstein equations and conformal flatness approximation for the 3-metric. Contrary to previous numerical approaches to this problem, we do not solve only the constraint equations but rather a set of five equations for the lapse function, the conformal factor and the shift vector. The orbital velocity is unambiguously determined by imposing that, at infinity, the metric behaves like the Schwarzschild one, a requirement which is equivalent to the virial theorem. The numerical scheme has been implemented using multi-domain spectral methods and passed numerous tests. A sequence of corotating black holes of equal mass is calculated. Defining the sequence by requiring that the ADM mass decrease is equal to the angular momentum decrease multiplied by the orbital angular velocity, it is found that the area of the apparent horizons is constant along the sequence. We also find a turning point in the ADM mass and angular momentum curves, which may be interpreted as an innermost stable circular orbit (ISCO). The values of the global quantities at the ISCO, especially the orbital velocity, are in much better agreement with those from third post-Newtonian calculations than with those resulting from previous numerical approaches.

We study equilibrium sequences of close binary systems composed of identical polytropic stars in Newtonian gravity. The solving method is a multi-domain spectral method which we have recently developed. An improvement is introduced here for accurate computations of binary systems with stiff equation of state ($\gamma > 2$). The computations are performed for both cases of synchronized and irrotational binary systems with adiabatic indices $\gamma=3,~2.5,~2.25,~2$ and 1.8. It is found that the turning points of total energy along a constant-mass sequence appear only for $\gamma \ge 1.8$ for synchronized binary systems and $\gamma \ge 2.3$ for irrotational ones. In the synchronized case, the equilibrium sequences terminate by the contact between the two stars. On the other hand, for irrotational binaries, it is found that the sequences terminate at a mass shedding limit which corresponds to a detached configuration.

We present a spectral method for solving elliptic equations which arise in general relativity, namely three-dimensional scalar Poisson equations, as well as generalized vectorial Poisson equations of the type $\Delta \vec{N} + \lambda \vec{\nabla}(\vec{\nabla}\cdot \vec{N}) = \vec{S}$ with $\lambda \not= -1$. The source can extend in all the Euclidean space ${\bf R}^3$, provided it decays at least as $r^{-3}$. A multi-domain approach is used, along with spherical coordinates $(r,\theta,\phi)$. In each domain, Chebyshev polynomials (in $r$ or $1/r$) and spherical harmonics (in $\theta$ and $\phi$) expansions are used. If the source decays as $r^{-k}$ the error of the numerical solution is shown to decrease at least as $N^{-2(k-2)}$, where $N$ is the number of Chebyshev coefficients. The error is even evanescent, i.e. decreases as $\exp(-N)$, if the source does not contain any spherical harmonics of index $l\geq k -3$ (scalar case) or $l\geq k-5$ (vectorial case).

We present a numerical method to compute quasiequilibrium configurations of close binary neutron stars in the pre-coalescing stage. A hydrodynamical treatment is performed under the assumption that the flow is either rigidly rotating or irrotational. The latter state is technically more complicated to treat than the former one (synchronized binary), but is expected to represent fairly well the late evolutionary stages of a binary neutron star system. As regards the gravitational field, an approximation of general relativity is used, which amounts to solving five of the ten Einstein equations (conformally flat spatial metric). The obtained system of partial differential equations is solved by means of a multi-domain spectral method. Two spherical coordinate systems are introduced, one centered on each star; this results in a precise description of the stellar interiors. Thanks to the multi-domain approach, this high precision is extended to the strong field regions. The computational domain covers the whole space so that exact boundary conditions are set to infinity. Extensive tests of the numerical code are performed, including comparisons with recent analytical solutions. Finally a constant baryon number sequence (evolutionary sequence) is presented in details for a polytropic equation of state with gamma=2.

Exact models of uniformly rotating strange stars, built of self bound quark matter, are calculated within the framework of general relativity. This is made possible thanks to a new numerical technique capable to handle the strong density discontinuity at the surface of these stars. Numerical calculations are done for a simple MIT bag model equation of state of strange quark matter. Evolutionary sequences of models of rotating strange stars at constant baryon mass are calculated. Maximally rotating configurations of strange stars are determined, assuming that the rotation frequency is limited by the mass shedding and the secular instability with respect to axisymmetric perturbations. Exact formulae which give the dependence of the maximum rotation frequency, and of the maximum mass and corresponding radius of rotating configurations, on the value of the bag constant, are obtained. The values of T/W for rapidly rotating massive strange stars are significantly higher than those for ordinary neutron stars. This might indicate particular susceptibility of rapidly rotating strange stars to triaxial instabilities.

We present results of numerical computations of quasiequilibrium sequences of binary neutron stars with zero vorticity, in the general relativistic framework. The Einstein equations are solved under the assumption of a conformally flat spatial 3-metric (Wilson-Mathews approximation). The evolution of the central density of each star is monitored as the orbit shrinks in response to gravitational wave emission. For a compactification ratio M/R=0.14, the central density remains rather constant (with a slight increase, below 0.1%) before decreasing. For a higher compactification ratio M/R=0.17 (i.e. stars closer to the maximum mass configuration), a very small density increase (at most 0.3%) is observed before the decrease. This effect remains within the error induced by the conformally flat approximation. It can be thus concluded that no substantial compression of the stars is found, which would have indicated a tendency to individually collapse to black hole prior to merger. Moreover, no turning point has been found in the binding energy or angular momentum along evolutionary sequences, which may indicate that these systems do not have any innermost stable circular orbit (ISCO).

We report on general relativistic calculations of quasiequilibrium configurations of binary neutron stars in circular orbits with zero vorticity. These configurations are expected to represent realistic situations as opposed to corotating configurations. The Einstein equations are solved under the assumption of a conformally flat spatial 3-metric (Wilson-Mathews approximation). The velocity field inside the stars is computed by solving an elliptical equation for the velocity scalar potential. Results are presented for sequences of constant baryon number (evolutionary sequences). Although the central density decreases much less with the binary separation than in the corotating case, it still decreases. Thus, no tendency is found for the stars to individually collapse to black hole prior to merger.

We present spectral methods developed in our group to solve three-dimensional partial differential equations. The emphasis is put on equations arising from astrophysical problems in the framework of general relativity.

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.

Both observational and theoretical rates of binary neutron star coalescence give low prospects for detection of a single event by the initial LIGO/VIRGO interferometers. However, by utilizing at the best all the a priori information on the expected signal, a positive detection can be achieved. This relies on the hypothesis that $\gamma$-ray bursts are the electromagnetic signature of neutron star coalescences. The information about the direction of the source can then be used to add in phase the signals from different detectors in order (i) to increase the signal-to-noise ratio and (ii) to make the noise more Gaussian. Besides, the information about the time of arrival can be used to drastically decrease the observation time and thereby the false alarm rate. Moreover the fluence of the $\gamma$-ray emission gives some information about the amplitude of the gravitational signal. One can then add the signals from $\sim 10^4$ observation boxes ($\sim$ number of $\gamma$-ray bursts during 10 years) to yield a positive detection. Such a detection, based on the Maximum a Posteriori Probability Criterium, is a minimal one, in the sense that no information on the position and time of the events, nor on any parameter of the model, is collected. The advantage is that this detection requires an improvement of the detector sensitivity by a factor of only $\sim 1.5$ with respect to the initial LIGO/VIRGO interferometers, and that, if positive, it will confirm the $\gamma$-ray burst model.

An analytical scheme and a numerical method in order to study the effects of general relativity on the viscosity driven secular bar mode instability of rapidly rotating stars are presented. The approach consists in perturbing an axisymmetric and stationary configuration and studying its evolution by constructing a series of triaxial quasi-equilibrium configurations. These are obtained by solution of an approximate set of field equations where only the dominant non-axisymmetric terms are taken into account. The progress with respect to our former investigation consists in a higher relativistic order of the non-axisymmetric terms included into the computation, namely the fully three-dimensional treatment of the vector part of the space-time metric tensor as opposed to the scalar part, solely, in the former case. The scheme is applied to rotating stars built on a polytropic equation of state and compared to our previous results. The 3D-vector part turns out to inhibit the symmetry breaking efficiently. Nevertheless, the bar mode instability is still possible for an astrophysically relevant mass of M_ns=1.4 M_sun when a stiff polytropic equation of state with an adiabatic index of gamma=2.5 is employed. Triaxial neutron stars may be efficient emitters of gravitational waves and are thus potentially interesting sources for the forthcoming laser interferometric gravitational wave detectors such as LIGO, VIRGO and GEO600. From a numerical point of view, the solution of the three-dimensional minimal-distortion shift vector equation in spherical coordinates is an important achievement of our code.

A general relativistic version of the Euler equation for perfect fluid hydrodynamics is applied to a system of two neutron stars orbiting each other. In the quasi-equilibrium phase of the evolution of this system, a first integral of motion can be derived for certain velocity fields of the neutron star fluid including the (academic) case of co-rotation with respect to the orbital motion (synchronized binaries) and the realistic case of counter-rotation with respect to the orbital motion. The velocity field leading to this latter configuration can be computed by solving three-dimensional vector and scalar Poisson equations.

We study the possibility to detect the gravitational wave background generated by all the neutron stars in the Galaxy with only one gravitational wave interferometric detector. The proposed strategy consists in squaring the detector's output and searching for a sidereal modulation. The shape of the squared signal is computed for a disk and a halo distribution of neutron stars. The required noise stability of the interferometric detector is discussed. We argue that a possible population of old neutron stars, originating from a high stellar formation rate at the birth of the Galaxy and not emitting as radio pulsars, could be detected by the proposed technique in the low frequency range of interferometric experiments.

Triaxial neutron stars may be important sources of gravitational radiation for the forthcoming generation of interferometric gravitational wave detectors such as LIGO, VIRGO, and GEO600. We investigate the viscosity triggered bar mode secular instability of rapidly rotating neutron stars by means of a perturbation analysis of numerically constructed ``exact'' general relativistic axisymmetric star models. In the theoretical approach, only the dominant parts of the nonaxisymmetric terms of the 3D-Einstein equations are taken into account. A comparison of our results with previous studies of Newtonian polytropic stars confirms James' classical result gamma_crit=2.238 for the critical polytropic index. Beyond the Newtonian regime, gamma_crit reveals a slight increase toward highly relativistic configurations. Six out of twelve employed realistic dense matter equations of state admit the spontaneous symmetry breaking for masses above 1.6 M_sol.

We show that accreting neutron stars in binary systems or in Landau-Thorne-Zytkow objects are good candidates for continuous gravitational wave emission. Their gravitational radiation is strong enough to be detected by the next generation of detectors having a typical noise of 10^{-23} Hz^{-1/2}.

Contents: 1. Introduction 2. Spontaneous symmetry breaking 2.1 Review of classical results about Maclaurin/Jacobi ellipsoids 2.2 Spontaneous breaking of symmetry: a general phenomenom 2.3 Previous results for compressible Newtonian stars 2.4 Generation of gravitational waves 2.5 Finding the equilibrium configurations of a rotating star in the Newtonian regime 2.6 Extension to general relativity 2.7 First integral of fluid motion in general relativity 2.8 Gravitational field equations 2.9 Numerical results 3. CW emission from pulsars 4. Discussion and conclusion Appendix: Relativistic hydrodynamics in an accelerated frame

Continuous wave gravitational radiation from isolated rotating neutron stars is discussed. The general waveform and orders of magnitude for the amplitude are presented for various known pulsars. The specific case of gravitational radiation resulting from the distortion induced by the stellar magnetic field is presented. Finally some preliminary results about the signal from the whole population of neutron stars in the Galaxy are discussed.

Possible consequences of ferromagnetic transition in dense matter suggested recently by Kutschera and W{\'o}jcik, for the magnetic properties of neutron stars, are studied. Specific model of dense matter, in which a small admixture of protons is completely polarized due to their interaction with neutrons, is considered. Magnetic field of neutron stars with a ferromagnetic core is calculated within the framework of general relativity. Two types of boundary conditions at the ferromagnetic core edge are considered, corresponding to normal and superconducting liquid envelope, respectively. Numerical results for the neutron star magnetic dipole moment are confronted with pulsar timing. To be consistent with observations, ferromagnetic cores surrounded by a non-superconducting envelope, should consist of weakly ordered ferromagnetic domains. If domains are highly ordered, ferromagnetic core should be screened by a superconducting envelope.

Several applications of spectral methods to problems related to the relativistic astrophysics of compact objects are presented. Based on a proper definition of the analytical properties of regular tensorial functions we have developed a spectral method in a general sphericallike coordinate system. The applications include the investigation of spherically symmetric neutron star collapse as well as the solution of the coupled 2D-Einstein-Maxwell equations for magnetized, rapidly rotating neutron stars. In both cases the resulting codes are efficient and give results typically several orders of magnitude more accurate than equivalent codes based on finite difference schemes. We further report the current status of a 3D-code aiming at the simulation of non-axisymmetric neutron star collapse where we have chosen a tensor based numerical scheme.

The gravitational wave emission by a distorted rotating fluid star is computed. The distortion is supposed to be symmetric around some axis inclined with respect to the rotation axis. In the general case, the gravitational radiation is emitted at two frequencies: $\Omega$ and $2\Omega$, where $\Omega$ is the rotation frequency. The obtained formul\ae\ are applied to the specific case of a neutron star distorted by its own magnetic field. Assuming that the period derivative $\dot P$ of pulsars is a measure of their magnetic dipole moment, the gravitational wave amplitude can be related to the observable parameters $P$ and $\dot P$ and to a factor $\beta$ which measures the efficiency of a given magnetic dipole moment in distorting the star. $\beta$ depends on the nuclear matter equation of state and on the magnetic field distribution. The amplitude at the frequency $2\Omega$, expressed in terms of $P$, $\dot P$ and $\beta$, is independent of the angle $\alpha$ between the magnetic axis and the rotation axis, whereas at the frequency $\Omega$, the amplitude increases as $\alpha$ decreases. The value of $\beta$ for specific models of magnetic field distributions has been computed by means of a numerical code giving self-consistent models of magnetized neutron stars within general relativity. It is found that the distortion at fixed magnetic dipole moment is very dependent of the magnetic field distribution; a stochastic magnetic field or a superconductor stellar interior greatly increases $\beta$ with respect to the uniformly magnetized perfect conductor case and might lead to gravitational waves detectable by the VIRGO or LIGO interferometers. The amplitude modulation of the signal induced by the daily rotation of the Earth has been computed and specified to the case of the Crab pulsar and VIRGO

We investigate the effects of general relativity upon the non-axisymmetric ``bar'' mode secular instability of rapidly rotating stars, i.e. the relativistic and compressible analog of the transition from Maclaurin spheroids to Jacobi ellipsoids. Our method consists in perturbing a stationary axisymmetric configuration, constructed by a 2-D general relativistic numerical code, and taking into account only the dominant terms in the non-axisymmetric part of the 3-D relativistic equations. For a polytropic equation of state, we have determined, as a function of the degree of relativity, the critical adiabatic index $\gamma_{\rm crit}$ above which rapidly rotating stars can break their axial symmetry. A by-product of the present study is the confirmation of the Newtonian value $\gamma_{\rm crit} = 2.238$ obtained by James (1964). We have also considered neutron star models contructed upon twelve nuclear matter equations of state taken from the literature. We found that five equations of state from this sample allow the symmetry breaking for sufficiently high rotation velocities. For the others, the Keplerian velocity (mass-shedding from the equator) is reached before the axisymmetry is broken. Rotating neutron stars that break their axial symmetry can be an important source of gravitational waves for the LIGO/VIRGO interferometric detectors.

We present a formalism well adapted to the numerical study of the encounter of an ordinary main sequence star with a massive black hole. Symmetry considerations, the use of a well adapted moving grid and a well adapted moving frame along with integration of the partial differential equations by means of pseudo-spectral methods result in a very powerful and accurate tool. The hydrodynamical equations are written in a moving frame which mimics the bulk of the movement of the fluid, resulting in very small relative velocities and a well suited spatial resolution throughout the calculations. Therefore, the numerical calculations are considerably simplified, smoothing in particular the Courant condition. Typical runs are performed within a few hours on a workstation with the high accuracy linked to the spectral methods. Predictions of the so-called affine are tested against this full numerical simulation.

We present the first numerical solutions of the coupled Einstein-Maxwell equations describing rapidly rotating neutron stars endowed with a magnetic field. These solutions are fully relativistic and self-consistent, all the effects of the electromagnetic field on the star's equilibrium (Lorentz force, spacetime curvature generated by the electromagnetic stress-energy) being taken into account. The magnetic field is axisymmetric and poloidal. Five dense matter equations of state are employed. The partial differential equation system is integrated by means of a pseudo-spectral method. Various tests passed by the numerical code are presented. The effects of the magnetic field on neutron stars structure are then investigated, especially by comparing magnetized and non-magnetized configurations with the same baryon number. The deformation of the star induced by the magnetic field is important only for huge values of B (B>10^{10} T). The maximum mass as well as the maximum rotational velocity are found to increase with the magnetic field. The maximum allowable poloidal magnetic field is of the order of 10^{14} T (10^{18} G) and is reached when the magnetic pressure is comparable to the fluid pressure at the centre of the star. For such values, the maximum mass of neutron stars is found to increase by 13 to 29 % (depending upon the EOS) with respect to the maximum mass of non-magnetized stars.