The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

Weak lensing, which is the deflection of light by matter along the line of sight, has proven to be an efficient method for constraining models of structure formation and reveal the nature of dark energy. So far, most weak-lensing studies have focused on the shear field that can be measured directly from the ellipticity of background galaxies. However, within the context of forthcoming full-sky weak-lensing surveys such as Euclid, convergence maps (mass maps) offer an important advantage over shear fields in terms of cosmological exploitation. While it carry the same information, the lensing signal is more compressed in the convergence maps than in the shear field. This simplifies otherwise computationally expensive analyses, for instance, non-Gaussianity studies. However, the inversion of the non-local shear field requires accurate control of systematic effects caused by holes in the data field, field borders, shape noise, and the fact that the shear is not a direct observable (reduced shear). We present the two mass-inversion methods that are included in the official Euclid data-processing pipeline: the standard Kaiser & Squires method (KS), and a new mass-inversion method (KS+) that aims to reduce the information loss during the mass inversion. This new method is based on the KS method and includes corrections for mass-mapping systematic effects. The results of the KS+ method are compared to the original implementation of the KS method in its simplest form, using the Euclid Flagship mock galaxy catalogue. In particular, we estimate the quality of the reconstruction by comparing the two-point correlation functions and third- and fourth- order moments obtained from shear and convergence maps, and we analyse each systematic effect independently and simultaneously...

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

In modern weak-lensing surveys, the common approach to correct for residual systematic biases in the shear is to calibrate shape measurement algorithms using simulations. These simulations must fully capture the complexity of the observations to avoid introducing any additional bias. In this paper we study the importance of faint galaxies below the observational detection limit of a survey. We simulate simplified Euclid VIS images including and excluding this faint population, and measure the shift in the multiplicative shear bias between the two sets of simulations. We measure the shear with three different algorithms: a moment-based approach, model fitting, and machine learning. We find that for all methods, a spatially uniform random distribution of faint galaxies introduces a shear multiplicative bias of the order of a few times $10^{-3}$. This value increases to the order of $10^{-2}$ when including the clustering of the faint galaxies, as measured in the Hubble Space Telescope Ultra-Deep Field. The magnification of the faint background galaxies due to the brighter galaxies along the line of sight is found to have a negligible impact on the multiplicative bias. We conclude that the undetected galaxies must be included in the calibration simulations with proper clustering properties down to magnitude 28 in order to reach a residual uncertainty on the multiplicative shear bias calibration of a few times $10^{-4}$, in line with the $2\times10^{-3}$ total accuracy budget required by the scientific objectives of the Euclid survey. We propose two complementary methods for including faint galaxy clustering in the calibration simulations.

We present the first joint inference of standard sirens and gravitational wave weak lensing by filtering of the same dataset. We imagine a post-LISA scenario emerging around the late 2030s when LISA will have accumulated a number of detections at high redshift; LIGO-VIRGO will have finished observing at low redshift, and Einstein Telescope will have started making new observations out to redshifts possibly overlapping with LISA. Euclid and other cosmological probes will have provided constraints at percent level by then, but mostly exhausted their ability to improve any further. We derive forecasts assuming $\sim1\,\text{deg}^{-2}$ detected sources, in conjunction with a spectroscopic follow-up (e.g. by Euclid, DESI, or ATHENA). Thanks to the statistical power of standard sirens as a geometry probe -- lifting key degeneracies in the gravitational wave weak lensing -- and no external priors assumed, the constraints on dark matter and its clustering, namely $\Omega_m$ and $\sigma_8$, could be achieved to $2\%$ and $3\%$. The Hubble constant could be constrained to better than $1\%$ in all cases; the dark energy density, $\Omega_\Lambda$, to $2\%$, and curvature, $\Omega_K$, to $0.02$; the amplitude and spectral tilt of the scalar fluctuations, $\ln(10^{10}A_s)$ and $n_s$, to $2\%$ and $7\%$. As a completely independent cosmological probe, with less calibration requirements, the joint inference of standard sirens and gravitational wave weak lensing might help solve the tensions currently observed between other cosmological probes, such as CMB, galaxy lensing, and Type Ia SNs, and distinguish between residual systematics and new physics.

With the first two detections in late 2015, astrophysics has officially entered into the new era of gravitational wave observations. Since then, much has been going on in the field with a lot of work focussing on the observations and implications for astrophysics and tests of general relativity in the strong regime. However much less is understood about how gravitational detectors really work at their fundamental level. For decades, the response to incoming signals has been customarily calculated using the very same physical principle, which has proved so successful in the first detections. In this paper we review the physical principle that is behind such a detection at the very fundamental level, and we try to highlight the peculiar subtleties that make it so hard in practice. We will then mention how detectors are built starting from this fundamental measurement element.

General relativity is supported by great experimental evidence. Yet there is a lot of interest in precisely setting its limits with on going and future experiments. A question to answer is about the validity of the Strong Equivalence Principle. Ground experiments and Lunar Laser Ranging have provided the best upper limit on the Nordtvedt parameter $\sigma[\eta]=4.4\times 10^{-4}$. With the future planetary mission BepiColombo, this parameter will be further improved by at least an order of magnitude. In this paper we envisage yet another possible testing environment with spacecraft ranging towards the nearby Sun-Earth collinear Lagrangian points. Neglecting errors in planetary masses and ephemerides, we forecast $\sigma[\eta]=6.4\,(2.0)\times10^{-4}$ (5 yr integration time) via ranging towards $L_1$ in a realistic (optimistic) scenario depending on current (future) range capabilities and knowledge of the Earth's ephemerides. A combined measurement, $L_1$+$L_2$, gives instead $4.8\,(1.7)\times10^{-4}$. In the optimistic scenario a single measurement of one year would be enough to reach $\approx3\times10^{-4}$. All figures are comparable with Lunar Laser Ranging, but worse than BepiColombo. Performances could be much improved if data were integrated over time and over the number of satellites flying around either of the two Lagrangian points. We point out that some systematics (gravitational perturbations of other planets or figure effects) are much more in control compared to other experiments. We do not advocate a specific mission to constrain the Strong Equivalence Principle, but we do suggest analysing ranging data of present and future spacecrafts flying around $L_1$/$L_2$ (one key mission is, for instance, LISA Pathfinder). This spacecraft ranging would be a new and complementary probe to constrain the Strong Equivalence Principle in space.

The measurement of frequency shifts for light beams exchanged between two test masses nearly in free fall is at the heart of gravitational wave detection. It is envisaged that the derivative of the frequency shift is in fact limited by differential forces acting on those test masses. We calculate the derivative of the frequency shift with a fully covariant, gauge-independent and coordinate-free method. This method is general and does not require a congruence of nearby beams' null geodesics as done in previous work. We show that the derivative of the parallel transport is the only means by which gravitational effects shows up in the frequency shift. This contribution is given as an integral of the Riemann tensor --the only physical observable of curvature-- along the beam's geodesic. The remaining contributions are: the difference of velocities, the difference of non-gravitational forces, and finally fictitious forces, either locally at the test masses or non-locally integrated along the beam's geodesic. As an application relevant to gravitational wave detection, we work out the frequency shift in the local Lorentz frame of nearby geodesics.

The basic constituent of interferometric gravitational wave detectors -- the test mass to test mass interferometric link -- behaves as a differential dynamometer measuring effective differential forces, comprising an integrated measure of gravity curvature, inertial effects, as well as non-gravitational spurious forces. This last contribution is going to be characterised by the LISA Pathfinder mission, a technology precursor of future space-borne detectors like eLISA. Changing the perspective from displacement to acceleration can benefit the data analysis of LISA Pathfinder and future detectors. The response in differential acceleration to gravitational waves is derived for a space-based detector's interferometric link. The acceleration formalism can also be integrated into time delay interferometry by building up the unequal-arm Michelson differential acceleration combination. The differential acceleration is nominally insensitive to the system free evolution dominating the slow displacement dynamics of low-frequency detectors. Working with acceleration also provides an effective way to subtract measured signals acting as systematics, including the actuation forces. Because of the strong similarity with the equations of motion, the optimal subtraction of systematic signals, known within some amplitude and time shift, with the focus on measuring the noise provides an effective way to solve the problem and marginalise over nuisance parameters. The $\mathcal{F}$-statistic, in widespread use throughout the gravitation waves community, is included in the method and suitably generalised to marginalise over linear parameters and noise at the same time. The method is applied to LPF simulator data and, thanks to its generality, can also be applied to the data reduction and analysis of future gravitational wave detectors.

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.

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.

LISA Pathfinder is a mission planned by the European Space Agency to test the key technologies that will allow the detection of gravitational waves in space. The instrument on-board, the LISA Technology package, will undergo an exhaustive campaign of calibrations and noise characterisation campaigns in order to fully describe the noise model. Data analysis plays an important role in the mission and for that reason the data analysis team has been developing a toolbox which contains all the functionalities required during operations. In this contribution we give an overview of recent activities, focusing on the improvements in the modelling of the instrument and in the data analysis campaigns performed both with real and simulated data.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

The LISA Pathfinder data analysis team has been developing in the last years the infrastructure and methods required to run the mission during flight operations. These are gathered in the LTPDA toolbox, an object oriented MATLAB toolbox that allows all the data analysis functionalities for the mission, while storing the history of all operations performed to the data, thus easing traceability and reproducibility of the analysis. The parameter estimation methods in the toolbox have been applied recently to data sets generated with the OSE (Off-line Simulations Environment), a detailed LISA Pathfinder non-linear simulator that will serve as a reference simulator during mission operations. These operational exercises aim at testing the on-orbit experiments in a realistic environment in terms of software and time constraints. These simulations, so called operational exercises, are the last verification step before translating these experiments into tele-command sequences for the spacecraft, producing therefore very relevant datasets to test our data analysis methods. In this contribution we report the results obtained with three different parameter estimation methods during one of these operational exercises.

LISA is the upcoming space-based Gravitational Wave telescope. LISA Pathfinder, to be launched in the coming years, will prove and verify the detection principle of the fundamental Doppler link of LISA on a flight hardware identical in design to that of LISA. LISA Pathfinder will collect a picture of all noise disturbances possibly affecting LISA, achieving the unprecedented pureness of geodesic motion necessary for the detection of gravitational waves. The first steps of both missions will crucially depend on a very precise calibration of the key system parameters. Moreover, robust parameters estimation is of fundamental importance in the correct assessment of the residual force noise, an essential part of the data processing for LISA. In this paper we present a maximum likelihood parameter estimation technique in time domain being devised for this calibration and show its proficiency on simulated data and validation through Monte Carlo realizations of independent noise runs. We discuss its robustness to non-standard scenarios possibly arising during the real-life mission, as well as its independence to the initial guess and non-gaussianities. Furthermore, we apply the same technique to data produced in mission-like fashion during operational exercises with a realistic simulator provided by ESA.

The precursor ESA mission LISA-Pathfinder, to be flown in 2013, aims at demonstrating the feasibility of the free-fall, necessary for LISA, the upcoming space-born gravitational wave observatory. LISA Technology Package (LTP) is planned to carry out a number of experiments, whose main targets are to identify and measure the disturbances on each test-mass, in order to reach an unprecedented low-level residual force noise. To fulfill this plan, it is then necessary to correctly design, set-up and optimize the experiments to be performed on-flight and do a full system parameter estimation. Here we describe the progress on the non-linear analysis using the methods developed in the framework of the \textit{LTPDA Toolbox}, an object-oriented MATLAB Data Analysis environment: the effort is to identify the critical parameters and remove the degeneracy by properly combining the results of different experiments coming from a closed-loop system like LTP.

This paper presents a quantitative assessment of the performance of the upcoming LISA Pathfinder geodesic explorer mission. The findings are based on the results of extensive ground testing and simulation campaigns using flight hardware and flight control and operations algorithms. The results show that, for the central experiment of measuring the stray differential acceleration between the LISA test masses, LISA Pathfinder will be able to verify the overall acceleration noise to within a factor two of the LISA requirement at 1 mHz and within a factor 6 at 0.1 mHz. We also discuss the key elements of the physical model of disturbances, coming from LISA Pathfinder and ground measurement, that will guarantee the LISA performance.

Binary systems of compact objects are strong emitters of gravitational waves whose amplitude depends on the binary orbital parameters as the component mass, the orbital semi-major axis and eccentricity. Here, in addition to the famous Hulse-Taylor binary system, we have studied the possibility to detect the gravitational wave signal emitted by binary systems at the center of our galaxy. In particular, recent infrared observation of the galactic center have revealed the existence of a cluster of stars each of which appears to orbit the central black hole in $SgrA^*$. For the stars labelled as S2 and S14, we have studied the emitted spectrum of gravitational wave and compare it with the sensitivity threshold of space-based interferometers like Lisa and Astrod. Furthermore, following recent observations, we have considered the possibility that $SgrA^*$ is actually a binary system of massive black holes and calculated the emission spectrum as a function of the system parameters. The diffraction pattern of gravitational waves emitted by a binary system by a cluster of stars has been also analyzed. We remark that this is only a preliminary-theoretical work than can acquire more interest in view of the next-coming gravitational wave astronomy era.

The diffraction effects on gravitational waves propagating through a stellar cluster are analyzed in the relevant approximation of Fresnel diffraction limit. We find that a gravitational wave scintillation effect - similar to the radio source scintillation effect - comes out naturally, implying that the gravitational wave intensity changes in a characteristic way as the observer moves.