Showing 1–50 of 219 results for author: Dietrich, T

Tests of General Relativity with GW230529: a neutron star merging with a lower mass-gap compact object

Authors: Elise M. Sänger, Soumen Roy, Michalis Agathos, Ofek Birnholtz, Alessandra Buonanno, Tim Dietrich, Maria Haney, Félix-Louis Julié, Geraint Pratten, Jan Steinhoff, Chris Van Den Broeck, Sylvia Biscoveanu, Prasanta Char, Anna Heffernan, Prathamesh Joshi, Atul Kedia, R. M. S. Schofield, M. Trevor, Michael Zevin

Abstract: On 29 May 2023, the LIGO Livingston observatory detected the gravitational-wave signal GW230529_181500 from the merger of a neutron star with a lower mass-gap compact object. Its long inspiral signal provides a unique opportunity to test General Relativity (GR) in a parameter space previously unexplored by strong-field tests. In this work, we performed parameterized inspiral tests of GR with GW230… ▽ More On 29 May 2023, the LIGO Livingston observatory detected the gravitational-wave signal GW230529_181500 from the merger of a neutron star with a lower mass-gap compact object. Its long inspiral signal provides a unique opportunity to test General Relativity (GR) in a parameter space previously unexplored by strong-field tests. In this work, we performed parameterized inspiral tests of GR with GW230529_181500. Specifically, we search for deviations in the frequency-domain GW phase by allowing for agnostic corrections to the post-Newtonian coefficients. We performed tests with the Flexible Theory Independent (FTI) and Test Infrastructure for General Relativity (TIGER) frameworks using several quasi-circular waveform models that capture different physical effects (higher modes, spins, tides). We find that the signal is consistent with GR for all deviation parameters. Assuming the primary object is a black hole, we obtain particularly tight constraints on the dipole radiation at $-1$PN order of $|δ\hat{\varphi}_{-2}| \lesssim 8 \times 10^{-5}$, which is a factor $\sim17$ times more stringent than previous bounds from the neutron star--black hole merger GW200115_042309, as well as on the 0.5PN and 1PN deviation parameters. We discuss some challenges that arise when analyzing this signal, namely biases due to correlations with tidal effects and the degeneracy between the 0PN deviation parameter and the chirp mass. To illustrate the importance of GW230529_181500 for tests of GR, we mapped the agnostic $-1$PN results to a class of Einstein-scalar-Gauss-Bonnet (ESGB) theories of gravity. We also conducted an analysis probing the specific phase deviation expected in ESGB theory and obtain an upper bound on the Gauss-Bonnet coupling of $\ell_{\rm GB} \lesssim 0.51~\rm{M}_\odot$ ($\sqrt{α_{\rm GB}} \lesssim 0.28$ km), which is better than any previously reported constraint. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 19 pages, 9 figures

Report number: LIGO-P2400200

What no one has seen before: gravitational waveforms from warp drive collapse

Authors: Katy Clough, Tim Dietrich, Sebastian Khan

Abstract: Despite originating in science fiction, warp drives have a concrete description in general relativity, with Alcubierre first proposing a spacetime metric that supported faster-than-light travel. Whilst there are numerous practical barriers to their implementation in real life, including a requirement for negative energy, computationally, one can simulate their evolution in time given an equation o… ▽ More Despite originating in science fiction, warp drives have a concrete description in general relativity, with Alcubierre first proposing a spacetime metric that supported faster-than-light travel. Whilst there are numerous practical barriers to their implementation in real life, including a requirement for negative energy, computationally, one can simulate their evolution in time given an equation of state describing the matter. In this work, we study the signatures arising from a warp drive "containment failure", assuming a stiff equation of state for the fluid. We compute the emitted gravitational-wave signal and track the energy fluxes of the fluid. Apart from its rather speculative application to the search for extraterrestrial life in gravitational-wave detector data, this work is interesting as a study of the dynamical evolution and stability of spacetimes that violate the null energy condition. Our work highlights the importance of exploring strange new spacetimes, to (boldly) simulate what no one has seen before. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 12 pages, 6 figures, plus appendix. Comments welcome!

The impact of gravitational waveform model systematics on the measurement of the Hubble Constant

Authors: Nina Kunert, Jonathan Gair, Peter T. H. Pang, Tim Dietrich

Abstract: Matching gravitational-wave observations of binary neutron stars with theoretical model predictions reveals important information about the sources, such as the masses and the distance to the stars. The latter can be used to determine the Hubble Constant, the rate at which the Universe expands. One general problem of all astrophysical measurements is that theoretical models only approximate the re… ▽ More Matching gravitational-wave observations of binary neutron stars with theoretical model predictions reveals important information about the sources, such as the masses and the distance to the stars. The latter can be used to determine the Hubble Constant, the rate at which the Universe expands. One general problem of all astrophysical measurements is that theoretical models only approximate the real underlying physics, which can lead to systematic uncertainties introducing biases. However, the extent of this bias for the distance measurement due to uncertainties of gravitational waveform models is unknown. In this study, we analyze a synthetic population of 38 binary neutron star sources measured with Advanced LIGO and Advanced Virgo at design sensitivity. We employ a set of four different waveform models and estimate model-dependent systematic biases on the extraction of the Hubble Constant using the bright siren method. Our results indicate that systematic biases are below statistical uncertainties for the current generation of gravitational wave detectors. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 9 pages, 3 figures

An Online Framework for Fitting Fast Transient Lightcurves

Authors: Tyler Barna, Brandon Reed, Igor Andreoni, Michael W. Coughlin, Tim Dietrich, Steven L. Groom, Theophile Jegou du Laz, Peter T. H Pang, Josiah N. Purdum, Ben Rusholme

Abstract: The identification of extragalactic fast optical transients (eFOTs) as potential multi-messenger sources is one of the main challenges in time-domain astronomy. However, recent developments have allowed for probes of rapidly-evolving transients. With the increasing number of alert streams from optical time-domain surveys, the next paradigm is building technologies to rapidly identify the most inte… ▽ More The identification of extragalactic fast optical transients (eFOTs) as potential multi-messenger sources is one of the main challenges in time-domain astronomy. However, recent developments have allowed for probes of rapidly-evolving transients. With the increasing number of alert streams from optical time-domain surveys, the next paradigm is building technologies to rapidly identify the most interesting transients for follow-up. One effort to make this possible is the fitting of objects to a variety of eFOT lightcurve models such as kilonovae and $γ$-ray burst (GRB) afterglows. In this work, we describe a new framework designed to efficiently fit transients to light curve models and flag them for further follow-up. We describe the pipeline's workflow and a handful of performance metrics, including the nominal sampling time for each model. We highlight as examples ZTF20abwysqy, the shortest long gamma ray burst discovered to date, and ZTF21abotose, a core-collapse supernova initially identified as a potential kilonova candidate. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: Accepted to MNRAS, 11 pages, 9 figures

Robust parameter estimation within minutes on gravitational wave signals from binary neutron star inspirals

Authors: Thibeau Wouters, Peter T. H. Pang, Tim Dietrich, Chris Van Den Broeck

Abstract: The gravitational waves emitted by binary neutron star inspirals contain information on nuclear matter above saturation density. However, extracting this information and conducting parameter estimation remains a computationally challenging and expensive task. Wong et al. introduced Jim arXiv:2302.05333, a parameter estimation pipeline that combines relative binning and jax features such as hardwar… ▽ More The gravitational waves emitted by binary neutron star inspirals contain information on nuclear matter above saturation density. However, extracting this information and conducting parameter estimation remains a computationally challenging and expensive task. Wong et al. introduced Jim arXiv:2302.05333, a parameter estimation pipeline that combines relative binning and jax features such as hardware acceleration and automatic differentiation into a normalizing flow-enhanced sampler for gravitational waves from binary black hole (BBH) mergers. In this work, we extend the Jim framework to analyze gravitational wave signals from binary neutron stars (BNS) mergers with tidal effects included. We demonstrate that Jim can be used for full Bayesian parameter estimation of gravitational waves from BNS mergers within a few tens of minutes, which includes the training of the normalizing flow and computing the reference parameters for relative binning. For instance, Jim can analyze GW170817 in 26 minutes (33 minutes) of total wall time using the TaylorF2 (IMRPhenomD_NRTidalv2) waveform, and GW190425 in around 21 minutes for both waveforms. We highlight the importance of such an efficient parameter estimation pipeline for several science cases as well as its ecologically friendly implementation of gravitational wave parameter estimation. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 17 pages, 7 figures. Code and data used to produce the figures in this paper is available at https://github.com/ThibeauWouters/TurboPE-BNS

Using Equation of State Constraints to Classify Low-Mass Compact Binary Mergers

Authors: Jacob Golomb, Isaac Legred, Katerina Chatziioannou, Adrian Abac, Tim Dietrich

Abstract: Compact objects observed via gravitational waves are classified as black holes or neutron stars primarily based on their inferred mass with respect to stellar evolution expectations. However, astrophysical expectations for the lowest mass range, $\lesssim 1.2 \,M_\odot$, are uncertain. If such low-mass compact objects exist, ground-based gravitational wave detectors may observe them in binary merg… ▽ More Compact objects observed via gravitational waves are classified as black holes or neutron stars primarily based on their inferred mass with respect to stellar evolution expectations. However, astrophysical expectations for the lowest mass range, $\lesssim 1.2 \,M_\odot$, are uncertain. If such low-mass compact objects exist, ground-based gravitational wave detectors may observe them in binary mergers. Lacking astrophysical expectations for classifying such observations, we go beyond the mass and explore the role of tidal effects. We evaluate how combined mass and tidal inference can inform whether each binary component is a black hole or a neutron star based on consistency with the supranuclear-density equation of state. Low-mass neutron stars experience a large tidal deformation; its observational identification (or lack thereof) can therefore aid in determining the nature of the binary components. Using simulated data, we find that the presence of a sub-solar mass neutron star (black hole) can be established with odds $\sim 100:1$ when two neutron stars (black holes) merge and emit gravitational waves at signal-to-noise ratio $\sim 20$. For the same systems, the absence of a black hole (neutron star) can be established with odds $\sim 10:1$. For mixed neutron star-black hole binaries, we can establish that the system contains a neutron star with odds $\gtrsim 5:1$. Establishing the presence of a black hole in mixed neutron star-black hole binaries is more challenging, except for the case of a $\lesssim 1\,M_{\odot}$ black hole with a $\gtrsim 1\,M_{\odot}$ neutron star companion. On the other hand, classifying each individual binary component suffers from an inherent labeling ambiguity. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: 16 pages, 9 figures

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

Black-hole formation in binary neutron star mergers: The impact of spin on the prompt-collapse scenario

Authors: Federico Schianchi, Maximiliano Ujevic, Anna Neuweiler, Henrique Gieg, Ivan Markin, Tim Dietrich

Abstract: Accurate modeling of the multi-messenger signatures connected to binary neutron star mergers requires proper knowledge on the final remnant's fate and the conditions under which black holes (BHs) can form in such mergers. In this article, we use a suite of 84 numerical-relativity simulations in 28 different physical setups to explore the impact of the individual stars' spin on the merger outcome a… ▽ More Accurate modeling of the multi-messenger signatures connected to binary neutron star mergers requires proper knowledge on the final remnant's fate and the conditions under which black holes (BHs) can form in such mergers. In this article, we use a suite of 84 numerical-relativity simulations in 28 different physical setups to explore the impact of the individual stars' spin on the merger outcome and on the early postmerger dynamics. We find that for setups close to the prompt-collapse threshold, the stars' intrinsic spin significantly changes the lifespan of the remnant before collapse and that the mass of the debris disk surrounding the BH is also altered. To enable a better understanding of BH formation, we check if there is at least a theoretical chance of observing densities that are above the maximum density allowed in a stable isolated neutron star, and we investigate the importance of different pressure contributions on the evolution of the postmerger remnant and BH formation. △ Less

Submitted 1 March, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

An overview of existing and new nuclear and astrophysical constraints on the equation of state of neutron-rich dense matter

Authors: Hauke Koehn, Henrik Rose, Peter T. H. Pang, Rahul Somasundaram, Brendan T. Reed, Ingo Tews, Adrian Abac, Oleg Komoltsev, Nina Kunert, Aleksi Kurkela, Michael W. Coughlin, Brian F. Healy, Tim Dietrich

Abstract: Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. Here, we apply these different pieces of data individually to a broad set of physics-agnostic candidate EOSs and analyze the resulting constraints. Specifically, we make use of… ▽ More Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. Here, we apply these different pieces of data individually to a broad set of physics-agnostic candidate EOSs and analyze the resulting constraints. Specifically, we make use of information from chiral effective field theory, perturbative quantum chromodynamics, as well as data from heavy-ion collisions and the PREX-II and CREX experiments. We also investigate the impact of current mass and radius measurements of neutron stars, such as radio timing measurements of heavy pulsars, NICER data, and other X-ray observations. We augment these by reanalyses of the gravitational-wave (GW) signal GW170817, its associated kilonova AT2017gfo and gamma-ray burst afterglow, the GW signal GW190425, and the GRB211211A afterglow, where we use improved models for the tidal waveform and kilonova light curves. Additionally, we consider the postmerger fate of GW170817 and its consequences for the EOS. This large and diverse set of constraints is eventually combined in numerous ways to explore limits on quantities such as the typical neutron-star radius, the maximum neutron-star mass, the nuclear symmetry-energy parameters, and the speed of sound. Based on the priors from our EOS candidate set, we find the radius of the canonical 1.4 M$_\odot$ neutron star to be $R_{1.4}= 12.27_{-0.94}^{+0.83}$ km and the TOV mass $M_{\rm TOV}= 2.26_{-0.22}^{+0.45}$ M$_\odot$ at 95% credibility, when including those constraints where systematic uncertainties are deemed small. A less conservative approach, combining all the presented constraints, similarly yields $R_{1.4}= 12.20_{-0.50}^{+0.53}$ km and $M_{\rm TOV}= 2.31_{-0.20}^{+0.08}$ M$_\odot$. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: 49 pages, 32 figures, webinterface for custom constraint combinations in https://enlil.gw.physik.uni-potsdam.de/eos_constraints

Report number: LA-UR-24-20420

IMRPhenomXP_NRTidalv2: An improved frequency-domain precessing binary neutron star waveform model

Authors: Marta Colleoni, Felip A. Ramis Vidal, Nathan K. Johnson-McDaniel, Tim Dietrich, Maria Haney, Geraint Pratten

Abstract: We present two new frequency-domain gravitational waveform models for the analysis of signals emitted by binary neutron star coalescences: IMRPhenomXAS_NRTidalv2 and IMRPhenomXP_NRTidalv2. Both models are available through the public algorithm library LALSuite and represent the first extensions of IMRPhenomX models including matter effects. We show here that these two models represent a significan… ▽ More We present two new frequency-domain gravitational waveform models for the analysis of signals emitted by binary neutron star coalescences: IMRPhenomXAS_NRTidalv2 and IMRPhenomXP_NRTidalv2. Both models are available through the public algorithm library LALSuite and represent the first extensions of IMRPhenomX models including matter effects. We show here that these two models represent a significant advancement in efficiency and accuracy with respect to their phenomenological predecessors, IMRPhenomD_NRTidalv2 and IMRPhenomPv2_NRTidalv2. The computational efficiency of the new models is achieved through the application of the same multibanding technique previously applied to binary black hole models. Furthermore, IMRPhenomXP_NRTidalv2 implements a more accurate description of the precession dynamics, including double-spin effects and, optionally, matter effects in the twisting-up construction. The latter are available through an option to use a numerical integration of the post-Newtonian precession equations. We show that the new precession descriptions allow the model to better reproduce the phenomenology observed in numerical-relativity simulations of precessing binary neutron stars. Finally, we present some applications of the new models to Bayesian parameter estimation studies, including a reanalysis of GW170817 and a study of simulated observations using numerical relativity waveforms for nonprecessing binary neutron stars with highly spinning components. We find that in these cases the new models make a negligible difference in the results. Nevertheless, by virtue of the aforementioned improvements, the new models represent valuable tools for the study of future detections of coalescing binary neutron stars. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 19 pages, 16 figures

NRTidalv3: A New and Robust Gravitational-Waveform Model for High-Mass-Ratio Binary Neutron Star Systems with Dynamical Tidal Effects

Authors: Adrian Abac, Tim Dietrich, Alessandra Buonanno, Jan Steinhoff, Maximiliano Ujevic

Abstract: For the analysis of gravitational-wave signals, fast and accurate gravitational-waveform models are required. These enable us to obtain information on the system properties from compact binary mergers. In this article, we introduce the NRTidalv3 model, which contains a closed-form expression that describes tidal effects, focusing on the description of binary neutron star systems. The model improve… ▽ More For the analysis of gravitational-wave signals, fast and accurate gravitational-waveform models are required. These enable us to obtain information on the system properties from compact binary mergers. In this article, we introduce the NRTidalv3 model, which contains a closed-form expression that describes tidal effects, focusing on the description of binary neutron star systems. The model improves upon previous versions by employing a larger set of numerical-relativity data for its calibration, by including high-mass ratio systems covering also a wider range of equations of state. It also takes into account dynamical tidal effects and the known post-Newtonian mass-ratio dependence of individual calibration parameters. We implemented the model in the publicly available LALSuite software library by augmenting different binary black hole waveform models (IMRPhenomD, IMRPhenomX, and SEOBNRv5_ROM). We test the validity of NRTidalv3 by comparing it with numerical-relativity waveforms, as well as other tidal models. Finally, we perform parameter estimation for GW170817 and GW190425 with the new tidal approximant and find overall consistent results with respect to previous studies. △ Less

Submitted 13 November, 2023; originally announced November 2023.

What can we learn about the unstable equation-of-state branch from neutron-star mergers?

Authors: Maximiliano Ujevic, Rahul Somasundaram, Tim Dietrich, Jerome Margueron, Ingo Tews

Abstract: The Equation of State (EOS) of dense strongly-interacting matter can be probed by astrophysical observations of neutron stars (NS), such as X-ray detections of pulsars or the measurement of the tidal deformability of NSs during the inspiral stage of NS mergers. These observations constrain the EOS at most up to the density of the maximum-mass configuration, $n_\textrm{TOV}$, which is the highest d… ▽ More The Equation of State (EOS) of dense strongly-interacting matter can be probed by astrophysical observations of neutron stars (NS), such as X-ray detections of pulsars or the measurement of the tidal deformability of NSs during the inspiral stage of NS mergers. These observations constrain the EOS at most up to the density of the maximum-mass configuration, $n_\textrm{TOV}$, which is the highest density that can be explored by stable NSs for a given EOS. However, under the right circumstances, binary neutron star (BNS) mergers can create a postmerger remnant that explores densities above $n_\textrm{TOV}$. In this work, we explore whether the EOS above $n_\textrm{TOV}$ can be measured from gravitational-wave or electromagnetic observations of the postmerger remnant. We perform a total of twenty-five numerical-relativity simulations of BNS mergers for a range of EOSs and find no case in which different descriptions of the matter above $n_{\rm TOV}$ have a detectable impact on postmerger observables. Hence, we conclude that the EOS above $n_\textrm{TOV}$ can likely not be probed through BNS merger observations for the current and next generation of detectors. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 8 pages, 4 figures

Report number: LA-UR-23-31683

Inference via Wild Bootstrap and Multiple Imputation under Fine-Gray Models with Incomplete Data

Authors: Marina T. Dietrich, Dennis Dobler, Mathisca C. M. de Gunst

Abstract: Fine-Gray models specify the subdistribution hazards for one out of multiple competing risks to be proportional. The estimators of parameters and cumulative incidence functions under Fine-Gray models have a simpler structure when data are censoring-complete than when they are more generally incomplete. This paper considers the case of incomplete data but it exploits the above-mentioned simpler est… ▽ More Fine-Gray models specify the subdistribution hazards for one out of multiple competing risks to be proportional. The estimators of parameters and cumulative incidence functions under Fine-Gray models have a simpler structure when data are censoring-complete than when they are more generally incomplete. This paper considers the case of incomplete data but it exploits the above-mentioned simpler estimator structure for which there exists a wild bootstrap approach for inferential purposes. The present idea is to link the methodology under censoring-completeness with the more general right-censoring regime with the help of multiple imputation. In a simulation study, this approach is compared to the estimation procedure proposed in the original paper by Fine and Gray when it is combined with a bootstrap approach. An application to a data set about hospital-acquired infections illustrates the method. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 32 pages, 2 figures, 1 table

MSC Class: 62N01; 62N02; 62J99; 62G09; 62G15

Wild Bootstrap for Counting Process-Based Statistics

Authors: Marina T. Dietrich, Dennis Dobler, Mathisca C. M. de Gunst

Abstract: The wild bootstrap is a popular resampling method in the context of time-to-event data analyses. Previous works established the large sample properties of it for applications to different estimators and test statistics. It can be used to justify the accuracy of inference procedures such as hypothesis tests or time-simultaneous confidence bands. This paper consists of two parts: in Part~I, a genera… ▽ More The wild bootstrap is a popular resampling method in the context of time-to-event data analyses. Previous works established the large sample properties of it for applications to different estimators and test statistics. It can be used to justify the accuracy of inference procedures such as hypothesis tests or time-simultaneous confidence bands. This paper consists of two parts: in Part~I, a general framework is developed in which the large sample properties are established in a unified way by using martingale structures. The framework includes most of the well-known non- and semiparametric statistical methods in time-to-event analysis and parametric approaches. In Part II, the Fine-Gray proportional sub-hazards model exemplifies the theory for inference on cumulative incidence functions given the covariates. The model falls within the framework if the data are censoring-complete. A simulation study demonstrates the reliability of the method and an application to a data set about hospital-acquired infections illustrates the statistical procedure. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 2 parts, 115 pages, 2 figures, 13 tables

MSC Class: 62G09; 62G20; 62N02; 62N03; 62F03; 62F12; 62F25; 62F40

Multi-band analyses of the bright GRB 230812B and the associated SN2023pel

Authors: T. Hussenot-Desenonges, T. Wouters, N. Guessoum, I. Abdi, A. Abulwfa, C. Adami, J. F. Agüí Fernández, T. Ahumada, V. Aivazyan, D. Akl, S. Anand, C. M. Andrade, S. Antier, S. A. Ata, P. D'Avanzo, Y. A. Azzam, A. Baransky, S. Basa, M. Blazek, P. Bendjoya, S. Beradze, P. Boumis, M. Bremer, R. Brivio, V. Buat , et al. (87 additional authors not shown)

Abstract: GRB~230812B is a bright and relatively nearby ($z =0.36$) long gamma-ray burst (GRB) that has generated significant interest in the community and has thus been observed over the entire electromagnetic spectrum. We report over 80 observations in X-ray, ultraviolet, optical, infrared, and sub-millimeter bands from the GRANDMA (Global Rapid Advanced Network for Multi-messenger Addicts) network of obs… ▽ More GRB~230812B is a bright and relatively nearby ($z =0.36$) long gamma-ray burst (GRB) that has generated significant interest in the community and has thus been observed over the entire electromagnetic spectrum. We report over 80 observations in X-ray, ultraviolet, optical, infrared, and sub-millimeter bands from the GRANDMA (Global Rapid Advanced Network for Multi-messenger Addicts) network of observatories and from observational partners. Adding complementary data from the literature, we then derive essential physical parameters associated with the ejecta and external properties (i.e. the geometry and environment) of the GRB and compare with other analyses of this event. We spectroscopically confirm the presence of an associated supernova, SN2023pel, and we derive a photospheric expansion velocity of v $\sim$ 17$\times10^3$ km s$^{-1}$. We analyze the photometric data first using empirical fits of the flux and then with full Bayesian Inference. We again strongly establish the presence of a supernova in the data, with a maximum (pseudo-)bolometric luminosity of $5.75 \times 10^{42}$ erg/s, at $15.76^{+0.81}_{-1.21}$ days (in the observer frame) after the trigger, with a half-max time width of 22.0 days. We compare these values with those of SN1998bw, SN2006aj, and SN2013dx. Our best-fit model favours a very low density environment ($\log_{10}({n_{\rm ISM}/{\rm cm}^{-3}}) = -2.38^{+1.45}_{-1.60}$) and small values for the jet's core angle $θ_{\rm core} = 1.54^{+1.02}_{-0.81} \ \rm{deg}$ and viewing angle $θ_{\rm obs} = 0.76^{+1.29}_{-0.76} \ \rm{deg}$. GRB 230812B is thus one of the best observed afterglows with a distinctive supernova bump. △ Less

Submitted 17 February, 2024; v1 submitted 22 October, 2023; originally announced October 2023.

Comparing gravitational waveform models for binary black hole mergers through a hypermodels approach

Authors: Anna Puecher, Anuradha Samajdar, Gregory Ashton, Chris Van Den Broeck, Tim Dietrich

Abstract: The inference of source parameters from gravitational-wave signals relies on theoretical models that describe the emitted waveform. Different model assumptions on which the computation of these models is based could lead to biases in the analysis of gravitational-wave data. In this work, we sample directly on four state-of-the-art binary black hole waveform models from different families, in order… ▽ More The inference of source parameters from gravitational-wave signals relies on theoretical models that describe the emitted waveform. Different model assumptions on which the computation of these models is based could lead to biases in the analysis of gravitational-wave data. In this work, we sample directly on four state-of-the-art binary black hole waveform models from different families, in order to investigate these systematic biases from the 13 heaviest gravitational-wave sources with moderate to high signal-to-noise ratios in the third Gravitational-Wave Transient Catalog (GWTC- 3). All models include spin-precession as well as higher-order modes. Using the "hypermodels" technique, we treat the waveform models as one of the sampled parameters, therefore directly getting the odds ratio of one waveform model over another from a single parameter estimation run. From the joint odds ratio over all 13 sources, we find the model NRSur7dq4 to be favoured over SEOBNRv4PHM, with an odds ratio of 29.43; IMRPhenomXPHM and IMRPhenomTPHM have an odds ratio, respectively, of 4.70 and 5.09 over SEOBNRv4PHM. However, this result is mainly determined by three events that show a strong preference for some of the models and that are all affected by possible data quality issues. If we do not consider these potentially problematic events, the odds ratio do not exhibit a significant preference for any of the models. Although further work studying a larger set of signals will be needed for robust quantitative results, the presented method highlights one possible avenue for future waveform model development. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Probing Quarkyonic Matter in Neutron Stars with the Bayesian Nuclear-Physics Multi-Messenger Astrophysics Framework

Authors: Peter T. H. Pang, Lars Sivertsen, Rahul Somasundaram, Tim Dietrich, Srimoyee Sen, Ingo Tews, Michael Coughlin, Chris Van Den Broeck

Abstract: The interior of neutron stars contains matter at the highest densities realized in our Universe. Interestingly, theoretical studies of dense matter, in combination with the existence of two solar mass neutron stars, indicate that the speed of sound $c_s$ has to increase to values well above the conformal limit ($c_s^2\sim 1/3$) before decreasing again at higher densities. The decrease could be exp… ▽ More The interior of neutron stars contains matter at the highest densities realized in our Universe. Interestingly, theoretical studies of dense matter, in combination with the existence of two solar mass neutron stars, indicate that the speed of sound $c_s$ has to increase to values well above the conformal limit ($c_s^2\sim 1/3$) before decreasing again at higher densities. The decrease could be explained by either a strong first-order phase transition or a cross-over transition from hadronic to quark matter. The latter scenario leads to a pronounced peak in the speed of sound reaching values above the conformal limit, naturally explaining the inferred behavior. In this work, we use the Nuclear-Physics Multi-Messenger Astrophysics framework \textsc{NMMA} to compare predictions of the quarkyonic matter model with astrophysical observations of neutron stars, with the goal of constraining model parameters. Assuming quarkyonic matter to be realized within neutron stars, we find that there can be a significant amount of quarks inside the core of neutron stars with masses in the two solar mass range, amounting to up to $\sim 0.13M_\odot$, contributing $\sim 5.9\%$ of the total mass. Furthermore, for the quarkyonic matter model investigated here, the radius of a $1.4M_\odot$ neutron star would be $13.44^{+1.69}_{-1.54} (13.54^{+1.02}_{-1.04})$ km, at $95\%$ credibility, without (with) the inclusion of AT2017gfo. △ Less

Submitted 28 February, 2024; v1 submitted 29 August, 2023; originally announced August 2023.

Comments: 14 pages, 10 figures

Report number: LA-UR-23-29118, NP3M-P2300017

Journal ref: Phys. Rev. C 109, 025807 (2024)

A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run

Authors: C. Fletcher, J. Wood, R. Hamburg, P. Veres, C. M. Hui, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, A. Goldstein, B. A. Hristov, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, S. Poolakkil, A. von Kienlin, C. A. Wilson-Hodge, The Fermi Gamma-ray Burst Monitor Team, M. Crnogorčević, J. DeLaunay, A. Tohuvavohu, R. Caputo, S. B. Cenko , et al. (1674 additional authors not shown)

Abstract: We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,… ▽ More We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Nuclear physics constraints from binary neutron star mergers in the Einstein Telescope era

Authors: Francesco Iacovelli, Michele Mancarella, Chiranjib Mondal, Anna Puecher, Tim Dietrich, Francesca Gulminelli, Michele Maggiore, Micaela Oertel

Abstract: The next generation of ground-based gravitational-wave detectors, Einstein Telescope (ET) and Cosmic Explorer (CE), present a unique opportunity to put constraints on dense matter, among many other groundbreaking scientific goals. In a recent study the science case of ET was further strengthened, studying in particular the performances of different detector designs. In this paper we present a more… ▽ More The next generation of ground-based gravitational-wave detectors, Einstein Telescope (ET) and Cosmic Explorer (CE), present a unique opportunity to put constraints on dense matter, among many other groundbreaking scientific goals. In a recent study the science case of ET was further strengthened, studying in particular the performances of different detector designs. In this paper we present a more detailed study of the nuclear physics section of that work. In particular, focusing on two different detector configurations (the single-site triangular-shaped design and a design consisting of two widely separated "L-shaped" interferometers), we study the detection prospects of binary neutron star (BNS) mergers, and how they can reshape our understanding of the underlying equation of state (EoS) of dense matter. We employ several state-of-the-art EoS models and state-of-the-art synthetic BNS merger catalogs, and we make use of the Fisher information formalism (FIM) to quantify statistical errors on the astrophysical parameters describing individual BNS events. To check the reliability of the FIM method, we further perform a full parameter estimation for a few simulated events. Based on the uncertainties on the tidal deformabilities associated to these events, we outline a mechanism to extract the underlying injected EoS using a recently developed meta-modelling approach within a Bayesian framework. Our results suggest that with $\gtrsim 500$ events with signal-to-noise ratio greater than $12$, we will be able to pin down very precisely the underlying EoS governing the neutron star matter. △ Less

Submitted 10 November, 2023; v1 submitted 22 August, 2023; originally announced August 2023.

Comments: 18 pages, 10 figures. v2 matches the version accepted for publication in PRD

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Chemical Distribution of the Dynamical Ejecta in the Neutron Star Merger GW170817

Authors: Shreya Anand, Peter T. H. Pang, Mattia Bulla, Michael W. Coughlin, Tim Dietrich, Brian Healy, Thomas Hussenot-Desenonges, Theophile Jegou du Laz, Mansi M. Kasliwal, Nina Kunert, Ivan Markin, Kunal Mooley, Vsevolod Nedora, Anna Neuweiler

Abstract: GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the i… ▽ More GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the inclusion of precise inclination information can reveal new insights about the ejecta properties, in particular, about the composition of the dynamical ejecta of AT2017gfo. Our analysis relies on updated kilonova modeling, which includes state-of-the-art heating rates, thermalization efficiencies, and opacities and is parameterized by $\bar{Y}_{\rm e,dyn}$, the average electron fraction of the dynamical ejecta component. Using this model, we incorporate the latest inclination angle constraint of AT2017gfo into a light curve fitting framework to derive updated parameter estimates. Our results suggest that the viewing angle of the observer is pointed towards the lanthanide-poor ($Y_{\rm e,dyn}\gtrsim0.25$), squeezed polar dynamical ejecta component, which can explain the early blue emission observed in the light curve of AT2017gfo. In contrast to a recent claim of spherical ejecta powering AT2017gfo, our study indicates that the composition of the dynamical ejecta has a strong angular dependence, with a lanthanide-rich ($Y_{\rm e,dyn}\lesssim0.25$), tidal component distributed around the merger plane with a half-opening angle of $35^\circ$. The inclination angle constraint reduces $\bar{Y}_{\rm e,dyn}$ from $0.24$ to $0.22$, with values $0.17\lesssim Y_{\rm e, dyn} \lesssim0.41$ enabling the robust production of $r$-process elements up to the $3^{\rm rd}$ peak in the tidal dynamical ejecta. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 8 pages, 6 figures

M1 neutrino transport within the numerical-relativistic code BAM with application to low mass binary neutron star mergers

Authors: Federico Schianchi, Henrique Gieg, Vsevolod Nedora, Anna Neuweiler, Maximiliano Ujevic, Mattia Bulla, Tim Dietrich

Abstract: Neutrino interactions are essential for an accurate understanding of the binary neutron star merger process. In this article, we extend the code infrastructure of the well-established numerical-relativity code BAM that until recently neglected neutrino-driven interactions. In fact, while previous work allowed already the usage of nuclear-tabulated equations of state and employing a neutrino leakag… ▽ More Neutrino interactions are essential for an accurate understanding of the binary neutron star merger process. In this article, we extend the code infrastructure of the well-established numerical-relativity code BAM that until recently neglected neutrino-driven interactions. In fact, while previous work allowed already the usage of nuclear-tabulated equations of state and employing a neutrino leakage scheme, we are moving forward by implementing a first-order multipolar radiation transport scheme (M1) for the advection of neutrinos. After testing our implementation on a set of standard scenarios, we apply it to the evolution of four low-mass binary systems, and we perform an analysis of ejecta properties. We also show that our new ejecta analysis infrastructure is able to provide numerical relativity-informed inputs for the codes $\texttt{POSSIS}$ and $\texttt{Skynet}$, for the computation of kilonova lightcurves and nucleosynthesis yields, respectively. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 29 pages, 27 figures

What is the nature of the HESS J1731-347 compact object?

Authors: Violetta Sagun, Edoardo Giangrandi, Tim Dietrich, Oleksii Ivanytskyi, Rodrigo Negreiros, Constança Providência

Abstract: Once further confirmed in future analyses, the radius and mass measurement of HESS J1731-347 with $M=0.77^{+0.20}_{-0.17}~M_{\odot}$ and $R=10.4^{+0.86}_{-0.78}~\rm km$ will be among the lightest and smallest compact objects ever detected. This raises many questions about its nature and opens up the window for different theories to explain such a measurement. In this article, we use the informatio… ▽ More Once further confirmed in future analyses, the radius and mass measurement of HESS J1731-347 with $M=0.77^{+0.20}_{-0.17}~M_{\odot}$ and $R=10.4^{+0.86}_{-0.78}~\rm km$ will be among the lightest and smallest compact objects ever detected. This raises many questions about its nature and opens up the window for different theories to explain such a measurement. In this article, we use the information from Doroshenko et al. (2022) on the mass, radius, and surface temperature together with the multimessenger observations of neutron stars to investigate the possibility that HESS J1731-347 is one of the lightest observed neutron star, a strange quark star, a hybrid star with an early deconfinement phase transition, or a dark matter-admixed neutron star. The nucleonic and quark matter are modeled within realistic equation of states (EOSs) with a self-consistent calculation of the pairing gaps in quark matter. By performing the joint analysis of the thermal evolution and mass-radius constraint, we find evidence that within a 1$σ$ confidence level, HESS J1731-347 is consistent with the neutron star scenario with the soft EOS as well as with a strange and hybrid star with the early deconfinement phase transition with a strong quark pairing and neutron star admixed with dark matter. △ Less

Submitted 12 November, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 9 pages, 3 figures

Journal ref: Astrophys. J. 958, 49 (2023)

Updated observing scenarios and multi-messenger implications for the International Gravitational-wave Network's O4 and O5

Authors: R. Weizmann Kiendrebeogo, Amanda M. Farah, Emily M. Foley, Abigail Gray, Nina Kunert, Anna Puecher, Andrew Toivonen, R. Oliver VandenBerg, Shreya Anand, Tomás Ahumada, Viraj Karambelkar, Michael W. Coughlin, Tim Dietrich, S. Zacharie Kam, Peter T. H. Pang, Leo P. Singer, Niharika Sravan

Abstract: An advanced LIGO and Virgo's third observing run brought another binary neutron star merger (BNS) and the first neutron-star black hole mergers. While no confirmed kilonovae were identified in conjunction with any of these events, continued improvements of analyses surrounding GW170817 allow us to project constraints on the Hubble Constant ($H_0$), the Galactic enrichment from $r$-process nucleosy… ▽ More An advanced LIGO and Virgo's third observing run brought another binary neutron star merger (BNS) and the first neutron-star black hole mergers. While no confirmed kilonovae were identified in conjunction with any of these events, continued improvements of analyses surrounding GW170817 allow us to project constraints on the Hubble Constant ($H_0$), the Galactic enrichment from $r$-process nucleosynthesis, and ultra-dense matter possible from forthcoming events. Here, we describe the expected constraints based on the latest expected event rates from the international gravitational-wave network (IGWN) and analyses of GW170817. We show the expected detection rate of gravitational waves and their counterparts, as well as how sensitive potential constraints are to the observed numbers of counterparts. We intend this analysis as support for the community when creating scientifically driven electromagnetic follow-up proposals. During the next observing run O4, we predict an annual detection rate of electromagnetic counterparts from BNS of $0.43^{+0.58}_{-0.26}$ ($1.97^{+2.68}_{-1.2}$) for the Zwicky Transient Facility (Rubin Observatory). △ Less

Submitted 12 December, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Journal ref: The Astrophysical Journal , 2023, 958, 158

General-Relativistic Hydrodynamics Simulation of a Neutron Star - Sub-Solar-Mass Black Hole Merger

Authors: Ivan Markin, Anna Neuweiler, Adrian Abac, Swami Vivekanandji Chaurasia, Maximiliano Ujevic, Mattia Bulla, Tim Dietrich

Abstract: Over the last few years, there has been an increasing interest in sub-solar mass black holes due to their potential to provide valuable information about cosmology or the black hole population. Motivated by this, we study observable phenomena connected to the merger of a sub-solar mass black hole with a neutron star. For this purpose, we perform new numerical-relativity simulations of a binary sys… ▽ More Over the last few years, there has been an increasing interest in sub-solar mass black holes due to their potential to provide valuable information about cosmology or the black hole population. Motivated by this, we study observable phenomena connected to the merger of a sub-solar mass black hole with a neutron star. For this purpose, we perform new numerical-relativity simulations of a binary system composed of a black hole with mass $0.5M_\odot$ and a neutron star with mass $1.4 M_\odot$. We investigate the merger dynamics of this exotic system and provide information about the connected gravitational-wave and kilonova signals. Our study indicates that current gravitational-waveform models are unable to adequately describe such systems and that phenomenological relations connecting the binary parameters with the ejecta and remnant properties are not applicable to our system. Furthermore, we find a dependence of the kilonova signal on the azimuthal viewing angle due to the asymmetric mass ejection. This first-of-its-kind simulation opens the door for the study of sub-solar mass black hole - neutron star mergers and could serve as a testing ground for future model development. △ Less

Submitted 23 April, 2023; originally announced April 2023.

Comments: 16 pages, 14 figures, to be submitted to PRD, comments welcome

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated… ▽ More Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 28 pages, 11 figures

Report number: LIGO-P2200031

Measuring tidal effects with the Einstein Telescope: A design study

Authors: Anna Puecher, Anuradha Samajdar, Tim Dietrich

Abstract: Over the last few years, there has been a large momentum to ensure that the third-generation era of gravitational wave detectors will find its realisation in the next decades, and numerous design studies have been ongoing for some time. Some of the main factors determining the cost of the Einstein Telescope lie in the length of the interferometer arms and its shape: L-shaped detectors versus a sin… ▽ More Over the last few years, there has been a large momentum to ensure that the third-generation era of gravitational wave detectors will find its realisation in the next decades, and numerous design studies have been ongoing for some time. Some of the main factors determining the cost of the Einstein Telescope lie in the length of the interferometer arms and its shape: L-shaped detectors versus a single triangular configuration. Both designs are further expected to include a xylophone configuration for improvement on both ends of the frequency bandwidth of the detector. We consider binary neutron star sources in our study, as examples of sources already observed with the current generation detectors and ones which hold most promise given the broader frequency band and higher sensitivity of the third-generation detectors. We estimate parameters of the sources, with different kinds of configurations of the Einstein Telescope detector, varying arm-lengths as well as shapes and alignments. Overall, we find little improvement with respect to changing the shape, or alignment. However, there are noticeable differences in the estimates of some parameters, including tidal deformability, when varying the arm-length of the detectors. In addition, we also study the effect of changing the laser power, and the lower limit of the frequency band in which we perform the analysis. △ Less

Submitted 9 April, 2023; originally announced April 2023.

Comments: 11 pages, 7 figures, 4 tables

Science with the Einstein Telescope: a comparison of different designs

Authors: Marica Branchesi, Michele Maggiore, David Alonso, Charles Badger, Biswajit Banerjee, Freija Beirnaert, Enis Belgacem, Swetha Bhagwat, Guillaume Boileau, Ssohrab Borhanian, Daniel David Brown, Man Leong Chan, Giulia Cusin, Stefan L. Danilishin, Jerome Degallaix, Valerio De Luca, Arnab Dhani, Tim Dietrich, Ulyana Dupletsa, Stefano Foffa, Gabriele Franciolini, Andreas Freise, Gianluca Gemme, Boris Goncharov, Archisman Ghosh , et al. (51 additional authors not shown)

Abstract: The Einstein Telescope (ET), the European project for a third-generation gravitational-wave detector, has a reference configuration based on a triangular shape consisting of three nested detectors with 10 km arms, where in each arm there is a `xylophone' configuration made of an interferometer tuned toward high frequencies, and an interferometer tuned toward low frequencies and working at cryogeni… ▽ More The Einstein Telescope (ET), the European project for a third-generation gravitational-wave detector, has a reference configuration based on a triangular shape consisting of three nested detectors with 10 km arms, where in each arm there is a `xylophone' configuration made of an interferometer tuned toward high frequencies, and an interferometer tuned toward low frequencies and working at cryogenic temperature. Here, we examine the scientific perspectives under possible variations of this reference design. We perform a detailed evaluation of the science case for a single triangular geometry observatory, and we compare it with the results obtained for a network of two L-shaped detectors (either parallel or misaligned) located in Europe, considering different choices of arm-length for both the triangle and the 2L geometries. We also study how the science output changes in the absence of the low-frequency instrument, both for the triangle and the 2L configurations. We examine a broad class of simple `metrics' that quantify the science output, related to compact binary coalescences, multi-messenger astronomy and stochastic backgrounds, and we then examine the impact of different detector designs on a more specific set of scientific objectives. △ Less

Submitted 17 June, 2023; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 197 pages, 73 figures. v2: corrections in the part on the sensitivity to stochastic backgrounds. Accepted in JCAP

Report number: ET-0084A-23

Impact of jets on kilonova photometric and polarimetric emission from binary neutron star mergers

Authors: Manisha Shrestha, Mattia Bulla, Lorenzo Nativi, Ivan Markin, Stephan Rosswog, Tim Dietrich

Abstract: A merger of binary neutron stars creates heavy unstable elements whose radioactive decay produces a thermal emission known as a kilonova. In this paper, we predict the photometric and polarimetric behaviour of this emission by performing 3-D Monte Carlo radiative transfer simulations. In particular, we choose three hydrodynamical models for merger ejecta, two including jets with different luminosi… ▽ More A merger of binary neutron stars creates heavy unstable elements whose radioactive decay produces a thermal emission known as a kilonova. In this paper, we predict the photometric and polarimetric behaviour of this emission by performing 3-D Monte Carlo radiative transfer simulations. In particular, we choose three hydrodynamical models for merger ejecta, two including jets with different luminosities and one without a jet structure, to help decipher the impact of jets on the light curve and polarimetric behaviour. In terms of photometry, we find distinct color evolutions across the three models. Models without a jet show the highest variation in light curves for different viewing angles. In contrast, to previous studies, we find models with a jet to produce fainter kilonovae when viewed from orientations close to the jet axis, compared to a model without a jet. In terms of polarimetry, we predict relatively low levels (<~0.3-0.4%) at all orientations that, however, remain non-negligible until a few days after the merger and longer than previously found. Despite the low levels, we find that the presence of a jet enhances the degree of polarization at wavelengths ranging from 0.25 to 2.5\micron, an effect that is found to increase with the jet luminosity. Thus, future photometric and polarimetric campaigns should observe kilonovae in blue and red filters for a few days after the merger to help constrain the properties of the ejecta (e.g. composition) and jet. △ Less

Submitted 24 May, 2023; v1 submitted 24 March, 2023; originally announced March 2023.

Comments: 11 pages, 7 figures, accepted for publication and in press

Revealing the strength of three-nucleon interactions with the Einstein Telescope

Authors: Henrik Rose, Nina Kunert, Tim Dietrich, Peter T. H. Pang, Rory Smith, Chris Van Den Broeck, Stefano Gandolfi, Ingo Tews

Abstract: Three-nucleon forces are crucial for the accurate description of nuclear systems, including dense matter probed in neutron stars. We explore nuclear Hamiltonians that reproduce two-nucleon scattering data and properties of light nuclei, but differ in the three-nucleon interactions among neutrons. While no significantly improved constraints can be obtained from current astrophysical data, we show t… ▽ More Three-nucleon forces are crucial for the accurate description of nuclear systems, including dense matter probed in neutron stars. We explore nuclear Hamiltonians that reproduce two-nucleon scattering data and properties of light nuclei, but differ in the three-nucleon interactions among neutrons. While no significantly improved constraints can be obtained from current astrophysical data, we show that observations of neutron star mergers by next-generation detectors like the proposed Einstein Telescope could provide strong evidence to distinguish between these Hamiltonians. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Report number: LA-UR-22-21458

Synthetic radio images of structured GRB and kilonova afterglows

Authors: Vsevolod Nedora, Tim Dietrich, Masaru Shibata

Abstract: In this paper, we compute and analyze synthetic radio images of gamma-ray bursts and kilonova afterglows. For modeling the former, we consider GRB170817A-inspired set of parameters, while for the latter, we employ ejecta profiles from numerical-relativity simulations. We find that the kilonova afterglow sky map has a doughnut-like structure at early times that becomes more ring-like at late times.… ▽ More In this paper, we compute and analyze synthetic radio images of gamma-ray bursts and kilonova afterglows. For modeling the former, we consider GRB170817A-inspired set of parameters, while for the latter, we employ ejecta profiles from numerical-relativity simulations. We find that the kilonova afterglow sky map has a doughnut-like structure at early times that becomes more ring-like at late times. This is caused by the fact that the synchrotron emission from electrons following Maxwellian distribution function dominates the early, beamed, emission while emissions from electrons following power-law distribution is important at late times. For an on-axis observer, the image flux centroid moves on the image plane initially away from the observer. Flux centroid displacement The image sizes, we find, are the largest for equal mass merger simulations with the soft equation of state. The presence of a kilonova afterglow affects the properties inferred from the source sky map even if the gamma-ray burst afterglow dominates the total flux density. The main effect is the reduction of the mean apparent velocity of the source, and an increase in the source size. Thus, neglecting the presence of the kilonova afterglow may lead to systematic errors in the inference of gamma-ray burst properties from the sky map observations. Notably, at the observing angle inferred for GRB170817A the presence of kilonova afterglow would affect the sky map properties only at very late times $t\gtrsim1500\,$days. △ Less

Submitted 24 February, 2023; originally announced February 2023.

GRANDMA and HXMT Observations of GRB 221009A -- the Standard-Luminosity Afterglow of a Hyper-Luminous Gamma-Ray Burst

Authors: D. A. Kann, S. Agayeva, V. Aivazyan, S. Alishov, C. M. Andrade, S. Antier, A. Baransky, P. Bendjoya, Z. Benkhaldoun, S. Beradze, D. Berezin, M. Boër, E. Broens, S. Brunier, M. Bulla, O. Burkhonov, E. Burns, Y. Chen, Y. P. Chen, M. Conti, M. W. Coughlin, W. W. Cui, F. Daigne, B. Delaveau, H. A. R. Devillepoix , et al. (91 additional authors not shown)

Abstract: GRB 221009A is the brightest Gamma-Ray Burst (GRB) detected in more than 50 years of study. In this paper, we present observations in the X-ray and optical domains after the GRB obtained by the GRANDMA Collaboration (which includes observations from more than 30 professional and amateur telescopes) and the Insight-HXMT Collaboration. We study the optical afterglow with empirical fitting from GRAND… ▽ More GRB 221009A is the brightest Gamma-Ray Burst (GRB) detected in more than 50 years of study. In this paper, we present observations in the X-ray and optical domains after the GRB obtained by the GRANDMA Collaboration (which includes observations from more than 30 professional and amateur telescopes) and the Insight-HXMT Collaboration. We study the optical afterglow with empirical fitting from GRANDMA+HXMT data, augmented with data from the literature up to 60 days. We then model numerically, using a Bayesian approach, the GRANDMA and HXMT-LE afterglow observations, that we augment with Swift-XRT and additional optical/NIR observations reported in the literature. We find that the GRB afterglow, extinguished by a large dust column, is most likely behind a combination of a large Milky-Way dust column combined with moderate low-metallicity dust in the host galaxy. Using the GRANDMA+HXMT-LE+XRT dataset, we find that the simplest model, where the observed afterglow is produced by synchrotron radiation at the forward external shock during the deceleration of a top-hat relativistic jet by a uniform medium, fits the multi-wavelength observations only moderately well, with a tension between the observed temporal and spectral evolution. This tension is confirmed when using the extended dataset. We find that the consideration of a jet structure (Gaussian or power-law), the inclusion of synchrotron self-Compton emission, or the presence of an underlying supernova do not improve the predictions, showing that the modelling of GRB22109A will require going beyond the most standard GRB afterglow model. Placed in the global context of GRB optical afterglows, we find the afterglow of GRB 221009A is luminous but not extraordinarily so, highlighting that some aspects of this GRB do not deviate from the global known sample despite its extreme energetics and the peculiar afterglow evolution. △ Less

Submitted 27 March, 2023; v1 submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted to ApJL for the special issue, 37 pages, 23 pages main text, 6 tables, 13 figures

Open data from the third observing run of LIGO, Virgo, KAGRA and GEO

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Al-Jodah, C. Alléné, A. Allocca , et al. (1719 additional authors not shown)

Abstract: The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in April of 2019 and lasting six months, O3b starting in November of 2019 and lasting five months, and O3GK starting in April of 2020 and lasti… ▽ More The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in April of 2019 and lasting six months, O3b starting in November of 2019 and lasting five months, and O3GK starting in April of 2020 and lasting 2 weeks. In this paper we describe these data and various other science products that can be freely accessed through the Gravitational Wave Open Science Center at https://gwosc.org. The main dataset, consisting of the gravitational-wave strain time series that contains the astrophysical signals, is released together with supporting data useful for their analysis and documentation, tutorials, as well as analysis software packages. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 27 pages, 3 figures

Report number: LIGO-P2200316

Quasi-equilibrium configurations of binary systems of dark matter admixed neutron stars

Authors: Hannes R. Rüter, Violetta Sagun, Wolfgang Tichy, Tim Dietrich

Abstract: Using an adapted version of the SGRID code, we construct for the first time consistent quasi-equilibrium configurations for a binary system consisting of two neutron stars in which each is admixed with dark matter. The stars are modelled as a system of two non-interacting fluids minimally coupled to gravity. For the fluid representing baryonic matter the SLy equation of state is used, whereas the… ▽ More Using an adapted version of the SGRID code, we construct for the first time consistent quasi-equilibrium configurations for a binary system consisting of two neutron stars in which each is admixed with dark matter. The stars are modelled as a system of two non-interacting fluids minimally coupled to gravity. For the fluid representing baryonic matter the SLy equation of state is used, whereas the second fluid, which corresponds to dark matter, is described using the equation of state of a degenerate Fermi gas. We consider two different scenarios for the distribution of the dark matter. In the first scenario the dark matter is confined to the core of the star, whereas in the second scenario the dark matter extends beyond the surface of the baryonic matter, forming a halo around the baryonic star. The presence of dark matter alters the star's reaction to the companion's tidal forces, which we investigate in terms of the coordinate deformation and mass shedding parameters. The constructed quasi-equilibrium configurations mark the first step towards consistent numerical-relativity simulations of dark matter admixed neutron star binaries. △ Less

Submitted 20 September, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: 13 pages, 11 figures, v2: Added investigation of dependence on dark matter particle mass

Bayesian model selection for GRB 211211A through multi-wavelength analyses

Authors: Nina Kunert, Sarah Antier, Vsevolod Nedora, Mattia Bulla, Peter T. H. Pang, Shreya Anand, Michael Coughlin, Ingo Tews, Jennifer Barnes, Thomas Hussenot-Desenonges, Brian Healy, Theophile Jegou du Laz, Meili Pilloix, Weizmann Kiendrebeogo, Tim Dietrich

Abstract: Although GRB 211211A is one of the closest gamma-ray bursts (GRBs), its classification is challenging because of its partially inconclusive electromagnetic signatures. In this paper, we investigate four different astrophysical scenarios as possible progenitors for GRB~211211A: a binary neutron-star merger, a black-hole--neutron-star merger, a core-collapse supernova, and an r-process enriched core… ▽ More Although GRB 211211A is one of the closest gamma-ray bursts (GRBs), its classification is challenging because of its partially inconclusive electromagnetic signatures. In this paper, we investigate four different astrophysical scenarios as possible progenitors for GRB~211211A: a binary neutron-star merger, a black-hole--neutron-star merger, a core-collapse supernova, and an r-process enriched core collapse of a rapidly rotating massive star (a collapsar). We perform a large set of Bayesian multi-wavelength analyses based on different models describing these scenarios and priors to investigate which astrophysical scenarios and processes might be related to GRB~211211A. Our analysis supports previous studies in which the presence of an additional component, likely related to $r$-process nucleosynthesis, is required to explain the observed light curves of GRB~211211A, as it can not solely be explained as a GRB afterglow. Fixing the distance to about $350~\rm Mpc$, namely the distance of the possible host galaxy SDSS J140910.47+275320.8, we find a statistical preference for a binary neutron-star merger scenario. △ Less

Submitted 15 December, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 14 pages, 6 figures

Report number: LA-UR-22-32364

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 2, January 2024, Pages 3900-3911

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO's and Advanced Virgo's third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1680 additional authors not shown)

Abstract: We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate t… ▽ More We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO's and Advanced Virgo's third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim 0.6$ (at 90% confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions we are unable to rule out $f_\mathrm{PBH} = 1$. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound $f_{\mathrm{DBH}} < 10^{-5}$ on the fraction of atomic dark matter collapsed into black holes. △ Less

Submitted 26 January, 2024; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: https://dcc.ligo.org/P2200139

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

Authors: Violetta Sagun, Edoardo Giangrandi, Oleksii Ivanytskyi, Costança Providência, Tim Dietrich

Abstract: 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… ▽ More 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. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 8 pages, 3 figures, Proceeding of the 15th Quark Confinement and the Hadron Spectrum Conference

Reverse phase transitions in binary neutron-star systems with exotic-matter cores

Authors: Maximiliano Ujevic, Henrique Gieg, Federico Schianchi, Swami Vivekanandji Chaurasia, Ingo Tews, Tim Dietrich

Abstract: Multi-messenger observations of binary neutron star mergers provide a unique opportunity to constrain the dense-matter equation of state. Although it is known from quantum chromodynamics that hadronic matter will undergo a phase transition to exotic forms of matter, e.g., quark matter, the onset density of such a phase transition cannot be computed from first principles. Hence, it remains an open… ▽ More Multi-messenger observations of binary neutron star mergers provide a unique opportunity to constrain the dense-matter equation of state. Although it is known from quantum chromodynamics that hadronic matter will undergo a phase transition to exotic forms of matter, e.g., quark matter, the onset density of such a phase transition cannot be computed from first principles. Hence, it remains an open question if such phase transitions occur inside isolated neutron stars or during binary neutron star mergers, or if they appear at even higher densities that are not realized in the Cosmos. In this article, we perform numerical-relativity simulations of neutron-star mergers and investigate scenarios in which the onset density of such a phase transition is exceeded in at least one inspiralling binary component. Our simulations reveal that shortly before the merger it is possible that such stars undergo a "reverse phase transition", i.e., densities decrease and the quark core inside the star disappears leaving a purely hadronic star at merger. After the merger, when densities increase once more, the phase transition occurs again and leads, in the cases considered in this work, to a rapid formation of a black hole. We compute the gravitational-wave signal and the mass ejection for our simulations of such scenarios and find clear signatures that are related to the post-merger phase transition, e.g., smaller ejecta masses due to the softening of the equation of state through the quark core formation. Unfortunately, we do not find measurable imprints of the reverse phase transition. △ Less

Submitted 5 January, 2023; v1 submitted 8 November, 2022; originally announced November 2022.

Comments: 12 pages, 9 figures

Report number: Report-no: LA-UR-22-31740

Parameter estimation methods for analyzing overlapping gravitational wave signals in the third-generation detector era

Authors: Justin Janquart, Tomasz Baka, Anuradha Samajdar, Tim Dietrich, Chris Van Den Broeck

Abstract: In the coming years, third-generation detectors such as Einstein Telescope and Cosmic Explorer will enter the network of ground-based gravitational-wave detectors. Their current design predicts a significantly improved sensitivity band with a lower minimum frequency than existing detectors. This, combined with the increased arm length, leads to two major effects: the detection of more signals and… ▽ More In the coming years, third-generation detectors such as Einstein Telescope and Cosmic Explorer will enter the network of ground-based gravitational-wave detectors. Their current design predicts a significantly improved sensitivity band with a lower minimum frequency than existing detectors. This, combined with the increased arm length, leads to two major effects: the detection of more signals and the detection of longer signals. Both will result in a large number of overlapping signals. It has been shown that such overlapping signals can lead to biases in the recovered parameters, which would adversely affect the science extracted from the observed binary merger signals. In this work, we analyze overlapping binary black hole coalescences with two methods to analyze multi-signal observations: \textit{hierarchical subtraction} and \textit{joint parameter estimation}. We find that these methods enable a reliable parameter extraction in most cases and that joint parameter estimation is usually more precise but comes with higher computational costs. △ Less

Submitted 7 March, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 14 pages, 17 figures

Second release of the CoRe database of binary neutron star merger waveforms

Authors: Alejandra Gonzalez, Francesco Zappa, Matteo Breschi, Sebastiano Bernuzzi, David Radice, Ananya Adhikari, Alessandro Camilletti, Swami Vivekanandji Chaurasia, Georgios Doulis, Surendra Padamata, Alireza Rashti, Maximiliano Ujevic, Bernd Brügmann, William Cook, Tim Dietrich, Albino Perego, Amit Poudel, Wolfgang Tichy

Abstract: We present the second data release of gravitational waveforms from binary neutron star merger simulations performed by the Computational Relativity (CoRe) collaboration. The current database consists of 254 different binary neutron star configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the… ▽ More We present the second data release of gravitational waveforms from binary neutron star merger simulations performed by the Computational Relativity (CoRe) collaboration. The current database consists of 254 different binary neutron star configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to $\ell=m=4$. They span a significant portion of the mass, mass-ratio,spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum. △ Less

Submitted 28 March, 2023; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: 33 pages, 14 figures

Journal ref: Class. Quantum Grav. 40 085011 (2023)

Search for gravitational-wave transients associated with magnetar bursts in Advanced LIGO and Advanced Virgo data from the third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1645 additional authors not shown)

Abstract: Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bu… ▽ More Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bursts come from two magnetars, SGR 1935$+$2154 and Swift J1818.0$-$1607. We also include three other electromagnetic burst events detected by Fermi GBM which were identified as likely coming from one or more magnetars, but they have no association with a known magnetar. No magnetar giant flares were detected during the analysis period. We find no evidence of gravitational waves associated with any of these 16 bursts. We place upper bounds on the root-sum-square of the integrated gravitational-wave strain that reach $2.2 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at 100 Hz for the short-duration search and $8.7 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at $450$ Hz for the long-duration search, given a detection efficiency of 50%. For a ringdown signal at 1590 Hz targeted by the short-duration search the limit is set to $1.8 \times 10^{-22}$ $/\sqrt{\text{Hz}}$. Using the estimated distance to each magnetar, we derive upper bounds on the emitted gravitational-wave energy of $3.2 \times 10^{43}$ erg ($7.3 \times 10^{43}$ erg) for SGR 1935$+$2154 and $8.2 \times 10^{42}$ erg ($2.8 \times 10^{43}$ erg) for Swift J1818.0$-$1607, for the short-duration (long-duration) search. Assuming isotropic emission of electromagnetic radiation of the burst fluences, we constrain the ratio of gravitational-wave energy to electromagnetic energy for bursts from SGR 1935$+$2154 with available fluence information. The lowest of these ratios is $3 \times 10^3$. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 30 pages with appendices, 5 figures, 10 tables

Report number: LIGO-P2100387

Unraveling information about supranuclear-dense matter from the complete binary neutron star coalescence process using future gravitational-wave detector networks

Authors: Anna Puecher, Tim Dietrich, Ka Wa Tsang, Chinmay Kalaghatgi, Soumen Roy, Yoshinta Setyawati, Chris Van Den Broeck

Abstract: Gravitational waves provide us with an extraordinary tool to study the matter inside neutron stars. In particular, the postmerger signal probes an extreme temperature and density regime and will help reveal information about the equation of state of supranuclear-dense matter. Although current detectors are most sensitive to the signal emitted by binary neutron stars before the merger, the upgrades… ▽ More Gravitational waves provide us with an extraordinary tool to study the matter inside neutron stars. In particular, the postmerger signal probes an extreme temperature and density regime and will help reveal information about the equation of state of supranuclear-dense matter. Although current detectors are most sensitive to the signal emitted by binary neutron stars before the merger, the upgrades of existing detectors and the construction of the next generation of detectors will make postmerger detections feasible. For this purpose, we present a new analytical, frequency-domain model for the inspiral-merger-postmerger signal emitted by binary neutron stars systems. The inspiral and merger part of the signals are modeled with IMRPhenomD_NRTidalv2, and we describe the main emission peak of postmerger with a three-parameter Lorentzian, using two different approaches: one in which the Lorentzian parameters are kept free, and one in which we model them via quasi-universal relations. We test the performance of our new complete waveform model in parameter estimation analyses, studying simulated signals obtained from both our developed model and by injecting numerical relativity waveforms. We investigate the performance of different detector networks to determine the improvement that future detectors will bring to our analysis. We consider Advanced LIGO+ and Advanced Virgo+, KAGRA, and LIGO-India. We also study the possible impact of a detector with high sensitivity in the kilohertz band like NEMO, and finally we compare these results to the ones we obtain with third-generation detectors, the Einstein Telescope and the Cosmic Explorer. △ Less

Submitted 20 June, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: Published version

The effects of self-interacting bosonic dark matter on neutron star properties

Authors: Edoardo Giangrandi, Violetta Sagun, Oleksii Ivanytskyi, Constança Providência, Tim Dietrich

Abstract: We propose a model of asymmetric bosonic dark matter (DM) with self-repulsion mediated by the vector field coupled to the complex scalar particles. By adopting the two-fluid formalism, we study different DM distribution regimes, either, fully condensed inside the core of a star or, otherwise, distributed in a dilute halo around a neutron star (NS). We show that DM condensed in a core leads to a de… ▽ More We propose a model of asymmetric bosonic dark matter (DM) with self-repulsion mediated by the vector field coupled to the complex scalar particles. By adopting the two-fluid formalism, we study different DM distribution regimes, either, fully condensed inside the core of a star or, otherwise, distributed in a dilute halo around a neutron star (NS). We show that DM condensed in a core leads to a decrease of the total gravitational mass, radius and tidal deformability compared to a pure baryonic star with the same central density, which we will perceive as an effective softening of the equation of state (EoS). On the other hand, the presence of a DM halo increases the tidal deformability and total gravitational mass. As a result, an accumulated DM inside compact stars could mimic an apparent stiffening of strongly interacting matter equation of state and constraints we impose on it at high densities. From the performed analysis of the effect of DM particles in a MeV-GeV mass-scale, interaction strength, and relative DM fractions inside NSs we obtained a rigorous constraint on model parameters. Finally, we discuss several smoking guns of the presence of DM that are free from the above mentioned apparent modification of the strongly interacting matter equation of state. With this we could be probed with the future astrophysical and gravitational wave (GW) surveys. △ Less

Submitted 14 August, 2023; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 14 pages, 5 figures

Journal ref: ApJ 953 1, 115 (2023)

Snowmass 2021 Cosmic Frontier White Paper: The Dense Matter Equation of State and QCD Phase Transitions

Authors: Slavko Bogdanov, Emmanuel Fonseca, Rahul Kashyap, Aleksi Kurkela, James M. Lattimer, Jocelyn S. Read, Bangalore S. Sathyaprakash, H. Thankful Cromartie, Tim Dietrich, Arnab Dhani, Timothy Dolch, Tyler Gorda, Sebastien Guillot, Wynn C. G. Ho, Rachael Huxford, Frederick K. Lamb, Philippe Landry, Bradley W. Meyers, M. Coleman Miller, Joonas Nättilä, Risto Paatelainen, Chanda Prescod-Weinstein, Saga Säppi, Ingrid H. Stairs, Nikolaos Stergioulas , et al. (4 additional authors not shown)

Abstract: Our limited understanding of the physical properties of matter at ultra-high density, high proton/neutron number asymmetry, and low temperature is presently one of the major outstanding problems in physics. As matter in this extreme state is known to only exist stably in the cores of neutron stars (NSs), complementary measurements from electromagnetic and gravitational wave astrophysical observati… ▽ More Our limited understanding of the physical properties of matter at ultra-high density, high proton/neutron number asymmetry, and low temperature is presently one of the major outstanding problems in physics. As matter in this extreme state is known to only exist stably in the cores of neutron stars (NSs), complementary measurements from electromagnetic and gravitational wave astrophysical observations of NSs, combined with terrestrial laboratory constraints and further theoretical investigations, hold the promise to provide important insight into the properties of matter in a region of the quantum chromodynamics phase space that is otherwise inaccessible. This multidisciplinary endeavor imposes the following requirements for facilities and resources in the upcoming decade and beyond: * A next generation of gravitational wave detectors to uncover more double NS and neutron star-black hole mergers; * Sensitive radio telescopes to find the most massive and fastest spinning NSs; * Large-area, high-time-resolution and/or high angular resolution X-ray telescopes to constrain the NS mass-radius relation; * Suitable laboratory facilities for nuclear physics experiments to constrain the dense matter equation of state; * Funding resources for theoretical studies of matter in this regime; * The availability of modern large-scale high performance computing infrastructure. The same facilities and resources would also enable significant advances in other high-profile fields of inquiry in modern physics such as the nature of dark matter, alternative theories of gravity, nucleon superfluidity and superconductivity, as well as an array of astrophysics, including but not limited to stellar evolution, nucleosynthesis, and primordial black holes. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021) under Cosmic Frontier (CF07: Cosmic probes of fundamental physics); 30 pages, 8 figures

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to bala… ▽ More We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more. △ Less

Submitted 2 January, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, Open Access Journal PDF

Report number: LIGO-P2100110-v13

Journal ref: The Astrophysical Journal Letters, 941, L30 (2022)

Long-Term Simulations of Dynamical Ejecta: Homologous Expansion and Kilonova Properties

Authors: Anna Neuweiler, Tim Dietrich, Mattia Bulla, Swami Vivekanandji Chaurasia, Stephan Rosswog, Maximiliano Ujevic

Abstract: Accurate numerical-relativity simulations are essential to study the rich phenomenology of binary neutron star systems. In this work, we focus on the material that is dynamically ejected during the merger process and on the kilonova transient it produces. Typically, radiative transfer simulations of kilonova light curves from ejecta make the assumption of homologous expansion, but this condition m… ▽ More Accurate numerical-relativity simulations are essential to study the rich phenomenology of binary neutron star systems. In this work, we focus on the material that is dynamically ejected during the merger process and on the kilonova transient it produces. Typically, radiative transfer simulations of kilonova light curves from ejecta make the assumption of homologous expansion, but this condition might not always be met at the end of usually very short numerical-relativity simulations. In this article, we adjust the infrastructure of the BAM code to enable longer simulations of the dynamical ejecta with the aim of investigating when the condition of homologous expansion is satisfied. In fact, we observe that the deviations from a perfect homologous expansion are about 30% at roughly 100ms after the merger. To determine how these deviations might affect the calculation of kilonova light curves, we extract the ejecta data for different reference times and use them as input for radiative transfer simulations. Our results show that the light curves for extraction times later than 80ms after the merger deviate by less than 0.4mag and are mostly consistent with numerical noise. Accordingly, deviations from the homologous expansion for the dynamical ejecta component are negligible for the purpose of kilonova modelling. △ Less

Submitted 31 August, 2022; v1 submitted 29 August, 2022; originally announced August 2022.

Modeling kilonova afterglows: Effects of the thermal electron population and interaction with GRB outflows

Authors: Vsevolod Nedora, Tim Dietrich, Masaru Shibata, Martin Pohl, Ludovica Crosato Menegazzi

Abstract: Given an increasing number of gamma-ray bursts accompanied by potential kilonovae there is a growing importance to advance modelling of kilonova afterglows. In this work, we investigate how the presence of two electron populations that follow a Maxwellian (thermal) and a power-law (non-thermal) distributions affect kilonova afterglow light curves. We employ semi-analytic afterglow model,… ▽ More Given an increasing number of gamma-ray bursts accompanied by potential kilonovae there is a growing importance to advance modelling of kilonova afterglows. In this work, we investigate how the presence of two electron populations that follow a Maxwellian (thermal) and a power-law (non-thermal) distributions affect kilonova afterglow light curves. We employ semi-analytic afterglow model, $\texttt{PyBlastAfterglow}$. We consider kilonova ejecta profiles from ab-initio numerical relativity binary neutron star merger simulations, targeted to GW170817. We do not perform model selection. We find that the emission from thermal electrons dominates at early times. If the interstellar medium density is high (${\simeq}0.1\,\ccm$) it adds an early time peak to the light curve. As ejecta decelerates the spectral and temporal indexes change in a characteristic way that, if observed, can be used to reconstruct the ejecta velocity distribution. For the low interstellar medium density, inferred for GRB 170817A, the emission from the non-thermal electron population generally dominates. We also assess how kilonova afterglow light curves change if the interstellar medium has been partially removed and pre-accelerated by laterally expanding gamma-ray burst ejecta. We find that the main effect is the emission suppression at early time ${\lesssim}10^{3}\,$days, and at its maximum it reaches ${\sim}40\%$ when the fast tail of the kilonova ejecta moves subsonically through the wake of laterally spreading gamma-ray burst ejecta. The subsequent rebrightening, when these ejecta break through and shocks form, is very mild (${\lesssim}10\%$), and may not be observable. △ Less

Submitted 13 January, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 20 pages, 14 figures

The GRANDMA network in preparation for the fourth gravitational-wave observing run

Authors: S. Agayeva, V. Aivazyan, S. Alishov, M. Almualla, C. Andrade, S. Antier, J. -M. Bai, A. Baransky, S. Basa, P. Bendjoya, Z. Benkhaldoun, S. Beradze, D. Berezin, U. Bhardwaj, M. Blazek, O. Burkhonov, E. Burns, S. Caudill, N. Christensen, F. Colas, A. Coleiro, W. Corradi, M. W. Coughlin, T. Culino, D. Darson , et al. (76 additional authors not shown)

Abstract: GRANDMA is a world-wide collaboration with the primary scientific goal of studying gravitational-wave sources, discovering their electromagnetic counterparts and characterizing their emission. GRANDMA involves astronomers, astrophysicists, gravitational-wave physicists, and theorists. GRANDMA is now a truly global network of telescopes, with (so far) 30 telescopes in both hemispheres. It incorpora… ▽ More GRANDMA is a world-wide collaboration with the primary scientific goal of studying gravitational-wave sources, discovering their electromagnetic counterparts and characterizing their emission. GRANDMA involves astronomers, astrophysicists, gravitational-wave physicists, and theorists. GRANDMA is now a truly global network of telescopes, with (so far) 30 telescopes in both hemispheres. It incorporates a citizen science programme (Kilonova-Catcher) which constitutes an opportunity to spread the interest in time-domain astronomy. The telescope network is an heterogeneous set of already-existing observing facilities that operate coordinated as a single observatory. Within the network there are wide-field imagers that can observe large areas of the sky to search for optical counterparts, narrow-field instruments that do targeted searches within a predefined list of host-galaxy candidates, and larger telescopes that are devoted to characterization and follow-up of the identified counterparts. Here we present an overview of GRANDMA after the third observing run of the LIGO/VIRGO gravitational-wave observatories in $2019-2020$ and its ongoing preparation for the forthcoming fourth observational campaign (O4). Additionally, we review the potential of GRANDMA for the discovery and follow-up of other types of astronomical transients. △ Less

Submitted 27 July, 2022; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: Submitted to the Proceedings of the SPIE, Astronomical Telescopes and Instrumentation 2022

Numerical Simulations of Dark Matter Admixed Neutron Star Binaries

Authors: Mattia Emma, Federico Schianchi, Francesco Pannarale, Violetta Sagun, Tim Dietrich

Abstract: 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 f… ▽ More 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. △ Less

Submitted 22 June, 2022; originally announced June 2022.

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

Authors: Henrique Gieg, Federico Schianchi, Tim Dietrich, Maximiliano Ujevic

Abstract: 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… ▽ More 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. △ Less

Submitted 10 July, 2022; v1 submitted 2 June, 2022; originally announced June 2022.

Comments: 28 pages, 14 figures

Journal ref: Universe 2022, 8(7), 370