We present a new method for the numerical solution of the radiative-transfer equation (RTE) in multidimensional scenarios commonly encountered in computational astrophysics. The method is based on the direct solution of the Boltzmann equation via an extension of the Lattice Boltzmann (LB) methods and allows to model the evolution of the radiation field as it interacts with a background fluid, via absorption, emission, and scattering. As a first application of this method, we restrict our attention to a frequency independent (``grey'') formulation within a special-relativistic framework, which can be employed also for classical computational astrophysics. For a number of standard tests that consider the performance of the method in optically thin, optically thick and intermediate regimes with a static fluid, we show the ability of the LB method to produce accurate and convergent results matching the analytic solutions. We also contrast the LB method with commonly employed moment-based schemes for the solution of the RTE, such as the M1 scheme. In this way, we are able to highlight that the LB method provides the correct solution for both non-trivial free-streaming scenarios and the intermediate optical-depth regime, for which the M1 method either fails or provides inaccurate solutions. When coupling to a dynamical fluid, on the other hand, we present the first self-consistent solution of the RTE with LB methods within a relativistic-hydrodynamic scenario. Finally, we show that besides providing more accurate results in all regimes, the LB method features smaller or comparable computational costs compared to the M1 scheme. We conclude that LB methods represent a competitive and promising avenue to the solution of radiative transport, one of the most common and yet important problems in computational astrophysics.

The recent detection of GW190814 featured the merger of a binary with a primary having a mass of $\sim 23\,M_{\odot}$ and a secondary with a mass of $\sim 2.6\,M_{\odot}$. While the primary was most likely a black hole, the secondary could be interpreted as either the lightest black hole or the most massive neutron star ever observed, but also as the indication of a novel class of exotic compact objects. We here argue that the secondary in GW190814 needs not be an ab-initio black hole nor an exotic object; rather, based on our current understanding of the nuclear-matter equation of state, it can be a rapidly rotating neutron star that collapsed to a rotating black hole at some point before merger. Using universal relations connecting the masses and spins of uniformly rotating neutron stars, we estimate the spin, $0.49 \lesssim \chi \lesssim 0.68$, of the secondary -- a quantity not constrained so far by the detection -- and a novel strict lower bound on the maximum mass, $M_{\rm TOV} > 2.08^{+0.04}_{-0.04}\, \,M_{\odot}$, of nonrotating neutron stars, consistent with recent observations of a very massive pulsar. The new lower bound also remains valid even in the less likely scenario in which the secondary neutron star never collapsed to a black hole.

We provide a systematic description of the steps necessary -- and of the potential pitfalls to be encountered -- when implementing a two-moment scheme within an Implicit-Explicit (IMEX) scheme to include radiative-transfer contributions in numerical simulations of general-relativistic (magneto-)hydrodynamics. We make use of the M1 closure, which provides an exact solution for the optically thin and thick limit, and an interpolation between these limits. Special attention is paid to the efficient solution of the emerging set of implicit conservation equations. In particular, we present an efficient method for solving these equations via the inversion of a $4\times 4$-matrix within an IMEX scheme. While this method relies on a few approximations, it offers a very good compromise between accuracy and computational efficiency. After a large number of tests in special relativity, we couple our new radiation code, \texttt{FRAC}, with the general-relativistic magnetohydrodynamics code \texttt{BHAC} to investigate the radiative Michel solution, namely, the problem of spherical accretion onto a black hole in the presence of a radiative field. By performing the most extensive exploration of the parameter space for this problem, we find that the accretion's efficiency can be expressed in terms of physical quantities such as temperature, $T$, luminosity, $L$, and black-hole mass, $M$, via the expression $\varepsilon=(L/L_{\rm Edd})/(\dot{M}/\dot{M}_{\rm Edd})= 7.41\times 10^{-7}\left(T/10^6\,\mathrm{K}\right)^{0.22} \left(L/L_\odot\right)^{0.48} \left(M/M_\odot\right)^{0.48}$, where $L_{\mathrm{Edd}}$ and $\dot{M}_{\mathrm{Edd}}$ are the Eddington luminosity and accretion rate, respectively. Finally, we also consider the accretion problem away from spherical symmetry, finding that the solution is stable under perturbations in the radiation field.

The first binary neutron-star merger event, GW170817, and its bright electromagnetic counterpart have provided a remarkable amount of information. By contrast, the second event, GW190425, with $M_{\rm tot}=3.4^{+0.3}_{-0.1}\,M_{\odot}$ and the lack of an electromagnetic counterpart, has hardly improved our understanding of neutron-star physics. While GW190425 is compatible with a scenario in which the merger has lead to a prompt collapse to a black hole and little ejected matter to power a counterpart, determining the mass ratio and the effective spin $\tilde{\chi}$ of the binary remains difficult. This is because gravitational waveforms cannot yet well constrain the component spins of the binary. However, since the mass of GW190425 is significantly larger than the maximum mass for nonrotating neutron stars, $M_{_{\rm TOV}}$, the mass ratio $q$ cannot be too small, as the heavier star would not be gravitationally stable. Making use of universal relations and a large number of equations of state, we provide limits in the $(\tilde{\chi},q)$ plane for GW190425, namely: $q_{\rm min} \geq 0.38$ and $ \tilde{\chi}_{\rm max}\leq 0.20$, assuming $M_\mathrm{tot} \simeq 3.4\, M_\odot$. Finally, we show how future observations of high-mass binaries can provide a lower bound on $M_{_{\rm TOV}}$.

The evolution of the gravitational-wave phase in the signal produced by inspiralling binaries of compact stars is modified by the nonzero deformability of the two stars. Hence, the measurement of these corrections has the potential of providing important information on the equation of state of nuclear matter. Extensive work has been carried out over the last decade to quantify these corrections, but it has so far been restricted to stars with zero intrinsic magnetic fields. While the corrections introduced by the magnetic tension and magnetic pressure are expected to be subdominant, it is nevertheless useful to determine the precise conditions under which these corrections become important. To address this question, we have carried out a second-order perturbative analysis of the tidal deformability of magnetised compact stars under a variety of magnetic-field strengths and equations of state describing either neutron stars or quark stars. Overall, we find that magnetically induced corrections to the tidal deformability will produce changes in the gravitational-wave phase evolution that are unlikely to be detected for realistic magnetic field i.e., $B\sim 10^{10} - 10^{12}\,{\rm G}$. At the same time, if the magnetic field is unrealistically large, i.e., $B\sim 10^{16}\,{\rm G}$, these corrections would produce a sizeable contribution to the phase evolution, especially for quark stars. In the latter case, the induced phase differences would represent a unique tool to measure the properties of the magnetic fields, providing information that is otherwise hard to quantify.

Very Long Baseline Interferometry (VLBI) at sub-millimeter waves has the potential to image the shadow of the black hole in the Galactic Center, Sagittarius A* (Sgr A*), and thereby test basic predictions of the theory of general relativity. We investigate the imaging prospects of a new Space VLBI mission concept. The setup consists of two satellites in polar or equatorial circular Medium-Earth Orbits with slightly different radii, resulting in a dense spiral-shaped uv-coverage with long baselines, allowing for extremely high-resolution and high-fidelity imaging of radio sources. We simulate observations of a general relativistic magnetohydrodynamics model of Sgr A* for this configuration with noise calculated from model system parameters. After gridding the $uv$-plane and averaging visibilities accumulated over multiple months of integration, images of Sgr A* with a resolution of up to 4 $\mu$as could be reconstructed, allowing for stronger tests of general relativity and accretion models than with ground-based VLBI.

With the first detection of gravitational waves from a binary system of neutron stars, GW170817, a new window was opened to study the properties of matter at and above nuclear-saturation density. Reaching densities a few times that of nuclear matter and temperatures up to $100\,\rm{MeV}$, such mergers also represent potential sites for a phase transition (PT) from confined hadronic matter to deconfined quark matter. While the lack of a postmerger signal in GW170817 has prevented us from assessing experimentally this scenario, two theoretical studies have explored the postmerger gravitational-wave signatures of PTs in mergers of binary systems of neutron stars. We here extend and complete the picture by presenting a novel signature of the occurrence of a PT. More specifically, using fully general-relativistic hydrodynamic simulations and employing a suitably constructed equation of state that includes a PT, we present the occurrence of a "delayed PT", i.e. a PT that develops only some time after the merger and produces a metastable object with a quark-matter core, i.e. a hypermassive hybrid star. Because in this scenario, the postmerger signal exhibits two distinct fundamental gravitational-wave frequencies -- before and after the PT -- the associated signature promises to be the strongest and cleanest among those considered so far, and one of the best signatures of the production of quark matter in the present Universe.

With the goal of providing more accurate and realistic estimates of the secular behavior of the mass accretion and drag rates in the "common-envelope" scenario encountered when a black hole or a neutron star moves in the stellar envelope of a red supergiant star, we have carried out the first general relativistic simulations of the accretion flow onto a nonrotating black hole moving supersonically in a medium with regular but different density gradients. The simulations reveal that the supersonic motion always rapidly reaches a stationary state and it produces a shock cone in the downstream part of the flow. In the absence of density gradients we recover the phenomenology already observed in the well-known Bondi--Hoyle--Lyttleton accretion problem, with super-Eddington mass accretion rate and a shock cone whose axis is stably aligned with the direction of motion. However, as the density gradient is made stronger, the accretion rate also increases and the shock cone is progressively and stably dragged toward the direction of motion. With sufficiently large gradients, the shock-cone axis can become orthogonal to the direction, or even move in the upstream region of the flow in the case of the largest density gradient. Together with the phenomenological aspects of the accretion flow, we have also quantified the rates of accretion of mass and momentum onto the black hole. Simple analytic expressions have been found for the rates of accretion of mass, momentum, drag force, and bremsstrahlung luminosity, all of which have been employed in the astrophysical modelling of the secular evolution of a binary system experiencing a common-envelope evolution...

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We report on the first detection of pulsed radio emission from a radio pulsar with the ALMA telescope. The detection was made in the Band-3 frequency range (85-101 GHz) using ALMA in the phased-array mode developed for VLBI observations. A software pipeline has been implemented to enable a regular pulsar observing mode in the future. We describe the pipeline and demonstrate the capability of ALMA to perform pulsar timing and searching. We also measure the flux density and polarization properties of the Vela pulsar (PSR J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons learned from the Vela observations for future pulsar studies with ALMA, particularly for searches near the supermassive black hole in the Galactic Center, and the potential of using pulsars for polarization calibration of ALMA.

We investigate the impact of using high-order numerical methods to study the merger of magnetised neutron stars with finite-temperature microphysics and neutrino cooling in full general relativity. By implementing a fourth-order accurate conservative finite-difference scheme we model the inspiral together with the early post-merger and highlight the differences to traditional second-order approaches at the various stages of the simulation. We find that even for finite-temperature equations of state, convergence orders higher than second order can be achieved in the inspiral and post-merger for the gravitational-wave phase. We further demonstrate that the second-order scheme overestimates the amount of proton-rich shock-heated ejecta, which can have an impact on the modelling of the dynamical part of the kilonova emission. Finally, we show that already at low resolution the growth rate of the magnetic energy is consistently resolved by using a fourth-order scheme.

The observations of GW170817/GRB170817A have confirmed that the coalescence of a neutron-star binary is the progenitor of a short gamma-ray burst. In the standard picture of a short gamma-ray burst, a collimated highly relativistic outflow is launched after merger and it successfully breaks out from the surrounding ejected matter. Using initial conditions inspired from numerical-relativity binary neutron-star merger simulations, we have performed general-relativistic hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations in which the jet is launched and propagates self-consistently. The complete set of simulations suggests that: (i) MHD jets have an intrinsic energy and velocity polar structure with a ``hollow core'' subtending an angle $\theta_{\rm core}\approx4^{\circ}-5^{\circ}$ and an opening angle of $\theta_{\rm jet}\gtrsim10^{\circ}$; (ii) MHD jets eject significant amounts of matter and two orders of magnitude more than HD jets; (iii) the energy stratification in MHD jets naturally yields the power-law energy scaling $E(>\Gamma\beta)\propto(\Gamma\beta)^{-4.5}$; (iv) MHD jets provide fits to the afterglow data from GRB170817A that are comparatively better than those of the HD jets and without free parameters; (v) finally, both of the best-fit HD/MHD models suggest an observation angle $\theta_{\rm obs} \simeq 21^{\circ}$ for GRB170817A.

One of the main dissipation processes acting on all scales in relativistic jets is thought to be governed by magnetic reconnection. Such dissipation processes have been studied in idealized environments, such as reconnection layers, which evolve in merging islands and lead to the production of plasmoids, ultimately resulting in efficient particle acceleration. In accretion flows onto black holes, reconnection layers can be developed and destroyed rapidly during the turbulent evolution of the flow. We present a series of two-dimensional general-relativistic magnetohydrodynamic simulations of tori accreting onto rotating black holes focusing our attention on the formation and evolution of current sheets. Initially, the tori are endowed with a poloidal magnetic field having a multi-loop structure along the radial direction and with an alternating polarity. During reconnection processes, plasmoids and plasmoid chains are developed leading to a flaring activity and hence to a variable electromagnetic luminosity. We describe the methods developed to track automatically the plasmoids that are generated and ejected during the simulation, contrasting the behaviour of multi-loop initial data with that encountered in typical simulations of accreting black holes having initial dipolar field composed of one loop only. Finally, we discuss the implications that our results have on the variability to be expected in accreting supermassive black holes.

High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.

We study in detail the nuclear aspects of a neutron-star merger in which deconfinement to quark matter takes place. For this purpose, we make use of the Chiral Mean Field (CMF) model, an effective relativistic model that includes self-consistent chiral symmetry restoration and deconfinement to quark matter and, for this reason, predicts the existence of different degrees of freedom depending on the local density/chemical potential and temperature. We then use the out-of-chemical-equilibrium finite-temperature CMF equation of state in full general-relativistic simulations to analyze which regions of different QCD phase diagrams are probed and which conditions, such as strangeness and entropy, are generated when a strong first-order phase transition appears. We also investigate the amount of electrons present in different stages of the merger and discuss how far from chemical equilibrium they can be and, finally, draw some comparisons with matter created in supernova explosions and heavy-ion collisions.

In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime.

Following the detection of GW170817 and the accompanying kilonova AT2017gfo, it has become crucial to model and understand the various channels through which mass is ejected in neutron-star binary mergers. We discuss the impact that high stellar spins prior to merger have on the ejection of mass focussing, in particular, on the dynamically ejected mass by performing general-relativistic magnetohydrodynamic simulations employing finite-temperature equations of state and neutrino-cooling effects. Using eight different models with dimensionless spins ranging from $\chi\simeq-0.14$ to $\chi\simeq0.29$ we discuss how the presence of different spins affects the angular distribution and composition of the ejected matter. Most importantly, we find that the dynamical component of the ejected mass can be strongly suppressed in the case of high spins aligned with the orbital angular momentum. In this case, in fact, the merger remnant has an excess angular momentum yielding a more extended and "colder" object, with reduced ability to shed mass dynamically. We discuss how this result impacts the analysis of the recent merger event GW170817 and its kilonova afterglow.

Binary neutron-star mergers have long been associated with short-duration gamma-ray bursts (GRBs). This connection was confirmed with the first coincident detection of gravitational waves together with electromagnetic radiation from GW170817. The basic paradigm for short-duration GRBs includes an ultra-relativistic jet, but the low-luminosity prompt emission together with follow-up radio and X-ray observations have hinted that this picture may be different in the case of GW170817. In particular, it has been proposed that large amounts of the magnetic energy that is amplified after the merger, can be released when the remnant collapses to a black hole, giving rise to a quasi-spherical explosion impacting on the merger ejecta. Through numerical simulations we investigate this scenario for a range of viewing angles, injected energies and matter densities at the time of the collapse. Depending on the magnitude of the energy injection and the remnant density, we find two types of outflows: one with a narrow relativistic core and one with a wide-angle, but mildly relativistic outflow. Furthermore, very wide outflows are possible, but require energy releases in excess of 10^52 erg.

The astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

Exploiting a very large library of physically plausible equations of state (EOSs) containing more than $10^{7}$ members and yielding more than $10^{9}$ stellar models, we conduct a survey of the impact that a neutron-star radius measurement via electromagnetic observations can have on the EOS of nuclear matter. Such measurements are soon to be expected from the ongoing \textit{Neutron Star Interior Composition Explorer} mission and will complement the constraints on the EOS from gravitational-wave detections. Thanks to the large statistical range of our EOS library, we can obtain a first quantitative estimate of the commonly made assumption that the high-density part of the EOS is best constrained when measuring the radius of the most massive, albeit rare, neutron stars with masses $M\gtrsim2.1\,M_\odot$. At the same time, we find that radius measurements of neutron stars with masses $M\simeq1.7-1.85\,M_\odot$ can provide the strongest constraints on the low-density part of the EOS. Finally, we quantify how radius measurements by future missions can further improve our understanding of the EOS of matter at nuclear densities.

Recent advances in black hole astrophysics, particularly the first visual evidence of a supermassive black hole at the center of the galaxy M87 by the Event Horizon Telescope (EHT), and the detection of an orbiting "hot spot" nearby the event horizon of Sgr A* in the Galactic center by the Gravity Collaboration, require the development of novel numerical methods to understand the underlying plasma microphysics. Non-thermal emission related to such hot spots is conjectured to originate from plasmoids that form due to magnetic reconnection in thin current layers in the innermost accretion zone. Resistivity plays a crucial role in current sheet formation, magnetic reconnection, and plasmoid growth in black hole accretion disks and jets. We included resistivity in the three-dimensional general-relativistic magnetohydrodynamics (GRMHD) code BHAC and present the implementation of an Implicit-Explicit scheme to treat the stiff resistive source terms of the GRMHD equations. The algorithm is tested in combination with adaptive mesh refinement to resolve the resistive scales and a constrained transport method to keep the magnetic field solenoidal. Several novel methods for primitive variable recovery, a key part in relativistic magnetohydrodynamics codes, are presented and compared for accuracy, robustness, and efficiency. We propose a new inversion strategy that allows for resistive-GRMHD simulations of low gas-to-magnetic pressure ratio and highly magnetized regimes as applicable for black hole accretion disks, jets, and neutron star magnetospheres. We apply the new scheme to study the effect of resistivity on accreting black holes, accounting for dissipative effects as reconnection.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

Worldwide very long baseline radio interferometry arrays are expected to obtain horizon-scale images of supermassive black hole candidates as well as of relativistic jets in several nearby active galactic nuclei. This motivates the development of models for magnetohydrodynamic flows in strong gravitational fields. The Black Hole Accretion Code (BHAC) intends to aid with the modelling of such sources by means of general relativistic magnetohydrodynamical (GRMHD) simulations in arbitrary stationary spacetimes. New additions were required to guarantee an accurate evolution of the magnetic field when small and large scales are captured simultaneously. We discuss the adaptive mesh refinement (AMR) techniques employed in BHAC, essential to keep several problems computationally tractable, as well as staggered-mesh-based constrained transport (CT) algorithms to preserve the divergence-free constraint of the magnetic field, including a general class of prolongation operators for face-allocated variables compatible with them. Through several standard tests, we show that the choice of divergence-control method can produce qualitative differences in simulations of scientifically relevant accretion problems. We demonstrate the ability of AMR to reduce the computational costs of accretion simulations while sufficiently resolving turbulence from the magnetorotational instability. In particular, we describe a simulation of an accreting Kerr black hole in Cartesian coordinates using AMR to follow the propagation of a relativistic jet while self-consistently including the jet engine, a problem set up-for which the new AMR implementation is particularly advantageous. The CT methods and AMR strategies discussed here are being employed in the simulations performed with BHAC used in the generation of theoretical models for the Event Horizon Telescope Collaboration.

The galaxy M 87 harbors a kiloparsec-scale relativistic jet, whose origin coincides with a supermassive black hole. Observational mm-VLBI campaigns are capable of resolving the jet-launching region at the scale of the event horizon. In order to provide a context for interpreting these observations, realistic general-relativistic magnetohydrodynamical (GRMHD) models of the accretion flow are constructed. The characteristics of the observed spectral-energy distribution (SED) depend on the shape of the electrons' energy-distribution function (eDF). The dependency on the eDF is omitted in the modeling of the first Event Horizon Telescope results. In this work, we aim to model the M 87 SED from radio up to NIR/optical frequencies using a thermal-relativistic Maxwell- J\"uttner distribution, as well as a relativistic $\kappa$-distribution function. The electrons are injected based on sub-grid, particle-in-cell parametrizations for sub-relativistic reconnection. A GRMHD simulation in Cartesian-Kerr-Schild coordinates, using eight levels of adaptive mesh refinement (AMR), forms the basis of our model. To obtain spectra and images, the GRMHD data is post-processed with the ray-tracing code RAPTOR, which is capable of ray tracing through AMR GRMHD simulation data. We obtain radio spectra in both the thermal-jet and $\kappa$-jet models consistent with radio observations. Additionally, the $\kappa$-jet models also recover the NIR/optical emission. The models recover the observed source sizes and core shifts and obtain a jet power of $\approx 10^{43}$ ergs/s. In the $\kappa$-jet models, both the accretion rates and jet powers are approximately two times lower than the thermal-jet model. The frequency cut-off observed at $\nu \approx 10^{15}$ Hz is recovered when the accelerator size is $10^6$ - $10^8$ cm, this could potentially point to an upper limit for plasmoid sizes in the jet of M 87.

If a phase transition is allowed to take place in the core of a compact star, a new stable branch of equilibrium configurations can appear, providing solutions with the same mass as the purely hadronic branch and hence giving rise to twin-star configurations. We perform an extensive analysis of the features of the phase transition leading twin-star configurations and, at the same time, fulfilling the constraints coming from the maximum mass of $2M_\odot$ and the information following gravitational-wave event GW170817. In particular, we use a general equation of state for the neutron-star matter that parametrizes the hadron-quark phase transition between the model describing the hadronic phase and a constant speed of sound for the quark phase. We find that the largest number of twin-star solutions has masses in the neutron-star branch in the range $1-2M_\odot$ and twin-branch masses $\gtrsim 2M_\odot$. The analysis of the masses, radii and tidal deformabilities also reveals that when twin stars appear, the tidal deformability shows two distinct branches with the same mass, thus differing considerably from the behaviour expected for neutron stars. In addition, we find that the data from GW170817 is compatible with the existence of hybrid stars and could also be interpreted as produced by the merger of a binary system of hybrid stars or of a hybrid star with a neutron star. The presence of a hybrid star in the inspiral phase can be established clearly if future gravitational-wave detections measure chirp masses $\mathcal{M}\lesssim 1.2M_\odot$ and tidal deformabilities of $\Lambda_{1.4}\lesssim 400$ for $1.4M_\odot$ stars. Finally, combining all observational information available, we set constraints on the parameters that characterise the phase transition, the maximum masses, and the radii of $1.4M_\odot$ stars described by equations of state leading to twin-star configurations.

In this work, we discuss the dense matter equation of state (EOS) for the extreme range of conditions encountered in neutron stars and their mergers. The calculation of the properties of such an EOS involves modeling different degrees of freedom (such as nuclei, nucleons, hyperons, and quarks), taking into account different symmetries, and including finite density and temperature effects in a thermodynamically consistent manner. We begin by addressing subnuclear matter consisting of nucleons and a small admixture of light nuclei in the context of the excluded volume approach. We then turn our attention to supranuclear homogeneous matter as described by the Chiral Mean Field (CMF) formalism. Finally, we present results from realistic neutron-star-merger simulations performed using the CMF model that predict signatures for deconfinement to quark matter in gravitational wave signals.

The main hard pulse of prompt gamma-ray emission in GRB$\,$170817A had a duration of $\sim0.5\,{\rm s}$ and its onset was delayed with respect to the gravitational-wave chirp signal by $t_{\rm del} \approx 1.74\,{\rm s}$. Detailed follow-up of the subsequent broadband kilonova emission revealed a two-component ejecta -- a lanthanide-poor ejecta with mass $M_{\rm ej,blue}\approx0.025\,M_\odot$ that powered the early but rapidly fading blue emission and a lanthanide-rich ejecta with mass $M_{\rm ej,red}\approx 0.04\,M_\odot$ that powered the longer lasting redder emission. Both the prompt gamma-ray onset delay and the existence of the blue ejecta with modest electron fraction, $0.2\lesssim Y_e\lesssim0.3$, can be explained if the collapse to a black hole was delayed by the formation of a hypermassive neutron star (HMNS). Here, we determine the survival time of the merger remnant by combining two different constraints, namely, the time needed to produce the requisite blue-ejecta mass and that necessary for the relativistic jet to bore its way out of the expanding ejecta. In this way, we determine that the remnant of GW170817 must have collapsed to a black hole after $t_{\rm coll}=0.98_{-0.26}^{+0.31}\,{\rm s}$. We also discuss how future detections and the delays between the gravitational and electromagnetic emissions can be used to constrain the properties of the merged object.

It has been proposed that Very Long Baseline Interferometry (VLBI) at sub-millimeter waves will allow us to image the shadow of the black hole in the center of our Milky Way, Sagittarius A* (Sgr A*), and thereby test basic predictions of general relativity. This paper presents imaging simulations of a new Space VLBI mission concept. An initial design study of the concept has been presented as the Event Horizon Imager (EHI). The EHI may be suitable for imaging Sgr A* at high frequencies (up to ~690 GHz), which has significant advantages over performing ground-based VLBI at 230 GHz. The concept EHI design consists of two or three satellites in polar or equatorial circular Medium-Earth Orbits with slightly different radii. Due to the relative drift of the satellites along the individual orbits, this setup will result in a dense spiral-shaped uv-coverage with long baselines (up to ~60 Glambda), allowing for extremely high-resolution and high-fidelity imaging of radio sources. We simulate EHI observations of general relativistic magnetohydrodynamics models of Sgr A* and calculate the expected noise based on preliminary system parameters. On long baselines, where the signal-to-noise ratio may be low, fringes could be detected if the system is sufficiently phase stable and the satellite orbits can be reconstructed with sufficient accuracy. Averaging visibilities accumulated over multiple epochs of observations could then help improving the image quality. With three satellites, closure phases could be used for imaging. Our simulations show that the EHI could be capable of imaging the black hole shadow of Sgr A* with a resolution of 4 uas (about 8% of the shadow diameter) within several months of observing time. The EHI concept could thus be used to measure black hole shadows much more precisely than with ground-based VLBI, allowing for stronger tests of general relativity and accretion models.

Radio-astronomical observations of the supermassive black-hole candidate in the galactic center will soon offer the possibility to study gravity in its strongest regimes and to test different models for these compact objects. Studies based on semi-analytic models and strong-field images of stationary plasma configurations around boson stars have stressed the difficulty to distinguish them from black holes. For the first time, we present general-relativistic magnetohydrodynamic simulations followed by consistent general-relativistic radiative-transfer calculations in a boson-star spacetime to revisit the appearance of these objects in a realistic accretion scenario. We find that the absence of an event horizon in a boson star leads to important differences in the dynamics of the accretion and results in both the formation of a small torus in the interior of the boson star and in the absence of an evacuated high-magnetization funnel in the polar regions. Synthetic reconstructed images considering realistic astronomical observing conditions show that differences in the appearance of the two compact objects are large enough to be detectable. These results, which also apply to other horizonless compact objects, strengthen confidence in the ability to determine the presence of an event horizon via radio observations and highlight the importance of self-consistent multidimensional simulations to study the compact object at the galactic center.

High-resolution Very-Long-Baseline Interferometry observations of NGC 1052 show a two sided jet with several regions of enhanced emission and a clear emission gap between the two jets.This gap shrinks with increasing frequency and vanishes around $\nu\sim43$ GHz. The observed structures are due to both the macroscopic fluid dynamics interacting with the surrounding ambient medium including an obscuring torus and the radiation microphysics. In this paper we investigate the possible physical conditions in relativistic jets of NGC 1052 by directly modelling the observed emission and spectra via state-of-the-art special-relativistic hydrodynamic (SRHD) simulations and radiative transfer calculations. To investigate the physical conditions in the relativistic jet we coupled our radiative transfer code to evolutionary algorithms and performed simultaneous modelling of the observed jet structure and the broadband radio spectrum. During the calculation of the radiation we consider both thermal and non-thermal emission. In order to compare our model to VLBI observations we take into account the sparse sampling of the u-v plane, the array properties and the imaging algorithm. We present for the first time an end-to-end pipeline for fitting numerical simulations to VLBI observations of relativistic jets taking into account the macrophysics including fluid dynamics and ambient medium configurations together with thermal/non-thermal emission and the properties of the observing array. The detailed analysis of our simulations shows that the structure and properties of the observed relativistic jets in NGC 1052 can be reconstructed by a slightly over-pressured jet ($d_k\sim1.5$) embedded in a decreasing pressure ambient medium

We study the lifetimes of the remnant produced by the merger of two neutron stars and revisit the determination of the threshold mass to prompt collapse, $M_{\rm th}$. Using a fully general-relativistic numerical approach and a novel method for a rigorous determination of $M_{\rm th}$, we show that a nonlinear universal relation exists between the threshold mass and the maximum compactness. For the temperature-dependent equations of state considered here, our results improve a similar linear relation found recently with methods that are less accurate but yield quantitatively similar results. Furthermore, exploiting the information from GW170817, we use the universal relation to set lower limits on the stellar radii for any mass.

Merging binaries of neutron stars are not only strong sources of gravitational waves, but also have the potential of revealing states of matter at densities and temperatures not accessible in laboratories. A crucial and long-standing question in this context is whether quarks are deconfined as a result of the dramatic increase in density and temperature following the merger. We present the first fully general-relativistic simulations of merging neutron stars including quarks at finite temperatures that can be switched off consistently in the equation of state. Within our approach, we can determine clearly what signatures a quark-hadron phase transition would leave in the gravitational-wave signal. In particular, we show that if the conditions are met for a phase transition to take place at several times nuclear saturation density, they would lead to a post-merger signal considerably different from the one expected from the inspiral, that can only probe the hadronic part of the equations of state, and to an anticipated collapse of the merged object. We also show that the phase transition leads to a very hot and dense quark core that, when it collapses to a black hole, produces a ringdown signal different from the hadronic one. Finally, in analogy with what is done in heavy-ion collisions, we use the evolution of the temperature and density in the merger remnant to illustrate the properties of the phase transition in a QCD phase diagram.

With the recent advent of multi-messenger gravitational-wave astronomy and in anticipation of more sensitive, next-generation gravitational-wave detectors, we investigate the dynamics, gravitational-wave emission, and nucleosynthetic yields of numerous eccentric binary neutron-star mergers having different equations of state. For each equation of state we vary the orbital properties around the threshold of immediate merger, as well as the binary mass ratio. In addition to a study of the gravitational-wave emission including $f$-mode oscillations before and after merger, we couple the dynamical ejecta output from the simulations to the nuclear-reaction network code \texttt{SkyNet} to compute nucleosynthetic yields and compare to the corresponding results in the case of a quasi-circular merger. We find that the amount and velocity of dynamically ejected material is always much larger than in the quasi-circular case, reaching maximal values of $M_{\rm ej, max} \sim 0.1 \, M_{\odot}$ and $v_{\rm max}/c \sim 0.75$. At the same time, the properties of this material are rather insensitive to the details of the orbit, such as pericenter distance or post-encounter apoastron distance. Furthermore, while the composition of the ejected matter depends on the orbital parameters and on the equation of state, the relative nucleosynthetic yields do not, thus indicating that kilonova signatures could provide information on the orbital properties of dynamically captured neutron-star binaries.

Gravitational waves have been detected from the inspiral of a binary neutron-star, GW170817, which allowed constraints to be placed on the neutron star equation of state. The equation of state can be further constrained if gravitational waves from a post-merger remnant are detected. Post-merger waveforms are currently generated by numerical-relativity simulations, which are computationally expensive. Here we introduce a hierarchical model trained on numerical-relativity simulations, which can generate reliable post-merger spectra in a fraction of a second. Our spectra have mean fitting factors of 0.95, which compares to fitting factors of 0.76 and 0.85 between different numerical-relativity codes that simulate the same physical system. This method is the first step towards generating large template banks of spectra for use in post-merger detection and parameter estimation.

We examine the range of rest-mass densities, temperatures and magnetic fields involved in simulations of binary neutron-star mergers and identify the conditions under which the ideal-magnetohydrodynamics approximation breaks down and hence the magnetic-field decay should be accounted for. We use recent calculations of the conductivities of warm correlated plasma in envelopes of compact stars and find that the magnetic-field decay timescales are much larger than the characteristic timescales of the merger process for lengthscales down to a meter. Because these are smaller than the currently available resolution in numerical simulations, the ideal-magnetohydrodynamics approximation is effectively valid for all realistic simulations. At the same time, we find that the Hall effect can be important at low densities and low temperatures, where it can induce a non-dissipative rearrangement of the magnetic field. Finally, we mark the region in temperature and density where the hydrodynamic description breaks down.

It has been suggested that a non-repeating fast radio burst (FRB) represents the final signal of a magnetized neutron star collapsing to a black hole. In this model, a supramassive neutron star supported by rapid rotation, will collapse to a black hole several thousand to million years after its birth as a result of spin down. The collapse violently snaps the magnetic-field lines anchored on the stellar surface, thus producing an electromagnetic pulse that will propagate outwards and accelerate electrons producing a massive radio burst, i.e. a "blitzar". We present a systematic study of the gravitational collapse of rotating and magnetised neutron stars with special attention to far-field evolution at late times after the collapse. By considering a series of neutron stars with rotation ranging from zero to millisecond periods and different magnetic-field strengths, we show that the blitzar emission is very robust and always characterised by a series sub-millisecond pulses decaying exponentially in amplitude. The luminosity and energy released when the magnetosphere is destroyed are well reproduced by a simple expression in terms of the stellar magnetic field and radius. Finally, we assess the occurrence of pair production during a blitzar scenario, concluding that for typical magnetic-field strengths of $10^{12}\,{\rm G}$ and spin frequencies of a few Hz, pair production is suppressed. Overall, the very good match between the results of the simulations and the luminosities normally observed for FRBs lends credibility to the blitzar model as a simple and yet plausible explanation for the phenomenology of non-repeating FRBs.

We explore in a parameterized manner a very large range of physically plausible equations of state (EOSs) for compact stars for matter that is either purely hadronic or that exhibits a phase transition. In particular, we produce two classes of EOSs with and without phase transitions, each containing one million EOSs. We then impose constraints on the maximum mass, ($M < 2.16 M_{\odot}$), and on the dimensionless tidal deformability ($\tilde{\Lambda} <800$) deduced from GW170817, together with recent suggestions of lower limits on $\tilde{\Lambda}$. Exploiting more than $10^9$ equilibrium models for each class of EOSs, we produce distribution functions of all the stellar properties and determine, among other quantities, the radius that is statistically most probable for any value of the stellar mass. In this way, we deduce that the radius of a purely hadronic neutron star with a representative mass of $1.4\,M_{\odot}$ is constrained to be $12.00\!<\!R_{1.4}/{\rm km}\!<\!13.45$ at a $2$-$\sigma$ confidence level, with a most likely value of $\bar{R}_{1.4}=12.39\,{\rm km}$; similarly, the smallest dimensionless tidal deformability is $\tilde{\Lambda}_{1.4}\!>\!375$, again at a $2$-$\sigma$ level. On the other hand, because EOSs with a phase transition allow for very compact stars on the so-called `twin-star' branch, small radii are possible with such EOSs although not probable, i.e. $8.53\!<\!R_{1.4}/{\rm km}\!<\!13.74$ and $\bar{R}_{1.4}=13.06\,{\rm km}$ at a $2$-$\sigma$ level, with $\tilde{\Lambda}_{1.4}\!>\!35.5$ at a $3$-$\sigma$ level. Finally, since these EOSs exhibit upper limits on $\tilde{\Lambda}$, the detection of a binary with total mass of $3.4\,M_{\odot}$ and $\tilde{\Lambda}_{1.7}\!>\!461$ can rule out twin-star solutions.

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.

Our Galactic Center, Sagittarius A* (Sgr A*), is believed to harbour a supermassive black hole (BH), as suggested by observations tracking individual orbiting stars. Upcoming sub-millimetre very-long-baseline-interferometry (VLBI) images of Sgr A* carried out by the Event-Horizon-Telescope Collaboration (EHTC) are expected to provide critical evidence for the existence of this supermassive BH. We assess our present ability to use EHTC images to determine if they correspond to a Kerr BH as predicted by Einstein's theory of general relativity (GR) or to a BH in alternative theories of gravity. To this end, we perform general-relativistic magnetohydrodynamical (GRMHD) simulations and use general-relativistic radiative transfer (GRRT) calculations to generate synthetic shadow images of a magnetised accretion flow onto a Kerr BH. In addition, and for the first time, we perform GRMHD simulations and GRRT calculations for a dilaton BH, which we take as a representative solution of an alternative theory of gravity. Adopting the VLBI configuration from the 2017 EHTC campaign, we find that it could be extremely difficult to distinguish between BHs from different theories of gravity, thus highlighting that great caution is needed when interpreting BH images as tests of GR.

The recent discovery of the electromagnetic counterpart of the gravitational wave source GW170817, has demonstrated the huge informative power of multi-messenger observations. During the next decade the nascent field of multi-messenger astronomy will mature significantly. Around 2030, third generation gravitational wave detectors will be roughly ten times more sensitive than the current ones. At the same time, neutrino detectors currently upgrading to multi km^3 telescopes, will include a 10 km^3 facility in the Southern hemisphere that is expected to be operational around 2030. In this review, we describe the most promising high frequency gravitational wave and neutrino sources that will be detected in the next two decades. In this context, we show the important role of the Transient High Energy Sky and Early Universe Surveyor (THESEUS), a mission concept proposed to ESA by a large international collaboration in response to the call for the Cosmic Vision Programme M5 missions. THESEUS aims at providing a substantial advancement in early Universe science as well as playing a fundamental role in multi-messenger and time-domain astrophysics, operating in strong synergy with future gravitational wave and neutrino detectors as well as major ground- and space-based telescopes. This review is an extension of the THESEUS white paper (Amati et al. 2017), also in light of the discovery of GW170817/GRB170817A that was announced on October 16th, 2017.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

At second order in perturbation theory, the unstable r-mode of a rotating star includes growing differential rotation whose form and growth rate are determined by gravitational-radiation reaction. With no magnetic field, the angular velocity of a fluid element grows exponentially until the mode reaches its nonlinear saturation amplitude and remains nonzero after saturation. With a background magnetic field, the differential rotation winds up and amplifies the field, and previous work where large mode amplitudes were considered suggests that the amplification may damp out the instability. A background magnetic field, however, turns the saturated time-independent perturbations corresponding to adding differential rotation into perturbations whose characteristic frequencies are of order the Alfv\'en frequency. As found in previous studies, we argue that magnetic- field growth is sharply limited by the saturation amplitude of an unstable mode. In contrast to previous work, however, we show that if the amplitude is small, i.e., of order 10^(-4), then the limit on the magnetic-field growth is stringent enough to prevent the loss of energy to the magnetic field from damping or significantly altering an unstable r-mode in nascent neutron stars with normal interiors and in cold stars whose interiors are type II superconductors. We show this result first for a toy model, and we then obtain an analogous upper limit on magnetic field growth using a more realistic model of a rotating neutron star. Our analysis depends on the assumption that there are no marginally unstable perturbations, and this may not hold when differential rotation leads to a magnetorotational instability.

Observational efforts to image the immediate environment of a black hole at the scale of the event horizon benefit from the development of efficient imaging codes that are capable of producing synthetic data, which may be compared with observational data. We aim to present RAPTOR, a new public code that produces accurate images, animations, and spectra of relativistic plasmas in strong gravity by numerically integrating the equations of motion of light rays and performing time-dependent radiative transfer calculations along the rays. The code is compatible with any analytical or numerical spacetime. It is hardware-agnostic and may be compiled and run both on GPUs and CPUs. We describe the algorithms used in RAPTOR and test the code's performance. We have performed a detailed comparison of RAPTOR output with that of other radiative-transfer codes and demonstrate convergence of the results. We then applied RAPTOR to study accretion models of supermassive black holes, performing time-dependent radiative transfer through general relativistic magneto-hydrodynamical (GRMHD) simulations and investigating the expected observational differences between the so-called fast-light and slow-light paradigms. Using RAPTOR to produce synthetic images and light curves of a GRMHD model of an accreting black hole, we find that the relative difference between fast-light and slow-light light curves is less than 5%. Using two distinct radiative-transfer codes to process the same data, we find integrated flux densities with a relative difference less than 0.01%. For two-dimensional GRMHD models, such as those examined in this paper, the fast-light approximation suffices as long as errors of a few percent are acceptable. The convergence of the results of two different codes demonstrates that they are, at a minimum, consistent.

Combining the GW observations of merging systems of binary neutron stars and quasi-universal relations, we set constraints on the maximum mass that can be attained by nonrotating stellar models of neutron stars. More specifically, exploiting the recent observation of the GW event GW 170817 and drawing from basic arguments on kilonova modeling of GRB 170817A, together with the quasi-universal relation between the maximum mass of nonrotating stellar models $M_{\rm TOV}$ and the maximum mass supported through uniform rotation $M_{\rm max}=\left(1.20^{+0.02}_{-0.05}\right) M_{\rm TOV}$ we set limits for the maximum mass to be $ 2.01^{+0.04}_{-0.04}\leq M_{\rm TOV}/M_{\odot}\lesssim 2.16^{+0.17}_{-0.15}$, where the lower limit in this range comes from pulsar observations. Our estimate, which follows a very simple line of arguments and does not rely on the modeling of the electromagnetic signal in terms of numerical simulations, can be further refined as new detections become available. We briefly discuss the impact that our conclusions have on the equation of state of nuclear matter.

We discuss magnetic field enhancement by unstable r-modes (driven by the gravitational radiation reaction force) in rotating stars. In the absence of a magnetic field, gravitational radiation exponentially increases the r-mode amplitude $\alpha$, and accelerates differential rotation (secular motion of fluid elements). For a magnetized star, differential rotation enhances the magnetic field energy. Rezzolla et al. (2000--2001) argued that if the magnetic energy grows faster than the gravitational radiation reaction force pumps energy into the r-modes, then the r-mode instability is suppressed. Chugunov (2015) demonstrated that without gravitational radiation, differential rotation can be treated as a degree of freedom decoupled from the r-modes and controlled by the back reaction of the magnetic field. In particular, the magnetic field windup does not damp r-modes. Here we discuss the effect of the back reaction of the magnetic field on differential rotation of unstable r-modes, and show that it limits the generated magnetic field and the magnetic energy growth rate preventing suppression of the r-mode instability by magnetic windup at low saturation amplitudes, $\alpha \ll 1$, predicted by current models.

Quasi-equilibrium models of uniformly rotating axisymmetric and triaxial quark stars are computed in general relativistic gravity scenario. The Isenberg-Wilson-Mathews (IWM) formulation is employed and the Compact Object CALculator (COCAL) code is extended to treat rotating stars with finite surface density and new equations of state (EOSs). Besides the MIT bag model for quark matter which is composed of de-confined quarks, we examine a new EOS proposed by Lai and Xu that is based on quark clustering and results in a stiff EOS that can support masses up to $3.3M_\odot$ in the case we considered. We perform convergence tests for our new code to evaluate the effect of finite surface density in the accuracy of our solutions and construct sequences of solutions for both small and high compactness. The onset of secular instability due to viscous dissipation is identified and possible implications are discussed. An estimate of the gravitational wave amplitude and luminosity based on quadrupole formulas is presented and comparison with neutron stars is discussed.

To date, the most precise tests of general relativity have been achieved through pulsar timing, albeit in the weak-field regime. Since pulsars are some of the most precise and stable "clocks" in the Universe, present observational efforts are focused on detecting pulsars in the vicinity of supermassive black holes (most notably in our Galactic Centre), enabling pulsar timing to be used as an extremely precise probe of strong-field gravity. In this paper a mathematical framework to describe test-particle dynamics in general black hole spacetimes is presented, and subsequently used to study a binary system comprising a pulsar orbiting a black hole. In particular, taking into account the parameterization of a general spherically symmetric black hole metric, general analytic expressions for both the advance of the periastron and for the orbital period of a massive test particle are derived. Furthermore, these expressions are applied to four representative cases of solutions arising in both general relativity and in alternative theories of gravity. Finally, this framework is applied to the Galactic Centre $S$-stars and four distinct pulsar toy models. It is shown that by adopting a fully general-relativistic description of test-particle motion which is independent of any particular theory of gravity, observations of pulsars can help impose better constraints on alternative theories of gravity than is presently possible.

The stability properties of rotating relativistic stars against prompt gravitational collapse to a black hole are rather well understood for uniformly rotating models. This is not the case for differentially rotating neutron stars, which are expected to be produced in catastrophic events such as the merger of binary system of neutron stars or the collapse of a massive stellar core. We consider sequences of differentially rotating equilibrium models using the $j$-constant law and by combining them with their dynamical evolution, we show that a sufficient stability criterion for differentially rotating neutron stars exists similar to the one of their uniformly rotating counterparts. Namely: along a sequence of constant angular momentum, a dynamical instability sets in for central rest-mass densities slightly below the one of the equilibrium solution at the turning point. In addition, following Breu & Rezzolla (2016), we show that "quasi-universal" relations can be found when calculating the turning-point mass. In turn, this allows us to compute the maximum mass allowed by differential rotation, $M_{\rm max,dr}$, in terms of the maximum mass of the nonrotating configuration, $M_{_{\rm TOV}}$, finding that $M_{\rm max, dr} \simeq \left(1.54 \pm 0.05\right) M_{_{\rm TOV}}$ for all the equations of state we have considered.

When binary systems of neutron stars merge, a very small fraction of their rest mass is ejected, either dynamically or secularly. This material is neutron-rich and its nucleosynthesis could provide the astrophysical site for the production of heavy elements in the universe, together with a kilonova signal confirming neutron-star mergers as the origin of short gamma-ray bursts. We perform full general-relativistic simulations of binary neutron-star mergers employing three different nuclear-physics EOSs, considering both equal- and unequal-mass configurations, and adopting a leakage scheme to account for neutrino radiative losses. Using a combination of techniques, we carry out an extensive and systematic study of the hydrodynamical, thermodynamical, and geometrical properties of the matter ejected dynamically, employing the WinNet nuclear-reaction network to recover the relative abundances of heavy elements produced by each configurations. Among the results obtained, three are particularly important. First, we find that both the properties of the dynamical ejecta and the nucleosynthesis yields are robust against variations of the EOS and masses, and match very well the observed chemical abundances. Second, using a conservative but robust criterion for unbound matter, we find that the amount of ejected mass is $\lesssim 10^{-3}\,M_{\odot}$, hence at least one order of magnitude smaller than what normally assumed in modelling kilonova signals. Finally, using a simplified and gray-opacity model we assess the observability of the infrared kilonova emission finding, that for all binaries the luminosity peaks around $\sim1/2$ day in the $H$-band, reaching a maximum magnitude of $-13$, and decreasing rapidly after one day. These rather low luminosities make the prospects for detecting kilonovae less promising than what assumed so far.

There have been many efforts to explain the dynamical mechanisms behind the phenomenology of quasi-periodic oscillations (QPOs) seen in the X-ray light curves of low-mass X-ray binaries. Up to now, none of the models can successfully explain all the frequencies observed in the power density spectrum of the light curve. Here we perform several hydrodynamic simulations of non-self-gravitating relativistic axisymmetric thick tori applied to the neutron star in the low-mass X-ray binary 4U 1636-53 and show how the observed oscillation modes triggered by different velocity perturbations give rise to a set of variability features similar to what we see in the observational X-ray data. When we match pairs of frequencies from the simulations of constant angular momentum ($l$) tori with the observed kilo Hertz QPOs, we find that certain combinations of frequencies lie on top of the observed relation, provided we assume a smaller mass for the neutron star than is generally assumed. However, constant-$l$ tori cannot match the entire range of frequencies observed for 4U 1636-53 due to physical constraints set by the torus size. We show that our model is consistent with the observed shift in QPO frequency that accompany state transitions of the accretion disk.

Inferring the properties of dense matter is one of the most exciting prospects from the measurement of gravitational waves from neutron star mergers. However, it will require reliable numerical simulations that incorporate viscous dissipation and energy transport if these can play a significant role within the survival time of the post-merger object. We calculate timescales for typical forms of dissipation and find that thermal transport and shear viscosity will not be important unless neutrino trapping occurs, which requires temperatures above about 10 MeV and gradients over lengthscales of 0.1 km or less. On the other hand, if direct-Urca processes remain suppressed, leaving modified-Urca processes to establish flavor equilibrium, then bulk viscous dissipation could provide significant damping to density oscillations observed right after the merger. When comparing with data from a state-of-the-art merger simulation we find that the bulk viscosity takes values close to its resonant maximum in a typical neutron-star merger, motivating a more careful assessment of the role of bulk viscous dissipation in the gravitational-wave signal from merging neutron stars.

Determining the differential-rotation law of compact stellar objects produced in binary neutron stars mergers or core-collapse supernovae is an old problem in relativistic astrophysics. Addressing this problem is important because it impacts directly on the maximum mass these objects can attain and hence on the threshold to black-hole formation under realistic conditions. Using the results from a large number of numerical simulations in full general relativity of binary neutron star mergers described with various equations of state and masses, we study the rotational properties of the resulting hypermassive neutron stars. We find that the angular-velocity distribution shows only a modest dependence on the equation of state, thus exhibiting the traits of "quasi-universality" found in other aspects of compact stars, both isolated and in binary systems. The distributions are characterized by an almost uniformly rotating core and a "disk". Such a configuration is significantly different from the $j-{\rm constant}$ differential-rotation law that is commonly adopted in equilibrium models of differentially rotating stars. Furthermore, the rest-mass contained in such a disk can be quite large, ranging from $\simeq 0.03\,M_{\odot}$ in the case of high-mass binaries with stiff equations of state, up to $\simeq 0.2\,M_{\odot}$ for low-mass binaries with soft equations of state. We comment on the astrophysical implications of our findings and on the long-term evolutionary scenarios that can be conjectured on the basis of our simulations.

We present the black hole accretion code (BHAC), a new multidimensional general-relativistic magnetohydrodynamics module for the MPI-AMRVAC framework. BHAC has been designed to solve the equations of ideal general-relativistic magnetohydrodynamics in arbitrary spacetimes and exploits adaptive mesh refinement techniques with an efficient block-based approach. Several spacetimes have already been implemented and tested. We demonstrate the validity of BHAC by means of various one-, two-, and three-dimensional test problems, as well as through a close comparison with the HARM3D code in the case of a torus accreting onto a black hole. The convergence of a turbulent accretion scenario is investigated with several diagnostics and we find accretion rates and horizon-penetrating fluxes to be convergent to within a few percent when the problem is run in three dimensions. Our analysis also involves the study of the corresponding thermal synchrotron emission, which is performed by means of a new general-relativistic radiative transfer code, BHOSS. The resulting synthetic intensity maps of accretion onto black holes are found to be convergent with increasing resolution and are anticipated to play a crucial role in the interpretation of horizon-scale images resulting from upcoming radio observations of the source at the Galactic Center.

We show how gravitational-wave observations with advanced detectors of tens to several tens of neutron-star binaries can measure the neutron-star radius with an accuracy of several to a few percent, for mass and spatial distributions that are realistic, and with none of the sources located within 100 Mpc. We achieve such an accuracy by combining measurements of the total mass from the inspiral phase with those of the compactness from the postmerger oscillation frequencies. For estimating the measurement errors of these frequencies we utilize analytical fits to postmerger numerical-relativity waveforms in the time domain, obtained here for the first time, for four nuclear-physics equations of state and a couple of values for the mass. We further exploit quasi-universal relations to derive errors in compactness from those frequencies. Measuring the average radius to well within 10% is possible for a sample of 100 binaries distributed uniformly in volume between 100 and 300 Mpc, so long as the equation of state is not too soft or the binaries are not too heavy.

High-resolution Very-Long-Baseline Interferometry observations of active galactic nuclei have revealed asymmetric structures in the jets of radio galaxies. These asymmetric structures may be due to internal asymmetries in the jet, could be induced by the different conditions in the surrounding ambient medium including the obscuring torus, or a combination of the two. In this paper we investigate the influence of the ambient medium (including the obscuring torus) on the observed properties of jets from radio galaxies. We performed special-relativistic hydrodynamic (RHD) simulations of over-pressured and pressure-matched jets using the special-relativistic hydrodynamics code \texttt{Ratpenat}, which is based on a second-order accurate finite-volume method and an approximate Riemann solver. Using a newly developed emission code to compute the electromagnetic emission, we have investigated the influence of different ambient medium and torus configurations on the jet structure and subsequently computed the non-thermal emission produced by the jet and the thermal absorption due to the torus. To better compare the emission simulations with observations we produced synthetic radio maps, taking into account the properties of the observatory. The detailed analysis of our simulations shows that the observed asymmetries can be produced by the interaction of the jet with the ambient medium and by the absorption properties of the obscuring torus.

The merger of binary neutron-stars systems combines in a single process: extreme gravity, copious emission of gravitational waves, complex microphysics, and electromagnetic processes that can lead to astrophysical signatures observable at the largest redshifts. We review here the recent progress in understanding what could be considered Einstein's richest laboratory, highlighting in particular the numerous significant advances of the last decade. Although special attention is paid to the status of models, techniques, and results for fully general-relativistic dynamical simulations, a review is also offered on initial data and advanced simulations with approximate treatments of gravity. Finally, we review the considerable amount of work carried out on the post-merger phase, including: black-hole formation, torus accretion onto the merged compact object, connection with gamma-ray burst engines, ejected material, and its nucleosynthesis.

We present the first systematic study of the gravitational collapse of rotating and magnetised neutron stars to charged and rotating (Kerr-Newman) black holes. In particular, we consider the collapse of magnetised and rotating neutron stars assuming that no pair-creation takes place and that the charge density in the magnetosphere is so low that the stellar exterior can be described as an electrovacuum. Under these assumptions, which are rather reasonable for a pulsar that has crossed the 'death line', we show that when the star is rotating, it acquires a net initial electrical charge, which is then trapped inside the apparent horizon of the newly formed back hole. We analyse a number of different quantities to validate that the black hole produced is indeed a Kerr-Newman one and show that, in the absence of rotation or magnetic field, the end result of the collapse is a Schwarzschild or Kerr black hole, respectively.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

(Abridged) We here continue our effort to model the behaviour of matter when orbiting or accreting onto a generic black hole by developing a new numerical code employing advanced techniques geared solve the equations of in general-relativistic hydrodynamics. The new code employs a number of high-resolution shock-capturing Riemann-solvers and reconstruction algorithms, exploiting the enhanced accuracy and the reduced computational cost of AMR techniques. In addition, the code makes use of sophisticated ray-tracing libraries that, coupled with general-relativistic radiation-transfer calculations, allow us to compute accurately the electromagnetic emissions from such accretion flows. We validate the new code by presenting an extensive series of stationary accretion flows either in spherical or axial symmetry and performed either in 2D or 3D. In addition, we consider the highly nonlinear scenario of a recoiling black hole produced in the merger of a supermassive black hole binary interacting with the surrounding circumbinary disc. In this way we can present, for the first time, ray-traced images of the shocked fluid and the light-curve resulting from consistent general-relativistic radiation-transport calculations from this process. The work presented here lays the ground for the development of a generic computational infrastructure employing AMR techniques to deal accurately and self-consistently with accretion flows onto compact objects. In addition to the accurate handling of the matter, we provide a self-consistent electromagnetic emission from these scenarios by solving the associated radiative-transfer problem. While magnetic fields are presently excluded from our analysis, the tools presented here can have a number of applications to study accretion flows onto black holes or neutron stars.

Collaborative international efforts under the name of the Event Horizon Telescope project, using sub- mm very long baseline interferometry, are soon expected to provide the first images of the shadow cast by the candidate supermassive black hole in our Galactic center, Sagittarius A*. Observations of this shadow would provide direct evidence of the existence of astrophysical black holes. Although it is expected that astrophysical black holes are described by the axisymmetric Kerr solution, there also exist many other black hole solutions, both in general relativity and in other theories of gravity, which cannot presently be ruled out. To this end, we present calculations of black hole shadow images from various metric theories of gravity as described by our recent work on a general parameterisation of axisymmetric black holes [R. Konoplya, L. Rezzolla and A. Zhidenko, Phys. Rev. D 93, 064015 (2016)]. An algorithm to perform general ray-tracing calculations for any metric theory of gravity is first outlined and then employed to demonstrate that even for extremal metric deformation parameters of various black hole spacetimes, this parameterisation is both robust and rapidly convergent to the correct solution.

The interferometric LIGO detectors have recently measured the first direct gravitational-wave signal from what has been interpreted as the inspiral, merger and ringdown of a binary system of black holes. The signal-to-noise ratio of the measured signal is large enough to leave little doubt that it does refer to the inspiral of two massive and ultracompact objects, whose merger yields a rotating black hole. Yet, the quality of the data is such that some room is left for alternative interpretations that do not involve black holes, but other objects that, within classical general relativity, can be equally massive and compact, namely, gravastars. We here consider the hypothesis that the merging objects were indeed gravastars and explore whether the merged object could therefore be not a black hole but a rotating gravastar. After comparing the real and imaginary parts of the ringdown signal of GW150914 with the corresponding quantities for a variety of gravastars, and notwithstanding the very limited knowledge of the perturbative response of rotating gravastars, we conclude it is not possible to model the measured ringdown of GW150914 as due to a rotating gravastar.

The most popular model for short gamma-ray bursts (sGRBs) involves the coalescence of binary neutron stars. Because the progenitor is actually hidden from view, we must consider under which circumstances such merging systems are capable of producing a successful sGRB. Soon after coalescence, winds are launched from the merger remnant. In this paper, we use realistic wind profiles derived from global merger simulations in order to investigate the interaction of sGRB jets with these winds using numerical simulations. We analyze the conditions for which these axisymmetric winds permit relativistic jets to breakout and produce a sGRB. We find that jets with luminosities comparable to those observed in sGRBs are only successful when their half-opening angles are below ~20{\deg}. This jet collimation mechanism leads to a simple physical interpretation of the luminosities and opening angles inferred for sGRBs. If wide, low luminosity jets are observed, they might be indicative of a different progenitor avenue such as the merger of a neutron star with a black hole. We also use the observed durations of sGRB to place constraints on the lifetime of the wind phase, which is determined by the time it takes the jet to breakout. In all cases we find that the derived limits argue against completely stable remnants for binary neutron star mergers that produce sGRBs.

We revisit the problem of predicting the spin magnitude and direction of the black hole resulting from the merger of two black holes with arbitrary masses and spins inspiralling in quasi-circular orbits. We do this by analyzing a catalog of 619 recent numerical-relativity simulations collected from the literature and spanning a large variety of initial conditions. By combining information from the post-Newtonian approximation, the extreme mass-ratio limit and perturbative calculations, we improve our previously proposed phenomenological formulae for the final remnant spin. In contrast with alternative suggestions in the literature, and in analogy with our previous expressions, the new formula is a simple algebraic function of the initial system parameters and is not restricted to binaries with spins aligned/anti-aligned with the orbital angular momentum, but can be employed for fully generic binaries. The accuracy of the new expression is significantly improved, especially for almost extremal progenitor spins and for small mass ratios, yielding a root-mean-square error $\sigma\approx0.002$ for aligned/anti-aligned binaries and $\sigma\approx0.006$ for generic binaries. Our new formula is suitable for cosmological applications and can be employed robustly in the analysis of the gravitational waveforms from advanced interferometric detectors.

A number of works have shown that important information on the equation of state of matter at nuclear density can be extracted from the gravitational waves emitted by merging neutron-star binaries. We present a comprehensive analysis of the gravitational-wave signal emitted during the inspiral, merger and post-merger of 56 neutron-star binaries. This sample of binaries, arguably the largest studied to date with realistic equations of state, spans across six different nuclear-physics equations of state and ten masses, allowing us to sharpen a number of results recently obtained on the spectral properties of the gravitational-wave signal. Overall we find that: (i) for binaries with masses differing no more than $20\%$, the frequency at gravitational-wave amplitude's maximum is related quasi-universally with the tidal deformability of the two stars; (ii) the spectral properties vary during the post-merger phase, with a transient phase lasting a few millisecond after the merger and followed by a quasi-stationary phase; (iii) when distinguishing the spectral peaks between these two phases, a number of ambiguities in the identification of the peaks disappear, leaving a simple and robust picture; (iv) using properly identified frequencies, quasi-universal relations are found between the spectral features and the properties of the neutron stars; (v) for the most salient peaks analytic fitting functions can be obtained in terms of the stellar tidal deformability or compactness. Altogether, these results support the idea that the equation of state of nuclear matter can be constrained tightly when a signal in gravitational waves from binary neutron stars is detected.

We present fully general-relativistic simulations of binary neutron star mergers with a temperature and composition dependent nuclear equation of state. We study the dynamical mass ejection from both quasi-circular and dynamical-capture eccentric mergers. We systematically vary the level of our treatment of the microphysics to isolate the effects of neutrino cooling and heating and we compute the nucleosynthetic yields of the ejecta. We find that eccentric binaries can eject significantly more material than quasi-circular binaries and generate bright infrared and radio emission. In all our simulations the outflow is composed of a combination of tidally- and shock-driven ejecta, mostly distributed over a broad $\sim 60^\circ$ angle from the orbital plane, and, to a lesser extent, by thermally driven winds at high latitudes. Ejecta from eccentric mergers are typically more neutron rich than those of quasi-circular mergers. We find neutrino cooling and heating to affect, quantitatively and qualitatively, composition, morphology, and total mass of the outflows. This is also reflected in the infrared and radio signatures of the binary. The final nucleosynthetic yields of the ejecta are robust and insensitive to input physics or merger type in the regions of the second and third r-process peaks. The yields for elements on the first peak vary between our simulations, but none of our models is able to explain the Solar abundances of first-peak elements without invoking additional first-peak contributions from either neutrino and viscously-driven winds operating on longer timescales after the mergers, or from core-collapse supernovae.

An important issue in the asteroseismology of compact and magnetized stars is the determination of the dissipation mechanism which is most efficient in damping the oscillations when these are produced. In a linear regime and for low-multipolarity modes these mechanisms are confined to either gravitational-wave or electromagnetic losses. We here consider the latter and compute the energy losses in the form of Poynting fluxes, Joule heating and Ohmic dissipation in a relativistic oscillating spherical star with a dipolar magnetic field in vacuum. While this approach is not particularly realistic for rapidly rotating stars, it has the advantage that it is fully analytic and that it provides expressions for the electric and magnetic fields produced by the most common modes of oscillation both in the vicinity of the star and far away from it. In this way we revisit and extend to a relativistic context the classical estimates of McDermott et al. Overall, we find that general-relativistic corrections lead to electromagnetic damping time-scales that are at least one order of magnitude smaller than in Newtonian gravity. Furthermore, with the only exception of $g$ (gravity) modes, we find that $f$ (fundamental), $p$ (pressure), $i$ (interface) and $s$ (shear) modes are suppressed more efficiently by gravitational losses than by electromagnetic ones.

We study the one-dimensional, longitudinally boost-invariant motion of an ideal fluid with infinite conductivity in the presence of a transverse magnetic field, i.e., in the ideal transverse magnetohydrodynamical limit. In an extension of our previous work Roy et al., [Phys. Lett. B 750, 45 (2015)], we consider the fluid to have a non-zero magnetization. First, we assume a constant magnetic susceptibility $\chi_{m}$ and consider an ultrarelativistic ideal gas equation of state. For a paramagnetic fluid (i.e., with $\chi_{m}>0$), the decay of the energy density slows down since the fluid gains energy from the magnetic field. For a diamagnetic fluid (i.e., with $\chi_{m}<0$), the energy density decays faster because it feeds energy into the magnetic field. Furthermore, when the magnetic field is taken to be external and to decay in proper time $\tau$ with a power law $\sim\tau^{-a}$, two distinct solutions can be found depending on the values of $a$ and $\chi_m$. Finally, we also solve the ideal magnetohydrodynamical equations for one-dimensional Bjorken flow with a temperature-dependent magnetic susceptibility and a realistic equation of state given by lattice-QCD data. We find that the temperature and energy density decay more slowly because of the non-vanishing magnetization. For values of the magnetic field typical for heavy-ion collisions, this effect is, however, rather small. Only for magnetic fields which are about an order of magnitude larger than expected for heavy-ion collisions, the system is substantially reheated and the lifetime of the quark phase might be extended.

A number of recent works have highlighted that it is possible to express the properties of general-relativistic stellar equilibrium configurations in terms of functions that do not depend on the specific equation of state employed to describe matter at nuclear densities. These functions are normally referred to as "universal relations" and have been found to apply, within limits, both to static or stationary isolated stars, as well as to fully dynamical and merging binary systems. Further extending the idea that universal relations can be valid also away from stability, we show that a universal relation is exhibited also by equilibrium solutions that are not stable. In particular, the mass of rotating configurations on the turning-point line shows a universal behaviour when expressed in terms of the normalised Keplerian angular momentum. In turn, this allows us to compute the maximum mass allowed by uniform rotation, M_{max}, simply in terms of the maximum mass of the nonrotating configuration, M_{TOV}, finding that M_{max} ~ (1.203 +- 0.022) M_{TOV} for all the equations of state we have considered. We further show that a universal relation can be found between the dimensionless moment of inertia and the stellar compactness. Although this relation is not surprising as it involves two quantities that have been shown to exhibit universal behaviour with other stellar properties, our parameterisation represents a refinement over a similar relation by Lattimer and Schutz (2005), where a different normalisation was used, and could provide an accurate tool to constrain the equation of state of nuclear matter when measurements of the moment of inertia become available.

Following previous work of ours in spherical symmetry, we here propose a new parametric framework to describe the spacetime of axisymmetric black holes in generic metric theories of gravity. In this case, the metric components are functions of both the radial and the polar angular coordinates, forcing a double expansion to obtain a generic axisymmetric metric expression. In particular, we use a continued-fraction expansion in terms of a compactified radial coordinate to express the radial dependence, while we exploit a Taylor expansion in terms of the cosine of the polar angle for the polar dependence. These choices lead to a superior convergence in the radial direction and to an exact limit on the equatorial plane. As a validation of our approach, we build parametrized representations of Kerr, rotating dilaton, and Einstein-dilaton-Gauss-Bonnet black holes. The match is already very good at lowest order in the expansion and improves as new orders are added. We expect a similar behavior for any stationary and axisymmetric black-hole metric.

We have studied the dynamics of an equal-mass magnetized neutron-star binary within a resistive magnetohydrodynamic (RMHD) approach in which the highly conducting stellar interior is matched to an electrovacuum exterior. Because our analysis is aimed at assessing the modifications introduced by resistive effects on the dynamics of the binary after the merger and through to collapse, we have carried out a close comparison with an equivalent simulation performed within the traditional ideal magnetohydrodynamic approximation. We have found that there are many similarities between the two evolutions but also one important difference: the survival time of the hyper massive neutron star increases in a RMHD simulation. This difference is due to a less efficient magnetic-braking mechanism in the resistive regime, in which matter can move across magnetic-field lines, thus reducing the outward transport of angular momentum. Both the RMHD and the ideal magnetohydrodynamic simulations carried here have been performed at higher resolutions and with a different grid structure than those in previous work of ours [L. Rezzolla, B. Giacomazzo, L. Baiotti, J. Granot, C. Kouveliotou, and M. A. Aloy, Astrophys. J. Letters 732, L6 (2011)], but confirm the formation of a low-density funnel with an ordered magnetic field produced by the black hole--torus system. In both regimes the magnetic field is predominantly toroidal in the highly conducting torus and predominantly poloidal in the nearly evacuated funnel. Reconnection processes or neutrino annihilation occurring in the funnel, none of which we model, could potentially increase the internal energy in the funnel and launch a relativistic outflow, which, however, is not produced in these simulations.

A large international effort is under way to assess the presence of a shadow in the radio emission from the compact source at the centre of our Galaxy, Sagittarius A$^*$ (Sgr A$^*$). If detected, this shadow would provide the first direct evidence of the existence of black holes and that Sgr A$^*$ is a supermassive black hole. In addition, the shape of the shadow could be used to learn about extreme gravity near the event horizon and to determine which theory of gravity better describes the observations. The mathematical description of the shadow has so far used a number of simplifying assumptions that are unlikely to be met by the real observational data. We here provide a general formalism to describe the shadow as an arbitrary polar curve expressed in terms of a Legendre expansion. Our formalism does not presume any knowledge of the properties of the shadow, e.g. the location of its centre, and offers a number of routes to characterize the distortions of the curve with respect to reference circles. These distortions can be implemented in a coordinate independent manner by different teams analysing the same data. We show that the new formalism provides an accurate and robust description of noisy observational data, with smaller error variances when compared to previous approaches for the measurement of the distortion.

In the initial stage of relativistic heavy-ion collisions, strong magnetic fields appear due to the large velocity of the colliding charges. The evolution of these fields appears as a novel and intriguing feature in the fluid-dynamical description of heavy-ion collisions. In this work, we study analytically the one-dimensional, longitudinally boost-invariant motion of an ideal fluid in the presence of a transverse magnetic field. Interestingly, we find that, in the limit of ideal magnetohydrodynamics, i.e., for infinite conductivity, and irrespective of the strength of the initial magnetization, the decay of the fluid energy density $e$ with proper time $\tau$ is the same as for the time-honored "Bjorken flow" without magnetic field. Furthermore, when the magnetic field is assumed to decay $\sim \tau^{-a}$, where $a$ is an arbitrary number, two classes of analytic solutions can be found depending on whether $a$ is larger or smaller than one. In summary, the analytic solutions presented here highlight that the Bjorken flow is far more general than formerly thought. These solutions can serve both to gain insight on the dynamics of heavy-ion collisions in the presence of strong magnetic fields and as testbeds for numerical codes.

We have performed two-dimensional special-relativistic magnetohydrodynamic simulations of non-equilibrium over-pressured relativistic jets in cylindrical geometry. Multiple stationary recollimation shock and rarefaction structures are produced along the jet by the nonlinear interaction of shocks and rarefaction waves excited at the interface between the jet and the surrounding ambient medium. Although initially the jet is kinematically dominated, we have considered axial, toroidal and helical magnetic fields to investigate the effects of different magnetic-field topologies and strengths on the recollimation structures. We find that an axial field introduces a larger effective gas-pressure and leads to stronger recollimation shocks and rarefactions, resulting in larger flow variations. The jet boost grows quadratically with the initial magnetic field. On the other hand, a toroidal field leads to weaker recollimation shocks and rarefactions, modifying significantly the jet structure after the first recollimation rarefaction and shock. The jet boost decreases systematically. For a helical field, instead, the behaviour depends on the magnetic pitch, with a phenomenology that ranges between the one seen for axial and toroidal magnetic fields, respectively. In general, however, a helical magnetic field yields a more complex shock and rarefaction substructure close to the inlet that significantly modifies the jet structure. The differences in shock structure resulting from different field configurations and strengths may have observable consequences for disturbances propagating through a stationary recollimation shock.

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.

Extending previous work by a number of authors, we have recently presented a new approach in which the detection of gravitational waves from merging neutron star binaries can be used to determine the equation of state of matter at nuclear density and hence the structure of neutron stars. In particular, after performing a large number of numerical-relativity simulations of binaries with nuclear equations of state, we have found that the post-merger emission is characterized by two distinct and robust spectral features. While the high-frequency peak was already shown to be associated with the oscillations of the hypermassive neutron star produced by the merger and to depend on the equation of state, we have highlighted that the low-frequency peak is related to the merger process and to the total compactness of the stars in the binary. This relation is essentially universal and provides a powerful tool to set tight constraints on the equation of state. We here provide additional information on the extensive analysis performed, illustrating the methods used, the tests considered, as well as the robustness of the results. We also discuss additional relations that can be deduced when exploring the data and how these correlate with various properties of the binary. Finally, we present a simple mechanical toy model that explains the main spectral features of the post-merger signal and can even reproduce analytically the complex waveforms emitted right after the merger.

We report simulations of the inspiral and merger of binary neutron stars performed with \texttt{WhiskyTHC}, the first of a new generation of numerical relativity codes employing higher than second-order methods for both the spacetime and the hydrodynamic evolution. We find that the use of higher-order schemes improves substantially the quality of the gravitational waveforms extracted from the simulations when compared to those computed using traditional second-order schemes. The reduced de-phasing and the faster convergence rate allow us to estimate the phase evolution of the gravitational waves emitted, as well as the magnitude of finite-resolution effects, without the need of phase- or time-alignments or rescalings of the waves, as sometimes done in other works. Furthermore, by using an additional unpublished simulation at very high resolution, we confirm the robustness of our high convergence order of $3.2$.

The merger of a binary of neutron stars provides natural explanations for many of the features of short gamma-ray bursts (SGRBs), such as the generation of a hot torus orbiting a rapidly rotating black hole, which can then build a magnetic jet and provide the energy reservoir to launch a relativistic outflow. Yet, this scenario has problems explaining the recently discovered long-term and sustained X-ray emission associated with the afterglows of a subclass of SGRBs. We propose a new model that explains how an X-ray afterglow can be sustained by the product of the merger and how the X-ray emission is produced before the corresponding emission in the gamma-band, although it is observed to follow it. Overall, our paradigm combines in a novel manner a number of well-established features of the emission in SGRBs and results from simulations. Because it involves the propagation of an ultra-relativistic outflow and its interaction with a confining medium, the paradigm also highlights a unifying phenomenology between short and long GRBs.

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.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

Determining the equation of state of matter at nuclear density and hence the structure of neutron stars has been a riddle for decades. We show how the imminent detection of gravitational waves from merging neutron star binaries can be used to solve this riddle. Using a large number of accurate numerical-relativity simulations of binaries with nuclear equations of state, we find that the postmerger emission is characterized by two distinct and robust spectral features. While the high-frequency peak has already been associated with the oscillations of the hypermassive neutron star produced by the merger and depends on the equation of state, a new correlation emerges between the low-frequency peak, related to the merger process, and the total compactness of the stars in the binary. More importantly, such a correlation is essentially universal, thus providing a powerful tool to set tight constraints on the equation of state. If the mass of the binary is known from the inspiral signal, the combined use of the two frequency peaks sets four simultaneous constraints to be satisfied. Ideally, even a single detection would be sufficient to select one equation of state over the others. We test our approach with simulated data and verify it works well for all the equations of state considered.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

We derive an analytic phenomenological expression that predicts the final mass of the black-hole remnant resulting from the merger of a generic binary system of black holes on quasi-circular orbits. Besides recovering the correct test-particle limit for extreme mass-ratio binaries, our formula reproduces well the results of all the numerical-relativity simulations published so far, both when applied at separations of a few gravitational radii, and when applied at separations of tens of thousands of gravitational radii. These validations make our formula a useful tool in a variety of contexts ranging from gravitational-wave physics to cosmology. As representative examples, we first illustrate how it can be used to decrease the phase error of the effective-one-body waveforms during the ringdown phase. Second, we show that, when combined with the recently computed self-force correction to the binding energy of nonspinning black-hole binaries, it provides an estimate of the energy emitted during the merger and ringdown. Finally, we use it to calculate the energy radiated in gravitational waves by massive black-hole binaries as a function of redshift, using different models for the seeds of the black-hole population.

Besides being among the most promising sources of gravitational waves, merging neutron star binaries also represent a leading scenario to explain the phenomenology of short gamma-ray bursts (SGRBs). Recent observations have revealed a large subclass of SGRBs with roughly constant luminosity in their X-ray afterglows, lasting $10\!-\!10^4$ s. These features are generally taken as evidence of a long-lived central engine powered by the magnetic spin-down of a uniformly rotating, magnetized object. We propose a different scenario in which the central engine powering the X-ray emission is a differentially rotating hypermassive neutron star (HMNS) that launches a quasi-isotropic and baryon-loaded wind driven by the magnetic field, which is built-up through differential rotation. Our model is supported by long-term, three-dimensional, general-relativistic, and ideal magnetohydrodynamic simulations, showing that this isotropic emission is a very robust feature. For a given HMNS, the presence of a collimated component depends sensitively on the initial magnetic field geometry, while the stationary electromagnetic luminosity depends only on the magnetic energy initially stored in the system. We show that our model is compatible with the observed timescales and luminosities and express the latter in terms of a simple scaling relation.

We present a new approach for achieving high-order convergence in fully general-relativistic hydrodynamic simulations. The approach is implemented in WhiskyTHC, a new code that makes use of state-of-the-art numerical schemes and was key in achieving, for the first time, higher than second-order convergence in the calculation of the gravitational radiation from inspiraling binary neutron stars Radice et al. (2013). Here, we give a detailed description of the algorithms employed and present results obtained for a series of classical tests involving isolated neutron stars. In addition, using the gravitational-wave emission from the late inspiral and merger of binary neutron stars, we make a detailed comparison between the results obtained with the new code and those obtained when using standard second-order schemes commonly employed for matter simulations in numerical relativity. We find that even at moderate resolutions and for binaries with large compactness, the phase accuracy is improved by a factor 50 or more.

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.

Recently, general relations among the quadrupole moment (Q), the moment of inertia (I), and the tidal deformability (Love number) of a neutron star were shown to exist. They are nearly independent of the nuclear matter equation of state and would be of great aid in extracting parameters from observed gravitational-waves and in testing general relativity. These relations, however, do not account for strong magnetic fields. We consider this problem by studying the effect of a strong magnetic field on slowly rotating relativistic neutron stars and show that, for simple magnetic field configurations that are purely poloidal or purely toroidal, the relation between Q and I is again nearly universal. However, different magnetic field geometries lead to different I-Q relations, and, in the case of a more realistic twisted-torus magnetic field configuration, the relation depends significantly on the equation of state, losing its universality. I-Love-Q relations must thus be used with very great care, since universality is lost for stars with long spin periods, i.e. P > 10 s, and strong magnetic fields, i.e. B > 10^12 G.

Several fast radio bursts have been discovered recently, showing a bright, highly dispersed millisecond radio pulse. The pulses do not repeat and are not associated with a known pulsar or gamma-ray burst. The high dispersion suggests sources at cosmological distances, hence implying an extremely high radio luminosity, far larger than the power of single pulses from a pulsar. We suggest that a fast radio burst represents the final signal of a supramassive rotating neutron star that collapses to a black hole due to magnetic braking. The neutron star is initially above the critical mass for non-rotating models and is supported by rapid rotation. As magnetic braking constantly reduces the spin, the neutron star will suddenly collapse to a black hole several thousand to million years after its birth. We discuss several formation scenarios for supramassive neutron stars and estimate the possible observational signatures {making use of the results of recent numerical general-relativistic calculations. While the collapse will hide the stellar surface behind an event horizon, the magnetic-field lines will snap violently. This can turn an almost ordinary pulsar into a bright radio "blitzar": Accelerated electrons from the travelling magnetic shock dissipate a significant fraction of the magnetosphere and produce a massive radio burst that is observable out to z>0.7. Only a few percent of the neutron stars needs to be supramassive in order to explain the observed rate. We suggest that fast radio bursts might trace the solitary formation of stellar mass black holes at high redshifts. These bursts could be an electromagnetic complement to gravitational-wave emission and reveal a new formation and evolutionary channel for black holes that are not seen as gamma-ray bursts. Radio observations of these bursts could trace the core-collapse supernova rate throughout the universe.

Inspiralling compact binaries as standard sirens will soon become an invaluable tool for cosmology when advanced interferometric gravitational-wave detectors begin their observations in the coming years. However, a degeneracy in the information carried by gravitational waves between the total rest-frame mass $M$ and the redshift $z$ of the source implies that neither can be directly extracted from the signal, but only the combination $M(1+z)$, the redshifted mass. Recent work has shown that for binary neutron star systems, a tidal correction to the gravitational-wave phase in the late-inspiral signal that depends on the rest-frame source mass could be used to break the mass-redshift degeneracy. We propose here to use the signature encoded in the post-merger signal to deduce the redshift to the source. This will allow an accurate extraction of the intrinsic rest-frame mass of the source, in turn permitting the determination of source redshift and luminosity distance solely from gravitational-wave observations. This will herald a new era in precision cosmography and astrophysics. Using numerical simulations of binary neutron star mergers of very slightly different mass, we model gravitational-wave signals at different redshifts and use Bayesian parameter estimation to determine the accuracy with which the redshift can be extracted for a source of known mass. We find that the Einstein Telescope can determine the source redshift to $sim 10$--$20%$ at redshifts of $z<0.04$.

Despite the recent rapid progress in numerical relativity, a convergence order less than the second has so far plagued codes solving the Einstein-Euler system of equations. We report simulations of the inspiral of binary neutron stars in quasi-circular orbits computed with a new code employing high-order, high-resolution shock-capturing, finite-differencing schemes that, for the first time, go beyond the second-order barrier. In particular, without any tuning or alignment, we measure a convergence order above three both in the phase and in the amplitude of the gravitational waves. Because the new code is able to calculate waveforms with very small phase errors already at modest resolutions, we are able to obtain accurate estimates of tidal effects in the inspiral that are essentially free from the large numerical viscosity typical of lower-order methods, and even for the challenging large compactness and small-deformability binary considered here. We find a remarkable agreement between our Richardson-extrapolated waveform and the one from the tidally corrected post-Newtonian (PN) Taylor-T4 model, with a de-phasing smaller than 0.2 radians during the seven orbits of the inspiral and up to the contact point. Because our results can be used reliably to assess the validity of the PN or other approximations at frequencies significantly larger than those considered so far in the literature, they seem to exclude at these compactnesses significant tidal amplifications from next-to-next-to-leading--order terms in the PN expansion.

We describe in detail the implementation of a simplified approach to radiative transfer in general relativity by means of the well-known neutrino leakage scheme (NLS). In particular, we carry out an extensive investigation of the properties and limitations of the NLS for isolated relativistic stars to a level of detail that has not been discussed before in a general-relativistic context. Although the numerous tests considered here are rather idealized, they provide a well-controlled environment in which to understand the relationship between the matter dynamics and the neutrino emission, which is important in order to model the neutrino signals from more complicated scenarios, such as binary neutron-star mergers. When considering nonrotating hot neutron stars we confirm earlier results of one-dimensional simulations, but also present novel results about the equilibrium properties and on how the cooling affects the stability of these configurations. In our idealized but controlled setup, we can then show that deviations from the thermal and weak-interaction equilibrium affect the stability of these models to radial perturbations, leading models that are stable in the absence of radiative losses, to a gravitational collapse to a black hole when neutrinos are instead radiated.

We show that magnetic fields stronger than about $10^{15}$ G are able to suppress the development of the hydrodynamical bar-mode instability in relativistic stars. The suppression is due to a change in the rest-mass density and angular velocity profiles due to the formation and to the linear growth of a toroidal component that rapidly overcomes the original poloidal one, leading to an amplification of the total magnetic energy. The study is carried out performing three-dimensional ideal-magnetohydrodynamics simulations in full general relativity, superimposing to the initial (matter) equilibrium configurations a purely poloidal magnetic field in the range $10^{14}-10^{16}$ G. When the seed field is a few parts in $10^{15}$ G or above, all the evolved models show the formation of a low-density envelope surrounding the star. For much weaker fields, no effect on the matter evolution is observed, while magnetic fields which are just below the suppression threshold are observed to slow down the growth-rate of the instability.

Understanding the properties of the internal magnetic field of neutron stars remains a theoretical challenge. Over the last years, twisted-torus geometries have been considered both in Newtonian and general-relativistic equilibrium models, as they represent a potentially good description of neutron star interiors. All of these works have found an apparent intrinsic limitation to geometries that are poloidal-field-dominated, with a toroidal-to-poloidal energy ratio inside the star that are <10%, unless surface currents are included and magnetic fields are allowed to be discontinuous. This limitation is in stark contrast with the general expectation that much higher toroidal fields should be present in the stellar interior and casts doubt about the stability and hence realism of these configurations. We here discuss how to overcome this limitation by adopting a new prescription for the azimuthal currents that leads to magnetized equilibria where the toroidal-to-total magnetic-field energy ratio can be as high as 90%, thus including geometries that are toroidal-field-dominated. Moreover, our results show that for a fixed exterior magnetic-field strength, a higher toroidal-field energy implies a much higher total magnetic energy stored in the star, with a potentially strong impact on the expected electromagnetic and gravitational-wave emission from highly magnetized neutron stars.

We present three-dimensional simulations of the dynamical bar-mode instability in magnetized and differentially rotating stars in full general relativity. Our focus is on the effects that magnetic fields have on the dynamics and the onset of the instability. In particular, we perform ideal-magnetohydrodynamics simulations of neutron stars that are known to be either stable or unstable against the purely hydrodynamical instability, but to which a poloidal magnetic field in the range of $10^{14}$--$10^{16}$ G is superimposed initially. As expected, the differential rotation is responsible for the shearing of the poloidal field and the consequent linear growth in time of the toroidal magnetic field. The latter rapidly exceeds in strength the original poloidal one, leading to a magnetic-field amplification in the the stars. Weak initial magnetic fields, i.e. $ \lesssim 10^{15}$ G, have negligible effects on the development of the dynamical bar-mode instability, simply braking the stellar configuration via magnetic-field shearing, and over a timescale for which we derived a simple algebraic expression. On the other hand, strong magnetic fields, i.e. $\gtrsim 10^{16}$ G, can suppress the instability completely, with the precise threshold being dependent also on the amount of rotation. As a result, it is unlikely that very highly magnetized neutron stars can be considered as sources of gravitational waves via the dynamical bar-mode instability.

A differentially rotating hypermassive neutron star (HMNS) is a metastable object which can be formed in the merger of neutron-star binaries. The eventual collapse of the HMNS into a black hole is a key element in generating the physical conditions expected to accompany the launch of a short gamma-ray burst. We investigate the influence of magnetic fields on HMNSs by performing three-dimensional simulations in general-relativistic magnetohydrodynamics. In particular, we provide direct evidence for the occurrence of the magnetorotational instability (MRI) in HMNS interiors. For the first time in simulations of these systems, rapidly-growing and spatially-periodic structures are observed to form with features like those of the channel flows produced by the MRI in other systems. Moreover, the growth time and wavelength of the fastest-growing mode are extracted and compared successfully with analytical predictions. The MRI emerges as an important mechanism to amplify magnetic fields over the lifetime of the HMNS, whose collapse to a black hole is accelerated. The evidence provided here that the MRI can actually develop in HMNSs could have a profound impact on the outcome of the merger of neutron-star binaries and on its connection to short gamma-ray bursts.

Using an extended set of equations of state and a multiple-group multiple-code collaborative effort to generate waveforms, we improve numerical-relativity-based data-analysis estimates of the measurability of matter effects in neutron-star binaries. We vary two parameters of a parameterized piecewise-polytropic equation of state (EOS) to analyze the measurability of EOS properties, via a parameter {\Lambda} that characterizes the quadrupole deformability of an isolated neutron star. We find that, to within the accuracy of the simulations, the departure of the waveform from point-particle (or spinless double black-hole binary) inspiral increases monotonically with {\Lambda}, and changes in the EOS that did not change {\Lambda} are not measurable. We estimate with two methods the minimal and expected measurability of {\Lambda} in second- and third- generation gravitational-wave detectors. The first estimate, using numerical waveforms alone, shows two EOS which vary in radius by 1.3km are distinguishable in mergers at 100Mpc. The second estimate relies on the construction of hybrid waveforms by matching to post-Newtonian inspiral, and estimates that the same EOS are distinguishable in mergers at 300Mpc. We calculate systematic errors arising from numerical uncertainties and hybrid construction, and we estimate the frequency at which such effects would interfere with template-based searches.

Recent work by McClarren & Hauck [29] suggests that the filtered spherical harmonics method represents an efficient, robust, and accurate method for radiation transport, at least in the two-dimensional (2D) case. We extend their work to the three-dimensional (3D) case and find that all of the advantages of the filtering approach identified in 2D are present also in the 3D case. We reformulate the filter operation in a way that is independent of the timestep and of the spatial discretization. We also explore different second- and fourth-order filters and find that the second-order ones yield significantly better results. Overall, our findings suggest that the filtered spherical harmonics approach represents a very promising method for 3D radiation transport calculations.

Numerical relativity has made big strides over the last decade. A number of problems that have plagued the field for years have now been mostly solved. This progress has transformed numerical relativity into a powerful tool to explore fundamental problems in physics and astrophysics, and I present here three representative examples. These "three little pieces" reflect a personal choice and describe work that I am particularly familiar with. However, many more examples could be made.

With the aim of determining the statistical properties of relativistic turbulence and unveiling novel and non-classical features, we resent the results of direct numerical simulations of driven turbulence in an ultrarelativistic hot plasma using high-order numerical schemes. We study the statistical properties of flows with average Mach number ranging from $\sim 0.4$ to $\sim 1.7$ and with average Lorentz factors up to $\sim 1.7$. We find that flow quantities, such as the energy density or the local Lorentz factor, show large spatial variance even in the subsonic case as compressibility is enhanced by relativistic effects. The velocity field is highly intermittent, but its power-spectrum is found to be in good agreement with the predictions of the classical theory of Kolmogorov. Overall, our results indicate that relativistic effects are able to significantly enhance the intermittency of the flow and affect the high-order statistics of the velocity field, while leaving unchanged the low-order statistics, which instead appear to be universal and in good agreement with the classical Kolmogorov theory. To the best of our knowledge, these are the most accurate simulations of driven relativistic turbulence to date.

In recent years, detailed observations and accurate numerical simulations have provided support to the idea that mergers of compact binaries containing either two neutron stars (NSs) or an NS and a black hole (BH) may constitute the central engine of short gamma-ray bursts (SGRBs). The merger of such compact binaries is expected to lead to the production of a spinning BH surrounded by an accreting torus. Several mechanisms can extract energy from this system and power the SGRBs. Here we connect observations and numerical simulations of compact binary mergers, and use the current sample of SGRBs with measured energies to constrain the mass of their powering tori. By comparing the masses of the tori with the results of fully general-relativistic simulations, we are able to infer the properties of the binary progenitors which yield SGRBs. By assuming a constant efficiency in converting torus mass into jet energy, epsilon_{jet}=10%, we find that most of the tori have masses smaller than 0.01M_{sun}, favoring "high-mass" binary NSs mergers, i.e., binaries with total masses >~1.5 the maximum mass of an isolated NS. This has important consequences for the gravitational-wave signals that may be detected in association with SGRBs, since "high-mass" systems do not form a long-lived hypermassive NS after the merger. While NS-BH systems cannot be excluded to be the engine of at least some of the SGRBs, the BH would need to have an initial spin of ~0.9, or higher.

Determining the conditions under which a black hole can be produced is a long-standing and fundamental problem in general relativity. We use numerical simulations of colliding selfgravitating fluid objects to study the conditions of black-hole formation when the objects are boosted to ultrarelativistic speeds. Expanding on previous work, we show that the collision is characterized by a type-I critical behaviour, with a black hole being produced for masses above a critical value, M_c, and a partially bound object for masses below the critical one. More importantly, we show for the first time that the critical mass varies with the initial effective Lorentz factor <\gamma> following a simple scaling of the type M_c ~ K <\gamma>^{-1.0}, thus indicating that a black hole of infinitesimal mass is produced in the limit of a diverging Lorentz factor. Furthermore, because a scaling is present also in terms of the initial stellar compactness, we provide a condition for black-hole formation in the spirit of the hoop conjecture.

We investigate the instability of purely poloidal magnetic fields in nonrotating neutron stars by means of three-dimensional general-relativistic magnetohydrodynamics simulations, extending the work presented in Ciolfi et al. (2011). Our aim is to draw a clear picture of the dynamics associated with the instability and to study the final configuration reached by the system, thus obtaining indications on possible equilibria in a magnetized neutron star. Furthermore, since the internal rearrangement of magnetic fields is a highly dynamical process, which has been suggested to be behind magnetar giant flares, our simulations can provide a realistic estimate of the electromagnetic and gravitational-wave emission which should accompany the flare event. Our main findings are the following: (i) the initial development of the instability meets all the expectations of perturbative studies in terms of the location of the seed of the instability, the timescale for its growth and the generation of a toroidal component; (ii) in the subsequent nonlinear reorganization of the system, ~90% of magnetic energy is lost in few Alfven timescales mainly through electromagnetic emission, and further decreases on a much longer timescale; (iii) all stellar models tend to achieve a significant amount of magnetic helicity and the equipartition of energy between poloidal and toroidal magnetic fields, and evolve to a new configuration which does not show a subsequent instability on dynamical or Alfven timescales; (iv) the electromagnetic emission matches the duration of the initial burst in luminosity observed in giant flares, giving support to the internal rearrangement scenario; (v) only a small fraction of the energy released during the process is converted into f-mode oscillations and in the consequent gravitational-wave emission, thus resulting in very low chances of detecting this signal with present and..

With the aim of determining the statistical properties of relativistic turbulence and unveiling novel and non-classical features, we present the results of direct numerical simulations of driven turbulence in an ultrarelativistic hot plasma using high-order numerical schemes. We study the statistical properties of flows with average Mach number ranging from $\sim 0.4$ to $\sim 1.7$ and with average Lorentz factors up to $\sim 1.7$. We find that flow quantities, such as the energy density or the local Lorentz factor, show large spatial variance even in the subsonic case as compressibility is enhanced by relativistic effects. The velocity field is highly intermittent, but its power-spectrum is found to be in good agreement with the predictions of the classical theory of Kolmogorov.

We present THC: a new high-order flux-vector-splitting code for Newtonian and special-relativistic hydrodynamics designed for direct numerical simulations of turbulent flows. Our code implements a variety of different reconstruction algorithms, such as the popular weighted essentially non oscillatory and monotonicity-preserving schemes, or the more specialised bandwidth-optimised WENO scheme that has been specifically designed for the study of compressible turbulence. We show the first systematic comparison of these schemes in Newtonian physics as well as for special-relativistic flows. In particular we will present the results obtained in simulations of grid-aligned and oblique shock waves and nonlinear, large-amplitude, smooth adiabatic waves. We will also discuss the results obtained in classical benchmarks such as the double-Mach shock reflection test in Newtonian physics or the linear and nonlinear development of the relativistic Kelvin-Helmholtz instability in two and three dimensions. Finally, we study the turbulent flow induced by the Kelvin-Helmholtz instability and we show that our code is able to obtain well-converged velocity spectra, from which we benchmark the effective resolution of the different schemes.

Einstein gravitational-wave Telescope (ET) is a design study funded by the European Commission to explore the technological challenges of and scientific benefits from building a third generation gravitational wave detector. The three-year study, which concluded earlier this year, has formulated the conceptual design of an observatory that can support the implementation of new technology for the next two to three decades. The goal of this talk is to introduce the audience to the overall aims and objectives of the project and to enumerate ET's potential to influence our understanding of fundamental physics, astrophysics and cosmology.

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 present new results on the dynamics and gravitational-wave emission from the collapse of differentially rotating neutron stars. We have considered a number of polytropic stellar models having different values of the dimensionless angular momentum J/M^2, where J and M are the asymptotic angular momentum and mass of the star, respectively. For neutron stars with J/M^2<1, i.e., "sub-Kerr" models, we were able to find models that are dynamically unstable and that collapse promptly to a rotating black hole. Both the dynamics of the collapse and the consequent emission of gravitational waves resemble the one seen for uniformly rotating stars, although with an overall decrease in the efficiency of gravitational-wave emission. For stellar models with J/M^2>1, i.e., "supra-Kerr" models, on the other hand, we were not able to find models that are dynamically unstable and all of the computed supra-Kerr models were found to be far from the stability threshold. For these models a gravitational collapse is possible only after a very severe and artificial reduction of the pressure, which then leads to a torus developing nonaxisymmetric instabilities and eventually contracting to a stable axisymmetric stellar configuration. While this does not exclude the possibility that a naked singularity can be produced by the collapse of a differentially rotating star, it also suggests that cosmic censorship is not violated and that generic conditions for a supra-Kerr progenitor do not lead to a naked singularity.

The stability properties of relativistic stars against gravitational collapse to black hole is a classical problem in general relativity. A sufficient criterion for secular instability was established by Friedman, Ipser and Sorkin (1988), who proved that a sequence of uniformly rotating barotropic stars is secularly unstable on one side of a turning point and then argued that a stronger result should hold: that the sequence should be stable on the opposite side, with the turning point marking the onset of secular instability. We show here that this expectation is not met. By computing in full general relativity the $F$-mode frequency for a large number of rotating stars, we show that the neutral-stability point, i.e., where the frequency becomes zero, differs from the turning point for rotating stars. Using numerical simulations we validate that the new criterion can be used to assess the dynamical stability of relativistic rotating stars.

We have developed the formalism necessary to employ the discontinuous-Galerkin approach in general-relativistic hydrodynamics. The formalism is firstly presented in a general 4-dimensional setting and then specialized to the case of spherical symmetry within a 3+1 splitting of spacetime. As a direct application, we have constructed a one-dimensional code, EDGES, which has been used to asses the viability of these methods via a series of tests involving highly relativistic flows in strong gravity. Our results show that discontinuous Galerkin methods are able not only to handle strong relativistic shock waves but, at the same time, to attain very high orders of accuracy and exponential convergence rates in smooth regions of the flow. Given these promising prospects and their affinity with a pseudospectral solution of the Einstein equations, discontinuous Galerkin methods could represent a new paradigm for the accurate numerical modelling in relativistic astrophysics.

The problem of the stability of magnetic fields in stars has a long history and has been investigated in detail in perturbation theory. Here we consider the nonlinear evolution of a non-rotating neutron star with a purely poloidal magnetic field, in general relativity. We find that an instability develops in the region of the closed magnetic field lines and over an Alfven timescale, as predicted by perturbation theory. After the initial unstable growth, our evolutions show that a toroidal magnetic field component is generated, which increases until it is locally comparable in strength with the poloidal one. On longer timescales the system relaxes to a new non-axisymmetric configuration with a reorganization of the stellar structure and large-amplitude oscillations, mostly in the fundamental mode. We discuss the energies involved in the instability and the impact they may have on the phenomenology of magnetar flares and on their detectability through gravitational-wave emission.

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.

Short Gamma-Ray Bursts (SGRBs) are among the most luminous explosions in the universe, releasing in less than one second the energy emitted by our Galaxy over one year. Despite decades of observations, the nature of their "central-engine" remains unknown. Considering a binary of magnetized neutron stars and solving Einstein equations, we show that their merger results in a rapidly spinning black hole surrounded by a hot and highly magnetized torus. Lasting over 35 ms and much longer than previous simulations, our study reveals that magnetohydrodynamical instabilities amplify an initially turbulent magnetic field of ~ 10^{12} G to produce an ordered poloidal field of ~ 10^{15} G along the black-hole spin-axis, within a half-opening angle of ~ 30 deg, which may naturally launch a relativistic jet. The broad consistency of our ab-initio calculations with SGRB observations shows that the merger of magnetized neutron stars can provide the basic physical conditions for the central-engine of SGRBs.

The strong tidal forces that arise during the last stages of the life of a black hole-neutron star binary may severely distort, and possibly disrupt, the star. Both phenomena will imprint signatures about the stellar structure in the emitted gravitational radiation. The information from the disruption, however, is confined to very high frequencies, where detectors are not very sensitive. We thus assess whether the lack of tidal distortion corrections in data-analysis pipelines will affect the detection of the inspiral part of the signal and whether these may yield information on the equation of state of matter at nuclear densities. Using recent post-Newtonian expressions and realistic equations of state to model these scenarios, we find that point-particle templates are sufficient for the detection of black hole-neutron star inspiralling binaries, with a loss of signals below 1% for both second and third-generation detectors. Such detections may be able to constrain particularly stiff equations of state, but will be unable to reveal the presence of a neutron star with a soft equation of state.

By performing new, long and numerically accurate general-relativistic simulations of magnetized, equal-mass neutron-star binaries, we investigate the role that realistic magnetic fields may have in the evolution of these systems. In particular, we study the evolution of the magnetic fields and show that they can influence the survival of the hypermassive-neutron star produced at the merger by accelerating its collapse to a black hole. We also provide evidence that even if purely poloidal initially, the magnetic fields produced in the tori surrounding the black hole have toroidal and poloidal components of equivalent strength. When estimating the possibility that magnetic fields could have an impact on the gravitational-wave signals emitted by these systems either during the inspiral or after the merger we conclude that for realistic magnetic-field strengths B<~1e12 G such effects could be detected, but only marginally, by detectors such as advanced LIGO or advanced Virgo. However, magnetically induced modifications could become detectable in the case of small-mass binaries and with the development of gravitational-wave detectors, such as the Einstein Telescope, with much higher sensitivities at frequencies larger than ~2 kHz.

The merger of a binary system composed of a black hole and a neutron star may leave behind a torus of hot, dense matter orbiting around the black hole. While numerical-relativity simulations are necessary to simulate this process accurately, they are also computationally expensive and unable at present to cover the large space of possible parameters, which include the relative mass ratio, the stellar compactness, and the black hole spin. To mitigate this and provide a first reasonable coverage of the space of parameters, we have developed a method for estimating the mass of the remnant torus from black hole-neutron star mergers. The toy model makes use of an improved relativistic affine model to describe the tidal deformations of an extended tri-axial ellipsoid orbiting around a Kerr black hole and measures the mass of the remnant torus by considering which of the fluid particles composing the star are on bound orbits at the time of the tidal disruption. We tune the toy model by using the results of fully general-relativistic simulations obtaining relative precisions of a few percent and use it to extensively investigate the space of parameters. In this way we find that the torus mass is largest for systems with highly spinning black holes, small stellar compactnesses, and large mass ratios. As an example, tori as massive as ~1.33 solar masses can be produced for a very extended star with compactness of ~0.1 inspiralling around a black hole with dimensionless spin equal to 0.85 and mass ratio of about 0.3. However, for a more astrophysically reasonable mass ratio of ~0.14 and a canonical value of the stellar compactness of ~0.145, the toy model sets a considerably smaller upper limit to the torus mass of less than ~0.34 solar masses.

The existence and detection of scalar fields could provide solutions to long-standing puzzles about the nature of dark matter, the dark compact objects at the centre of most galaxies, and other phenomena. Yet, self-interacting scalar fields are very poorly constrained by astronomical observations, leading to great uncertainties in estimates of the mass $m_\phi$ and the self-interacting coupling constant $\lambda$ of these fields. To counter this, we have systematically employed available astronomical observations to develop new constraints, considerably restricting this parameter space. In particular, by exploiting precise observations of stellar dynamics at the centre of our Galaxy and assuming that these dynamics can be explained by a single boson star, we determine an upper limit for the boson star compactness and impose significant limits on the values of the properties of possible scalar fields. Requiring the scalar field particle to follow a collisional dark matter model further narrows these constraints. Most importantly, we find that if a scalar dark matter particle does exist, then it cannot account for both the dark-matter halos and the existence of dark compact objects in galactic nuclei

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.

The generation of a large recoil velocity from the inspiral and merger of binary black holes represents one of the most exciting results of numerical-relativity calculations. While many aspects of this process have been investigated and explained, the "antikick", namely the sudden deceleration after the merger, has not yet found a simple explanation. We show that the antikick can be understood in terms of the radiation from a deformed black hole where the anisotropic curvature distribution on the horizon correlates with the direction and intensity of the recoil. Our analysis is focussed on Robinson-Trautman spacetimes and allows us to measure both the energies and momenta radiated in a gauge-invariant manner. At the same time, this simpler setup provides the qualitative and quantitative features of merging black holes, opening the way to a deeper understanding of the nonlinear dynamics of black-hole spacetimes.

We present new results on instabilities in rapidly and differentially rotating neutron stars. We model the stars in full general relativity and describe the stellar matter adopting a cold realistic equation of state based on the unified SLy prescription. We provide evidence that rapidly and differentially rotating stars that are below the expected threshold for the dynamical bar-mode instability, beta_c = T/|W| ~ 0.25, do nevertheless develop a shear instability on a dynamical timescale and for a wide range of values of beta. This class of instability, which has so far been found only for small values of beta and with very small growth rates, is therefore more generic than previously found and potentially more effective in producing strong sources of gravitational waves. Overall, our findings support the phenomenological predictions made by Watts, Andersson and Jones on the nature of the low-T/|W|.

We present new results from accurate and fully general-relativistic simulations of the coalescence of unmagnetized binary neutron stars with various mass ratios. The evolution of the stars is followed through the inspiral phase, the merger and prompt collapse to a black hole, up until the appearance of a thick accretion disk, which is studied as it enters and remains in a regime of quasi-steady accretion. Although a simple ideal-fluid equation of state with \Gamma=2 is used, this work presents a systematic study within a fully general relativistic framework of the properties of the resulting black-hole--torus system produced by the merger of unequal-mass binaries. More specifically, we show that: (1) The mass of the torus increases considerably with the mass asymmetry and equal-mass binaries do not produce significant tori if they have a total baryonic mass M_tot >~ 3.7 M_sun; (2) Tori with masses M_tor ~ 0.2 M_sun are measured for binaries with M_tot ~ 3.4 M_sun and mass ratios q ~ 0.75-0.85; (3) The mass of the torus can be estimated by the simple expression M_tor(q, M_tot) = [c_1 (1-q) + c_2](M_max-M_tot), involving the maximum mass for the binaries and coefficients constrained from the simulations, and suggesting that the tori can have masses as large as M_tor ~ 0.35 M_sun for M_tot ~ 2.8 M_sun and q ~ 0.75-0.85; (4) Using a novel technique to analyze the evolution of the tori we find no evidence for the onset of non-axisymmetric instabilities and that very little, if any, of their mass is unbound; (5) Finally, for all the binaries considered we compute the complete gravitational waveforms and the recoils imparted to the black holes, discussing the prospects of detection of these sources for a number of present and future detectors.

Binary black-hole systems with spins aligned or anti-aligned to the orbital angular momentum provide the natural ground to start detailed studies of the influence of strong-field spin effects on gravitational wave observations of coalescing binaries. Furthermore, such systems may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins parallel to the orbital angular momentum but otherwise arbitrary. Our attention is particularly focused on the gravitational-wave emission so as to quantify how much spin effects contribute to the signal-to-noise ratio, to the horizon distances, and to the relative event rates for the representative ranges in masses and detectors. As expected, the signal-to-noise ratio increases with the projection of the total black hole spin in the direction of the orbital momentum. We find that equal-spin binaries with maximum spin aligned with the orbital angular momentum are more than "three times as loud" as the corresponding binaries with anti-aligned spins, thus corresponding to event rates up to 30 times larger. We also consider the waveform mismatch between the different spinning configurations and find that, within our numerical accuracy, binaries with opposite spins S_1=-S_2 cannot be distinguished whereas binaries with spin S_1=S_2 have clearly distinct gravitational-wave emissions. Finally, we derive a simple expression for the energy radiated in gravitational waves and find that the binaries always have efficiencies E_rad/M > 3.6%, which can become as large as E_rad/M = 10% for maximally spinning binaries with spins aligned with the orbital angular momentum.

As one step towards a systematic modeling of the electromagnetic (EM) emission from an inspiralling black hole binary we consider a simple scenario in which the binary moves in a uniform magnetic field anchored to a distant circumbinary disc. We study this system by solving the Einstein-Maxwell equations in which the EM fields are chosen with astrophysically consistent strengths. We consider binaries with spins aligned or anti-aligned with the orbital angular momentum and study the dependence of gravitational and EM signals with these spin configurations. Overall we find that the EM radiation in the lowest l=2, m=2 multipole accurately reflects the gravitational one, with identical phase evolutions and amplitudes that differ only by a scaling factor. We also compute the efficiency of the energy emission in EM waves and find that it is given by E^rad_EM/M ~ 10^-15 (M/10^8 M_Sun)^2 (B/10^4 G)^2, hence 13 orders of magnitude smaller than the gravitational energy for realistic magnetic fields. The corresponding luminosity is much smaller than the accretion luminosity if the system is accreting at near the Eddington rate. Most importantly, this EM emission is at frequencies of 10^-4 (10^8 M_Sun/M) Hz, well outside those accessible to astronomical radio observations. As a result, it is unlikely that the EM emission discussed here can be detected directly and simultaneously with the gravitational-wave one. However, indirect processes, driven by changes in the EM fields behavior could yield observable events. In particular if the accretion rate of the circumbinary disc is small and sufficiently stable over the timescale of the final inspiral, then the EM emission may be observable indirectly as it will alter the accretion rate through the magnetic torques exerted by the distorted magnetic field lines.

We discuss different ways that neutron stars can generate gravitational waves, describe recent improvements in modelling the relevant scenarios in the context of improving detector sensitivity, and show how observations are beginning to test our understanding of fundamental physics. The main purpose of the discussion is to establish promising science goals for third-generation ground-based detectors, like the Einstein Telescope, and identify the various challenges that need to be met if we want to use gravitational-wave data to probe neutron star physics.

(Abridged.) The accretion-induced collapse (AIC) of a white dwarf (WD) may lead to the formation of a protoneutron star and a collapse-driven supernova explosion. This process represents a path alternative to thermonuclear disruption of accreting white dwarfs in Type Ia supernovae. Neutrino and gravitational-wave (GW) observations may provide crucial information necessary to reveal a potential AIC. Motivated by the need for systematic predictions of the GW signature of AIC, we present results from an extensive set of general-relativistic AIC simulations using a microphysical finite-temperature equation of state and an approximate treatment of deleptonization during collapse. Investigating a set of 114 progenitor models in rotational equilibrium, with a wide range of rotational configurations, temperatures and central densities, we extend previous Newtonian studies and find that the GW signal has a generic shape akin to what is known as a "Type III" signal in the literature. We discuss the detectability of the emitted GWs, showing that the signal-to-noise ratio for current or next-generation interferometer detectors could be high enough to detect such events in our Galaxy. Some of our AIC models form massive quasi-Keplerian accretion disks after bounce. In rapidly differentially rotating models, the disk mass can be as large as ~0.8-Msun. Slowly and/or uniformly rotating models produce much smaller disks. Finally, we find that the postbounce cores of rapidly spinning white dwarfs can reach sufficiently rapid rotation to develop a nonaxisymmetric rotational instability.

Knowledge of the spin of the black hole resulting from the merger of a generic black-hole binary is of great importance for studying the cosmological evolution of supermassive black holes. Several attempts have been made to model the spin via simple expressions exploiting the results of numerical-relativity simulations. While these expressions are in reasonable agreement with the simulations, they neglect the precession of the binary's orbital plane, and cannot therefore be applied directly -- i.e., without evolving the system to small separations using post-Newtonian theory -- to binaries with separations larger than a few hundred gravitational radii. While not a problem in principle, this may be impractical if the formulas are employed in cosmological merger-trees or N-body simulations, which provide the spins and angular momentum of the two black holes when their separation is of hundreds or thousands of gravitational radii. The formula that we propose is instead built on improved assumptions and gives, for any separation, a very accurate prediction both for the norm of the final spin and for its direction. By comparing with the numerical data, we also show that the final-spin direction is very accurately aligned with the binary's total angular momentum at large separation. Hence, observations of the final-spin direction (e.g. via a jet) can provide information on the binary's orbital plane at large separations and could be relevant, for instance, for studying X-shaped radio sources.

Using accurate and fully general-relativistic simulations we assess the effect that magnetic fields have on the gravitational-wave emission produced during the inspiral and merger of magnetized neutron stars. In particular, we show that magnetic fields have an impact after the merger, because they are amplified by a Kelvin-Helmholtz instability, but also during the inspiral, most likely because the magnetic tension reduces the stellar tidal deformation for extremely large initial magnetic fields, B_0>~10^{17}G. We quantify the influence of magnetic fields by computing the overlap, O, between the waveforms produced during the inspiral by magnetized and unmagnetized binaries. We find that for any realistic magnetic field strength B_0<~10^{14}G the overlap during the inspiral is O>~0.999 and is quite insensitive to the mass of the neutron stars. Only for unrealistically large magnetic fields like B_0~10^{17}G the overlap does decrease noticeably, becoming at our resolutions O<~0.76/0.67 for stars with baryon masses M_b~1.4/1.6 Msun, respectively. Because neutron stars are expected to merge with magnetic fields ~10^{8}-10^{10}G and because present detectors are sensitive to O<~0.995, we conclude that it is very unlikely that the present detectors will be able to discern the presence of magnetic fields during the inspiral of neutron stars.

We have recently presented an investigation in full general relativity of the dynamics and gravitational-wave emission from binary neutron stars which inspiral and merge, producing a black hole surrounded by a torus (see arXiv:0804.0594). We here discuss in more detail the convergence properties of the results presented in arXiv:0804.0594 and, in particular, the deterioration of the convergence rate at the merger and during the survival of the merged object, when strong shocks are formed and turbulence develops. We also show that physically reasonable and numerically convergent results obtained at low-resolution suffer however from large truncation errors and hence are of little physical use. We summarize our findings in an "error budget", which includes the different sources of possible inaccuracies we have investigated and provides a first quantitative assessment of the precision in the modelling of compact fluid binaries.

Gravitational waveforms from the inspiral and ring-down stages of the binary black hole coalescences can be modelled accurately by approximation/perturbation techniques in general relativity. Recent progress in numerical relativity has enabled us to model also the non-perturbative merger phase of the binary black-hole coalescence problem. This enables us to \emph{coherently} search for all three stages of the coalescence of non-spinning binary black holes using a single template bank. Taking our motivation from these results, we propose a family of template waveforms which can model the inspiral, merger, and ring-down stages of the coalescence of non-spinning binary black holes that follow quasi-circular inspiral. This two-dimensional template family is explicitly parametrized by the physical parameters of the binary. We show that the template family is not only \emph{effectual} in detecting the signals from black hole coalescences, but also \emph{faithful} in estimating the parameters of the binary. We compare the sensitivity of a search (in the context of different ground-based interferometers) using all three stages of the black hole coalescence with other template-based searches which look for individual stages separately. We find that the proposed search is significantly more sensitive than other template-based searches for a substantial mass-range, potentially bringing about remarkable improvement in the event-rate of ground-based interferometers. As part of this work, we also prescribe a general procedure to construct interpolated template banks using non-spinning black hole waveforms produced by numerical relativity.

An increase in the central density of a neutron star may trigger a phase transition from hadronic matter to deconfined quark matter in the core, causing it to collapse to a more compact hybrid-star configuration. We present a study of this, building on previous work by Lin et al. (2006). We follow them in considering a supersonic phase transition and using a simplified equation of state, but our calculations are general relativistic (using 2D simulations in the conformally flat approximation) as compared with their 3D Newtonian treatment. We also improved the treatment of the initial phase transformation, avoiding the introduction of artificial convection. As before, we find that the emitted gravitational-wave spectrum is dominated by the fundamental quasi-radial and quadrupolar pulsation modes but the strain amplitudes are much smaller than suggested previously, which is disappointing for the detection prospects. However, we see significantly smaller damping and observe a nonlinear mode resonance which substantially enhances the emission in some cases. We explain the damping mechanisms operating, giving a different view from the previous work. Finally, we discuss the detectability of the gravitational waves, showing that the signal-to-noise ratio for current or second generation interferometers could be high enough to detect such events in our Galaxy, although third generation detectors would be needed to observe them out to the Virgo cluster, which would be necessary for having a reasonable event rate.

Over the last few years enormous progress has been made in the numerical description of the inspiral and merger of binary black holes. A particular effort has gone into the modelling of the physical properties of the final black hole, namely its spin and recoil velocity, as these quantities have direct impact in astrophysics, cosmology and, of course, general relativity. As numerical-relativity calculations still remain computationally very expensive and cannot be used to investigate the complete space of possible parameters, semi-analytic approaches have been developed and shown to reproduce with very high precision the numerical results. I here collect and review these efforts, pointing out the relative strengths and weaknesses, and discuss which directions are more promising to further improve them.

Binary neutron-star (BNS) systems represent primary sources for the gravitational-wave (GW) detectors. We present a systematic investigation in full GR of the dynamics and GW emission from BNS which inspiral and merge, producing a black hole (BH) surrounded by a torus. Our results represent the state of the art from several points of view: (i) We use HRSC methods for the hydrodynamics equations and high-order finite-differencing techniques for the Einstein equations; (ii) We employ AMR techniques with "moving boxes"; (iii) We use as initial data BNSs in irrotational quasi-circular orbits; (iv) We exploit the isolated-horizon formalism to measure the properties of the BHs produced in the merger; (v) Finally, we use two approaches, based either on gauge-invariant perturbations or on Weyl scalars, to calculate the GWs. These techniques allow us to perform accurate evolutions on timescales never reported before (ie ~30 ms) and to provide the first complete description of the inspiral and merger of a BNS leading to the prompt or delayed formation of a BH and to its ringdown. We consider either a polytropic or an ideal fluid EOS and show that already with this idealized EOSs a very interesting phenomenology emerges. In particular, we show that while high-mass binaries lead to the prompt formation of a rapidly rotating BH surrounded by a dense torus, lower-mass binaries give rise to a differentially rotating NS, which undergoes large oscillations and emits large amounts of GWs. Eventually, also the NS collapses to a rotating BH surrounded by a torus. Finally, we also show that the use of a non-isentropic EOS leads to significantly different evolutions, giving rise to a delayed collapse also with high-mass binaries, as well as to a more intense emission of GWs and to a geometrically thicker torus.

Many astrophysical processes involving magnetic fields and quasi-stationary processes are well described when assuming the fluid as a perfect conductor. For these systems, the ideal-magnetohydrodynamics (MHD) description captures the dynamics effectively and a number of well-tested techniques exist for its numerical solution. Yet, there are several astrophysical processes involving magnetic fields which are highly dynamical and for which resistive effects can play an important role. The numerical modeling of such non-ideal MHD flows is significantly more challenging as the resistivity is expected to change of several orders of magnitude across the flow and the equations are then either of hyperbolic-parabolic nature or hyperbolic with stiff terms. We here present a novel approach for the solution of these relativistic resistive MHD equations exploiting the properties of implicit-explicit (IMEX) Runge Kutta methods. By examining a number of tests we illustrate the accuracy of our approach under a variety of conditions and highlight its robustness when compared with alternative methods, such as the Strang-splitting. Most importantly, we show that our approach allows one to treat, within a unified framework, both those regions of the flow which are fluid-pressure dominated (such as in the interior of compact objects) and those which are instead magnetic-pressure dominated (such as in their magnetospheres)

We provide a compact analytic formula to compute the spin of the black hole produced by the coalescence of two black holes. The expression, which uses an analytic fit of numerical-relativity data and relies on four assumptions, aims at modelling generic initial spin configurations and mass ratios. A comparison with numerical-relativity simulations already shows very accurate agreements with all of the numerical data available to date, but we also suggest a number of ways in which our predictions can be further improved.

We have studied extreme mass-ratio inspirals (EMRIs) in spacetimes containing a rotating black hole and a non self-gravitating torus with constant specific angular momentum. We have found that the effect of the hydrodynamic drag exerted by the torus on the satellite is much smaller than the corresponding one due to radiation reaction, for systems such as those generically expected in AGNs and at distances from the SMBH which can be probed with LISA. However, given the uncertainty on the parameters of these systems, there exist configurations in which the effect of the hydrodynamic drag can be comparable to the radiation-reaction one in phases of the inspiral which are detectable by LISA. This is the case, for instance, for a 10^6 M_sun SMBH surrounded by a corotating torus of comparable mass and with radius of 10^3-10^4 gravitational radii, or for a 10^5 M_sun SMBH surrounded by a corotating 10^4 M_sun torus with radius of 10^5 gravitational radii. Should these conditions be met in astrophysical systems, EMRI-gravitational waves could provide a characteristic signature of the presence of the torus. In fact, while radiation reaction always increases the inclination of the orbit with respect to the equatorial plane, the hydrodynamic drag from a torus corotating with the SMBH always decreases it. However, even when initially dominating over radiation reaction, the influence of the hydrodynamic drag decays very rapidly as the satellite moves into the very strong-field region of the SMBH (i.e., p <~ 5M). Although our results have been obtained for a specific class of tori, we argue that they will be qualitatively valid also for more generic distributions of the specific angular momentum.

During the first 40 s after their birth, proto-neutron stars are expected to be subject to at least two types of instability: the convective instability and the neutron-finger one. Both instabilities involve convective motions and hence can trigger dynamo actions which may be responsible for the large magnetic fields in neutron stars and magnetars. We have solved the mean-field induction equation in a simplified one-dimensional model of both the convective and the neutron-finger instability zones. Although very idealized, the model includes the nonlinearities introduced by the feedback processes which tend to saturate the growth of the magnetic field (alpha-quenching) and suppress its turbulent diffusion (eta-quenching). The possibility of a dynamo action is studied within a dynamical model of turbulent diffusivity where the boundary of the unstable zone is allowed to move. We show that the dynamo action can be operative and that the amplification of the magnetic field can still be very effective. Furthermore, we confirm the existence of a critical spin-period, below which the dynamo is always excited independently of the degree of differential rotation, and whose value is related to the size of the neutron-finger instability zone. Finally we provide a relation for the intensity of the final field as a function of the spin of the star and of its differential rotation. Although they were obtained by using a toy model, we expect that our results are able to capture the qualitative and asymptotic behaviour of a mean-field dynamo action developing in the neutron-finger instability zone. Overall, we find that such a dynamo is very efficient in producing magnetic fields well above equipartition and thus that it could represent a possible explanation for the large surface magnetic fields observed in neutron stars.

Binary black-hole systems with spins aligned with the orbital angular momentum are of special interest as they may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins aligned with the orbital angular momentum but otherwise arbitrary. By least-square fitting the results of these simulations we have constructed two "spin diagrams" which provide straightforward information about the recoil velocity |v_kick| and the final black-hole spin a_fin in terms of the dimensionless spins a_1 and a_2 of the two initial black holes. Overall they suggest a maximum recoil velocity of |v_kick|=441.94 km/s, and minimum and maximum final spins a_fin=0.3471 and a_fin=0.9591, respectively.

Determining the final spin of a black-hole (BH) binary is a question of key importance in astrophysics. Modelling this quantity in general is made difficult by the fact that it depends on the 7-dimensional space of parameters characterizing the two initial black holes. However, in special cases, when symmetries can be exploited, the description can become simpler. For black-hole binaries with unequal masses but with equal spins which are aligned with the orbital angular momentum, we show that the use of recent simulations and basic but exact constraints derived from the extreme mass-ratio limit allow to model this quantity with a simple analytic expression. Despite the simple dependence, the expression models very accurately all of the available estimates, with errors of a couple of percent at most. We also discuss how to use the fit to predict when a Schwarzschild BH is produced by the merger of two spinning BHs, when the total angular momentum of the spacetime ``flips'' sign, or under what conditions the final BH is ``spun-up'' by the merger. Finally, suggest an extension of the fit to include unequal-spin binaries, thus potentially providing a complete description of the final spin from the coalescence of generic black-hole binaries with spins aligned to the orbital angular momentum.

We have performed a detailed analysis of orbital motion in the vicinity of a nearly extremal Kerr black hole. For very rapidly rotating black holes (spin a=J/M>0.9524M) we have found a class of very strong field eccentric orbits whose angular momentum L_z increases with the orbit's inclination with respect to the equatorial plane, while keeping latus rectum and eccentricity fixed. This behavior is in contrast with Newtonian intuition, and is in fact opposite to the "normal" behavior of black hole orbits. Such behavior was noted previously for circular orbits; since it only applies to orbits very close to the black hole, they were named "nearly horizon-skimming orbits". Our analysis generalizes this result, mapping out the full generic (inclined and eccentric) family of nearly horizon-skimming orbits. The earlier work on circular orbits reported that, under gravitational radiation emission, nearly horizon-skimming orbits tend to evolve to smaller orbit inclination, toward prograde equatorial configuration. Normal orbits, by contrast, always demonstrate slowly growing orbit inclination (orbits evolve toward the retrograde equatorial configuration). Using up-to-date Teukolsky-fluxes, we have concluded that the earlier result was incorrect: all circular orbits, including nearly horizon-skimming ones, exhibit growing orbit inclination. Using kludge fluxes based on a Post-Newtonian expansion corrected with fits to circular and to equatorial Teukolsky-fluxes, we argue that the inclination grows also for eccentric nearly horizon-skimming orbits. We also find that the inclination change is, in any case, very small. As such, we conclude that these orbits are not likely to have a clear and peculiar imprint on the gravitational waveforms expected to be measured by the space-based detector LISA.

Binary black-hole systems with spins aligned with the orbital angular momentum are of special interest, as studies indicate that this configuration is preferred in nature. If the spins of the two bodies differ, there can be a prominent beaming of the gravitational radiation during the late plunge, causing a recoil of the final merged black hole. We perform an accurate and systematic study of recoil velocities from a sequence of equal-mass black holes whose spins are aligned with the orbital angular momentum, and whose individual spins range from a = +0.584 to -0.584. In this way we extend and refine the results of a previous study and arrive at a consistent maximum recoil of 448 +- 5 km/s for anti-aligned models as well as to a phenomenological expression for the recoil velocity as a function of spin ratio. This relation highlights a nonlinear behavior, not predicted by the PN estimates, and can be readily employed in astrophysical studies on the evolution of binary black holes in massive galaxies. An essential result of our analysis is the identification of different stages in the waveform, including a transient due to lack of an initial linear momentum in the initial data. Furthermore we are able to identify a pair of terms which are largely responsible for the kick, indicating that an accurate computation can be obtained from modes up to l=3. Finally, we provide accurate measures of the radiated energy and angular momentum, finding these to increase linearly with the spin ratio, and derive simple expressions for the final spin and the radiated angular momentum which can be easily implemented in N-body simulations of compact stellar systems. Our code is calibrated with strict convergence tests and we verify the correctness of our measurements by using multiple independent methods whenever possible.

Recent progress in numerical relativity has enabled us to model the non-perturbative merger phase of the binary black-hole coalescence problem. Based on these results, we propose a phenomenological family of waveforms which can model the inspiral, merger, and ring-down stages of black hole coalescence. We also construct a template bank using this family of waveforms and discuss its implementation in the search for signatures of gravitational waves produced by black-hole coalescences in the data of ground-based interferometers. This template bank might enable us to extend the present inspiral searches to higher-mass binary black-hole systems, i.e., systems with total mass greater than about 80 solar masses, thereby increasing the reach of the current generation of ground-based detectors.

Gravastars have been recently proposed as potential alternatives to explain the astrophysical phenomenology traditionally associated to black holes, raising the question of whether the two objects can be distinguished at all. Leaving aside the debate about the processes that would lead to the formation of a gravastar and the astronomical evidence in their support, we here address two basic questions: Is a gravastar stable against generic perturbations? If stable, can an observer distinguish it from a black hole of the same mass? To answer these questions we construct a general class of gravastars and determine the conditions they must satisfy in order to exist as equilibrium solutions of the Einstein equations. For such models we perform a systematic stability analysis against axial-perturbations, computing the real and imaginary parts of the eigenfrequencies. Overall, we find that gravastars are stable to axial perturbations, but also that their quasi-normal modes differ from those of a black hole of the same mass and thus can be used to discern, beyond dispute, a gravastar from a black hole.

The accurate modelling of astrophysical scenarios involving compact objects and magnetic fields, such as the collapse of rotating magnetized stars to black holes or the phenomenology of gamma-ray bursts, requires the solution of the Einstein equations together with those of general-relativistic magnetohydrodynamics. We present a new numerical code developed to solve the full set of general-relativistic magnetohydrodynamics equations in a dynamical and arbitrary spacetime with high-resolution shock-capturing techniques on domains with adaptive mesh refinements. After a discussion of the equations solved and of the techniques employed, we present a series of testbeds carried out to validate the code and assess its accuracy. Such tests range from the solution of relativistic Riemann problems in flat spacetime, over to the stationary accretion onto a Schwarzschild black hole and up to the evolution of oscillating magnetized stars in equilibrium and constructed as consistent solutions of the coupled Einstein-Maxwell equations.

We present new results on dynamical instabilities in rapidly rotating neutron-stars. In particular, using numerical simulations in full General Relativity, we analyse the effects that the stellar compactness has on the threshold for the onset of the dynamical bar-mode instability, as well as on the appearance of other dynamical instabilities. By using an extrapolation technique developed and tested in our previous study [1], we explicitly determine the threshold for a wide range of compactnesses using four sequences of models of constant baryonic mass comprising a total of 59 stellar models. Our calculation of the threshold is in good agreement with the Newtonian prediction and improves the previous post-Newtonian estimates. In addition, we find that for stars with sufficiently large mass and compactness, the m=3 deformation is the fastest growing one. For all of the models considered, the non-axisymmetric instability is suppressed on a dynamical timescale with an m=1 deformation dominating the final stages of the instability. These results, together with those presented in [1], suggest that an m=1 deformation represents a general and late-time feature of non-axisymmetric dynamical instabilities both in full General Relativity and in Newtonian gravity.

We present accurate simulations of the dynamical bar-mode instability in full General Relativity focussing on two aspects which have not been investigated in detail in the past. Namely, on the persistence of the bar deformation once the instability has reached its saturation and on the precise determination of the threshold for the onset of the instability in terms of the parameter $\beta={T}/{|W|}$. We find that generic nonlinear mode-coupling effects appear during the development of the instability and these can severely limit the persistence of the bar deformation and eventually suppress the instability. In addition, we observe the dynamics of the instability to be strongly influenced by the value $\beta$ and on its separation from the critical value $\beta_c$ marking the onset of the instability. We discuss the impact these results have on the detection of gravitational waves from this process and provide evidence that the classical perturbative analysis of the bar-mode instability for Newtonian and incompressible Maclaurin spheroids remains qualitatively valid and accurate also in full General Relativity.

To investigate the imprint on the gravitational-wave emission from extreme mass-ratio inspirals in non-pure Kerr spacetimes, we have studied the ``kludge'' waveforms generated in highly-accurate, numerically-generated spacetimes containing a black hole and a self-gravitating, homogeneous torus with comparable mass and spin. In order to maximize their impact on the produced waveforms, we have considered tori that are compact, massive and close to the central black hole, investigating under what conditions the LISA experiment could detect their presence. Our results show that for a large portion of the space of parameters the waveforms produced by EMRIs in these black hole-torus systems are indistinguishable from pure-Kerr waveforms. Hence, a ``confusion problem'' will be present for observations carried out over a timescale below or comparable to the dephasing time.

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.

A paradigm deeply rooted in modern numerical relativity calculations prescribes the removal of those regions of the computational domain where a physical singularity may develop. We here challenge this paradigm by performing three-dimensional simulations of the collapse of uniformly rotating stars to black holes without excision. We show that this choice, combined with suitable gauge conditions and the use of minute numerical dissipation, improves dramatically the long-term stability of the evolutions. In turn, this allows for the calculation of the waveforms well beyond what previously possible, providing information on the black-hole ringing and setting a new mark on the present knowledge of the gravitational-wave emission from the stellar collapse to a rotating black hole.

We discuss the procedure for the exact solution of the Riemann problem in special relativistic magnetohydrodynamics (MHD). We consider both initial states leading to a set of only three waves analogous to the ones in relativistic hydrodynamics, as well as generic initial states leading to the full set of seven MHD waves. Because of its generality, the solution presented here could serve as an important test for those numerical codes solving the MHD equations in relativistic regimes.

We present a new method for the calculation of black hole perturbations induced by extended sources in which the solution of the nonlinear hydrodynamics equations is coupled to a perturbative method based on Regge-Wheeler/Zerilli and Bardeen-Press-Teukolsky equations when these are solved in the frequency domain. In contrast to alternative methods in the time domain which may be unstable for rotating black-hole spacetimes, this approach is expected to be stable as long as an accurate evolution of the matter sources is possible. Hence, it could be used under generic conditions and also with sources coming from three-dimensional numerical relativity codes. As an application of this method we compute the gravitational radiation from an oscillating high-density torus orbiting around a Schwarzschild black hole and show that our method is remarkably accurate, capturing both the basic quadrupolar emission of the torus and the excited emission of the black hole.

Velocities close to the speed of light are a robust observational property of the jets observed in microquasars and AGNs, and are expected to be behind much of the phenomenology of GRBs. Yet, the mechanism boosting relativistic jets to such large Lorentz factors is still essentially unknown. Building on recent general-relativistic, multidimensional simulations of progenitors of short GRBs, we discuss a new effect in relativistic hydrodynamics which can act as an efficient booster in jets. This effect is purely hydrodynamical and occurs when large velocities tangential to a discontinuity are present in the flow, yielding Lorentz factors $\Gamma \sim 10^2-10^3$ or larger in flows with moderate initial Lorentz factors. Although without a Newtonian counterpart, this effect can be explained easily through the most elementary hydrodynamical flow: i.e., a relativistic Riemann problem.

We use a relativistic ray-tracing code to analyze the X-ray emission from a pressure-supported oscillating relativistic torus around a black hole. We show that a strong correlation exists between the {\it intrinsic} frequencies of the torus normal modes and the {\it extrinsic} frequencies seen in the observed light curve power spectrum. This correlation demonstrates the feasibility of the oscillating-torus model to explain the multiple peaks seen in black hole high-frequency quasi-periodic oscillations. Using an optically thin, monochromatic emission model, we also determine how a relativistically broadened emission line and the amplitude of the X-ray modulations are dependent on the observer's inclination angle and on the torus oscillation amplitudes. Observations of these features can provide important information about the torus as well as the black hole.

We present the first calculation of gravitational wave emission produced in the gravitational collapse of uniformly rotating neutron stars to black holes in fully three-dimensional simulations. The initial stellar models are relativistic polytropes which are dynamically unstable and with angular velocities ranging from slow rotation to the mass-shedding limit. An essential aspect of these simulations is the use of progressive mesh-refinement techniques which allow to move the outer boundaries of the computational domain to regions where gravitational radiation attains its asymptotic form. The waveforms have been extracted using a gauge-invariant approach in which the numerical spacetime is matched with the non-spherical perturbations of a Schwarzschild spacetime. Overall, the results indicate that the waveforms have features related to the properties of the initial stellar models (in terms of their w-mode oscillations) and of the newly produced rotating black holes (in terms of their quasi-normal modes). While our waveforms are in good qualitative agreement with those computed by Stark and Piran in two-dimensional simulations, our amplitudes are about one order of magnitude smaller and this difference is mostly likely due to our less severe pressure reduction. For a neutron star rotating uniformly near mass-shedding and collapsing at 10 kpc, the signal-to-noise ratio computed uniquely from the burst is S/N ~ 0.25, but this grows to be S/N <~ 4 in the case of LIGO II.

Results are presented from general relativistic numerical computations of primordial black-hole formation during the radiation-dominated era of the universe. Growing-mode perturbations are specified within the linear regime and their subsequent evolution is followed as they become nonlinear. We use a spherically symmetric Lagrangian code and study both super-critical perturbations, which go on to produce black holes, and sub-critical perturbations, for which the overdensity eventually disperses into the background medium. For super-critical perturbations, we confirm the results of previous work concerning scaling-laws but note that the threshold amplitude for a perturbation to lead to black-hole formation is substantially reduced when the initial conditions are taken to represent purely growing modes. For sub-critical cases, where an initial collapse is followed by a subsequent re-expansion, strong compressions and rarefactions are seen for perturbation amplitudes near to the threshold. We have also investigated the effect of including a significant component of vacuum energy and have calculated the resulting changes in the threshold and in the slope of the scaling law.

We present a new three-dimensional fully general-relativistic hydrodynamics code using high-resolution shock-capturing techniques and a conformal traceless formulation of the Einstein equations. Besides presenting a thorough set of tests which the code passes with very high accuracy, we discuss its application to the study of the gravitational collapse of uniformly rotating neutron stars to Kerr black holes. The initial stellar models are modelled as relativistic polytropes which are either secularly or dynamically unstable and with angular velocities which range from slow rotation to the mass-shedding limit. We investigate the gravitational collapse by carefully studying not only the dynamics of the matter, but also that of the trapped surfaces, i.e. of both the apparent and event horizons formed during the collapse. The use of these surfaces, together with the dynamical horizon framework, allows for a precise measurement of the black-hole mass and spin. The ability to successfully perform these simulations for sufficiently long times relies on excising a region of the computational domain which includes the singularity and is within the apparent horizon. The dynamics of the collapsing matter is strongly influenced by the initial amount of angular momentum in the progenitor star and, for initial models with sufficiently high angular velocities, the collapse can lead to the formation of an unstable disc in differential rotation. All the simulations performed with uniformly rotating initial data and a polytropic or ideal-fluid equation of state show no evidence of shocks or of the presence of matter on stable orbits outside the black hole.

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.

This is the second of a series of papers investigating the oscillation properties of relativistic, non-selfgravitating tori orbiting around black holes. Extending the work done in a Schwarzschild background, we here consider the axisymmetric oscillations of vertically integrated tori in a Kerr spacetime. The tori are modeled with a number of different non-Keplerian distributions of specific angular momentum and we discuss how the oscillation properties depend on these and on the rotation of the central black hole. We first consider a local analysis to highlight the relations between acoustic and epicyclic oscillations in a Kerr spacetime and subsequently perform a global eigenmode analysis to compute the axisymmetric p modes. In analogy with what found in a Schwarzschild background, these modes behave as sound waves that are modified by rotation and are globally trapped in the torus. For constant distributions of specific angular momentum, the eigenfrequencies appear in a sequence 2:3:4:... which is essentially independent of the size of the disc and of the black hole rotation. For non-constant distributions of angular momentum, on the other hand, the sequence depends on the properties of the disc and on the spin of the black hole, becoming harmonic for sufficiently large tori. We also comment on how p modes could explain the high frequency quasi-periodic oscillations observed in low-mass X-ray binaries with a black hole candidate and the properties of an equivalent model in Newtonian physics.

Relativistic stars are endowed with intense electromagnetic fields but are also subject to oscillations of various types. We here investigate the impact that oscillations have on the electric and magnetic fields external to a relativistic star in vacuum. In particular, modelling the star as a relativistic polytrope with infinite conductivity, we consider the solution of the general relativistic Maxwell equations both in the vicinity of the stellar surface and far from it, once a perturbative velocity field is specified. In this first paper we present general analytic expressions that are not specialized to a particular magnetic field topology or velocity field. However, as a validating example and an astrophysically important application, we consider a dipolar magnetic field and the velocity field corresponding to the rotation of the misaligned dipole. Besides providing analytic expressions for the electromagnetic fields produced by this configuration, we calculate, for the first time, the general relativistic energy loss through dipolar electromagnetic radiation. We find that the widely used Newtonian expression under-estimates this loss by a factor of 2-6 depending on the stellar compactness. This correction could have important consequences in the study of the spin evolution of pulsars.

We present analytic solutions of Maxwell equations in the internal and external background spacetime of a slowly rotating magnetized neutron star. The star is considered isolated and in vacuum, with a dipolar magnetic field not aligned with the axis of rotation. With respect to a flat spacetime solution, general relativity introduces corrections related both to the monopolar and the dipolar parts of the gravitational field. In particular, we show that in the case of infinite electrical conductivity general relativistic corrections due to the dragging of reference frames are present, but only in the expression for the electric field. In the case of finite electrical conductivity, however, corrections due both to the spacetime curvature and to the dragging of reference frames are shown to be present in the induction equation. These corrections could be relevant for the evolution of the magnetic fields of pulsars and magnetars. The solutions found, while obtained through some simplifying assumption, reflect a rather general physical configuration and could therefore be used in a variety of astrophysical situations.

Observations of X-ray emissions from binary systems have always been considered important tools to test the validity of General Relativity in strong-field regimes. The pairs and triplets of high frequency quasi-periodic oscillations 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 over the years to explain the QPOs 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 QPOs result from basic p-mode oscillations of a non-Keplerian disc of finite size. We show that within this new model all of the key properties of the QPOs: a) harmonic ratios of frequencies even as the frequencies change; b) variations in the relative strength of the frequencies with spectral energy distribution and with photon energy; c) small and systematic changes in the frequencies, can all be explained simply given a single reasonable assumption.

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.

We have investigated the turbulent mean-field dynamo action in protoneutron stars that are subject to convective and neutron finger instabilities. While the first one develops mostly in the inner regions of the star, the second one is favoured in the outer regions, where the Rossby number is much smaller and a mean-field dynamo action is more efficient. By solving the mean-field induction equation we have computed the critical spin period below which no dynamo action is possible and found it to be ~1 s for a wide range of stellar models. Because this critical period is substantially longer than the characteristic spin period of very young pulsars, we expect that a mean-field dynamo will be effective for most protoneutron stars.

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.

I present results of two-dimensional general relativistic hydrodynamical simulations of constant specific angular momentum tori orbiting around a Schwarzschild black hole. After introducing axisymmetric perturbations, these objects either become unstable to the runaway instability or respond with regular oscillations. The latter, in particular, are responsible for quasi-periodic bursts of accretion onto the black hole as well as for the emission of intense gravitational radiation, with signal-to-noise ratios at the detector which are comparable or even larger than the typical ones expected in stellar-core collapse.

These lectures aim at providing an introduction to the properties of gravitational waves and in particular to those gravitational waves that are expected as a consequence of perturbations of black holes and neutron stars. Imprinted in the gravitational radiation emitted by these objects is, in fact, a wealth of physical information. In the case of black holes, a detailed knowledge of the gravitational radiation emitted as a response to perturbations will reveal us important details about their mass and spin, but also about the fundamental properties of the event horizon. In the case of neutron stars, on the other hand, this information can provide a detailed map of their internal structure and tell us about the equation of state of matter at very high density, thus filling-in a gap in energies and densities that cannot be investigated by experiments in terrestrial laboratories.

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 computed the eigenfrequencies of f-modes for a constant-rest-mass sequences of rapidly rotating relativistic inviscid stars in differential rotation. The frequencies have been calculated neglecting the metric perturbations (the relativistic Cowling approximation) and expressed as a function of the ratio between the rotational kinetic energy and the absolute value of the gravitational energy of the stellar model beta=T/|W|. The zeros and the end-points of these sequences mark respectively the onset of the secular instability driven by gravitational radiation-reaction and the maximum value of beta at which an equilibrium model exists. In differentially rotating stars the secular stability limits appear at a beta larger than those found for uniformly rotating stars. Differential rotation, on the other hand, also allows for the existence of equilibrium models at values of beta larger than those for uniformly rotating stars, moving the end-point of the sequences to larger beta. As a result, for some degrees of differential rotation, the onset of the secular instability for f-modes is generally favoured by the presence of differential rotation.

We show that, at first order in the angular velocity, the general relativistic description of Rossby-Haurwitz waves (the analogues of r-waves on a thin shell) can be obtained from the corresponding Newtonian one after a coordinate transformation. As an application, we show that the results recently obtained by Rezzolla and Yoshida (2001) in the analysis of Newtonian Rossby-Haurwitz waves of a slowly and differentially rotating, fluid shell apply also in General Relativity, at first order in the angular velocity.

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.

This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the application of this code to problems in general relativistic astrophysics. In particular, we report on the accuracy of our code in the long-term dynamical evolution of relativistic stars and on some new physics results obtained in the process of code testing. The tests involve single non-rotating stars in stable equilibrium, non-rotating stars undergoing radial and quadrupolar oscillations, non-rotating stars on the unstable branch of the equilibrium configurations migrating to the stable branch, non-rotating stars undergoing gravitational collapse to a black hole, and rapidly rotating stars in stable equilibrium and undergoing quasi-radial oscillations. The numerical evolutions have been carried out in full general relativity using different types of polytropic equations of state using either the rest-mass density only, or the rest-mass density and the internal energy as independent variables. New variants of the spacetime evolution and new high resolution shock capturing (HRSC) treatments based on Riemann solvers and slope limiters have been implemented and the results compared with those obtained from previous methods. Finally, we have obtained the first eigenfrequencies of rotating stars in full general relativity and rapid rotation. A long standing problem, such frequencies have not been obtained by other methods. Overall, and to the best of our knowledge, the results presented in this paper represent the most accurate long-term three-dimensional evolutions of relativistic stars available to date.

Following the general formalism presented by Rezzolla, Ahmedov and Miller (MNRAS, 322, 723 2001), we here derive analytic solutions of the electromagnetic fields equations in the internal and external background spacetime of a slowly rotating highly conducting magnetized neutron star. The star is assumed to be isolated and in vacuum, with a dipolar magnetic field not aligned with the axis of rotation. Our results indicate that the electromagnetic fields of a slowly rotating neutron star are modified by general relativistic effects arising from both the monopolar and the dipolar parts of the gravitational field. The results presented here differ from the ones discussed by Rezzolla, Ahmedov and Miller (MNRAS, 322, 723 2001) mainly in that we here consider the interior magnetic field to be dipolar with the same radial dependence as the external one. While this assumption might not be a realistic one, it should be seen as the application of our formalism to a case often discussed in the literature.

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.

The evolution of the r-mode instability is likely to be accompanied by secular kinematic effects which will produce differential rotation with large scale drifts of fluid elements, mostly in the azimuthal direction. As first discussed by Rezzolla, Lamb and Shapiro 2000, the interaction of these secular velocity fields with a pre-existing neutron star magnetic field could result in the generation of intense and large scale toroidal fields. Following their derivation in the companion paper, we here discuss the numerical solution of the evolution equations for the magnetic field. The values of the magnetic fields obtained in this way are used to estimate the conditions under which the r-mode instability might be prevented or suppressed. We also assess the impact of the generation of large magnetic fields on the gravitational wave detectability of r-mode unstable neutron stars. Our results indicate that the signal to noise ratio in the detection of gravitational waves from the r-mode instability might be considerably decreased if the latter develops in neutron stars with initial magnetic fields larger than 10^10 G.

The instability of r-mode oscillations in rapidly rotating neutron stars has attracted attention as a potential mechanism for producing high frequency, almost periodic gravitational waves. The analyses carried so far have shown the existence of these modes and have considered damping by shear and bulk viscosity. However, the magnetohydrodynamic coupling of the modes with a stellar magnetic field and its role in the damping of the instability has not been fully investigated yet. Following our introductory paper (Rezzolla, Lamb and Shapiro 2000), we here discuss in more detail the existence of secular higher-order kinematical effects which will produce toroidal fluid drifts. We also define the sets of equations that account for the time evolution of the magnetic fields produced by these secular velocity fields and show that the magnetic fields produced can reach equipartition in less than a year. The full numerical calculations as well as the evaluation of the impact of strong magnetic fields on the onset and evolution of the r-mode instability will be presented in a companion paper.

Recent studies have raised doubts about the occurrence of r modes in Newtonian stars with a large degree of differential rotation. To assess the validity of this conjecture we have solved the eigenvalue problem for Rossby-Haurwitz waves (the analogues of r waves on a thin-shell) in the presence of differential rotation. The results obtained indicate that the eigenvalue problem is never singular and that, at least for the case of a thin-shell, the analogues of r modes can be found for arbitrarily large degrees of differential rotation. This work clarifies the puzzling results obtained in calculations of differentially rotating axi-symmetric Newtonian stars.

Gravitational wave emission is considered to be the driving force for the evolution of short-period cataclysmic binary stars, making them a potential test for the validity of General Relativity. In spite of continuous refinements of the physical description, a 10% mismatch exists between the theoretical minimum period ($P_{\rm turn} \simeq 70$ min) and the short-period cut-off ($P_{\rm min} \simeq 80$ min) observed in the period distribution for cataclysmic variable binaries. A possible explanation for this mismatch was associated with the use of the Roche model. We here present a systematic comparison between self-consistent, numerically constructed sequences of hydrostatic models of binary stars and Roche models of semi-detached binaries. On the basis of our approach, we also derive a value for the minimum period of cataclysmic variable binaries. The results obtained through the comparison indicate that the Roche model is indeed very good, with deviations from the numerical solution which are of a few percent at most. Our results therefore suggest that additional sources of angular momentum loss or alternative explanations need to be considered in order to justify the mismatch.

We show that r-mode oscillations distort the magnetic fields of neutron stars and that their occurrence is likely to be limited by this interaction. If the field is gtrsim 10^{16} (Omega/Omega_B) G, where Omega and Omega_B are the angular velocities of the star and at which mass shedding occurs, r-mode oscillations cannot occur. Much weaker fields will prevent gravitational radiation from exciting r-mode oscillations or damp them on a relatively short timescale by extracting energy from the modes faster than gravitational wave emission can pump energy into them. For example, a 10^{10} G poloidal magnetic field that threads the star's superconducting core is likely to prevent the ell=2 mode from being excited unless Omega exceeds 0.35 Omega_B. If Omega is larger than 0.35 Omega_B initially, the ell=2 mode may be excited but is likely to decay rapidly once Omega falls below 0.35 Omega_B, which happens in lesssim 15^d if the saturation amplitude is gtrsim 0.1. The r-mode oscillations may play an important role in determining the structure of neutron star magnetic fields.

We report on our numerical implementation of fully relativistic hydrodynamics coupled to Einstein's field equations in three spatial dimensions. We briefly review several steps in our code development, including our recasting of Einstein's equations and several tests which demonstrate its advantages for numerical integrations. We outline our implementation of relativistic hydrodynamics, and present numerical results for the evolution of both stable and unstable Oppenheimer-Volkov equilibrium stars, which represent a very promising first test of our code.

We present a new set of 3.5 Post-Newtonian equations in which Newtonian hydrodynamics is coupled to the nonconservative effects of gravitational radiation emission. Our formalism differs in two significant ways from a similar 3.5 Post-Newtonian approach proposed by Blanchet (1993, 1997). Firstly we concentrate only on the radiation-reaction effects produced by a time-varying mass-current quadrupole $S_{ij}$. Secondly, we adopt a gauge in which the radiation-reaction force densities depend on the fourth time derivative of $S_{ij}$, rather than on the fifth, as in Blanchet's approach. This difference makes our formalism particularly well-suited to numerical implementation and could prove useful in performing fully numerical simulations of the recently discovered $r$-mode instability for rotating neutron stars subject to axial perturbations.

One of the most interesting questions regarding a possible first order cosmological quark--hadron phase transition concerns the final fate of the baryon number contained within the disconnected quark regions at the end of the transition. We here present a detailed investigation of the hydrodynamical evolution of an evaporating quark drop, using a multi-component fluid description to follow the mechanisms of baryon number segregation. With this approach, we are able to take account of the simultaneous effects of baryon number flux suppression at the phase interface, entropy extraction by means of particles having long mean-free-paths, and baryon number diffusion. A range of computations has been performed to investigate the permitted parameter-space and this has shown that significant baryon number concentrations, perhaps even up to densities above that of nuclear matter, represent an inevitable outcome within this scenario.

We present results of a linear stability analysis of relativistic detonation fronts, which have been considered as representing phase interfaces in cosmological first order phase transitions. After discussing general stability conditions for detonation fronts, we concentrate on the properties of the fronts with respect to corrugation instabilities and discuss separately the cases of Chapman-Jouguet and strong detonation waves. Contrarily to what recently claimed, we find that strong detonations are both evolutionary and stable with respect to corrugations of the front. Moreover, Chapman-Jouguet detonations appear to be unconditionally linearly stable. The implications of the stability results for first order cosmological phase transitions are presented and a discussion of the causal structure of reaction fronts is also given.

We discuss the results of a full relativistic treatment of the hydrodynamics of disconnected quark regions during the final stages of the cosmological quark-hadron transition. In this study, which represents a further development of a previous analysis of the evaporation of cosmological quark drops, the effects of long range energy and momentum transfer via electromagnetically interacting particles are consistently taken into account. For this purpose, a set of relativistic Lagrangian equations describing the evolution of the strongly interacting fluids is coupled to a system of equivalent equations accounting for the hydrodynamics of the fluid of electromagnetically interacting particles. The complete set of equations has then been solved numerically and results are presented from this. The inclusion of relativistic radiative transfer produces significantly different results, with the formation of high density regions at the end of the drop evaporation being particularly relevant. A comparison is made with results obtained for the previous radiation-free model and the cosmological implications concerning baryon number concentrations are briefly discussed.

We have carried out a study of the hydrodynamics of disconnected quark regions during the final stages of the cosmological quark-hadron transition. A set of relativistic Lagrangian equations is presented for following the evaporation of a single quark drop and results from the numerical solution of this are discussed. A self-similar solution is shown to exist and the formation of baryon number density inhomogeneities at the end of the drop contraction is discussed.

Two previous papers in this series have presented a study of the growth of hadronic bubbles during the cosmological Quark--Hadron transition, treating the material within each phase as a single perfect fluid. Here, we extend the analysis to include the effects of long-range energy and momentum transfer by weakly and electromagnetically interacting particles. After a short review of the formalism adopted, we discuss the numerical strategies used in the computer code which has been constructed in order to solve this system of equations. Results for the growth of single hadronic bubbles are also presented.

We present a new complete set of Lagrangian relativistic hydrodynamical equations describing the transfer of energy and momentum between a standard fluid and a radiation fluid in a general non-stationary spherical flow. The new set of equations has been derived for a particular application to the study of the cosmological Quark--Hadron transition but can also be used in other contexts.