No Arabic abstract
For a uniform population of neutron stars whose spin-down is dominated by the emission of gravitational radiation, an old argument of Blandford states that the expected gravitational-wave amplitude of the nearest source is independent of the deformation and rotation frequency of the objects. Recent work has improved and extended this argument to set upper limits on the expected amplitude from neutron stars that also emit electromagnetic radiation. We restate these arguments in a more general framework, and simulate the evolution of such a population of stars in the gravitational potential of our galaxy. The simulations allow us to test the assumptions of Blandfords argument on a realistic model of our galaxy. We show that the two key assumptions of the argument (two dimensionality of the spatial distribution and a steady-state frequency distribution) are in general not fulfilled. The effective scaling dimension D of the spatial distribution of neutron stars is significantly larger than two, and for frequencies detectable by terrestrial instruments the frequency distribution is not in a steady state unless the ellipticity is unrealistically large. Thus, in the cases of most interest, the maximum expected gravitational-wave amplitude does have a strong dependence on the deformation and rotation frequency of the population. The results strengthen the previous upper limits on the expected gravitational-wave amplitude from neutron stars by a factor of 6 for realistic values of ellipticity.
The direct detection of continuous gravitational waves from pulsars is a much anticipated discovery in the emerging field of multi-messenger gravitational wave (GW) astronomy. Because putative pulsar signals are exceedingly weak large amounts of data need to be integrated to achieve desired sensitivity. Contemporary searches use ingenious ad-hoc methods to reduce computational complexity. In this paper we provide analytical expressions for the Fourier transform of realistic pulsar signals. This provides description of the manifold of pulsar signals in the Fourier domain, used by many search methods. We analyze the shape of the Fourier transform and provide explicit formulas for location and size of peaks resulting from stationary frequencies. We apply our formulas to analysis of recently identified outlier at 1891.76 Hz.
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.
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.
Observations of gravitational waves from compact binary mergers have enabled unique tests of general relativity in the dynamical and non-linear regimes. One of the most important such tests are constraints on the post-Newtonian (PN) corrections to the phase of the gravitational wave signal. The values of these PN coefficients can be calculated within standard general relativity, and these values are different in many alternate theories of gravity. It is clearly of great interest to constrain these deviations based on gravitational wave observations. In the majority of such tests which have been carried out, and which yield by far the most stringent constraints, it is common to vary these PN coefficients individually. While this might in principle be useful for detecting certain deviations from standard general relativity, it is a serious limitation. For example, we would expect alternate theories of gravity to generically have additional parameters. The corrections to the PN coefficients would be expected to depend on these additional non-GR parameters whence, we expect that the various PN coefficients to be highly correlated. We present an alternate analysis here using data from the binary neutron star coalescence GW170817. Our analysis uses an appropriate linear combination of non-GR parameters that represent absolute deviations from the corresponding post-Newtonian inspiral coefficients in the TaylorF2 approximant phase. These combinations represent uncorrelated non-GR parameters which correspond to principal directions of their covariance matrix in the parameter subspace. Our results illustrate good agreement with GR. In particular, the integral non-GR phase is $Psi_{mbox{non-GR}} = (0.447pm253)times10^{-1}$ and the deviation from GR percentile is $p^{mbox{Dev-GR}}_{n}=25.85%$.
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.