Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userWhen is a gravitational-wave signal stochastic?
133
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
Roughly every 2-10 minutes, a pair of stellar mass black holes merge somewhere in the Universe. A small fraction of these mergers are detected as individually resolvable gravitational-wave events by advanced detectors such as LIGO and Virgo. The rest contribute to a stochastic background. We derive the statistically optimal search strategy for a background of unresolved binaries. Our method applies Bayesian parameter estimation to all available data. Using Monte Carlo simulations, we demonstrate that the search is both safe and effective: it is not fooled by instrumental artefacts such as glitches, and it recovers simulated stochastic signals without bias. Given realistic assumptions, we estimate that the search can detect the binary black hole background with about one day of design sensitivity data versus $approx 40$ months using the traditional cross-correlation search. This framework independently constrains the merger rate and black hole mass distribution, breaking a degeneracy present in the cross-correlation approach. The search provides a unified framework for population studies of compact binaries, which is cast in terms of hyper-parameter estimation. We discuss a number of extensions and generalizations including: application to other sources (such as binary neutron stars and continuous-wave sources), simultaneous estimation of a continuous Gaussian background, and applications to pulsar timing.
We present a novel method for sampling iso-likelihood contours in nested sampling using a type of machine learning algorithm known as normalising flows and incorporate it into our sampler nessai. Nessai is designed for problems where computing the likelihood is computationally expensive and therefore the cost of training a normalising flow is offset by the overall reduction in the number of likelihood evaluations. We validate our sampler on 128 simulated gravitational wave signals from compact binary coalescence and show that it produces unbiased estimates of the system parameters. Subsequently, we compare our results to those obtained with dynesty and find good agreement between the computed log-evidences whilst requiring 2.07 times fewer likelihood evaluations. We also highlight how the likelihood evaluation can be parallelised in nessai without any modifications to the algorithm. Finally, we outline diagnostics included in nessai and how these can be used to tune the samplers settings.
A central challenge in Gravitational Wave Astronomy is identifying weak signals in the presence of non-stationary and non-Gaussian noise. The separation of gravitational wave signals from noise requires good models for both. When accurate signal models are available, such as for binary Neutron star systems, it is possible to make robust detection statements even when the noise is poorly understood. In contrast, searches for un-modeled transient signals are strongly impacted by the methods used to characterize the noise. Here we take a Bayesian approach and introduce a multi-component, variable dimension, parameterized noise model that explicitly accounts for non-stationarity and non-Gaussianity in data from interferometric gravitational wave detectors. Instrumental transients (glitches) and burst sources of gravitational waves are modeled using a Morlet-Gabor continuous wavelet frame. The number and placement of the wavelets is determined by a trans-dimensional Reversible Jump Markov Chain Monte Carlo algorithm. The Gaussian component of the noise and sharp line features in the noise spectrum are modeled using the BayesLine algorithm, which operates in concert with the wavelet model.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
