No Arabic abstract
Galactic ultra compact binaries are expected to be the dominant source of gravitational waves in the milli-Hertz frequency band. Of the tens of millions of galactic binaries with periods shorter than an hour, it is estimated that a few tens of thousand will be resolved by the future Laser Interferometer Space Antenna (LISA). The unresolved remainder will be the main source of ``noise between 1-3 milli-Hertz. Typical galactic binaries are millions of years from merger, and consequently their signals will persist for the the duration of the LISA mission. Extracting tens of thousands of overlapping galactic signals and characterizing the unresolved component is a central challenge in LISA data analysis, and a key contribution to arriving at a global solution that simultaneously fits for all signals in the band. Here we present an end-to-end analysis pipeline for galactic binaries that uses trans-dimensional Bayesian inference to develop a time-evolving catalog of sources as data arrive from the LISA constellation.
We discuss the detection of gravitational-wave backgrounds in the context of Bayesian inference and suggest a practical definition of what it means for a signal to be considered stochastic---namely, that the Bayesian evidence favors a stochastic signal model over a deterministic signal model. A signal can further be classified as Gaussian-stochastic if a Gaussian signal model is favored. In our analysis we use Bayesian model selection to choose between several signal and noise models for simulated data consisting of uncorrelated Gaussian detector noise plus a superposition of sinusoidal signals from an astrophysical population of gravitational-wave sources. For simplicity, we consider co-located and co-aligned detectors with white detector noise, but the method can be extended to more realistic detector configurations and power spectra. The general trend we observe is that a deterministic model is favored for small source numbers, a non-Gaussian stochastic model is preferred for intermediate source numbers, and a Gaussian stochastic model is preferred for large source numbers. However, there is very large variation between individual signal realizations, leading to fuzzy boundaries between the three regimes. We find that a hybrid, trans-dimensional model comprised of a deterministic signal model for individual bright sources and a Gaussian-stochastic signal model for the remaining confusion background outperforms all other models in most instances.
A number of works have shown that important information on the equation of state of matter at nuclear density can be extracted from the gravitational waves emitted by merging neutron-star binaries. We present a comprehensive analysis of the gravitational-wave signal emitted during the inspiral, merger and post-merger of 56 neutron-star binaries. This sample of binaries, arguably the largest studied to date with realistic equations of state, spans across six different nuclear-physics equations of state and ten masses, allowing us to sharpen a number of results recently obtained on the spectral properties of the gravitational-wave signal. Overall we find that: (i) for binaries with masses differing no more than $20%$, the frequency at gravitational-wave amplitudes maximum is related quasi-universally with the tidal deformability of the two stars; (ii) the spectral properties vary during the post-merger phase, with a transient phase lasting a few millisecond after the merger and followed by a quasi-stationary phase; (iii) when distinguishing the spectral peaks between these two phases, a number of ambiguities in the identification of the peaks disappear, leaving a simple and robust picture; (iv) using properly identified frequencies, quasi-universal relations are found between the spectral features and the properties of the neutron stars; (v) for the most salient peaks analytic fitting functions can be obtained in terms of the stellar tidal deformability or compactness. Altogether, these results support the idea that the equation of state of nuclear matter can be constrained tightly when a signal in gravitational waves from binary neutron stars is detected.
Ultralight bosons can form large clouds around stellar-mass black holes via the superradiance instability. Through processes such as annihilation, these bosons can source continuous gravitational wave signals with frequencies within the range of LIGO and Virgo. If boson annihilation occurs, then the Galactic black hole population will give rise to many gravitational signals; we refer to this as the ensemble signal. We characterize the ensemble signal as observed by the gravitational-wave detectors; this is important because the ensemble signal carries the primary signature that a continuous wave signal has a boson annihilation origin. We explore how a broad set of black hole population parameters affects the resulting spin-0 boson annihilation signal and consider its detectability by recent searches for continuous gravitational waves. A population of $10^8$ black holes with masses up to $30mathrm{M}_odot$ and a flat dimensionless initial spin distribution between zero and unity produces up to a thousand signals loud enough to be in principle detected by these searches. For a more moderately spinning population the number of signals drops by about an order of magnitude, still yielding up to a hundred detectable signals for some boson masses. A non-detection of annihilation signals at frequencies between 100 and 1200 Hz disfavors the existence of scalar bosons with rest energies between $2times10^{-13}$ and $2.5times10^{-12}$ eV. Finally we show that, depending on the black hole population parameters, care must be taken in assuming that the continuous wave upper limits from searches for isolated signals are still valid for signals that are part of a dense ensemble: Between 200 and 300 Hz, we urge caution when interpreting a null result for bosons between 4 and $6times10^{-13}$ eV.
Extending previous work by a number of authors, we have recently presented a new approach in which the detection of gravitational waves from merging neutron star binaries can be used to determine the equation of state of matter at nuclear density and hence the structure of neutron stars. In particular, after performing a large number of numerical-relativity simulations of binaries with nuclear equations of state, we have found that the post-merger emission is characterized by two distinct and robust spectral features. While the high-frequency peak was already shown to be associated with the oscillations of the hypermassive neutron star produced by the merger and to depend on the equation of state, we have highlighted that the low-frequency peak is related to the merger process and to the total compactness of the stars in the binary. This relation is essentially universal and provides a powerful tool to set tight constraints on the equation of state. We here provide additional information on the extensive analysis performed, illustrating the methods used, the tests considered, as well as the robustness of the results. We also discuss additional relations that can be deduced when exploring the data and how these correlate with various properties of the binary. Finally, we present a simple mechanical toy model that explains the main spectral features of the post-merger signal and can even reproduce analytically the complex waveforms emitted right after the merger.
The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star, a neutron star black hole binary and a binary black hole, where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with non-spinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence parameter space.