Related papers: EMRI_MC: A GPU-based Python code for Bayesian infe…
We introduce $\texttt{GPUmonty}$, a CUDA/C-based Monte Carlo radiative transfer code accelerated using graphics processing units (GPUs). $\texttt{GPUmonty}$ derives from the CPU-based code $\texttt{grmonty}$ and offloads the most…
Gravitational waves from extreme mass-ratio inspirals (EMRIs), the inspirals of stellar-mass compact objects into massive black holes, are predicted to be observed by the Laser Interferometer Space Antenna (LISA). A sufficiently large…
A novel multi-task Gaussian process (GP) framework is proposed, by using a common mean process for sharing information across tasks. In particular, we investigate the problem of time series forecasting, with the objective to improve…
We present a novel algorithm that is based on a Bayesian Markov Chain Monte Carlo (MCMC) technique for performing robust profile analysis of a data cube from either single-dish or interferometric radio telescopes. It fits a set of models…
The PE for GW merger events relies on a waveform model calibrated using numerical simulations. Within the Bayesian framework, this waveform model represents the GW signal produced during the merger and is crucial for estimating the…
Bayesian spectral deconvolution provides a data-driven framework for mathematical model selection and parameter estimation from spectral data. Although highly versatile, it becomes computationally expensive as the number of model…
We present a forecast for the upcoming Einstein Telescope (ET) interferometer with two new methods to infer cosmological parameters. We consider the emission of Gravitational Waves (GWs) from compact binary coalescences, whose…
We introduce an algorithm to marginalize the likelihood for a gravitational wave signal from a quasi-circular binary merger over its extrinsic parameters, accounting for the effects of higher harmonics and spin-induced precession. The…
We introduce $\texttt{GWFAST}$, a Fisher information matrix $\texttt{Python}$ code that allows easy and efficient estimation of signal-to-noise ratios and parameter measurement errors for large catalogs of resolved sources observed by…
Extreme Mass Ratio Inspirals (EMRIs) are key observational targets for the Laser Interferometer Space Antenna (LISA) mission. Unresolvable EMRI signals contribute to the formation of a gravitational wave background (GWB). Characterizing the…
Extreme mass ratio inspirals (EMRIs) provide unique probes of near-horizon dissipation through the tidal heating. We present a full Bayesian analysis of tidal heating in equatorial eccentric EMRIs by performing injection-recovery studies…
The Bayesian approach to Inverse Problems relies predominantly on Markov Chain Monte Carlo methods for posterior inference. The typical nonlinear concentration of posterior measure observed in many such Inverse Problems presents severe…
Bayesian graphical models are a useful tool for understanding dependence relationships among many variables, particularly in situations with external prior information. In high-dimensional settings, the space of possible graphs becomes…
This document serves to complement our website which was developed with the aim of exposing the students to Gaussian Processes (GPs). GPs are non-parametric Bayesian regression models that are largely used by statisticians and geospatial…
Geodesic Monte Carlo (gMC) is a powerful algorithm for Bayesian inference on non-Euclidean manifolds. The original gMC algorithm was cleverly derived in terms of its progenitor, the Riemannian manifold Hamiltonian Monte Carlo (RMHMC). Here,…
Folding uncertainty in theoretical models into Bayesian parameter estimation is necessary in order to make reliable inferences. A general means of achieving this is by marginalizing over model uncertainty using a prior distribution…
An important step in the planning of future gravitational-wave (GW) detectors and of the networks they will form is the estimation of their detection and parameter-estimation capabilities, which is the basis of science-case studies. Several…
MRChem is a code for molecular electronic structure calculations, based on a multiwavelet adaptive basis representation. We provide a description of our implementation strategy and several benchmark calculations. Systems comprising more…
We take a new look at parameter estimation for Gaussian Mixture Models (GMMs). In particular, we propose using \emph{Riemannian manifold optimization} as a powerful counterpart to Expectation Maximization (EM). An out-of-the-box invocation…
In this paper we discuss the development of a fast and accurate waveform model for the quasi-circular orbital evolution of extreme-mass-ratio-inspirals (EMRIs). This model simply employs the data of a few numerical Teukoulsky-based energy…