Researcher profile

Luigi Ferraioli

Luigi Ferraioli contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2014arXiv

Bayesian Model Selection for LISA Pathfinder

The main goal of the LISA Pathfinder (LPF) mission is to fully characterize the acceleration noise models and to test key technologies for future space-based gravitational-wave observatories similar to the eLISA concept. The data analysis team has developed complex three-dimensional models of the LISA Technology Package (LTP) experiment on-board LPF. These models are used for simulations, but more importantly, they will be used for parameter estimation purposes during flight operations. One of the tasks of the data analysis team is to identify the physical effects that contribute significantly to the properties of the instrument noise. A way of approaching this problem is to recover the essential parameters of a LTP model fitting the data. Thus, we want to define the simplest model that efficiently explains the observations. To do so, adopting a Bayesian framework, one has to estimate the so-called Bayes Factor between two competing models. In our analysis, we use three main different methods to estimate it: The Reversible Jump Markov Chain Monte Carlo method, the Schwarz criterion, and the Laplace approximation. They are applied to simulated LPF experiments where the most probable LTP model that explains the observations is recovered. The same type of analysis presented in this paper is expected to be followed during flight operations. Moreover, the correlation of the output of the aforementioned methods with the design of the experiment is explored.

preprint2014arXiv

Data series subtraction with unknown and unmodeled background noise

LISA Pathfinder (LPF), ESA's precursor mission to a gravitational wave observatory, will measure the degree to which two test-masses can be put into free-fall, aiming to demonstrate a residual relative acceleration with a power spectral density (PSD) below 30 fm/s$^2$/Hz$^{1/2}$ around 1 mHz. In LPF data analysis, the measured relative acceleration data series must be fit to other various measured time series data. This fitting is required in different experiments, from system identification of the test mass and satellite dynamics to the subtraction of noise contributions from measured known disturbances. In all cases, the background noise, described by the PSD of the fit residuals, is expected to be coloured, requiring that we perform such fits in the frequency domain. This PSD is unknown {\it a priori}, and a high accuracy estimate of this residual acceleration noise is an essential output of our analysis. In this paper we present a fitting method based on Bayesian parameter estimation with an unknown frequency-dependent background noise. The method uses noise marginalisation in connection with averaged Welch's periodograms to achieve unbiased parameter estimation, together with a consistent, non-parametric estimate of the residual PSD. Additionally, we find that the method is equivalent to some implementations of iteratively re-weighted least-squares fitting. We have tested the method both on simulated data of known PSD, and to analyze differential acceleration from several experiments with the LISA Pathfinder end-to-end mission simulator.

preprint2012arXiv

Bayesian Inference for LISA Pathfinder using Markov Chain Monte Carlo Methods

We present a parameter estimation procedure based on a Bayesian framework by applying a Markov Chain Monte Carlo algorithm to the calibration of the dynamical parameters of a space based gravitational wave detector. The method is based on the Metropolis-Hastings algorithm and a two-stage annealing treatment in order to ensure an effective exploration of the parameter space at the beginning of the chain. We compare two versions of the algorithm with an application to a LISA Pathfinder data analysis problem. The two algorithms share the same heating strategy but with one moving in coordinate directions using proposals from a multivariate Gaussian distribution, while the other uses the natural logarithm of some parameters and proposes jumps in the eigen-space of the Fisher Information matrix. The algorithm proposing jumps in the eigen-space of the Fisher Information matrix demonstrates a higher acceptance rate and a slightly better convergence towards the equilibrium parameter distributions in the application to LISA Pathfinder data . For this experiment, we return parameter values that are all within $\sim1σ$ of the injected values. When we analyse the accuracy of our parameter estimation in terms of the effect they have on the force-per-unit test mass noise estimate, we find that the induced errors are three orders of magnitude less than the expected experimental uncertainty in the power spectral density.

preprint2011arXiv

Analysis of the residual force noise for the LISA Technology Package

The analysis of the noise sources perturbing a test mass (TM) geodesic motion is the main scientific objective of the LISA Technology Package experiment (LTP) on board of the LISA Pathfinder space mission. Information on force noise acting on TMs are obtained with a data reduction procedure involving system parameters. Such parameters can be estimated from dedicated experimental runs. Therefore the final estimation of force noise is affected by two sources of uncertainty. One is statistical and connected to the random nature of noisy signals. The other is connected to the uncertainties on the system parameters. The analysis of simulated LTP data is indicating that the major contribution to the force noise power spectral density uncertainties is coming from the statistical properties of the spectrum estimator.

preprint2010arXiv

Calibrating spectral estimation for the LISA Technology Package with multichannel synthetic noise generation

The scientific objectives of the Lisa Technology Package (LTP) experiment, on board of the LISA Pathfinder mission, demand for an accurate calibration and validation of the data analysis tools in advance of the mission launch. The levels of confidence required on the mission outcomes can be reached only with an intense activity on synthetically generated data. A flexible procedure allowing the generation of cross-correlated stationary noise time series was set-up. Multi-channel time series with the desired cross correlation behavior can be generated once a model for a multichannel cross-spectral matrix is provided. The core of the procedure is the synthesis of a noise coloring multichannel filter through a frequency-by-frequency eigendecomposition of the model cross-spectral matrix and a Z-domain fit. The common problem of initial transients in noise time series is solved with a proper initialization of the filter recursive equations. The noise generator performances were tested in a two dimensional case study of the LTP dynamics along the two principal channels of the sensing interferometer.

preprint2009arXiv

Discrete derivative estimation in LISA Pathfinder data reduction

Data analysis for the LISA Technology package (LTP) experiment to be flown aboard the LISA Pathfinder mission requires the solution of the system dynamics for the calculation of the force acting on the test masses (TMs) starting from interferometer position data. The need for a solution to this problem has prompted us to implement a discrete time domain derivative estimator suited for the LTP experiment requirements. We first report on the mathematical procedures for the definition of two methods; the first based on a parabolic fit approximation and the second based on a Taylor series expansion. These two methods are then generalized and incorporated in a more general class of five point discrete derivative estimators. The same procedure employed for the second derivative can be applied to the estimation of the first derivative and of a data smoother allowing defining a class of simple five points estimators for both. The performances of three particular realization of the five point second derivative estimator are analyzed with simulated noisy data. This analysis pointed out that those estimators introducing large amount of high frequency noise can determine systematic errors in the estimation of low frequencies noise levels.