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

Pablo Laguna

Pablo Laguna 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.HEastro-ph.COastro-phhep-phastro-ph.GAastro-ph.SRcond-mathep-thphysics.comp-phRobotics

Catalog footprint

What is connected

43works
11topics
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

43 published item(s)

preprint2026arXiv

GW200105: A detailed study of eccentricity in the neutron star-black hole binary

GW200105_162426 is the first neutron star-black hole merger to be confidently confirmed through either gravitational-wave or electromagnetic observations. Although initially analyzed after detection, the event has recently gained renewed attention following a study [Morras et al. arXiv:2503.15393] that employed a post-Newtonian inspiral-only waveform model and reported strong evidence for orbital eccentricity. In this work, we perform a detailed analysis of GW200105 using state-of-the-art effective-one-body waveform models. Importantly, we present the first study of this event utilizing a physically complete model that incorporates both orbital eccentricity and spin precession across the full inspiral, merger, and ringdown stages, along with higher-order gravitational wave modes. Our results support the presence of eccentricity in the signal, with zero eccentricity excluded from the 99% credible interval, but yielding a mass ratio closer to the original LIGO-Virgo-KAGRA analysis, differing from the findings of [Morras et al. arXiv:2503.15393]. Additionally, similar to a previous eccentric-only analysis [de Lluc Planas et al. Astrophys. J. 995, 47 (2025).], we observe a multimodal structure in the eccentricity posterior distribution. We conduct targeted investigations to understand the origin of this multimodality and complement our analysis with numerical relativity simulations to examine how the inclusion of eccentricity impacts the merger dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Calculating Quasi-Normal Modes of Schwarzschild Black Holes with Physics Informed Neural Networks

Machine learning, particularly neural networks, has rapidly permeated most activities and work where data has a story to tell. Recently, deep learning has started to be used for solving differential equations with input from physics, also known as Physics Informed Neural Networks (PINNs). We present a study showing the efficacy of PINNs for solving the Zerilli and the Regge-Wheeler equations in the time domain to calculate the quasi-normal oscillation modes of a Schwarzschild black hole. We compare the extracted modes with those obtained with finite difference methods. Although the PINN results are competitive, with a few percent differences in the quasi-normal modes estimates relative to those computed with finite difference methods, the real power of PINNs will emerge when applied to large dimensionality problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Numerical relativity for next-generation gravitational-wave probes of fundamental physics

The next generation of gravitational-wave detectors, conceived to begin operations in the 2030s, will probe fundamental physics with exquisite sensitivity. These observations will measure the equation of state of dense nuclear matter in the most extreme environments in the universe, reveal with exquisite fidelity the nonlinear dynamics of warped spacetime, put general relativity to the strictest test, and perhaps use black holes as cosmic particle detectors. Achieving each of these goals will require a new generation of numerical relativity simulations that will run at scale on the supercomputers of the 2030s to achieve the necessary accuracy, which far exceeds the capabilities of numerical relativity and high-performance computing infrastructures available today.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Field-mediated locomotor dynamics on highly deformable surfaces

In many systems motion occurs on deformed and deformable surfaces, setting up the possibility for dynamical interactions solely mediated by the coupling of the entities with their environment. Here we study the "two-body" dynamics of robot locomotion on a highly deformable spandex membrane in two scenarios: one in which a robot orbits a large central depression and the other where the two robots affect each other's motion solely through mutual environmental deformations. Inspired by the resemblance of the orbits of the single robot with those of general relativistic orbits around black holes, we recast the vehicle plus membrane dynamics in physical space into the geodesic motion of a "test particle" in a fiducial curved space-time and demonstrate how this framework facilitates understanding the observed dynamics. The two-robot problem also exhibits a resemblance with Einstein's general relativistic view of gravity, which in the words of Wheeler: "spacetime tells matter how to move; matter tells spacetime how to curve." We generalize this case the mapping to include a reciprocal coupling that translates into robotic curvature-based control schemes which modify interaction (promoting avoidance or aggregation) without long-range sensing. Our work provides a starting point for developing a mechanical analog gravity system as well as develops a framework that can provide insights into active matter in deformable environments and robot exploration in complex landscapes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inside the Final Black Hole: Puncture and Trapped Surface Dynamics

A popular approach in numerical simulations of black hole binaries is to model black holes as punctures in the fabric of spacetime. The location and the properties of the black hole punctures are tracked with apparent horizons, namely outermost marginally outer trapped surfaces (MOTSs). As the holes approach each other, a common apparent horizon suddenly appears, engulfing the two black holes and signaling the merger. The evolution of common apparent horizons and their connection with gravitational wave emission have been studied in detail with the framework of dynamical horizons. We present a study of the dynamics of the MOTSs and their punctures in the interior of the final black hole. The study focuses on head-on mergers for various initial separations and mass ratios. We find that MOTSs intersect for most of the parameter space. We show that for those situations in which they do not, it is because of the singularity avoidance property of the moving puncture gauge condition used in the study. Although we are unable to carry out evolutions that last long enough to show the ultimate fate of the punctures, our results suggest that MOTSs always intersect and that at late times their overlap is only partial. As a consequence, the punctures inside the MOTSs, although close enough to each other to act effectively as a single puncture, do not merge.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bowen-York Type Initial Data for Binaries with Neutron Stars

A new approach to construct initial data for binary systems with neutron star components is introduced. The approach is a generalization of the puncture initial data method for binary black holes based on Bowen-York solutions to the momentum constraint. As with binary black holes, the method allows setting orbital configurations with direct input from post-Newtonian approximations and involves solving only the Hamiltonian constraint. The effectiveness of the method is demonstrated with evolutions of double neutron star and black hole -- neutron star binaries in quasi-circular orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Relativistic tidal effects in non standard Kerr space-time

Astrophysical phenomena involving massive black holes (BHs) in close binaries are expected to leave detectable signatures in the electromagnetic and gravitational-wave spectrum. Such imprints may provide precious information to probe the space-time around rotating BHs and to reveal new insights on the nature of gravity in the strong-field regime. To support this observational window, it is crucial to develop suitable tests to verify the predictions of General Relativity. In this framework, the metric recently proposed by Johannsen and Psaltis parametrizes strong-field deviations from a Kerr space-time in a theory-independent way. In the following, we make use of this approach to describe the tidal field produced by spinning BHs. We compute the gravito-magnetic and gravito-electric tidal tensors for particles moving on equatorial circular geodesics, comparing our results with those obtained in the standard General Relativity scenario. Our calculations show significant differences even for distances far form the last stable orbit, which may affect the evolution of the binary and leave detectable signatures. We test our framework computing quasiequilibrium sequences of BH-white dwarf systems by means of the affine model, for different binary configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Busting Up Binaries: Encounters Between Compact Binaries and a Supermassive Black Hole

Given the stellar density near the galactic center, close encounters between compact object binaries and the supermassive black hole are a plausible occurrence. We present results from a numerical study of close to 13 million such encounters. Consistent with previous studies, we corroborate that, for binary systems tidally disrupted by the black hole, the component of the binary remaining bound to the hole has eccentricity ~ 0.97 and circularizes dramatically by the time it enters the classical LISA band. Our results also show that the population of surviving binaries merits attention. These binary systems experience perturbations to their internal orbital parameters with potentially interesting observational consequences. We investigated the regions of parameter space for survival and estimated the distribution of orbital parameters post-encounter. We found that surviving binaries harden and their eccentricity increases, thus accelerating their merger due gravitational radiation emission and increasing the predicted merger rates by up to 1%.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Testing General Relativity with Present and Future Astrophysical Observations

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ultra-Close Encounters of Stars With Massive Black Holes: Tidal Disruption Events With Prompt Hyperaccretion

A bright flare from a galactic nucleus followed at late times by a $t^{-5/3}$ decay in luminosity is often considered the signature of the complete tidal disruption of a star by a massive black hole. The flare and power-law decay are produced when the stream of bound debris returns to the black hole, self-intersects, and eventually forms an accretion disk or torus. In the canonical scenario of a solar-type star disrupted by a $10^{6}\; M_\odot$ black hole, the time between the disruption of the star and the formation of the accretion torus could be years. We present fully general relativistic simulations of a new class of tidal disruption events involving ultra-close encounters of solar-type stars with intermediate mass black holes. In these encounters, a thick disk forms promptly after disruption, on timescales of hours. After a brief initial flare, the accretion rate remains steady and highly super-Eddington for a few days at $\sim 10^2\,M_\odot\,{\rm yr}^{-1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Critical Collapse of Scalar Fields Beyond Axisymmetry

We investigate non-spherically symmetric, scalar field collapse of a family of initial data consisting of a spherically symmetric profile with a deformation proportional to the real part of the spherical harmonic $Y_{21}(θ,φ)$. Independent of the strength of the anisotropy in the data, we find that supercritical collapse yields a black hole mass scaling $M_h \propto (p-p^*)^γ$ with $γ\approx 0.37$, a value remarkably close to the critical exponent obtained by Choptuik in his pioneering study in spherical symmetry. We also find hints of discrete self-similarity. However, the collapse experiments are not sufficiently close to the critical solution to unequivocally claim that the detected periodicity is from critical collapse echoing.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electromagnetic waves from neutron stars and black holes driven by polar gravitational perturbations

Neutron stars and black holes are the most compact astrophysical objects we can think of and as a consequence they are the main sources of gravitational waves. There are many astrophysically relevant scenarios in which these objects are immersed in or endowed with strong magnetic fields, in such a way that gravitational perturbations can couple to electromagnetic ones and can potentially trigger synergistic electromagnetic signatures. In a recent paper we derived the main equations for gravito-electromagnetic perturbations and studied in detail the case of polar electromagnetic perturbations driven by axial gravitational perturbations. In this paper we deal with the case of axial electromagnetic perturbations driven by polar black-hole or neutron stars oscillations, in which the energy emitted in case is considerably larger than in the previous case. In the case of neutron stars the phenomenon lasts considerably longer since the fluid acts as an energy reservoir that shakes the magnetic field for a timescale of the order of secs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration

The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravitationally Driven Electromagnetic Perturbations of Neutron Stars and Black Holes

Gravitational perturbations of neutron stars and black holes are well known sources of gravitational radiation. If the compact object is immersed in or endowed with a magnetic field, the gravitational perturbations would couple to electromagnetic perturbations and potentially trigger synergistic electromagnetic signatures. We present a detailed analytic calculation of the dynamics of coupled gravitational and electromagnetic perturbations for both neutron stars and black holes. We discuss the prospects for detecting the electromagnetic waves in these scenarios and the potential that these waves have for providing information about their source.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

GRB060218 as a Tidal Disruption of a White Dwarf by an Intermediate Mass Black Hole

A highly unusual pair of a gamma-ray burst (GRB) GRB060218 and an associated supernova SN2006aj has puzzled theorists for years. A supernova shock breakout and a jet from a newborn stellar mass compact object were put forward to explain its multiwavelength signature. We propose that the source is naturally explained by another channel, a tidal disruption of a white dwarf (WD) by an intermediate mass black hole (IMBH). The tidal disruption is accompanied by a tidal pinching, which leads to the ignition of a WD and a supernova. Some debris falls back onto the IMBH, forms a disk, which quickly amplifies the magnetic field, and launches a jet. We successfully fit soft X-ray spectrum with the Comptonized blackbody emission from a jet photosphere. The optical/UV emission is consistent with self-absorbed synchrotron from the expanding jet front. The accretion rate temporal dependence Mdot(t) in a tidal disruption provides a good fit to soft X-ray lightcurve. The IMBH mass is found to be about 10^4Msun in three independent estimates: (1) fitting tidal disruption Mdot(t) to soft X-ray lightcurve; (2) computing the jet base radius in a jet photospheric emission model; (3) inferring the central BH mass based on a host dwarf galaxy stellar mass. The supernova position is consistent with the center of the host galaxy, while low supernova ejecta mass is consistent with a WD mass. High expected rate of tidal disruptions in dwarf galaxies is consistent with one source observed by Swift satellite over several years at GRB060218 distance of 150Mpc. The encounters with the WDs provide a lot of fuel for IMBH growth.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Template Mode Hierarchies for Binary Black Hole Mergers

Matched filtering is a popular data analysis framework used to search for gravitational wave signals emitted by compact object binaries. The templates used in matched filtering searches are constructed predominantly from the quadrupolar mode because this mode is the energetically most dominant channel. However, for highly precessing binaries or binaries with moderately large mass ratios, significant power is also carried by higher-order modes. We investigate how the inclusion of higher modes in the templates increases the prospects for detecting gravitational waves. Specifically, we use numerical relativity waveforms from the late inspiral and coalescence of binary black holes to identify mode hierarchies that cover the sky of binary orientations. We show that the ordering in these hierarchies depends on the characteristics of the binary system and the mode strengths. Our study demonstrates that detecting moderately high precessing or unequal mass binaries requires the inclusion of higher modes in the templates banks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Transient Gravitational-Wave Sky

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Impact of Higher-order Modes on the Detection of Binary Black Hole Coalescences

The inspiral and merger of black-hole binary systems are a promising source of gravitational waves. The most effective method to look for a signal with a well understood waveform, such as the binary black hole signal, is matched filtering against a library of model waveforms. Current model waveforms are comprised solely of the dominant radiation mode, the quadrupole mode, although it is known that there can be significant power in the higher-order modes. The binary black hole waveforms produced by numerical relativity are accurate through late inspiral, merger, and ringdown and include the higher-order modes. The available numerical-relativity waveforms span an increasing portion of the physical parameter space of unequal mass, spin and precession. In this paper, we investigate the degree to which gravitational-wave searches could be improved by the inclusion of higher modes in the model waveforms, for signals with a variety of initial mass ratios and generic spins. Our investigation studies how well the quadrupole-only waveform model matches the signal as a function of the inclination and orientation of the source and how the modes contribute to the distance reach into the Universe of Advanced LIGO for a fixed set of internal source parameters. The mismatch between signals and quadrupole-only waveform can be large, dropping below 0.97 for up to 65% of the source-sky for the non-precessing cases we studied, and over a larger area in one precessing case. There is a corresponding 30% increase in detection volume that could be achieved by adding higher modes to the search; however, this is mitigated by the fact that the mismatch is largest for signals which radiate the least energy and to which the search is therefore least sensitive. Likewise, the mismatch is largest in directions from the source along which the least energy is radiated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Prompt Emission from Tidal Disruptions of White Dwarfs by Intermediate Mass Black Holes

We present a qualitative picture of prompt emission from tidal disruptions of white dwarfs (WD) by intermediate mass black holes (IMBH). The smaller size of an IMBH compared to a supermassive black hole and a smaller tidal radius of a WD disruption lead to a very fast event with high peak luminosity. Magnetic field is generated in situ following the tidal disruption, which leads to effective accretion. Since large-scale magnetic field is also produced, geometrically thick super-Eddington inflow leads to a relativistic jet. The dense jet possesses a photosphere, which emits quasi-thermal radiation in soft X-rays. The source can be classified as a long low-luminosity gamma-ray burst (ll-GRB). Tidal compression of a WD causes nuclear ignition, which is observable as an accompanying supernova. We suggest that GRB060218 and SN2006aj is such a pair of ll-GRB and supernova. We argue that in a flux-limited sample the disruptions of WDs by IMBHs are more frequent then the disruptions of other stars by IMBHs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries

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 and parameter-estimation algorithms using numerically generated waveforms, and to foster closer collaboration between the numerical relativity and data analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian - numerical relativity hybrid waveforms, large numbers of injections, and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical relativity groups have contributed 63 hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the $(\ell,m)=(2,2)$ mode. Based on these studies we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tidal Disruptions of White Dwarfs from Ultra-Close Encounters with Intermediate Mass Spinning Black Holes

We present numerical relativity results of tidal disruptions of white dwarfs from ultra-close encounters with a spinning, intermediate mass black hole. These encounters require a full general relativistic treatment of gravity. We show that the disruption process and prompt accretion of the debris strongly depend on the magnitude and orientation of the black hole spin. However, the late-time accretion onto the black hole follows the same decay, $\dot{M}$ ~ t^{-5/3}, estimated from Newtonian gravity disruption studies. We compute the spectrum of the disk formed from the fallback material using a slim disk model. The disk spectrum peaks in the soft X-rays and sustains Eddington luminosity for 1-3 yrs after the disruption. For arbitrary black hole spin orientations, the disrupted material is scattered away from the orbital plane by relativistic frame dragging, which often leads to obscuration of the inner fallback disk by the outflowing debris. The disruption events also yield bursts of gravitational radiation with characteristic frequencies of ~3.2 Hz and strain amplitudes of ~10^{-18} for galactic intermediate mass black holes. The optimistic rate of considered ultra-close disruptions is consistent with no sources found in ROSAT all-sky survey. The future missions like Wide-Field X-ray Telescope (WFXT) could observe dozens of events.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Late Inspiral and Merger of Binary Black Holes in Scalar-Tensor Theories of Gravity

Gravitational wave observations will probe non-linear gravitational interactions and thus enable strong tests of Einstein's theory of general relativity. We present a numerical relativity study of the late inspiral and merger of binary black holes in scalar-tensor theories of gravity. We consider black hole binaries in an inhomogeneous scalar field, specifically binaries inside a scalar field bubble, in some cases with a potential. We calculate the emission of dipole radiation. We also show how these configurations trigger detectable differences between gravitational waves in scalar-tensor gravity and the corresponding waves in general relativity. We conclude that, barring an external mechanism to induce dynamics in the scalar field, scalar-tensor gravity binary black holes alone are not capable of awaking a dormant scalar field, and are thus observationally indistinguishable from their general relativistic counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mergers of Supermassive Black Holes in Astrophysical Environments

Modeling the late inspiral and merger of supermassive black holes is central to understanding accretion processes and the conditions under which electromagnetic emission accompanies gravitational waves. We use fully general relativistic, hydrodynamics simulations to investigate how electromagnetic signatures correlate with black hole spins, mass ratios, and the gaseous environment in this final phase of binary evolution. In all scenarios, we find some form of characteristic electromagnetic variability whose pattern depends on the spins and binary mass ratios. Binaries in hot accretion flows exhibit a flare followed by a sudden drop in luminosity associated with the plunge and merger, as well as quasi-periodic oscillations correlated with the gravitational waves during the inspiral. Conversely, circumbinary disk systems are characterized by a low luminosity of variable emission, suggesting challenging prospects for their detection.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Properties of Accretion Flows Around Coalescing Supermassive Black Holes

What are the properties of accretion flows in the vicinity of coalescing supermassive black holes (SBHs)? The answer to this question has direct implications for the feasibility of coincident detections of electromagnetic (EM) and gravitational wave (GW) signals from coalescences. Such detections are considered to be the next observational grand challenge that will enable testing general relativity in the strong, nonlinear regime and improve our understanding of evolution and growth of these massive compact objects. In this paper we review the properties of the environment of coalescing binaries in the context of the circumbinary disk and hot, radiatively inefficient accretion flow models and use them to mark the extent of the parameter space spanned by this problem. We report the results from an ongoing, general relativistic, hydrodynamical study of the inspiral and merger of black holes, motivated by the latter scenario. We find that correlated EM+GW oscillations can arise during the inspiral phase followed by the gradual rise and subsequent drop-off in the light curve at the time of coalescence. While there are indications that the latter EM signature is a more robust one, a detection of either signal coincidentally with GWs would be a convincing evidence for an impending SBH binary coalescence. The observability of an EM counterpart in the hot accretion flow scenario depends on the details of a model. In the case of the most massive binaries observable by the Laser Interferometer Space Antenna, upper limits on luminosity imply that they may be identified by EM searches out to z~0.1-1. However, given the radiatively inefficient nature of the gas flow, we speculate that a majority of massive binaries may appear as low luminosity AGN in the local universe.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Super-Extremal Spinning Black Holes via Accretion

A Kerr black hole with mass $M$ and angular momentum $J$ satisfies the extremality inequality $|J| \le M^2$. In the presence of matter and/or gravitational radiation, this bound needs to be reformulated in terms of local measurements of the mass and the angular momentum directly associated with the black hole. The isolated and dynamical horizon framework provides such quasi-local characterization of black hole mass and angular momentum. With this framework, it is possible in axisymmetry to reformulate the extremality limit as $|J| \le 2\,M_H^2$, with $M_H$ the irreducible mass of the black hole computed from its apparent horizon area and $J$ obtained using approximate rotational Killing vectors on the apparent horizon. The $|J| \le 2\,M_H^2$ condition is also equivalent to requiring a non-negative black hole surface gravity. We present numerical experiments of an accreting black hole that temporarily violates this extremality inequality. The initial configuration consists of a single, rotating black hole surrounded by a thick, shell cloud of negative energy density. For these numerical experiments, we introduce a new matter-without-matter evolution method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Einstein Toolkit: A Community Computational Infrastructure for Relativistic Astrophysics

We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, compact objects, and collapsing stars, as well as a full suite of analysis tools. The Einstein Toolkit is currently based on the Cactus Framework for high-performance computing and the Carpet adaptive mesh refinement driver. It implements spacetime evolution via the BSSN evolution system and general-relativistic hydrodynamics in a finite-volume discretization. The toolkit is under continuous development and contains many new code components that have been publicly released for the first time and are described in this article. We discuss the motivation behind the release of the toolkit, the philosophy underlying its development, and the goals of the project. A summary of the implemented numerical techniques is included, as are results of numerical test covering a variety of sample astrophysical problems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Comparisons of eccentric binary black hole simulations with post-Newtonian models

We present the first comparison between numerical relativity (NR) simulations of an eccentric binary black hole system with corresponding post-Newtonian (PN) results. We evolve an equal-mass, non-spinning configuration with an initial eccentricity e ~ 0.1 for 21 gravitational wave cycles before merger, and find agreement in the gravitational wave phase with an adiabatic eccentric PN model with 2 PN radiation reaction within 0.1 radians for 10 cycles. The NR and PN phase difference grows to 0.7 radians by 5 cycles before merger. We find that these results can be obtained by expanding the eccentric PN expressions in terms of the frequency-related variable x = (omega M)^{2/3} with M the total mass of the binary. When using instead the mean motion n = 2 π/P, where P is the orbital period, the comparison leads to significant disagreements with NR.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Final Mass and Spin of Merged Black Holes and the Golden Black Hole

We present results on the mass and spin of the final black hole from mergers of equal mass, spinning black holes. The study extends over a broad range of initial orbital configurations, from direct plunges to quasi-circular inspirals to more energetic orbits (generalizations of Newtonian elliptical orbits). It provides a comprehensive search of those configurations that maximize the final spin of the remnant black hole. We estimate that the final spin can reach a maximum spin $a/M_h \approx 0.99\pm 0.01$ for extremal black hole mergers. In addition, we find that, as one increases the orbital angular momentum from small values, the mergers produce black holes with mass and spin parameters $\lbrace M_h/M, a/M_h \rbrace$ ~spiraling around the values $\lbrace \hat M_h/M, \hat a/M_h \rbrace$ of a {\it golden} black hole. Specifically, $(M_h-\hat M_h)/M \propto e^{\pm B\,ϕ}\cosϕ$ and $(a-\hat a)/M_h \propto e^{\pm C\,ϕ}\sinϕ$, with $ϕ$ a monotonically growing function of the initial orbital angular momentum. We find that the values of the parameters for the \emph{golden} black hole are those of the final black hole obtained from the merger of a binary with the corresponding spinning black holes in a quasi-circular inspiral.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Integrated Sachs-Wolfe Effect for Gravitational Radiation

Gravitational waves are messengers carrying valuable information about their sources. For sources at cosmological distances, the waves will contain also the imprint left by the intervening matter. The situation is in close analogy with cosmic microwave photons, for which the large-scale structures the photons traverse contribute to the observed temperature anisotropies, in a process known as the integrated Sachs-Wolfe effect. We derive the gravitational wave counterpart of this effect for waves propagating on a Friedman-Robertson-Walker background with scalar perturbations. We find that the phase, frequency and amplitude of the gravitational waves experience Sachs-Wolfe type integrated effects, this in addition to the magnification effects on the amplitude from gravitational lensing. We show that for supermassive black hole binaries, the integrated effects could account for measurable changes on the frequency, chirp mass and luminosity distance of the binary, thus unveiling the presence of inhomogeneities, and potentially dark energy, in the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Relativistic Mergers of Supermassive Black Holes and their Electromagnetic Signatures

Coincident detections of electromagnetic (EM) and gravitational wave (GW) signatures from coalescence events of supermassive black holes are the next observational grand challenge. Such detections will provide the means to study cosmological evolution and accretion processes associated with these gargantuan compact objects. More generally, the observations will enable testing general relativity in the strong, nonlinear regime and will provide independent cosmological measurements to high precision. Understanding the conditions under which coincidences of EM and GW signatures arise during supermassive black hole mergers is therefore of paramount importance. As an essential step towards this goal, we present results from the first fully general relativistic, hydrodynamical study of the late inspiral and merger of equal-mass, spinning supermassive black hole binaries in a gas cloud. We find that variable EM signatures correlated with GWs can arise in merging systems as a consequence of shocks and accretion combined with the effect of relativistic beaming. The most striking EM variability is observed for systems where spins are aligned with the orbital axis and where orbiting black holes form a stable set of density wakes, but all systems exhibit some characteristic signatures that can be utilized in searches for EM counterparts. In the case of the most massive binaries observable by the Laser Interferometer Space Antenna, calculated luminosities imply that they may be identified by EM searches to z = 1, while lower mass systems and binaries immersed in low density ambient gas can only be detected in the local universe.

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
preprint2009arXiv

The Samurai Project: verifying the consistency of black-hole-binary waveforms for gravitational-wave detection

We quantify the consistency of numerical-relativity black-hole-binary waveforms for use in gravitational-wave (GW) searches with current and planned ground-based detectors. We compare previously published results for the $(\ell=2,| m | =2)$ mode of the gravitational waves from an equal-mass nonspinning binary, calculated by five numerical codes. We focus on the 1000M (about six orbits, or 12 GW cycles) before the peak of the GW amplitude and the subsequent ringdown. We find that the phase and amplitude agree within each code's uncertainty estimates. The mismatch between the $(\ell=2,| m| =2)$ modes is better than $10^{-3}$ for binary masses above $60 M_{\odot}$ with respect to the Enhanced LIGO detector noise curve, and for masses above $180 M_{\odot}$ with respect to Advanced LIGO, Virgo and Advanced Virgo. Between the waveforms with the best agreement, the mismatch is below $2 \times 10^{-4}$. We find that the waveforms would be indistinguishable in all ground-based detectors (and for the masses we consider) if detected with a signal-to-noise ratio of less than $\approx14$, or less than $\approx25$ in the best cases.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Gravitational Recoil from Binary Black Hole Mergers: the Close-Limit Approximation

[abridged] The coalescence of a binary black hole system is one of the main sources of gravitational waves that present and future detectors will study. Apart from the energy and angular momentum that these waves carry, for unequal-mass binaries there is also a net flux of linear momentum that implies a recoil velocity of the resulting final black hole in the opposite direction. We present a computation of the recoil velocity based on the close-limit approximation scheme, which gives excellent results for head-on and grazing collisions of black holes when compared to full numerical relativistic calculations. We obtain a maximum recoil velocity of ~ 64 km/s for a symmetric mass ratio eta = M_1 M_2/(M_1+M_2)^2 ~ 0.19 and an initial proper separation of 4 M, where M is the total ADM mass of the system. This separation is the maximum at which the close-limit approximation is expected to provide accurate results. If we supplement this estimate with PN calculations up to the innermost stable circular orbit, we obtain a lower bound for the recoil velocity, with a maximum around 84 km/s. This is a lower bound because it neglects the initial merger phase. We can however obtain a rough estimate by using PN methods or the close-limit approximation. Since both methods are known to overestimate the amount of radiation, we obtain in this way an upper bound for the recoil with maxima in the range of 220-265 km/s. We also provide non-linear fits to these estimated upper and lower bounds. These estimates are subject to uncertainties related to issues such as the choice of initial data and higher effects in perturbation theory. Nonetheless, our estimates are consistent with previous results in the literature and suggest a narrower range of possible recoil velocities.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Gravitational recoil from spinning binary black hole mergers

The inspiral and merger of binary black holes will likely involve black holes with both unequal masses and arbitrary spins. The gravitational radiation emitted by these binaries will carry angular as well as linear momentum. A net flux of emitted linear momentum implies that the black hole produced by the merger will experience a recoil or kick. Previous studies have focused on the recoil velocity from unequal mass, non-spinning binaries. We present results from simulations of equal mass but spinning black hole binaries and show how a significant gravitational recoil can also be obtained in these situations. We consider the case of black holes with opposite spins of magnitude $a$ aligned/anti-aligned with the orbital angular momentum, with $a$ the dimensionless spin parameters of the individual holes. For the initial setups under consideration, we find a recoil velocity of $V = 475 \KMS a$. Supermassive black hole mergers producing kicks of this magnitude could result in the ejection from the cores of dwarf galaxies of the final hole produced by the collision.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Constant Crunch Coordinates for Black Hole Simulations

We reinvestigate the utility of time-independent constant mean curvature foliations for the numerical simulation of a single spherically-symmetric black hole. Each spacelike hypersurface of such a foliation is endowed with the same constant value of the trace of the extrinsic curvature tensor, $K$. Of the three families of $K$-constant surfaces possible (classified according to their asymptotic behaviors), we single out a sub-family of singularity-avoiding surfaces that may be particularly useful, and provide an analytic expression for the closest approach such surfaces make to the singularity. We then utilize a non-zero shift to yield families of $K$-constant surfaces which (1) avoid the black hole singularity, and thus the need to excise the singularity, (2) are asymptotically null, aiding in gravity wave extraction, (3) cover the physically relevant part of the spacetime, (4) are well behaved (regular) across the horizon, and (5) are static under evolution, and therefore have no ``grid stretching/sucking'' pathologies. Preliminary numerical runs demonstrate that we can stably evolve a single spherically-symmetric static black hole using this foliation. We wish to emphasize that this coordinatization produces $K$-constant surfaces for a single black hole spacetime that are regular, static and stable throughout their evolution.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Construction of Astrophysical Initial Data for Perturbations of Relativistic Stars

We develop a framework for constructing initial data sets for perturbations about spherically symmetric matter distributions. This framework facilitates setting initial data representing astrophysical sources of gravitational radiation involving relativistic stars. The procedure is based on Lichnerowicz-York's conformal approach to solve the constraints in Einstein's equations. The correspondence of these initial data sets in terms of the standard gauge perturbation variables in the Regge-Wheeler perturbation variables is established, and examples of initial data sets of merging neutron stars under the close-limit approximation are presented.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Dynamics of scalar fields in the background of rotating black holes II: A note on superradiance

We analyze the amplification due to so-called superradiance from the scattering of pulses off rotating black holes as a numerical time evolution problem. We consider the "worst possible case" of scalar field pulses for which superradiance effects yield amplifications $< 1%$. We show that this small effect can be isolated by numerically evolving quasi-monochromatic, modulated pulses with a recently developed Teukolsky code. The results show that it is possible to study superradiance in the time domain, but only if the initial data is carefully tuned. This illustrates the intrinsic difficulties of detecting superradiance in more general evolution scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Critical Dynamics of Symmetry Breaking: Quenches, Dissipation and Cosmology

Symmetry-breaking phase transitions may leave behind topological defects \cite{Kibble} with a density dependent on the quench rate \cite{Zurek}. We investigate the dynamics of such quenches for the one-dimensional, Landau-Ginzburg case and show that the density of kinks, $n$, scales differently with the quench timescale, $τ_Q$, depending on whether the dynamics in the vicinity of the critical point is overdamped ($n \propto τ_Q^{-1/4}$) or underdamped ($n \propto τ_Q^{-1/3}$). Either of these cases may be relevant to the early Universe, and we derive bounds on the initial density of topological defects in cosmological phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Density of kinks after a quench: When symmetry breaks, how big are the pieces?

Numerical study of order parameter dynamics in the course of second order (Landau-Ginzburg) symmetry breaking transitions shows that the density of topological defects, kinks, is proportional to the fourth root of the rate of the quench. This confirms the more general theory of domain-size evolution in the course of symmetry breaking transformations proposed by one of us \cite{Zurek85}. Using these ideas, it is possible to compute the density of topological defects from the quench timescale and from the equilibrium scaling of the correlation length and relaxation time near the critical point.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

The Imprint of Proper Motion of Nonlinear Structures on the Cosmic Microwave Background

We investigate the imprint of nonlinear matter condensations on the Cosmic Microwave Background (CMB) in an $Ω=1$, Cold Dark Matter (CDM) model universe. Temperature anisotropies are obtained by numerically evolving matter inhomogeneities and CMB photons from the beginning of decoupling until the present epoch. The underlying density field produced by the inhomogeneities is followed from the linear, through the weakly clustered, into the fully nonlinear regime. We concentrate on CMB temperature distortions arising from variations in the gravitational potentials of nonlinear structures. We find two sources of temperature fluctuations produced by time-varying potentials: (1) anisotropies due to intrinsic changes in the gravitational potentials of the inhomogeneities and (2) anisotropies generated by the peculiar, bulk motion of the structures across the microwave sky. Both effects generate CMB anisotropies in the range of $10^{-7} \siml ΔT/T \siml 10^{-6}$ on scales of $ \sim 1^{\circ}$. For isolated structures, anisotropies due to proper motion exhibit a dipole-like signature in the CMB sky that in principle could yield information on the transverse velocity of the structures.

0 citations0 reviewsNo code linkedOpen access