BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Frans Pretorius

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEhep-thastro-ph.COhep-phastro-ph.GAastro-ph.SRmath-phmath.MPnucl-thastro-ph.EPphysics.comp-ph

Catalog footprint

What is connected

49works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

49 published item(s)

preprint2024arXiv

Nonlinear Effects In Black Hole Ringdown From Scattering Experiments I: spin and initial data dependence of quadratic mode coupling

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.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Causal, stable first-order viscous relativistic hydrodynamics with ideal gas microphysics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Conservative finite volume scheme for first-order viscous relativistic hydrodynamics

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).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamics of a Z2 symmetric EdGB gravity in spherical symmetry

We report on a numerical investigation of black hole evolution in an Einstein dilaton Gauss-Bonnet (EdGB) gravity theory where the Gauss-Bonnet coupling and scalar (dilaton) field potential are symmetric under a global change in sign of the scalar field (a 'Z2' symmetry). We find that for sufficiently small Gauss-Bonnet couplings Schwarzschild black holes are stable to radial scalar field perturbations, and are unstable to such perturbations for sufficiently large couplings. For the latter case, we provide numerical evidence that there is a band of coupling parameters and black hole masses where the end states are stable scalarized black hole solutions, in general agreement with the results of Macedo et al (2019 Phys. Rev. D 99 104041). For Gauss-Bonnet couplings larger than those in the stable band, we find that an 'elliptic region' forms outside of the black hole horizon, indicating the theory does not possess a well-posed initial value formulation in that regime.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Robustness of slow contraction to cosmic initial conditions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Scalarized Black Hole dynamics in Einstein dilaton Gauss-Bonnet Gravity

We report on a numerical investigation of the stability of scalarized black holes in Einstein dilaton Gauss-Bonnet (EdGB) gravity in the full dynamical theory, though restricted to spherical symmetry. We find evidence that for sufficiently small curvature-couplings the resulting scalarized black hole solutions are nonlinearly stable. For such small couplings, we show that an elliptic region forms inside these EdGB black hole spacetimes (prior to any curvature singularity), and give evidence that this region remains censored from asymptotic view. However, for coupling values "superextremal" relative to a given black hole mass, an elliptic region forms exterior to the horizon, implying the exterior Cauchy problem is ill-posed in this regime.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The rotating black hole interior: Insights from gravitational collapse in $AdS_3$ spacetime

We present results from a numerical study of rotating black hole formation in 3-dimensional asymptotically anti-de Sitter (AdS) spacetime, focusing on the structure of the black hole interior. While black holes in $AdS_3$ are of theoretical interest for a wide variety of reasons, we choose to study this system primarily as a toy model for astrophysical (4-dimensional) black holes formed from gravitational collapse. We investigate the effect of angular momentum on the geometry inside the event horizon, and see qualitative changes in the interior structure as a function of the spin parameter. For low spins, we find that a central spacelike curvature singularity forms, connecting to a singular, null Cauchy horizon. For spins above a threshold consistent with the linear analysis of Dias, Reall and Santos, curvature on the Cauchy horizon remains bounded, signaling a violation of the strong cosmic censorship conjecture. Further increasing the spin leads to a decrease in the relative size of the spacelike branch of the singularity, which vanishes completely above a second threshold. In these high-spin cases, the interior evolution is bounded by a regular Cauchy horizon, which extends all the way inward to a regular, timelike origin. We further explore the geodesic focusing ("gravitational shock-wave") effect predicted to occur along the outgoing branch of the inner horizon, first described by Marolf and Ori. Remarkably, we observe the effect at late times in all of the black holes we form, even those in which the inner apparent horizon collapses to zero radius early in their evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Equation of state effects and one-arm spiral instability in hypermassive neutron stars formed in eccentric neutron star mergers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Relativistic Simulations of Eccentric Binary Neutron Star Mergers: One-arm Spiral Instability and Effects of Neutron Star Spin

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spontaneous Scalarization with Massive Fields

We study the effect of a mass term in the spontaneous scalarization of neutron stars, for a wide range of scalar field parameters and neutron star equations of state. Even though massless scalars have been the focus of interest in spontaneous scalarization so far, recent observations of binary systems rule out most of their interesting parameter space. We point out that adding a mass term to the scalar field potential is a natural extension to the model that avoids these observational bounds if the Compton wavelength of the scalar is small compared to the binary separation. Our model is formally similar to the asymmetron scenario recently introduced in application to cosmology, though here we are interested in consequences for neutron stars and thus consider a mass term that does not modify the geometry on cosmological scales. We review the allowed values for the mass and scalarization parameters in the theory given current binary system observations and black hole spin measurements. We show that within the allowed ranges, spontaneous scalarization can have nonperturbative, strong effects that may lead to observable signatures in binary neutron star or black hole-neutron star mergers, or even in isolated neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Theoretical Physics Implications of the Binary Black-Hole Mergers GW150914 and GW151226

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Eccentric mergers of black holes with spinning neutron stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gravity-dominated unequal-mass black hole collisions

We continue our series of studies of high-energy collisions of black holes investigating unequal-mass, boosted head-on collisions in four dimensions. We show that the fraction of the center-of-mass energy radiated as gravitational waves becomes independent of mass ratio and approximately equal to $13\%$ at large energies. We support this conclusion with calculations using black hole perturbation theory and Smarr's zero-frequency limit approximation. These results lend strong support to the conjecture that the detailed structure of the colliding objects is irrelevant at high energies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

One-arm Spiral Instability in Hypermassive Neutron Stars Formed by Dynamical-Capture Binary Neutron Star Mergers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Black Hole Superradiance in Dynamical Spacetime

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Detecting gravitational waves from highly eccentric compact binaries

In dense stellar regions, highly eccentric binaries of black holes and neutron stars can form through various n-body interactions. Such a binary could emit a significant fraction of its binding energy in a sequence of largely isolated gravitational wave bursts prior to merger. Given expected black hole and neutron star masses, many such systems will emit these repeated bursts at frequencies within the sensitive band of contemporary ground-based gravitational wave detectors. Unfortunately, existing gravitational wave searches are ill-suited to detect these signals. In this work, we adapt a "power stacking" method to the detection of gravitational wave signals from highly eccentric binaries. We implement this method as an extension of the Q-transform, a projection onto a multiresolution basis of windowed complex exponentials that has previously been used to analyze data from the network of LIGO/Virgo detectors. Our method searches for excess power over an ensemble of time-frequency tiles. We characterize the performance of our method using Monte Carlo experiments with signals injected in simulated detector noise. Our results indicate that the power stacking method achieves substantially better sensitivity to eccentric binary signals than existing localized burst searches.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical Relativity and Astrophysics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Eccentric binary neutron star mergers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Linear Stability Analysis and the Speed of Gravitational Waves in Dynamical Chern-Simons Modified Gravity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonperturbative analysis of the evolution of cosmological perturbations through a nonsingular bounce

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Observing complete gravitational wave signals from dynamical capture binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulating extreme-mass-ratio systems in full general relativity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ultrarelativistic black hole formation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Universality, maximum radiation and absorption in high-energy collisions of black holes with spin

We explore the impact of black hole spins on the dynamics of high-energy black hole collisions. We report results from numerical simulations with gamma-factors up to 2.49 and dimensionless spin parameter +0.85, +0.6, 0, -0.6, -0.85. We find that the scattering threshold becomes independent of spin at large center-of-mass energies, confirming previous conjectures that structure does not matter in ultrarelativistic collisions. It has further been argued that in this limit all of the kinetic energy of the system may be radiated by fine tuning the impact parameter to threshold. On the contrary, we find that only about 60% of the kinetic energy is radiated for gamma=2.49. By monitoring apparent horizons before and after scattering events we show that the "missing energy" is absorbed by the individual black holes in the encounter, and moreover the individual black-hole spins change significantly. We support this conclusion with perturbative calculations. An extrapolation of our results to the limit gamma->infinity suggests that about half of the center-of-mass energy of the system can be emitted in gravitational radiation, while the rest must be converted into rest-mass and spin energy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Conformal Thin-Sandwich Solver for Generic Initial Data

We present a new scheme for constructing initial data for the Einstein field equations using the conformal thin-sandwich formulation that does not assume conformal flatness or approximate Killing vectors. This includes a method for determining free data based on superposition, as well as a way to handle black hole singularities without excision. We numerically solve the constraint equations using a multigrid algorithm with mesh refinement. We demonstrate the efficacy of the method with initial data solutions for several applications: a quasicircular binary black hole merger, a dynamical capture black hole-neutron star merger, and an ultrarelativistic collision.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Detectable seismic consequences of the interaction of a primordial black hole with Earth

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamical Capture Binary Neutron Star Mergers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Eccentric black hole-neutron star mergers: effects of black hole spin and equation of state

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hydrodynamics in full general relativity with conservative AMR

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

NR/HEP: roadmap for the future

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulation of Asymptotically AdS5 Spacetimes with a Generalized Harmonic Evolution Scheme

Motivated by the gauge/gravity duality, we introduce a numerical scheme based on generalized harmonic evolution to solve the Einstein field equations on asymptotically anti-de Sitter (AdS) spacetimes. We work in global AdS5, which can be described by the (t,r,χ,θ,ϕ) spherical coordinates adapted to the R{\times}S3 boundary. We focus on solutions that preserve an SO(3) symmetry that acts to rotate the 2-spheres parametrized by θ,ϕ. In the boundary conformal field theory (CFT), the way in which this symmetry manifests itself hinges on the way we choose to embed Minkowski space in R{\times}S3. We present results from an ongoing study of prompt black hole formation via scalar field collapse, and explore the subsequent quasi-normal ringdown. Beginning with initial data characterized by highly distorted apparent horizon geometries, the metrics quickly evolve, via quasi-normal ringdown, to equilibrium static black hole solutions at late times. The lowest angular number quasi-normal modes are consistent with the linear modes previously found in perturbative studies, whereas the higher angular modes are a combination of linear modes and of harmonics arising from non-linear mode-coupling. We extract the stress energy tensor of the dual CFT on the boundary, and find that despite being highly inhomogeneous initially, it nevertheless evolves from the outset in a manner that is consistent with a thermalized N=4 SYM fluid. As a first step towards closer contact with relativistic heavy ion collision physics, we map this solution to a Minkowski piece of the R{\times}S3 boundary, and obtain a corresponding fluid flow in Minkowski space.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The imminent detection of gravitational waves from massive black-hole binaries with pulsar timing arrays

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Dynamical Systems Approach to Schwarzschild Null Geodesics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Eccentric Black Hole-Neutron Star Mergers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Evaporation of 2-Dimensional Black Holes

We present a detailed analysis of results from a new study of the quantum evaporation of Callan-Giddings-Harvey-Strominger (CGHS) black holes within the mean-field approximation. This semi-classical theory incorporates back reaction. Our analytical and numerical calculations show that, while some of the assumptions underlying the standard evaporation paradigm are borne out, several are not. One of the anticipated properties we confirm is that the semi-classical space-time is asymptotically flat at right future null infinity, $\spr$, yet incomplete in the sense that null observers reach a future Cauchy horizon in finite affine time. Unexpected behavior includes that the Bondi mass traditionally used in the literature can become negative even when the area of the horizon is macroscopic; an improved Bondi mass remains positive until the end of semi-classical evaporation, yet the final value can be arbitrarily large relative to the Planck mass; and the flux of the quantum radiation at $\spr$ is non-thermal even when the horizon area is large compared to the Planck scale. Furthermore, if the black hole is initially macroscopic, the evaporation process exhibits remarkable universal properties. Although the literature on CGHS black holes is quite rich, these features had escaped previous analyses, in part because of lack of required numerical precision, and in part due to misinterpretation of certain properties and symmetries of the model. Finally, our results provide support for the full quantum scenario recently developed by Ashtekar, Taveras and Varadarajan, and also offer a number of interesting problems to the mathematical relativity and geometric analysis communities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Final State of Gregory-Laflamme Instability

We describe the behavior of a perturbed 5-dimensional black string subject to the Gregory-Laflamme instability. We show that the horizon evolves in a self-similar manner, where at any moment in the late-time development of the instability the horizon can be described as a sequence of 3-dimensional spherical black holes of varying size, joined by black string segments of similar radius. As with the initial black string, each local string segment is itself unstable, and this fuels the self-similar cascade to (classically) arbitrarily small scales; in the process the horizon develops a fractal structure. In finite asymptotic time, the remaining string segments shrink to zero-size, yielding a naked singularity. Since no fine-tuning is required to excite the instability, this constitutes a generic violation of cosmic censorship. We further discuss how this behavior is related to satellite formation in low-viscosity fluid streams subject to the Rayleigh-Plateau instability, and estimate the fractal dimension of the horizon prior to formation of the naked singularity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational Wave Tests of General Relativity with the Parameterized Post-Einsteinian Framework

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Surprises in the Evaporation of 2-Dimensional Black Holes

Quantum evaporation of Callan-Giddings-Harvey-Strominger (CGHS) black holes is analyzed in the mean field approximation. This semi-classical theory incorporates back reaction. Detailed analytical and numerical calculations show that, while some of the assumptions underlying the standard evaporation paradigm are borne out, several are not. Furthermore, if the black hole is initially macroscopic, the evaporation process exhibits remarkable universal properties (which are distinct from the features observed in the simplified, exactly soluble models). Although the literature on CGHS black holes is quite rich, these features had escaped previous analyses, in part because of lack of required numerical precision, and in part because certain properties and symmetries of the model were not fully recognized. Finally, our results provide support for the full quantum gravity scenario recently developed by Ashtekar, Taveras and Varadarajan.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Black Strings, Low Viscosity Fluids, and Violation of Cosmic Censorship

We describe the behavior of 5-dimensional black strings, subject to the Gregory-Laflamme instability. Beyond the linear level, the evolving strings exhibit a rich dynamics, where at intermediate stages the horizon can be described as a sequence of 3-dimensional spherical black holes joined by black string segments. These segments are themselves subject to a Gregory-Laflamme instability, resulting in a self-similar cascade, where ever-smaller satellite black holes form connected by ever-thinner string segments. This behavior is akin to satellite formation in low-viscosity fluid streams subject to the Rayleigh-Plateau instability. The simulation results imply that the string segments will reach zero radius in finite asymptotic time, whence the classical space-time terminates in a naked singularity. Since no fine-tuning is required to excite the instability, this constitutes a generic violation of cosmic censorship.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraining the evolutionary history of Newton's constant with gravitational wave observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Semianalytical estimates of scattering thresholds and gravitational radiation in ultrarelativistic black hole encounters

Ultrarelativistic collisions of black holes are ideal gedanken experiments to study the nonlinearities of general relativity. In this paper we use semianalytical tools to better understand the nature of these collisions and the emitted gravitational radiation. We explain many features of the energy spectra extracted from numerical relativity simulations using two complementary semianalytical calculations. In the first calculation we estimate the radiation by a "zero-frequency limit" analysis of the collision of two point particles with finite impact parameter. In the second calculation we replace one of the black holes by a point particle plunging with arbitrary energy and impact parameter into a Schwarzschild black hole, and we explore the multipolar structure of the radiation paying particular attention to the near-critical regime. We also use a geodesic analogy to provide qualitative estimates of the dependence of the scattering threshold on the black hole spin and on the dimensionality of the spacetime.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Superkicks in ultrarelativistic encounters of spinning black holes

We study ultrarelativistic encounters of two spinning, equal-mass black holes through simulations in full numerical relativity. Two initial data sequences are studied in detail: one that leads to scattering and one that leads to a grazing collision and merger. In all cases, the initial black hole spins lie in the orbital plane, a configuration that leads to the so-called "superkicks". In astrophysical, quasicircular inspirals, such kicks can be as large as ~3,000 km/s; here, we find configurations that exceed ~15,000 km/s. We find that the maximum recoil is to a good approximation proportional to the total amount of energy radiated in gravitational waves, but largely independent of whether a merger occurs or not. This shows that the mechanism predominantly responsible for the superkick is not related to merger dynamics. Rather, a consistent explanation is that the "bobbing" motion of the orbit causes an asymmetric beaming of the radiation produced by the in-plane orbital motion of the binary, and the net asymmetry is balanced by a recoil. We use our results to formulate some conjectures on the ultimate kick achievable in any black hole encounter.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Two-dimensional quantum black holes: Numerical methods

We present details of a new numerical code designed to study the formation and evaporation of 2-dimensional black holes within the CGHS model. We explain several elements of the scheme that are crucial to resolve the late-time behavior of the spacetime, including regularization of the field variables, compactification of the coordinates, the algebraic form of the discretized equations of motion, and the use of a modified Richardson extrapolation scheme to achieve high-order convergence. Physical interpretation of our results will be discussed in detail elsewhere.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ultra Relativistic Particle Collisions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Status of NINJA: the Numerical INJection Analysis project

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.

0 citations0 reviewsNo code linkedOpen access