Cosmology using advanced gravitational-wave detectors alone

The advent of gravitational wave (GW) astronomy has provided a new window through which to view and understand the Universe. To fully exploit the potential of GW astronomy, an understanding of all the potential electromagnetic counterparts to a gravitational wave detected source will help maximise the science returns. Here I present a study of the electromagnetic emission from relativistic jets that accompany the merger of binary neutron stars or black hole-neutron star systems. These counterparts provide a probe for the structure and dynamics of these relativistic outflows. Binary neutron star, or neutron star-black hole, mergers are thought to be the dominant progenitor of short gamma-ray bursts (GRBs). Here we investigate the possibility that there is a hidden population of low-Lorentz factor jets resulting in failed GRBs, on-axis orphan afterglows, and what kind of counterparts can be expected given a merger-jet population dominated by these failed-GRB jets. I find that for GW detected mergers, ∼ 80% of the population of on-axis events may result in a failed GRB afterglow. The afterglow of a failed GRB is characterised by the lack of any prompt emission; where the γ-rays are emitted within an optically thick region of the low-Lorentz factor (Γ) outflow and significant suppression via pair production and a high opacity results in the photons coupled to the pair plasma. This plasma will undergo adiabatic expansion, and the photons will decouple at the photospheric radius. The energy in the prompt photons, for a sufficiently low-Γ outflow, will have been significantly suppressed. GW detected mergers have a Malmquist bias towards on-axis events (i.e. the rotational axis of the system), where the peak of the probability distribution is an inclination ∼ 30⁰. If the jets from these mergers have an intrinsic structure out to wider angles, then the majority of mergers will be accompanied by electromagnetic counterparts from these various jet structures. By making some simple assumptions about the energetic structure of a jet outside of a bright core region, the various temporal features that result from a given jet structure can be predicted. Where the population of merger jets is dominated by a single structure model, I show the expected fractions of optical counterparts brighter than m_AB = 21. On 17 August 2017, the Light Interferometer Gravitational Wave Observatory (LIGO) in collaboration with Virgo detected the merger of a binary neutron star system. Various electromagnetic counterparts were detected: the GRB 170817A by Fermi/GBM and INTEGRAL; an optical, blue to red, macro/kilo-nova from ∼ 1/2 day post merger to ∼ 5 − 10 days; and a brightening radio, and X-ray counterpart from ∼ 10 days. Optical detection of this counterpart at a magnitude ∼ 26 was made at ∼ 100 days post-merger. Analysis of this counterpart is consistent with the afterglow of a Gaussian structured jet viewed at the system inclination, ∼ 18 ± 8⁰. If all short GRB jets have a similar jet structure, then the rates of orphan afterglows in deep drilling blind surveys e.g. the Large Synoptic Survey Telescope (LSST), will be higher than those expected from a homogeneous, or ‘top-hat’ jet, population. The rates for the various jet structures for orphan afterglows from mergers are discussed, showing that for a population of failed GRBs, or an intrinsic Gaussian structure, an excess in the orphan rate may be apparent. Understanding the dynamics and structure for the jets from black-hole systems born at the merger of a compact binary can help give clues as to the nature of jets from black holes on all scales. As an aside, I show empirically that regardless of black hole mass or system phenomenology, the relativistic jets from such systems share a universal scaling for the jet power and emitted γ-ray luminosity. This scaling could be due to the similar efficiencies of various processes, or alternatively, the scaling may be able to give insights into the emission and physical processes that are responsible for high-energy photons from these outflows. GW astronomy offers a probe of the most extreme relativistic outflows in the Universe, GRBs. The predicted electromagnetic counterparts from these outflows, in association with GW detections, provides a way to probe the Lorentz-factor distribution for merger-jets. Additionally, the phenomenological shape of the afterglows, at various inclinations, gives an indication of the intrinsic structure of these jets. An understanding of these dynamical and structural qualities can be used to constrain the parent population, merger rates, and binary evolution models for compact binary systems.

We develop a method for implementing a proposal on utilizing knowledge of neutron star (NS) equation of state (EoS) for inferring the Hubble constant from a population of binary neutron star (BNS) mergers. This method is useful in exploiting BNSs as standard sirens when their redshifts are not available. Gravitational wave (GW) signals from compact object binaries provide a direct measurement of their luminosity distances, but not their redshifts. Unlike in the past, here we employ a realistic EoS parametrization in a Bayesian framework to simultaneously measure the Hubble constant and refine the constraints on the EoS parameters. The uncertainty in the redshift depends on the uncertainties in the EoS and the mass parameters estimated from GW data. Combining the inferred BNS redshifts with the corresponding luminosity distances, one constructs a redshift-distance relation and deduces the Hubble constant from it. Here, we show that in the Cosmic Explorer era, one can measure the Hubble constant to a precision of $\lesssim 5\%$ (with a $90\%$ credible interval) with a realistic distribution of a thousand BNSs, while allowing for uncertainties in their EoS parameters. Such a measurement can potentially resolve the current tension in the measurements of the Hubble constant from the early- and late-time universe. The methodology implemented in this work demonstrates a comprehensive prescription for inferring the NS EoS and the Hubble constant by simultaneously combining GW observations from merging NSs, while employing a simple population model for NS masses and keeping the merger rate of NSs constant in redshift. This method can be immediately extended to incorporate merger rate, population properties, and additional cosmological parameters.

By reaching through shrouding blastwaves, efficiently discovering off-axis events and probing the central engine at work, gravitational wave (GW) observations will soon revolutionize the study of gamma-ray bursts. Already, analyses of GW data targeting gamma-ray bursts have helped constrain the central engines of selected events. Advanced GW detectors with significantly improved sensitivities are under construction. After outlining the GW emission mechanisms from gamma-ray burst progenitors (binary coalescences, stellar core collapses, magnetars and others) that may be detectable with advanced detectors, we review how GWs will improve our understanding of gamma-ray burst central engines, their astrophysical formation channels and the prospects and methods for different search strategies. We place special emphasis on multimessenger searches. To achieve the most scientific benefit, GW, electromagnetic and neutrino observations should be combined to provide greater discriminating power and science reach.

In testing gravity a model-independent way, one of crucial tests is measuring the propagation speed of a gravitational wave (GW). In general relativity, a GW propagates with the speed of light, while in the alternative theories of gravity the propagation speed could deviate from the speed of light due to the modification of gravity or spacetime structure at a quantum level. Previously we proposed the method measuring the GW speed by directly comparing the arrival times between a GW and a photon from the binary merger of neutron stars or neutron star and black hole, assuming that it is associated with a short gamma-ray burst. The sensitivity is limited by the intrinsic time delay between a GW and a photon at the source. In this paper, we extend the method to distinguish the intrinsic time delay from the true signal caused by anomalous GW speed with multiple events at cosmological distances, also considering the redshift distribution of GW sources, redshift-dependent GW propagation speed, and the statistics of intrinsic time delays. We show that an advanced GW detector such as Einstein Telescope will constrain the GW propagation speed at the precision of ~10^{-16}. We also discuss the optimal statistic to measure the GW speed, performing numerical simulations.

Context. The last decade has seen the emergence of two new fields within astrophysics: gamma-ray polarimetry and gravitational wave (GW) astronomy. The former, which aims to measure the polarization of gamma rays in the energy range of tens to hundreds of keV, from astrophysical sources, saw the launch of the first dedicated polarimeters such as GAP and POLAR. Due to both a large scientific interest as well as their large signal-to-noise ratios, gamma-ray bursts (GRBs) are the primary source of interest of the first generation of polarimeters. Polarization measurements are theorized to provide a unique probe of the mechanisms at play in these extreme phenomena. On the other hand, GW astronomy started with the detection of the first black hole mergers by LIGO in 2015, followed by the first multi-messenger detection in 2017. Aims. While the potential of the two individual fields has been discussed in detail in the literature, the potential for joint observations has thus far been ignored. In this article, we aim to define how GW observations can best contribute to gamma-ray polarimetry and study the scientific potential of joint analyses. In addition, we aim to provide predictions on feasibility of such joint measurements in the near future. Methods. We study which GW observables can be combined with measurements from gamma-ray polarimetry to improve the discriminating power regarding GRB emission models. We then provide forecasts for the joint detection capabilities of current and future GW detectors and polarimeters. Results. Our results show that by adding GW data to polarimetry, a single precise joint detection would allow for the majority of emission models to be ruled out. We show that in the coming years, joint detections between GW and gamma-ray polarimeters might already be possible. Although these would allow one to constrain part of the model space, the probability of highly constraining joint detections will remain small in the near future. However, the scientific merit held by even a single such measurement makes it important to pursue such an endeavour. Furthermore, we show that using the next generation of GW detectors, such as the Einstein Telescope, joint detections for which GW data can better complement the polarization data become possible.

Determining the masses of neutron stars (NSs) accurately improves our understanding of the NS interior and complicated binary evolution. However, the masses of the systems are degenerate with the orbital inclination angle when using solely gravitational waves (GWs) or electromagnetic measurements, especially for face-on binaries. Taking advantage of both GWs and optical observations for LISA NS–white dwarf (WD) binaries, we propose a mass determination method utilising multimessenger observational information. By combining the binary mass function obtained from optical observations and a GW mass function, which we introduce, derived from GW observations, we demonstrate how we can set improved constraints on the NS mass and break the degeneracy in the mass and viewing inclination determination. We further comment on the universal relation of the error bar of the GW mass function versus the GW signal-to-noise ratio and propose a simple method for estimating the capability of using GW observations for mass determination with LISA. We show that for ultracompact NS–WD binaries within our Galaxy, the mass of the NS can be constrained to within an accuracy of ±0.2 M ⊙ with the proposed method.

After the successful detection of a gravitational-wave (GW) signal and its associated electromagnetic (EM) counterparts from GW170817, neutron star–black hole (NSBH) mergers have been highly expected to be the next type of multimessenger source. However, despite the detection of several NSBH merger candidates during the GW third observation run, no confirmed EM counterparts from these sources have been identified. The most plausible explanation is that these NSBH merger candidates were plunging events mainly because the primary black holes (BHs) had near-zero projected aligned spins based on GW observations. In view of the fact that neutron stars (NSs) can be easily tidally disrupted by BHs with high projected aligned spins, we study an evolution channel to form NSBH binaries with fast-spinning BHs, the properties of BH mass and spin, and their associated tidal disruption probability. We find that if the NSs are born first, the companion helium stars would be tidally spun up efficiently, and would thus finally form fast-spinning BHs. If BHs do not receive significant natal kicks at birth, these NSBH binaries that can merge within Hubble time would have BHs with projected aligned spins χ z ≳ 0.8 and, hence, can certainly allow tidal disruption to happen. Even if significant BH kicks are considered for a small fraction of NSBH binaries, the projected aligned spins of BHs are χ z ≳ 0.2. These systems can still be disrupted events unless the NSs are very massive. Thus, NS-first-born NSBH mergers would be promising multimessenger sources. We discuss various potential EM counterparts associated with these systems and their detectability in the upcoming fourth observation run.

It has been shown that time delay of strong gravitational lensing is not only an effective cosmological probe to constrain the Hubble constant, the matter density and the curvature of the universe but also useful for breaking the degeneracy of the dark energy equation of state and thus provide complementarity with other popular probes, such as type Ia supernovae, baryon acoustic oscillations, and cosmic microwave background radiation. Interestingly, compared to traditional strong lensing systems where quasars act as sources, strongly lensed gravitational waves (GWs) from compact binary coalescence and their electromagnetic (EM) counterparts systems have been recently proposed to be more powerful for studying cosmology since GWs and their EM counterparts are transients. Short durations of GWs and their EM counterparts are very advantageous to time delay measurement and lens profile modeling. Here, in the framework of Chevalier-Polarski-Linder (CPL) parametrization, we investigate improvement of constraining power on dark energy equation of state due to including time delay measurements of strong lensed GW systems. It is suggested that, on the basis of the third generation ground-based detector, e.g. Einstein Telescope, adding time delay of only 30 strong lensed GW systems to type Ia supernovae and cosmic microwave background radiation can improve the dark energy figure of merit by a factor 2. For the traditional standard siren method where the uncertainties of luminosity distances of GWs are $\sim10\%$, a few $\times~10^4$ events are expected to present similar constraints. In the precision cosmology era, this progress is of great significance for studying the nature of dark energy.

Gravitational waves (GWs) from compact binary stars at cosmological distances are promising and powerful cosmological probes, referred to as the GW standard sirens. With future GW detectors, we will be able to precisely measure source luminosity distances out to a redshift z∼5. To extract cosmological information, previously proposed cosmological studies using the GW standard sirens rely on source redshift information obtained through an extensive electromagnetic follow-up campaign. However, the redshift identification is typically time consuming and rather challenging. Here, we propose a novel method for cosmology with the GW standard sirens free from the redshift measurements. Utilizing the anisotropies of the number density and luminosity distances of compact binaries originated from the large-scale structure, we show that, once GW observations will be well established in the future, (i)these anisotropies can be measured even at very high redshifts (z≥2), where the identification of the electromagnetic counterpart is difficult, (ii)the expected constraints on the primordial non-Gaussianity with the Einstein Telescope would be comparable to or even better than the other large-scale structure probes at the same epoch, and (iii)the cross-correlation with other cosmological observations is found to have high-statistical significance, providing additional cosmological information at very high redshifts.

We investigate the distribution of neutron star masses in different populations of binaries, employing Bayesian statistical techniques. In particular, we explore the differences in neutron star masses between sources that have experienced distinct evolutionary paths and accretion episodes. We find that the distribution of neutron star masses in non-recycled eclipsing high-mass binaries as well as of slow pulsars, which are all believed to be near their birth masses, has a mean of 1.28 M_solar and a dispersion of 0.24 M_solar. These values are consistent with expectations for neutron star formation in core-collapse supernovae. On the other hand, double neutron stars, which are also believed to be near their birth masses, have a much narrower mass distribution, peaking at 1.33 M_solar but with a dispersion of only 0.05 M_solar. Such a small dispersion cannot easily be understood and perhaps points to a particular and rare formation channel. The mass distribution of neutron stars that have been recycled has a mean of 1.48 M_solar and a dispersion of 0.2 M_solar, consistent with the expectation that they have experienced extended mass accretion episodes. The fact that only a very small fraction of recycled neutron stars in the inferred distribution have masses that exceed ~2 M_solar suggests that only a few of these neutron stars cross the mass threshold to form low mass black holes.

We report on the population properties of compact binary mergers inferred from gravitational-wave observations of these systems during the first three LIGO-Virgo observing runs. The Gravitational-Wave Transient Catalog 3 (GWTC-3) contains signals consistent with three classes of binary mergers: binary black hole, binary neutron star, and neutron star–black hole mergers. We infer the binary neutron star merger rate to be between 10 and 1700 Gpc−3 yr−1 and the neutron star–black hole merger rate to be between 7.8 and 140 Gpc−3 yr−1, assuming a constant rate density in the comoving frame and taking the union of 90% credible intervals for methods used in this work. We infer the binary black hole merger rate, allowing for evolution with redshift, to be between 17.9 and 44 Gpc−3 yr−1 at a fiducial redshift (z=0.2). The rate of binary black hole mergers is observed to increase with redshift at a rate proportional to (1+z)κ with κ=2.9−1.8+1.7 for z≲1. Using both binary neutron star and neutron star–black hole binaries, we obtain a broad, relatively flat neutron star mass distribution extending from 1.2−0.2+0.1 to 2.0−0.3+0.3M⊙. We confidently determine that the merger rate as a function of mass sharply declines after the expected maximum neutron star mass, but cannot yet confirm or rule out the existence of a lower mass gap between neutron stars and black holes. We also find the binary black hole mass distribution has localized over- and underdensities relative to a power-law distribution, with peaks emerging at chirp masses of 8.3−0.5+0.3 and 27.9−1.8+1.9M⊙. While we continue to find that the mass distribution of a binary’s more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above approximately 60M⊙, which would indicate the presence of a upper mass gap. Observed black hole spins are small, with half of spin magnitudes below χi≈0.25. While the majority of spins are preferentially aligned with the orbital angular momentum, we infer evidence of antialigned spins among the binary population. We observe an increase in spin magnitude for systems with more unequal-mass ratio. We also observe evidence of misalignment of spins relative to the orbital angular momentum. Published by the American Physical Society 2023

Gravitational wave (GW) observations of coalescing compact binaries will be unique probes of strong-field, dynamical aspects of relativistic gravity. We present a short review of various schemes proposed in the literature to test general relativity (GR) and alternative theories of gravity using inspiral waveforms. Broadly these schemes may be classified into two types: model dependent and model independent. In the model dependent category, GW observations are compared against a specific waveform model representative of a particular theory or a class of theories such as scalar-tensor theories, dynamical Chern–Simons theory and massive graviton theories. Model independent tests are attempts to write down a parametrized gravitational waveform where the free parameters take different values for different theories and (at least some of) which can be constrained by GW observations. We revisit some of the proposed bounds in the case of downscaled LISA configuration (eLISA) and compare them with the original LISA configuration. We also compare the expected bounds on alternative theories of gravity from ground-based and space-based detectors and find that space-based GW detectors can test GR and other theories of gravity with unprecedented accuracies. We then focus on a recent proposal to use singular value decomposition of the Fisher information matrix to improve the accuracies with which post-Newtonian theory can be tested. We extend those results to the case of space-based detector eLISA and discuss its implications.

It has been predicted before the first detection of such events by LIGO, that the gravitational wave (GW) emission from compact binary coalescence can be used to constrain the expansion rate of the universe. Here we show how H0 can be derived purely from the GW of neutron star-black hole (NSBH) mergers using a Bayesian Hierarchical inference framework. This method provides an estimate of H0 spanning the redshift range, z < 0.25 with current GW sensitivity and without the need for any afterglow detection. We utilise the inherent distribution of neutron star masses together with the NSBH waveform amplitude and frequency to estimate distance and redshift respectively, thereby measuring H0 for the NSBH events up to systematics and observational uncertainties. With the most updated observations, our first estimate is $H_0 = 96^{+44}_{-50}$ km s−1 Mpc−1 for the likely NSBH events GW190426 and GW200115 using the proof-of-concept analysis outlined in this work. Taking into account of realistic detectors response and sensitivity, we forecast that soon, with 10 more such NSBH events we can reach competitive precision of δH0/H0 ≲ 20 %. It is valuable to have a completely independent, both in terms observations and relevant physics, method of determining H0 as is emphasized by the current controversy on H0 and attempts to explain it invoking new physics.

Gravitational waves were predicted by Albert Einstein more than a century ago, but only in 2015, the Laser Interferometer Gravitational-Wave Observatory (LIGO) was able to detect them. The gravitational wave phenomenon can be compared to spreading water from a lake after a stone has been thrown into it. Here, gravitational wave generation comes from an astronomical binary system formed by black holes or neutron stars. However, unlike the water, the amplitude of those gravitational waves is on a scale smaller than a proton’s size. Despite this, we can describe it with simple equations in a phenomenological way. We can model those waves on a regular computer using post Newtonian physics. Here we were able to generate gravitational waves from computational simulations and make its data analyze. When a gravitational wave is detected in the real-world problem, there is a great interest in establishing the physical features of the astronomical bodies involved in the process. In this way, we propose applying a simple neural network to receive the gravitational wave data and infer information about the astronomical bodies’ mass. The experimental results show that a simple neural network can extract mass information from the gravitational wave data. The recognition process proposed is much more straightforward than the complex computation based on numerical relativity for gravitational wave data analysis.

GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of state. For the binary neutron star hypothesis, we find that we cannot rule out the majority of theoretical models considered. In addition, the gravitational-wave data alone does not rule out the possibility that one or both objects were low-mass black holes. We discuss the possible outcomes in the case of a binary neutron star merger, finding that all scenarios from prompt collapse to long-lived or even stable remnants are possible. For long-lived remnants, we place an upper limit of 1.9 kHz on the rotation rate. If a black hole was formed any time after merger and the coalescing stars were slowly rotating, then the maximum baryonic mass of non-rotating neutron stars is at most , and three equations of state considered here can be ruled out. We obtain a tighter limit of for the case that the merger results in a hypermassive neutron star.