Recent Gravitational Wave Observations Research Articles

Analysis of the Principle of Gravitational Wave Detection and State- of-art Results

Over one hundred years after Einstein’s first theoretical prediction of gravitational waves, the recent observations of gravitational waves by interferometric detectors have risen excitement in the science community with the light from a new era of directly observing perturbations in spacetime itself. With the large quantity of new research in the gravitational wave detection area, this study provides a brief summary of the current achievements of the field. This paper first gives a short introduction of the field, then offers a summarize of the theories and principles of gravitational wave detection. Afterwards, it demonstrates modern detectors using the advanced LIGO detectors as an example, and presents its groundbreaking results, signal GW150914 and GW151226. Finally, the limitations and prospects of current observations are proposed. By summarizing fundamental concepts and major achievements in the field of gravitational wave detection, it is hoped to provide future students and researchers with a basic idea of the field and thus promote new ideas and advancements in the rising era of gravitational astronomy.

Black hole-neutron star mergers in Einstein-scalar-Gauss-Bonnet gravity

Gravitational wave observations of black hole-neutron star binaries, particularly those where the black hole has a lower mass compared to other observed systems, have the potential to place strong constraints on modifications to general relativity that arise at small curvature length scales. Here we study the dynamics of black hole-neutron star mergers in shift-symmetric Einstein-scalar-Gauss-Bonnet gravity, a representative example of such a theory, by numerically evolving the full equations of motion. We consider quasi-circular binaries with different mass-ratios that are consistent with recent gravitational wave observations, including cases with and without tidal disruption of the star, and quantify the impact of varying the coupling controlling deviations from general relativity on the gravitational wave signal and scalar radiation. We find that the main effect on the late inspiral is the accelerated frequency evolution compared to general relativity, and that—even considering Gauss-Bonnet coupling values approaching those where the theory breaks down—the impact on the merger gravitational wave signal is mild, predominately manifesting as a small change in the amplitude of the ringdown. We compare our results to current post-Newtonian calculations and find consistency throughout the inspiral. Published by the American Physical Society 2024

Which Black Hole Is Spinning? Probing the Origin of Black Hole Spin with Gravitational Waves

Theoretical studies of angular momentum transport suggest that isolated stellar-mass black holes are born with negligible dimensionless spin magnitudes χ ≲ 0.01. However, recent gravitational-wave observations indicate ≳40% of binary black hole systems contain at least one black hole with a nonnegligible spin magnitude. One explanation is that the firstborn black hole spins up the stellar core of what will become the second-born black hole through tidal interactions. Typically, the second-born black hole is the “secondary” (less massive) black hole though it may become the “primary” (more massive) black hole through a process known as mass-ratio reversal. We investigate this hypothesis by analyzing data from the third gravitational-wave transient catalog using a “single-spin” framework in which only one black hole may spin in any given binary. Given this assumption, we show that at least 28% (90% credibility) of the LIGO–Virgo–KAGRA binaries contain a primary with significant spin, possibly indicative of mass-ratio reversal. We find no evidence for binaries that contain a secondary with significant spin. However, the single-spin framework is moderately disfavored (natural log Bayes factor lnB=3.1 ) when compared to a model that allows both black holes to spin. If future studies can firmly establish that most merging binaries contain two spinning black holes, it may call into question our understanding of formation mechanisms for binary black holes or the efficiency of angular momentum transport in black hole progenitors.

General relativistic simulations of collapsing binary neutron star mergers with Monte Carlo neutrino transport

Recent gravitational wave observations of neutron-star-neutron-star and neutron-star-black-hole binaries appear to indicate that massive neutron stars may not be too uncommon in merging systems. These discoveries have led to an increased interest in the simulation of merging compact binaries involving massive stars. In this paper, we present a first set of evolution of massive neutron star binaries using Monte Carlo radiation transport for the evolution of neutrinos. We study a range of systems, from nearly symmetric binaries that collapse to a black hole before forming a disk or ejecting material, to more asymmetric binaries in which tidal disruption of the lower mass star leads to the production of more interesting postmerger remnants. For the latter type of systems, we additionally study the impact of viscosity on the properties of the outflows, and compare our results to two recent simulations of identical binaries performed with the whiskythc code. We find agreement on the black hole properties, disk mass, and mass and velocity of the outflows within expected numerical uncertainties, and some minor but noticeable differences in the evolution of the electron fraction when using a subgrid viscosity model, with viscosity playing a more minor role in our simulations. The method used to account for r-process heating in the determination of the outflow properties appears to have a larger impact on our result than those differences between numerical codes. We also use the simulation with the most ejected material to verify that our newly implemented Lagrangian tracers provide a reasonable sampling of the matter outflows as they leave the computational grid. We note that, given the lack of production of hot outflows in these mergers, the main role of neutrinos in these systems is to set the composition of the postmerger remnant. One of the main potential uses of our simulations is, thus, as improved initial conditions for longer evolutions of such remnants.

Structure and stability of differentially rotating compact stellar objects

AbstractDepending on the dynamics of a binary neutron star merger, the collision may result in a differentially rotating compact object. Differentially rotating stars can sustain a total mass considerably higher than that of a uniformly rotating star, giving rise to “hypermassive” objects like hypermassive neutron stars. These stars are likely to exhibit extreme structural deformation along the radial axis due to their high masses. Both the increased mass and structural deformations supported by differential rotation allow the compact remnant to remain stable in otherwise unstable configurations on short, dynamical timescales. In this work, we numerically simulate differentially rotating neutron stars to explore an increase in mass and structural deformation for three relativistic mean‐field equations of state models. Results are used to predict outcomes for recent gravitational wave observations of binary neutron star mergers.

Anisotropic solution for polytropic stars in 4D Einstein–Gauss–Bonnet gravity

In the present work we have investigated a new anisotropic solution for polytropic stars in the framework of 4D Einstein–Gauss–Bonnet (EGB) gravity. The possibility of determining the masses and radii of compact stars which puts some limitations on equation of state (EoS) above the nuclear saturation density. For this purpose, the 4D EGB field equations are solved by taking a generalized polytropic equation of state (EoS) with Finch–Skea ansatz. The generalized solution for anisotropic model has been tested for different values of Gauss–Bonnet constant alpha which satisfies all the physical criteria including causality with static stability via mass vs central mass density (M-rho _c), Bondi and Abreu criterion. The adiabatic index shows a minor influence of the GB coupling constant whereas the central and surface redshifts in the EGB gravity always remain lower than the GR. We present the possibility of fitting the mass and radius for some known compact star via M-R curve which satisfies the recent gravitational wave observations from GW 170817 event.

Chaos in quadratic gravity

While recent gravitational wave observations by LIGO and Virgo allow for tests of general relativity in the extreme gravity regime, these observations are still blind to a large swath of phenomena outside these instruments' sensitivity curves. Future gravitational-wave detectors, such as LISA, will enable probes of longer-duration and lower-frequency events. In particular, LISA will enable the characterization of the non-linear dynamics of extreme mass-ratio inspirals, when a small compact object falls into a supermassive black hole. In this paper, we study the motion of test particles around spinning black holes in two quadratic gravity theories: scalar Gauss-Bonnet and dynamical Chern-Simons gravity. We show that geodesic trajectories around rotating black holes in these theories are likely to not have a fourth constant of the motion. In particular, we show that Poincar\'e sections of the orbital phase space present chaotic features that will affect the inspiral of small compact objects into supermassive black holes in these theories. Nevertheless, the characteristic size of these chaotic features is tiny and their location in parameter space is very close to the event horizon of the supermassive black hole. Therefore, the detection of such chaotic features with LISA is likely very challenging, at best.

Theory-agnostic framework for inspiral tests of general relativity with higher-harmonic gravitational waves

Recent gravitational wave observations show evidence for the presence of higher harmonics, thus possibly indicating that these waves were generated in the inspiral of compact objects with asymmetric mass ratios. Signals with higher harmonics contain a trove of information that can lead to a better estimation of system parameters and possibly to more stringent tests of general relativity. Gravitational wave model that include higher harmonics, however, have only been developed within general relativity, while models to test theory-agnostic deviations from general relativity have been purely based on the signal's dominant mode. We here extend the parameterized post-Einsteinian framework to include the $\ell=2, 3$ and $4$ higher harmonics to first post-Newtonian order, therefore providing a ready-to-use Fourier-domain waveform model for tests of general relativity with higher harmonics. We find that the deformations to the higher harmonics of the Fourier phase can be easily mapped to the deformation of the dominant harmonic, while the deformations to the higher-harmonics of the Fourier amplitude in general cannot in a theory-agnostic way. Nonetheless, we develop a simple ansatz for the deformations of the waveform amplitude (through a re-scaling deformation of the time-domain amplitude) that both minimizes the number of independent amplitude deformations parameters and captures the predictions of all known modified theories to date.

Constraints on Einstein-dilation-Gauss-Bonnet gravity from black hole-neutron star gravitational wave events

Recent gravitational wave observations allow us to probe gravity in the strong and dynamical field regime. In this paper, we focus on testing Einstein-dilation Gauss-Bonnet gravity which is motivated by string theory. In particular, we use two new neutron star black hole binaries (GW200105 and GW200115). We also consider GW190814 which is consistent with both a binary black hole and a neutron star black hole binary. Adopting the leading post-Newtonian correction and carrying out a Bayesian Markov-chain Monte Carlo analyses, we derive the 90\% credible upper bound on the coupling constant of the theory as $\sqrt{\alpha_{GB}} \lesssim 1.33\,\rm km$, whose consistency is checked with an independent Fisher analysis. This bound is stronger than the bound obtained in previous literature by combining selected binary black hole events in GWTC-1 and GWTC-2 catalogs. We also derive a combined bound of $\sqrt{\alpha_{GB}} \lesssim 1.18\,\rm km$ by stacking GW200105, GW200115, GW190814, and selected binary black hole events. In order to check the validity of the effect of higher post-Newtonian terms, we derive corrections to the waveform phase up to second post Newtonian order by mapping results in scalar-tensor theories to Einstein-dilation Gauss-Bonnet gravity. We find that such higher-order terms improve the bounds by $14.5\%$ for GW200105 and $6.9\%$ for GW200115 respectively.

A template-free approach for waveform extraction of gravitational wave events

We develop a general data-driven and template-free method for the extraction of event waveforms in the presence of background noise. Recent gravitational-wave observations provide one of the significant scientific areas requiring data analysis and waveform extraction capability. We use our method to find the waveforms for the reported events from the first, second, and third LIGO observation runs (O1, O2, and O3). Using the instantaneous frequencies derived by the Hilbert transform of the extracted waveforms, we provide the physical time delays between the arrivals of gravitational waves to the detectors.

Failure of the Fisher matrix when including tidal terms: Considering construction of template banks of tidally deformed binary neutron stars

Recent gravitational-wave observations have begun to constrain the internal physics of neutron stars. However, current detection searches for neutron star systems assume that potential neutron stars are low-mass black holes, ignoring any affect on the gravitational-wave signal due to the internal neutron-star physics. We wish to create a template bank of binary neutron star waveforms including the effect of tidal deformability. However, we find that the Fisher matrix, which is commonly used to approximate match calculations when placing template banks, is unsuitable to predict the match between two binary neutron star waveforms. We find that the Fisher matrix can predict errors on the mismatch that are larger than $100\%$ when attempting to identify waveforms with a match of $0.97$. We explore the regime in which the Fisher matrix cannot be trusted and examine why it breaks down. We demonstrate that including higher-order terms in the Taylor series expansion of the match can reliably compute matches for these examples, but that it is prohibitively computationally expensive to do so. Finally, we demonstrate that stochastic placement can still be used to construct a template bank of tidally deformed neutron-star waveforms.

Complexity analysis of gravitational waves

The recent observations of gravitational waves during the coalescence of black holes has confirmed various predictions of general relativity. However, some concerns have been raised that these putative waves are artifacts. We have analyzed the ten confident gravitational waves detected in Observation Run 1 (O1) and Observation Run 2 (O2) and the gravitational wave emitted by the single known example of a neutron star coalescence, all of which have been verified by the Advanced Laser Interferometer Gravitational Wave Observatory and/or the Advanced Virgo Collaboration. On three occasions both Advanced LIGO and Advanced Virgo detected the same event (Abbott et al., 2019).Whitening of the Hanford data (Gravitational Wave Open Science Center) was employed in this study to reduce colored non-transient noise, but template matching was not employed. In each of the ten events mentioned above we verified the presence of such waves by the application of non-stationary data analysis via the Morlet wavelet. We examined the strain data before, during and after each wave. In each case we employed seven measures of complexity or information before, and during each wave and analyzed these data using doubly multivariate MANOVA, the appropriate statistical test for repeated measures (pre-wave vs wave). In addition we compared the strain data before and after the putative event to see if there was persistence of any effect of the event. In each case the information content was significantly elevated during the wave compared to similar time samples taken before the wave. No persistent effects on the strain data were noted after the event. We conclude our analysis is consistent with the actual detection of gravitational waves by the LIGO and Virgo Collaboration. Our results suggest black hole coalescence generates mathematically complex gravitational waves which may elucidate aspects of such coalescence.

Recent Observations of Gravitational Waves by LIGO and Virgo Detectors

In this paper we present the most recent observations of gravitational waves (GWs) by LIGO and Virgo detectors. We also discuss contributions of the recent Nobel prize winner, Sir Roger Penrose to understanding gravitational radiation and black holes (BHs). We make a short introduction to GW phenomenon in general relativity (GR) and we present main sources of detectable GW signals. We describe the laser interferometric detectors that made the first observations of GWs. We briefly discuss the first direct detection of GW signal that originated from a merger of two BHs and the first detection of GW signal form merger of two neutron stars (NSs). Finally we present in more detail the observations of GW signals made during the first half of the most recent observing run of the LIGO and Virgo projects. Finally we present prospects for future GW observations.

Superfluid Neutron Matter with a Twist

Superfluid neutron matter is a key ingredient in the composition of neutron stars. The physics of the inner crust are largely dependent on those of its S-wave neutron superfluid, which has made its presence known through pulsar glitches and modifications in neutron star cooling. Moreover, with recent gravitational-wave observations of neutron star mergers, the need for an equation of state for the matter of these compact stars is further accentuated and a model-independent treatment of neutron superfluidity is important. Ab initio techniques developed for finite systems can be guided to perform extrapolations to the thermodynamic limit and attain this model-independent extraction of various quantities of infinite superfluid neutron matter. To inform such an extrapolation scheme, we performed calculations of the neutron 1S0 pairing gap using model-independent odd–even staggering in the context of the particle-conserving, projected Bardeen–Cooper–Schrieffer (BCS) theory under twisted boundary conditions. While the practice of twisted boundary conditions is standard in solid-state physics and has been used repeatedly in the past to reduce finite-size effects, this is the first time that it has been employed in the context of pairing. We find that a twist-averaging approach results in a substantial reduction of the finite-size effects, bringing systems with N⪆50 within a 2% error margin from the infinite system. This can significantly reduce extrapolation-related errors in the extraction of superfluid neutron matter quantities.

Accessing the axion via compact object binaries

Black holes in binaries with other compact objects can provide natural venues for indirect detection of axions or other ultralight fields. The superradiant instability associated with a rapidly spinning black hole leads to the creation of an axion cloud which carries energy and angular momentum from the black hole. This cloud will then decay via gravitational wave emission. We show that the energy lost as a result of this process tends toward an outspiraling of the binary orbit. A given binary system is sensitive to a narrow range of axion masses, determined by the mass of the black hole. This proposal provides a complementary alternative to other approaches for detecting or constraining light particles created by superradiance, such as directly measuring the black hole spin or detecting the resulting gravitational wave signal. Pulsar-black hole binaries, once detected in the electromagnetic band, will allow high-precision measurements of black hole mass loss via timing measurements of the companion pulsar. This avenue of investigation is particularly promising in light of the recent preliminary announcements of two candidate black hole-neutron star mergers by LIGO/Virgo (#S190814bv and #S190426c). We demonstrate that for such a binary system with a typical millisecond pulsar and a 3M⊙ black hole, axions with masses between 2.7 × 10−12 eV and 3.2 × 10−12 eV are detectable. Recent gravitational wave observations by LIGO/Virgo of binary black hole mergers imply that, for these binaries, gravitational radiation from the rotating quadrupole moment is a dominant effect, causing an inspiraling orbit. With some reasonable assumptions about the period of the binary when it formed and the spins of the black holes, these observations rule out possible axion masses between 3 × 10−13 eV and 6 × 10−13 eV . Future binary black hole observations, for example by LISA, are expected to provide more robust bounds. In some circumstances, neutron stars may also undergo superradiant instabilities, and binary pulsars could be used to exclude axions with certain masses and matter couplings.

Parameter estimation for tests of general relativity with the astrophysical stochastic gravitational wave background

Recent observations of gravitational waves from binary black holes and neutron stars allow us to probe the strong and dynamical field regime of gravity. On the other hand, a collective signal from many individual, unresolved sources results in what is known as a stochastic background. We here consider probing gravity with such a background from stellar-mass binary black hole mergers. We adopt a simple power-law spectrum and carry out a parameter estimation study with a network of current and future ground-based detectors by including both general relativistic and beyond general relativistic variables. For a network of second-generation detectors, we find that one can place meaningful bounds on the deviation parameter in the gravitational-wave amplitude if it enters at a sufficiently negative post-Newtonian order. However, such future bounds from a stochastic background are weaker than existing bounds from individual sources, such as GW150914 and GW151226. We also find that systematic errors due to mismodeling of the spectrum is much smaller than statistical errors, which justifies our use of the power-law model. Regarding a network of third-generation detectors, we find that the bounds on the deviation parameter from statistical errors improve upon the second-generation case, though systematic errors now dominate the error budget and thus one needs to use a more realistic spectrum model. We conclude that individual sources seem to be more powerful in probing general relativity than the astrophysical stochastic background.

Estimates for disk and ejecta masses produced in compact binary mergers

There is irresistible observational evidence that binary systems of compact objects with at least one neutron star are progenitors of short gamma-ray bursts, as well as a production site for r-process elements, at least when some matter is ejected by the merger and an accretion disk is formed. The recent observations of gravitational waves in conjunction with electromagnetic counterparts fuel the need for models predicting the outcome of a given merger and the properties of the associated matter outflows as a function of the initial parameters of the binary. In this manuscript, we provide updated fitting formulae that estimate the disk mass for double neutron star binaries and ejecta masses for black hole-neutron star and double neutron star binaries, fitted to the results of numerical simulations. Our proposed fitting formulae improve on existing models by aiming for analytical simplicity, by covering a larger region of parameter space, and by accounting for regions of parameter space not covered by numerical simulations but with physically manifest merger outcomes.

Reviving Horndeski theory using teleparallel gravity after GW170817

Horndeski gravity was highly constrained from the recent gravitational wave observations by the LIGO Collaboration down to $|c_{g}/c-1|\gtrsim 10^{-15}$. In this Letter we study the tensorial perturbations in a flat cosmological background for an analogue version of Horndenki gravity which is based in Teleparallel Gravity constructed from a flat manifold with a nonvanishing torsion tensor. It is found that in this approach, one can construct a more general Horndeski theory satisfying $c_T=c_g/c=1$ without eliminating the coupling functions $G_5(\phi,X)$ and $G_4(\phi,X)$ that were highly constrained in standard Horndeski theory. Hence, in the Teleparallel approach one is able to restore these terms, creating an interesting way to revive Horndeski gravity.

GW170817 constraints analyzed with Gogny forces and momentum-dependent interactions

A set of equations of state obtained from finite-range Gogny forces and momentum-dependent interactions is used to investigate the recent observation of gravitational waves from the binary neutron star merger GW170817 event. For this set of interactions, we have calculated the neutron star tidal deformabilities (related to the second Love number), the mass-radius diagram, and the moment of inertia (I). The I-Love relation has been verified. We also have found strong correlations among the tidal deformability of the canonical neutron star, its radius, and the derivatives of the nuclear symmetry energy at the saturation density. Most of the obtained results are located within the constraints of the tidal deformabilities extracted from the GW170817 detection.

Neutron and quark stars: constraining the parameters for simple EoS using the GW170817

It is well known that the equation of state (EoS) of compact objects like neutron and quark stars is not determined despite there are several sophisticated models to describe it. From the electromagnetic observations, summarized in \cite{Lattimer01}, and the recent observation of gravitational waves from binary neutron star inspiral GW170817 \cite{Abbott2017_etal} and GW190425 \cite{Abbott2019}, it is possible to make an estimation of the range of masses and so constraint the mass of the neutron and quark stars, determining not only the best approximation for the EoS, but which kind of stars we would be observing. In this paper we explore several configurations of neutron stars assuming a simple polytropic equation of state, using a single layer model without crust. In particular, when the EoS depends on the mass rest density, $p=K \rho_{0}^{\Gamma}$, and when it depends on the energy density $p=K \rho^{\Gamma}$, considerable differences in the mass-radius relationships are found. On the other hand, we also explore quark stars models using the MIT bag EoS for different values of the vacuum energy density $B$.