We present a new class of high-order accurate numerical algorithms for solving the equations of general-relativistic ideal magnetohydrodynamics in curved spacetimes. In this paper we assume the background spacetime to be given and static, i.e., we make use of the Cowling approximation. The governing partial differential equations are solved via a new family of fully-discrete and arbitrary high-order accurate path-conservative discontinuous Galerkin (DG) finite-element methods combined with adaptive mesh refinement and time accurate local timestepping. In order to deal with shock waves and other discontinuities, the highorder DG schemes are supplemented with a novel a-posteriori subcell finite-volume limiter, which makes the new algorithms as robust as classical second-order total-variation diminishing finite-volume methods at shocks and discontinuities, but also as accurate as unlimited high-order DG schemes in smooth regions of the flow. We show the advantages of this new approach by means of various classical two- and three-dimensional benchmark problems on fixed spacetimes. Finally, we present a performance and accuracy comparisons between Runge-Kutta DG schemes and ADER high-order finite-volume schemes, showing the higher efficiency of DG schemes.

We present a strongly hyperbolic first-order formulation of the Einstein equations based on the conformal and covariant Z4 system (CCZ4) with constraint-violation damping, which we refer to as FO-CCZ4. As CCZ4, this formulation combines the advantages of a conformal and traceless formulation, with the suppression of constraint violations given by the damping terms, but being first order in time and space, it is particularly suited for a discontinuous Galerkin (DG) implementation. The strongly hyperbolic first-order formulation has been obtained by making careful use of first and second-order ordering constraints. A proof of strong hyperbolicity is given for a selected choice of standard gauges via an analytical computation of the entire eigenstructure of the FO-CCZ4 system. The resulting governing partial differential equations system is written in non-conservative form and requires the evolution of 58 unknowns. A key feature of our formulation is that the first-order CCZ4 system decouples into a set of pure ordinary differential equations and a reduced hyperbolic system of partial differential equations that contains only linearly degenerate fields. We implement FO-CCZ4 in a high-order path-conservative arbitrary-high-order-method-using-derivatives (ADER)-DG scheme with adaptive mesh refinement and local time-stepping, supplemented with a third-order ADER-WENO subcell finite-volume limiter in order to deal with singularities arising with black holes. We validate the correctness of the formulation through a series of standard tests in vacuum, performed in one, two and three spatial dimensions, and also present preliminary results on the evolution of binary black-hole systems. To the best of our knowledge, these are the first successful three-dimensional simulations of moving punctures carried out with high-order DG schemes using a first-order formulation of the Einstein equations.

This paper presents an arbitrary h.o. accurate ADER DG method on space-time adaptive meshes (AMR) for the solution of two important families of non-linear time dependent PDE for compr. dissipative flows: the compr. Navier-Stokes equations and the equations of visc. and res. MHD in 2 and 3 space-dimensions. The work continues a recent series of papers concerning the development and application of a proper a posteriori subcell FV limiting procedure suitable for DG methods. It is a well known fact that a major weakness of h.o. DG methods lies in the difficulty of limiting discontinuous solutions, which generate spurious oscillations, namely the so-called 'Gibbs phenomenon'. In the present work the main benefits of the MOOD paradigm, i.e. the computational robustness even in the presence of strong shocks, are preserved and the numerical diffusion is considerably reduced also for the limited cells by resorting to a proper sub-grid. An important feature of our new scheme is its ability to cure even floating point errors that may occur during a simulation, for example when taking real roots of negative numbers or after divisions by zero. We apply the whole approach for the first time to the equations of compr. gas dynamics and MHD in the presence of viscosity, thermal conductivity and magnetic resistivity, therefore extending our family of adaptive ADER-DG schemes to cases for which the numerical fluxes also depend on the gradient of the state vector. The distinguished high-resolution properties of the presented numerical scheme stands out against a wide number of non-trivial test cases both for the compr. Navier-Stokes and the viscous and resistive MHD equations. The present results show clearly that the shock-capturing capability of the news schemes are significantly enhanced within a cell-by-cell Adaptive Mesh Refinement implementation together with time accurate local time stepping (LTS).

We present a new version of conservative ADER-WENO finite volume schemes, in which both the high order spatial reconstruction as well as the time evolution of the reconstruction polynomials in the local space-time predictor stage are performed in primitive variables, rather than in conserved ones. Since the underlying finite volume scheme is still written in terms of cell averages of the conserved quantities, our new approach performs the spatial WENO reconstruction twice: the first WENO reconstruction is carried out on the known cell averages of the conservative variables. The WENO polynomials are then used at the cell centers to compute point values of the conserved variables, which are converted into point values of the primitive variables. A second WENO reconstruction is performed on the point values of the primitive variables to obtain piecewise high order reconstruction polynomials of the primitive variables. The reconstruction polynomials are subsequently evolved in time with a novel space-time finite element predictor that is directly applied to the governing PDE written in primitive form. We have verified the validity of the new approach over the classical Euler equations of gas dynamics, the special relativistic hydrodynamics (RHD) and ideal magnetohydrodynamics (RMHD) equations, as well as the Baer-Nunziato model for compressible two-phase flows. In all cases we have noticed that the new ADER schemes provide less oscillatory solutions when compared to ADER finite volume schemes based on the reconstruction in conserved variables, especially for the RMHD and the Baer-Nunziato equations. For the RHD and RMHD equations, the accuracy is improved and the CPU time is reduced by about 25%. We recommend to use this version of ADER as the standard one in the relativistic framework. The new approach can be extended to ADER-DG schemes on space-time adaptive grids.

We study the dynamics of radiation pressure supported tori around Schwarzschild black holes, focusing on their oscillatory response to an external perturbation. Using KORAL, a general relativistic radiation hydrodynamics code capable of modeling all radiative regimes from the optically thick to the optically thin, we monitor a sample of models at different initial temperatures and opacities, evolving them in two spatial dimensions for $\sim 165$ orbital periods. The dynamics of models with high opacity is very similar to that of purely hydrodynamics models, and it is characterized by regular oscillations which are visible also in the light curves. As the opacity is decreased, the tori quickly and violently migrate towards the gas-pressure dominated regime, collapsing towards the equatorial plane. When the spectra of the $L_2$ norm of the mass density are considered, high frequency inertial-acoustic modes of oscillations are detected (with the fundamental mode at a frequency $68 M_{\rm BH}^{-1}\,\rm Hz$), in close analogy to the phenomenology of purely hydrodynamic models. An additional mode of oscillation, at a frequency $129 M_{\rm BH}^{-1}\,\rm Hz$, is also found, which can be unambiguously attributed to the radiation. The spectra extracted from the light curves are typically more noisy, indicating that in a real observation such modes would not be easily detected.

In this paper we present a novel arbitrary high order accurate discontinuous Galerkin (DG) finite element method on space-time adaptive Cartesian meshes (AMR) for hyperbolic conservation laws in multiple space dimensions, using a high order \aposteriori sub-cell ADER-WENO finite volume \emph{limiter}. Notoriously, the original DG method produces strong oscillations in the presence of discontinuous solutions and several types of limiters have been introduced over the years to cope with this problem. Following the innovative idea recently proposed in \cite{Dumbser2014}, the discrete solution within the troubled cells is \textit{recomputed} by scattering the DG polynomial at the previous time step onto a suitable number of sub-cells along each direction. Relying on the robustness of classical finite volume WENO schemes, the sub-cell averages are recomputed and then gathered back into the DG polynomials over the main grid. In this paper this approach is implemented for the first time within a space-time adaptive AMR framework in two and three space dimensions, after assuring the proper averaging and projection between sub-cells that belong to different levels of refinement. The combination of the sub-cell resolution with the advantages of AMR allows for an unprecedented ability in resolving even the finest details in the dynamics of the fluid. The spectacular resolution properties of the new scheme have been shown through a wide number of test cases performed in two and in three space dimensions, both for the Euler equations of compressible gas dynamics and for the magnetohydrodynamics (MHD) equations.

We present a new numerical tool for solving the special relativistic ideal MHD equations that is based on the combination of the following three key features: (i) a one-step ADER discontinuous Galerkin (DG) scheme that allows for an arbitrary order of accuracy in both space and time, (ii) an a posteriori subcell finite volume limiter that is activated to avoid spurious oscillations at discontinuities without destroying the natural subcell resolution capabilities of the DG finite element framework and finally (iii) a space-time adaptive mesh refinement (AMR) framework with time-accurate local time-stepping. The divergence-free character of the magnetic field is instead taken into account through the so-called "divergence-cleaning" approach. The convergence of the new scheme is verified up to 5th order in space and time and the results for a set of significant numerical tests including shock tube problems, the RMHD rotor and blast wave problems, as well as the Orszag-Tang vortex system are shown. We also consider a simple case of the relativistic Kelvin-Helmholtz instability with a magnetic field, emphasizing the potential of the new method for studying turbulent RMHD flows. We discuss the advantages of our new approach when the equations of relativistic MHD need to be solved with high accuracy within various astrophysical systems.

We present a new code for performing general-relativistic radiation-hydrodynamics simulations of accretion flows onto black holes. The radiation field is treated in the optically-thick approximation, with the opacity contributed by Thomson scattering and thermal bremsstrahlung. Our analysis is concentrated on a detailed numerical investigation of hot two-dimensional, Bondi-Hoyle accretion flows with various Mach numbers. We find significant differences with respect to purely hydrodynamical evolutions. In particular, once the system relaxes to a radiation-pressure dominated regime, the accretion rates become about two orders of magnitude smaller than in the purely hydrodynamical case, remaining however super-Eddington as are the luminosities. Furthermore, when increasing the Mach number of the inflowing gas, the accretion rates become smaller because of the smaller cross section of the black hole, but the luminosities increase as a result a stronger emission in the shocked regions. Overall, our approach provides the first self-consistent calculation of the Bondi-Hoyle luminosity, most of which is emitted within r~100 M from the black hole, with typical values L/L_Edd ~ 1-7, and corresponding energy efficiencies eta_BH ~ 0.09-0.5. The possibility of computing luminosities self-consistently has also allowed us to compare with the bremsstrahlung luminosity often used in modelling the electromagnetic counterparts to supermassive black-hole binaries, to find that in the optically-thick regime these more crude estimates are about 20 times larger than our radiation-hydrodynamics results.

We propose a simple model for an optically thick radiative torus in local thermodynamic equilibrium around a Kerr black hole. The hydrodynamics structure, which is not affected by the radiation field, is the same as for the so--called polish doughnuts. Under the assumption of isentropic fluid and polytropic equation of state, a simple stationary and axisymmetric solution to the relativistic radiation hydrodynamics equations is possible, for which the temperature of the torus scales like the specific enthalpy. The astrophysical relevance of the model is briefly discussed.

We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method. Due to the one-step nature of the underlying scheme, the resulting algorithm is particularly well suited for an AMR strategy on space-time adaptive meshes, i.e.with time-accurate local time stepping. The AMR property has been implemented 'cell-by-cell', with a standard tree-type algorithm, while the scheme has been parallelized via the Message Passing Interface (MPI) paradigm. The new scheme has been tested over a wide range of examples for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations of compressible gas dynamics and the equations of magnetohydrodynamics (MHD). High order in space and time have been confirmed via a numerical convergence study and a detailed analysis of the computational speed-up with respect to highly refined uniform meshes is also presented. We also show test problems where the presented high order AMR scheme behaves clearly better than traditional second order AMR methods. The proposed scheme that combines for the first time high order ADER methods with space--time adaptive grids in two and three space dimensions is likely to become a useful tool in several fields of computational physics, applied mathematics and mechanics.

We revisit particle acceleration in the polar cap region of a neutron star by taking into account both general relativistic effects and the presence of toroidal oscillations at the star surface. In particular, we address the question of whether toroidal oscillations at the stellar surface can affect the acceleration properties in the polar cap. We solve numerically the relativistic electrodynamics equations in the stationary regime, focusing on the computation of the Lorentz factor of a space-charge-limited electron flow accelerated in the polar cap region of a rotating and oscillating pulsar. To this extent, we adopt the correct expression of the general relativistic Goldreich-Julian charge density in the presence of toroidal oscillations. Depending on the ratio of the actual charge density of the pulsar magnetosphere to the Goldreich-Julian charge density, we distinguish two different regimes of the Lorentz factor of the particle flow, namely an oscillatory regime produced for sub-GJ current density configurations, which does not produce an efficient acceleration, and a true accelerating regime for super-GJ current density configurations. We find that star oscillations may be responsible for a significant asymmetry in the pulse profile that depends on the orientation of the oscillations with respect to the pulsar magnetic field. In particular, significant enhancements of the Lorentz factor are produced by stellar oscillations in the super-GJ current density regime.

We analyze a class of physical properties, forming the content of the so-called von Zeipel theorem, which characterizes stationary, axisymmetric, non-selfgravitating perfect fluids in circular motion in the gravitational field of a compact object. We consider the extension of the theorem to the magnetohydrodynamic regime, under the assumption of an infinitely conductive fluid, both in the Newtonian and in the relativistic framework. When the magnetic field is toroidal, the conditions required by the theorem are equivalent to integrability conditions, as it is the case for purely hydrodynamic flows. When the magnetic field is poloidal, the analysis for the relativistic regime is substantially different with respect to the Newtonian case and additional constraints, in the form of PDEs, must be imposed on the magnetic field in order to guarantee that the angular velocity $\Omega$ depends only on the specific angular momentum $\ell$. In order to deduce such physical constraints, it is crucial to adopt special coordinates, which are adapted to the $\Omega={\rm const}$ surfaces. The physical significance of these results is briefly discussed.

We present a high order one-step ADER-WENO finite volume scheme with space-time adaptive mesh refinement (AMR) for the solution of the special relativistic hydrodynamic and magnetohydrodynamic equations. By adopting a local discontinuous Galerkin predictor method, a high order one-step time discretization is obtained, with no need for Runge--Kutta sub-steps. This turns out to be particularly advantageous in combination with space-time adaptive mesh refinement, which has been implemented following a "cell-by-cell" approach. As in existing second order AMR methods, also the present higher order AMR algorithm features time-accurate local time stepping (LTS), where grids on different spatial refinement levels are allowed to use different time steps. We also compare two different Riemann solvers for the computation of the numerical fluxes at the cell interfaces. The new scheme has been validated over a sample of numerical test problems in one, two and three spatial dimensions, exploring its ability in resolving the propagation of relativistic hydrodynamical and magnetohydrodynamical waves in different physical regimes. The astrophysical relevance of the new code for the study of the Richtmyer--Meshkov instability is briefly discussed in view of future applications.

We study the Richtmyer--Meshkov (RM) instability of a relativistic perfect fluid by means of high order numerical simulations with adaptive mesh refinement (AMR). The numerical scheme adopts a finite volume Weighted Essentially Non-Oscillatory (WENO) reconstruction to increase accuracy in space, a local space-time discontinuous Galerkin predictor method to obtain high order of accuracy in time and a high order one-step time update scheme together with a "cell-by-cell" space-time AMR strategy with time-accurate local time stepping. In this way, third order accurate (both in space and in time) numerical simulations of the RM instability are performed, spanning a wide parameter space. We present results both for the case in which a light fluid penetrates into a higher density one (Atwood number $A>0$), and for the case in which a heavy fluid penetrates into a lower density one (Atwood number $A<0$). We find that, for large Lorentz factors \gamma_s of the incident shock wave, the relativistic RM instability is substantially weakened and ultimately suppressed. More specifically, the growth rate of the RM instability in the linear phase has a local maximum which occurs at a critical value of \gamma_s ~ [1.2,2]. Moreover, we have also revealed a genuine relativistic effect, absent in Newtonian hydrodynamics, which arises in three dimensional configurations with a non-zero velocity component tangent to the incident shock front. In this case, the RM instability is strongly affected, typically resulting in less efficient mixing of the fluid.

We try to explain the subpulse drift phenomena adopting the space-charge limited flow (SCLF) model and comparing the plasma drift velocity in the inner region of pulsar magnetospheres with the observed velocity of drifting subpulses. We apply the approach described in a recent paper of van Leeuwen & Timokhin (2012), where it was shown that the standard estimation of the subpulse drift velocity through the total value of the scalar potential drop in the inner gap gives inaccurate results, while the exact expression relating the drift velocity to the gradient of the scalar potential should be used instead. After considering a selected sample of sources taken from the catalog of Weltevrede, Edwards & Stappers (2006) with coherently drifting subpulses and reasonably known observing geometry, we show that their subpulse drift velocities would correspond to the drift of the plasma located very close or above the pair formation front. Moreover, a detailed analysis of PSR B0826-34 and PSR B0818-41 reveals that the variation of the subpulse separation with the pulse longitude can be successfully explained by the dependence of the plasma drift velocity on the angular coordinates.

A high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions is presented. A high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method, while a high order spatial accuracy is obtained through a WENO reconstruction. Thanks to the one-step nature of the underlying scheme, the resulting algorithm can be efficiently imported within an AMR framework on space-time adaptive meshes. We provide convincing evidence that the presented high order AMR scheme behaves better than traditional second order AMR methods. Tests are shown of the new scheme for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations and the equations of ideal magnetohydrodynamics. The proposed scheme is likely to become a useful tool in several astrophysical scenarios.

We present the implementation of an implicit-explicit (IMEX) Runge-Kutta numerical scheme for general relativistic hydrodynamics coupled to an optically thick radiation field in two existing GR-hydrodynamics codes. We argue that the necessity of such an improvement arises naturally in astrophysically relevant regimes where the optical thickness is high as the equations become stiff. By performing several 1D tests we verify the codes' new ability to deal with this stiffness and show consistency. Then, still in 1D, we compute a luminosity versus accretion rate diagram for the setup of spherical accretion onto a Schwarzschild black hole and find good agreement with previous work. Lastly, we revisit the supersonic Bondi Hoyle Lyttleton (BHL) accretion in 2D where we can now present simulations of realistic temperatures, down to T~10^6 K. Here we find that radiation pressure plays an important role, but also that these highly dynamical set-ups push our approximate treatment towards the limit of physical applicability. The main features of radiation hydrodynamics BHL flows manifest as (i) an effective adiabatic index approaching gamma_effective ~ 4/3; (ii) accretion rates two orders of magnitude lower than without radiation pressure; (iii) luminosity estimates around the Eddington limit, hence with an overall radiative efficiency as small as eta ~ 10^{-2}; (iv) strong departures from thermal equilibrium in shocked regions; (v) no appearance of the flip-flop instability. We conclude that the current optically thick approximation to the radiation transfer does give physically substantial improvements over the pure hydro also in set-ups departing from equilibrium, and, once accompanied by an optically thin treatment, is likely to provide a fundamental tool for investigating accretion flows in a large variety of astrophysical systems.

We revisit the suggestion that dual jets can be produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. By performing independent calculations and by computing the electromagnetic emission in a way which is consistent with estimates using the Poynting flux, we show that a dual-jet structure is present but energetically subdominant with respect to a non-collimated and predominantly quadrupolar emission, which is similar to the one computed when the binary is in electrovacuum. While our findings set serious restrictions on the detectability of dual jets from coalescing binaries, they also increase the chances of detecting an EM counterpart from these systems.

We investigate the conditions for radio emission in rotating and oscillating magnetars, by focusing on the main physical processes determining the position of their death-lines in the P-\dot{P} diagram, i.e. of those lines that separate the regions where the neutron star may be radio-loud or radio-quiet. After using the general relativistic expression for the electromagnetic scalar potential in the magnetar magnetosphere, we find that larger compactness parameters of the star as well as larger inclination angles between the rotation axis and the magnetic moment produce death-lines well above the majority of known magnetars. This is consistent with the observational evidence of no regular radio emission from the magnetars in the frequency range typical for the ordinary pulsars. On the contrary, when oscillations of the magnetar are taken into account, the death-lines shift downward and the conditions necessary for the generation of radio emission in the magnetosphere are met. Present observations showing a close connection between the burst activity of magnetars and the generation of the radio emission in the magnetar magnetosphere are naturally accounted for within our interpretation.

We perform a global linear perturbative analysis, and investigate the oscillation properties of relativistic, non-selfgravitating tori orbiting around neutron stars in the slow rotation limit approximation. Extending the work done in Schwarzschild and Kerr backgrounds, we consider the axisymmetric oscillations of vertically integrated tori in the Hartle-Thorne spacetime. The equilibrium models are constructed by selecting a number of different non-Keplerian distributions of specific angular momentum, allowing for disc sizes $L \sim 0.5 - 600$ gravitational radii. Our results, obtained after solving a global eigenvalue problem to compute the xisymmetric $p$-modes, indicate that such oscillation modes could account with most observed lower ($\nu_L$) and upper ($\nu_U$) high frequency quasi-periodic oscillations for Sco X-1, and for some Z sources and Atoll sources with $\nu_L\ gtrsim 500$ Hz. However, when $\nu_L \lesssim 500$ Hz, $p$-modes oscillations do not account for the linear relation $\nu_U = A \nu_L + B$, $B \neq 0$ between the upper and the lower high frequency quasi-periodic oscillations that are observed in neutron star low-mass X-ray binaries.

We present the results of two-dimensional magnetohydrodynamical numerical simulations of relativistic magnetic reconnection, with particular emphasis on the dynamics of Petschek-type configurations with high Lundquist numbers, S ~ 10^5-10^8. The numerical scheme adopted, allowing for unprecedented accuracy for this type of calculations, is based on high order finite volume and discontinuous Galerkin methods as recently proposed by Dumbser & Zanotti (2009). The possibility of producing high Lorentz factors is discussed, by studying the effects produced on the dynamics by different magnetization and resistivity regimes. We show that Lorentz factors close to ~4 can be produced for a plasma magnetization parameter sigma=20. Moreover, we find that the Sweet-Parker layers are unstable, generating secondary magnetic islands, but only for S>S_c~10^8, much larger than what is reported in the Newtonian regime.

We show the results of two dimensional general relativistic inviscid and isothermal hydrodynamical simulations comparing the behavior of co-rotating (with respect to the black hole rotation) and counter-rotating circumbinary quasi-Keplerian discs in the post merger phase of a supermassive binary black hole system. While confirming the spiral shock generation within the disc due to the combined effects of mass loss and recoil velocity of the black hole, we find that the maximum luminosity of counter-rotating discs is a factor ~(2-12) higher than in the co-rotating case, depending on the spin of the black hole. On the other hand, the luminosity peak happens ~10 days later with respect to the co-rotating case, for a binary with a total mass M~10^6 M_\odot. Although the global dynamics of counter-rotating discs in the post merger phase of a merging event is very similar to that for co-rotating discs, an important difference has been found. In fact, increasing the spin of the central black hole produces more luminous co-rotating discs while less luminous counter-rotating ones.

We present the results of two-dimensional and three-dimensional magnetohydrodynamical numerical simulations of relativistic magnetic reconnection, with particular emphasis on the dynamics of the plasma in a Petschek-type configuration with high Lundquist numbers, S\sim 10^5-10^8. The numerical scheme adopted, allowing for unprecedented accuracy for this type of calculations, is based on high order finite volume and discontinuous Galerkin methods as recently proposed by \citet{Dumbser2009}. The possibility of producing high Lorentz factors is discussed, showing that Lorentz factors close to \sim 4 can be produced for a plasma parameter \sigma_m=20. Moreover, we find that the Sweet-Parker layers are unstable, generating secondary magnetic islands, but only for S > S_c = 10^8, much larger than what is reported in the Newtonian regime. Finally, the effects of a mildly anisotropic Ohm law are considered in a configuration with a guide magnetic field. Such effects produce only slightly faster reconnection rates and Lorentz factors of about 1% larger with respect to the perfectly isotropic Ohm law.

The numerical investigation of Bondi-Hoyle accretion onto a moving black hole has a long history, both in Newtonian and in general-relativistic physics. By performing new two-dimensional and general-relativistic simulations onto a rotating black hole, we point out a novel feature, namely, that quasi-periodic oscillations (QPOs) are naturally produced in the shock cone that develops in the downstream part of the flow. Because the shock cone in the downstream part of the flow acts as a cavity trapping pressure perturbations, modes with frequencies in the integer ratios 2:1 and 3:1 are easily produced. The frequencies of these modes depend on the black-hole spin and on the properties of the flow, and scale linearly with the inverse of the black-hole mass. Our results may be relevant for explaining the detection of QPOs in Sagittarius A*, once such detection is confirmed by further observations. Finally, we report on the development of the flip-flop instability, which can affect the shock cone under suitable conditions; such an instability has been discussed before in Newtonian simulations but was never found in a relativistic regime.

We investigate the dynamics of a circumbinary disc that responds to the loss of mass and to the recoil velocity of the black hole produced by the merger of a binary system of supermassive black holes. We perform the first two-dimensional general relativistic hydrodynamics simulations of \textit{extended} non-Keplerian discs and employ a new technique to construct a "shock detector", thus determining the precise location of the shocks produced in the accreting disc by the recoiling black hole. In this way we can study how the properties of the system, such as the spin, mass and recoil velocity of the black hole, affect the mass accretion rate and are imprinted on the electromagnetic emission from these sources. We argue that the estimates of the bremsstrahlung luminosity computed without properly taking into account the radiation transfer yield cooling times that are unrealistically short. At the same time we show, through an approximation based on the relativistic isothermal evolution, that the luminosity produced can reach a peak value above $L \simeq 10^{43} \ {\rm erg/s} $ at about $\sim 30\,{\rm d}$ after the merger of a binary with total mass $M\simeq 10^6 M_\odot$ and persist for several days at values which are a factor of a few smaller. If confirmed by more sophisticated calculations such a signal could indeed lead to an electromagnetic counterpart of the merger of binary black-hole system.

We study the magnetosphere of a slowly rotating magnetized neutron star subject to toroidal oscillations in the relativistic regime. Under the assumption of a zero inclination angle between the magnetic moment and the angular momentum of the star, we analyze the Goldreich-Julian charge density and derive a second-order differential equation for the electrostatic potential. The analytical solution of this equation in the polar cap region of the magnetosphere shows the modification induced by stellar toroidal oscillations on the accelerating electric field and on the charge density. We also find that, after decomposing the oscillation velocity in terms of spherical harmonics, the first few modes with $m=0,1$ are responsible for energy losses that are almost linearly dependent on the amplitude of the oscillation and that, for the mode $(l,m)=(2,1)$, can be a factor $\sim8$ larger than the rotational energy losses, even for a velocity oscillation amplitude at the star surface as small as $\eta=0.05 \ \Omega \ R$. The results obtained in this paper clarify the extent to which stellar oscillations are reflected in the time variation of the physical properties at the surface of the rotating neutron star, mainly by showing the existence of a relation between $P\dot{P}$ and the oscillation amplitude. Finally, we propose a qualitative model for the explanation of the phenomenology of intermittent pulsars in terms of stellar oscillations that are periodically excited by star glitches.

In this paper we propose the first better than second order accurate method in space and time for the numerical solution of the resistive relativistic magnetohydrodynamics (RRMHD) equations on unstructured meshes in multiple space dimensions. The nonlinear system under consideration is purely hyperbolic and contains a source term, the one for the evolution of the electric field, that becomes stiff for low values of the resistivity. For the spatial discretization we propose to use high order $\PNM$ schemes as introduced in \cite{Dumbser2008} for hyperbolic conservation laws and a high order accurate unsplit time discretization is achieved using the element-local space-time discontinuous Galerkin approach proposed in \cite{DumbserEnauxToro} for one-dimensional balance laws with stiff source terms. The divergence free character of the magnetic field is accounted for through the divergence cleaning procedure of Dedner et al. \cite{Dedneretal}. To validate our high order method we first solve some numerical test cases for which exact analytical reference solutions are known and we also show numerical convergence studies in the stiff limit of the RRMHD equations using $\PNM$ schemes from third to fifth order of accuracy in space and time. We also present some applications with shock waves such as a classical shock tube problem with different values for the conductivity as well as a relativistic MHD rotor problem and the relativistic equivalent of the Orszag-Tang vortex problem. We have verified that the proposed method can handle equally well the resistive regime and the stiff limit of ideal relativistic MHD. For these reasons it provides a powerful tool for relativistic astrophysical simulations involving the appearance of magnetic reconnection.

We present a new numerical code, ECHO, based on an Eulerian Conservative High Order scheme for time dependent three-dimensional general relativistic magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at providing a shock-capturing conservative method able to work at an arbitrary level of formal accuracy (for smooth flows), where the other existing GRMHD and GRMD schemes yield an overall second order at most. Moreover, our goal is to present a general framework, based on the 3+1 Eulerian formalism, allowing for different sets of equations, different algorithms, and working in a generic space-time metric, so that ECHO may be easily coupled to any solver for Einstein's equations. Various high order reconstruction methods are implemented and a two-wave approximate Riemann solver is used. The induction equation is treated by adopting the Upwind Constrained Transport (UCT) procedures, appropriate to preserve the divergence-free condition of the magnetic field in shock-capturing methods. The limiting case of magnetodynamics (also known as force-free degenerate electrodynamics) is implemented by simply replacing the fluid velocity with the electromagnetic drift velocity and by neglecting the matter contribution to the stress tensor. ECHO is particularly accurate, efficient, versatile, and robust. It has been tested against several astrophysical applications, including a novel test on the propagation of large amplitude circularly polarized Alfven waves. In particular, we show that reconstruction based on a Monotonicity Preserving filter applied to a fixed 5-point stencil gives highly accurate results for smooth solutions, both in flat and curved metric (up to the nominal fifth order), while at the same time providing sharp profiles in tests involving discontinuities.

We present a numerical study of the dynamics of magnetized, relativistic, non-self-gravitating, axisymmetric tori orbiting in the background spacetimes of Schwarzschild and Kerr black holes. The initial models have a constant specific angular momentum and are built with a non-zero toroidal magnetic field component, for which equilibrium configurations have recently been obtained. In this work we extend our previous investigations which dealt with purely hydrodynamical thick discs, and study the dynamics of magnetized tori subject to perturbations which, for the values of the magnetic field strength considered here, trigger quasi-periodic oscillations lasting for tens of orbital periods. Overall, we have found that the dynamics of the magnetized tori analyzed is very similar to that found in the corresponding unmagnetized models. The spectral distribution of the eigenfrequencies of oscillation shows the presence of a fundamental p mode and of a series of overtones in a harmonic ratio 2:3:.... These simulations, therefore, extend the validity of the model of Rezzolla et al.(2003a) for explaining the high-frequency QPOs observed in the spectra of LMXBs containing a black-hole candidate also to the case of magnetized discs with purely toroidal magnetic field distribution. If sufficiently compact and massive, these oscillations can also lead to the emission of intense gravitational radiation which is potentially detectable for sources within the Galaxy.

By combining the numerical solution of the nonlinear hydrodynamics equations with the solution of the linear inhomogeneous Zerilli-Moncrief and Regge-Wheeler equations we investigate the properties of the gravitational radiation emitted during the axisymmetric accretion of matter onto a Schwarzschild black hole. The matter models considered include quadrupolar dust shells and thick accretion disks, permitting us to simulate situations which may be encountered at the end stages of stellar gravitational collapse or binary neutron star merger. We focus on the interference pattern appearing in the energy spectra of the emitted gravitational waves and on the amount of excitation of the quasi-normal modes of the accreting black hole. We show that, quite generically in the presence of accretion, the black hole ringdown is not a simple superposition of quasi-normal modes, although the fundamental mode is usually present and often dominates the gravitational-wave signal. We interpret this as due to backscattering of waves off the non-exponentially decaying part of the black-hole potential and to the finite spatial extension of the accreting matter. Our results suggest that the black-hole QNM contributions to the full gravitational-wave signal should be extremely small and possibly not detectable in generic astrophysical scenarios involving the accretion of extended distributions of matter.

Quasi-periodic oscillations of high density thick accretion disks orbiting a Schwarzschild black hole have been recently addressed as interesting sources of gravitational waves. The aim of this paper is to compare the gravitational waveforms emitted from these sources when computed using (variations of) the standard quadrupole formula and gauge-invariant metric perturbation theory. To this goal we evolve representative disk models using an existing general relativistic hydrodynamics code which has been previously employed in investigations of such astrophysical systems. Two are the main results of this work: First, for stable and marginally stable disks, no excitation of the black hole quasi-normal modes is found. Secondly, we provide a simple, relativistic modification of the Newtonian quadrupole formula which, in certain regimes, yields excellent agreement with the perturbative approach. This holds true as long as back-scattering of GWs is negligible. Otherwise, any functional form of the quadrupole formula yields systematic errors of the order of 10%.

We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbolic character in the solution procedure. All theoretical ingredients necessary to build up high-resolution shock-capturing schemes based on the solution of local Riemann problems (i.e. Godunov-type schemes) are described. In particular, we use a renormalized set of regular eigenvectors of the flux Jacobians of the relativistic magnetohydrodynamics equations. In addition, the paper describes a procedure based on the equivalence principle of general relativity that allows the use of Riemann solvers designed for special relativistic magnetohydrodynamics in GRMHD. Our formulation and numerical methodology are assessed by performing various test simulations recently considered by different authors. These include magnetized shock tubes, spherical accretion onto a Schwarzschild black hole, equatorial accretion onto a Kerr black hole, and magnetized thick accretion disks around a black hole prone to the magnetorotational instability.

We present a comprehensive numerical study of the dynamics of relativistic axisymmetric accretion tori with a power-law distribution of specific angular momentum orbiting in the background spacetime of a Kerr black hole. By combining general relativistic hydrodynamics simulations with a linear perturbative approach we investigate the main dynamical properties of these objects over a large parameter space. The astrophysical implications of our results extend and improve two interesting results that have been recently reported in the literature. Firstly, the induced quasi-periodic variation of the mass quadrupole moment makes relativistic tori of nuclear matter densities, as those formed during the last stages of binary neutron star mergers, promising sources of gravitational radiation, potentially detectable by interferometric instruments. Secondly, $p$-mode oscillations in relativistic tori of low rest-mass densities could be used to explain high frequency quasi-periodic oscillations observed in X-ray binaries containing a black hole candidate under conditions more generic than those considered so far.

We consider the effects of a cosmological constant on the dynamics of constant angular momentum discs orbiting Schwarzschild-de Sitter black holes. The motivation behind this study is to investigate whether the presence of a radial force contrasting the black hole's gravitational attraction can influence the occurrence of the runaway instability, a robust feature of the dynamics of constant angular momentum tori in Schwarzschild and Kerr spacetimes. In addition to the inner cusp near the black hole horizon through which matter can accrete onto the black hole, in fact, a positive cosmological constant introduces also an outer cusp through which matter can leave the torus without accreting onto the black hole. To assess the impact of this outflow on the development of the instability we have performed time-dependent and axisymmetric hydrodynamical simulations of equilibrium initial configurations in a sequence of background spacetimes of Schwarzschild-de Sitter black holes with increasing masses. The simulations have been performed with an unrealistic value for the cosmological constant which, however, yields sufficiently small discs to be resolved accurately on numerical grids and thus provides a first qualitative picture of the dynamics. The calculations, carried out for a wide range of initial conditions, show that the mass-loss from the outer cusp can have a considerable impact on the instability, with the latter being rapidly suppressed if the outflow is large enough.

This is the first of a series of papers investigating the oscillation properties of relativistic, non-selfgravitating tori orbiting around a black hole. In this initial paper we consider the axisymmetric oscillation modes of a torus constructed in a Schwarzschild spacetime. To simplify the treatment and make it as analytical as possible, we build our tori with vertically integrated and vertically averaged quantities. Our investigation progresses by steps. We first consider a local analysis in Newtonian gravity and determine the properties of acoustic wave propagation within these objects, as well as the relations between acoustic and epicyclic oscillations. Next, we extend the local analysis to a general relativistic framework. Finally, we perform a global analysis and determine both the eigenfunctions and the eigenfrequencies of the axisymmetric oscillations corresponding to the p modes of relativistic tori. These behave as sound waves globally trapped in the torus and possess eigenfrequencies appearing in the simple sequence 2:3:4:..., independently of the distribution of angular momentum considered. The properties of the modes investigated here are in good agreement with those observed in recent numerical simulations and could have a number of different applications. In X-ray binary systems containing a black hole candidate, for instance, p-mode oscillations could be used to explain the harmonic relations in the high frequency quasi-periodic oscillations observed. In systems comprising a massive torus orbiting a black hole, on the other hand, p-mode oscillations could be used to explain the development or the suppression of the runaway instability.

Observations of X-ray emissions from binary systems have long since been considered important tools to test General Relativity in strong-field regimes. The high frequency quasi-periodic oscillations (HFQPOs) observed in binaries containing a black hole candidate, in particular, have been proposed as a means to measure more directly the black hole properties such as its mass and spin. Numerous models have been suggested to explain the HFQPOs and the rich phenomenology accompanying them. Many of these models rest on a number of assumptions and are at times in conflict with the most recent observations. We here propose a new, simple model in which the HFQPOs result from basic $p$-mode oscillations of a small accretion torus orbiting close to the black hole. We show that within this model the key properties of the HFQPOs can be explained simply, given a single reasonable assumption. We also discuss observational tests that can falsify the model.

We present general relativistic hydrodynamics simulations of constant specific angular momentum tori orbiting a Schwarzschild black hole. These tori are expected to form as a result of stellar gravitational collapse, binary neutron star merger or disruption, can reach very high rest-mass densities and behave effectively as neutron stars but with a toroidal topology (i.e. ``toroidal neutron stars''). Our attention is here focussed on the dynamical response of these objects to axisymmetric perturbations. We show that, upon the introduction of perturbations, these systems either become unstable to the runaway instability or exhibit a regular oscillatory behaviour resulting in a quasi-periodic variation of the accretion rate as well as of the mass quadrupole. The latter, in particular, is responsible for the emission of intense gravitational radiation whose signal-to-noise ratio at the detector is comparable or larger than the typical one expected in stellar-core collapse, making these new sources of gravitational waves potentially detectable. We discuss a systematic investigation of the parameter space both in the linear and nonlinear regimes, providing estimates of how the gravitational radiation emitted depends on the mass of the torus and on the strength of the perturbation.

We extend our approach for the exact solution of the Riemann problem in relativistic hydrodynamics to the case in which the fluid velocity has components tangential to the initial discontinuity. As in one-dimensional flows, we here show that the wave-pattern produced in a multidimensional relativistic Riemann problem can be predicted entirely by examining the initial conditions. Our method is logically very simple and allows for a numerical implementation of an exact Riemann solver which is both straightforward and computationally efficient. The simplicity of the approach is also important for revealing special relativistic effects responsible for a smooth transition from one wave-pattern to another when the tangential velocities in the initial states are suitably varied. While the content of this paper is focussed on a flat spacetime, the local Lorentz invariance allows its use also in fully general relativistic calculations.

In Newtonian and relativistic hydrodynamics the Riemann problem consists of calculating the evolution of a fluid which is initially characterized by two states having different values of uniform rest-mass density, pressure and velocity. When the fluid is allowed to relax, one of three possible wave-patterns is produced, corresponding to the propagation in opposite directions of two nonlinear hydrodynamical waves. New effects emerge in a special relativistic Riemann problem when velocities tangential to the initial discontinuity surface are present. We show that a smooth transition from one wave-pattern to another can be produced by varying the initial tangential velocities while otherwise maintaining the initial states unmodified. These special relativistic effects are produced by the coupling through the relativistic Lorentz factors and do not have a Newtonian counterpart.

We have solved numerically the general relativistic induction equations in the interior background spacetime of a slowly rotating magnetized neutron star. The analytic form of these equations was discussed in a recent paper (Rezzolla et al 2001a), where corrections due both to the spacetime curvature and to the dragging of reference frames were shown to be present. Through a number of calculations we have investigated the evolution of the magnetic field with different rates of stellar rotation, different inclination angles between the magnetic moment and the rotation axis, as well as different values of the electrical conductivity. All of these calculations have been performed for a constant temperature relativistic polytropic star and make use of a consistent solution of the initial value problem which avoids the use of artificial analytic functions. Our results show that there exist general relativistic effects introduced by the rotation of the spacetime which tend to decrease the decay rate of the magnetic field. The rotation-induced corrections are however generally hidden by the high electrical conductivity of the neutron star matter and when realistic values for the electrical conductivity are considered, these corrections become negligible even for the fastest known pulsar.

A Riemann problem with prescribed initial conditions will produce one of three possible wave patterns corresponding to the propagation of the different discontinuities that will be produced once the system is allowed to relax. In general, when solving the Riemann problem numerically, the determination of the specific wave pattern produced is obtained through some initial guess which can be successively discarded or improved. We here discuss a new procedure, suitable for implementation in an exact Riemann solver in one dimension, which removes the initial ambiguity in the wave pattern. In particular we focus our attention on the relativistic velocity jump between the two initial states and use this to determine, through some analytic conditions, the wave pattern produced by the decay of the initial discontinuity. The exact Riemann problem is then solved by means of calculating the root of a nonlinear equation. Interestingly, in the case of two rarefaction waves, this root can even be found analytically. Our procedure is straightforward to implement numerically and improves the efficiency of numerical codes based on exact Riemann solvers.

We study the processes of particle acceleration which take place in the field of a rotating black hole as part of a mechanism of formation of galactic jets within the first parsec from the central source, where gravitation is supposed to be dominant. We find the Lorentz factor that a stream of particles acquires as function of distance, when the orbital parameters vary slightly due to a local electromagnetic field or a pressure gradient.