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

Tim Dietrich

Tim Dietrich contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEgr-qcastro-ph.COastro-ph.SRnucl-thastro-ph.IMhep-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
8topics
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

33 published item(s)

preprint2026arXiv

Black hole spectroscopy of collapsing and merging neutron stars

Black hole spectroscopy is an important pillar when studying gravitational waves from black holes and enables tests of general relativity. Most of the gravitational-wave signals observed over the last decade originate from binary black hole systems. Binary neutron star or black hole-neutron star systems are rarer but of particular interest for the next-generation ground-based gravitational-wave detectors. These events offer the exciting possibility of studying matter effects on the ringdown of "dirty black holes". In this work, we ask the question: Does matter matter? Using numerical-relativity, we simulate a wide range of collapsing neutron stars producing matter environments, both in isolated scenarios and in binary mergers. Qualitatively, the resulting ringdown signals can be classified into "clean", "modified", and "distorted" cases, depending on the amount of matter that is present. We apply standard strategies for extracting quasinormal modes of clean signals, using both theory-agnostic and theory-specific assumptions. Even in the presence of matter, possible modifications of quasinormal modes seem to be dominated by ringdown modeling systematics. We find that incorporating multiple quasinormal modes allows one to drastically reduce mismatches and errors in estimating the final black hole mass at early times. If not treated carefully, deviations in the fundamental quasinormal mode might artificially be overestimated and falsely attributed to the presence of matter or violations of general relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Challenging a binary neutron star merger interpretation of GW230529

GW230529_181500 represented the first gravitational-wave detection with one of the component objects' mass inferred to lie in the previously hypothesized mass gap between the heaviest neutron stars and the lightest observed black holes. Given the expected maximum mass values for neutron stars, this object was identified as a black hole, and, with the secondary component being a neutron star, the detection was classified as a neutron star-black hole merger. However, due to the low signal-to-noise ratio and the known waveform degeneracy between the spin and mass ratio in the employed gravitational-wave models, GW230529_181500 could also be interpreted as a merger of two heavy ($\gtrsim 2 \mathrm{M}_\odot$) neutron stars with high spins. We investigate the distinguishability of these scenarios by performing parameter estimation on simulated signals obtained from numerical-relativity waveforms for both neutron star-black hole and binary neutron star systems, with parameters consistent with GW230529_181500, and comparing them to the analysis of the real event data. We find that GW230529_181500 is more likely to have originated from a neutron star-black hole merger, though the possibility of a binary neutron star origin can not be ruled out. Moreover, we use the simulation data to estimate the signatures of potential electromagnetic counterparts emitted by the systems. We find them to be too dim to be located by current wide-field surveys if only the dynamical ejecta is considered, and detectable by the Vera C. Rubin Observatory during the first two days after merger if one accounts for additional disk wind ejecta.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Detecting electromagnetic counterparts to LIGO/Virgo/KAGRA gravitational wave events with DECam: Neutron Star Mergers

With GW170817 being the only multimessenger gravitational wave (GW) event with an associated kilonova detected so far, there exists a pressing need for realistic estimation of the GW localization uncertainties and rates, as well as optimization of available telescope time to enable the detection of new kilonovae. We simulate GW events assuming a data-driven distribution of binary parameters for the LIGO/Virgo/KAGRA (LVK) fourth and fifth observing runs (O4 and O5). We map the binary neutron star (BNS) and neutron star-black hole (NSBH) properties to the kilonova optical light curves. We use the simulated population of kilonovae to generate follow-up observing plans, with the primary goal of optimizing detection with the Gravitational Wave Multi-Messenger Astronomy DECam Survey (GW-MMADS). We explore the dependence of kilonova detectability on the mass, distance, inclination, and spin of the binaries. Assuming that no BNS was detected during O4 until the end of 2024, we present updated GW BNS (NSBH) merger detection rates. We expect to detect BNS (NSBH) kilonovae with DECam at a per-year rate of: $0$-$2.0$ ($0$) in O4, and $2.0$-$19$ ($0$-$1.0$) in O5. We expect the majority of BNS detections and also those accompanied by a detectable kilonova to produce a hypermassive NS remnant, with a significant fraction of the remaining BNSs promptly collapsing to a BH. We release GW simulations and depths required to detect kilonovae based on our predictions to support the astronomical community in their multimessenger follow-up campaigns and analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

ESO Expanding Horizon White Paper: Revealing the properties of matter at supranuclear densities with gravitational waves

Understanding dense matter under extreme conditions is one of the most fundamental puzzles in modern physics. Complex interactions give rise to emergent, collective phenomena. While nuclear experiments and Earth - based colliders provide valuable insights, much of the quantum chromodynamics phase diagram at high density and low temperature remains accessible only through astrophysical observations of neutron stars, neutron star mergers, and stellar collapse. Astronomical observations thus offer a direct window to the physics on subatomic scales with gravitational waves presenting an especially clean channel. Next-generation gravitational - wave observatories, such as the Einstein Telescope, would serve as unparalleled instruments to transform our understanding of neutron star matter. They will enable the detection of up to tens of thousands of binary neutron star and neutron star - black hole mergers per year, a dramatic increase over the few events accessible with current detectors. They will provide an unprecedented precision in probing cold, dense matter during the binary inspiral, exceeding by at least an order of magnitude what current facilities can achieve. Moreover, these observatories will allow us to explore uncharted regimes of dense matter at finite temperatures produced in a subset of neutron star mergers, areas that remain entirely inaccessible to current instruments. Together with multimessenger observations, these measurements will significantly deepen our knowledge of dense nuclear matter.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Efficient Bayesian analysis of kilonovae and gamma ray burst afterglows with fiesta

Gamma-ray burst (GRB) afterglows and kilonovae (KNe) are electromagnetic transients that can accompany binary neutron star (BNS) mergers. Therefore, studying their emission processes is of general interest for constraining cosmological parameters or the behavior of ultra-dense matter. One common method to analyze electromagnetic data from BNS mergers is to sample a Bayesian posterior over the parameters of a physical model for the transient. However, sampling the posterior is computationally costly and because of the many likelihood evaluations required in this process, detailed models are too expensive to be used directly in Bayesian inference. In this paper, we address the problem by introducing fiesta, a python package to train machine learning (ML) surrogates for GRB afterglow and kilonova models that have the capacity to accelerate likelihood evaluations. Specifically, we introduce extensive ML surrogates for the state-of-the-art GRB afterglow models afterglowpy and pyblastafterglow, along with a new surrogate for KN emission based on the possis code. Our surrogates enable evaluation of the light-curve posterior within minutes. We also provide built-in posterior sampling capabilities in fiesta that rely on the flowMC package, which efficiently scale to higher dimensions when adding up to tens of nuisance sampling parameters. Because of its use of the JAX framework, fiesta also allows for GPU acceleration during both surrogate training and posterior sampling. We applied our framework to reanalyze AT2017gfo/GRB170817A and GRB211211A with our surrogates, thus employing the new pyblastafterglow model for the first time in Bayesian inference.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a $1.4M_\odot$ neutron star to be $R=11.98^{+0.35}_{-0.40}$km.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Advancing the Landscape of Multimessenger Science in the Next Decade

The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Elliptica: a new pseudo-spectral code for the construction of initial data

Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudo-spectral code Elliptica, an infrastructure for construction of initial data for various binary and single gravitational systems of all kinds. The elliptic equations under consideration are solved on a single spatial hypersurface of the spacetime manifold. Using coordinate maps, the hypersurface is covered by patches whose boundaries can adapt to the surface of the compact objects. To solve elliptic equations with arbitrary boundary condition, Elliptica deploys a Schur complement domain decomposition method with a direct solver. In this version, we use cubed sphere coordinate maps and the fields are expanded using Chebyshev polynomials of the first kind. Here, we explain the building blocks of Elliptica and the initial data construction algorithm for a black hole-neutron star binary system. We perform convergence tests and evolve the data to validate our results. Within our framework, the neutron star can reach spin values close to breakup with arbitrary direction, while the black hole can have arbitrary spin with dimensionless spin magnitude $\sim 0.8$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

High-accuracy high-mass ratio simulations for binary neutron stars and their comparison to existing waveform models

The subsequent observing runs of the advanced gravitational-wave detector network will likely provide us with various gravitational-wave observations of binary neutron star systems. For an accurate interpretation of these detections, we need reliable gravitational-wave models. To test and to point out how existing models could be improved, we perform a set of high-resolution numerical-relativity simulations for four different physical setups with mass ratios $q$ = $1.25$, $1.50$, $1.75$, $2.00$, and total gravitational mass $M = 2.7M_\odot$ . Each configuration is simulated with five different resolutions to allow a proper error assessment. Overall, we find approximately 2nd order converging results for the dominant $(2,2)$, but also subdominant $(2,1)$, $(3,3)$, $(4,4)$ modes, while, generally, the convergence order reduces slightly for an increasing mass ratio. Our simulations allow us to validate waveform models, where we find generally good agreement between state-of-the-art models and our data, and to prove that scaling relations for higher modes currently employed for binary black hole waveform modeling also apply for the tidal contribution. Finally, we also test if the current NRTidal model to describe tidal effects is a valid description for high-mass ratio systems. We hope that our simulation results can be used to further improve and test waveform models in preparation for the next observing runs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

How does dark matter affect compact star properties and high density constraints of strongly interacting matter

We study the impact of asymmetric bosonic dark matter on neutron star properties, including possible changes of tidal deformability, maximum mass, radius, and matter distribution inside the star. The conditions at which dark matter particles tend to condensate in the star's core or create an extended halo are presented. We show that dark matter condensed in a core leads to a decrease of the total gravitational mass and tidal deformability compared to a pure baryonic star, which we will perceive as an effective softening of the equation of state. On the other hand, the presence of a dark matter halo increases those observable quantities. Thus, observational data on compact stars could be affected by accumulated dark matter and, consequently, constraints we put on strongly interacting matter at high densities. To confirm the presence of dark matter in the compact star's interior, and to break the degeneracy between the effect of accumulated dark matter and strongly interacting matter properties at high densities, several astrophysical and GW tests are proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Incorporating a radiative hydrodynamics scheme in the numerical-relativity code BAM

To study binary neutron star systems and to interpret observational data such as gravitational-wave and kilonova signals, one needs an accurate description of the processes that take place during the final stages of the coalescence, e.g., through numerical-relativity simulations. In this work, we present an updated version of the numerical-relativity code BAM in order to incorporate nuclear-theory based Equations of State and a simple description of neutrino interactions through a Neutrino Leakage Scheme. Different test simulations, for stars undergoing a neutrino-induced gravitational collapse and for binary neutron stars systems, validate our new implementation. For the binary neutron stars systems, we show that we can evolve stably and accurately distinct microphysical models employing the different equations of state: SFHo, DD2, and the hyperonic BHB$Λϕ$. Overall, our test simulations have good agreement with those reported in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multi-messenger constraints on the Hubble constant through combination of gravitational waves, gamma-ray bursts and kilonovae from neutron star mergers

The simultaneous detection of gravitational waves and light from the binary neutron star merger GW170817 led to independent measurements of distance and redshift, providing a direct estimate of the Hubble constant $H_0$ that does not rely on a cosmic distance ladder nor assumes a specific cosmological model. By using gravitational waves as ''standard sirens'', this approach holds promise to arbitrate the existing tension between the $H_0$ value inferred from the cosmic microwave background and those obtained from local measurements. However, the known degeneracy in the gravitational-wave analysis between distance and inclination of the source lead to a $H_0$ value from GW170817 that was not precise enough to resolve the existing tension. In this review, we summarize recent works exploiting the viewing-angle dependence of the electromagnetic signal, namely the associated short gamma-ray burst and kilonova, to constrain the system inclination and improve on $H_0$. We outline the key ingredients of the different methods, summarize the results obtained in the aftermath of GW170817 and discuss the possible systematics introduced by each of these methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical Simulations of Dark Matter Admixed Neutron Star Binaries

Multi-messenger observations of compact binary mergers provide a new way to constrain the nature of dark matter that may accumulate in and around neutron stars. In this article, we extend the infrastructure of our numerical-relativity code BAM to enable the simulation of neutron stars that contain an additional mirror dark matter component. We perform single star tests to verify our code and the first binary neutron star simulations of this kind. We find that the presence of dark matter reduces the lifetime of the merger remnant and favors a prompt collapse to a black hole. Furthermore, we find differences in the merger time for systems with the same total mass and mass ratio, but different amounts of dark matter. Finally, we find that electromagnetic signals produced by the merger of binary neutron stars admixed with dark matter are very unlikely to be as bright as their dark matter free counterparts. Given the increasing sensitivity of multi-messenger facilities, our analysis gives a new perspective on how to probe the presence of dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantifying modeling uncertainties when combining multiple gravitational-wave detections from binary neutron star sources

With the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure that uncertainties inherent in the gravitational-wave models will not lead to systematic biases when information from multiple detections are combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyse them with a set of four different waveform models. Based on our analysis with about 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to $\pm 250\rm m$ ($2\%$ at $90\%$ credible interval) but that systematic differences between currently employed waveform models can be twice as large. Hence, it will be essential to ensure that systematic biases will not become dominant in inferences of the neutron-star equation of state when capitalizing on future developments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An aligned-spin neutron-star--black-hole waveform model based on the effective-one-body approach and numerical-relativity simulations

After the discovery of gravitational waves from binary black holes (BBHs) and binary neutron stars (BNSs) with the LIGO and Virgo detectors, neutron-star--black-holes (NSBHs) are the natural next class of binary systems to be observed. In this work, we develop a waveform model for aligned-spin neutron-star--black-holes (NSBHs) combining a BBH baseline waveform (available in the effective-one-body approach) with a phenomenological description of tidal effects (extracted from numerical-relativity simulations), and correcting the amplitude during the late inspiral, merger and ringdown to account for the NS tidal disruption. We calibrate the amplitude corrections using NSBH waveforms obtained with the SpEC and the SACRA codes. The model was calibrated using simulations with NS masses in the range $1.2-1.4 M_\odot$, tidal deformabilities up to $4200$ (for a 1.2 $M_\odot$ NS), and dimensionless BH spin magnitude up to 0.9. Based on the simulations used, and on checking that sensible waveforms are produced, we recommend our model to be employed with NS mass in the range $1\mbox{--}3 M_\odot$, tidal deformability $0\mbox{--}5000$, and BH spin magnitude up to $0.9$. We also validate our model against two new, highly accurate NSBH waveforms with BH spin 0.9 and mass ratios 3 and 4, characterized by tidal disruption, produced with SpEC, and find very good agreement. We find that it will be challenging for the advanced LIGO-Virgo--detector network at design sensitivity to distinguish different source classes. We perform parameter-estimation on a synthetic numerical-relativity signal in zero noise to study parameter biases. Finally, we reanalyze GW170817, with the hypothesis that it is a NSBH. We do not find evidence to distinguish the BNS and NSBH hypotheses, however the posterior for the mass ratio is shifted to less equal masses under the NSBH hypothesis. [Abstract abridged for arxiv].

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Data-driven expectations for electromagnetic counterpart searches based on LIGO/Virgo public alerts

Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in 2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared kilonova, and panchromatic afterglow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo's third observing run, O3, were much larger than predicted. We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly localized events have increased the overall number of detections, of which well-localized, gold-plated events make up a smaller proportion overall. We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics of O3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and timely forecasting of gravitational-wave localization accuracy is paramount given the large commitments of telescope time and the need to prioritize which events are followed up. We include a data release of our simulated localizations as a public proposal planning resource for astronomers.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Discovery and confirmation of the shortest gamma ray burst from a collapsar

Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the universe. The duration and hardness distribution of GRBs has two clusters, now understood to reflect (at least) two different progenitors. Short-hard GRBs (SGRBs; T90 <2 s) arise from compact binary mergers, while long-soft GRBs (LGRBs; T90 >2 s) have been attributed to the collapse of peculiar massive stars (collapsars). The discovery of SN 1998bw/GRB 980425 marked the first association of a LGRB with a collapsar and AT 2017gfo/GRB 170817A/GW170817 marked the first association of a SGRB with a binary neutron star merger, producing also gravitational wave (GW). Here, we present the discovery of ZTF20abwysqy (AT2020scz), a fast-fading optical transient in the Fermi Satellite and the InterPlanetary Network (IPN) localization regions of GRB 200826A; X-ray and radio emission further confirm that this is the afterglow. Follow-up imaging (at rest-frame 16.5 days) reveals excess emission above the afterglow that cannot be explained as an underlying kilonova (KN), but is consistent with being the supernova (SN). Despite the GRB duration being short (rest-frame T90 of 0.65 s), our panchromatic follow-up data confirms a collapsar origin. GRB 200826A is the shortest LGRB found with an associated collapsar; it appears to sit on the brink between a successful and a failed collapsar. Our discovery is consistent with the hypothesis that most collapsars fail to produce ultra-relativistic jets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

High-accuracy simulations of highly spinning binary neutron star systems

With an increasing number of expected gravitational-wave detections of binary neutron star mergers, it is essential that gravitational-wave models employed for the analysis of observational data are able to describe generic compact binary systems. This includes systems in which the individual neutron stars are millisecond pulsars for which spin effects become essential. In this work, we perform numerical-relativity simulations of binary neutron stars with aligned and anti-aligned spins within a range of dimensionless spins of $χ\sim [-0.28,0.58]$. The simulations are performed with multiple resolutions, show a clear convergence order and, consequently, can be used to test existing waveform approximants. We find that for very high spins gravitational-wave models that have been employed for the interpretation of GW170817 and GW190425 are not capable of describing our numerical-relativity dataset. We verify through a full parameter estimation study in which clear biases in the estimate of the tidal deformability and effective spin are present. We hope that in preparation of the next gravitational-wave observing run of the Advanced LIGO and Advanced Virgo detectors our new set of numerical-relativity data can be used to support future developments of new gravitational-wave models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigating the mass-ratio dependence of the prompt-collapse threshold with numerical-relativity simulations

The next observing runs of advanced gravitational-wave detectors will lead to a variety of binary neutron star detections and numerous possibilities for multi-messenger observations of binary neutron star systems. In this context a clear understanding of the merger process and the possibility of prompt black hole formation after merger is important, as the amount of ejected material strongly depends on the merger dynamics. These dynamics are primarily affected by the total mass of the binary, however, the mass ratio also influences the postmerger evolution. To determine the effect of the mass ratio, we investigate the parameter space around the prompt-collapse threshold with a new set of fully relativistic simulations. The simulations cover three equations of state and seven mass ratios in the range of $1.0 \leq q \leq 1.75$, with five to seven simulations of binary systems of different total mass in each case. The threshold mass is determined through an empirical relation based on the collapse-time, which allows us to investigate effects of the mass-ratio on the threshold mass and also on the properties of the remnant system. Furthermore, we model effects of mass ratio and equation of state on tidal parameters of threshold configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Multi-messenger Bayesian parameter inference of a binary neutron-star merger

The combined detection of a binary neutron-star merger in both gravitational waves (GWs) and electromagnetic (EM) radiation spanning the entire spectrum -- GW170817 / AT2017gfo / GRB170817A -- marks a breakthrough in the field of multi-messenger astronomy. Between the plethora of modeling and observations, the rich synergy that exists among the available data sets creates a unique opportunity to constrain the binary parameters, the equation of state of supranuclear density matter, and the physical processes at work during the kilonova and gamma-ray burst. We report, for the first time, Bayesian parameter estimation combining information from GW170817, AT2017gfo, GRB170817 to obtain truly multi-messenger constraints on the tidal deformability $\tildeΛ \in [302,860]$, total binary mass $M \in [2.722,2.751] M_\odot$, the radius of a $1.4$ solar mass neutron star $R \in [11.3,13.5] \rm km$ (with additional $0.2\ \rm km$ systematic uncertainty), and an upper bound on the mass ratio of $q \leq 1.27$, all at 90% confidence. Our joint novel analysis makes use of new phenomenological descriptions of the dynamical ejecta, debris disk mass, and remnant black hole properties, all derived from a large suite of numerical relativity simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the nature of GW190814 and its impact on the understanding of supranuclear matter

The observation of a compact object with a mass of $2.50-2.67M_{\odot}$ on August 14, 2019, by the LIGO Scientific and Virgo collaborations (LVC) has the potential to improve our understanding of the supranuclear equation of state. While the gravitational-wave analysis of the LVC suggests that GW190814 likely was a binary black hole system, the secondary component could also have been the heaviest neutron star observed to date. We use our previously derived nuclear-physics-multimessenger astrophysics framework to address the nature of this object. Based on our findings, we determine GW190814 to be a binary black hole merger with a probability of $>99.9\%$. Even if we weaken previously employed constraints on the maximum mass of neutron stars, the probability of a binary black hole origin is still $\sim 81\%$. Furthermore, we study the impact that this observation has on our understanding of the nuclear equation of state by analyzing the allowed region in the mass-radius diagram of neutron stars for both a binary black hole or neutron star--black hole scenario. We find that the unlikely scenario in which the secondary object was a neutron star requires rather stiff equations of state with a maximum speed of sound $c_s\geq \sqrt{0.6}$ times the speed of light, while the binary black hole scenario does not offer any new insight.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Optimizing Cadences with Realistic Light Curve Filtering for Serendipitous Kilonova Discovery with Vera Rubin Observatory

Current and future optical and near-infrared wide-field surveys have the potential of finding kilonovae, the optical and infrared counterparts to neutron star mergers, independently of gravitational-wave or high-energy gamma-ray burst triggers. The ability to discover fast and faint transients such as kilonovae largely depends on the area observed, the depth of those observations, the number of re-visits per field in a given time frame, and the filters adopted by the survey; it also depends on the ability to perform rapid follow-up observations to confirm the nature of the transients. In this work, we assess kilonova detectability in existing simulations of the LSST strategy for the Vera C. Rubin Wide Fast Deep survey, with focus on comparing rolling to baseline cadences. Although currently available cadences can enable the detection of more than 300 kilonovae out to 1400 Mpc over the ten-year survey, we can expect only 3-32 kilonovae similar to GW170817 to be recognizable as fast-evolving transients. We also explore the detectability of kilonovae over the plausible parameter space, focusing on viewing angle and ejecta masses. We find that observations in redder izy bands are crucial for identification of nearby (within 300 Mpc) kilonovae that could be spectroscopically classified more easily than more distant sources. Rubin&#39;s potential for serendipitous kilonova discovery could be increased by gain of efficiency with the employment of individual 30s exposures (as opposed to 2x15s snap pairs), with the addition of red-band observations coupled with same-night observations in g- or r-bands, and possibly with further development of a new rolling-cadence strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constructing Love-Q-Relations with Gravitational Wave Detections

Quasi-universal relations connecting the tidal deformability and the quadrupole moment of individual neutron stars are predicted by theoretical computations, but have not been measured experimentally. However, such relations are employed during the interpretation of gravitational waves and, therefore, have a direct impact on the interpretation of real data. In this work, we study how quasi-universal relations can be tested and measured from gravitational wave signals connected to binary neutron star coalescences. We study a population of $120$ binary neutron star systems and find that Advanced LIGO and Advanced Virgo at design sensitivity could find possible deviations of predicted relations if the observed neutron stars are highly spinning. In the future, a network of third generation (3G) detectors will be able to even allow a measurement of quasi-universal relations. Thus, the outlined approach provides a new test of general relativity and nuclear physics predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Distinguishing high-mass binary neutron stars from binary black holes with second- and third-generation gravitational wave observatories

(Abridged) While the gravitational-wave (GW) signal GW170817 was accompanied by a variety of electromagnetic (EM) counterparts, sufficiently high-mass binary neutron star (BNS) mergers are expected to be unable to power bright EM counterparts. The putative high-mass binary BNS merger GW190425, for which no confirmed EM counterpart has been identified, may be an example of such a system. It is thus important to understand how well we will be able to distinguish high-mass BNSs and low-mass binary black holes (BBHs) solely from their GW signals. To do this, we consider the imprint of the tidal deformability of the neutron stars on the GW signal for systems undergoing prompt black hole formation after merger. We model the BNS signals using hybrid numerical relativity -- tidal effective-one-body waveforms. Specifically, we consider a set of five nonspinning equal-mass BNS signals with masses of 2.7, 3.0, 3.2 Msun and with three different equations of state, as well as the analogous BBH signals. We perform parameter estimation on these signals in three networks: Advanced LIGO-Advanced Virgo and Advanced LIGO-Advanced Virgo-KAGRA with sensitivities similar to O3 and O4, respectively, and a 3G network of two Cosmic Explorers (CEs) and one Einstein Telescope, with a CE sensitivity similar to Stage 2. Our analysis suggests that we cannot distinguish the signals from high-mass BNSs and BBHs at a 90% credible level with the O3-like network even at 40 Mpc. However, we can distinguish all but the most compact BNSs that we consider in our study from BBHs at 40 Mpc at a >= 95% credible level using the O4-like network, and can even distinguish them at a > 99.2% (>= 97%) credible level at 369 (835) Mpc using the 3G network. Additionally, we present a simple method to compute the leading effect of the Earth&#39;s rotation on the response of a gravitational wave detector in the frequency domain.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gravitational waves and mass ejecta from binary neutron star mergers: Effect of the spin orientation

We continue our study of the binary neutron star parameter space by investigating the effect of the spin orientation on the dynamics, gravitational wave emission, and mass ejection during the binary neutron star coalescence. We simulate seven different configurations using multiple resolutions to allow a reasonable error assessment. Due to the particular choice of the setups, five configurations show precession effects, from which two show a precession (&#34;wobbling&#34;) of the orbital plane, while three show a &#34;bobbing&#34; motion, i.e., the orbital angular momentum does not precess, while the orbital plane moves along the orbital angular momentum axis. Considering the ejection of mass, we find that precessing systems can have an anisotropic mass ejection, which could lead to a final remnant kick of $\sim 40 \rm km/s$ for the studied systems. Furthermore, for the chosen configurations, antialigned spins lead to larger mass ejecta than aligned spins, so that brighter electromagnetic counterparts could be expected for these configurations. Finally, we compare our simulations with the precessing, tidal waveform approximant IMRPhenomPv2_NRTidalv2 and find good agreement between the approximant and our numerical relativity waveforms with phase differences below 1.2 rad accumulated over the last $\sim$ 16 gravitational wave cycles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Implications of the search for optical counterparts during the second part of the Advanced LIGO&#39;s and Advanced Virgo&#39;s third observing run: lessons learned for future follow-up observations

Joint multi-messenger observations with gravitational waves and electromagnetic data offer new insights into the astrophysical studies of compact objects. The third Advanced LIGO and Advanced Virgo observing run began on April 1, 2019; during the eleven months of observation, there have been 14 compact binary systems candidates for which at least one component is potentially a neutron star. Although intensive follow-up campaigns involving tens of ground and space-based observatories searched for counterparts, no electromagnetic counterpart has been detected. Following on a previous study of the first six months of the campaign, we present in this paper the next five months of the campaign from October 2019 to March 2020. We highlight two neutron star - black hole candidates (S191205ah, S200105ae), two binary neutron star candidates (S191213g and S200213t) and a binary merger with a possible neutron star and a &#34;MassGap&#34; component, S200115j. Assuming that the gravitational-wave candidates are of astrophysical origin and their location was covered by optical telescopes, we derive possible constraints on the matter ejected during the events based on the non-detection of counterparts. We find that the follow-up observations during the second half of the third observing run did not meet the necessary sensitivity to constrain the source properties of the potential gravitational-wave candidate. Consequently, we suggest that different strategies have to be used to allow a better usage of the available telescope time. We examine different choices for follow-up surveys to optimize sky localization coverage vs.\ observational depth to understand the likelihood of counterpart detection.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Interpreting Binary Neutron Star Mergers: Describing the Binary Neutron Star Dynamics, Modelling Gravitational Waveforms, and Analyzing Detections

Gravitational waves emitted from the coalescence of neutron star binaries open a new window to probe matter and fundamental physics in unexplored, extreme regimes. To extract information about the supranuclear matter inside neutron stars and the properties of the compact binary systems, robust theoretical prescriptions are required. We give an overview about general features of the dynamics and the gravitational wave signal during the binary neutron star coalescence. We briefly describe existing analytical and numerical approaches to investigate the highly dynamical, strong-field region during the merger. We review existing waveform approximants and discuss properties and possible advantages and shortcomings of individual waveform models, and their application for real gravitational-wave data analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lensed or not lensed: Determining lensing magnifications for binary neutron star mergers from a single detection

Advanced LIGO and Advanced Virgo could observe the first lensed gravitational wave sources in the coming years, while the future Einstein Telescope could observe hundreds of lensed events. It is, therefore, crucial to develop methodologies to distinguish between lensed from unlensed gravitational-wave observations. A lensed signal not identified as such will lead to biases during the interpretation of the source. In particular, sources will appear to have intrinsically higher masses. No robust method currently exists to distinguish between the magnification bias caused by lensing and intrinsically high-mass sources. In this work, we show how to recognize lensed and unlensed binary neutron star systems through the measurement of their tidal effects for highly magnified sources as a proof-of-principle. The proposed method could be used to identify lensed binary neutron stars, which are the chief candidate for lensing cosmography studies. We apply our method on GW190425, finding no evidence in favor of lensing, mainly due to the poor measurement of the event&#39;s tidal effects. However, we expect that future detections with better tidal measurements can yield better constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring the Hubble Constant with a sample of kilonovae

Kilonovae produced by the coalescence of compact binaries with at least one neutron star are promising standard sirens for an independent measurement of the Hubble constant ($H_0$). Through their detection via follow-up of gravitational-wave (GW), short gamma-ray bursts (sGRBs) or optical surveys, a large sample of kilonovae (even without GW data) can be used for $H_0$ contraints. Here, we show measurement of $H_0$ using light curves associated with four sGRBs, assuming these are attributable to kilonovae, combined with GW170817. Including a systematic uncertainty on the models that is as large as the statistical ones, we find $H_0 = 73.8^{+6.3}_{-5.8}$\,$\mathrm{km}$ $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ and $H_0 = 71.2^{+3.2}_{-3.1}$\,$\mathrm{km}$ $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for two different kilonova models that are consistent with the local and inverse-distance ladder measurements. For a given model, this measurement is about a factor of 2-3 more precise than the standard-siren measurement for GW170817 using only GWs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modeling the gravitational wave signature of neutron star black hole coalescences: PhenomNSBH

Accurate gravitational-wave (GW) signal models exist for black hole binary (BBH) and neutron-star binary (BNS) systems, which are consistent with all of the published GW observations to date. Detections of a third class of compact-binary systems, neutron-star black hole (NSBH) binaries, have not yet been confirmed, but are eagerly awaited in the near future. For NSBH systems, GW models do not exist across the viable parameter space of signals. In this work we present the frequency-domain phenomenological model, PhenomNSBH, for GWs produced by NSBH systems with mass ratios from equal-mass up to 15, spin on the black hole up to a dimensionless spin of $|χ|=0.5$, and tidal deformabilities ranging from 0 (the BBH limit) to 5000. We extend previous work on a phenomenological amplitude model for NSBH systems to produce an amplitude model that is parameterized by a single tidal deformability parameter. This amplitude model is combined with an analytic phase model describing tidal corrections. The resulting approximant is compared to publicly-available NSBH numerical-relativity simulations and hybrid waveforms constructed from numerical-relativity simulations and tidal inspiral approximants. For most signals observed by second-generation ground-based detectors, it will be difficult to use the GW signal alone to distinguish single NSBH systems from either BNSs or BBHs, and therefore to unambiguously identify an NSBH system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multipolar Effective-One-Body Waveforms for Precessing Binary Black Holes: Construction and Validation

As gravitational-wave detectors become more sensitive, we will access a greater variety of signals emitted by compact binary systems, shedding light on their astrophysical origin and environment. A key physical effect that can distinguish among formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binary&#39;s orbital plane to precess. To accurately model such systems, it is crucial to include multipoles beyond the dominant quadrupole. Here, we develop the first multipolar precessing waveform model in the effective-one-body (EOB) formalism for the inspiral, merger and ringdown (IMR) of binary black holes: SEOBNRv4PHM. In the nonprecessing limit, the model reduces to SEOBNRv4HM, which was calibrated to numerical-relativity (NR) simulations, and waveforms from perturbation theory. We validate SEOBNRv4PHM by comparing it to the public catalog of 1405 precessing NR waveforms of the Simulating eXtreme Spacetimes (SXS) collaboration, and also to new 118 precessing NR waveforms, which span mass ratios 1-4 and spins up to 0.9. We stress that SEOBNRv4PHM is not calibrated to NR simulations in the precessing sector. We compute the unfaithfulness against the 1523 SXS precessing NR waveforms, and find that, for $94\%$ ($57\%$) of the cases, the maximum value, in the total mass range $20-200 M_\odot$, is below $3\%$ ($1\%$). Those numbers become $83\%$ ($20\%$) when using the IMR, multipolar, precessing phenomenological model IMRPhenomPv3HM. We investigate the impact of such unfaithfulness values with two parameter-estimation studies on synthetic signals. We also compute the unfaithfulness between those waveform models and identify in which part of the parameter space they differ the most. We validate them also against the multipolar, precessing NR surrogate model NRSur7dq4, and find that the SEOBNRv4PHM model outperforms IMRPhenomPv3HM.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On standardizing kilonovae and their use as standard candles to measure the Hubble constant

The detection of GW170817 is revolutionizing many areas of astrophysics with the joint observation of gravitational waves and electromagnetic emissions. These multi-messenger events provide a new approach to determine the Hubble constant, thus, they are a promising candidate for mitigating the tension between measurements of Type Ia supernovae via the local distance ladder and the Cosmic Microwave Background. In addition to the &#34;standard siren&#34; provided by the gravitational-wave measurement, the kilonova itself has characteristics that allow to improve existing measurements or to perform yet another, independent measurement of the Hubble constant without gravitational-wave information. Here, we employ standardization techniques borrowed from the type-Ia community and apply them to kilonovae, not using any information from the gravitational-wave signal. We use two versions of this technique, one derived from direct observables measured from the lightcurve, and the other based on inferred ejecta parameters, e.g., mass, velocity, and composition, for two different models. These lead to Hubble Constant measurements of $H_0 = 109^{+49}_{-35}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for the measured analysis, and $H_0 = 85^{+22}_{-17}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ and $H_0 = 79^{+23}_{-15}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for the inferred analyses. This measurement has error bars within ~$\sim$\,2 to the gravitational-wave measurements ($H_0=74^{+16}_{-8}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$), showing its promise as an independent constraint on $H_0$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Implications of the search for optical counterparts during the first six months of the Advanced LIGO&#39;s and Advanced Virgo&#39;s third observing run: possible limits on the ejecta mass and binary properties

GW170817 showed that neutron star mergers not only emit gravitational waves but also can release electromagnetic signatures in multiple wavelengths. Within the first half of the third observing run of the Advanced LIGO and Virgo detectors, there have been a number of gravitational wave candidates of compact binary systems for which at least one component is potentially a neutron star. In this article, we look at the candidates S190425z, S190426c, S190510g, S190901ap, and S190910h, predicted to have potentially a non-zero remnant mass, in more detail. All these triggers have been followed up with extensive campaigns by the astronomical community doing electromagnetic searches for their optical counterparts; however, according to the released classification, there is a high probability that some of these events might not be of extraterrestrial origin. Assuming that the triggers are caused by a compact binary coalescence and that the individual source locations have been covered during the EM follow-up campaigns, we employ three different kilonova models and apply them to derive possible constraints on the matter ejection consistent with the publicly available gravitational-wave trigger information and the lack of a kilonova detection. These upper bounds on the ejecta mass can be related to limits on the maximum mass of the binary neutron star candidate S190425z and to constraints on the mass-ratio, spin, and NS compactness for the potential black hole-neutron star candidate S190426c. Our results show that deeper electromagnetic observations for future gravitational wave events near the horizon limit of the advanced detectors are essential.

0 citations0 reviewsNo code linkedOpen access