BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Luciano Rezzolla appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
In this short note, we briefly comment on the analytical bounds that must be imposed on the parameter space of the Rezzolla-Zhidenko (RZ) metric-parametrization approach introduced in Ref. [1]. We hope this will clarify some of the confusion recently emerged on this issue [2].
Imaging a supermassive black hole and extracting physical information requires good knowledge of both the gravitational and the astrophysical conditions near the black hole. When the geometrical properties of the black hole are well understood, extracting information on the emission properties is possible. Similarly, when the emission properties are well understood, extracting information on the black-hole geometry is possible. At present however, uncertainties are present both in the geometry and in the emission, and this inevitably leads to degeneracies in the interpretation of the observations. We explore here the impact of varying geometry and emission coefficient when modelling the imaging of a spherically-accreting black hole. Adopting the Rezzolla-Zhidenko parametric metric to model arbitrary static black-holes, we first demonstrate how shadow-size measurements leave degeneracies in the multidimensional space of metric-deviation parameters, even in the limit of infinite-precision measurements. Then, at finite precision, we show that these degenerate regions can be constrained when multiple pieces of information, such as the shadow-size and the peak image intensity contrast, are combined. Such degeneracies can potentially be eliminated with measurements at increased angular-resolution and flux-sensitivity. While our approach is restricted to spherical symmetry and hence idealised, we expect our results to hold also when more complex geometries and emission processes are considered.
Determining the phase structure of Quantum Chromodynamics (QCD) and its Equation of State (EOS) at densities and temperatures realized inside neutron stars and their mergers is a long-standing open problem. The holographic V-QCD framework provides a model for the EOS of dense and hot QCD, which describes the deconfinement phase transition between a dense baryonic and a quark matter phase. We use this model in fully general relativistic hydrodynamic (GRHD) simulations to study the formation of quark matter and the emitted gravitational wave signal of binary systems that are similar to the first ever observed neutron star merger event GW170817.
The gravitational-wave events GW170817 and GW190425 have led to a number of important insights on the equation of state of dense matter and the properties of neutron stars, such as their radii and the maximum mass. Some of these conclusions have been drawn on the basis of numerical-relativity simulations of binary neutron-star mergers with vanishing initial spins. While this may be a reasonable assumption in equal-mass systems, it may be violated in the presence of large mass asymmetries accompanied by the presence of high spins. To quantify the impact of high spins on multi-messenger gravitational-wave events, we have carried out a series of high-mass binary neutron-star mergers with a highly spinning primary star and large mass asymmetries that have been modelled self-consistently using two temperature-dependent equations of state. We show that, when compared with equal-mass, irrotational binaries, these systems can lead to significant differences in the remnant lifetime, in the dynamical ejecta, in the remnant disc masses, in the secular ejecta, and on the bulk kilonova properties. These differences could be exploited to remove the degeneracy between low- and high-spin priors in the detection of gravitational waves from binary neutron-star mergers.
The maximum mass of a nonrotating neutron star, $M_{\rm TOV}$, plays a very important role in deciphering the structure and composition of neutron stars and in revealing the equation of state (EOS) of nuclear matter. Although with a large-error bar, the recent mass estimate for the black-widow binary pulsar PSR J0952-0607, i.e. $M=2.35\pm0.17~M_\odot$, provides the strongest lower bound on $M_{\rm TOV}$ and suggests that neutron stars with very large masses can in principle be observed. Adopting an agnostic modelling of the EOS, we study the impact that large masses have on the neutron-star properties. In particular, we show that assuming $M_{\rm TOV}\gtrsim 2.35\,M_\odot$ constrains tightly the behaviour of the pressure as a function of the energy density and moves the lower bounds for the stellar radii to values that are significantly larger than those constrained by the NICER measurements, rendering the latter ineffective in constraining the EOS. We also provide updated analytic expressions for the lower bound on the binary tidal deformability in terms of the chirp mass and show how larger bounds on $M_{\rm TOV}$ lead to tighter constraints for this quantity. In addition, we point out a novel quasi-universal relation for the pressure profile inside neutron stars that is only weakly dependent from the EOS and the maximum-mass constraint. Finally, we study how the sound speed and the conformal anomaly are distributed inside neutron stars and show how these quantities depend on the imposed maximum-mass constraints.
Magnetic reconnection is thought to be one of the main energy-dissipation mechanisms fueling energy to the plasma in the vicinity of a black hole. Indeed, plasmoids formed through magnetic reconnection may play a key role in $γ$-ray, X-ray and near-infrared flares from the black hole at the center of our galaxy, SgrA*. We report the results of three-dimensional general-relativistic ideal and resistive magnetohydrodynamics simulations modelling magnetic reconnection in accretion flows around astrophysical black holes. As an important difference with similar works, our accretion discs have an initial dipolar magnetic-field configuration with loops of alternating polarity. We show that current sheets are formed and destroyed rapidly in the turbulent environment of black-hole accretion. Plasmoids are formed from current sheets close to the event horizon, in a region of $\sim2-15$ gravitational radii. We further quantify the magnetic dissipation and the process of energy transfer to the plasmoids, reporting the reconnection rate, the relative current density with respect to the local magnetic field, and the size of the plasmoids. We find that plasmoids gain energy through reconnection and heat up to relativistic temperatures, with the largest ones being sufficiently energetic to leave the black hole near the polar regions. During their evolution, plasmoids are stretched and elongated, becoming disrupted when the shear is sufficiently large, although some plasmoids survive as well-distinguished structures at distances of $\sim30-40$ gravitational radii from the black hole. Finally, we find that in some cases the plasmoids acquire a super-Keplerian azimuthal velocity, as suggested by recent observations of flares from Sgr~A*.
We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.
Ultra-high angular resolution in astronomy has always been an important vehicle for making fundamental discoveries. Recent results in direct imaging of the vicinity of the supermassive black hole in the nucleus of the radio galaxy M87 by the millimeter VLBI system Event Horizon Telescope and various pioneering results of the Space VLBI mission RadioAstron provided new momentum in high angular resolution astrophysics. In both mentioned cases, the angular resolution reached the values of about 10-20 microrcseconds. Further developments toward at least an order of magnitude "sharper" values are dictated by the needs of astrophysical studies and can only be achieved by placing millimeter and submillimeter wavelength interferometric systems in space. A concept of such the system, called Terahertz Exploration and Zooming-in for Astrophysics (THEZA), has been proposed in the framework of the ESA Call for White Papers for the Voayage 2050 long term plan in 2019. In the current paper we discuss several approaches for addressing technological challenges of the THEZA concept. In particular, we consider a novel configuration of a space-borne millimeter/sub-millimeter antenna which might resolve several bottlenecks in creating large precise mechanical structures. The paper also presents an overview of prospective space-qualified technologies of low-noise analogue front-end instrumentation for millimeter/sub-millimeter telescopes, data handling and processing. The paper briefly discusses approaches to the interferometric baseline state vector determination and synchronisation and heterodyning system. In combination with the original ESA Voyage 2050 White Paper, the current work sharpens the case for the next generation microarcsceond-level imaging instruments and provides starting points for further in-depth technology trade-off studies.
We perform three-dimensional (3D) general-relativistic magnetohydrodynamic simulations to model the jet break-out from the ejecta expected to be produced in a binary neutron-star merger. The structure of the relativistic outflow from the 3D simulation confirms our previous results from 2D simulations, namely, that a relativistic magnetized outflow breaking out from the merger ejecta exhibits a hollow core of $θ_{\rm core}\approx4^{\circ}$, an opening angle of $θ_{\rm jet}\gtrsim10^{\circ}$, and is accompanied by a wind of ejected matter that will contribute to the kilonova emission. We also compute the non-thermal afterglow emission of the relativistic outflow and fit it to the panchromatic afterglow from GRB170817A, together with the superluminal motion reported from VLBI observations. In this way, we deduce an observer angle of $θ_{\rm obs}= 35.7^{\circ \,\,+1.8}_{\phantom{\circ \,\,}-2.2}$. We further compute the afterglow emission from the ejected matter and constrain the parameter space for a scenario in which the matter responsible for the thermal kilonova emission will also lead to a non-thermal emission yet to be observed.
The detection of the binary events GW170817 and GW190814 has provided invaluable constraints on the maximum mass of nonrotating configurations of neutron stars, $M_{_{\rm TOV}}$. However, the large differences in the neutron-star masses measured in GW170817 and GW190814 has also lead to a significant tension between the predictions for such maximum masses, with GW170817 suggesting that $M_{_{\rm TOV}} \lesssim 2.3\,M_{\odot}$, and GW190814 requiring $M_{_{\rm TOV}} \gtrsim 2.5\,M_{\odot}$ if the secondary was a (non- or slowly rotating) neutron star at merger. Using a genetic algorithm, we sample the multidimensional space of parameters spanned by gravitational-wave and astronomical observations associated with GW170817. Consistent with previous estimates, we find that all of the physical quantities are in agreement with the observations if the maximum mass is in the range $M_{_{\rm TOV}} = 2.210^{+0.116}_{-0.123} \,M_{\odot}$ within a $2\textrm{-}σ$ confidence level. By contrast, maximum masses with $M_{_{\rm TOV}} \gtrsim 2.5\,M_{\odot}$, not only require efficiencies in the gravitational-wave emission that are well above the numerical-relativity estimates, but they also lead to a significant under-production of the ejected mass. Hence, the tension can be released by assuming that the secondary in GW190814 was a black hole at merger, although it could have been a rotating neutron star before.
The recent 230 GHz observations of the Event Horizon Telescope (EHT) are able to image the innermost structure of the M87 and show a ring-like structure which is in agreement with thermal synchrotron emission generated in a torus surrounding a supermassive black hole. However, at lower frequencies M87 is characterised by a large-scale and edge-brightened jet with clear signatures of non-thermal emission. In order to bridge the gap between these scales and to provide a theoretical interpretation of these observations we perform general relativistic magnetohydrodynamic simulations of accretion on to black holes and jet launching. M87 has been the target for multiple observations across the entire electromagnetic spectrum. Among these VLBI observations provide unique details on the collimation profile of the jet down to several gravitational radii. In this work we aim to model the observed broad-band spectrum of M87 from the radio to the NIR regime and at the same time fit the jet structure as observed with Global mm-VLBI at 86 GHz. We use general relativistic magnetohydrodynamics and simulate the accretion of the magnetised plasma onto Kerr-black holes in 3D. The radiative signatures of these simulations are computed taking different electron distribution functions into account and a detailed parameter survey is performed in order to match the observations. The results of our simulations show that magnetically arrested disks around fast spinning black holes ($a_\star\geq0.5$) together with a mixture of thermal and non-thermal particle distributions are able to model simultaneously the broad-band spectrum and the innermost jet structure of M87
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.
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...
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.
ExaHyPE ("An Exascale Hyperbolic PDE Engine") is a software engine for solving systems of first-order hyperbolic partial differential equations (PDEs). Hyperbolic PDEs are typically derived from the conservation laws of physics and are useful in a wide range of application areas. Applications powered by ExaHyPE can be run on a student's laptop, but are also able to exploit thousands of processor cores on state-of-the-art supercomputers. The engine is able to dynamically increase the accuracy of the simulation using adaptive mesh refinement where required. Due to the robustness and shock capturing abilities of ExaHyPE's numerical methods, users of the engine can simulate linear and non-linear hyperbolic PDEs with very high accuracy. Users can tailor the engine to their particular PDE by specifying evolved quantities, fluxes, and source terms. A complete simulation code for a new hyperbolic PDE can often be realised within a few hours - a task that, traditionally, can take weeks, months, often years for researchers starting from scratch. In this paper, we showcase ExaHyPE's workflow and capabilities through real-world scenarios from our two main application areas: seismology and astrophysics.
The theory of gauge-invariant non-spherical metric perturbations of Schwarzschild black hole spacetimes is now well established. Yet, as different notations and conventions have been used throughout the years, the literature on the subject is often confusing and sometimes confused. The purpose of this paper is to review and collect the relevant expressions related to the Regge-Wheeler and Zerilli equations for the odd and even-parity perturbations of a Schwarzschild spacetime. Special attention is paid to the form they assume in the presence of matter-sources and, for the two most popular conventions in the literature, to the asymptotic expressions and gravitational-wave amplitudes. Besides pointing out some inconsistencies in the literature, the expressions collected here could serve as a quick reference for the calculation of the perturbations of Schwarzschild black hole spacetimes driven by generic sources and for those approaches in which gravitational waves are extracted from numerically generated spacetimes.
The capability of the Event Horizon Telescope (EHT) to image the nearest supermassive black hole candidates at horizon-scale resolutions offers a novel means to study gravity in its strongest regimes and to test different models for these objects. Here, we study the observational appearance at 230 GHz of a surfaceless black hole mimicker, namely a non-rotating boson star, in a scenario consistent with the properties of the accretion flow onto Sgr A*. To this end, we perform general relativistic magnetohydrodynamic simulations followed by general relativistic radiative transfer calculations in the boson star space-time. Synthetic reconstructed images considering realistic astronomical observing conditions show that, despite qualitative similarities, the differences in the appearance of a black hole -- either rotating or not -- and a boson star of the type considered here are large enough to be detectable. These differences arise from dynamical effects directly related to the absence of an event horizon, in particular, the accumulation of matter in the form of a small torus or a spheroidal cloud in the interior of the boson star, and the absence of an evacuated high-magnetization funnel in the polar regions. The mechanism behind these effects is general enough to apply to other horizonless and surfaceless black hole mimickers, strengthening confidence in the ability of the EHT to identify such objects via radio observations.
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 $θ_{\rm core}\approx4^{\circ}-5^{\circ}$ and an opening angle of $θ_{\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(>Γβ)\propto(Γβ)^{-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 $θ_{\rm obs} \simeq 21^{\circ}$ for GRB170817A.
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 $μ$as could be reconstructed, allowing for stronger tests of general relativity and accretion models than with ground-based VLBI.
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.
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.
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χ$ 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χ,q)$ plane for GW190425, namely: $q_{\rm min} \geq 0.38$ and $ \tildeχ_{\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}}$.
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.
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.
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.
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 $χ_{m}$ and consider an ultrarelativistic ideal gas equation of state. For a paramagnetic fluid (i.e., with $χ_{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 $χ_{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 $τ$ with a power law $\simτ^{-a}$, two distinct solutions can be found depending on the values of $a$ and $χ_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.
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.
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.
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.
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.
As numerical calculations of inspiralling neutron-star binaries reach values of accuracy that are comparable with those of binary black holes, a fine budgeting of the various sources of error becomes increasingly important. Among such sources, the initial data is normally not accounted for, the rationale being that the error on the initial spacelike hypersurface is always far smaller than the one gained during the evolution. We here consider critically this assumption and perform a comparative analysis of the gravitational waveforms relative to essentially the same physical binary configuration when computed with two different initial-data codes, and then evolved with the same evolution code. More specifically, we consider the evolution of irrotational neutron-star binaries computed either with the pseudo-spectral code \lorene{}, or with the newly developed finite-difference code \cocal{}; both sets of initial data are subsequently evolved with the high-order evolution code \whiskythc{}. In this way we find that despite the initial data shows global (local) differences that are $\lesssim 0.02\%\ (1\%)$, the gravitational-wave phase at the merger time differs by $\sim 0.5$ radians after $\sim 3$ orbits, a surprisingly large value. Our results highlight the highly nonlinear impact that errors in the initial data can have on the subsequent evolution and the importance of using exactly the same initial data when comparative studies are done.
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.
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 $σ\approx0.002$ for aligned/anti-aligned binaries and $σ\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.
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.
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 $τ$ is the same as for the time-honored "Bjorken flow" without magnetic field. Furthermore, when the magnetic field is assumed to decay $\sim τ^{-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 present a new code for performing general-relativistic radiation-hydrodynamics simulations of accretion flows onto black holes. The radiation field is treated in the optically-thick approximation, with the opacity contributed by Thomson scattering and thermal bremsstrahlung. Our analysis is concentrated on a detailed numerical investigation of hot two-dimensional, Bondi-Hoyle accretion flows with various Mach numbers. We find significant differences with respect to purely hydrodynamical evolutions. In particular, once the system relaxes to a radiation-pressure dominated regime, the accretion rates become about two orders of magnitude smaller than in the purely hydrodynamical case, remaining however super-Eddington as are the luminosities. Furthermore, when increasing the Mach number of the inflowing gas, the accretion rates become smaller because of the smaller cross section of the black hole, but the luminosities increase as a result a stronger emission in the shocked regions. Overall, our approach provides the first self-consistent calculation of the Bondi-Hoyle luminosity, most of which is emitted within r~100 M from the black hole, with typical values L/L_Edd ~ 1-7, and corresponding energy efficiencies eta_BH ~ 0.09-0.5. The possibility of computing luminosities self-consistently has also allowed us to compare with the bremsstrahlung luminosity often used in modelling the electromagnetic counterparts to supermassive black-hole binaries, to find that in the optically-thick regime these more crude estimates are about 20 times larger than our radiation-hydrodynamics results.
We 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$.
We present the extension of our \cocal~- Compact Object CALculator - code to compute general-relativistic initial data for binary compact-star systems. In particular, we construct quasiequilibrium initial data for equal-mass binaries with spins that are either aligned or antialigned with the orbital angular momentum. The Isenberg-Wilson-Mathews formalism is adopted and the constraint equations are solved using the representation formula with a suitable choice of a Green's function. We validate the new code with solutions for equal-mass binaries and explore its capabilities for a wide range of compactnesses, from a white dwarf binary with compactness $\sim 10^{-4}$, up to a highly relativistic neutron-star binary with compactness $\sim 0.22$. We also present a comparison with corotating and irrotational quasiequilibrium sequences from the spectral code \lorene [Taniguchi and Gourgoulhon, Phys. Rev. D {\bf 66}, 104019 (2002)] and with different compactness, showing that the results from the two codes agree to a precision of the order of $0.05\%$. Finally, we present equilibria for spinning configurations with a nuclear-physics equation of state in a piecewise polytropic representation.
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.
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.
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.
Models of rotating relativistic stars with a toroidal magnetic field have been computed for a sample of eight equations of state of cold dense matter. Non-rotating models admit important levels of magnetization and quadrupole distortion accompanied by a seemingly unlimited growth in size. Rotating models reach the mass-shedding limit at smaller angular velocities than in the non-magnetized case according to the larger circumferential equatorial radius induced by the magnetic field. Moreover, they can be classified as prolate-prolate, oblate-prolate, or oblate-oblate with respect to surface deformation and quadrupole distortion. Simple expressions for surface and quadrupole deformation are provided that are valid up to magnetar field strengths and rapid rotation.
We have computed models of rotating relativistic stars with a toroidal magnetic field and investigated the combined effects of magnetic field and rotation on the apparent shape (i.e. the surface deformation), which could be relevant for the electromagnetic emission, and on the internal matter distribution (i.e. the quadrupole distortion), which could be relevant for the emission of gravitational waves. Using a sample of eight different cold nuclear physics equations of state, we have computed models of maximum field strength, as well as the distortion coefficients for the surface and the quadrupolar deformations. Surprisingly, we find that non-rotating models admit arbitrary levels of magnetization, accompanied by a growth of size and quadrupole distortion to which we could not find a limit. Rotating models, on the other hand, are subject to a mass-shedding limit at frequencies well below the corresponding ones for unmagnetized stars. Overall, the space of solutions can be split into three distinct classes for which the surface deformation and the quadrupole distortion are either prolate and prolate, oblate and prolate, or oblate and oblate, respectively. We also derive a simple formula expressing the relativistic distortion coefficients, which allows one to compute the surface deformation and the quadrupole distortion up to significant levels of rotation and magnetization, essentially covering all known magnetars. Such a formula replaces Newtonian equivalent expressions that overestimate the magnetic quadrupole distortion by about a factor of 6 and are inadequate for strongly relativistic objects like neutron stars.
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.
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.
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 propose a new parametric framework to describe in generic metric theories of gravity the spacetime of spherically symmetric and slowly rotating black holes. In contrast to similar approaches proposed so far, we do not use a Taylor expansion in powers of M/r, where M and r are the mass of the black hole and a generic radial coordinate, respectively. Rather, we use a continued-fraction expansion in terms of a compactified radial coordinate. This choice leads to superior convergence properties and allows us to approximate a number of known metric theories with a much smaller set of coefficients. The measure of these coefficients via observations of near-horizon processes can be used to effectively constrain and compare arbitrary metric theories of gravity. Although our attention is here focussed on spherically symmetric black holes, we also discuss how our approach could be extended to rotating black holes.
We extend the Wald solution to a black hole that is also boosted. More specifically, we derive analytic solutions for the Maxwell equations for a rotating black hole moving at constant speed in an asymptotically uniform magnetic test field. By adopting Kerr-Schild coordinates we avoid singular behaviours at the horizon and obtain a complete description of the charge and current distributions in terms of the black-hole spin and boost. Using this solution, we compute the energy losses expected when charged particles are accelerated along the magnetic field lines, improving previous estimates that had to cope with singular electromagnetic fields on the horizon. When used to approximate the emission from binary black holes in a uniform magnetic field, our estimates match reasonably well those from numerical-relativity calculations in the force-free approximation.
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.
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.
We probe the gravitational interaction of two black holes in the strong-field regime by computing the scattering angle $χ$ of hyperbolic-like, close binary-black-hole encounters as a function of the impact parameter. The fully general-relativistic result from numerical relativity is compared to two analytic approximations: post-Newtonian theory and the effective-one-body formalism. As the impact parameter decreases, so that black holes pass within a few times their Schwarzschild radii, we find that the post-Newtonian prediction becomes quite inaccurate, while the effective-one-body one keeps showing a good agreement with numerical results. Because we have explored a regime which is very different from the one considered so far with binaries in quasi-circular orbits, our results open a new avenue to improve analytic representations of the general-relativistic two-body Hamiltonian.
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.
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.
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.
Following previous work in vacuum spacetimes, we investigate the constraint-damping properties in the presence of matter of the recently developed traceless, conformal and covariant Z4 (CCZ4) formulation of the Einstein equations. First, we evolve an isolated neutron star with an ideal gas equation of state and subject to a constraint-violating perturbation. We compare the evolution of the constraints using the CCZ4 and Baumgarte-Shibata-Shapiro-Nakamura-Oohara-Kojima (BSSNOK) systems. Second, we study the collapse of an unstable spherical star to a black hole. Finally, we evolve binary neutron star systems over several orbits until the merger, the formation of a black hole, and up to the ringdown. We show that the CCZ4 formulation is stable in the presence of matter and that the constraint violations are one or more orders of magnitude smaller than for the BSSNOK formulation. Furthermore, by comparing the CCZ4 and the BSSNOK formulations also for neutron star binaries with large initial constraint violations, we investigate their influence on the errors on physical quantities. We also give a new, simple and robust prescription for the damping parameter that removes the instabilities found when using the fully covariant version of CCZ4 in the evolution of black holes. Overall, we find that at essentially the same computational costs the CCZ4 formulation provides solutions that are stable and with a considerably smaller violation of the Hamiltonian constraint than the BSSNOK formulation. We also find that the performance of the CCZ4 formulation is very similar to another conformal and traceless, but noncovariant formulation of the Z4 system, i.e. the Z4c formulation.
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.
The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters.
We present a new numerical implementation of the general-relativistic resistive magnetohydrodynamics (MHD) equations within the Whisky code. The numerical method adopted exploits the properties of implicit-explicit Runge-Kutta numerical schemes to treat the stiff terms that appear in the equations for large electrical conductivities. Using tests in one, two, and three dimensions, we show that our implementation is robust and recovers the ideal-MHD limit in regimes of very high conductivity. Moreover, the results illustrate that the code is capable of describing scenarios in a very wide range of conductivities. In addition to tests in flat spacetime, we report simulations of magnetized nonrotating relativistic stars, both in the Cowling approximation and in dynamical spacetimes. Finally, because of its astrophysical relevance and because it provides a severe testbed for general-relativistic codes with dynamical electromagnetic fields, we study the collapse of a nonrotating star to a black hole. We show that also in this case our results on the quasinormal mode frequencies of the excited electromagnetic fields in the Schwarzschild background agree with the perturbative studies within 0.7% and 5.6% for the real and the imaginary part of the l=1 mode eigenfrequency, respectively. Finally we provide an estimate of the electromagnetic efficiency of this process.
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$.
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.
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 Λ 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 Λ, and changes in the EOS that did not change Λ are not measurable. We estimate with two methods the minimal and expected measurability of Λ 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.
The merger of two neutron stars will in general lead to the formation of a torus surrounding a black hole whose rotational energy can be tapped to potentially power a short gamma-ray burst. We have studied the merger of equal-mass binaries with spins aligned with the orbital angular momentum to determine the maximum spin the black hole can reach. Our initial data consists of irrotational binaries to which we add various amounts of rotation to increase the total angular momentum. Although the initial data violates the constraint equations, the use of the constraint-damping CCZ4 formulation yields evolutions with violations smaller than those with irrotational initial data and standard formulations. Interestingly, we find that a limit of $J/M^2 \simeq 0.89$ exists for the dimensionless spin and that any additional angular momentum given to the binary ends up in the torus rather than in the black hole, thus providing another nontrivial example supporting the cosmic censorship hypothesis.
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.
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.
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.
We provide additional information on our recent study of the electromagnetic emission 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. As anticipated in a recent letter, our results show that although a dual-jet structure is present, the associated luminosity is ~ 100 times smaller than the total one, which is predominantly quadrupolar. We here discuss the details of our implementation of the equations in which the force-free condition is not implemented at a discrete level, but rather obtained via a damping scheme which drives the solution to satisfy the correct condition. We show that this is important for a correct and accurate description of the current sheets that can develop in the course of the simulation. We also study in greater detail the three-dimensional charge distribution produced as a consequence of the inspiral and show that during the inspiral it possesses a complex but ordered structure which traces the motion of the two black holes. Finally, we provide quantitative estimates of the scaling of the electromagnetic emission with frequency, with the diffused part having a dependence that is the same as the gravitational-wave one and that scales as L ~ Omega^{10/3-8/3}, while the collimated one scales as L ~ Omega^{5/3-6/3}, thus with a steeper dependence than previously estimated. We discuss the impact of these results on the potential detectability of dual jets from supermassive black holes and the steps necessary for more accurate estimates.
The understanding of strong-field dynamics near black-hole horizons is a long-standing and challenging prob- lem in general relativity. Recent advances in numerical relativity and in the geometric characterization of black- hole horizons open new avenues into the problem. In this first paper in a series of two, we focus on the analysis of the recoil occurring in the merger of binary black holes, extending the analysis initiated in [1] with Robinson- Trautman spacetimes. More specifically, we probe spacetime dynamics through the correlation of quantities defined at the black-hole horizon and at null infinity. The geometry of these hypersurfaces responds to bulk gravitational fields acting as test screens in a scattering perspective of spacetime dynamics. Within a 3 + 1 approach we build an effective-curvature vector from the intrinsic geometry of dynamical-horizon sections and correlate its evolution with the flux of Bondi linear momentum at large distances. We employ this setup to study numerically the head-on collision of nonspinning black holes and demonstrate its validity to track the qualita- tive aspects of recoil dynamics at infinity. We also make contact with the suggestion that the antikick can be described in terms of a "slowness parameter" and how this can be computed from the local properties of the horizon. In a companion paper [2] we will further elaborate on the geometric aspects of this approach and on its relation with other approaches to characterize dynamical properties of black-hole horizons.
In a companion paper [1], we have presented a cross-correlation approach to near-horizon physics in which bulk dynamics is probed through the correlation of quantities defined at inner and outer spacetime hypersurfaces acting as test screens. More specifically, dynamical horizons provide appropriate inner screens in a 3+1 setting and, in this context, we have shown that an effective-curvature vector measured at the common horizon produced in a head-on collision merger can be correlated with the flux of linear Bondi-momentum at null infinity. In this paper we provide a more sound geometric basis to this picture. First, we show that a rigidity property of dynamical horizons, namely foliation uniqueness, leads to a preferred class of null tetrads and Weyl scalars on these hypersurfaces. Second, we identify a heuristic horizon news-like function, depending only on the geometry of spatial sections of the horizon. Fluxes constructed from this function offer refined geometric quantities to be correlated with Bondi fluxes at infinity, as well as a contact with the discussion of quasi-local 4-momentum on dynamical horizons. Third, we highlight the importance of tracking the internal horizon dual to the apparent horizon in spatial 3-slices when integrating fluxes along the horizon. Finally, we discuss the link between the dissipation of the non-stationary part of the horizon's geometry with the viscous-fluid analogy for black holes, introducing a geometric prescription for a "slowness parameter" in black-hole recoil dynamics.
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 <γ> following a simple scaling of the type M_c ~ K <γ>^{-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.
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.
We present a new formulation of the Einstein equations based on a conformal and traceless decomposition of the covariant form of the Z4 system. This formulation combines the advantages of a conformal decomposition, such as the one used in the BSSNOK formulation (i.e. well-tested hyperbolic gauges, no need for excision, robustness to imperfect boundary conditions) with the advantages of a constraint-damped formulation, such as the generalized harmonic one (i.e. exponential decay of constraint violations when these are produced). We validate the new set of equations through standard tests and by evolving binary black hole systems. Overall, the new conformal formulation leads to a better behavior of the constraint equations and a rapid suppression of the violations when they occur. The changes necessary to implement the new conformal formulation in standard BSSNOK codes are very small as are the additional computational costs.
We revisit the suggestion that dual jets can be produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. By performing independent calculations and by computing the electromagnetic emission in a way which is consistent with estimates using the Poynting flux, we show that a dual-jet structure is present but energetically subdominant with respect to a non-collimated and predominantly quadrupolar emission, which is similar to the one computed when the binary is in electrovacuum. While our findings set serious restrictions on the detectability of dual jets from coalescing binaries, they also increase the chances of detecting an EM counterpart from these systems.
We investigate the 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..
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.
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.
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.
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.
Binary neutron-star systems represent one of the most promising sources of gravitational waves. In order to be able to extract important information, notably about the equation of state of matter at nuclear density, it is necessary to have in hands an accurate analytical model of the expected waveforms. Following our recent work, we here analyze more in detail two general-relativistic simulations spanning about 20 gravitational-wave cycles of the inspiral of equal-mass binary neutron stars with different compactnesses, and compare them with a tidal extension of the effective-one-body (EOB) analytical model. The latter tidally extended EOB model is analytically complete up to the 1.5 post-Newtonian level, and contains an analytically undetermined parameter representing a higher-order amplification of tidal effects. We find that, by calibrating this single parameter, the EOB model can reproduce, within the numerical error, the two numerical waveforms essentially up to the merger. By contrast, analytical models (either EOB, or Taylor-T4) that do not incorporate such a higher-order amplification of tidal effects, build a dephasing with respect to the numerical waveforms of several radians.
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.
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.
We investigate the stationary end-state obtained by evolving a collapsing spherical star with the gauges routinely adopted to study puncture black holes. We compare the end-state of the collapse with the trumpet solution found in the evolution of a single wormhole slice and show that the two solutions closely agree. We demonstrate that the agreement is caused by the use of the Gamma-driver shift condition, which allows the matter to fall inwards into a region of spacetime that is not resolved by the numerical grid, and which simultaneously finds the stationary coordinates of the trumpet outside the matter.
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.
We present a new three-dimensional general relativistic hydrodynamics code, the Whisky code. This code incorporates the expertise developed over the past years in the numerical solution of Einstein equations and of the hydrodynamics equations in a curved spacetime, and is the result of a collaboration of several European Institutes. We here discuss the ability of the code to carry out long-term accurate evolutions of the linear and nonlinear dynamics of isolated relativistic stars.
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 Γ=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.
To detect the gravitational-wave (GW) signal from binary neutron stars and extract information about the equation of state of matter at nuclear density, it is necessary to match the signal with a bank of accurate templates. We present the two longest (to date) general-relativistic simulations of equal-mass binary neutron stars with different compactnesses, C=0.12 and C=0.14, and compare them with a tidal extension of the effective-one-body (EOB)model. The typical numerical phasing errors over the $\simeq 22$ GW cycles are $Δϕ\simeq \pm 0.24$ rad. By calibrating only one parameter (representing a higher-order amplification of tidal effects), the EOB model can reproduce, within the numerical error, the two numerical waveforms essentially up to the merger. By contrast, the third post-Newtonian Taylor-T4 approximant with leading-order tidal corrections dephases with respect to the numerical waveforms by several radians.
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_ϕ$ and the self-interacting coupling constant $λ$ 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 consider the evolution in full general relativity of a family of linearly unstable isolated spherical neutron stars under the effects of very small, perturbations as induced by the truncation error. Using a simple ideal-fluid equation of state we find that this system exhibits a type-I critical behaviour, thus confirming the conclusions reached by Liebling et al. [1] for rotating magnetized stars. Exploiting the relative simplicity of our system, we are able carry out a more in-depth study providing solid evidences of the criticality of this phenomenon and also to give a simple interpretation of the putative critical solution as a spherical solution with the unstable mode being the fundamental F-mode. Hence for any choice of the polytropic constant, the critical solution will distinguish the set of subcritical models migrating to the stable branch of the models of equilibrium from the set of subcritical models collapsing to a black hole. Finally, we study how the dynamics changes when the numerically perturbation is replaced by a finite-size, resolution independent velocity perturbation and show that in such cases a nearly-critical solution can be changed into either a sub or supercritical. The work reported here also lays the basis for the analysis carried in a companion paper, where the critical behaviour in the the head-on collision of two neutron stars is instead considered [2].
We consider the head-on collision of equal-mass neutron stars boosted towards each other and we study the behavior of such systems near the threshold of black-hole formation. In particular, we confirm the previous findings by [1] that a type-I critical phenomenon can be observed by fine-tuning the initial mass of the two neutron stars. At the same time, we argue against the interpretation that the critical solution is not a perturbed spherical star and show instead that the metastable star corresponds to a (perturbed) equilibrium solution on the unstable branch of the equilibrium configurations. As a result, the head-on collision of two neutron stars near the critical threshold can be seen as a transition in the space of configurations from an initial stable solution over to a critical metastable one which can either migrate to a stable solution or collapse to a black hole. The critical exponent for this process shows a fine structure which was already observed in the case of the critical collapse of scalar fields but never before for perfect fluids.
Suitable gauge conditions are fundamental for stable and accurate numerical-relativity simulations of inspiralling compact binaries. A number of well-studied conditions have been developed over the last decade for both the lapse and the shift and these have been successfully used both in vacuum and non-vacuum spacetimes when simulating binaries with comparable masses. At the same time, recent evidence has emerged that the standard "Gamma-driver" shift condition requires a careful and non-trivial tuning of its parameters to ensure long-term stable evolutions of unequal-mass binaries. We present a novel gauge condition in which the damping constant is promoted to be a dynamical variable and the solution of an evolution equation. We show that this choice removes the need for special tuning and provides a shift damping term which is free of instabilities in our simulations and dynamically adapts to the individual positions and masses of the binary black-hole system. Our gauge condition also reduces the variations in the coordinate size of the apparent horizon of the larger black hole and could therefore be useful when simulating binaries with very small mass ratios.
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.
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.
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.
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.
The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.
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.
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)
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.
The final evolution of a binary black-hole system gives rise to a recoil velocity if an asymmetry is present in the emitted gravitational radiation. Measurements of this effect for non-spinning binaries with unequal masses have pointed out that kick velocities $\sim~175$ km/s can be reached for a mass ratio $\simeq 0.36$. However, a larger recoil can be obtained for equal-mass binaries if the asymmetry is provided by the spins. Using two independent methods we show that the merger of such binaries yields velocities as large as $\sim 440$ km/s for black holes having unequal spins that are antialigned and parallel to the orbital angular momentum.
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.
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.