Related papers: Efficient Gravitational Wave Searches with Pulsar …
Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the stochastic background of gravitational waves (GWB)…
The goal of the MeerKAT radio telescope's pulsar timing array programme (MPTA) is the detection of gravitational waves (GWs) of nanohertz frequencies. Evidence for such a signal was recently announced by the MPTA and several other pulsar…
Pulsar timing arrays (PTAs) are being used to search for very low frequency gravitational waves. A gravitational wave signal appears in pulsar timing residuals through two components: one independent of and one dependent on the pulsar's…
In this paper, we propose Barrier Hamiltonian Monte Carlo (BHMC), a version of the HMC algorithm which aims at sampling from a Gibbs distribution $\pi$ on a manifold $\mathrm{M}$, endowed with a Hessian metric $\mathfrak{g}$ derived from a…
Gravitational waves (GWs) at ultra-low frequencies (${\lesssim 100\,\mathrm{nHz}}$) are key to understanding the assembly and evolution of astrophysical black hole (BH) binaries with masses $\sim 10^{6}-10^{9}\,M_\odot$ at low redshifts.…
In this paper, we discuss an extension of the Split Hamiltonian Monte Carlo (Split HMC) method for Gaussian process model (GPM). This method is based on splitting the Hamiltonian in a way that allows much of the movement around the state…
Pulsar timing arrays (PTA) are a promising probe to the cosmologically novel nanohertz gravitational wave (GW) regime through the stochastic GW background. In this work, we consider subluminal GW modes as a possible source of correlations…
Accurate and efficient estimation of rare events probabilities is of significant importance, since often the occurrences of such events have widespread impacts. The focus in this work is on precisely quantifying these probabilities, often…
Markov Chain Monte Carlo (MCMC) algorithms are essential tools in computational statistics for sampling from unnormalised probability distributions, but can be fragile when targeting high-dimensional, multimodal, or complex target…
Recent research show that the cosmological components of the Universe should influence on the propagation of Gravitational Waves (GWs) and even it has been proposed a new way to measure the cosmological constant using Pulsar Timing Arrays…
Gravitational wave (GW) searches using pulsar timing arrays (PTAs) are commonly assumed to be limited to a GW frequency of $\lesssim 4\times 10^{-7}$Hz given by the Nyquist rate associated with the average observational cadence of $2$ weeks…
Hamiltonian Monte Carlo is a widely used algorithm for sampling from posterior distributions of complex Bayesian models. It can efficiently explore high-dimensional parameter spaces guided by simulated Hamiltonian flows. However, the…
Pulsar timing arrays (PTAs) are anticipated to detect continuous gravitational waves (GWs) from individual supermassive black hole binaries (SMBHBs) in the near future. To identify the host galaxy of a GW source, PTAs require significantly…
Time-lapse images carry out important information about dynamic changes in Earth's interior which can be inferred using different Full Waveform Inversion (FWI) schemes. The estimation process is performed by manipulating more than one…
We review and expand on a Bayesian model selection technique for the detection of gravitational waves from neutron star ring-downs associated with pulsar glitches. The algorithm works with power spectral densities constructed from…
Particle Markov Chain Monte Carlo (PMCMC) is a general computational approach to Bayesian inference for general state space models. Our article scales up PMCMC in terms of the number of observations and parameters by generating the…
Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to…
High-precision pulsar timing is central to a wide range of astrophysics and fundamental physics applications. When timing an ensemble of millisecond pulsars in different sky positions, known as a pulsar timing array (PTA), one can search…
Pulsar timing arrays (PTAs) detect gravitational waves (GWs) via the correlations they create in the arrival times of pulses from different pulsars. The mean correlation, a function of the angle $\gamma$ between the directions to two…
Markov Chain Monte Carlo (MCMC) is a well-established family of algorithms primarily used in Bayesian statistics to sample from a target distribution when direct sampling is challenging. Existing work on Bayesian decision trees uses MCMC.…