LISA Pathfinder (LPF) has been a space-based mission designed to test new technologies that will be required for a gravitational wave observatory in space. Magnetically driven forces play a key role in the instrument sensitivity in the low-frequency regime (mHz and below), the measurement band of interest for a space-based observatory. The magnetic field can couple to the magnetic susceptibility and remanent magnetic moment from the test masses and disturb them from their geodesic movement. LISA Pathfinder carried on-board a dedicated magnetic measurement subsystem with noise levels of 10 $ \rm nT \ Hz^{-1/2}$ from 1 Hz down to 1 mHz. In this paper we report on the magnetic measurements throughout LISA Pathfinder operations. We characterise the magnetic environment within the spacecraft, study the time evolution of the magnetic field and its stability down to 20 $\mu$Hz, where we measure values around 200 $ \rm nT \ Hz^{-1/2}$ and identify two different frequency regimes, one related to the interplanetary magnetic field and the other to the magnetic field originating inside the spacecraft. Finally, we characterise the non-stationary component of the fluctuations of the magnetic field below the mHz and relate them to the dynamics of the solar wind.

The zodiacal dust complex, a population of dust and small particles that pervades the Solar System, provides important insight into the formation and dynamics of planets, comets, asteroids, and other bodies. Here we present a new set of data obtained using a novel technique: direct measurements of momentum transfer to a spacecraft from individual particle impacts. This technique is made possible by the extreme precision of the instruments flown on the LISA Pathfinder spacecraft, a technology demonstrator for a future space-based gravitational wave observatory that operated near the first Sun-Earth Lagrange point from early 2016 through Summer of 2017. Using a simple model of the impacts and knowledge of the control system, we show that it is possible to detect impacts and measure properties such as the transferred momentum (related to the particle's mass and velocity), direction of travel, and location of impact on the spacecraft. In this paper, we present the results of a systematic search for impacts during 4348 hours of Pathfinder data. We report a total of 54 candidates with momenta ranging from 0.2$\,\mu\textrm{Ns}$ to 230$\,\mu\textrm{Ns}$. We furthermore make a comparison of these candidates with models of micrometeoroid populations in the inner solar system including those resulting from Jupiter-family comets, Oort-cloud comets, Hailey-type comets, and Asteroids. We find that our measured population is consistent with a population dominated by Jupiter-family comets with some evidence for a smaller contribution from Hailey-type comets. This is in agreement with consensus models of the zodiacal dust complex in the momentum range sampled by LISA Pathfinder.

We report on the results of the LISA Pathfinder (LPF) free-fall mode experiment, in which the control force needed to compensate the quasistatic differential force acting on two test masses is applied intermittently as a series of "impulse" forces lasting a few seconds and separated by roughly 350 s periods of true free fall. This represents an alternative to the normal LPF mode of operation in which this balancing force is applied continuously, with the advantage that the acceleration noise during free fall is measured in the absence of the actuation force, thus eliminating associated noise and force calibration errors. The differential acceleration noise measurement presented here with the free-fall mode agrees with noise measured with the continuous actuation scheme, representing an important and independent confirmation of the LPF result. An additional measurement with larger actuation forces also shows that the technique can be used to eliminate actuation noise when this is a dominant factor.

LISA Pathfinder (LPF) was a technology pioneering mission designed to test key technologies required for gravitational wave detection in space. In the low frequency regime (milli-Hertz and below), where space-based gravitational wave observatories will operate, temperature fluctuations play a crucial role since they can couple into the interferometric measurement and the test masses' free-fall accuracy in many ways. A dedicated temperature measurement subsystem, with noise levels in 10$\,\mu$K$\,$Hz$^{-1/2}$ down to $1\,$mHz was part of the diagnostics unit on board LPF. In this paper we report on the temperature measurements throughout mission operations, characterize the thermal environment, estimate transfer functions between different locations and report temperature stability (and its time evolution) at frequencies as low as 10$\,\mu$Hz, where typically values around $1\,$K$\,$Hz$^{-1/2}$ were measured.

Non-recurrent short term variations of the galactic cosmic-ray (GCR) flux above 70 MeV n$^{-1}$ were observed between 2016 February 18 and 2017 July 3 aboard the European Space Agency LISA Pathfinder (LPF) mission orbiting around the Lagrange point L1 at 1.5$\times$10$^6$ km from Earth. The energy dependence of three Forbush decreases (FDs) is studied and reported here. A comparison of these observations with others carried out in space down to the energy of a few tens of MeV n$^{-1}$ shows that the same GCR flux parameterization applies to events of different intensity during the main phase. FD observations in L1 with LPF and geomagnetic storm occurrence is also presented. Finally, the characteristics of GCR flux non-recurrent variations (peaks and depressions) of duration $<$ 2 days and their association with interplanetary structures are investigated. It is found that, most likely, plasma compression regions between subsequent corotating high-speed streams cause peaks, while heliospheric current sheet crossing cause the majority of the depressions.

LISA Pathfinder (LPF) was a European Space Agency mission with the aim to test key technologies for future space-borne gravitational-wave observatories like LISA. The main scientific goal of LPF was to demonstrate measurements of differential acceleration between free-falling test masses at the sub-femto-g level, and to understand the residual acceleration in terms of a physical model of stray forces, and displacement readout noise. A key step toward reaching the LPF goals was the correct calibration of the dynamics of LPF, which was a three-body system composed by two test-masses enclosed in a single spacecraft, and subject to control laws for system stability. In this work, we report on the calibration procedures adopted to calculate the residual differential stray force per unit mass acting on the two test-masses in their nominal positions. The physical parameters of the adopted dynamical model are presented, together with their role on LPF performance. The analysis and results of these experiments show that the dynamics of the system was accurately modeled and the dynamical parameters were stationary throughout the mission. Finally, the impact and importance of calibrating system dynamics for future space-based gravitational wave observatories is discussed.

Galactic cosmic-ray (GCR) energy spectra observed in the inner heliosphere are modulated by the solar activity, the solar polarity and structures of solar and interplanetary origin. A high counting rate particle detector (PD) aboard LISA Pathfinder (LPF), meant for subsystems diagnostics, was devoted to the measurement of galactic cosmic-ray and solar energetic particle integral fluxes above 70 MeV n$^{-1}$ up to 6500 counts s$^{-1}$. PD data were gathered with a sampling time of 15 s. Characteristics and energy-dependence of GCR flux recurrent depressions and of a Forbush decrease dated August 2, 2016 are reported here. The capability of interplanetary missions, carrying PDs for instrument performance purposes, in monitoring the passage of interplanetary coronal mass ejections (ICMEs) is also discussed.

Test mass charging caused by cosmic rays will be a significant source of acceleration noise for space-based gravitational wave detectors like LISA. Operating between December 2015 and July 2017, the technology demonstration mission LISA Pathfinder included a bespoke monitor to help characterise the relationship between test mass charging and the local radiation environment. The radiation monitor made in situ measurements of the cosmic ray flux while also providing information about its energy spectrum. We describe the monitor and present measurements which show a gradual 40% increase in count rate coinciding with the declining phase of the solar cycle. Modulations of up to 10% were also observed with periods of 13 and 26 days that are associated with co-rotating interaction regions and heliospheric current sheet crossings. These variations in the flux above the monitor detection threshold (approximately 70 MeV) are shown to be coherent with measurements made by the IREM monitor on-board the Earth orbiting INTEGRAL spacecraft. Finally we use the measured deposited energy spectra, in combination with a GEANT4 model, to estimate the galactic cosmic ray differential energy spectrum over the course of the mission.

We report on electrostatic measurements made on board the European Space Agency mission LISA Pathfinder. Detailed measurements of the charge-induced electrostatic forces exerted on free-falling test masses (TMs) inside the capacitive gravitational reference sensor are the first made in a relevant environment for a space-based gravitational wave detector. Employing a combination of charge control and electric-field compensation, we show that the level of charge-induced acceleration noise on a single TM can be maintained at a level close to 1.0 fm/s^2/sqrt(Hz) across the 0.1-100 mHz frequency band that is crucial to an observatory such as LISA. Using dedicated measurements that detect these effects in the differential acceleration between the two test masses, we resolve the stochastic nature of the TM charge build up due to interplanetary cosmic rays and the TM charge-to-force coupling through stray electric fields in the sensor. All our measurements are in good agreement with predictions based on a relatively simple electrostatic model of the LISA Pathfinder instrument.

Future space borne gravitational wave detectors will require a precise definition of calibration signals to ensure the achievement of their design sensitivity. The careful design of the test signals plays a key role in the correct understanding and characterisation of these instruments. In that sense, methods achieving optimal experiment designs must be considered as complementary to the parameter estimation methods being used to determine the parameters describing the system. The relevance of experiment design is particularly significant for the LISA Pathfinder mission, which will spend most of its operation time performing experiments to characterise key technologies for future space borne gravitational wave observatories. Here we propose a framework to derive the optimal signals ---in terms of minimum parameter uncertainty--- to be injected to these instruments during its calibration phase. We compare our results with an alternative numerical algorithm which achieves an optimal input signal by iteratively improving an initial guess. We show agreement of both approaches when applied to the LISA Pathfinder case.

The Solar System contains a population of dust and small particles originating from asteroids, comets, and other bodies. These particles have been studied using a number of techniques ranging from in-situ satellite detectors to analysis of lunar microcraters to ground-based observations of zodiacal light. In this paper, we describe an approach for using the LISA Pathfinder (LPF) mission as an instrument to detect and characterize the dynamics of dust particles in the vicinity of Earth-Sun L1. Launching in late 2015, LPF is a dedicated technology demonstrator mission that will validate several key technologies for a future space-based gravitational-wave observatory. The primary science instrument aboard LPF is a precision accelerometer which we show will be capable of sensing discrete momentum impulses as small as $4\times 10^{-8}\,\textrm{N}\cdot\textrm{s}$. We then estimate the rate of such impulses resulting from impacts of micrometeoroids based on standard models of the micrometeoroid environment in the inner solar system. We find that LPF may detect dozens to hundreds of individual events corresponding to impacts of particles with masses $> 10^{-9}\,$g during LPF's roughly six-month science operations phase in a $5\times 10^5\,\textrm{km}$ by $8\times 10^5\,\textrm{km}$ Lissajous orbit around L1. In addition, we estimate the ability of LPF to characterize individual impacts by measuring quantities such as total momentum transferred, direction of impact, and location of impact on the spacecraft. Information on flux and direction provided by LPF may provide insight as to the nature and origin of the individual impact and help constrain models of the interplanetary dust complex in general. Additionally, this direct in-situ measurement of micrometeoroid impacts will be valuable to designers of future spacecraft targeting the environment around L1.

Arm-locking is a technique for stabilizing the frequency of a laser in an inter-spacecraft interferometer by using the spacecraft separation as the frequency reference. A candidate technique for future space-based gravitational wave detectors such as the Laser Interferometer Space Antenna (LISA), arm-locking has been extensive studied in this context through analytic models, time-domain simulations, and hardware-in-the-loop laboratory demonstrations. In this paper we show the Laser Ranging Interferometer instrument flying aboard the upcoming Gravity Recovery and Climate Experiment Follow-On (GRACE-FO) mission provides an appropriate platform for an on-orbit demonstration of the arm-locking technique. We describe an arm-locking controller design for the GRACE-FO system and a series of time-domain simulations that demonstrate its feasibility. We conclude that it is possible to achieve laser frequency noise suppression of roughly two orders of magnitude around a Fourier frequency of 1Hz with conservative margins on the system's stability. We further demonstrate that `pulling' of the master laser frequency due to fluctuating Doppler shifts and lock acquisition transients is less than $100\,$MHz over several GRACE-FO orbits. These findings motivate further study of the implementation of such a demonstration.

The LISA Pathfinder mission will demonstrate the technology of drag-free test masses for use as inertial references in future space-based gravitational wave detectors. To accomplish this, the Pathfinder spacecraft will perform drag-free flight about a test mass while measuring the acceleration of this primary test mass relative to a second reference test mass. Because the reference test mass is contained within the same spacecraft, it is necessary to apply forces on it to maintain its position and attitude relative to the spacecraft. These forces are a potential source of acceleration noise in the LISA Pathfinder system that are not present in the full LISA configuration. While LISA Pathfinder has been designed to meet it's primary mission requirements in the presence of this noise, recent estimates suggest that the on-orbit performance may be limited by this `suspension noise'. The drift-mode or free-flight experiments provide an opportunity to mitigate this noise source and further characterize the underlying disturbances that are of interest to the designers of LISA-like instruments. This article provides a high-level overview of these experiments and the methods under development to analyze the resulting data.

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.

We investigate the capability of LISA to measure the sky position of equal-mass, nonspinning black hole binaries, combining for the first time the entire inspiral-merger-ringdown signal, the effect of the LISA orbits, and the complete three-channel LISA response. We consider an ensemble of systems near the peak of LISA's sensitivity band, with total rest mass of 2\times10^6 M\odot, a redshift of z = 1, and randomly chosen orientations and sky positions. We find median sky localization errors of approximately \sim3 arcminutes. This is comparable to the field of view of powerful electromagnetic telescopes, such as the James Webb Space Telescope, that could be used to search for electromagnetic signals associated with merging massive black holes. We investigate the way in which parameter errors decrease with measurement time, focusing specifically on the additional information provided during the merger-ringdown segment of the signal. We find that this information improves all parameter estimates directly, rather than through diminishing correlations with any subset of well- determined parameters. Although we have employed the baseline LISA design for this study, many of our conclusions regarding the information provided by mergers will be applicable to alternative mission designs as well.

Arm locking is a technique that has been proposed for reducing laser frequency fluctuations in the Laser Interferometer Space Antenna (LISA), a gravitational-wave observatory sensitive in the milliHertz frequency band. Arm locking takes advantage of the geometric stability of the triangular constellation of three spacecraft that comprise LISA to provide a frequency reference with a stability in the LISA measurement band that exceeds that available from a standard reference such as an optical cavity or molecular absorption line. We have implemented a time-domain simulation of arm locking including the expected limiting noise sources (shot noise, clock noise, spacecraft jitter noise, and residual laser frequency noise). The effect of imperfect a priori knowledge of the LISA heterodyne frequencies and the associated 'pulling' of an arm locked laser is included. We find that our implementation meets requirements both on the noise and dynamic range of the laser frequency.

We investigate the precision with which the parameters describing the characteristics and location of nonspinning black hole binaries can be measured with the Laser Interferometer Space Antenna (LISA). By using complete waveforms including the inspiral, merger and ringdown portions of the signals, we find that LISA will have far greater precision than previous estimates for nonspinning mergers that ignored the merger and ringdown. Our analysis covers nonspinning waveforms with moderate mass ratios, q >= 1/10, and total masses 10^5 < M/M_{Sun} < 10^7. We compare the parameter uncertainties using the Fisher matrix formalism, and establish the significance of mass asymmetry and higher-order content to the predicted parameter uncertainties resulting from inclusion of the merger. In real-time observations, the later parts of the signal lead to significant improvements in sky-position precision in the last hours and even the final minutes of observation. For comparable mass systems with total mass M/M_{Sun} = ~10^6, we find that the increased precision resulting from including the merger is comparable to the increase in signal-to-noise ratio. For the most precise systems under investigation, half can be localized to within O(10 arcmin), and 10% can be localized to within O(1 arcmin).

Recent advances in numerical relativity provide a detailed description of the waveforms of coalescing massive black hole binaries (MBHBs), expected to be the strongest detectable LISA sources. We present a preliminary study of LISA's sensitivity to MBHB parameters using a hybrid numerical/analytic waveform for equal-mass, non-spinning holes. The Synthetic LISA software package is used to simulate the instrument response and the Fisher information matrix method is used to estimate errors in the parameters. Initial results indicate that inclusion of the merger signal can significantly improve the precision of some parameter estimates. For example, the median parameter errors for an ensemble of systems with total redshifted mass of one million Solar masses at a redshift of one were found to decrease by a factor of slightly more than two for signals with merger as compared to signals truncated at the Schwarzchild ISCO.