X-ray pulse profile modeling of PSR J0740+6620, the most massive known pulsar, with data from the NICER and XMM-Newton observatories recently led to a measurement of its radius. We investigate this measurements implications for the neutron star equat
ion of state (EoS), employing a nonparametric EoS model based on Gaussian processes and combining information from other x-ray, radio and gravitational-wave observations of neutron stars. Our analysis mildly disfavors EoSs that support a disconnected hybrid star branch in the mass-radius relation, a proxy for strong phase transitions, with a Bayes factor of $6.9$. For such EoSs, the transition mass from the hadronic to the hybrid branch is constrained to lie outside ($1,2$) $M_{odot}$. We also find that the conformal sound-speed bound is violated inside neutron star cores, which implies that the core matter is strongly interacting. The squared sound speed reaches a maximum of $0.75^{+0.25}_{-0.24}, c^2$ at $3.60^{+2.25}_{-1.89}$ times nuclear saturation density at 90% credibility. Since all but the gravitational-wave observations prefer a relatively stiff EoS, PSR J0740+6620s central density is only $3.57^{+1.3}_{-1.3}$ times nuclear saturation, limiting the density range probed by observations of cold, nonrotating neutron stars in $beta$-equilibrium.
Gravitational waves in general relativity contain two polarization degrees of freedom, commonly labeled plus and cross. Besides those two tensor modes, generic theories of gravity predict up to four additional polarization modes: two scalar and two v
ector. Detection of nontensorial modes in gravitational wave data would constitute a clean signature of physics beyond general relativity. Previous measurements have pointed to the unambiguous presence of tensor modes in gravitational waves, but the presence of additional generic nontensorial modes has not been directly tested. We propose a model-independent analysis capable of detecting and characterizing mixed tensor and nontensor components in transient gravitational wave signals, including those from compact binary coalescences. This infrastructure can constrain the presence of scalar or vector polarization modes on top of the tensor modes predicted by general relativity. Our analysis is morphology-independent (as it does not rely on a waveform templates), phase-coherent, and agnostic about the source sky location. We apply our analysis to data from GW190521 and simulated data and demonstrate that it is capable of placing upper limits on the strength of nontensorial modes when none are present, or characterizing their morphology in the case of a positive detection. Tests of the polarization content of a transient gravitational wave signal hinge on an extended detector network, wherein each detector observes a different linear combination of polarization modes. We therefore anticipate that our analysis will yield precise polarization constraints in the coming years, as the current ground-based detectors LIGO Hanford, LIGO Livingston, and Virgo are joined by KAGRA and LIGO India.
Transient non-gaussian noise in gravitational wave detectors, commonly referred to as glitches, pose challenges for inference of the astrophysical properties of detected signals when the two are coincident in time. Current analyses aim towards modeli
ng and subtracting the glitches from the data using a flexible, morphology-independent model in terms of sine-gaussian wavelets before the signal source properties are inferred using templates for the compact binary signal. We present a new analysis of gravitational wave data that contain both a signal and glitches by simultaneously modeling the compact binary signal in terms of templates and the instrumental glitches using sine-gaussian wavelets. The model for the glitches is generic and can thus be applied to a wide range of glitch morphologies without any special tuning. The simultaneous modeling of the astrophysical signal with templates allows us to efficiently separate the signal from the glitches, as we demonstrate using simulated signals injected around real O2 glitches in the two LIGO detectors. We show that our new proposed analysis can separate overlapping glitches and signals, estimate the compact binary parameters, and provide ready-to-use glitch-subtracted data for downstream inference analyses.
Finite-size effects on the gravitational wave signal from a neutron star merger typically manifest at high frequencies where detector sensitivity decreases. Proposed sensitivity improvements can give us access both to stronger signals and to a myriad
of weak signals from cosmological distances. The latter will outnumber the former and the relevant part of signal will be redshifted towards the detector most sensitive band. We study the redshift dependence of information about neutron star matter and find that single-scale properties, such as the star radius or the post-merger frequency, are better measured from the distant weak sources from $zsim 1$.
We present a systematic comparison of the binary black hole (BBH) signal waveform reconstructed by two independent and complementary approaches used in LIGO and Virgo source inference: a template-based analysis, and a morphology-independent analysis.
We apply the two approaches to real events and to two sets of simulated observations made by adding simulated BBH signals to LIGO and Virgo detector noise. The first set is representative of the 10 BBH events in the first Gravitational Wave Transient Catalog (GWTC-1). The second set is constructed from a population of BBH systems with total mass and signal strength in the ranges that ground based detectors are typically sensitive. We find that the reconstruction quality of the GWTC-1 events is consistent with the results of both sets of simulated signals. We also demonstrate a simulated case where the presence of a mismodelled effect in the observed signal, namely higher order modes, can be identified through the morphology-independent analysis. This study is relevant for currently progressing and future observational runs by LIGO and Virgo.
Estimates of the source parameters of gravitational-wave (GW) events produced by compact binary mergers rely on theoretical models for the GW signal. We present the first frequency-domain model for inspiral, merger and ringdown of the GW signal from
precessing binary-black-hole systems that also includes multipoles beyond the leading-order quadrupole. Our model, {tt PhenomPv3HM}, is a combination of the higher-multipole non-precessing model {tt PhenomHM} and the spin-precessing model {tt PhenomPv3} that includes two-spin precession via a dynamical rotation of the GW multipoles. We validate the new model by comparing to a large set of precessing numerical-relativity simulations and find excellent agreement across the majority of the parameter space they cover. For mass ratios $<5$ the mismatch improves, on average, from $sim6%$ to $sim 2%$ compared to {tt PhenomPv3} when we include higher multipoles in the model. However, we find mismatches $sim8%$ for the mass-ratio $6$ and highly spinning simulation. As a first application of the new model we have analysed the binary black hole event GW170729. We find larger values for the primary black hole mass of $58.25^{+11.73}_{-12.53} , M_odot$ (90% credible interval). The lower limit ($sim 46 , M_odot$) is comparable to the proposed maximum black hole mass predicted by different stellar evolution models due to the pulsation pair-instability supernova (PPISN) mechanism. If we assume that the primary ac{BH} in GW170729 formed through a PPISN then out of the four PPISN models we considered only the model of Woosley (2017) is consistent with our mass measurements at the 90% level.
Estimating the parameters of gravitational wave signals detected by ground-based detectors requires an understanding of the properties of the detectors noise. In particular, the most commonly used likelihood function for gravitational wave data analy
sis assumes that the noise is Gaussian, stationary, and of known frequency-dependent variance. The variance of the colored Gaussian noise is used as a whitening filter on the data before computation of the likelihood function. In practice the noise variance is not known and it evolves over timescales of dozens of seconds to minutes. We study two methods for estimating this whitening filter for ground-based gravitational wave detectors with the goal of performing parameter estimation studies. The first method uses large amounts of data separated from the specific segment we wish to analyze and computes the power spectral density of the noise through the mean-median Welch method. The second method uses the same data segment as the parameter estimation analysis, which potentially includes a gravitational wave signal, and obtains the whitening filter through a fit of the power spectrum of the data in terms of a sum of splines and Lorentzians. We compare these two methods and argue that the latter is more reliable for gravitational wave parameter estimation.
We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the
effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusion of higher-order modes in the models results in an improved estimate for the mass ratio of $(0.3-0.8)$ at the 90% credible level. Our updated measurement excludes equal masses at that level. We also find that models with higher-order modes lead to the data being more consistent with a smaller effective spin, with the probability that the effective spin is greater than zero being reduced from $99%$ to $94%$. The 90% credible interval for the effective spin parameter is now $(-0.01-0.50)$. Additionally, the recovered signal-to-noise ratio increases by $sim0.3$ units compared to analyses without higher-order modes. We study the effect of common spin priors on the derived spin and mass measurements, and observe small shifts in the spins, while the masses remain unaffected. We argue that our conclusions are robust against systematic errors in the waveform models. We also compare the above waveform-based analysis which employs compact-binary waveform models to a more flexible wavelet- and chirplet-based analysis. We find consistency between the two, with overlaps of $sim 0.9$, typical of what is expected from simulations of signals similar to GW170729, confirming that the data are well-described by the existing waveform models. Finally, we study the possibility that the primary component of GW170729 was the remnant of a past merger of two black holes and find this scenario to be indistinguishable from the standard formation scenario.
In the coming years gravitational-wave detectors will undergo a series of improvements, with an increase in their detection rate by about an order of magnitude. Routine detections of gravitational-wave signals promote novel astrophysical and fundamen
tal theory studies, while simultaneously leading to an increase in the number of detections temporally overlapping with instrumentally- or environmentally-induced transients in the detectors (glitches), often of unknown origin. Indeed, this was the case for the very first detection by the LIGO and Virgo detectors of a gravitational-wave signal consistent with a binary neutron star coalescence, GW170817. A loud glitch in the LIGO-Livingston detector, about one second before the merger, hampered coincident detection (which was initially achieved solely with LIGO-Hanford data). Moreover, accurate source characterization depends on specific assumptions about the behavior of the detector noise that are rendered invalid by the presence of glitches. In this paper, we present the various techniques employed for the initial mitigation of the glitch to perform source characterization of GW170817 and study advantages and disadvantages of each mitigation method. We show that, despite the presence of instrumental noise transients louder than the one affecting GW170817, we are still able to produce unbiased measurements of the intrinsic parameters from simulated injections with properties similar to GW170817.
We present a robust method to characterize the gravitational wave emission from the remnant of a neutron star coalescence. Our approach makes only minimal assumptions about the morphology of the signal and provides a full posterior probability distri
bution of the underlying waveform. We apply our method on simulated data from a network of advanced ground-based detectors and demonstrate the gravitational wave signal reconstruction. We study the reconstruction quality for different binary configurations and equations of state for the colliding neutron stars. We show how our method can be used to constrain the yet-uncertain equation of state of neutron star matter. The constraints on the equation of state we derive are complimentary to measurements of the tidal deformation of the colliding neutron stars during the late inspiral phase. In the case of a non-detection of a post-merger signal following a binary neutron star inspiral we show that we can place upper limits on the energy emitted.
