Source author record

Bruno Giacomazzo

Bruno Giacomazzo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2026arXiv

ESO Expanding Horizon White Paper: Revealing the properties of matter at supranuclear densities with gravitational waves

Understanding dense matter under extreme conditions is one of the most fundamental puzzles in modern physics. Complex interactions give rise to emergent, collective phenomena. While nuclear experiments and Earth - based colliders provide valuable insights, much of the quantum chromodynamics phase diagram at high density and low temperature remains accessible only through astrophysical observations of neutron stars, neutron star mergers, and stellar collapse. Astronomical observations thus offer a direct window to the physics on subatomic scales with gravitational waves presenting an especially clean channel. Next-generation gravitational - wave observatories, such as the Einstein Telescope, would serve as unparalleled instruments to transform our understanding of neutron star matter. They will enable the detection of up to tens of thousands of binary neutron star and neutron star - black hole mergers per year, a dramatic increase over the few events accessible with current detectors. They will provide an unprecedented precision in probing cold, dense matter during the binary inspiral, exceeding by at least an order of magnitude what current facilities can achieve. Moreover, these observatories will allow us to explore uncharted regimes of dense matter at finite temperatures produced in a subset of neutron star mergers, areas that remain entirely inaccessible to current instruments. Together with multimessenger observations, these measurements will significantly deepen our knowledge of dense nuclear matter.

preprint2025arXiv

Multi-messenger observations in the Einstein Telescope era: binary neutron star and black hole - neutron star mergers

The Einstein Telescope (ET), a proposed next-generation gravitational wave (GW) observatory, will expand the reach of GW astronomy of stellar-mass compact object binaries to unprecedented distances, enhancing opportunities for multi-messenger observations. Here we investigate multi-messenger emission properties of binary neutron star (NSNS) and black hole-neutron star (BHNS) mergers detectable by ET, providing projections to optimize observational strategies and maximize scientific insights from these sources. Using a synthetic population of compact binary mergers, we model each source's GW signal-to-noise ratio, sky localization uncertainty, kilonova (KN) light curves in optical and near-infrared bands, fluence of the relativistic jet gamma-ray burst (GRB) prompt emission and afterglow light curves across radio, optical, X-ray and very high energy wavelengths. We analyze multi-messenger detectability prospects for ET as a standalone observatory with two different configurations and within a network of next-generation GW detectors. ET will detect over $10^4$ NSNS mergers annually, enabling potential observation of tens to hundreds of electromagnetic (EM) counterparts. BHNS mergers have more limited multi-messenger prospects, but joint GW-EM rates will increase by an order of magnitude compared to current-generation instruments. We quantify uncertainties due to the NS equation of state (EoS) and mass distribution of NSNSs, as well as the NS EoS and BH spin for BHNSs. While a single ET will achieve an impressive GW detection rate, the fraction of well-localized events is orders of magnitude lower than in a network with additional detectors. This significantly limits efficient EM follow-up and science cases requiring well-characterized counterparts or early observations. The challenge is even greater for BHNS mergers due to their low EM rate.

preprint2024arXiv

GRMHD study of accreting massive black hole binaries in astrophysical environment: a review

We present an overview of recent numerical advances in the theoretical characterization of massive binary black hole (MBBH) mergers in astrophysical environments. These systems are among the loudest sources of gravitational waves (GWs) in the universe and particularly promising candidates for multimessenger astronomy. Coincident detection of GWs and electromagnetic (EM) signals from merging MBBHs is at the frontier of contemporary astrophysics. One major challenge in observational efforts searching for these systems is the scarcity of strong predictions for EM signals arising before, during, and after merger. Therefore, a great effort in theoretical work to-date has been to characterize EM counterparts emerging from MBBHs concurrently to the GW signal, aiming to determine distinctive observational features that will guide and assist EM observations. To produce sharp EM predictions of MBBH mergers it is key to model the binary inspiral down to coalescence in a full general relativistic fashion by solving Einstein's field equations coupled with the magnetohydrodynamics equations that govern the evolution of the accreting plasma in strong-gravity. We review the general relativistic numerical investigations that have explored the astrophysical manifestations of MBBH mergers in different environments and focused on predicting potentially observable smoking-gun EM signatures that accompany the gravitational signal.

preprint2022arXiv

Accretion-to-jet energy conversion efficiency in GW170817

Based on previously published multi-wavelength modelling of the GRB 170817A jet afterglow, that includes information from the VLBI centroid motion, we construct the posterior probability density distribution on the total energy in the bipolar jets launched by the GW170817 merger remnant. By applying a new numerical-relativity-informed fitting formula for the accretion disk mass, we construct the posterior probability density distribution of the GW170817 remnant disk mass. By combining the two, we estimate the accretion-to-jet energy conversion efficiency in this system, carefully accounting for uncertainties. The accretion-to-jet energy conversion efficiency in GW170817 is $η\sim 10^{-3}$ with an uncertainty of slightly less than two orders of magnitude. This low efficiency is in good agreement with expectations from the $ν\barν$ mechanism, which therefore cannot be excluded by this measurement alone. Such an efficiency also agrees with that anticipated for the Blandford-Znajek mechanism, provided that the magnetic field in the disk right after the merger is predominantly toroidal (which is expected as a result of the merger dynamics).

preprint2022arXiv

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

preprint2022arXiv

Implementing a new recovery scheme for primitive variables in the general relativistic magnetohydrodynamic code Spritz

General relativistic magnetohydrodynamic (GRMHD) simulations represent a fundamental tool to probe various underlying mechanisms at play during binary neutron star (BNS) and neutron star (NS) - black hole (BH) mergers. Contemporary flux-conservative GRMHD codes numerically evolve a set of conservative equations based on `conserved' variables which then need to be converted back into the fundamental (`primitive') variables. The corresponding conservative-to-primitive variable recovery procedure, based on root-finding algorithms, constitutes one of the core elements of such GRMHD codes. Recently, a new robust, accurate and efficient recovery scheme called RePrimAnd was introduced, which has demonstrated the ability to always converge to a unique solution. The scheme provides fine-grained error policies to handle invalid states caused by evolution errors, and also provides analytical bounds for the error of all primitive variables. In this work, we describe the technical aspects of implementing the RePrimAnd scheme into the GRMHD code Spritz. To check our implementation as well as to assess the various features of the scheme, we perform a number of GRMHD tests in three dimensions. Our tests, which include critical cases such as a NS collapse to a BH as well as the early evolution (~50 ms) of a Fishbone-Moncrief BH-accrection disk system, show that RePrimAnd is able to support magnetized, low density environments with magnetic-to-fluid pressure ratios as high as 10^4, in situations where the previously used recovery scheme fails.

preprint2022arXiv

Misaligned Spinning Binary Black Hole Mergers in Hot Magnetized Plasma

We present general relativistic magneto-hydrodynamical simulations of equal-mass spinning black hole binary mergers embedded in a magnetized gas cloud. We focus on the effect of the spin orientation relative to the orbital angular momentum on the flow dynamics, mass accretion rate and Poynting luminosity. We find that, across the inspiral, the gas accreting onto the individual black holes concentrates into disk-like overdensities, whose angular momenta are oriented towards the spin axes and which persist until merger. We identify quasi-periodic modulations occurring in the mass accretion rate at the level of 1-20%, evolving in parallel with the gravitational wave chirp. The similarity between the accretion rate time-series and the gravitational strain is a consequence of the interplay between strong, dynamical gravitational fields and magnetic fields in the vicinity of the inspiralling black holes. This result suggests that quasi-periodicity in the pre-merger accretion rate of massive binaries is not exclusive of environments in which the black holes are embedded in a circumbinary accretion disk, and could provide an additional useful signature of electromagnetic emission concurrent to low-frequency gravitational wave detection.

preprint2020arXiv

Spritz: a new fully general-relativistic magnetohydrodynamic code

The new era of multimessenger astrophysics requires the capability of studying different aspects of the evolution of compact objects. In particular, the merger of neutron star binaries is a strong source of gravitational waves and electromagnetic radiation, from radio to $γ$-rays, as demonstrated by the detection of GW170817 and its electromagnetic counterparts. In order to understand the physical mechanisms involved in such systems, it is necessary to employ fully general relativistic magnetohydrodynamic (GRMHD) simulations able to include the effects of a composition and temperature dependent equation of state describing neutron star matter as well as neutrino emission and reabsorption. Here, we present our new code named $\texttt{Spritz}$ that solves the GRMHD equations in 3D Cartesian coordinates and on a dynamical spacetime. The code can support tabulated equations of state, taking into account finite temperature effects and allowing for the inclusion of neutrino radiation. In this first paper, we present the general features of the code and a series of tests performed in special and general relativity to assess the robustness of the basic GRMHD algorithms implemented. Among these tests, we also present the first comparison between a non-staggered and a staggered formulation of the vector potential evolution, which is used to guarantee the divergence-less character of the magnetic field. With respect to other publicly available GRMHD codes, $\texttt{Spritz}$ combines the robust approach of a staggered formulation of the vector potential together with the use of an equation of state driver ($\texttt{EOS_Omni}$) that can allow the code to use finite temperature equations of state. A next version of the code will fully test the $\texttt{EOS_Omni}$ driver by coupling it with a neutrino leakage scheme.

preprint2016arXiv

Binary Neutron Star Mergers and Short Gamma-Ray Bursts: Effects of Magnetic Field Orientation, Equation of State, and Mass Ratio

We present fully GRMHD simulations of the merger of binary neutron star (BNS) systems. We consider BNSs producing a hypermassive neutron star (HMNS) that collapses to a spinning black hole (BH) surrounded by a magnetized accretion disk in a few tens of ms. We investigate whether such systems may launch relativistic jets and power short gamma-ray bursts. We study the effects of different equations of state (EOSs), different mass ratios, and different magnetic field orientations. For all cases, we present a detailed investigation of the matter dynamics and of the magnetic field evolution, with particular attention to its global structure and possible emission of relativistic jets. The main result of this work is that we found the formation of an organized magnetic field structure. This happens independently of EOS, mass ratio, and initial magnetic field orientation. We also show that those models that produce a longer-lived HMNS lead to a stronger magnetic field before collapse to BH. Such larger fields make it possible, for at least one of our models, to resolve the MRI and hence further amplify the magnetic field. However, by the end of our simulations, we do not observe a magnetically dominated funnel and hence neither a relativistic outflow. With respect to the recent simulations of Ruiz et al 2016, we evolve models with lower and more realistic initial magnetic field strengths and, because of computational reasons, we do not evolve the accretion disk for the long timescales that seem to be required in order to see a relativistic outflow. Since all our models produce a similar ordered magnetic field structure, we expect that the results found in Ruiz et al 2016, where they only considered an equal-mass system with an ideal fluid EOS, should be general and, at least from a qualitative point of view, independent from mass-ratio, magnetic field orientation, and EOS.

preprint2016arXiv

General Relativistic Magnetohydrodynamic Simulations of Binary Neutron Star Mergers with the APR4 Equation of State

We present new results of fully general relativistic magnetohydrodynamic simulations of binary neutron star (BNS) mergers performed with the Whisky code. All the models use a piecewise polytropic approximation of the APR4 equation of state for cold matter, together with a 'hybrid' part to incorporate thermal effects during the evolution. We consider both equal and unequal-mass models, with total masses such that either a supramassive NS or a black hole is formed after merger. Each model is evolved with and without a magnetic field initially confined to the stellar interior. We present the different gravitational wave (GW) signals as well as a detailed description of the matter dynamics (magnetic field evolution, ejected mass, post-merger remnant/disk properties). Our simulations provide new insights into BNS mergers, the associated GW emission and the possible connection with the engine of short gamma-ray bursts (both in the 'standard' and in the 'time-reversal' scenarios) and other electromagnetic counterparts.

preprint2016arXiv

Short Gamma-Ray Bursts from the Merger of Two Black Holes

Short Gamma-Ray Bursts (GRBs) are explosions of cosmic origin believed to be associated with the merger of two compact objects, either two neutron stars, or a neutron star and a black hole. The presence of at least one neutron star has long been thought to be an essential element of the model: its tidal disruption provides the needed baryonic material whose rapid accretion onto the post-merger black hole powers the burst. The recent tentative detection by the Fermi satellite of a short GRB in association with the gravitational wave signal GW150914 produced by the merger of two black holes has challenged this standard paradigm. Here we show that the evolution of two high-mass, low-metallicity stars with main sequence rotational speeds a few tens of percent of the critical speed eventually undergoing a weak supernova explosion {\em can} produce a short gamma-ray burst. The outer layers of the envelope of the last exploding star remain bound and circularize at large radii. With time, the disk cools and becomes neutral, suppressing the magneto-rotational instability, and hence the viscosity. The disk remains 'long-lived dead' until tidal torques and shocks during the pre-merger phase heat it up and re-ignite accretion, rapidly consuming the disk and powering the short gamma-ray burst.

preprint2015arXiv

Producing Magnetar Magnetic Fields in the Merger of Binary Neutron Stars

The merger of binary neutron stars (BNSs) can lead to large amplifications of the magnetic field due to the development of turbulence and instabilities in the fluid, such as the Kelvin-Helmholtz shear instability, which drive small-scale dynamo activity. In order to properly resolve such instabilities and obtain the correct magnetic field amplification, one would need to employ resolutions that are currently unfeasible in global general relativistic magnetohydrodynamic (GRMHD) simulations of BNS mergers. Here, we present a subgrid model that allows global simulations to take into account the small-scale amplification of the magnetic field which is caused by the development of turbulence during BNS mergers. Assuming dynamo saturation, we show that magnetar-level fields ($\sim 10^{16}\,{\rm G}$) can be easily reached, and should therefore be expected from the merger of magnetized BNSs. The total magnetic energy can reach values up to $\sim 10^{51}\,{\rm erg}$ and the post-merger remnant can therefore emit strong electromagnetic signals and possibly produce short gamma-ray bursts.

preprint2014arXiv

Gravitational waves from massive magnetars formed in binary neutron star mergers

Binary neutron star (NS) mergers are among the most promising sources of gravitational waves (GWs), as well as candidate progenitors for short Gamma-Ray Bursts (SGRBs). Depending on the total initial mass of the system, and the NS equation of state (EOS), the post-merger phase can see a prompt collapse to a black hole, or the formation of a supramassive NS, or even a stable NS. In the case of post-merger NS (PMNS) formation, magnetic field amplification during the merger will produce a magnetar with a large induced mass quadrupole moment, and millisecond spin. If the timescale for orthogonalization of the magnetic symmetry axis with the spin axis is sufficiently short the NS will radiate its spin down energy primarily via GWs. Here we study this scenario for various outcomes of NS formation: we generalise the set of equilibrium states for a twisted torus magnetic configuration to include solutions that, at a fixed exterior dipole field, carry a larger magnetic energy reservoir; we hence compute their magnetic ellipticity and the strength of the expected GW signal as a function of the magnitude of the dipole and toroidal field. The relative number of GW detections from PMNSs and from binary NSs is a strong function of the NS equation of state (EOS), being higher (~ 1%) for the stiffest EOSs and negligibly small for the softest ones. For intermediate-stiffness EOSs, such as the n=4/7 polytrope recently used by Giacomazzo \& Perna or the GM1 used by Lasky et al., the relative fraction is ~0.3%; correspondingly we estimate a GW detection rate from stable PMNSs of ~ (0.1-1) yr$^{-1}$ with Advanced detectors, and of ~ (100-1000) yr$^{-1}$ with third generation detectors such as the Einstein Telescope. Measurement of such GW signal would provide strong constraints on the NS EOS and on the nature of the binary progenitors giving rise to SGRBs.

preprint2013arXiv

Formation of Stable Magnetars from Binary Neutron Star Mergers

By performing fully general relativistic magnetohydrodynamic simulations of binary neutron star mergers, we investigate the possibility that the end result of the merger is a stable magnetar. In particular, we show that, for a binary composed of two equal-mass neutron stars (NSs) of gravitational mass M~1.2 Msun and equation of state similar to Shen et al. at high densities, the merger product is a stable NS. Such NS is found to be differentially rotating and ultraspinning with spin parameter J/M^2~0.86, where J is its total angular momentum, and it is surrounded by a disk of ~0.1 Msun. While in our global simulations the magnetic field is amplified by about two orders of magnitude, local simulations have shown that hydrodynamic instabilities and the onset of the magnetorotational instability could further increase the magnetic field strength up to magnetar levels. This leads to the interesting possibility that, for some NS mergers, a stable and magnetized NS surrounded by an accretion disk could be formed. We discuss the impact of these new results for the emission of electromagnetic counterparts of gravitational wave signals and for the central engine of short gamma-ray bursts.

preprint2013arXiv

General-relativistic resistive magnetohydrodynamics in three dimensions: Formulation and tests

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.

preprint2013arXiv

Matter effects on binary neutron star waveforms

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.

preprint2013arXiv

The Transient Gravitational-Wave Sky

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

preprint2012arXiv

Compact Binary Progenitors of Short Gamma-Ray Bursts

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.

preprint2012arXiv

General Relativistic Simulations of Accretion Induced Collapse of Neutron Stars to Black Holes

Neutron stars (NSs) in the astrophysical universe are often surrounded by accretion disks. Accretion of matter onto an NS may increase its mass above the maximum value allowed by its equation of state, inducing its collapse to a black hole (BH). Here we study this process for the first time, in three-dimensions, and in full general relativity. By considering three initial NS configurations, each with and without a surrounding disk (of mass ~7% M_{NS}), we investigate the effect of the accretion disk on the dynamics of the collapse and its imprint on both the gravitational wave (GW) and electromagnetic (EM) signals that can be emitted by these sources. We show in particular that, even if the GW signal is similar for the accretion induced collapse (AIC) and the collapse of an NS in vacuum (and detectable only for Galactic sources), the EM counterpart could allow us to discriminate between these two types of events. In fact, our simulations show that, while the collapse of an NS in vacuum leaves no appreciable baryonic matter outside the event horizon, an AIC is followed by a phase of rapid accretion of the surviving disk onto the newly formed BH. The post-collapse accretion rates, on the order of ~10^{-2} M_{sun} s^{-1}, make these events tantalizing candidates as engines of short gamma-ray bursts.

preprint2012arXiv

General Relativistic Simulations of Magnetized Plasmas around Merging Supermassive Black Holes

Coalescing supermassive black hole binaries are produced by the mergers of galaxies and are the most powerful sources of gravitational waves accessible to space-based gravitational observatories. Some such mergers may occur in the presence of matter and magnetic fields and hence generate an electromagnetic counterpart. In this Letter, we present the first general relativistic simulations of magnetized plasma around merging supermassive black holes using the general relativistic magnetohydrodynamic code Whisky. By considering different magnetic field strengths, going from non-magnetically dominated to magnetically dominated regimes, we explore how magnetic fields affect the dynamics of the plasma and the possible emission of electromagnetic signals. In particular we observe a total amplification of the magnetic field of ~2 orders of magnitude which is driven by the accretion onto the binary and that leads to much stronger electromagnetic signals, more than a factor of 10^4 larger than comparable calculations done in the force-free regime where such amplifications are not possible.

preprint2011arXiv

Accurate evolutions of inspiralling and magnetized neutron-stars: equal-mass binaries

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.

preprint2011arXiv

Accurate numerical simulations of inspiralling binary neutron stars and their comparison with effective-one-body analytical models

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.

preprint2011arXiv

Collapse of differentially rotating neutron stars and cosmic censorship

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.

preprint2011arXiv

The missing link: Merging neutron stars naturally produce jet-like structures and can power short Gamma-Ray Bursts

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.

preprint2010arXiv

Accurate evolutions of unequal-mass neutron-star binaries: properties of the torus and short GRB engines

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.

preprint2010arXiv

Analytic modelling of tidal effects in the relativistic inspiral of binary neutron stars

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.