BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Bela Szilagyi

Bela Szilagyi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
gr-qcastro-ph.HEastro-phastro-ph.COphysics.geo-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 published item(s)

preprint2020arXiv

Constraining the parameters of GW150914 & GW170104 with numerical relativity surrogates

Gravitational-wave detectors have begun to observe coalescences of heavy black holes at a consistent pace for the past few years. Accurate models of gravitational waveforms are essential for unbiased and precise estimation of source parameters. Recently developed surrogate models based on high-accuracy numerical relativity (NR) simulations are ideal for constraining physical parameters of heavy black hole merger events. In this paper, we first demonstrate the viability of these multi-modal surrogates as reliable parameter estimation tools. We show that NR surrogates can extract additional information from GW data that is inaccessible to traditional models, by analyzing a set of synthetic signals with the NR surrogate and other approximants. We also consider the case of two of the earliest binary black holes detected by the LIGO observatories: GW150914 and GW170104. We reanalyze their data with fully-precessing NR-surrogate templates and freely provide the resulting posterior samples as supplemental material. We find that our refined analysis is able to extract information from sub-dominant GW harmonics in data, and therefore better resolve the degeneracy in measuring source luminosity distance and orbital inclination for both events. We estimate the sources of both events to be 20-25% further away than was previously estimated. Our analyses also constrain their orbital orientation more tightly around face-on or face-off configurations than before. Additionally, for GW150914 we constrain the effective inspiral spin more tightly around zero. This work is one of the first to unambiguously extract sub-dominant GW mode information from real events. It is also a first step toward eliminating the approximations used in semi-analytic waveform models from GW parameter estimation. It also motivates that NR surrogates be extended to cover more of the binary black hole parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GipsyX/RTGx, A New Tool Set for Space Geodetic Operations and Research

GipsyX/RTGx is the Jet Propulsion Laboratory&#39;s (JPL) next generation software package for positioning, navigation, timing, and Earth science using measurements from three geodetic techniques: Global Navigation Satellite Systems (GNSS), Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS); with Very Long Baseline Interferometry (VLBI) under development. The software facilitates combined estimation of geodetic and geophysical parameters using a Kalman filter approach on real or simulated data in both post-processing and in real-time. The estimated parameters include station coordinates and velocities, satellite orbits and clocks, Earth orientation, ionospheric and tropospheric delays. The software is also capable of full realization of a dynamic terrestrial reference through analysis and combination of time series of ground station coordinates. We present some key aspects of its new architecture, and describe some of its major applications, including Real-time orbit determination and ephemeris predictions in the U.S. Air Force Next Generation GPS Operational Control Segment (OCX), as well as in JPL&#39;s Global Differential GPS (GDGPS) System, supporting User Range Error (URE) of $<$ 5 cm RMS; precision post-processing GNSS orbit determination, including JPL&#39;s contributions to the International GNSS Service (IGS) with URE in the 2 cm RMS range; Precise point positioning (PPP) with ambiguity resolution, both statically and kinematically, for geodetic applications with 2 mm horizontal, and 6.5 mm vertical repeatability for static positioning; Operational orbit and clock determination for Low Earth Orbiting (LEO) satellites, such as NASA&#39;s Gravity Recovery and Climate Experiment (GRACE) mission with GRACE relative clock alignment at the 20 ps level.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Key Elements of Robustness in Binary Black Hole Evolutions using Spectral Methods

As a network of advanced-era gravitational wave detectors is nearing its design sensitivity, efficient and accurate waveform modeling becomes more and more relevant. Understanding of the nature of the signal being sought can have an order unity effect on the event rates seen in these instruments. The paper provides a description of key elements of the Spectral Einstein Code ({\tt SpEC}), with details of our spectral adaptive mesh refinement (AMR) algorithm that has been optimized for binary black hole (BBH) evolutions. We expect that the gravitational waveform catalog produced by our code will have a central importance in both the detection and parameter estimation of gravitational waves in these instruments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Neutron star-black hole mergers with a nuclear equation of state and neutrino cooling: Dependence in the binary parameters

We present a first exploration of the results of neutron star-black hole mergers using black hole masses in the most likely range of $7M_\odot-10M_\odot$, a neutrino leakage scheme, and a modeling of the neutron star material through a finite-temperature nuclear-theory based equation of state. In the range of black hole spins in which the neutron star is tidally disrupted ($χ_{\rm BH}\gtrsim 0.7$), we show that the merger consistently produces large amounts of cool ($T\lesssim 1\,{\rm MeV}$), unbound, neutron-rich material ($M_{\rm ej}\sim 0.05M_\odot-0.20M_\odot$). A comparable amount of bound matter is initially divided between a hot disk ($T_{\rm max}\sim 15\,{\rm MeV}$) with typical neutrino luminosity $L_ν\sim 10^{53}\,{\rm erg/s}$, and a cooler tidal tail. After a short period of rapid protonization of the disk lasting $\sim 10\,{\rm ms}$, the accretion disk cools down under the combined effects of the fall-back of cool material from the tail, continued accretion of the hottest material onto the black hole, and neutrino emission. As the temperature decreases, the disk progressively becomes more neutron-rich, with dimmer neutrino emission. This cooling process should stop once the viscous heating in the disk (not included in our simulations) balances the cooling. These mergers of neutron star-black hole binaries with black hole masses $M_{\rm BH}\sim 7M_\odot-10M_\odot$ and black hole spins high enough for the neutron star to disrupt provide promising candidates for the production of short gamma-ray bursts, of bright infrared post-merger signals due to the radioactive decay of unbound material, and of large amounts of r-process nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Black hole-neutron star mergers at realistic mass ratios: Equation of state and spin orientation effects

Black hole-neutron star mergers resulting in the disruption of the neutron star and the formation of an accretion disk and/or the ejection of unbound material are prime candidates for the joint detection of gravitational-wave and electromagnetic signals when the next generation of gravitational-wave detectors comes online. However, the disruption of the neutron star and the properties of the post-merger remnant are very sensitive to the parameters of the binary. In this paper, we study the impact of the radius of the neutron star and the alignment of the black hole spin for systems within the range of mass ratio currently deemed most likely for field binaries (M_BH ~ 7 M_NS) and for black hole spins large enough for the neutron star to disrupt (J/M^2=0.9). We find that: (i) In this regime, the merger is particularly sensitive to the radius of the neutron star, with remnant masses varying from 0.3M_NS to 0.1M_NS for changes of only 2 km in the NS radius; (ii) 0.01-0.05M_sun of unbound material can be ejected with kinetic energy >10^51 ergs, a significant increase compared to low mass ratio, low spin binaries. This ejecta could power detectable optical and radio afterglows. (iii) Only a small fraction (<3%) of the Advanced LIGO events in this parameter range have gravitational-wave signals which could offer constraints on the equation of state of the neutron star. (iv) A misaligned black hole spin works against disk formation, with less neutron star material remaining outside of the black hole after merger, and a larger fraction of that material remaining in the tidal tail instead of the forming accretion disk. (v) Large kicks (v>300 km/s) can be given to the final black hole as a result of a precessing BHNS merger, when the disruption of the neutron star occurs just outside or within the innermost stable spherical orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First direct comparison of non-disrupting neutron star-black hole and binary black hole merger simulations

We present the first direct comparison of numerical simulations of neutron star-black hole and black hole-black hole mergers in full general relativity. We focus on a configuration with non spinning objects and within the most likely range of mass ratio for neutron star-black hole systems (q=6). In this region of the parameter space, the neutron star is not tidally disrupted prior to merger, and we show that the two types of mergers appear remarkably similar. The effect of the presence of a neutron star on the gravitational wave signal is not only undetectable by the next generation of gravitational wave detectors, but also too small to be measured in the numerical simulations: even the plunge, merger and ringdown signals appear in perfect agreement for both types of binaries. The characteristics of the post-merger remnants are equally similar, with the masses of the final black holes agreeing within dM< 5 10^{-4}M_BH and their spins within da< 10^{-3}M_BH. The rate of periastron advance in the mixed binary agrees with previously published binary black hole results, and we use the inspiral waveforms to place constraints on the accuracy of our numerical simulations independent of algorithmic choices made for each type of binary. Overall, our results indicate that non-disrupting neutron star-black hole mergers are exceptionally well modeled by black hole-black hole mergers, and that given the absence of mass ejection, accretion disk formation, or differences in the gravitational wave signals, only electromagnetic precursors could prove the presence of a neutron star in low-spin systems of total mass ~10Msun, at least until the advent of gravitational wave detectors with a sensitivity comparable to that of the proposed Einstein Telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Massive disk formation in the tidal disruption of a neutron star by a nearly extremal black hole

Black hole-neutron star (BHNS) binaries are important sources of gravitational waves for second-generation interferometers, and BHNS mergers are also a proposed engine for short, hard gamma-ray bursts. The behavior of both the spacetime (and thus the emitted gravitational waves) and the neutron star matter in a BHNS merger depend strongly and nonlinearly on the black hole&#39;s spin. While there is a significant possibility that astrophysical black holes could have spins that are nearly extremal (i.e. near the theoretical maximum), to date fully relativistic simulations of BHNS binaries have included black-hole spins only up to $S/M^2$=0.9, which corresponds to the black hole having approximately half as much rotational energy as possible, given the black hole&#39;s mass. In this paper, we present a new simulation of a BHNS binary with a mass ratio $q=3$ and black-hole spin $S/M^2$=0.97, the highest simulated to date. We find that the black hole&#39;s large spin leads to the most massive accretion disk and the largest tidal tail outflow of any fully relativistic BHNS simulations to date, even exceeding the results implied by extrapolating results from simulations with lower black-hole spin. The disk appears to be remarkably stable. We also find that the high black-hole spin persists until shortly before the time of merger; afterwards, both merger and accretion spin down the black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Suitability of hybrid gravitational waveforms for unequal-mass binaries

This article studies sufficient accuracy criteria of hybrid post-Newtonian (PN) and numerical relativity (NR) waveforms for parameter estimation of strong binary black-hole sources in second- generation ground-based gravitational-wave detectors. We investigate equal-mass non-spinning binaries with a new 33-orbit NR waveform, as well as unequal-mass binaries with mass ratios 2, 3, 4 and 6. For equal masses, the 33-orbit NR waveform allows us to recover previous results and to extend the analysis toward matching at lower frequencies. For unequal masses, the errors between different PN approximants increase with mass ratio. Thus, at 3.5PN, hybrids for higher-mass-ratio systems would require NR waveforms with many more gravitational-wave (GW) cycles to guarantee no adverse impact on parameter estimation. Furthermore, we investigate the potential improvement in hybrid waveforms that can be expected from 4th order post-Newtonian waveforms, and find that knowledge of this 4th post-Newtonian order would significantly improve the accuracy of hybrid waveforms.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Are different approaches to constructing initial data for binary black hole simulations of the same astrophysical situation equivalent?

Initial data for numerical evolutions of binary-black holes have been dominated by &#34;conformally flat&#34; (CF) data (i.e., initial data where the conformal background metric is chosen to be flat) because they are easy to construct. However, CF initial data cannot simulate nearly extremal spins, while more complicated &#34;conformally curved&#34; initial data (i.e., initial data in which the background metric is \emph{not} explicitly chosen to be flat), such as initial data where the spatial metric is chosen to be proportional to a weighted superposition of two Kerr-Schild (SKS) black holes can. Here we establish the consistency between the astrophysical results of these two initial data schemes for nonspinning binary systems. We evolve the inspiral, merger, and ringdown of two equal-mass, nonspinning black holes using SKS initial data and compare with an analogous simulation using CF initial data. We find that the resultant gravitational-waveform phases agree to within $δϕ\lesssim 10^{-2}$ radians and the amplitudes agree to within $δA/A \lesssim 5 \times 10^{-3}$, which are within the numerical errors of the simulations. Furthermore, we find that the final mass and spin of the remnant black hole agree to one part in $10^{5}

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Black hole-neutron star mergers for 10 solar mass black holes

General relativistic simulations of black hole-neutron star mergers have currently been limited to low-mass black holes (less than 7 solar mass), even though population synthesis models indicate that a majority of mergers might involve more massive black holes (10 solar mass or more). We present the first general relativistic simulations of black hole-neutron star mergers with 10 solar mass black holes. For massive black holes, the tidal forces acting on the neutron star are usually too weak to disrupt the star before it reaches the innermost stable circular orbit of the black hole. Varying the spin of the black hole in the range a/M = 0.5-0.9, we find that mergers result in the disruption of the star and the formation of a massive accretion disk only for large spins a/M>0.7-0.9. From these results, we obtain updated constraints on the ability of BHNS mergers to be the progenitors of short gamma-ray bursts as a function of the mass and spin of the black hole. We also discuss the dependence of the gravitational wave signal on the black hole parameters, and provide waveforms and spectra from simulations beginning 7-8 orbits before merger.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations of non-equal mass black hole binaries with spectral methods

This paper presents techniques and results for simulations of unequal mass, non-spinning black hole binaries with pseudo-spectral methods. Specifically, we develop an efficient root-finding procedure to ensure the black hole initial data have the desired masses and spins, we extend the dual coordinate frame method and eccentricity removal to asymmetric binaries. Furthermore, we describe techniques to simulate mergers of unequal mass black holes. The second part of the paper presents numerical simulations of non-spinning black hole binaries with mass ratios 2, 3, 4 and 6, covering between 15 and 22 orbits, merger and ringdown. We discuss the accuracy of these simulations, the evolution of the (initially zero) black hole spins, and the remnant black hole properties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accurate gravitational waveforms for binary-black-hole mergers with nearly extremal spins

Motivated by the possibility of observing gravitational waves from merging black holes whose spins are nearly extremal (i.e., 1 in dimensionless units), we present numerical waveforms from simulations of merging black holes with the highest spins simulated to date: (1) a 25.5-orbit inspiral, merger, and ringdown of two holes with equal masses and spins of magnitude 0.97 aligned with the orbital angular momentum; and (2) a previously reported 12.5-orbit inspiral, merger, and ringdown of two holes with equal masses and spins of magnitude 0.95 anti-aligned with the orbital angular momentum. First, we consider the horizon mass and spin evolution of the new aligned-spin simulation. During the inspiral, the horizon area and spin evolve in remarkably close agreement with Alvi&#39;s analytic predictions, and the remnant hole&#39;s final spin agrees reasonably well with several analytic predictions. We also find that the total energy emitted by a real astrophysical system with these parameters---almost all of which is radiated during the time included in this simulation---would be 10.952% of the initial mass at infinite separation. Second, we consider the gravitational waveforms for both simulations. After estimating their uncertainties, we compare the waveforms to several post-Newtonian approximants, finding significant disagreement well before merger, although the phase of the TaylorT4 approximant happens to agree remarkably well with the numerical prediction in the aligned-spin case. We find that the post-Newtonian waveforms have sufficient uncertainty that hybridized waveforms will require far longer numerical simulations (in the absence of improved post-Newtonian waveforms) for accurate parameter estimation of low-mass binary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Spin Diagrams for Equal-Mass Black-Hole Binaries with Aligned Spins

Binary black-hole systems with spins aligned with the orbital angular momentum are of special interest as they may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins aligned with the orbital angular momentum but otherwise arbitrary. By least-square fitting the results of these simulations we have constructed two &#34;spin diagrams&#34; which provide straightforward information about the recoil velocity |v_kick| and the final black-hole spin a_fin in terms of the dimensionless spins a_1 and a_2 of the two initial black holes. Overall they suggest a maximum recoil velocity of |v_kick|=441.94 km/s, and minimum and maximum final spins a_fin=0.3471 and a_fin=0.9591, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Some mathematical problems in numerical relativity

The main goal of numerical relativity is the long time simulation of highly nonlinear spacetimes that cannot be treated by perturbation theory. This involves analytic, computational and physical issues. At present, the major impasses to achieving global simulations of physical usefulness are of an analytic/computational nature. We present here some examples of how analytic insight can lend useful guidance for the improvement of numerical approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Well-Posed Initial-Boundary Evolution in General Relativity

Maximally dissipative boundary conditions are applied to the initial-boundary value problem for Einstein&#39;s equations in harmonic coordinates to show that it is well-posed for homogeneous boundary data and for boundary data that is small in a linearized sense. The method is implemented as a nonlinear evolution code which satisfies convergence tests in the nonlinear regime and is robustly stable in the weak field regime. A linearized version has been stably matched to a characteristic code to compute the gravitational waveform radiated to infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Boundary conditions in linearized harmonic gravity

We investigate the initial-boundary value problem for linearized gravitational theory in harmonic coordinates. Rigorous techniques for hyperbolic systems are applied to establish well-posedness for various reductions of the system into a set of six wave equations. The results are used to formulate computational algorithms for Cauchy evolution in a 3-dimensional bounded domain. Numerical codes based upon these algorithms are shown to satisfy tests of robust stability for random constraint violating initial data and random boundary data; and shown to give excellent performance for the evolution of typical physical data. The results are obtained for plane boundaries as well as piecewise cubic spherical boundaries cut out of a Cartesian grid.

0 citations0 reviewsNo code linkedOpen access