-
Cosmic twins are likely rare: exploring pathways to hadron-quark stars with the same mass as neutron stars
Authors:
Mahdi Naseri,
Gabriele Bozzola,
Vasileios Paschalidis
Abstract:
A viable model for the dense matter equation of state above the nuclear saturation density is a hadron-to-quark phase transition at densities relevant to compact objects. In this case, stable hybrid hadron-quark stars can arise. An even more interesting scenario is one where the hadron-to-quark phase transition results in the emergence of a third branch of stable compact objects (in addition to wh…
▽ More
A viable model for the dense matter equation of state above the nuclear saturation density is a hadron-to-quark phase transition at densities relevant to compact objects. In this case, stable hybrid hadron-quark stars can arise. An even more interesting scenario is one where the hadron-to-quark phase transition results in the emergence of a third branch of stable compact objects (in addition to white dwarfs and neutron stars). Inherent to the presence of a third family of compact stars is the existence of twin stars -- hybrid stars with the same mass as the corresponding neutron stars, but with smaller radii. Interestingly, the neutron star-twin star scenario is consistent with GW170817. If twin stars exist in nature, it raises a question about the mechanism that leads to their formation. Here, we explore gravitational collapse as a pathway to the formation of low-mass twin stars. We perform fully general relativistic simulations of the collapse of a stellar iron core, modeled as a cold degenerate gas, to investigate whether the end product is a neutron star or a twin star. Our simulations show that even with unrealistically large perturbations in the initial conditions, the core bounces well below the hadron-to-quark phase transition density, if the initial total rest mass is in the twin star range. Following cooling, these configurations produce neutron stars. We find that twin stars can potentially form due to mass loss, e.g., through winds, from a slightly more massive hybrid star that was initially produced in the collapse of a more massive core or if the maximum NS mass is below the Chandrasekhar mass limit. The challenge in producing twin stars in gravitational collapse, in conjunction with the fine tuning required because of their narrow mass range, suggest the rarity of twin stars in nature.
△ Less
Submitted 21 June, 2024;
originally announced June 2024.
-
Waveform Modelling for the Laser Interferometer Space Antenna
Authors:
LISA Consortium Waveform Working Group,
Niayesh Afshordi,
Sarp Akçay,
Pau Amaro Seoane,
Andrea Antonelli,
Josu C. Aurrekoetxea,
Leor Barack,
Enrico Barausse,
Robert Benkel,
Laura Bernard,
Sebastiano Bernuzzi,
Emanuele Berti,
Matteo Bonetti,
Béatrice Bonga,
Gabriele Bozzola,
Richard Brito,
Alessandra Buonanno,
Alejandro Cárdenas-Avendaño,
Marc Casals,
David F. Chernoff,
Alvin J. K. Chua,
Katy Clough,
Marta Colleoni,
Mekhi Dhesi,
Adrien Druart
, et al. (121 additional authors not shown)
Abstract:
LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmologic…
▽ More
LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome.
△ Less
Submitted 20 December, 2023; v1 submitted 2 November, 2023;
originally announced November 2023.
-
Not all spacetime coordinates for general-relativistic ray tracing are created equal
Authors:
Gabriele Bozzola,
Chi-kwan Chan,
Vasileios Paschalidis
Abstract:
Models for the observational appearance of astrophysical black holes rely critically on accurate general-relativistic ray tracing and radiation transport to compute the intensity measured by a distant observer. In this paper, we illustrate how the choice of coordinates and initial conditions affect this process. In particular, we show that propagating rays from the camera to the source leads to di…
▽ More
Models for the observational appearance of astrophysical black holes rely critically on accurate general-relativistic ray tracing and radiation transport to compute the intensity measured by a distant observer. In this paper, we illustrate how the choice of coordinates and initial conditions affect this process. In particular, we show that propagating rays from the camera to the source leads to different solutions if the spatial part of the momentum of the photon points towards the horizon or away from it. In doing this, we also show that coordinates that are well suited for numerical General-Relativistic MagnetoHydroDynamic (GRMHD) simulations are typically not optimal for generic ray tracing. We discuss the implications for black-hole images and show that radiation transport in optimal and non-optimal spacetime coordinates lead to the same images up to numerical errors and algorithmic choices.
△ Less
Submitted 3 October, 2023;
originally announced October 2023.
-
Quantifying uncertainties in general relativistic magnetohydrodynamic codes
Authors:
Pedro Luis Espino,
Gabriele Bozzola,
Vasileios Paschalidis
Abstract:
In this paper, we show that similar open-source codes for general relativistic magnetohydrodynamic (GRMHD) produce different results for key features of binary neutron star mergers. First, we present a new open-source version of the publicly available IllinoisGRMHD code that provides support for realistic, finite temperature equations of state. After stringent tests of our upgraded code, we perfor…
▽ More
In this paper, we show that similar open-source codes for general relativistic magnetohydrodynamic (GRMHD) produce different results for key features of binary neutron star mergers. First, we present a new open-source version of the publicly available IllinoisGRMHD code that provides support for realistic, finite temperature equations of state. After stringent tests of our upgraded code, we perform a code comparison between GRHydro, IllinoisGRMHD, Spritz, and WhiskyTHC, which implement the same physics, but slightly different computational methods. The benefit of the comparison is that all codes are embedded in the EinsteinToolkit suite, hence their only difference is algorithmic. We find similar convergence properties, fluid dynamics, and gravitational waves, but different merger times, remnant lifetimes, and gravitational wave phases. Such differences must be resolved before the post-merger dynamics modeled with such simulations can be reliably used to infer the properties of nuclear matter especially in the era of precision gravitational wave astronomy.
△ Less
Submitted 26 October, 2022; v1 submitted 24 October, 2022;
originally announced October 2022.
-
Kicks in charged black hole binaries
Authors:
Raimon Luna,
Gabriele Bozzola,
Vitor Cardoso,
Vasileios Paschalidis,
Miguel Zilhão
Abstract:
We compute the emission of linear momentum (kicks) by both gravitational and electromagnetic radiation in fully general-relativistic numerical evolutions of quasi-circular charged black hole binaries. We derive analytical expressions for slowly moving bodies and explore numerically a variety of mass ratios and charge-to-mass ratios. We find that for the equal mass case our analytical expression is…
▽ More
We compute the emission of linear momentum (kicks) by both gravitational and electromagnetic radiation in fully general-relativistic numerical evolutions of quasi-circular charged black hole binaries. We derive analytical expressions for slowly moving bodies and explore numerically a variety of mass ratios and charge-to-mass ratios. We find that for the equal mass case our analytical expression is in excellent agreement with the observed values and, contrarily to what happens in the vacuum case, we find that in presence of electromagnetic fields there is emission of momentum by gravitational waves. We also find that the strong gravitational kicks of binaries with unequal masses affect the electromagnetic kicks, causing them to strongly deviate from Keplerian predictions. For the values of charge-to-mass ratio considered in this work, we observe that magnitudes of the electromagnetic kicks are always smaller than the gravitational ones.
△ Less
Submitted 9 November, 2022; v1 submitted 13 July, 2022;
originally announced July 2022.
-
Numerical-relativity simulations of the quasi-circular inspiral and merger of non-spinning, charged black holes: methods and comparison with approximate approaches
Authors:
Gabriele Bozzola,
Vasileios Paschalidis
Abstract:
We present fully general relativistic simulations of the quasi-circular inspiral and merger of charged, non-spinning, binary black holes with charge-to-mass ratio $λ\le 0.3$. We discuss the key features that enabled long term and stable evolutions of these binaries. We also present a formalism for computing the angular momentum carried away by electromagnetic waves, and the electromagnetic contrib…
▽ More
We present fully general relativistic simulations of the quasi-circular inspiral and merger of charged, non-spinning, binary black holes with charge-to-mass ratio $λ\le 0.3$. We discuss the key features that enabled long term and stable evolutions of these binaries. We also present a formalism for computing the angular momentum carried away by electromagnetic waves, and the electromagnetic contribution to black-hole horizon properties. We implement our formalism and present the results for the first time in numerical-relativity simulations. In addition, we compare our full non-linear solutions with existing approximate models for the inspiral and ringdown phases. We show that Newtonian models based on the quadrupole approximation have errors of 20 % - 100 % in key gauge-invariant quantities. On the other hand, for the systems considered, we find that estimates of the remnant black hole spin based on the motion of test particles in Kerr-Newman spacetimes agree with our non-linear calculations to within a few percent. Finally, we discuss the prospects for detecting black hole charge by future gravitational-wave detectors using either the inspiral-merger-ringdown signal or the ringdown signal alone.
△ Less
Submitted 11 August, 2022; v1 submitted 14 April, 2021;
originally announced April 2021.
-
kuibit: Analyzing Einstein Toolkit simulations with Python
Authors:
Gabriele Bozzola
Abstract:
In the era of gravitational-wave astronomy, general-relativistic simulations of compact objects play a role of paramount importance. These calculations can be performed with the Einstein Toolkit, an open-source and community-supported software for numerical-relativity and relativistic astrophysics. The code comes with multiple solvers for Einstein's equations and for the equations of general-relat…
▽ More
In the era of gravitational-wave astronomy, general-relativistic simulations of compact objects play a role of paramount importance. These calculations can be performed with the Einstein Toolkit, an open-source and community-supported software for numerical-relativity and relativistic astrophysics. The code comes with multiple solvers for Einstein's equations and for the equations of general-relativistic magneto-hydrodynamics, along with a series of useful diagnostics. However, analyzing the output of the Einstein Toolkit can be a challenging task. Usually, the process involves a series of technical obstacles, like combining data from different restarts or working with HDF5 files. Here, we present kuibit, a Python library that takes care of all these low-level details (and many other more) and that provides high-level, intuitive, representations of the data. Kuibit ships with a wide range of features that include full support for 1-3D ASCII and HDF5 grid data, time and frequency series, gravitational waves, and apparent horizons. With kuibit, users can inspect most of the content of a simulation with just a few lines of code. Importantly, kuibit is designed to be a code for the community: it is user-friendly and does not require any proprietary software to run, it has documentation and examples, and it is openly developed with emphasis on extensibility and maintainability.
△ Less
Submitted 13 April, 2021;
originally announced April 2021.
-
General relativistic simulations of the quasi-circular inspiral and merger of charged black holes: GW150914 and fundamental physics implications
Authors:
Gabriele Bozzola,
Vasileios Paschalidis
Abstract:
We perform general-relativistic simulations of charged black holes targeting GW150914. We show that the inspiral is most efficient for detecting black hole charge through gravitational waves and that GW150914 is compatible with having charge-to-mass ratio as high as 0.3. Our work applies to electric and magnetic charge, and to theories with black holes endowed with U(1) (hidden or dark) charges. U…
▽ More
We perform general-relativistic simulations of charged black holes targeting GW150914. We show that the inspiral is most efficient for detecting black hole charge through gravitational waves and that GW150914 is compatible with having charge-to-mass ratio as high as 0.3. Our work applies to electric and magnetic charge, and to theories with black holes endowed with U(1) (hidden or dark) charges. Using our results we place an upper bound on the deviation from general relativity in the dynamical, strong-field regime of the so-called theory of MOdified Gravity (MOG).
△ Less
Submitted 27 January, 2021; v1 submitted 28 June, 2020;
originally announced June 2020.
-
Maximum mass and universal relations of rotating relativistic hybrid hadron-quark stars
Authors:
Gabriele Bozzola,
Pedro Luis Espino,
Collin Davis Lewin,
Vasileios Paschalidis
Abstract:
We construct equilibrium models of uniformly and differentially rotating hybrid hadron-quark stars using equations of state (EOSs) with a first-order phase transition that gives rise to a third family of compact objects. We find that the ratio of the maximum possible mass of uniformly rotating configurations - the supramassive limit - to the Tolman-Oppenheimer-Volkoff (TOV) limit mass is not EOS-i…
▽ More
We construct equilibrium models of uniformly and differentially rotating hybrid hadron-quark stars using equations of state (EOSs) with a first-order phase transition that gives rise to a third family of compact objects. We find that the ratio of the maximum possible mass of uniformly rotating configurations - the supramassive limit - to the Tolman-Oppenheimer-Volkoff (TOV) limit mass is not EOS-independent, and is between 1.15 and 1.31,in contrast with the value of 1.20 previously found for hadronic EOSs. Therefore, some of the constraints placed on the EOS from the observation of the gravitational wave event GW170817 do not apply to hadron-quark EOSs. However, the supramassive limit mass for the family of EOSs we treat is consistent with limits set by GW170817, strengthening the possibility of interpreting GW170817 with a hybrid hadron-quark EOSs. We also find that along constant angular momentum sequences of uniformly rotating stars, the third family maximum and minimum mass models satisfy approximate EOS-independent relations, and the supramassive limit of the third family is approximately 16.5 % larger than the third family TOV limit. For differentially rotating spheroidal stars, we find that a lower-limit on the maximum supportable rest mass is 123 % more than the TOV limit rest mass. Finally, we verify that the recently discovered universal relations relating angular momentum, rest mass and gravitational mass for turning-point models hold for hybrid hadron-quark EOSs when uniform rotation is considered, but have a clear dependence on the degree of differential rotation.
△ Less
Submitted 30 April, 2019;
originally announced May 2019.
-
Initial data for general relativistic simulations of multiple electrically charged black holes with linear and angular momenta
Authors:
Gabriele Bozzola,
Vasileios Paschalidis
Abstract:
A general relativistic, stationary and axisymmetric black hole in a four-dimensional asymptotically-flat spacetime is fully determined by its mass, angular momentum and electric charge. The expectation that astrophysically relevant black holes do not posses charge has resulted in a limited number of investigations of moving and charged black holes in the dynamical, strong-field gravitational (and…
▽ More
A general relativistic, stationary and axisymmetric black hole in a four-dimensional asymptotically-flat spacetime is fully determined by its mass, angular momentum and electric charge. The expectation that astrophysically relevant black holes do not posses charge has resulted in a limited number of investigations of moving and charged black holes in the dynamical, strong-field gravitational (and electromagnetic) regime, where numerical studies are necessary. Apart from having a theoretical interest, the advent of multimessenger astronomy with gravitational waves offers new ways to think about charged black holes. In this work, we initiate an exploration of charged binary black holes by generating valid initial data for general relativistic simulations of black hole systems that have generic electric charge, linear and angular momenta. We develop our initial data formalism within the framework of the conformal transverse-traceless (Bowen-York) technique using the puncture approach, and apply the theory of isolated horizons to attribute physical parameters (mass, charge and angular momentum) to each hole. We implemented our formalism in the case of a binary system by modifying the publicly available TwoPunctures and QuasiLocalMeasures codes. We demonstrate that our code can recover existing solutions and that it has excellent self-convergence properties for a generic configuration of two black holes.
△ Less
Submitted 18 June, 2019; v1 submitted 3 March, 2019;
originally announced March 2019.
-
Universal relations for differentially rotating relativistic stars at the threshold to collapse
Authors:
Gabriele Bozzola,
Nikolaos Stergioulas,
Andreas Bauswein
Abstract:
A binary neutron star merger produces a rapidly and differentially rotating compact remnant whose lifespan heavily affects the electromagnetic and gravitational emissions. Its stability depends on both the equation of state (EOS) and the rotation law and it is usually investigated through numerical simulations. Nevertheless, by means of a sufficient criterion for secular instability, equilibrium s…
▽ More
A binary neutron star merger produces a rapidly and differentially rotating compact remnant whose lifespan heavily affects the electromagnetic and gravitational emissions. Its stability depends on both the equation of state (EOS) and the rotation law and it is usually investigated through numerical simulations. Nevertheless, by means of a sufficient criterion for secular instability, equilibrium sequences can be used as a computationally inexpensive way to estimate the onset of dynamical instability, which, in general, is close to the secular one. This method works well for uniform rotation and relies on the location of turning points: stellar models that are stationary points in a sequence of equilibrium solutions with constant rest mass or angular momentum. Here, we investigate differentially rotating models (using a large number of equations of state and different rotation laws) and find that several universal relations between properly scaled gravitational mass, rest mass and angular momentum of the turning-point models that are valid for uniform rotation, are insensitive to the degree of differential rotation, to high accuracy.
△ Less
Submitted 8 September, 2017;
originally announced September 2017.