-
Constraining the nuclear equation of state from orbits of primordial black holes inside neutron stars
Authors:
Anna Chen,
Juan Diego DelPrado,
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
Lacking terrestrial experimental data, our best constraints on the behavior of matter at high densities up to and above nuclear density arise from observations of neutron stars. Current constraints include those based on measurements of stellar masses, radii, and tidal deformabilities. Here we explore how orbits of primordial black holes - should they exist - inside neutron stars could provide com…
▽ More
Lacking terrestrial experimental data, our best constraints on the behavior of matter at high densities up to and above nuclear density arise from observations of neutron stars. Current constraints include those based on measurements of stellar masses, radii, and tidal deformabilities. Here we explore how orbits of primordial black holes - should they exist - inside neutron stars could provide complementary constraints on the nuclear equation of state (EOS). Specifically, we consider a sample of candidate EOSs, construct neutron star models for these EOSs, and compute orbits of primordial black holes inside these stars. We discuss how the pericenter advance of eccentric orbits, i.e. orbital precession, results in beat phenomena in the emitted gravitational wave signal. Observing this beat frequency could constrain the nuclear EOS and break possible degeneracies arising from other constraints, as well as provide information about the host star.
△ Less
Submitted 6 September, 2024;
originally announced September 2024.
-
Primordial black holes captured by neutron stars: simulations in general relativity
Authors:
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
We present self-consistent numerical simulations in general relativity of putative primordial black holes inside neutron stars. Complementing a companion paper in which we assumed the black hole mass $m$ to be much smaller than the mass $M_*$ of the neutron star, thereby justifying a point-mass treatment, we here consider black holes with masses large enough so that their effect on the neutron sta…
▽ More
We present self-consistent numerical simulations in general relativity of putative primordial black holes inside neutron stars. Complementing a companion paper in which we assumed the black hole mass $m$ to be much smaller than the mass $M_*$ of the neutron star, thereby justifying a point-mass treatment, we here consider black holes with masses large enough so that their effect on the neutron star cannot be neglected. We develop and employ several new numerical techniques, including initial data describing boosted black holes in neutron-star spacetimes, a relativistic determination of the escape speed, and a gauge condition that keeps the black hole hole at a fixed coordinate location. We then perform numerical simulations that highlight different aspects of the capture of primordial black holes by neutron stars. In particular, we simulate the initial passage of the black hole through the star, demonstrating that the neutron star remains dynamically stable provided the black-hole mass is sufficiently small, $m \lesssim 0.05 M_*$. We also model the late evolution of a black hole oscillating about the center of an initially stable neutron star while accreting stellar mass and ultimately triggering gravitational collapse.
△ Less
Submitted 16 May, 2024;
originally announced May 2024.
-
Jetlike structures in low-mass binary neutron star merger remnants
Authors:
Jamie Bamber,
Antonios Tsokaros,
Milton Ruiz,
Stuart L. Shapiro
Abstract:
GW170817 and GRB 170817A provided direct evidence that binary neutron star (NSNS) mergers can produce short gamma-ray bursts (sGRBs). However, questions remain about the nature of the central engine. Depending on the mass, the remnant from a NSNS merger may promptly collapse to a black hole (BH), form a hypermassive neutron star (HMNS) which undergoes a delayed collapse to a BH, a supramassive neu…
▽ More
GW170817 and GRB 170817A provided direct evidence that binary neutron star (NSNS) mergers can produce short gamma-ray bursts (sGRBs). However, questions remain about the nature of the central engine. Depending on the mass, the remnant from a NSNS merger may promptly collapse to a black hole (BH), form a hypermassive neutron star (HMNS) which undergoes a delayed collapse to a BH, a supramassive neutron star (SMNS) with a much longer lifetime, or an indefinitely stable NS. There is strong evidence that a BH with an accretion disk can launch a sGRB-compatible jet via the Blandford-Znajek mechanism, but whether a supramassive star can do the same is less clear. We have performed general relativistic magnetohydrodynamics simulations of the merger of both irrotational and spinning, equal-mass NSNSs constructed from a piecewise polytropic representation of the SLy equation of state, with a range of gravitational masses that yield remnants with mass above and below the supramassive limit. Each NS is endowed with a dipolar magnetic field extending from the interior into the exterior, as in a radio pulsar. We examine cases with different initial binary masses, including a case which produces a HMNS which collapses to a BH, and lower mass binaries that produce SMNS remnants. We find similar jetlike structures for both the SMNS and HMNS remnants that meet our basic critera for an incipient jet. The outflow for the HMNS case is consistent with a Blandford-Znajek (BZ) jet. There is sufficient evidence that such BZ-powered outflows can break out and produce ulrarelativistic jets so that we can describe the HMNS system as a sGRB progenitor. However, the incipient jets from the SMNS remnants have much more baryon pollution and we see indications of inefficient acceleration and mixing with the surrounding debris. Therefore, we cannot conclude that SMNS outflows are the progenitors of sGRBs.
△ Less
Submitted 24 August, 2024; v1 submitted 3 May, 2024;
originally announced May 2024.
-
Primordial black holes captured by neutron stars: relativistic point-mass treatment
Authors:
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
Primordial black holes (PBHs), if they exist, may collide with and be captured by neutron stars. We adopt a relativistic point-mass approximation to study this capture, the subsequent confinement of the PBH of mass $m$ inside the neutron star of mass $M_* \gg m$, and the PBH's growth by accretion of stellar material. Building on earlier treatments we systematically study the capture, confinement,…
▽ More
Primordial black holes (PBHs), if they exist, may collide with and be captured by neutron stars. We adopt a relativistic point-mass approximation to study this capture, the subsequent confinement of the PBH of mass $m$ inside the neutron star of mass $M_* \gg m$, and the PBH's growth by accretion of stellar material. Building on earlier treatments we systematically study the capture, confinement, and accretion process, characterize the emitted quasiperiodic continuous gravitational-wave signal, track the evolution of the PBH's orbital parameters, and compare the effects of different choices for the prescription of the dissipative forces. Our point-mass treatment here is applicable in the limit of small PBH masses, for which its effects on the neutron star can be ignored.
△ Less
Submitted 24 June, 2024; v1 submitted 12 April, 2024;
originally announced April 2024.
-
Could long-period transients be powered by primordial black hole capture?
Authors:
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
Long-period radio transients have unusual properties that challenge their interpretation as pulsars or magnetars. We examine whether they might instead be powered by primordial black holes (PBHs) making repeated passages through a host star, thereby providing a signature of elusive dark-matter candidates. We demonstrate that constraints derived from the transients' period and period derivative alo…
▽ More
Long-period radio transients have unusual properties that challenge their interpretation as pulsars or magnetars. We examine whether they might instead be powered by primordial black holes (PBHs) making repeated passages through a host star, thereby providing a signature of elusive dark-matter candidates. We demonstrate that constraints derived from the transients' period and period derivative alone already rule out this scenario for most potential host stars. While white dwarfs may satisfy these constraints, they are unlikely to capture PBHs in the required mass range.
△ Less
Submitted 16 February, 2024;
originally announced February 2024.
-
Primordial black hole capture, gravitational wave beats, and the nuclear equation of state
Authors:
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
Primordial black holes (PBHs), if captured by neutron stars (NSs), would emit a characteristic gravitational wave (GW) signal as they orbit inside the host star. We identify a specific and qualitatively new feature of these signals, namely quasi-periodic beats caused by the precession of noncircular PBH orbits. We demonstrate numerically and analytically that the beat frequency depends rather sens…
▽ More
Primordial black holes (PBHs), if captured by neutron stars (NSs), would emit a characteristic gravitational wave (GW) signal as they orbit inside the host star. We identify a specific and qualitatively new feature of these signals, namely quasi-periodic beats caused by the precession of noncircular PBH orbits. We demonstrate numerically and analytically that the beat frequency depends rather sensitively on the NS structure, so that hypothetical future observations with next-generation GW detectors would provide valuable constraints on the nuclear equation of state.
△ Less
Submitted 9 July, 2024; v1 submitted 2 February, 2024;
originally announced February 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.
-
Spikes and accretion of unbound, collisionless matter around black holes
Authors:
Stuart L. Shapiro
Abstract:
We consider the steady-state density and velocity dispersion profiles of collisionless matter around a Schwarzschild black hole (BH) and its associated rate of accretion onto the BH. We treat matter, which could be stars or dark matter particles, whose orbits are {\it unbound} to the BH, but still governed by its gravitational field. We consider two opposite spatial geometries for the matter distr…
▽ More
We consider the steady-state density and velocity dispersion profiles of collisionless matter around a Schwarzschild black hole (BH) and its associated rate of accretion onto the BH. We treat matter, which could be stars or dark matter particles, whose orbits are {\it unbound} to the BH, but still governed by its gravitational field. We consider two opposite spatial geometries for the matter distributions: an infinite, 3D cluster and a 2D razor-thin disk, both with zero net angular momentum. We demonstrate that the results depend critically on the adopted geometry, even in the absence of angular momentum. We adopt a simple monoenergetic, isotropic, phase-space distribution function for the matter as a convenient example to illustrate this dependence. The effect of breaking strict isotropy by incorporating an unreplenished loss cone due to BH capture of low-angular momentum matter is also considered. Calculations are all analytic and performed in full general relativity, though key results are also evaluated in the Newtonian limit. We consider one application to show that the rate of BH accretion from an ambient cluster can be significantly less than that from a thin disk to which it may collapse, although both rates are considerably smaller than Bondi accretion for a (collisional) fluid with a similar asymptotic particle density and velocity dispersion (i.e., sound speed).
△ Less
Submitted 20 October, 2023;
originally announced October 2023.
-
Effect of magnetic fields on the dynamics and gravitational wave emission of PPI-saturated self-gravitating accretion disks: simulations in full GR
Authors:
Erik Wessel,
Vasileios Paschalidis,
Antonios Tsokaros,
Milton Ruiz,
Stuart L. Shapiro
Abstract:
We explore the effect magnetic fields have on self-gravitating accretion disks around spinning black holes via numerical evolutions in full dynamical magnetohydrodynamic spacetimes. The configurations we study are unstable to the Papaloizou-Pringle Instability (PPI). PPI-saturated accretion tori have been shown to produce gravitational waves, detectable to cosmological distances by third-generatio…
▽ More
We explore the effect magnetic fields have on self-gravitating accretion disks around spinning black holes via numerical evolutions in full dynamical magnetohydrodynamic spacetimes. The configurations we study are unstable to the Papaloizou-Pringle Instability (PPI). PPI-saturated accretion tori have been shown to produce gravitational waves, detectable to cosmological distances by third-generation gravitational wave (GW) observatories. While the PPI operates strongly for purely hydrodynamic disks, the situation can be different for disks hosting initially small magnetic fields. Evolutions of disks without self-gravity in fixed BH spacetimes have shown that small seed fields can initiate the rapid growth of the magneto-rotational instability (MRI), which then strongly suppresses the PPI. Since realistic astrophysical disks are expected to be magnetized, PPI-generated GW signals may be suppressed as well. However, it is unclear what happens when the disk self-gravity is restored. Here, we study the impact of magnetic fields on the PPI-saturated state of a self-gravitating accretion disk around a spinning BH ($χ= 0.7$) aligned with the disk angular momentum, as well as one around a non-spinning BH. We find the MRI is effective at reducing the amplitude of PPI modes and their associated GWs, but the systems still generate GWs. Estimating the detectability of these systems accross a wide range of masses, we show that magnetic fields reduce the maximum detection distance by Cosmic Explorer from 300Mpc (in the pure hydrodynamic case) to 45Mpc for a $10 M_{\odot}$ system, by LISA from 11500Mpc to 2700Mpc for a $2 \times 10^{5} M_{\odot}$ system, and by DECIGO from $z \approx 5$ down to $z \approx 2$ for a $1000 M_{\odot}$ system.
△ Less
Submitted 26 June, 2023; v1 submitted 14 April, 2023;
originally announced April 2023.
-
General Relativistic Magnetohydrodynamic Simulations of Accretion Disks Around Tilted Binary Black Holes of Unequal Mass
Authors:
Milton Ruiz,
Antonios Tsokaros,
Stuart L. Shapiro
Abstract:
We perform general relativistic simulations of magnetized, accreting disks onto spinning binary black holes (BHBHs) with different mass ratios (MRs). The magnitude of the individual BH spins are all $χ= 0.26$ and lie either along the initial orbital plane or $45^\circ$ above it. We evolve these systems throughout the inspiral, merger and postmerger phases to identify the impact of the BH spins and…
▽ More
We perform general relativistic simulations of magnetized, accreting disks onto spinning binary black holes (BHBHs) with different mass ratios (MRs). The magnitude of the individual BH spins are all $χ= 0.26$ and lie either along the initial orbital plane or $45^\circ$ above it. We evolve these systems throughout the inspiral, merger and postmerger phases to identify the impact of the BH spins and the MR on any jet and their electromagnetic (EM) signatures. We find that incipient jets are launched from both BHs regardless of the MR and along the spin directions as long as the force-free parameter $B^2/(8\,πρ_0)$ in the funnel and above their poles is larger than one. At large distances the two jets merge into a single one what may prevent the EM detection of individual jets. As the accretion rate reaches a quasistationary state during predecoupling, we observe a sudden amplification of the outgoing Poynting luminosity that depends on the MR. Following merger, the sudden change in the direction of the spin of the BH remnant with respect to the spins of its progenitors causes a reorientation of the jet, which drives a single, high-velocity, outward narrow beam collimated by a tightly wound, helical Bfield which, in turn, boosts the Poynting luminosity. This effect is nearly MR independent. During this process, a kink is produced in the Bfield lines, confining the jet. The kink propagates along the jet but rapidly decays, leaving no memory of the spin-shift. Our results suggest that the merger of misaligned, low-spinning, BHBH mergers in low-mass disks may not provide a viable scenario to explain X-shaped radio galaxies if other features are not taken into account. However, the sudden changes in the outgoing luminosity at merger may help to identify mergers of BHs in active galactic nuclei, shedding light on BH growth mechanisms and the observed co-evolution of their host galaxies.
△ Less
Submitted 17 February, 2023;
originally announced February 2023.
-
The Effect of Stars on the Dark Matter Spike Around a Black Hole: A Tale of Two Treatments
Authors:
Stuart L. Shapiro,
Douglas C. Heggie
Abstract:
We revisit the role that gravitational scattering off stars plays in establishing the steady-state distribution of collisionless dark matter (DM) around a massive black hole (BH). This is a physically interesting problem that has potentially observable signatures, such as $γ-$rays from DM annihilation in a density spike. The system serves as a laboratory for comparing two different dynamical appro…
▽ More
We revisit the role that gravitational scattering off stars plays in establishing the steady-state distribution of collisionless dark matter (DM) around a massive black hole (BH). This is a physically interesting problem that has potentially observable signatures, such as $γ-$rays from DM annihilation in a density spike. The system serves as a laboratory for comparing two different dynamical approaches, both of which have been widely used: a Fokker-Planck treatment and a two-component conduction fluid treatment. In our Fokker-Planck analysis we extend a previous analytic model to account for a nonzero flux of DM particles into the BH, as well as a cut-off in the distribution function near the BH due to relativistic effects or, further out, possible DM annihilation. In our two-fluid analysis, following an approximate analytic treatment, we recast the equations as a "heated Bondi accretion" problem and solve the equations numerically without approximation. While both the Fokker-Planck and two-fluid methods yield basically the same DM density and velocity dispersion profiles away from the boundaries in the spike interior, there are other differences, especially the determination of the DM accretion rate. We discuss limitations of the two treatments, including the assumption of an isotropic velocity dispersion.
△ Less
Submitted 16 September, 2022;
originally announced September 2022.
-
Self-gravitating disks around rapidly spinning, tilted black holes: General relativistic simulations
Authors:
Antonios Tsokaros,
Milton Ruiz,
Stuart L. Shapiro,
Vasileios Paschalidis
Abstract:
We perform general relativistic simulations of self-gravitating black hole-disks in which the spin of the black hole is significantly tilted ($45^\circ$ and $90^\circ$) with respect to the angular momentum of the disk and the disk-to-black hole mass ratio is $16\%-28\%$. The black holes are rapidly spinning with dimensionless spins up to $\sim 0.97$. These are the first self-consistent hydrodynami…
▽ More
We perform general relativistic simulations of self-gravitating black hole-disks in which the spin of the black hole is significantly tilted ($45^\circ$ and $90^\circ$) with respect to the angular momentum of the disk and the disk-to-black hole mass ratio is $16\%-28\%$. The black holes are rapidly spinning with dimensionless spins up to $\sim 0.97$. These are the first self-consistent hydrodynamic simulations of such systems, which can be prime sources for multimessenger astronomy. In particular tilted black hole-disk systems lead to: i) black hole precession; ii) disk precession and warping around the black hole; iii) earlier saturation of the Papaloizou-Pringle instability compared to aligned/antialigned systems, although with a shorter mode growth timescale; iv) acquisition of a small black-hole kick velocity; v) significant gravitational wave emission via various modes beyond, but as strong as, the typical $(2,2)$ mode; and vi) the possibility of a broad alignment of the angular momentum of the disk with the black hole spin. This alignment is not related to the Bardeen-Petterson effect and resembles a solid body rotation. Our simulations suggest that any electromagnetic luminosity from our models may power relativistic jets, such as those characterizing short gamma-ray bursts. Depending on the black hole-disk system scale the gravitational waves may be detected by LIGO/Virgo, LISA and/or other laser interferometers.
△ Less
Submitted 9 September, 2022;
originally announced September 2022.
-
Jet Launching from Binary Neutron Star Mergers: Incorporating Neutrino Transport and Magnetic Fields
Authors:
Lunan Sun,
Milton Ruiz,
Stuart L. Shapiro,
Antonios Tsokaros
Abstract:
We perform general relativistic, magnetohydrodynamic (GRMHD) simulations of merging binary neutron stars incorporating neutrino transport and magnetic fields. Our new radiative transport module for neutrinos adopts a general relativistic, truncated-moment (M1) formalism. The binaries consist of two identical, irrotational stars modeled by the SLy nuclear equation of state (EOS). They are initially…
▽ More
We perform general relativistic, magnetohydrodynamic (GRMHD) simulations of merging binary neutron stars incorporating neutrino transport and magnetic fields. Our new radiative transport module for neutrinos adopts a general relativistic, truncated-moment (M1) formalism. The binaries consist of two identical, irrotational stars modeled by the SLy nuclear equation of state (EOS). They are initially in quasicircular orbit and threaded with a poloidal magnetic field that extends from the stellar interior into the exterior, as in typical pulsars. We insert neutrino processes shortly after the merger and focus on the role of neutrinos in launching a jet following the collapse of the hypermassive neutron star (HMNS) remnant to a spinning black hole (BH). We treat two microphysical versions: one (a "warm-up") evolving a single neutrino species and considering only charged-current processes, and the other evolving three species $(ν_e, \barν_e, ν_{\rm x})$ and related processes. We trace the evolution until the system reaches a quasiequilibrium state after BH formation. We find that the BH + disk remnant eventually launches an incipient jet. The electromagnetic Poynting luminosity is $\sim 10^{53} \rm \, erg\, s^{-1}$, consistent with that of typical short gamma-ray bursts (sGRBs). The effect of neutrino cooling shortens the lifetime of the HMNS, and lowers the amplitude of the major peak of the gravitational wave (GW) power spectrum somewhat. After BH formation, neutrinos help clear out the matter near the BH poles, resulting in lower baryon-loaded surrounding debris. The neutrino luminosity resides in the range $\sim 10^{52-53} \rm \,erg\,s^{-1}$ once quasiequilibrium is achieved. Comparing with the neutrino-free models, we observe that the inclusion of neutrinos yields similar ejecta masses and is inefficient in carrying off additional angular momentum.
△ Less
Submitted 16 May, 2022; v1 submitted 25 February, 2022;
originally announced February 2022.
-
Did a close tidal encounter cause the Great Dimming of Betelgeuse?
Authors:
Hailey Aronson,
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
We assess whether gravity darkening, induced by a tidal interaction during a stellar fly-by, might be sufficient to explain the Great Dimming of Betelgeuse. Adopting several simple approximations, we calculate the tidal deformation and associated gravity darkening in a close tidal encounter, as well as the reduction in the radiation flux as seen by a distant observer. We show that, in principle, t…
▽ More
We assess whether gravity darkening, induced by a tidal interaction during a stellar fly-by, might be sufficient to explain the Great Dimming of Betelgeuse. Adopting several simple approximations, we calculate the tidal deformation and associated gravity darkening in a close tidal encounter, as well as the reduction in the radiation flux as seen by a distant observer. We show that, in principle, the duration and degree of the resulting stellar dimming can be used to estimate the minimum pericenter separation and mass of a fly-by object, which, even if it remains undetected otherwise, might be a black hole, neutron star, or white dwarf. Our estimates show that, while such fly-by events may occur in other astrophysical scenarios, where our analysis should be applicable, they likely are not large enough to explain the Great Dimming of Betelgeuse by themselves.
△ Less
Submitted 21 October, 2022; v1 submitted 20 January, 2022;
originally announced January 2022.
-
Magnetohydrodynamic simulations of self-consistent rotating neutron stars with mixed poloidal and toroidal magnetic fields
Authors:
Antonios Tsokaros,
Milton Ruiz,
Stuart L. Shapiro,
Kōji Uryū
Abstract:
We perform the first magnetohydrodynamic simulations in full general relativity of self-consistent rotating neutron stars (NSs) with ultrastrong mixed poloidal and toroidal magnetic fields. The initial uniformly rotating NS models are computed assuming perfect conductivity, stationarity, and axisymmetry. Although the specific geometry of the mixed field configuration can delay or accelerate the de…
▽ More
We perform the first magnetohydrodynamic simulations in full general relativity of self-consistent rotating neutron stars (NSs) with ultrastrong mixed poloidal and toroidal magnetic fields. The initial uniformly rotating NS models are computed assuming perfect conductivity, stationarity, and axisymmetry. Although the specific geometry of the mixed field configuration can delay or accelerate the development of various instabilities known from analytic perturbative studies, all our models finally succumb to them. Differential rotation is developed spontaneously in the cores of our magnetars which, after sufficient time, is converted back to uniform rotation. The rapidly rotating magnetars show a significant amount of ejecta, which can be responsible for transient kilonova signatures. However no highly collimated, helical magnetic fields or incipient jets, which are necessary for gamma-ray bursts, arise at the poles of these magnetars by the time our simulations are terminated.
△ Less
Submitted 29 October, 2021;
originally announced November 2021.
-
Jet Launching from Merging Magnetized Binary Neutron Stars with Realistic Equations of State
Authors:
Milton Ruiz,
Antonios Tsokaros,
Stuart L. Shapiro
Abstract:
We perform general relativistic, magnetohydrodynamic (GRMHD) simulations of binary neutron stars in quasi-circular orbit that merge and undergo delayed or prompt collapse to a black hole (BH). The stars are irrotational and modeled using an SLy or an H4 nuclear equation of state. To assess the impact of the initial magnetic field configuration on jet launching, we endow the stars with a purely pol…
▽ More
We perform general relativistic, magnetohydrodynamic (GRMHD) simulations of binary neutron stars in quasi-circular orbit that merge and undergo delayed or prompt collapse to a black hole (BH). The stars are irrotational and modeled using an SLy or an H4 nuclear equation of state. To assess the impact of the initial magnetic field configuration on jet launching, we endow the stars with a purely poloidal magnetic field that is initially unimportant dynamically and is either confined to the stellar interior or extends from the interior into the exterior as in typical pulsars. Consistent with our previous results, we find that only the BH + disk remnants originating from binaries that form hypermassive neutron stars (HMNSs) and undergo delayed collapse can drive magnetically-powered jets. We find that the closer the total mass of the binary is to the threshold value for prompt collapse, the shorter is the time delay between the gravitational wave peak amplitude and jet launching. This time delay also strongly depends on the initial magnetic field configuration. We also find that seed magnetic fields confined to the stellar interior can launch a jet over $\sim 25\,\rm ms$ later than those with pulsar-like magnetic fields. The lifetime of the jet [$Δt\lesssim 150\,\rm ms$] and its outgoing Poynting luminosity [$L_{\rm EM}\sim 10^{52\pm 1}\rm erg/s$] are consistent with typical short gamma-ray burst central engine lifetimes, as well as with the Blandford--Znajek mechanism for launching jets and their associated Poynting luminosities. Our numerical results also suggest that the dynamical ejection of matter can be enhanced by the magnetic field. Therefore, GRMHD studies are required to fully understand kilonova signals from GW170818-like events.
△ Less
Submitted 16 December, 2021; v1 submitted 22 October, 2021;
originally announced October 2021.
-
Accretion onto black holes inside neutron stars with piecewise-polytropic equations of state: analytic and numerical treatments
Authors:
Sophia C. Schnauck,
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
We consider spherically symmetric accretion onto a small, possibly primordial, black hole residing at the center of a neutron star governed by a cold nuclear equation of state (EOS). We generalize the relativistic Bondi solution for such EOSs, approximated by piecewise polytropes, and thereby obtain analytical expressions for the steady-state matter profiles and accretion rates. We compare these r…
▽ More
We consider spherically symmetric accretion onto a small, possibly primordial, black hole residing at the center of a neutron star governed by a cold nuclear equation of state (EOS). We generalize the relativistic Bondi solution for such EOSs, approximated by piecewise polytropes, and thereby obtain analytical expressions for the steady-state matter profiles and accretion rates. We compare these rates with those found by time-dependent, general relativistic hydrodynamical simulations upon relaxation and find excellent agreement. We consider several different candidate EOSs, neutron star masses and central densities and find that the accretion rates vary only little, resulting in an accretion rate that depends primarily on the black hole mass, and only weakly on the properties of the neutron star.
△ Less
Submitted 7 December, 2021; v1 submitted 15 October, 2021;
originally announced October 2021.
-
Accretion onto a small black hole at the center of a neutron star
Authors:
Chloe B. Richards,
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
We revisit the system consisting of a neutron star that harbors a small, possibly primordial, black hole at its center, focusing on a nonspinning black hole embedded in a nonrotating neutron star. Extending earlier treatments, we provide an analytical treatment describing the rate of secular accretion of the neutron star matter onto the black hole, adopting the relativistic Bondi accretion formali…
▽ More
We revisit the system consisting of a neutron star that harbors a small, possibly primordial, black hole at its center, focusing on a nonspinning black hole embedded in a nonrotating neutron star. Extending earlier treatments, we provide an analytical treatment describing the rate of secular accretion of the neutron star matter onto the black hole, adopting the relativistic Bondi accretion formalism for stiff equations of state that we presented elsewhere. We use these accretion rates to sketch the evolution of the system analytically until the neutron star is completely consumed. We also perform numerical simulations in full general relativity for black holes with masses up to nine orders of magnitude smaller than the neutron star mass, including a simulation of the entire evolution through collapse for the largest black hole mass. We construct relativistic initial data for these simulations by generalizing the black hole puncture method to allow for the presence of matter, and evolve these data with a code that is optimally designed to resolve the vastly different length scales present in this problem. We compare our analytic and numerical results, and provide expressions for the lifetime of neutron stars harboring such endoparasitic black holes.
△ Less
Submitted 18 May, 2021; v1 submitted 18 February, 2021;
originally announced February 2021.
-
Multimessenger Binary Mergers Containing Neutron Stars: Gravitational Waves, Jets, and $\boldsymbolγ$-Ray Bursts
Authors:
Milton Ruiz,
Stuart L. Shapiro,
Antonios Tsokaros
Abstract:
Neutron stars (NSs) are extraordinary not only because they are the densest form of matter in the visible Universe but also because they can generate B-fields ten orders of magnitude larger than those currently constructed on Earth. The combination of extreme gravity with the enormous electromagnetic (EM) fields gives rise to spectacular phenomena like those observed on August 2017 with the merger…
▽ More
Neutron stars (NSs) are extraordinary not only because they are the densest form of matter in the visible Universe but also because they can generate B-fields ten orders of magnitude larger than those currently constructed on Earth. The combination of extreme gravity with the enormous electromagnetic (EM) fields gives rise to spectacular phenomena like those observed on August 2017 with the merger of a binary neutron star (NSNS) system, an event that generated a gravitational wave (GW) signal, a short $γ$-ray burst (sGRB), and a kilonova. This event serves as the highlight so far of the era of multimessenger astronomy. In this review, we present the current state of our theoretical understanding of compact binary mergers containing NSs as gleaned from the latest general relativistic magnetohydrodynamic simulations. Such mergers can lead to events like the one on August 2017, GW170817, and its EM counterparts, GRB 170817 and AT 2017gfo. In addition to exploring the GW emission from binary black hole-neutron star and NSNS mergers, we also focus on their counterpart EM signals. In particular, we are interested in identifying the conditions under which a relativistic jet can be launched following these mergers. Such a jet is an essential feature of most sGRB models and provides the main conduit of energy from the central object to the outer radiation regions. Jet properties, including their lifetimes and Poynting luminosities, the effects of the initial B-field geometries and spins of the coalescing NSs, as well as their governing equation of state, are discussed. Lastly, we present our current understanding of how the Blandford-Znajek mechanism arises from merger remnants as the trigger for launching jets, if, when and how a horizon is necessary for this mechanism, and the possibility that it can turn on in magnetized neutron ergostars, which contain ergoregions, but no horizons.
△ Less
Submitted 8 April, 2021; v1 submitted 5 February, 2021;
originally announced February 2021.
-
Neutron Stars Harboring a Primordial Black Hole: Maximum Survival Time
Authors:
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
We explore in general relativity the survival time of neutron stars that host an endoparasitic, possibly primordial, black hole at their center. Corresponding to the minimum steady-state Bondi accretion rate for adiabatic flow that we found earlier for stiff nuclear equations of state (EOSs), we derive analytically the maximum survival time after which the entire star will be consumed by the black…
▽ More
We explore in general relativity the survival time of neutron stars that host an endoparasitic, possibly primordial, black hole at their center. Corresponding to the minimum steady-state Bondi accretion rate for adiabatic flow that we found earlier for stiff nuclear equations of state (EOSs), we derive analytically the maximum survival time after which the entire star will be consumed by the black hole. We also show that this maximum survival time depends only weakly on the stiffness for polytropic EOSs with $Γ\geq 5/3$, so that this survival time assumes a nearly universal value that depends on the initial black hole mass alone. Establishing such a value is important for constraining the contribution of primordial black holes in the mass range $10^{-16} M_\odot \lesssim M \lesssim 10^{-10} M_\odot$ to the dark-matter content of the Universe.
△ Less
Submitted 13 April, 2021; v1 submitted 28 January, 2021;
originally announced January 2021.
-
Relativistic Bondi accretion for stiff equations of state
Authors:
Chloe B. Richards,
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
We revisit Bondi accretion - steady-state, adiabatic, spherical gas flow onto a Schwarzschild black hole at rest in an asymptotically homogeneous medium - for stiff polytropic equations of state (EOSs) with adiabatic indices $Γ> 5/3$. A general relativistic treatment is required to determine their accretion rates, for which we provide exact expressions. We discuss several qualitative differences b…
▽ More
We revisit Bondi accretion - steady-state, adiabatic, spherical gas flow onto a Schwarzschild black hole at rest in an asymptotically homogeneous medium - for stiff polytropic equations of state (EOSs) with adiabatic indices $Γ> 5/3$. A general relativistic treatment is required to determine their accretion rates, for which we provide exact expressions. We discuss several qualitative differences between results for soft and stiff EOSs - including the appearance of a minimum steady-state accretion rate for EOSs with $Γ\geq 5/3$ - and explore limiting cases in order to examine these differences. As an example we highlight results for $Γ= 2$, which is often used in numerical simulations to model the EOS of neutron stars. We also discuss a special case with this index, the ultra-relativistic `causal' EOS, $P = ρ$. The latter serves as a useful limit for the still undetermined neutron-star EOS above nuclear density. The results are useful, for example, to estimate the accretion rate onto a mini-black hole residing at the center of a neutron star.
△ Less
Submitted 7 July, 2021; v1 submitted 21 January, 2021;
originally announced January 2021.
-
Black hole-neutron star coalescence: effects of the neutron star spin on jet launching and dynamical ejecta mass
Authors:
Milton Ruiz,
Vasileios Paschalidis,
Antonios Tsokaros,
Stuart L. Shapiro
Abstract:
Black hole-neutron star (BHNS) mergers are thought to be sources of gravitational waves (GWs) with coincident electromagnetic (EM) counterparts. To further probe whether these systems are viable progenitors of short gamma-ray bursts (sGRBs) and kilonovae, and how one may use (the lack of) EM counterparts associated with LIGO/Virgo candidate BHNS GW events to sharpen parameter estimation, we study…
▽ More
Black hole-neutron star (BHNS) mergers are thought to be sources of gravitational waves (GWs) with coincident electromagnetic (EM) counterparts. To further probe whether these systems are viable progenitors of short gamma-ray bursts (sGRBs) and kilonovae, and how one may use (the lack of) EM counterparts associated with LIGO/Virgo candidate BHNS GW events to sharpen parameter estimation, we study the impact of neutron star spin in BHNS mergers. Using dynamical spacetime magnetohydrodynamic simulations of BHNSs initially on a quasicircular orbit, we survey configurations that differ in the BH spin ($a_{\rm BH}/M_{\rm BH}=0$ and $0.75$), the NS spin ($a_{\rm NS}/M_{\rm NS}=-0.17,\,0,\,0.23$ and $0.33$), and the binary mass ratio ($q\equiv M_{\rm BH}:M_{\rm NS}=3:1$ and $5:1$). The general trend we find is that increasing the NS prograde spin increases both the rest mass of the accretion disk onto the remnant black hole, and the rest mass of dynamically ejected matter. By a time~$Δt\sim 3500-5500M\sim 88-138(M_{\rm NS}/1.4M_\odot)\,\rm ms$ after the peak gravitational wave amplitude, a magnetically--driven jet is launched only for $q=3:1$ regardless of the initial NS spin. The lifetime of the jets [$Δt\sim 0.5-0.8(M_{\rm NS}/1.4 M_\odot)\,\rm s$] and their outgoing Poynting luminosity [$L_{\rm Poyn}\sim 10^{51.5\pm 0.5}\,\rm erg/s$] are consistent with typical sGRBs luminosities and expectations from the Blandford-Znajek mechanism. By the time we terminate our simulations, we do not observe either an outflow or a large-scale magnetic field collimation for the other systems we considered. The mass range of dynamically ejected matter is $10^{-4.5}-10^{-2}~(M_{\rm NS}/1.4M_\odot)M_\odot$, which can power kilonovae with peak bolometric luminosities $L_{\rm knova}\sim 10^{40}-10^{41.4}$ erg/s with rise times $\lesssim 6.5\,\rm h$ and potentially detectable by the LSST.
△ Less
Submitted 1 January, 2021; v1 submitted 17 November, 2020;
originally announced November 2020.
-
Gravitational Waves from Disks Around Spinning Black Holes: Simulations in Full General Relativity
Authors:
Erik Wessel,
Vasileios Paschalidis,
Antonios Tsokaros,
Milton Ruiz,
Stuart L. Shapiro
Abstract:
We present fully general-relativistic numerical evolutions of self-gravitating tori around spinning black holes with dimensionless spin $a/M = 0.7$ parallel or anti-parallel to the disk angular momentum. The initial disks are unstable to the hydrodynamic Papaloizou-Pringle Instability which causes them to grow persistent orbiting matter clumps. The effect of black hole spin on the growth and satur…
▽ More
We present fully general-relativistic numerical evolutions of self-gravitating tori around spinning black holes with dimensionless spin $a/M = 0.7$ parallel or anti-parallel to the disk angular momentum. The initial disks are unstable to the hydrodynamic Papaloizou-Pringle Instability which causes them to grow persistent orbiting matter clumps. The effect of black hole spin on the growth and saturation of the instability is assessed. We find that the instability behaves similarly to prior simulations with non-spinning black holes, with a shift in frequency due to spin-induced changes in disk orbital period. Copious gravitational waves are generated by these systems, and we analyze their detectability by current and future gravitational wave observatories for large range of masses. We find that systems of $10 M_\odot$ - relevant for black hole-neutron star mergers - are detectable by Cosmic Explorer out to $\sim300$ Mpc, while DECIGO (LISA) will be able to detect systems of $1000 M_\odot$ ($10^5M_\odot$) - relevant for disks forming in collapsing supermassive stars - out to cosmological redshift of $z\sim5$ ($z\sim 1$). Computing the accretion rate of these systems we find that these systems may also be promising sources of coincident electromagnetic signals.
△ Less
Submitted 8 November, 2020;
originally announced November 2020.
-
Relativistic radiation hydrodynamics in a reference-metric formulation
Authors:
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
We adopt a two-moment formalism, together with a reference-metric approach, to express the equations of relativistic radiation hydrodynamics in a form that is well-suited for numerical implementations in curvilinear coordinates. We illustrate the approach by employing a gray opacity in an optically thick medium. As numerical demonstrations we present results for two test problems, namely stationar…
▽ More
We adopt a two-moment formalism, together with a reference-metric approach, to express the equations of relativistic radiation hydrodynamics in a form that is well-suited for numerical implementations in curvilinear coordinates. We illustrate the approach by employing a gray opacity in an optically thick medium. As numerical demonstrations we present results for two test problems, namely stationary, slab-symmetric solutions in flat spacetimes, including shocks, and heated Oppenheimer-Snyder collapse to a black hole. For the latter, we carefully analyze the transition from an initial transient to a post-transient phase that is well described by an analytically-known diffusion solution. We discuss the properties of the numerical solution when rendered in moving-puncture coordinates.
△ Less
Submitted 18 September, 2020;
originally announced September 2020.
-
Magnetic Ergostars, Jet Formation and Gamma-Ray Bursts: Ergoregions versus Horizons
Authors:
Milton Ruiz,
Antonios Tsokaros,
Stuart L. Shapiro,
Kyle C. Nelli,
Sam Qunell
Abstract:
We perform the first fully general relativistic, magnetohydrodynamic simulations of dynamically stable hypermassive neutron stars with and without ergoregions to assess the impact of ergoregions on launching magnetically--driven outflows. The hypermassive neutron stars are modeled by a compressible and causal equation of state and are initially endowed with a dipolar magnetic field extending from…
▽ More
We perform the first fully general relativistic, magnetohydrodynamic simulations of dynamically stable hypermassive neutron stars with and without ergoregions to assess the impact of ergoregions on launching magnetically--driven outflows. The hypermassive neutron stars are modeled by a compressible and causal equation of state and are initially endowed with a dipolar magnetic field extending from the stellar interior into its exterior. We find that, after a few Alfvén times, magnetic field lines in the ergostar (star that contains ergoregions) and the normal star have been tightly wound in both cases into a helical funnel within which matter begins to flow outward. The maximum Lorentz factor in the outflow is $Γ_L \sim 2.5$, while the force-free parameter holds at $B^2/8πρ_0\lesssim 10$. These values are incompatible with highly relativistic, magnetically-driven outflows (jets) and short $γ$-ray bursts. We compare these results with those of a spinning black hole surrounded by a magnetized, massless accretion disk that launches a bona fide jet. Our simulations suggest that the Blandford-Znajek mechanism for launching relativistic jets only operates when a black hole is present, though the Poynting luminosity in all cases is comparable. Therefore, one cannot distinguish a magnetized, accreting black hole from a magnetized hypermassive neutron star in the so-called mass-gap based solely on the value of the observed Poynting luminosity. These results complement our previous studies of supramassive remnants and suggest that it would be challenging for either normal neutron stars or ergostars in a hypermassive state to be the progenitors of short $γ$-ray bursts.
△ Less
Submitted 9 November, 2020; v1 submitted 18 September, 2020;
originally announced September 2020.
-
GW190814: Spin and equation of state of a neutron star companion
Authors:
Antonios Tsokaros,
Milton Ruiz,
Stuart L. Shapiro
Abstract:
The recent discovery by LIGO/Virgo of a merging binary having a $\sim 23 M_\odot$ black hole and a $\sim 2.6 M_\odot$ compact companion has triggered a debate regarding the nature of the secondary, which falls into the so-called mass gap. Here we explore some consequences of the assumption that the secondary was a neutron star (NS). We show with concrete examples of heretofore viable equations of…
▽ More
The recent discovery by LIGO/Virgo of a merging binary having a $\sim 23 M_\odot$ black hole and a $\sim 2.6 M_\odot$ compact companion has triggered a debate regarding the nature of the secondary, which falls into the so-called mass gap. Here we explore some consequences of the assumption that the secondary was a neutron star (NS). We show with concrete examples of heretofore viable equations of state (EOSs) that rapid uniform rotation may neither be necessary for some EOSs nor sufficient for others to explain the presence of a NS. Absolute upper limits for the maximum mass of a spherical NS derived from GW170817 already suggest that this unknown compact companion might be a slowly or even a nonrotating NS. However several soft NS EOSs favored by GW170817 with maximum spherical masses $\lesssim 2.1 M_\odot$ cannot be invoked to explain this object, even allowing for maximum uniform rotation. By contrast, sufficiently stiff EOSs that yield $2.6 M_\odot$ NSs which are slowly rotating or, in some cases, nonrotating, and are compatible with GW170817 and the results of NICER, can account for the black hole companion.
△ Less
Submitted 14 December, 2020; v1 submitted 10 July, 2020;
originally announced July 2020.
-
Locating ergostar models in parameter space
Authors:
Antonios Tsokaros,
Milton Ruiz,
Stuart L. Shapiro
Abstract:
Recently, we have shown that dynamically stable ergostar solutions (equilibrium neutron stars that contain an ergoregion) with a compressible and causal equation of state exist [A. Tsokaros, M. Ruiz, L. Sun, S. L. Shapiro, and K. Uryū, Phys. Rev. Lett. 123, 231103 (2019)]. These stars are hypermassive, differentially rotating, and highly compact. In this work, we make a systematic study of equilib…
▽ More
Recently, we have shown that dynamically stable ergostar solutions (equilibrium neutron stars that contain an ergoregion) with a compressible and causal equation of state exist [A. Tsokaros, M. Ruiz, L. Sun, S. L. Shapiro, and K. Uryū, Phys. Rev. Lett. 123, 231103 (2019)]. These stars are hypermassive, differentially rotating, and highly compact. In this work, we make a systematic study of equilibrium models in order to locate the position of ergostars in parameter space. We adopt four equations of state that differ in the matching density of a maximally stiff core. By constructing a large number of models both with uniform and differential rotation of different degrees, we identify the parameters for which ergostars appear. We find that the most favorable conditions for the appearance of dynamically stable ergostars are a significant finite density close to the surface of the star (i.e., similar to self-bound quark stars) and a small degree of differential rotation.
△ Less
Submitted 30 March, 2020; v1 submitted 4 February, 2020;
originally announced February 2020.
-
Prospects for Fundamental Physics with LISA
Authors:
Enrico Barausse,
Emanuele Berti,
Thomas Hertog,
Scott A. Hughes,
Philippe Jetzer,
Paolo Pani,
Thomas P. Sotiriou,
Nicola Tamanini,
Helvi Witek,
Kent Yagi,
Nicolas Yunes,
T. Abdelsalhin,
A. Achucarro,
K. V. Aelst,
N. Afshordi,
S. Akcay,
L. Annulli,
K. G. Arun,
I. Ayuso,
V. Baibhav,
T. Baker,
H. Bantilan,
T. Barreiro,
C. Barrera-Hinojosa,
N. Bartolo
, et al. (296 additional authors not shown)
Abstract:
In this paper, which is of programmatic rather than quantitative nature, we aim to further delineate and sharpen the future potential of the LISA mission in the area of fundamental physics. Given the very broad range of topics that might be relevant to LISA, we present here a sample of what we view as particularly promising directions, based in part on the current research interests of the LISA sc…
▽ More
In this paper, which is of programmatic rather than quantitative nature, we aim to further delineate and sharpen the future potential of the LISA mission in the area of fundamental physics. Given the very broad range of topics that might be relevant to LISA, we present here a sample of what we view as particularly promising directions, based in part on the current research interests of the LISA scientific community in the area of fundamental physics. We organize these directions through a "science-first" approach that allows us to classify how LISA data can inform theoretical physics in a variety of areas. For each of these theoretical physics classes, we identify the sources that are currently expected to provide the principal contribution to our knowledge, and the areas that need further development. The classification presented here should not be thought of as cast in stone, but rather as a fluid framework that is amenable to change with the flow of new insights in theoretical physics.
△ Less
Submitted 27 April, 2020; v1 submitted 27 January, 2020;
originally announced January 2020.
-
Magnetohydrodynamic Simulations of Binary Neutron Star Mergers in General Relativity: Effects of Magnetic Field Orientation on Jet Launching
Authors:
Milton Ruiz,
Antonios Tsokaros,
Stuart L. Shapiro
Abstract:
Binary neutron star (NSNS) mergers can be sources of gravitational waves coincident with electromagnetic counterpart emission. To solidify their role as multimessenger sources, we present fully 3D, general relativistic, magnetohydrodynamic simulations of spinning NSNSs initially on quasicircular orbits that merge and undergo delayed collapse to a black hole (BH). The NSNSs consist of two identical…
▽ More
Binary neutron star (NSNS) mergers can be sources of gravitational waves coincident with electromagnetic counterpart emission. To solidify their role as multimessenger sources, we present fully 3D, general relativistic, magnetohydrodynamic simulations of spinning NSNSs initially on quasicircular orbits that merge and undergo delayed collapse to a black hole (BH). The NSNSs consist of two identical stars modeled as $Γ=2$ polytropes with spin $χ_{NS}= 0.36$ aligned along the direction of the total orbital angular momentum $L$. Each star is initially threaded by a dynamical unimportant interior dipole B-field. The field is extended into the exterior where a nearly force-free magnetosphere resembles that of a pulsar. The magnetic dipole moment $μ$ is either aligned or perpendicular to $L$ and has the same initial magnitude for each orientation. For comparison, we also impose symmetry across the orbital plane in one case where $μ$ in both stars is aligned along $L$. We find that the lifetime of the transient hypermassive neutron star remnant, the jet launching time, and the ejecta are very sensitive to the B-field orientation. By contrast, the physical properties of the BH + disk remnant, such as the mass and spin of the BH, the accretion rate, and the electromagnetic luminosity, are roughly independent of the initial B-field orientation. In addition, we find imposing symmetry across the orbital plane does not play a significant role in the final outcome of the mergers. Our results show that an incipient jet emerges only when the seed B-field has a sufficiently large-scale poloidal component aligned to $L$. The lifetime [$Δt\gtrsim 140(M_{NS}/1.625M_\odot)\rm ms$] and Poynting luminosities [$L_{EM}\simeq 10^{52}$erg/s] of the jet, when it forms, are consistent with typical short gamma-ray bursts, as well as with the Blandford--Znajek mechanism for launching jets.
△ Less
Submitted 19 March, 2020; v1 submitted 24 January, 2020;
originally announced January 2020.
-
Faint Dark Matter Annihilation Signals and the Milky Way's Supermassive Black Hole
Authors:
Barry T. Chiang,
Stuart L. Shapiro,
Jessie Shelton
Abstract:
A wide range of mechanisms predict present-day s-wave dark matter (DM) annihilation cross-sections that are orders of magnitude below current experimental sensitivity. We explore the capability of DM density spikes around the Milky Way's supermassive black hole to probe such faint signals of DM annihilations, considering a range of possible spike and halo distributions. As an exemplar of a theory…
▽ More
A wide range of mechanisms predict present-day s-wave dark matter (DM) annihilation cross-sections that are orders of magnitude below current experimental sensitivity. We explore the capability of DM density spikes around the Milky Way's supermassive black hole to probe such faint signals of DM annihilations, considering a range of possible spike and halo distributions. As an exemplar of a theory with a suppressed s-wave annihilation cross-section, we consider a hidden sector axion portal model of DM. In this model, the leading contribution to the annihilation cross-section in the early universe is p-wave, while s-wave annihilations occur at higher order in the coupling constant. We provide a unified treatment of DM freezeout in this model including both s- and p-wave annihilations and analytically determine the photon spectrum for the dominant DM annihilation process in the universe today. We find that Fermi and H.E.S.S. observations of the Galactic Center offer excellent sensitivity to this model over a wide range of parameter space, with prospects depending sensitively on the properties of the DM spike as well as the central halo.
△ Less
Submitted 9 July, 2020; v1 submitted 19 December, 2019;
originally announced December 2019.
-
Enabling real-time multi-messenger astrophysics discoveries with deep learning
Authors:
E. A. Huerta,
Gabrielle Allen,
Igor Andreoni,
Javier M. Antelis,
Etienne Bachelet,
Bruce Berriman,
Federica Bianco,
Rahul Biswas,
Matias Carrasco,
Kyle Chard,
Minsik Cho,
Philip S. Cowperthwaite,
Zachariah B. Etienne,
Maya Fishbach,
Francisco Förster,
Daniel George,
Tom Gibbs,
Matthew Graham,
William Gropp,
Robert Gruendl,
Anushri Gupta,
Roland Haas,
Sarah Habib,
Elise Jennings,
Margaret W. G. Johnson
, et al. (35 additional authors not shown)
Abstract:
Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravit…
▽ More
Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.
△ Less
Submitted 26 November, 2019;
originally announced November 2019.
-
Great Impostors: Extremely Compact, Merging Binary Neutron Stars in the Mass Gap Posing as Binary Black Holes
Authors:
Antonios Tsokaros,
Milton Ruiz,
Stuart L. Shapiro,
Lunan Sun,
Kōji Uryū
Abstract:
Can one distinguish a binary black hole undergoing a merger from a binary neutron star if the individual compact companions have masses that fall inside the so-called mass gap of $3-5\ M_\odot$? For neutron stars, achieving such masses typically requires extreme compactness and in this work we present initial data and evolutions of binary neutron stars initially in quasiequilibrium circular orbits…
▽ More
Can one distinguish a binary black hole undergoing a merger from a binary neutron star if the individual compact companions have masses that fall inside the so-called mass gap of $3-5\ M_\odot$? For neutron stars, achieving such masses typically requires extreme compactness and in this work we present initial data and evolutions of binary neutron stars initially in quasiequilibrium circular orbits having a compactness $C=0.336$. These are the most compact, nonvacuum, quasiequilibrium binary objects that have been constructed and evolved to date, including boson stars. The compactness achieved is only slightly smaller than the maximum possible imposed by causality, $C_{\rm max}=0.355$, which requires the sound speed to be less than the speed of light. By comparing the emitted gravitational waveforms from the late inspiral to merger and postmerger phases between such a binary neutron star vs a binary black hole of the same total mass we identify concrete measurements that serve to distinguish them. With that level of compactness, the binary neutron stars exhibit no tidal disruption up until merger, whereupon a prompt collapse is initiated even before a common core forms. Within the accuracy of our simulations the black hole remnants from both binaries exhibit ringdown radiation that is not distinguishable from a perturbed Kerr spacetime. However, their inspiral leads to phase differences of the order of $\sim 5$ rad over an $\sim 81$ km separation (1.7 orbits) while typical neutron stars exhibit phase differences of $\geq 20$ rad. Although a difference of $\sim 5$ rad can be measured by current gravitational wave laser interferometers (e.g., aLIGO/Virgo), uncertainties in the individual masses and spins will likely prevent distinguishing such compact, massive neutron stars from black holes.
△ Less
Submitted 26 February, 2020; v1 submitted 15 November, 2019;
originally announced November 2019.
-
Dynamically stable ergostars exist: General relativistic models and simulations
Authors:
Antonios Tsokaros,
Milton Ruiz,
Lunan Sun,
Stuart L. Shapiro,
Kōji Uryū
Abstract:
We construct the first dynamically stable ergostars (equilibrium neutron stars that contain an ergoregion) for a compressible, causal equation of state. We demonstrate their stability by evolving both strict and perturbed equilibrium configurations in full general relativity for over a hundred dynamical timescales ($\gtrsim 30$ rotational periods) and observing their stationary behavior. This stab…
▽ More
We construct the first dynamically stable ergostars (equilibrium neutron stars that contain an ergoregion) for a compressible, causal equation of state. We demonstrate their stability by evolving both strict and perturbed equilibrium configurations in full general relativity for over a hundred dynamical timescales ($\gtrsim 30$ rotational periods) and observing their stationary behavior. This stability is in contrast to earlier models which prove radially unstable to collapse. Our solutions are highly differentially rotating hypermassive neutron stars with a corresponding spherical compaction of $C=0.3$. Such ergostars can provide new insights into the geometry of spacetimes around highly compact, rotating objects and on the equation of state at supranuclear densities. Ergostars may form as remnants of extreme binary neutron star mergers and possibly provide another mechanism for powering short gamma-ray bursts.
△ Less
Submitted 3 November, 2019; v1 submitted 8 July, 2019;
originally announced July 2019.
-
Dark matter heating of gas accreting onto Sgr A$^*$
Authors:
Elizabeth R. Bennewitz,
Cristian Gaidau,
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
We study effects of heating by dark matter (DM) annihilation on black hole gas accretion. We observe that, for reasonable assumptions about DM densities in spikes around supermassive black holes, as well as DM masses and annihilation cross-sections within the standard WIMP model, heating by DM annihilation may have an appreciable effect on the accretion onto Sgr A$^*$ in the Galactic center. Motiv…
▽ More
We study effects of heating by dark matter (DM) annihilation on black hole gas accretion. We observe that, for reasonable assumptions about DM densities in spikes around supermassive black holes, as well as DM masses and annihilation cross-sections within the standard WIMP model, heating by DM annihilation may have an appreciable effect on the accretion onto Sgr A$^*$ in the Galactic center. Motivated by this observation we study the effects of such heating on Bondi accretion, i.e. spherically symmetric, steady-state Newtonian accretion onto a black hole. We consider different adiabatic indices for the gas, and different power-law exponents for the DM density profile. We find that typical transonic solutions with heating have a significantly reduced accretion rate. However, for many plausible parameters, transonic solutions do not exist, suggesting a breakdown of the underlying assumptions of steady-state Bondi accretion. Our findings indicate that heating by DM annihilation may play an important role in the accretion onto supermassive black holes at the center of galaxies, and may help explain the low accretion rate observed for Sgr A$^*$.
△ Less
Submitted 3 October, 2019; v1 submitted 28 June, 2019;
originally announced July 2019.
-
On the dynamical stability of quasi-toroidal differentially rotating neutron stars
Authors:
Pedro L. Espino,
Vasileios Paschalidis,
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
We investigate the dynamical stability of relativistic, differentially rotating, quasi-toroidal models of neutron stars through hydrodynamical simulations in full general relativity. We find that all quasi-toroidal configurations studied in this work are dynamically unstable against the growth of non-axisymmetric modes. Both one-arm and bar mode instabilities grow during their evolution. We find t…
▽ More
We investigate the dynamical stability of relativistic, differentially rotating, quasi-toroidal models of neutron stars through hydrodynamical simulations in full general relativity. We find that all quasi-toroidal configurations studied in this work are dynamically unstable against the growth of non-axisymmetric modes. Both one-arm and bar mode instabilities grow during their evolution. We find that very high rest mass configurations collapse to form black holes. Our calculations suggest that configurations whose rest mass is less than the binary neutron star threshold mass for prompt collapse to black hole transition dynamically to spheroidal, differentially rotating stars that are dynamically stable, but secularly unstable. Our study shows that the existence of extreme quasi-toroidal neutron star equilibrium solutions does not imply that long-lived binary neutron star merger remnants can be much more massive than previously found. Finally, we find models that are initially supra-Kerr ($J/M^2>1$) and undergo catastrophic collapse on a dynamical timescale, in contrast to what was found in earlier works. However, cosmic censorship is respected in all of our cases. Our work explicitly demonstrates that exceeding the Kerr bound in rotating neutron star models does not imply dynamical stability.
△ Less
Submitted 20 June, 2019;
originally announced June 2019.
-
Maximally Rotating Supermassive Stars at the Onset of Collapse: Effects of Gas Pressure
Authors:
Kenneth A. Dennison,
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
The "direct collapse" scenario has emerged as a promising evolutionary track for the formation of supermassive black holes early in the Universe. In an idealized version of such a scenario, a uniformly rotating supermassive star spinning at the mass-shedding (Keplerian) limit collapses gravitationally after it reaches a critical configuration. Under the assumption that the gas is dominated by radi…
▽ More
The "direct collapse" scenario has emerged as a promising evolutionary track for the formation of supermassive black holes early in the Universe. In an idealized version of such a scenario, a uniformly rotating supermassive star spinning at the mass-shedding (Keplerian) limit collapses gravitationally after it reaches a critical configuration. Under the assumption that the gas is dominated by radiation pressure, this critical configuration is characterized by unique values of the dimensionless parameters $J/M^2$ and $R_p/M$, where $J$ is the angular momentum, $R_p$ the polar radius, and $M$ the mass. Motivated by a previous perturbative treatment we adopt a fully nonlinear approach to evaluate the effects of gas pressure on these dimensionless parameters for a large range of masses. We find that gas pressure has a significant effect on the critical configuration even for stellar masses as large as $M \simeq 10^6 M_{\odot}$. We also calibrate two approximate treatments of the gas pressure perturbation in a comparison with the exact treatment, and find that one commonly used approximation in particular results in increasing deviations from the exact treatment as the mass decreases, and the effects of gas pressure increase. The other approximation, however, proves to be quite robust for all masses $M \gtrsim 10^4 M_{\odot}$.
△ Less
Submitted 12 July, 2019; v1 submitted 10 June, 2019;
originally announced June 2019.
-
Effect of spin on the inspiral of binary neutron stars
Authors:
Antonios Tsokaros,
Milton Ruiz,
Vasileios Paschalidis,
Stuart L. Shapiro,
Kōji Uryū
Abstract:
We perform long-term simulations of spinning binary neutron stars, with our highest dimensionless spin being $χ\sim 0.32$. To assess the importance of spin during the inspiral we vary the spin, and also use two equations of state, one that consists of plain nuclear matter and produces compact stars (SLy), and a hybrid one that contains both nuclear and quark matter and leads to larger stars (ALF2)…
▽ More
We perform long-term simulations of spinning binary neutron stars, with our highest dimensionless spin being $χ\sim 0.32$. To assess the importance of spin during the inspiral we vary the spin, and also use two equations of state, one that consists of plain nuclear matter and produces compact stars (SLy), and a hybrid one that contains both nuclear and quark matter and leads to larger stars (ALF2). Using high resolution that has grid spacing $Δx\sim 98$ m on the finest refinement level, we find that the effects of spin in the phase evolution of a binary system can be larger than the one that comes from tidal forces. Our calculations demonstrate explicitly that although tidal effects are dominant for small spins ($\lesssim 0.1$), this is no longer true when the spins are larger, but still much smaller than the Keplerian limit.
△ Less
Submitted 31 May, 2019;
originally announced June 2019.
-
Gravity and Light: Combining Gravitational Wave and Electromagnetic Observations in the 2020s
Authors:
R. J. Foley,
K. D. Alexander,
I. Andreoni,
I. Arcavi,
K. Auchettl,
J. Barnes,
G. Baym,
E. C. Bellm,
A. M. Beloborodov,
N. Blagorodnova,
J. P. Blakeslee,
P. R. Brady,
M. Branchesi,
J. S. Brown,
N. Butler,
M. Cantiello,
R. Chornock,
D. O. Cook,
J. Cooke,
D. L. Coppejans,
A. Corsi,
S. M. Couch,
M. W. Coughlin,
D. A. Coulter,
P. S. Cowperthwaite
, et al. (88 additional authors not shown)
Abstract:
As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse s…
▽ More
As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
Effects of spin on magnetized binary neutron star mergers and jet launching
Authors:
Milton Ruiz,
Antonios Tsokaros,
Vasileios Paschalidis,
Stuart L. Shapiro
Abstract:
Events GW170817 and GRB 170817A provide the best confirmation so far that compact binary mergers where at least one of the companions is a neutron star (NS) can be the progenitors of short gamma-ray bursts (sGRBs). An open question for GW170817 remains the values and impact of the initial NS spins. The initial spins could possibly affect the remnant black hole (BH) mass and spin, the remnant disk…
▽ More
Events GW170817 and GRB 170817A provide the best confirmation so far that compact binary mergers where at least one of the companions is a neutron star (NS) can be the progenitors of short gamma-ray bursts (sGRBs). An open question for GW170817 remains the values and impact of the initial NS spins. The initial spins could possibly affect the remnant black hole (BH) mass and spin, the remnant disk and the formation and lifetime of a jet and its luminosity. Here we summarize our general relativistic magnetohydrodynamic simulations of spinning, NS binaries undergoing merger and delayed collapse to a BH. The binaries consist of two identical NSs, modeled as $Γ=2$ polytropes, in quasicircular orbit, each with spins $χ_{\rm{NS}}=-0.053,\,0,\,0.24$, or $0.36$. The stars are endowed initially with a dipolar magnetic field extending from the interior into the exterior, as in a radio pulsar. Following merger, the redistribution of angular momentum by magnetic braking and magnetic turbulent viscosity in the hypermassive neutron star (HMNS) remnant, along with the loss of angular momentum due to gravitational radiation, induce the formation of a massive, nearly uniformly rotating inner core surrounded by a magnetized Keplerian disk-like envelope. The HMNS eventually collapses to a BH, with spin $a/M_{\rm BH} \simeq 0.78$ independent of the initial spin of the NSs, surrounded by a magnetized accretion disk. The larger the initial NS spin the heavier the disk. After $Δt\sim 3000-4000 M \sim 45-60(M_{\rm NS}/1.625M_\odot)\rm ms$ following merger, a mildly relativistic jet is launched. The lifetime of the jet [$Δt\sim 100-140(M_{\rm NS}/1.625M_\odot)\rm ms$] and its outgoing Poynting luminosity [$L_{\rm EM}\sim 10^{51.5\pm 1}\rm erg/s$] are consistent with typical sGRBs, as well as with the Blandford--Znajek mechanism for launching jets and their associated Poynting luminosities.
△ Less
Submitted 18 April, 2019; v1 submitted 22 February, 2019;
originally announced February 2019.
-
Deep Learning for Multi-Messenger Astrophysics: A Gateway for Discovery in the Big Data Era
Authors:
Gabrielle Allen,
Igor Andreoni,
Etienne Bachelet,
G. Bruce Berriman,
Federica B. Bianco,
Rahul Biswas,
Matias Carrasco Kind,
Kyle Chard,
Minsik Cho,
Philip S. Cowperthwaite,
Zachariah B. Etienne,
Daniel George,
Tom Gibbs,
Matthew Graham,
William Gropp,
Anushri Gupta,
Roland Haas,
E. A. Huerta,
Elise Jennings,
Daniel S. Katz,
Asad Khan,
Volodymyr Kindratenko,
William T. C. Kramer,
Xin Liu,
Ashish Mahabal
, et al. (23 additional authors not shown)
Abstract:
This report provides an overview of recent work that harnesses the Big Data Revolution and Large Scale Computing to address grand computational challenges in Multi-Messenger Astrophysics, with a particular emphasis on real-time discovery campaigns. Acknowledging the transdisciplinary nature of Multi-Messenger Astrophysics, this document has been prepared by members of the physics, astronomy, compu…
▽ More
This report provides an overview of recent work that harnesses the Big Data Revolution and Large Scale Computing to address grand computational challenges in Multi-Messenger Astrophysics, with a particular emphasis on real-time discovery campaigns. Acknowledging the transdisciplinary nature of Multi-Messenger Astrophysics, this document has been prepared by members of the physics, astronomy, computer science, data science, software and cyberinfrastructure communities who attended the NSF-, DOE- and NVIDIA-funded "Deep Learning for Multi-Messenger Astrophysics: Real-time Discovery at Scale" workshop, hosted at the National Center for Supercomputing Applications, October 17-19, 2018. Highlights of this report include unanimous agreement that it is critical to accelerate the development and deployment of novel, signal-processing algorithms that use the synergy between artificial intelligence (AI) and high performance computing to maximize the potential for scientific discovery with Multi-Messenger Astrophysics. We discuss key aspects to realize this endeavor, namely (i) the design and exploitation of scalable and computationally efficient AI algorithms for Multi-Messenger Astrophysics; (ii) cyberinfrastructure requirements to numerically simulate astrophysical sources, and to process and interpret Multi-Messenger Astrophysics data; (iii) management of gravitational wave detections and triggers to enable electromagnetic and astro-particle follow-ups; (iv) a vision to harness future developments of machine and deep learning and cyberinfrastructure resources to cope with the scale of discovery in the Big Data Era; (v) and the need to build a community that brings domain experts together with data scientists on equal footing to maximize and accelerate discovery in the nascent field of Multi-Messenger Astrophysics.
△ Less
Submitted 1 February, 2019;
originally announced February 2019.
-
Magnetic Braking and Damping of Differential Rotation in Massive Stars
Authors:
Lunan Sun,
Milton Ruiz,
Stuart L. Shapiro
Abstract:
Fragmentation of highly differentially rotating massive stars that undergo collapse has been suggested as a possible channel for binary black hole formation. Such a scenario could explain the formation of the new population of massive black holes detected by the LIGO/VIRGO gravitational wave laser interferometers. We probe that scenario by performing general relativistic magnetohydrodynamic simula…
▽ More
Fragmentation of highly differentially rotating massive stars that undergo collapse has been suggested as a possible channel for binary black hole formation. Such a scenario could explain the formation of the new population of massive black holes detected by the LIGO/VIRGO gravitational wave laser interferometers. We probe that scenario by performing general relativistic magnetohydrodynamic simulations of differentially rotating massive stars supported by thermal radiation pressure plus a gas pressure perturbation. The stars are initially threaded by a dynamically weak, poloidal magnetic field confined to the stellar interior. We find that magnetic braking and turbulent viscous damping via magnetic winding and the magnetorotational instability in the bulk of the star redistribute angular momentum, damp differential rotation and induce the formation of a massive and nearly uniformly rotating inner core surrounded by a Keplerian envelope. The core + disk configuration evolves on a secular timescale and remains in quasi-stationary equilibrium until the termination of our simulations. Our results suggest that the high degree of differential rotation required for $m=2$ seed density perturbations to trigger gas fragmentation and binary black hole formation is likely to be suppressed during the normal lifetime of the star prior to evolving to the point of dynamical instability to collapse. Other cataclysmic events, such as stellar mergers leading to collapse, may therefore be necessary to reestablish sufficient differential rotation and density perturbations to drive nonaxisymmetric modes leading to binary black hole formation.
△ Less
Submitted 7 April, 2019; v1 submitted 7 December, 2018;
originally announced December 2018.
-
Jet launching from binary black hole-neutron star mergers: Dependence on black hole spin, binary mass ratio and magnetic field orientation
Authors:
Milton Ruiz,
Stuart L. Shapiro,
Antonios Tsokaros
Abstract:
Black hole-neutron star (BHNS) mergers are one of the most promising targets for multimessenger astronomy. Using general relativistic magnetohydrodynamic simulations of BHNS undergoing merger we showed that a magnetically--driven jet can be launched by the remnant if the NS is endowed with a dipole B field extending from the interior into the exterior as in a radio pulsar. These self-consistent st…
▽ More
Black hole-neutron star (BHNS) mergers are one of the most promising targets for multimessenger astronomy. Using general relativistic magnetohydrodynamic simulations of BHNS undergoing merger we showed that a magnetically--driven jet can be launched by the remnant if the NS is endowed with a dipole B field extending from the interior into the exterior as in a radio pulsar. These self-consistent studies considered a BHNS system with mass ratio $q=3:1$, BH spin $a/M_{BH}=0.75$ aligned with the total orbital angular momentum (OAM), and a NS that is irrotational, threaded by an aligned B field, and modeled by an $Γ$--law equation of state with $Γ=2$. Here, as a crucial step in establishing BHNS systems as viable progenitors of central engines that power short gamma--ray bursts (sGRBs) and thereby solidify their role as multimessenger sources, we survey different BHNS configurations that differ in BH spin ($a/M_{BH} =-0.5,\,0,\,0.5,\,0.75$), in the mass ratio ($q=3:1$ and $q=5:1$), and in the orientation of the B field (aligned and tilted by $90^\circ$ with respect to the OAM). We find that by $Δt\sim 3500M-4000M \sim 88(M_{NS}/1.4M_\odot){\rm ms}-100(M_{NS}/1.4M_\odot)\rm ms$ after the peak gravitational wave signal a jet is launched in the cases where the initial BH spin is $a/M_{BH}= 0.5$ or $0.75$. The lifetime of the jets[$Δt\sim 0.5(M_{NS}/1.4M_\odot){\rm s-0.7}(M_{NS}/1.4M_\odot)\rm s$] and their Poynting luminosities [$L_{jet}\sim 10^{51\pm 1}\rm erg/s$] are consistent with sGRBs, as well as with the Blandford--Znajek mechanism. By the time we terminate our simulations, we do not observe either an outflow or a large-scale B field collimation in the other configurations we simulate. These results suggest that future multimessenger detections from BHNSs are more likely produced by binaries with highly spinning BH companions and small tilt-angle B fields.
△ Less
Submitted 22 December, 2018; v1 submitted 19 October, 2018;
originally announced October 2018.
-
Complete initial value spacetimes containing black holes in general relativity: Application to black hole-disk systems
Authors:
Antonios Tsokaros,
Kōji Uryū,
Stuart L. Shapiro
Abstract:
We present a new initial data formulation to solve the full set of Einstein equations for spacetimes that contain a black hole under general conditions. The method can be used to construct complete initial data for spacetimes (the full metric) that contain a black hole. Contrary to most current studies the formulation requires minimal assumptions. For example, rather than imposing the form of the…
▽ More
We present a new initial data formulation to solve the full set of Einstein equations for spacetimes that contain a black hole under general conditions. The method can be used to construct complete initial data for spacetimes (the full metric) that contain a black hole. Contrary to most current studies the formulation requires minimal assumptions. For example, rather than imposing the form of the spatial conformal metric we impose 3 gauge conditions adapted to the coordinates describing the system under consideration. For stationary, axisymmetric spacetimes our method yields Kerr-Schild black holes in vacuum and rotating equilibrium neutron stars. We demonstrate the power of our new method by solving for the first time the whole system of Einstein equations for a nonaxisymmetric, self-gravitating torus in the presence of a black hole. The black hole has dimensionless spin $J_{\rm bh}/M_{\rm bh}^2=0.9918$, a rotation axis tilted at a $30^\circ$ angle with respect to the angular momentum of the disk, and a mass of $\sim 1/5$ of the disk.
△ Less
Submitted 26 February, 2019; v1 submitted 5 October, 2018;
originally announced October 2018.
-
Constant circulation sequences of binary neutron stars and their spin characterization
Authors:
Antonios Tsokaros,
Koji Uryu,
Milton Ruiz,
Stuart L. Shapiro
Abstract:
For isentropic fluids, dynamical evolution of a binary system conserves the baryonic mass and circulation; therefore, sequences of constant rest mass and constant circulation are of particular importance. In this work, we present the extension of our Compact Object CALculator (\cocal{}) code to compute such quasiequilibria and compare them with the well-known corotating and irrotational sequences,…
▽ More
For isentropic fluids, dynamical evolution of a binary system conserves the baryonic mass and circulation; therefore, sequences of constant rest mass and constant circulation are of particular importance. In this work, we present the extension of our Compact Object CALculator (\cocal{}) code to compute such quasiequilibria and compare them with the well-known corotating and irrotational sequences, the latter being the simplest, zero-circulation case. The circulation as a measure of the spin for a neutron star in a binary system has the advantage of being exactly calculable since it is a local quantity. To assess the different measures of spin, such as the angular velocity of the star, the quasilocal, dimensionless spin parameter $J/M^2$, or the circulation $\mathcal{C}$, we first compute sequences of single, uniformly rotating stars and describe how the different spin diagnostics are related to each other. The connection to spinning binary systems is accomplished through the concept of circulation and the use of the constant rotational velocity formulation. Finally, we explore a modification of the latter formulation that naturally leads to differentially rotating binary systems.
△ Less
Submitted 31 December, 2018; v1 submitted 21 September, 2018;
originally announced September 2018.
-
Star Clusters, Self-Interacting Dark Matter Halos and Black Hole Cusps: The Fluid Conduction Model and its Extension to General Relativity
Authors:
Stuart L. Shapiro
Abstract:
We adopt the fluid conduction approximation to study the evolution of spherical star clusters and self-interacting dark matter (SIDM) halos. We also explore the formation and dynamical impact of density cusps that arise in both systems due to the presence of a massive, central black hole. The large N-body, self-gravitating systems we treat are "weakly-collisional": the mean free time between star…
▽ More
We adopt the fluid conduction approximation to study the evolution of spherical star clusters and self-interacting dark matter (SIDM) halos. We also explore the formation and dynamical impact of density cusps that arise in both systems due to the presence of a massive, central black hole. The large N-body, self-gravitating systems we treat are "weakly-collisional": the mean free time between star or SIDM particle collisions is much longer than their characteristic crossing (dynamical) time scale, but shorter than the system lifetime. The fluid conduction model reliably tracks the "gravothermal catastrophe" in star clusters and SIDM halos without black holes. For a star cluster with a massive, central black hole, this approximation reproduces the familiar Bahcall-Wolf quasistatic density cusp for the stars bound to the black hole and shows how the cusp halts the "gravothermal catastrophe" and causes the cluster to re-expand. An SIDM halo with an initial black hole central density spike that matches onto to an exterior NFW profile relaxes to a core-halo structure with a central density cusp determined by the velocity dependence of the SIDM interaction cross section. The success and relative simplicity of the fluid conduction approach in evolving such "weakly-collisional", quasiequilibrium Newtonian systems motivates its extension to relativistic systems. We present a general relativistic extension here.
△ Less
Submitted 7 September, 2018;
originally announced September 2018.
-
Simulating the Magnetorotational Collapse of Supermassive Stars: Incorporating Gas Pressure Perturbations and Different Rotation Profiles
Authors:
Lunan Sun,
Milton Ruiz,
Stuart L. Shapiro
Abstract:
Collapsing supermassive stars (SMSs) with masses $M \gtrsim 10^{4-6}M_\odot$ have long been speculated to be the seeds that can grow and become supermassive black holes (SMBHs). We previously performed GRMHD simulations of marginally stable magnetized $Γ= 4/3$ polytropes uniformly rotating at the mass-shedding limit to model the direct collapse of SMSs. These configurations are supported entirely…
▽ More
Collapsing supermassive stars (SMSs) with masses $M \gtrsim 10^{4-6}M_\odot$ have long been speculated to be the seeds that can grow and become supermassive black holes (SMBHs). We previously performed GRMHD simulations of marginally stable magnetized $Γ= 4/3$ polytropes uniformly rotating at the mass-shedding limit to model the direct collapse of SMSs. These configurations are supported entirely by thermal radiation pressure and model SMSs with $M \gtrsim 10^{6}M_\odot$. We found that around $90\%$ of the initial stellar mass forms a spinning black hole (BH) surrounded by a massive, hot, magnetized torus, which eventually launches an incipient jet. Here we perform GRMHD simulations of $Γ\gtrsim 4/3$, polytropes to account for the perturbative role of gas pressure in SMSs with $M \lesssim 10^{6}M_\odot$. We also consider different initial stellar rotation profiles. The stars are initially seeded with a dynamically weak dipole magnetic field that is either confined to the stellar interior or extended from its interior into the stellar exterior. We find that the mass of the BH remnant is $90\%-99\%$ of the initial stellar mass, depending sharply on $Γ-4/3$ as well as on the initial stellar rotation profile. After $t\sim 250-550M\approx 1-2\times 10^3(M/10^6M_\odot)$s following the BH formation, a jet is launched and it lasts for $\sim 10^4-10^5(M/10^6M_\odot)$s, consistent with the duration of long gamma-ray bursts. Our results suggest that the Blandford-Znajek mechanism powers the jet. They are also in agreement with our proposed universal model that estimates accretion rates and luminosities that characterize magnetized BH-disk remnant systems that launch a jet. This model helps explain why the outgoing luminosities for vastly different BH-disk formation scenarios all reside within a narrow range ($\sim 10^{52 \pm 1} \rm erg/s$), roughly independent of $M$.
△ Less
Submitted 13 November, 2018; v1 submitted 20 July, 2018;
originally announced July 2018.
-
Maximally Rotating Supermassive Stars at the Onset of Collapse: The Perturbative Effects of Gas Pressure, Magnetic Fields, Dark Matter and Dark Energy
Authors:
Satya P. Butler,
Alicia R. Lima,
Thomas W. Baumgarte,
Stuart L. Shapiro
Abstract:
The discovery of quasars at increasingly large cosmological redshifts may favor "direct collapse" as the most promising evolutionary route to the formation of supermassive black holes. In this scenario, supermassive black holes form when their progenitors - supermassive stars - become unstable to gravitational collapse. For uniformly rotating stars supported by pure radiation pressure and spinning…
▽ More
The discovery of quasars at increasingly large cosmological redshifts may favor "direct collapse" as the most promising evolutionary route to the formation of supermassive black holes. In this scenario, supermassive black holes form when their progenitors - supermassive stars - become unstable to gravitational collapse. For uniformly rotating stars supported by pure radiation pressure and spinning at the mass-shedding limit, the critical configuration at the onset of collapse is characterized by universal values of the dimensionless spin and radius parameters $J/M^2$ and $R/M$, independent of mass $M$. We consider perturbative effects of gas pressure, magnetic fields, dark matter and dark energy on these parameters, and thereby determine the domain of validity of this universality. We obtain leading-order corrections for the critical parameters and establish their scaling with the relevant physical parameters. We compare two different approaches to approximate the effects of gas pressure, which plays the most important role, find identical results for the above dimensionless parameters, and also find good agreement with recent numerical results.
△ Less
Submitted 26 July, 2018; v1 submitted 2 February, 2018;
originally announced February 2018.
-
Disks Around Merging Binary Black Holes: From GW150914 to Supermassive Black Holes
Authors:
Abid Khan,
Vasileios Paschalidis,
Milton Ruiz,
Stuart L. Shapiro
Abstract:
We perform magnetohydrodynamic simulations in full general relativity of disk accretion onto nonspinning black hole binaries with mass ratio 36:29. We survey different disk models which differ in their scale height, total size and magnetic field to quantify the robustness of previous simulations on the initial disk model. Scaling our simulations to LIGO GW150914 we find that such systems could exp…
▽ More
We perform magnetohydrodynamic simulations in full general relativity of disk accretion onto nonspinning black hole binaries with mass ratio 36:29. We survey different disk models which differ in their scale height, total size and magnetic field to quantify the robustness of previous simulations on the initial disk model. Scaling our simulations to LIGO GW150914 we find that such systems could explain possible gravitational wave and electromagnetic counterparts such as the Fermi GBM hard X-ray signal reported 0.4s after GW150915 ended. Scaling our simulations to supermassive binary black holes, we find that observable flow properties such as accretion rate periodicities, the emergence of jets throughout inspiral, merger and post-merger, disk temperatures, thermal frequencies, and the time-delay between merger and the boost in jet outflows that we reported in earlier studies display only modest dependence on the initial disk model we consider here.
△ Less
Submitted 24 February, 2018; v1 submitted 8 January, 2018;
originally announced January 2018.
-
GW170817, General Relativistic Magnetohydrodynamic Simulations, and the Neutron Star Maximum Mass
Authors:
Milton Ruiz,
Stuart L. Shapiro,
Antonios Tsokaros
Abstract:
Recent numerical simulations in general relativistic magnetohydrodynamics (GRMHD) provide useful constraints for the interpretation of the GW170817 discovery. Combining the observed data with these simulations leads to a bound on the maximum mass of a cold, spherical neutron star (the TOV limit): ${M_{\rm max}^{\rm sph}}\lesssim 2.74/β$, where $β$ is the ratio of the maximum mass of a uniformly ro…
▽ More
Recent numerical simulations in general relativistic magnetohydrodynamics (GRMHD) provide useful constraints for the interpretation of the GW170817 discovery. Combining the observed data with these simulations leads to a bound on the maximum mass of a cold, spherical neutron star (the TOV limit): ${M_{\rm max}^{\rm sph}}\lesssim 2.74/β$, where $β$ is the ratio of the maximum mass of a uniformly rotating neutron star (the supramassive limit) over the maximum mass of a nonrotating star. Causality arguments allow $β$ to be as high as $1.27$, while most realistic candidate equations of state predict $β$ to be closer to $1.2$, yielding ${M_{\rm max}^{\rm sph}}$ in the range $2.16-2.28 M_\odot$. A minimal set of assumptions based on these simulations distinguishes this analysis from previous ones, but leads to a similar estimate. There are caveats, however, and they are enumerated and discussed. The caveats can be removed by further simulations and analysis to firm up the basic argument.
△ Less
Submitted 4 May, 2018; v1 submitted 1 November, 2017;
originally announced November 2017.
-
GRMHD simulations of prompt-collapse neutron star mergers: the absence of jets
Authors:
Milton Ruiz,
Stuart L. Shapiro
Abstract:
Inspiraling and merging binary neutron stars are not only important source of gravitational waves, but also promising candidates for coincident electromagnetic counterparts. These systems are thought to be progenitors of short gamma-ray bursts (sGRBs). We have shown previously that binary neutron star mergers that undergo {\it delayed} collapse to a black hole surrounded by a {\it weighty} magneti…
▽ More
Inspiraling and merging binary neutron stars are not only important source of gravitational waves, but also promising candidates for coincident electromagnetic counterparts. These systems are thought to be progenitors of short gamma-ray bursts (sGRBs). We have shown previously that binary neutron star mergers that undergo {\it delayed} collapse to a black hole surrounded by a {\it weighty} magnetized accretion disk can drive magnetically-powered jets. We now perform magnetohydrodynamic simulations in full general relativity of binary neutron stars mergers that undergo {\it prompt} collapse to explore the possibility of jet formation from black hole-{\it light} accretion disk remnants. We find that after $t-t_{\rm BH}\sim 26(M_{\rm NS}/1.8M_\odot)$ms [$M_{\rm NS}$ is the ADM mass] following prompt black hole formation, there is no evidence of mass outflow or magnetic field collimation. The rapid formation of the black hole following merger prevents magnetic energy from approaching force-free values above the magnetic poles, which is required for the launching of a jet by the usual Blandford--Znajek mechanism. Detection of gravitational waves in coincidence with sGRBs may provide constraints on the nuclear equation of state (EOS): the fate of an NSNS merger--delayed or prompt collapse, and hence the appearance or nonappearance of an sGRB--depends on a critical value of the total mass of the binary, and this value is sensitive to the EOS.
△ Less
Submitted 1 November, 2017; v1 submitted 1 September, 2017;
originally announced September 2017.