Introduction to spectral methods

This proceeding is intended to be a first introduction to spectral methods. It is written around some simple problems that are solved explicitly and in details and that aim at demonstrating the power of those methods. The mathematical foundation of the spectral approximation is first introduced, based on the Gauss quadratures. The two usual basis of Legendre and Chebyshev polynomials are then presented. The next section is devoted to one dimensional equation solvers using only one domain. Three different methods are described. Techniques using several domains are shown in the last section of this paper and their various merits discussed

Circular geodesics and thick tori around rotating boson stars

Accretion disks play an important role in the evolution of their relativistic inner compact objects. The emergence of a new generation of interferometers will allow to resolve these accretion disks and provide more information about the properties of the central gravitating object. Due to this instrumental leap forward it is crucial to investigate the accretion disk physics near various types of inner compact objects now to deduce later constraints on the central objects from observations. A possible candidate for the inner object is the boson star. Here, we will try to analyze the differences between accretion structures surrounding boson stars and black holes. We aim at analysing the physics of circular geodesics around boson stars and study simple thick accretion tori (so-called Polish doughnuts) in the vicinity of these stars. We realize a detailed study of the properties of circular geodesics around boson stars. We then perform a parameter study of thick tori with constant angular momentum surrounding boson stars. This is done using the boson star models computed by a code constructed with the spectral solver library KADATH. We demonstrate that all the circular stable orbits are bound. In the case of a constant angular momentum torus, a cusp in the torus surface exists only for boson stars with a strong gravitational scalar field. Moreover, for each inner radius of the disk, the allowed specific angular momentum values lie within a constrained range which depends on the boson star considered. We show that the accretion tori around boson stars have different characteristics than in the vicinity of a black hole. With future instruments it could be possible to use these differences to constrain the nature of compact objects.Comment: Accepted for publication in CQ

Numerical simulation of oscillatons: extracting the radiating tail

Spherically symmetric, time-periodic oscillatons -- solutions of the Einstein-Klein-Gordon system (a massive scalar field coupled to gravity) with a spatially localized core -- are investigated by very precise numerical techniques based on spectral methods. In particular the amplitude of their standing-wave tail is determined. It is found that the amplitude of the oscillating tail is very small, but non-vanishing for the range of frequencies considered. It follows that exactly time-periodic oscillatons are not truly localized, and they can be pictured loosely as consisting of a well (exponentially) localized nonsingular core and an oscillating tail making the total mass infinite. Finite mass physical oscillatons with a well localized core -- solutions of the Cauchy-problem with suitable initial conditions -- are only approximately time-periodic. They are continuously losing their mass because the scalar field radiates to infinity. Their core and radiative tail is well approximated by that of time-periodic oscillatons. Moreover the mass loss rate of physical oscillatons is estimated from the numerical data and a semi-empirical formula is deduced. The numerical results are in agreement with those obtained analytically in the limit of small amplitude time-periodic oscillatons.Comment: 22 figures, accepted for publication in PR

Searching for Gravitational Waves from the Inspiral of Precessing Binary Systems: New Hierarchical Scheme using "Spiky" Templates

In a recent investigation of the effects of precession on the anticipated detection of gravitational-wave inspiral signals from compact object binaries with moderate total masses, we found that (i) if precession is ignored, the inspiral detection rate can decrease by almost a factor of 10, and (ii) previously proposed ``mimic'' templates cannot improve the detection rate significantly (by more than a factor of 2). In this paper we propose a new family of templates that can improve the detection rate by factors of 5--6 in cases where precession is most important. Our proposed method for these new ``mimic'' templates involves a hierarchical scheme of efficient, two-parameter template searches that can account for a sequence of spikes that appear in the residual inspiral phase, after one corrects for the any oscillatory modification in the phase. We present our results for two cases of compact object masses (10 and 1.4 solar masses and 7 and 3 solar masses) as a function of spin properties. Although further work is needed to fully assess the computational efficiency of this newly proposed template family, we conclude that these ``spiky templates'' are good candidates for a family of precession templates used in realistic searches, that can improve detection rates of inspiral events.Comment: 17 pages, 22 figures, version accepted by PRD. Minor revision

Searching for Gravitational Waves from the Inspiral of Precessing Binary Systems: Astrophysical Expectations and Detection Efficiency of "Spiky'' Templates

Relativistic spin-orbit and spin-spin couplings has been shown to modify the gravitational waveforms expected from inspiraling binaries with a black hole and a neutron star. As a result inspiral signals may be missed due to significant losses in signal-to-noise ratio, if precession effects are ignored in gravitational-wave searches. We examine the sensitivity of the anticipated loss of signal-to-noise ratio on two factors: the accuracy of the precessing waveforms adopted as the true signals and the expected distributions of spin-orbit tilt angles, given the current understanding of their physical origin. We find that the results obtained using signals generated by approximate techniques are in good agreement with the ones obtained by integrating the 2PN equations. This shows that a complete account of all high-order post-Newtonian effects is usually not necessary for the determination of detection efficiencies. Based on our current astrophysical expectations, large tilt angles are not favored and as a result the decrease in detection rate varies rather slowly with respect to the black hole spin magnitude and is within 20--30% of the maximum possible values.Comment: 7 fig., accepted by Phys. Rev. D Minor modification

A Fully Pseudospectral Scheme for Solving Singular Hyperbolic Equations

With the example of the spherically symmetric scalar wave equation on Minkowski space-time we demonstrate that a fully pseudospectral scheme (i.e. spectral with respect to both spatial and time directions) can be applied for solving hyperbolic equations. The calculations are carried out within the framework of conformally compactified space-times. In our formulation, the equation becomes singular at null infinity and yields regular boundary conditions there. In this manner it becomes possible to avoid "artificial" conditions at some numerical outer boundary at a finite distance. We obtain highly accurate numerical solutions possessing exponential spectral convergence, a feature known from solving elliptic PDEs with spectral methods. Our investigations are meant as a first step towards the goal of treating time evolution problems in General Relativity with spectral methods in space and time.Comment: 24 pages, 12 figure