Related papers: GWDALI: A Fisher-matrix based software for gravita…
We study Bayesian inference methods for solving linear inverse problems, focusing on hierarchical formulations where the prior or the likelihood function depend on unspecified hyperparameters. In practice, these hyperparameters are often…
In this paper we investigate the gravitational waves emission by stellar dynamical structures as complex systems in the quadrupole approximation considering bounded and unbounded orbits. Precisely, after deriving analytical expressions for…
This short note reviews so-called Natural Gradient Descent (NGD) for multivariate Gaussians. The Fisher Information Matrix (FIM) is derived for several different parameterizations of Gaussians. Careful attention is paid to the symmetric…
Identifying the source parameters from a gravitational-wave measurement alone is limited by our ability to discriminate signals from different sources and the accuracy of the waveform family employed in the search. Here we address both…
We introduce a new template for the detection of gravitational waves from compact binary systems which is based on Chebyshev polynomials of the first kind. As well as having excellent convergence properties, these polynomials are also very…
An isotropic stochastic background of nanohertz gravitational waves creates excess residual power in pulsar-timing-array datasets, with characteristic inter-pulsar correlations described by the Hellings-Downs function. These correlations…
Data assimilation leads naturally to a Bayesian formulation in which the posterior probability distribution of the system state, given the observations, plays a central conceptual role. The aim of this paper is to use this Bayesian…
Rapid identification, characterization, and localization of gravitational waves from binary compact object mergers can enable well-informed follow-on multimessenger observations. In this work, we investigate a small modification to the RIFT…
The two-dimensional Hubbard model is studied using the variational quantum Monte Carlo technique with Gutzwiller-type variational wave functions. In addition to the simple one-site correlated Gutzwiller wave function, we use a form with…
Extreme mass-ratio inspirals~(EMRIs) detectable by the Laser Inteferometric Space Antenna~(LISA) are unique probes of astrophysics and fundamental physics. Parameter estimation for these sources is challenging, especially because the…
The set up of matched filters for the detection of gravitational waves from in-spiraling compact binaries is usually carried out using the restricted post-Newtonian approximation: the filter phase is modelled including post-Newtonian…
Natural Gradient Descent, a second-degree optimization method motivated by the information geometry, makes use of the Fisher Information Matrix instead of the Hessian which is typically used. However, in many cases, the Fisher Information…
Gravitational wave (GW) astronomy has opened new frontiers in understanding the cosmos, while the integration of artificial intelligence (AI) in science promises to revolutionize data analysis methodologies. However, a significant gap…
We propose a nonparametric density estimator based on the Gaussian process (GP) and derive three novel closed form learning algorithms based on Fisher divergence (FD) score matching. The density estimator is formed by multiplying a base…
Long-term precise timing of Galactic millisecond pulsars holds great promise for measuring the long-period (months-to-years) astrophysical gravitational waves. Several gravitational-wave observational programs, called Pulsar Timing Arrays…
We introduce a reduced basis approach as a new paradigm for modeling, representing and searching for gravitational waves. We construct waveform catalogs for non-spinning compact binary coalescences, and we find that for accuracies of 99%…
The Fisher matrix formalism has in recent times become the standard method for predicting the precision with which various cosmological parameters can be extracted from future data. This approach is fast, and generally returns accurate…
Inference and simulation in the context of high-dimensional dynamical systems remain computationally challenging problems. Some form of dimensionality reduction is required to make the problem tractable in general. In this paper, we propose…
Due to the highly non-convex nature of large-scale robust parameter estimation, avoiding poor local minima is challenging in real-world applications where input data is contaminated by a large or unknown fraction of outliers. In this paper,…
We consider estimation of undirected Gaussian graphical models and inverse covariances in high-dimensional scenarios by penalizing the corresponding precision matrix. While single $L_1$ (Graphical Lasso) and $L_2$ (Graphical Ridge)…