-
Exploring Black Hole Mimickers: Electromagnetic and Gravitational Signatures of AdS Black Shells
Authors:
Suvendu Giri,
Ulf Danielsson,
Luis Lehner,
Frans Pretorius
Abstract:
We study electromagnetic and gravitational properties of AdS black shells (also referred to as AdS black bubbles) -- a class of quantum gravity motivated black hole mimickers, that in the classical limit are described as ultra compact shells of matter. We find that their electromagnetic properties are remarkably similar to black holes. We then discuss the extent to which these objects are distingu…
▽ More
We study electromagnetic and gravitational properties of AdS black shells (also referred to as AdS black bubbles) -- a class of quantum gravity motivated black hole mimickers, that in the classical limit are described as ultra compact shells of matter. We find that their electromagnetic properties are remarkably similar to black holes. We then discuss the extent to which these objects are distinguishable from black holes, both for intrinsic interest within the black shell model, and as a guide for similar efforts in other sub-classes of exotic compact objects (ECOs). We study photon rings and lensing band characteristics, relevant for very large baseline inteferometry (VLBI) observations, as well as gravitational wave observables -- quasinormal modes in the eikonal limit and the static tidal Love number for non-spinning shells -- relevant for ongoing and upcoming gravitational wave observations.
△ Less
Submitted 13 May, 2024;
originally announced May 2024.
-
Nonlinear Effects In Black Hole Ringdown From Scattering Experiments I: spin and initial data dependence of quadratic mode coupling
Authors:
Hengrui Zhu,
Justin L. Ripley,
Frans Pretorius,
Sizheng Ma,
Keefe Mitman,
Robert Owen,
Michael Boyle,
Yitian Chen,
Nils Deppe,
Lawrence E. Kidder,
Jordan Moxon,
Kyle C. Nelli,
Harald P. Pfeiffer,
Mark A. Scheel,
William Throwe,
Nils L. Vu
Abstract:
We investigate quadratic quasinormal mode coupling in black hole spacetime through numerical simulations of single perturbed black holes using both numerical relativity and second-order black hole perturbation theory. Focusing on the dominant $\ell=|m|=2$ quadrupolar modes, we find good agreement (within $\sim10\%$) between these approaches, with discrepancies attributed to truncation error and un…
▽ More
We investigate quadratic quasinormal mode coupling in black hole spacetime through numerical simulations of single perturbed black holes using both numerical relativity and second-order black hole perturbation theory. Focusing on the dominant $\ell=|m|=2$ quadrupolar modes, we find good agreement (within $\sim10\%$) between these approaches, with discrepancies attributed to truncation error and uncertainties from mode fitting. Our results align with earlier studies extracting the coupling coefficients from select binary black hole merger simulations, showing consistency for the same remnant spins. Notably, the coupling coefficient is insensitive to a diverse range of initial data, including configurations that led to a significant (up to $5\%$) increase in the remnant black hole mass. These findings present opportunities for testing the nonlinear dynamics of general relativity with ground-based gravitational wave observatories. Lastly, we provide evidence of a bifurcation in coupling coefficients between counter-rotating and co-rotating quasinormal modes as black hole spin increases.
△ Less
Submitted 1 January, 2024;
originally announced January 2024.
-
A Survey of Gravitational Waves
Authors:
Frans Pretorius
Abstract:
We review the state of the field of gravitational wave astrophysics, framing the challenges, current observations, and future prospects within the context of the predictions of Einstein's theory of general relativity.
We review the state of the field of gravitational wave astrophysics, framing the challenges, current observations, and future prospects within the context of the predictions of Einstein's theory of general relativity.
△ Less
Submitted 6 June, 2023;
originally announced June 2023.
-
Causal, stable first-order viscous relativistic hydrodynamics with ideal gas microphysics
Authors:
Alex Pandya,
Elias R. Most,
Frans Pretorius
Abstract:
We present the first numerical analysis of causal, stable first-order relativistic hydrodynamics with ideal gas microphysics, based in the formalism developed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK theory). The BDNK approach provides definitions for the conserved stress-energy tensor and baryon current, and rigorously proves causality, local well-posedness, strong hyperbolicity, and linea…
▽ More
We present the first numerical analysis of causal, stable first-order relativistic hydrodynamics with ideal gas microphysics, based in the formalism developed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK theory). The BDNK approach provides definitions for the conserved stress-energy tensor and baryon current, and rigorously proves causality, local well-posedness, strong hyperbolicity, and linear stability (about equilibrium) for the equations of motion, subject to a set of coupled nonlinear inequalities involving the undetermined model coefficients (the choice for which defines the "hydrodynamic frame"). We present a class of hydrodynamic frames derived from the relativistic ideal gas "gamma-law" equation of state which satisfy the BDNK constraints, and explore the properties of the resulting model for a series of (0+1)D and (1+1)D tests in 4D Minkowski spacetime. These tests include a comparison of the dissipation mechanisms in Eckart, BDNK, and Muller-Israel-Stewart theories, as well as investigations of the impact of hydrodynamic frame on the causality and stability properties of Bjorken flow, planar shockwave, and heat flow solutions.
△ Less
Submitted 31 December, 2022; v1 submitted 19 September, 2022;
originally announced September 2022.
-
Binary neutron star mergers in Einstein-scalar-Gauss-Bonnet gravity
Authors:
William E. East,
Frans Pretorius
Abstract:
Binary neutron star mergers, which can lead to less massive black holes relative to other known astrophysical channels, have the potential to probe modifications to general relativity that arise at smaller curvature scales compared to more massive compact object binaries. As a representative example of this, here we study binary neutron star mergers in shift-symmetric Einstein-scalar-Gauss-Bonnet…
▽ More
Binary neutron star mergers, which can lead to less massive black holes relative to other known astrophysical channels, have the potential to probe modifications to general relativity that arise at smaller curvature scales compared to more massive compact object binaries. As a representative example of this, here we study binary neutron star mergers in shift-symmetric Einstein-scalar-Gauss-Bonnet gravity using evolutions of the full, nonperturbative evolution equations. We find that the impact on the inspiral is small, even at large values of the modified gravity coupling (as expected, as neutron stars do not have scalar charge in this theory). However, postmerger there can be strong scalar effects, including radiation. When a black hole forms, it develops scalar charge, impacting the ringdown gravitational wave signal. In cases where a longer-lived remnant star persists postmerger, we find that the oscillations of the star source levels of scalar radiation similar to the black hole formation cases. In remnant stars, we further find that at coupling values comparable to the maximum value for which black hole solutions of the same mass exist, there is significant nonlinear enhancement in the scalar field, which if sufficiently large leads to a breakdown in the evolution, seemingly due to loss of hyperbolicity of the underlying equations.
△ Less
Submitted 28 November, 2022; v1 submitted 19 August, 2022;
originally announced August 2022.
-
Conservative finite volume scheme for first-order viscous relativistic hydrodynamics
Authors:
Alex Pandya,
Elias R. Most,
Frans Pretorius
Abstract:
We present the first conservative finite volume numerical scheme for the causal, stable relativistic Navier-Stokes equations developed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK). BDNK theory has arisen very recently as a promising means of incorporating entropy-generating effects (viscosity, heat conduction) into relativistic fluid models, appearing as a possible alternative to the so-called…
▽ More
We present the first conservative finite volume numerical scheme for the causal, stable relativistic Navier-Stokes equations developed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK). BDNK theory has arisen very recently as a promising means of incorporating entropy-generating effects (viscosity, heat conduction) into relativistic fluid models, appearing as a possible alternative to the so-called Müller-Israel-Stewart (MIS) theory successfully used to model quark-gluon plasma. The major difference between the two lies in the structure of the system of PDEs: BDNK theory only has a set of conservation laws, whereas MIS also includes a set of evolution equations for its dissipative degrees of freedom. The simpler structure of the BDNK PDEs in this respect allows for rigorous proofs of stability, causality, and hyperbolicity in full generality which have as yet been impossible for MIS. To capitalize on these advantages, we present the first fully conservative multi-dimensional fluid solver for the BDNK equations suitable for physical applications. The scheme includes a flux-conservative discretization, non-oscillatory reconstruction, and a central-upwind numerical flux, and is designed to smoothly transition to a high-resolution shock-capturing perfect fluid solver in the inviscid limit. We assess the robustness of our new method in a series of flat-spacetime tests for a conformal fluid, and provide a detailed comparison with previous approaches of Pandya & Pretorius (2021).
△ Less
Submitted 1 June, 2022; v1 submitted 28 January, 2022;
originally announced January 2022.
-
Dynamics and Observational Signatures of Shell-like Black Hole Mimickers
Authors:
Ulf Danielsson,
Luis Lehner,
Frans Pretorius
Abstract:
We undertake the task of studying the non-linear dynamics of quantum gravity motivated alternatives to black holes that in the classical limit appear as ultra-compact shells of matter. We develop a formalism that should be amenable to numerical solution in generic situations. For a concrete model we focus on the spherically symmetric AdS black bubble -- a shell of matter at the Buchdahl radius sep…
▽ More
We undertake the task of studying the non-linear dynamics of quantum gravity motivated alternatives to black holes that in the classical limit appear as ultra-compact shells of matter. We develop a formalism that should be amenable to numerical solution in generic situations. For a concrete model we focus on the spherically symmetric AdS black bubble -- a shell of matter at the Buchdahl radius separating a Schwarzschild exterior from an AdS interior. We construct a numerical code to study the radial dynamics of and accretion onto AdS black bubbles, with exterior matter provided by scalar fields. In doing so we develop numerical methods that could be extended to future studies beyond spherical symmetry. Regarding AdS black bubbles in particular, we find that the original prescription for the internal matter fluxes needed to stabilize the black bubble is inadequate in dynamical settings, and we propose a two parameter generalization of the flux model to fix this. To allow for more efficient surveys of parameter space, we develop a simpler numerical model adapted to spherically symmetric bubble dynamics. We identify regions of parameter space that do allow for stable black bubbles, and moreover allow control to a desired end-state after an accretion episode. Based on these results, and evolution of scalar fields on black bubble backgrounds, we speculate on some observational consequences if what are currently presumed to be black holes in the universe were actually black bubbles.
△ Less
Submitted 20 September, 2021;
originally announced September 2021.
-
Dynamical attractors in contracting spacetimes dominated by kinetically coupled scalar fields
Authors:
Anna Ijjas,
Frans Pretorius,
Paul J. Steinhardt,
David Garfinkle
Abstract:
We present non-perturbative numerical relativity simulations of slowly contracting spacetimes in which the scalar field driving slow contraction is coupled to a second scalar field through an exponential non-linear sigma model-type kinetic interaction. These models are important because they can generate a nearly scale-invariant spectrum of super-Hubble density fluctuations fully consistent with c…
▽ More
We present non-perturbative numerical relativity simulations of slowly contracting spacetimes in which the scalar field driving slow contraction is coupled to a second scalar field through an exponential non-linear sigma model-type kinetic interaction. These models are important because they can generate a nearly scale-invariant spectrum of super-Hubble density fluctuations fully consistent with cosmic microwave background observations. We show that the non-linear evolution rapidly approaches a homogeneous, isotropic and flat Friedmann- Robertson-Walker (FRW) geometry for a wide range of inhomogeneous and anisotropic initial conditions. Ultimately, we find, the kinetic coupling causes the evolution to deflect away from flat FRW and towards a novel Kasner-like stationary point, but in general this occurs on time scales that are too long to be observationally relevant.
△ Less
Submitted 20 September, 2021;
originally announced September 2021.
-
Projecting the likely importance of weak-interaction-driven bulk viscosity in neutron star mergers
Authors:
Elias R. Most,
Steven P. Harris,
Christopher Plumberg,
Mark G. Alford,
Jorge Noronha,
Jacquelyn Noronha-Hostler,
Frans Pretorius,
Helvi Witek,
Nicolás Yunes
Abstract:
In this work, we estimate how much bulk viscosity driven by Urca processes is likely to affect the gravitational wave signal of a neutron star coalescence. In the late inspiral, we show that bulk viscosity affects the binding energy at fourth post-Newtonian (PN) order. Even though this effect is enhanced by the square of the gravitational compactness, the coefficient of bulk viscosity is likely to…
▽ More
In this work, we estimate how much bulk viscosity driven by Urca processes is likely to affect the gravitational wave signal of a neutron star coalescence. In the late inspiral, we show that bulk viscosity affects the binding energy at fourth post-Newtonian (PN) order. Even though this effect is enhanced by the square of the gravitational compactness, the coefficient of bulk viscosity is likely too small to lead to observable effects in the waveform during the late inspiral, when only considering the orbital motion itself. In the post-merger, however, the characteristic time-scales and spatial scales are different, potentially leading to the opposite conclusion. We post-process data from a state-of-the-art equal-mass binary neutron star merger simulation to estimate the effects of bulk viscosity (which was not included in the simulation itself). In that scenario, we find that bulk viscosity can reach high values in regions of the merger. We compute several estimates of how much it might directly affect the global dynamics of the considered merger scenario, and find that it could become significant. Even larger effects could arise in different merger scenarios or in simulations that include non-linear effects. This assessment is reinforced by a quantitative comparison with relativistic heavy-ion collisions where such effects have been explored extensively.
△ Less
Submitted 11 July, 2021;
originally announced July 2021.
-
The Effects of Multiple Modes and Reduced Symmetry on the Rapidity and Robustness of Slow Contraction
Authors:
Anna Ijjas,
Frans Pretorius,
Paul J. Steinhardt,
Andrew P. Sullivan
Abstract:
We demonstrate that the rapidity and robustness of slow contraction in homogenizing and flattening the universe found in simulations in which the initial conditions were restricted to non-perturbative variations described by a single fourier mode along only a single spatial direction are in general enhanced if the initial variations are along two spatial directions, include multiple modes, and the…
▽ More
We demonstrate that the rapidity and robustness of slow contraction in homogenizing and flattening the universe found in simulations in which the initial conditions were restricted to non-perturbative variations described by a single fourier mode along only a single spatial direction are in general enhanced if the initial variations are along two spatial directions, include multiple modes, and thereby have reduced symmetry. Particularly significant are shear effects that only become possible when variations are allowed along two or more spatial dimensions. Based on the numerical results, we conjecture that the counterintuitive enhancement occurs because more degrees of freedom are activated which drive spacetime away from an unstable Kasner fixed point and towards the stable Friedmann-Robertson-Walker fixed point.
△ Less
Submitted 25 April, 2021;
originally announced April 2021.
-
A numerical exploration of first-order relativistic hydrodynamics
Authors:
Alex Pandya,
Frans Pretorius
Abstract:
We present the first numerical solutions of the causal, stable relativistic Navier-Stokes equations as formulated by Bemfica, Disconzi, Noronha, and Kovtun (BDNK). For this initial investigation we restrict to plane-symmetric configurations of a conformal fluid in Minkowski spacetime. We consider evolution of three classes of initial data: a smooth (initially) stationary concentration of energy, a…
▽ More
We present the first numerical solutions of the causal, stable relativistic Navier-Stokes equations as formulated by Bemfica, Disconzi, Noronha, and Kovtun (BDNK). For this initial investigation we restrict to plane-symmetric configurations of a conformal fluid in Minkowski spacetime. We consider evolution of three classes of initial data: a smooth (initially) stationary concentration of energy, a standard shock tube setup, and a smooth shockwave setup. We compare these solutions to those obtained with the Muller-Israel-Stewart (MIS) formalism, variants of which are the common tools used to model relativistic, viscous fluids. We find that for the two smooth initial data cases, simple finite difference methods are adequate to obtain stable, convergent solutions to the BDNK equations. For low viscosity, the MIS and BDNK evolutions show good agreement. At high viscosity the solutions begin to differ in regions with large gradients, and there the BDNK solutions can (as expected) exhibit violation of the weak energy condition. This behavior is transient, and the solutions evolve toward a hydrodynamic regime in a way reminiscent of an approach to a universal attractor. For the shockwave problem, we give evidence that if a hydrodynamic frame is chosen so that the maximum characteristic speed of the BDNK system is the speed of light (or larger), arbitrarily strong shockwaves are smoothly resolved. Regarding the shock tube problem, it is unclear whether discontinuous initial data is mathematically well-posed for the BDNK system, even in a weak sense. Nevertheless we attempt numerical solution, and then need to treat the perfect fluid terms using high-resolution shock-capturing methods. When such methods can successfully evolve the solution beyond the initial time, subsequent evolution agrees with corresponding MIS solutions, as well as the perfect fluid solution in the limit of zero viscosity.
△ Less
Submitted 23 June, 2021; v1 submitted 1 April, 2021;
originally announced April 2021.
-
Ultralocality and Slow Contraction
Authors:
Anna Ijjas,
Andrew P. Sullivan,
Frans Pretorius,
Paul J. Steinhardt,
William G. Cook
Abstract:
We study the detailed process by which slow contraction smooths and flattens the universe using an improved numerical relativity code that accepts initial conditions with non-perturbative deviations from homogeneity and isotropy along two independent spatial directions. Contrary to common descriptions of the early universe, we find that the geometry first rapidly converges to an inhomogeneous, spa…
▽ More
We study the detailed process by which slow contraction smooths and flattens the universe using an improved numerical relativity code that accepts initial conditions with non-perturbative deviations from homogeneity and isotropy along two independent spatial directions. Contrary to common descriptions of the early universe, we find that the geometry first rapidly converges to an inhomogeneous, spatially-curved and anisotropic ultralocal state in which all spatial gradient contributions to the equations of motion decrease as an exponential in time to negligible values. This is followed by a second stage in which the geometry converges to a homogeneous, spatially flat and isotropic spacetime. In particular, the decay appears to follow the same history whether the entire spacetime or only parts of it are smoothed by the end of slow contraction.
△ Less
Submitted 28 February, 2021;
originally announced March 2021.
-
Robustness of slow contraction to cosmic initial conditions
Authors:
Anna Ijjas,
William G. Cook,
Frans Pretorius,
Paul J. Steinhardt,
Elliot Y. Davies
Abstract:
We present numerical relativity simulations of cosmological scenarios in which the universe is smoothed and flattened by undergoing a phase of slow contraction and test their sensitivity to a wide range of initial conditions. Our numerical scheme enables the variation of all freely specifiable physical quantities that characterize the initial spatial hypersurface, such as the initial shear and spa…
▽ More
We present numerical relativity simulations of cosmological scenarios in which the universe is smoothed and flattened by undergoing a phase of slow contraction and test their sensitivity to a wide range of initial conditions. Our numerical scheme enables the variation of all freely specifiable physical quantities that characterize the initial spatial hypersurface, such as the initial shear and spatial curvature contributions as well as the initial field and velocity distributions of the scalar that drives the cosmological evolution. In particular, we include initial conditions that are far outside the perturbative regime of the well-known attractor scaling solution. We complement our numerical results by analytically performing a complete dynamical systems analysis and show that the two approaches yield consistent results.
△ Less
Submitted 8 July, 2020; v1 submitted 8 June, 2020;
originally announced June 2020.
-
Supersmoothing through Slow Contraction
Authors:
William G. Cook,
Iryna A. Glushchenko,
Anna Ijjas,
Frans Pretorius,
Paul J. Steinhardt
Abstract:
Performing a fully non-perturbative analysis using the tools of numerical general relativity, we demonstrate that a period of slow contraction is a `supersmoothing' cosmological phase that homogenizes, isotropizes and flattens the universe both classically and quantum mechanically and can do so far more robustly and rapidly than had been realized in earlier studies.
Performing a fully non-perturbative analysis using the tools of numerical general relativity, we demonstrate that a period of slow contraction is a `supersmoothing' cosmological phase that homogenizes, isotropizes and flattens the universe both classically and quantum mechanically and can do so far more robustly and rapidly than had been realized in earlier studies.
△ Less
Submitted 1 June, 2020;
originally announced June 2020.
-
Einstein-Vlasov Calculations of Structure Formation
Authors:
William E. East,
Radosław Wojtak,
Frans Pretorius
Abstract:
We study the dynamics of small inhomogeneities in an expanding universe collapsing to form bound structures using full solutions of the Einstein-Vlasov (N-body) equations. We compare these to standard Newtonian N-body solutions using quantities defined with respect to fiducial observers in order to bound relativistic effects. We focus on simplified initial conditions containing a limited range of…
▽ More
We study the dynamics of small inhomogeneities in an expanding universe collapsing to form bound structures using full solutions of the Einstein-Vlasov (N-body) equations. We compare these to standard Newtonian N-body solutions using quantities defined with respect to fiducial observers in order to bound relativistic effects. We focus on simplified initial conditions containing a limited range of length scales, but vary the inhomogeneities from small magnitude, where the Newtonian and general-relativistic calculations agree quite well, to large magnitude, where the background metric receives an order one correction. For large inhomogeneities, we find that the collapse of overdensities tends to happen faster in Newtonian calculations relative to fully general-relativistic ones. Even in this extreme regime, the differences in the spacetime evolution outside the regions of large gravitational potential and velocity are small. For standard cosmological values, we corroborate the robustness of Newtonian N-body simulations to model large scale perturbations and the related cosmic variance in the local expansion rate.
△ Less
Submitted 22 November, 2019; v1 submitted 15 August, 2019;
originally announced August 2019.
-
Binary neutron star mergers: Effects of spin and post-merger dynamics
Authors:
William E. East,
Vasileios Paschalidis,
Frans Pretorius,
Antonios Tsokaros
Abstract:
Spin can have significant effects on the electromagnetic transients accompanying binary neutron star mergers. The measurement of spin can provide important information about binary formation channels. In the absence of a strong neutron star spin prior, the degeneracy of spin with other parameters leads to significant uncertainties in their estimation, in particular limiting the power of gravitatio…
▽ More
Spin can have significant effects on the electromagnetic transients accompanying binary neutron star mergers. The measurement of spin can provide important information about binary formation channels. In the absence of a strong neutron star spin prior, the degeneracy of spin with other parameters leads to significant uncertainties in their estimation, in particular limiting the power of gravitational waves to place tight constraints on the nuclear equation of state. Thus detailed studies of highly spinning neutron star mergers are essential to understand all aspects of multimessenger observation of such events. We perform a systematic investigation of the impact of neutron star spin---considering dimensionless spin values up to $a_{\rm NS}=0.33$---on the merger of equal mass, quasicircular binary neutron stars using fully general-relativistic simulations. We find that the peak frequency of the post-merger gravitational wave signal is only weakly influenced by the neutron star spin, with cases where the spin is aligned (antialigned) with the orbital angular momentum giving slightly lower (higher) values compared to the irrotational case. We find that the one-arm instability arises in a number of cases, with some dependence on spin. Spin has a pronounced impact on the mass, velocity, and angular distribution of the dynamical ejecta, and the mass of the disk that remains outside the merger remnant. We discuss the implications of these findings on anticipated electromagnetic signals, and on constraints that have been placed on the equation of state based on multimessenger observations of GW170817.
△ Less
Submitted 12 December, 2019; v1 submitted 12 June, 2019;
originally announced June 2019.
-
Gravitational Collapse in Einstein dilaton Gauss-Bonnet Gravity
Authors:
Justin L Ripley,
Frans Pretorius
Abstract:
We present results from a numerical study of spherical gravitational collapse in shift symmetric Einstein dilaton Gauss-Bonnet (EdGB) gravity. This modified gravity theory has a single coupling parameter that when zero reduces to general relativity (GR) minimally coupled to a massless scalar field. We first show results from the weak EdGB coupling limit, where we obtain solutions that smoothly app…
▽ More
We present results from a numerical study of spherical gravitational collapse in shift symmetric Einstein dilaton Gauss-Bonnet (EdGB) gravity. This modified gravity theory has a single coupling parameter that when zero reduces to general relativity (GR) minimally coupled to a massless scalar field. We first show results from the weak EdGB coupling limit, where we obtain solutions that smoothly approach those of the Einstein-Klein-Gordon system of GR. Here, in the strong field case, though our code does not utilize horizon penetrating coordinates, we nevertheless find tentative evidence that approaching black hole formation the EdGB modifications cause the growth of scalar field "hair", consistent with known static black hole solutions in EdGB gravity. For the strong EdGB coupling regime, in a companion paper we first showed results that even in the weak field (i.e. far from black hole formation), the EdGB equations are of mixed type: evolution of the initially hyperbolic system of partial differential equations lead to formation of a region where their character changes to elliptic. Here, we present more details about this regime. In particular, we show that an effective energy density based on the Misner-Sharp mass is negative near these elliptic regions, and similarly the null convergence condition is violated then.
△ Less
Submitted 15 October, 2019; v1 submitted 18 March, 2019;
originally announced March 2019.
-
Hyperbolicity in Spherical Gravitational Collapse in a Horndeski Theory
Authors:
Justin L Ripley,
Frans Pretorius
Abstract:
We numerically study spherical gravitational collapse in shift symmetric Einstein dilaton Gauss Bonnet (EdGB) gravity. We find evidence that there are open sets of initial data for which the character of the system of equations changes from hyperbolic to elliptic type in a compact region of the spacetime. In these cases evolution of the system, treated as a hyperbolic initial boundary value proble…
▽ More
We numerically study spherical gravitational collapse in shift symmetric Einstein dilaton Gauss Bonnet (EdGB) gravity. We find evidence that there are open sets of initial data for which the character of the system of equations changes from hyperbolic to elliptic type in a compact region of the spacetime. In these cases evolution of the system, treated as a hyperbolic initial boundary value problem, leads to the equations of motion becoming ill-posed when the elliptic region forms. No singularities or discontinuities are encountered on the corresponding effective "Cauchy horizon". Therefore it is conceivable that a well-posed formulation of EdGB gravity (at least within spherical symmetry) may be possible if the equations are appropriately treated as mixed-type.
△ Less
Submitted 15 October, 2019; v1 submitted 4 February, 2019;
originally announced February 2019.
-
Stability and the Gauge Problem in Non-Perturbative Cosmology
Authors:
Anna Ijjas,
Frans Pretorius,
Paul J. Steinhardt
Abstract:
In this paper, we describe the first steps towards fully non-perturbative cosmology. We explain why the conventional methods used by cosmologists based on the ADM formulation are generally inadequate for this purpose and why it is advantageous instead to adapt the harmonic formulation pioneered and utilized in mathematical and numerical relativity. Here we focus on using this approach to evaluatin…
▽ More
In this paper, we describe the first steps towards fully non-perturbative cosmology. We explain why the conventional methods used by cosmologists based on the ADM formulation are generally inadequate for this purpose and why it is advantageous instead to adapt the harmonic formulation pioneered and utilized in mathematical and numerical relativity. Here we focus on using this approach to evaluating the linear mode stability in homogeneous and nearly homogeneous backgrounds and devising a valid scheme and diagnostics for numerical computation. We also briefly touch on the relevance of these methods for extracting cosmological observables from non-perturbative simulations.
△ Less
Submitted 19 September, 2018;
originally announced September 2018.
-
Evolution of Highly Eccentric Binary Neutron Stars Including Tidal Effects
Authors:
Huan Yang,
William E. East,
Vasileios Paschalidis,
Frans Pretorius,
Raissa F. P. Mendes
Abstract:
This work is the first in a series of studies aimed at understanding the dynamics of highly eccentric binary neutron stars, and constructing an appropriate gravitational-waveform model for detection. Such binaries are possible sources for ground-based gravitational wave detectors, and are expected to form through dynamical scattering and multi-body interactions in globular clusters and galactic nu…
▽ More
This work is the first in a series of studies aimed at understanding the dynamics of highly eccentric binary neutron stars, and constructing an appropriate gravitational-waveform model for detection. Such binaries are possible sources for ground-based gravitational wave detectors, and are expected to form through dynamical scattering and multi-body interactions in globular clusters and galactic nuclei. In contrast to black holes, oscillations of neutron stars are generically excited by tidal effects after close pericenter passage. Depending on the equation of state, this can enhance the loss of orbital energy by up to tens of percent over that radiated away by gravitational waves during an orbit. Under the same interaction mechanism, part of the orbital angular momentum is also transferred to the star. We calculate the impact of the neutron star oscillations on the orbital evolution of such systems, and compare these results to full numerical simulations. Utilizing a Post-Newtonian flux description we propose a preliminary model to predict the timing of different pericenter passages. A refined version of this model (taking into account Post-Newtonian corrections to the tidal coupling and the oscillations of the stars) may serve as a waveform model for such highly eccentric systems.
△ Less
Submitted 31 May, 2018;
originally announced June 2018.
-
Superradiant Instability and Backreaction of Massive Vector Fields around Kerr Black Holes
Authors:
William E. East,
Frans Pretorius
Abstract:
We study the growth and saturation of the superradiant instability of a complex, massive vector (Proca) field as it extracts energy and angular momentum from a spinning black hole, using numerical solutions of the full Einstein-Proca equations. We concentrate on a rapidly spinning black hole ($a=0.99$) and the dominant $m=1$ azimuthal mode of the Proca field, with real and imaginary components of…
▽ More
We study the growth and saturation of the superradiant instability of a complex, massive vector (Proca) field as it extracts energy and angular momentum from a spinning black hole, using numerical solutions of the full Einstein-Proca equations. We concentrate on a rapidly spinning black hole ($a=0.99$) and the dominant $m=1$ azimuthal mode of the Proca field, with real and imaginary components of the field chosen to yield an axisymmetric stress-energy tensor and, hence, spacetime. We find that in excess of $9\%$ of the black hole's mass can be transferred into the field. In all cases studied, the superradiant instability smoothly saturates when the black hole's horizon frequency decreases to match the frequency of the Proca cloud that spontaneously forms around the black hole.
△ Less
Submitted 14 June, 2017; v1 submitted 16 April, 2017;
originally announced April 2017.
-
Equation of state effects and one-arm spiral instability in hypermassive neutron stars formed in eccentric neutron star mergers
Authors:
William E. East,
Vasileios Paschalidis,
Frans Pretorius
Abstract:
We continue our investigations of the development and importance of the one-arm spiral instability in long-lived hypermassive neutron stars (HMNSs) formed in dynamical capture binary neutron star mergers. Employing hydrodynamic simulations in full general relativity, we find that the one-arm instability is generic in that it can develop in HMNSs within a few tens of milliseconds after merger for a…
▽ More
We continue our investigations of the development and importance of the one-arm spiral instability in long-lived hypermassive neutron stars (HMNSs) formed in dynamical capture binary neutron star mergers. Employing hydrodynamic simulations in full general relativity, we find that the one-arm instability is generic in that it can develop in HMNSs within a few tens of milliseconds after merger for all equations of state in our survey. We find that mergers with stiffer equations of state tend to produce HMNSs with stronger $m=2$ azimuthal mode density deformations, and weaker $m=1$ components, relative to softer equations of state. We also find that for equations of state that can give rise to double-core HMNSs, large $m=1$ density modes can already be present due to asymmetries in the two cores. This results in the generation of $l=2$, $m=1$ gravitational wave modes even before the dominance of a one-arm mode that ultimately arises following merger of the two cores. Our results further suggest that stiffer equations of state give rise to HMNSs generating lower $m=1$ gravitational wave frequencies. Thus, if gravitational waves from the one-arm instability are detected, they could in principle constrain the neutron star equation of state. We estimate that, depending on the equation of state, the one-arm mode could potentially be detectable by second generation gravitational wave detectors at $\sim 10$ Mpc and by third generation ones at $\sim 100$ Mpc. Finally, we provide estimates of the properties of dynamical ejecta, as well as the accompanying kilonovae signatures.
△ Less
Submitted 2 September, 2016;
originally announced September 2016.
-
Theoretical Physics Implications of the Binary Black-Hole Mergers GW150914 and GW151226
Authors:
Nicolas Yunes,
Kent Yagi,
Frans Pretorius
Abstract:
The gravitational wave observations GW150914 and GW151226 by Advanced LIGO provide the first opportunity to learn about physics in the extreme gravity environment of coalescing binary black holes. The LIGO Scientific Collaboration and the Virgo Collaboration have verified that this observation is consistent with General Relativity. This paper expands their analysis to a larger class of anomalies,…
▽ More
The gravitational wave observations GW150914 and GW151226 by Advanced LIGO provide the first opportunity to learn about physics in the extreme gravity environment of coalescing binary black holes. The LIGO Scientific Collaboration and the Virgo Collaboration have verified that this observation is consistent with General Relativity. This paper expands their analysis to a larger class of anomalies, highlighting the inferences that can be drawn on non-standard theoretical physics mechanisms. We find that these events constrain a plethora of mechanisms associated with the generation and propagation of gravitational waves, including the activation of scalar fields, gravitational leakage into large extra dimensions, the variability of Newton's constant, a modified dispersion relation, gravitational Lorentz violation and the strong equivalence principle. Though other observations limit many of these mechanisms already, GW150914 and GW151226 are unique in that they are direct probes of dynamical strong-field gravity and of gravitational wave propagation. We also show that GW150914 constrains inferred properties of exotic compact object alternatives to Kerr black holes. We argue, however, that the true potential for GW150914 to both rule out exotic objects and constrain physics beyond General Relativity is severely limited by the lack of understanding of the merger regime in almost all relevant modified gravity theories. This event thus significantly raises the bar that these theories have to pass, both in terms of having a sound theoretical underpinning, and being able to solve the equations of motion for binary merger events. We conclude with a discussion of the additional inferences that can be drawn if the lower-confidence observation of an electromagnetic counterpart to GW150914 holds true; this would provide dramatic constraints on the speed of gravity and gravitational Lorentz violation.
△ Less
Submitted 8 September, 2016; v1 submitted 29 March, 2016;
originally announced March 2016.
-
Relativistic Simulations of Eccentric Binary Neutron Star Mergers: One-arm Spiral Instability and Effects of Neutron Star Spin
Authors:
William E. East,
Vasileios Paschalidis,
Frans Pretorius,
Stuart L. Shapiro
Abstract:
We perform general-relativistic hydrodynamical simulations of dynamical capture binary neutron star mergers, emphasizing the role played by the neutron star spin. Dynamical capture mergers may take place in globular clusters, as well as other dense stellar systems, where most neutron stars have large spins. We find significant variability in the merger outcome as a function of initial neutron star…
▽ More
We perform general-relativistic hydrodynamical simulations of dynamical capture binary neutron star mergers, emphasizing the role played by the neutron star spin. Dynamical capture mergers may take place in globular clusters, as well as other dense stellar systems, where most neutron stars have large spins. We find significant variability in the merger outcome as a function of initial neutron star spin. For cases where the spin is aligned with the orbital angular momentum, the additional centrifugal support in the remnant hypermassive neutron star can prevent the prompt collapse to a black hole, while for antialigned cases the decreased total angular momentum can facilitate the collapse to a black hole. We show that even moderate spins can significantly increase the amount of ejected material, including the amount unbound with velocities greater than half the speed of light, leading to brighter electromagnetic signatures associated with kilonovae and interaction of the ejecta with the interstellar medium. Furthermore, we find that the initial neutron star spin can strongly affect the already rich phenomenology in the postmerger gravitational wave signatures that arise from the oscillation modes of the hypermassive neutron star. In several of our simulations, the resulting hypermassive neutron star develops the one-arm ($m=1$) spiral instability, the most pronounced cases being those with small but non-negligible neutron star spins. For long-lived hypermassive neutron stars, the presence of this instability leads to improved prospects for detecting these events through gravitational waves, and thus may give information about the neutron star equation of state.
△ Less
Submitted 13 January, 2016; v1 submitted 3 November, 2015;
originally announced November 2015.
-
One-arm Spiral Instability in Hypermassive Neutron Stars Formed by Dynamical-Capture Binary Neutron Star Mergers
Authors:
Vasileios Paschalidis,
William E. East,
Frans Pretorius,
Stuart L. Shapiro
Abstract:
Using general-relativistic hydrodynamical simulations, we show that merging binary neutron stars can form hypermassive neutrons stars that undergo the one-arm spiral instability. We study the particular case of a dynamical capture merger where the stars have a small spin, as may arise in globular clusters, and focus on an equal-mass scenario where the spins are aligned with the orbital angular mom…
▽ More
Using general-relativistic hydrodynamical simulations, we show that merging binary neutron stars can form hypermassive neutrons stars that undergo the one-arm spiral instability. We study the particular case of a dynamical capture merger where the stars have a small spin, as may arise in globular clusters, and focus on an equal-mass scenario where the spins are aligned with the orbital angular momentum. We find that this instability develops when post-merger fluid vortices lead to the generation of a toroidal remnant - a configuration whose maximum density occurs in a ring around the center-of-mass - with high vorticity along its rotation axis. The instability quickly saturates on a timescale of $\sim 10$ ms, with the $m=1$ azimuthal density multipole mode dominating over higher modes. The instability also leaves a characteristic imprint on the post-merger gravitational wave signal that could be detectable if the instability persists in long-lived remnants.
△ Less
Submitted 29 October, 2015; v1 submitted 12 October, 2015;
originally announced October 2015.
-
Eccentric mergers of black holes with spinning neutron stars
Authors:
William E. East,
Vasileios Paschalidis,
Frans Pretorius
Abstract:
We study dynamical capture binary black hole-neutron star (BH-NS) mergers focusing on the effects of the neutron star spin. These events may arise in dense stellar regions, such as globular clusters, where the majority of neutron stars are expected to be rapidly rotating. We initialize the BH-NS systems with positions and velocities corresponding to marginally unbound Newtonian orbits, and evolve…
▽ More
We study dynamical capture binary black hole-neutron star (BH-NS) mergers focusing on the effects of the neutron star spin. These events may arise in dense stellar regions, such as globular clusters, where the majority of neutron stars are expected to be rapidly rotating. We initialize the BH-NS systems with positions and velocities corresponding to marginally unbound Newtonian orbits, and evolve them using general-relativistic hydrodynamical simulations. We find that even moderate spins can significantly increase the amount of mass in unbound material. In some of the more extreme cases, there can be up to a third of a solar mass in unbound matter. Similarly, large amounts of tidally stripped material can remain bound and eventually accrete onto the BH---as much as a tenth of a solar mass in some cases. These simulations demonstrate that it is important to treat neutron star spin in order to make reliable predictions of the gravitational wave and electromagnetic transient signals accompanying these sources.
△ Less
Submitted 25 June, 2015; v1 submitted 24 March, 2015;
originally announced March 2015.
-
Probing Strong Field Gravity Through Numerical Simulations
Authors:
Matthew W. Choptuik,
Luis Lehner,
Frans Pretorius
Abstract:
This article is an overview of the contributions numerical relativity has made to our understanding of strong field gravity, to be published in the book "General Relativity and Gravitation: A Centennial Perspective", commemorating the 100th anniversary of general relativity.
This article is an overview of the contributions numerical relativity has made to our understanding of strong field gravity, to be published in the book "General Relativity and Gravitation: A Centennial Perspective", commemorating the 100th anniversary of general relativity.
△ Less
Submitted 24 February, 2015;
originally announced February 2015.
-
Numerical Relativity and Astrophysics
Authors:
Luis Lehner,
Frans Pretorius
Abstract:
Throughout the Universe many powerful events are driven by strong gravitational effects that require general relativity to fully describe them. These include compact binary mergers, black hole accretion and stellar collapse, where velocities can approach the speed of light, and extreme gravitational fields --$Φ_{\rm Newt}/c^2 \simeq 1$-- mediate the interactions. Many of these processes trigger em…
▽ More
Throughout the Universe many powerful events are driven by strong gravitational effects that require general relativity to fully describe them. These include compact binary mergers, black hole accretion and stellar collapse, where velocities can approach the speed of light, and extreme gravitational fields --$Φ_{\rm Newt}/c^2 \simeq 1$-- mediate the interactions. Many of these processes trigger emission across a broad range of the electromagnetic spectrum. Compact binaries further source strong gravitational wave emission that could directly be detected in the near future. This feat will open up a gravitational wave window into our Universe and revolutionize its understanding. Describing these phenomena requires general relativity, and --where dynamical effects strongly modify gravitational fields-- the full Einstein equations coupled to matter sources. Numerical relativity is a field within general relativity concerned with studying such scenarios that cannot be accurately modeled via perturbative or analytical calculations. In this review, we examine results obtained within this discipline, with a focus on its impact in astrophysics.
△ Less
Submitted 19 May, 2014;
originally announced May 2014.
-
A Parametrized post-Einsteinian Framework for Gravitational Wave Bursts
Authors:
Nicholas Loutrel,
Nicolas Yunes,
Frans Pretorius
Abstract:
The population of stellar-mass, compact object binaries that merge with non-negligible eccentricity may be large enough to motivate searches with ground-based gravitational wave detectors. Such events could be exceptional laboratories to test General Relativity in the dynamical, strong-field regime, as a larger fraction of the energy is emitted at high-velocities, compared to quasi-circular inspir…
▽ More
The population of stellar-mass, compact object binaries that merge with non-negligible eccentricity may be large enough to motivate searches with ground-based gravitational wave detectors. Such events could be exceptional laboratories to test General Relativity in the dynamical, strong-field regime, as a larger fraction of the energy is emitted at high-velocities, compared to quasi-circular inspirals. A serious obstacle here, however, is the challenge of computing theoretical waveforms for eccentric systems with the requisite accuracy for use in a matched-filter search. The corresponding waveforms are more a sequence of concentrated bursts of energy emitted near periapse than a continuous waveform. Based on this, an alternative approach, stacking excess power over the set of time-frequency tiles coincident with the bursts, was recently suggested as a more practical (though sub-optimal) detection strategy. The leading-order "observable" that would be inferred from such a detection would be a sequence of discrete numbers characterizing the position and size of each time-frequency tile. In General Relativity, this (possibly large) sequence of numbers is uniquely determined by the small set of parameters describing the binary at formation. In this work, following the spirit of the parameterized post-Einsteinian framework developed for quasi-circular inspiral, we propose a simple, parameterized deformation of the baseline general relativistic burst algorithm for eccentric inspiral events that would allow for model-independent tests of Einstein's theory in this high-velocity, strong-field regime.
△ Less
Submitted 2 October, 2017; v1 submitted 31 March, 2014;
originally announced April 2014.
-
Black Hole Superradiance in Dynamical Spacetime
Authors:
William E. East,
Fethi M. Ramazanoğlu,
Frans Pretorius
Abstract:
We study the superradiant scattering of gravitational waves by a nearly extremal black hole (dimensionless spin $a=0.99$) by numerically solving the full Einstein field equations, thus including backreaction effects. This allows us to study the dynamics of the black hole as it loses energy and angular momentum during the scattering process. To explore the nonlinear phase of the interaction, we con…
▽ More
We study the superradiant scattering of gravitational waves by a nearly extremal black hole (dimensionless spin $a=0.99$) by numerically solving the full Einstein field equations, thus including backreaction effects. This allows us to study the dynamics of the black hole as it loses energy and angular momentum during the scattering process. To explore the nonlinear phase of the interaction, we consider gravitational wave packets with initial energies up to $10%$ of the mass of the black hole. We find that as the incident wave energy increases, the amplification of the scattered waves, as well as the energy extraction efficiency from the black hole, is reduced. During the interaction the apparent horizon geometry undergoes sizable nonaxisymmetric oscillations. The largest amplitude excitations occur when the peak frequency of the incident wave packet is above where superradiance occurs, but close to the dominant quasinormal mode frequency of the black hole.
△ Less
Submitted 14 March, 2014; v1 submitted 16 December, 2013;
originally announced December 2013.
-
Nonperturbative analysis of the evolution of cosmological perturbations through a nonsingular bounce
Authors:
BingKan Xue,
David Garfinkle,
Frans Pretorius,
Paul J. Steinhardt
Abstract:
In bouncing cosmology, the primordial fluctuations are generated in a cosmic contraction phase before the bounce into the current expansion phase. For a nonsingular bounce, curvature and anisotropy grow rapidly during the bouncing phase, raising questions about the reliability of perturbative analysis. In this paper, we study the evolution of adiabatic perturbations in a nonsingular bounce by nonp…
▽ More
In bouncing cosmology, the primordial fluctuations are generated in a cosmic contraction phase before the bounce into the current expansion phase. For a nonsingular bounce, curvature and anisotropy grow rapidly during the bouncing phase, raising questions about the reliability of perturbative analysis. In this paper, we study the evolution of adiabatic perturbations in a nonsingular bounce by nonperturbative methods including numerical simulations of the nonsingular bounce and the covariant formalism for calculating nonlinear perturbations. We show that the bounce is disrupted in regions of the universe with significant inhomogeneity and anisotropy over the background energy density, but is achieved in regions that are relatively homogeneous and isotropic. Sufficiently small perturbations, consistent with observational constraints, can pass through the nonsingular bounce with negligible alteration from nonlinearity. We follow scale invariant perturbations generated in a matter-like contraction phase through the bounce. Their amplitude in the expansion phase is determined by the growing mode in the contraction phase, and the scale invariance is well preserved across the bounce.
△ Less
Submitted 14 August, 2013;
originally announced August 2013.
-
Simulating extreme-mass-ratio systems in full general relativity
Authors:
William E. East,
Frans Pretorius
Abstract:
We introduce a new method for numerically evolving the full Einstein field equations in situations where the spacetime is dominated by a known background solution. The technique leverages the knowledge of the background solution to subtract off its contribution to the truncation error, thereby more efficiently achieving a desired level of accuracy. We demonstrate the method by applying it to the r…
▽ More
We introduce a new method for numerically evolving the full Einstein field equations in situations where the spacetime is dominated by a known background solution. The technique leverages the knowledge of the background solution to subtract off its contribution to the truncation error, thereby more efficiently achieving a desired level of accuracy. We demonstrate the method by applying it to the radial infall of a solar-type star into supermassive black holes with mass ratios $\geq 10^6$. The self-gravity of the star is thus consistently modeled within the context of general relativity, and the star's interaction with the black hole computed with moderate computational cost, despite the over five orders of magnitude difference in gravitational potential (as defined by the ratio of mass to radius). We compute the tidal deformation of the star during infall, and the gravitational wave emission, finding the latter is close to the prediction of the point-particle limit.
△ Less
Submitted 12 May, 2013; v1 submitted 6 March, 2013;
originally announced March 2013.
-
Observing complete gravitational wave signals from dynamical capture binaries
Authors:
William E. East,
Sean T. McWilliams,
Janna Levin,
Frans Pretorius
Abstract:
We assess the detectability of the gravitational wave signals from highly eccentric compact binaries. We use a simple model for the inspiral, merger, and ringdown of these systems. The model is based on mapping the binary to an effective single black hole system described by a Kerr metric, thereby including certain relativistic effects such as zoom-whirl-type behavior. The resultant geodesics sour…
▽ More
We assess the detectability of the gravitational wave signals from highly eccentric compact binaries. We use a simple model for the inspiral, merger, and ringdown of these systems. The model is based on mapping the binary to an effective single black hole system described by a Kerr metric, thereby including certain relativistic effects such as zoom-whirl-type behavior. The resultant geodesics source quadrupolar radiation and, in turn, are evolved under its dissipative effects. At the light ring, we attach a merger model that was previously developed for quasicircular mergers but also performs well for eccentric mergers with little modification. We apply this model to determine the detectability of these sources for initial, Enhanced, and Advanced LIGO across the parameter space of nonspinning close capture compact binaries. We conclude that, should these systems exist in nature, the vast majority will be missed by conventional burst searches or by quasicircular waveform templates in the advanced detector era. Other methods, such as eccentric templates or, more practically, a stacked excess power search, must be developed to avoid losing these sources. These systems would also have been missed frequently in the initial LIGO data analysis. Thus, previous null coincidence results with detected gamma-ray bursts cannot exclude the possibility of coincident gravitational wave signals from eccentric binaries.
△ Less
Submitted 11 February, 2013; v1 submitted 4 December, 2012;
originally announced December 2012.
-
Gravitational waves and stalled satellites from massive galaxy mergers at z <= 1
Authors:
Sean T. McWilliams,
Jeremiah P. Ostriker,
Frans Pretorius
Abstract:
We present a model for merger-driven evolution of the mass function for massive galaxies and their central supermassive black holes at late times. We discuss the current observational evidence in favor of merger-driven massive galaxy evolution during this epoch, and demonstrate that the observed evolution of the mass function can be reproduced by evolving an initial mass function under the assumpt…
▽ More
We present a model for merger-driven evolution of the mass function for massive galaxies and their central supermassive black holes at late times. We discuss the current observational evidence in favor of merger-driven massive galaxy evolution during this epoch, and demonstrate that the observed evolution of the mass function can be reproduced by evolving an initial mass function under the assumption of negligible star formation. We calculate the stochastic gravitational wave signal from the resulting black-hole binary mergers in the low redshift universe (z <= 1) implied by this model, and find that this population has a signal-to-noise ratio as much as ~5x larger than previous estimates for pulsar timing arrays, with an expectation value for the characteristic strain h_c (f=1 yr^{-1}) = 4.1 x 10^{-15} that may already be in tension with observational constraints, and a {2-sigma, 3-sigma} lower limit within this model of h_c (f=1 yr^{-1}) = {1.1 x 10^{-15}, 6.8 x 10^{-16}}. The strength of this signal is sufficient to make it detectable with high probability under conservative assumptions within the next several years, if the principle assumption of merger-driven galaxy evolution since z = 1 holds true. For cases where a galaxy merger fails to lead to a black hole merger, we estimate the probability for a given number of satellite unmerged black holes to remain within a massive host galaxy, and interpret the result in light of ULX observations. In particular, we find that the brightest cluster galaxies should have 1-2 such sources with luminosities above 10^{39} erg/s, which is consistent with the statistics of observed ULXs.
△ Less
Submitted 6 February, 2013; v1 submitted 22 November, 2012;
originally announced November 2012.
-
The imminent detection of gravitational waves from massive black-hole binaries with pulsar timing arrays
Authors:
Sean T. McWilliams,
Jeremiah P. Ostriker,
Frans Pretorius
Abstract:
Recent observations of massive galaxies indicate that they double in mass and quintuple in size between redshift z = 1 and the present, despite undergoing very little star formation, suggesting that galaxy mergers drive the evolution. Since these galaxies will contain supermassive black holes, this suggests a larger black hole merger rate, and therefore a larger gravitational-wave signal, than pre…
▽ More
Recent observations of massive galaxies indicate that they double in mass and quintuple in size between redshift z = 1 and the present, despite undergoing very little star formation, suggesting that galaxy mergers drive the evolution. Since these galaxies will contain supermassive black holes, this suggests a larger black hole merger rate, and therefore a larger gravitational-wave signal, than previously expected. We calculate the merger-driven evolution of the mass function, and find that merger rates are 10 to 30 times higher and gravitational waves are 3 to 5 times stronger than previously estimated, so that the gravitational-wave signal may already be detectable with existing data from pulsar timing arrays. We also provide an explanation for the disagreement with past estimates that were based on dark matter halo simulations.
△ Less
Submitted 19 November, 2012;
originally announced November 2012.
-
Ultrarelativistic black hole formation
Authors:
William E. East,
Frans Pretorius
Abstract:
We study the ultrarelativistic head-on collision of equal mass particles, modeled as self-gravitating fluid spheres, by numerically solving the coupled Einstein-hydrodynamic equations. We focus on cases well within the kinetic energy dominated regime, where between 88-92% ($γ=8$ to 12) of the initial net energy of the spacetime resides in the translation kinetic energy of the particles. We find th…
▽ More
We study the ultrarelativistic head-on collision of equal mass particles, modeled as self-gravitating fluid spheres, by numerically solving the coupled Einstein-hydrodynamic equations. We focus on cases well within the kinetic energy dominated regime, where between 88-92% ($γ=8$ to 12) of the initial net energy of the spacetime resides in the translation kinetic energy of the particles. We find that for sufficiently large boosts, black hole formation occurs. Moreover, near yet above the threshold of black hole formation, the collision initially leads to the formation of two distinct apparent horizons that subsequently merge. We argue that this can be understood in terms of a focusing effect, where one boosted particle acts as a gravitational lens on the other and vice versa, and that this is further responsible for the threshold being lower (by a factor of a few) compared to simple hoop conjecture estimates. Cases slightly below threshold result in complete disruption of the model particles. The gravitational radiation emitted when black holes form reaches luminosities of 0.014 $c^5/G$, carrying $16\pm2%$ of the total energy.
△ Less
Submitted 20 February, 2013; v1 submitted 1 October, 2012;
originally announced October 2012.
-
Dynamical Capture Binary Neutron Star Mergers
Authors:
William E. East,
Frans Pretorius
Abstract:
We study dynamical capture binary neutron star mergers as may arise in dense stellar regions such as globular clusters. Using general-relativistic hydrodynamics, we find that these mergers can result in the prompt collapse to a black hole or in the formation of a hypermassive neutron star, depending not only on the neutron star equation of state but also on impact parameter. We also find that thes…
▽ More
We study dynamical capture binary neutron star mergers as may arise in dense stellar regions such as globular clusters. Using general-relativistic hydrodynamics, we find that these mergers can result in the prompt collapse to a black hole or in the formation of a hypermassive neutron star, depending not only on the neutron star equation of state but also on impact parameter. We also find that these mergers can produce accretion disks of up to a tenth of a solar mass and unbound ejected material of up to a few percent of a solar mass. We comment on the gravitational radiation and electromagnetic transients that these sources may produce.
△ Less
Submitted 31 October, 2012; v1 submitted 26 August, 2012;
originally announced August 2012.
-
Detectable seismic consequences of the interaction of a primordial black hole with Earth
Authors:
Yang Luo,
Shravan Hanasoge,
Jeroen Tromp,
Frans Pretorius
Abstract:
Galaxies observed today are likely to have evolved from density perturbations in the early universe. Perturbations that exceeded some critical threshold are conjectured to have undergone gravitational collapse to form primordial black holes (PBHs) at a range of masses. Such PBHs serve as candidates for cold dark matter and their detection would shed light on conditions in the early universe. Here…
▽ More
Galaxies observed today are likely to have evolved from density perturbations in the early universe. Perturbations that exceeded some critical threshold are conjectured to have undergone gravitational collapse to form primordial black holes (PBHs) at a range of masses. Such PBHs serve as candidates for cold dark matter and their detection would shed light on conditions in the early universe. Here we propose a mechanism to search for transits of PBHs through/nearby Earth by studying the associated seismic waves. Using a spectral-element method, we simulate and visualize this seismic wave field in Earth's interior. We predict the emergence of two unique signatures, namely, a wave that would arrive almost simultaneously everywhere on Earth's free surface and the excitation of unusual spheroidal modes with a characteristic frequency-spacing in free oscillation spectra. These qualitative characteristics are unaffected by the speed or proximity of the PBH trajectory. The seismic energy deposited by a proximal ${M^{PBH} = 10^{15}}$ g PBH is comparable to a magnitude $M_w=4$ earthquake. The non-seismic collateral damage due to the actual impact of such small PBHs with Earth would be negligible. Unfortunately, the expected collision rate is very low even if PBHs constituted all of dark matter, at ${\sim 10^{-7} {yr}^{-1}}$, and since the rate scales as ${1/M^{PBH}}$, fortunately encounters with larger, Earth-threatening PBHs are exceedingly unlikely. However, the rate at which non-colliding close encounters of PBHs could be detected by seismic activity alone is roughly two orders of magnitude larger --- that is once every hundred thousand years --- than the direct collision rate.
△ Less
Submitted 16 March, 2012;
originally announced March 2012.
-
NR/HEP: roadmap for the future
Authors:
Vitor Cardoso,
Leonardo Gualtieri,
Carlos Herdeiro,
Ulrich Sperhake,
Paul M. Chesler,
Luis Lehner,
Seong Chan Park,
Harvey S. Reall,
Carlos F. Sopuerta,
Daniela Alic,
Oscar J. C. Dias,
Roberto Emparan,
Valeria Ferrari,
Steven B. Giddings,
Mahdi Godazgar,
Ruth Gregory,
Veronika E. Hubeny,
Akihiro Ishibashi,
Greg Landsberg,
Carlos O. Lousto,
David Mateos,
Vicki Moeller,
Hirotada Okawa,
Paolo Pani,
M. Andy Parker
, et al. (7 additional authors not shown)
Abstract:
Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong fi…
▽ More
Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.
△ Less
Submitted 24 January, 2012;
originally announced January 2012.
-
Hydrodynamics in full general relativity with conservative AMR
Authors:
William E. East,
Frans Pretorius,
Branson C. Stephens
Abstract:
There is great interest in numerical relativity simulations involving matter due to the likelihood that binary compact objects involving neutron stars will be detected by gravitational wave observatories in the coming years, as well as to the possibility that binary compact object mergers could explain short-duration gamma-ray bursts. We present a code designed for simulations of hydrodynamics cou…
▽ More
There is great interest in numerical relativity simulations involving matter due to the likelihood that binary compact objects involving neutron stars will be detected by gravitational wave observatories in the coming years, as well as to the possibility that binary compact object mergers could explain short-duration gamma-ray bursts. We present a code designed for simulations of hydrodynamics coupled to the Einstein field equations targeted toward such applications. This code has recently been used to study eccentric mergers of black hole-neutron star binaries. We evolve the fluid conservatively using high-resolution shock-capturing methods, while the field equations are solved in the generalized-harmonic formulation with finite differences. In order to resolve the various scales that may arise, we use adaptive mesh refinement (AMR) with grid hierarchies based on truncation error estimates. A noteworthy feature of this code is the implementation of the flux correction algorithm of Berger and Colella to ensure that the conservative nature of fluid advection is respected across AMR boundaries. We present various tests to compare the performance of different limiters and flux calculation methods, as well as to demonstrate the utility of AMR flux corrections.
△ Less
Submitted 29 May, 2012; v1 submitted 13 December, 2011;
originally announced December 2011.
-
Eccentric black hole-neutron star mergers: effects of black hole spin and equation of state
Authors:
William E. East,
Frans Pretorius,
Branson C. Stephens
Abstract:
There is a high level of interest in black hole-neutron star binaries, not only because their mergers may be detected by gravitational wave observatories in the coming years, but also because of the possibility that they could explain a class of short duration gamma-ray bursts. We study black hole-neutron star mergers that occur with high eccentricity as may arise from dynamical capture in dense s…
▽ More
There is a high level of interest in black hole-neutron star binaries, not only because their mergers may be detected by gravitational wave observatories in the coming years, but also because of the possibility that they could explain a class of short duration gamma-ray bursts. We study black hole-neutron star mergers that occur with high eccentricity as may arise from dynamical capture in dense stellar regions such as nuclear or globular clusters. We perform general relativistic simulations of binaries with a range of impact parameters, three different initial black hole spins (zero, aligned and anti-aligned with the orbital angular momentum), and neutron stars with three different equations of state. We find a rich diversity across these parameters in the resulting gravitational wave signals and matter dynamics, which should also be reflected in the consequent electromagnetic emission. Before tidal disruption, the gravitational wave emission is significantly larger than perturbative predictions suggest for periapsis distances close to effective innermost stable separations, exhibiting features reflecting the zoom-whirl dynamics of the orbit there. Guided by the simulations, we develop a simple model for the change in orbital parameters of the binary during close encounters. Depending upon the initial parameters of the system, we find that mass transfer during non-merging close encounters can range from essentially zero to a sizable fraction of the initial neutron star mass. The same holds for the amount of material outside the black hole post-merger, and in some cases roughly half of this material is estimated to be unbound. We also see that non-merging close encounters generically excite large oscillations in the neutron star that are qualitatively consistent with f-modes.
△ Less
Submitted 18 April, 2012; v1 submitted 13 November, 2011;
originally announced November 2011.
-
Eccentric binary neutron star mergers
Authors:
Roman Gold,
Sebastiano Bernuzzi,
Marcus Thierfelder,
Bernd Bruegmann,
Frans Pretorius
Abstract:
Neutron star binaries offer a rich phenomenology in terms of gravitational waves and merger remnants. However, most general relativistic studies have been performed for nearly circular binaries, with the exception of head-on collisions. We present the first numerical relativity investigation of mergers of eccentric equal-mass neutron-star binaries that probes the regime between head-on and circula…
▽ More
Neutron star binaries offer a rich phenomenology in terms of gravitational waves and merger remnants. However, most general relativistic studies have been performed for nearly circular binaries, with the exception of head-on collisions. We present the first numerical relativity investigation of mergers of eccentric equal-mass neutron-star binaries that probes the regime between head-on and circular. In addition to gravitational waves generated by the orbital motion, we find that the signal also contains a strong component due to stellar oscillations (f-modes) induced by tidal forces, extending a classical result for Newtonian binaries. The merger can lead to rather massive disks on the order of 10% of the total initial mass.
△ Less
Submitted 8 January, 2013; v1 submitted 23 September, 2011;
originally announced September 2011.
-
Eccentric Black Hole-Neutron Star Mergers
Authors:
Branson C. Stephens,
William E. East,
Frans Pretorius
Abstract:
Within the next few years gravitational waves (GWs) from merging black holes (BHs) and neutron stars (NSs) may be directly detected, making a thorough theoretical understanding of these systems a high priority. As an additional motivation, these systems may represent a subset of short-duration gamma-ray burst (sGRB) progenitors. BH-NS mergers are expected to result from primordial, quasi-circular…
▽ More
Within the next few years gravitational waves (GWs) from merging black holes (BHs) and neutron stars (NSs) may be directly detected, making a thorough theoretical understanding of these systems a high priority. As an additional motivation, these systems may represent a subset of short-duration gamma-ray burst (sGRB) progenitors. BH-NS mergers are expected to result from primordial, quasi-circular inspiral as well as dynamically formed capture binaries. The latter channel allows mergers with high eccentricity, resulting in a richer variety of outcomes. We perform general relativistic simulations of BH-NS interactions with a range of impact parameters, and find significant variation in the properties of these events that have potentially observable consequences, namely the GW signature, remnant accretion disk mass, and amount of unbound material.
△ Less
Submitted 22 July, 2011; v1 submitted 16 May, 2011;
originally announced May 2011.
-
Gravitational Wave Tests of General Relativity with the Parameterized Post-Einsteinian Framework
Authors:
Neil Cornish,
Laura Sampson,
Nico Yunes,
Frans Pretorius
Abstract:
Gravitational wave astronomy has tremendous potential for studying extreme astrophysical phenomena and exploring fundamental physics. The waves produced by binary black hole mergers will provide a pristine environment in which to study strong field, dynamical gravity. Extracting detailed information about these systems requires accurate theoretical models of the gravitational wave signals. If grav…
▽ More
Gravitational wave astronomy has tremendous potential for studying extreme astrophysical phenomena and exploring fundamental physics. The waves produced by binary black hole mergers will provide a pristine environment in which to study strong field, dynamical gravity. Extracting detailed information about these systems requires accurate theoretical models of the gravitational wave signals. If gravity is not described by General Relativity, analyses that are based on waveforms derived from Einstein's field equations could result in parameter biases and a loss of detection efficiency. A new class of "parameterized post-Einsteinian" (ppE) waveforms has been proposed to cover this eventuality. Here we apply the ppE approach to simulated data from a network of advanced ground based interferometers (aLIGO/aVirgo) and from a future spaced based interferometer (LISA). Bayesian inference and model selection are used to investigate parameter biases, and to determine the level at which departures from general relativity can be detected. We find that in some cases the parameter biases from assuming the wrong theory can be severe. We also find that gravitational wave observations will beat the existing bounds on deviations from general relativity derived from the orbital decay of binary pulsars by a large margin across a wide swath of parameter space.
△ Less
Submitted 23 August, 2011; v1 submitted 10 May, 2011;
originally announced May 2011.
-
A Dynamical Systems Approach to Schwarzschild Null Geodesics
Authors:
Edward Belbruno,
Frans Pretorius
Abstract:
The null geodesics of a Schwarzschild black hole are studied from a dynamical systems perspective. Written in terms of Kerr-Schild coordinates, the null geodesic equation takes on the simple form of a particle moving under the influence of a Newtonian central force with an inverse-cubic potential. We apply a McGehee transformation to these equations, which clearly elucidates the full phase space o…
▽ More
The null geodesics of a Schwarzschild black hole are studied from a dynamical systems perspective. Written in terms of Kerr-Schild coordinates, the null geodesic equation takes on the simple form of a particle moving under the influence of a Newtonian central force with an inverse-cubic potential. We apply a McGehee transformation to these equations, which clearly elucidates the full phase space of solutions. All the null geodesics belong to one of four families of invariant manifolds and their limiting cases, further characterized by the angular momentum L of the orbit: for |L|>|L_c|, (1) the set that flow outward from the white hole, turn around, then fall into the black hole, (2) the set that fall inward from past null infinity, turn around outside the black hole to continue to future null infinity, and for |L|<|L_c|, (3) the set that flow outward from the white hole and continue to future null infinity, (4) the set that flow inward from past null infinity and into the black hole. The critical angular momentum Lc corresponds to the unstable circular orbit at r=3M, and the homoclinic orbits associated with it. There are two additional critical points of the flow at the singularity at r=0. Though the solutions of geodesic motion and Hamiltonian flow we describe here are well known, what we believe is a novel aspect of this work is the mapping between the two equivalent descriptions, and the different insights each approach can give to the problem. For example, the McGehee picture points to a particularly interesting limiting case of the class (1) that move from the white to black hole: in the limit as L goes to infinity, as described in Schwarzschild coordinates, these geodesics begin at r=0, flow along t=constant lines, turn around at r=2M, then continue to r=0. During this motion they circle in azimuth exactly once, and complete the journey in zero affine time.
△ Less
Submitted 2 March, 2011;
originally announced March 2011.
-
Linear Stability Analysis and the Speed of Gravitational Waves in Dynamical Chern-Simons Modified Gravity
Authors:
David Garfinkle,
Frans Pretorius,
Nicolas Yunes
Abstract:
We perform a linear stability analysis of dynamical Chern-Simons modified gravity in the geometric optics approximation and find that it is linearly stable on the backgrounds considered. Our analysis also reveals that gravitational waves in the modified theory travel at the speed of light in Minkowski spacetime. However, on a Schwarzschild background the characteristic speed of propagation along a…
▽ More
We perform a linear stability analysis of dynamical Chern-Simons modified gravity in the geometric optics approximation and find that it is linearly stable on the backgrounds considered. Our analysis also reveals that gravitational waves in the modified theory travel at the speed of light in Minkowski spacetime. However, on a Schwarzschild background the characteristic speed of propagation along a given direction splits into two modes, one subluminal and one superluminal. The width of the splitting depends on the azimuthal components of the propagation vector, is linearly proportional to the mass of the black hole, and decreases with the third inverse power of the distance from the black hole. Radial propagation is unaffected, implying that as probed by gravitational waves the location of the event horizon of the spacetime is unaltered. The analysis further reveals that when a high frequency, pure gravitational wave is scattered from a black hole, a scalar wave of comparable amplitude is excited, and vice-versa.
△ Less
Submitted 8 October, 2013; v1 submitted 14 July, 2010;
originally announced July 2010.
-
Constraining the evolutionary history of Newton's constant with gravitational wave observations
Authors:
Nicolas Yunes,
Frans Pretorius,
David Spergel
Abstract:
Space-borne gravitational wave detectors, such as the proposed Laser Interferometer Space Antenna, are expected to observe black hole coalescences to high redshift and with large signal-to-noise ratios, rendering their gravitational waves ideal probes of fundamental physics. The promotion of Newton's constant to a time-function introduces modifications to the binary's binding energy and the grav…
▽ More
Space-borne gravitational wave detectors, such as the proposed Laser Interferometer Space Antenna, are expected to observe black hole coalescences to high redshift and with large signal-to-noise ratios, rendering their gravitational waves ideal probes of fundamental physics. The promotion of Newton's constant to a time-function introduces modifications to the binary's binding energy and the gravitational wave luminosity, leading to corrections in the chirping frequency. Such corrections propagate into the response function and, given a gravitational wave observation, they allow for constraints on the first time-derivative of Newton's constant at the time of merger. We find that space-borne detectors could indeed place interesting constraints on this quantity as a function of sky position and redshift, providing a {\emph{constraint map}} over the entire range of redshifts where binary black hole mergers are expected to occur. A LISA observation of an equal-mass inspiral event with total redshifted mass of 10^5 solar masses for three years should be able to measure $\dot{G}/G$ at the time of merger to better than 10^(-11)/yr.
△ Less
Submitted 13 April, 2010; v1 submitted 15 December, 2009;
originally announced December 2009.
-
Comment on "Kerr Black Holes as Particle Accelerators to Arbitrarily High Energy"
Authors:
Emanuele Berti,
Vitor Cardoso,
Leonardo Gualtieri,
Frans Pretorius,
Ulrich Sperhake
Abstract:
It has been suggested that rotating black holes could serve as particle colliders with arbitrarily high center-of-mass energy. Astrophysical limitations on the maximal spin, back-reaction effects and sensitivity to the initial conditions impose severe limits on the likelihood of such collisions.
It has been suggested that rotating black holes could serve as particle colliders with arbitrarily high center-of-mass energy. Astrophysical limitations on the maximal spin, back-reaction effects and sensitivity to the initial conditions impose severe limits on the likelihood of such collisions.
△ Less
Submitted 2 December, 2009; v1 submitted 11 November, 2009;
originally announced November 2009.
-
Fundamental Theoretical Bias in Gravitational Wave Astrophysics and the Parameterized Post-Einsteinian Framework
Authors:
Nicolas Yunes,
Frans Pretorius
Abstract:
We consider the concept of fundamental bias in gravitational wave astrophysics as the assumption that general relativity is the correct theory of gravity during the entire wave-generation and propagation regime. Such an assumption is valid in the weak field, as verified by precision experiments and observations, but it need not hold in the dynamical strong-field regime where tests are lacking. F…
▽ More
We consider the concept of fundamental bias in gravitational wave astrophysics as the assumption that general relativity is the correct theory of gravity during the entire wave-generation and propagation regime. Such an assumption is valid in the weak field, as verified by precision experiments and observations, but it need not hold in the dynamical strong-field regime where tests are lacking. Fundamental bias can cause systematic errors in the detection and parameter estimation of signals, which can lead to a mischaracterization of the universe through incorrect inferences about source event rates and populations. We propose a remedy through the introduction of the parameterized post-Einsteinian framework, which consists of the enhancement of waveform templates via the inclusion of post-Einsteinian parameters. These parameters would ostensibly be designed to interpolate between templates constructed in general relativity and well-motivated alternative theories of gravity, and also include extrapolations that follow sound theoretical principles, such as consistency with conservation laws and symmetries. As an example, we construct parameterized post-Einsteinian templates for the binary coalescence of equal-mass, non-spinning compact objects in a quasi-circular inspiral. The parametrized post-Einsteinian framework should allow matched filtered data to select a specific set of post-Einsteinian parameters without a priori assuming the validity of the former, thus either verifying general relativity or pointing to possible dynamical strong-field deviations.
△ Less
Submitted 7 December, 2009; v1 submitted 17 September, 2009;
originally announced September 2009.
-
Ultra Relativistic Particle Collisions
Authors:
Matthew W. Choptuik,
Frans Pretorius
Abstract:
We present results from numerical solution of the Einstein field equations describing the head-on collision of two solitons boosted to ultra relativistic energies. We show, for the first time, that at sufficiently high energies the collision leads to black hole formation, consistent with hoop conjecture arguments. This implies that the non-linear gravitational interaction between the kinetic energ…
▽ More
We present results from numerical solution of the Einstein field equations describing the head-on collision of two solitons boosted to ultra relativistic energies. We show, for the first time, that at sufficiently high energies the collision leads to black hole formation, consistent with hoop conjecture arguments. This implies that the non-linear gravitational interaction between the kinetic energy of the solitons causes gravitational collapse, and that arguments for black hole formation in super-Planck scale particle collisions are robust.
△ Less
Submitted 11 July, 2010; v1 submitted 12 August, 2009;
originally announced August 2009.