Models of rotating relativistic stars with a toroidal magnetic field have been computed for a sample of eight equations of state of cold dense matter. Non-rotating models admit important levels of magnetization and quadrupole distortion accompanied by a seemingly unlimited growth in size. Rotating models reach the mass-shedding limit at smaller angular velocities than in the non-magnetized case according to the larger circumferential equatorial radius induced by the magnetic field. Moreover, they can be classified as prolate-prolate, oblate-prolate, or oblate-oblate with respect to surface deformation and quadrupole distortion. Simple expressions for surface and quadrupole deformation are provided that are valid up to magnetar field strengths and rapid rotation.

We have computed models of rotating relativistic stars with a toroidal magnetic field and investigated the combined effects of magnetic field and rotation on the apparent shape (i.e. the surface deformation), which could be relevant for the electromagnetic emission, and on the internal matter distribution (i.e. the quadrupole distortion), which could be relevant for the emission of gravitational waves. Using a sample of eight different cold nuclear physics equations of state, we have computed models of maximum field strength, as well as the distortion coefficients for the surface and the quadrupolar deformations. Surprisingly, we find that non-rotating models admit arbitrary levels of magnetization, accompanied by a growth of size and quadrupole distortion to which we could not find a limit. Rotating models, on the other hand, are subject to a mass-shedding limit at frequencies well below the corresponding ones for unmagnetized stars. Overall, the space of solutions can be split into three distinct classes for which the surface deformation and the quadrupole distortion are either prolate and prolate, oblate and prolate, or oblate and oblate, respectively. We also derive a simple formula expressing the relativistic distortion coefficients, which allows one to compute the surface deformation and the quadrupole distortion up to significant levels of rotation and magnetization, essentially covering all known magnetars. Such a formula replaces Newtonian equivalent expressions that overestimate the magnetic quadrupole distortion by about a factor of 6 and are inadequate for strongly relativistic objects like neutron stars.

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.

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.

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.

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.