Search

The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase., Comment: 155 pages, with executive summary and table of contents

Interspacecraft ranging is crucial for the suppression of laser frequency noise via time-delay interferometry (TDI). So far, the effects of on-board delays and ambiguities on the LISA ranging observables were neglected in LISA modelling and data processing investigations. In reality, on-board delays cause offsets and timestamping delays in the LISA measurements, and pseudo-random noise (PRN) ranging is ambiguous, as it only determines the range up to an integer multiple of the PRN code length. In this article, we identify the four LISA ranging observables: PRN ranging, the sideband beatnotes at the interspacecraft interferometer, TDI ranging, and ground-based observations. We derive their observation equations in the presence of on-board delays, noise, and ambiguities. We then propose a three-stage ranging sensor fusion to combine these observables in order to gain accurate and precise ranging estimates. We propose to calibrate the on-board delays on ground and to compensate the associated offsets and timestamping delays in an initial data treatment (stage 1). We identify the ranging-related routines, which need to run continuously during operation (stage 2), and implement them numerically. Essentially, this involves the reduction of ranging noise, for which we develop a Kalman filter combining the PRN ranging and the sideband beatnotes. We further implement crosschecks for the PRN ranging ambiguities and offsets (stage 3). We show that both ground-based observations and TDI ranging can be used to resolve the PRN ranging ambiguities. Moreover, we apply TDI ranging to estimate the PRN ranging offsets.

Finding a stochastic gravitational-wave background (SGWB) of astrophysical or primordial origin is one of the quests of current and future gravitational-wave observatories. While detector networks such as LIGO-Virgo-Kagra or pulsar timing arrays can use cross-correlations to tell instrumental noise and SGWB apart, LISA is likely to be the only flying detector of its kind in 2035. This particularity poses a challenge for data analysis. To tackle it, we present a strategy based on Bayesian model selection. We use a flexible noise power spectral density~(PSD) model and the knowledge of noise and signal transfer functions to allow SGWBs detection when the noise PSD is unknown. With this technique, we then probe the parameter space accessible by LISA for power-law SGWB shapes., Comment: 5 pages, 1 figure, contribution to the 2023 Gravitation session of the 57th Rencontres de Moriond

Time-delay interferometry (TDI) is a crucial step in the on-ground data processing pipeline of the Laser Interferometer Space Antenna (LISA), as it reduces otherwise overwhelming laser noise and allows for the detection of gravitational waves. This being said, several laser noise couplings have been identified that limit the performance of TDI. First, on-board processing, which is used to decimate the sampling rate from tens of MHz down to a few Hz, requires careful design of the anti-aliasing filters to mitigate folding of laser noise power into the observation band. Furthermore, the flatness of those filters is important to limit the effect of the flexing-filtering coupling. Secondly, the post-processing delays applied in TDI are subject to ranging and interpolation errors. All of these effects are partially described in the literature. In this paper, we present them in a unified framework and give a more complete description of aliased laser noise and the coupling of interpolation errors. Furthermore, for the first time, we discuss the impact of laser locking on laser noise residuals in the final TDI output. To verify the validity of the analytic PSD models we derive, we run numerical simulations using LISA Instrument and calculate second-generation TDI variables with PyTDI. We consider a setup with six independent lasers and with locked lasers (locking configuration N1-12). We find that laser locking indeed affects the laser noise residual in the TDI combination as it introduces correlations among the six lasers inducing slight modulations of the PSD compared to the case of six independent lasers. This implies further studies on laser noise residuals should consider the various locking configurations to produce accurate results., Comment: 17 pages, 8 figures

The data produced by the future space-based millihertz gravitational-wave detector LISA will require nontrivial pre-processing, which might affect the science results. It is crucial to demonstrate the feasibility of such processing algorithms and assess their performance and impact on the science. We are building an end-to-end pipeline that includes state-of-the-art simulations and noise reduction algorithms. The simulations must include a detailed model of the full measurement chain, capturing the main features that affect the instrument performance and processing algorithms. In particular, we include in these simulations, for the first time, proper relativistic treatment of reference frames with realistic numerically-optimized orbits; a model for onboard clocks and clock synchronization measurements; proper modeling of total laser frequencies, including laser locking, frequency planning and Doppler shifts; and a better treatment of onboard processing. Using these simulated data, we show that our pipeline is able to reduce the most critical noises and form synchronized observables. By injecting signals from a verification binary, we demonstrate that good parameter estimation can be obtained on this more realistic setup, extending existing results from previous LISA Data Challenges., Comment: Moriond Gravitation 2023 Proceedings, 7 pages, 5 figures

Within its Voyage 2050 planning cycle, the European Space Agency (ESA) is considering long-term large class science mission themes. Gravitational-wave astronomy is among the topics under study. Building on previous work by other authors, this paper studies a gravitational-wave interferometer concept, dubbed "LISAmax", consisting of three spacecraft, each located close to one of the Sun-Earth libration points L3, L4 and L5, forming a triangular constellation with an arm length of 259 million kilometers (to be compared to LISA's 2.5 million kilometer arms). We argue that this is the largest triangular formation that can be reached from Earth without a major leap in mission complexity and cost (hence the name). The sensitivity curve of such a detector is at least two orders of magnitude lower in amplitude than that of LISA, at frequencies below 1 mHz. This makes the observatory sensitive to gravitational waves in the {\mu}Hz range and opens a new window for gravitational-wave astronomy, not covered by any other planned detector concept. We analyze in detail the constellation stability for a 10-year mission in the full numerical model including insertion, dispersion, and self-gravity-induced accelerations. We compute the orbit transfers using a European launcher and chemical propulsion. Different orbit options, such as precessing, inclined orbits, the use of flybys for the transfer, and the launch strategy, are discussed. The payload design parameters are assessed, and the expected sensitivity curve is compared with a number of potential gravitational-wave sources. No show stoppers are identified at this point of the analysis., Comment: 20 pages, 11 figures, 3 tables

Detecting stochastic background radiation of cosmological origin is an exciting possibility for current and future gravitational-wave (GW) detectors. However, distinguishing it from other stochastic processes, such as instrumental noise and astrophysical backgrounds, is challenging. It is even more delicate for the space-based GW observatory LISA since it cannot correlate its observations with other detectors, unlike today's terrestrial network. Nonetheless, with multiple measurements across the constellation and high accuracy in the noise level, detection is still possible. In the context of GW background detection, previous studies have assumed that instrumental noise has a known, possibly parameterized, spectral shape. To make our analysis robust against imperfect knowledge of the instrumental noise, we challenge this crucial assumption and assume that the single-link interferometric noises have an arbitrary and unknown spectrum. We investigate possible ways of separating instrumental and GW contributions by using realistic LISA data simulations with time-varying arms and second-generation time-delay interferometry. By fitting a generic spline model to the interferometer noise and a power-law template to the signal, we can detect GW stochastic backgrounds up to energy density levels comparable with fixed-shape models. We also demonstrate that we can probe a region of the GW background parameter space that today's detectors cannot access., Comment: 15 pages, 10 figures. Added references and minor changes

Motivated by multi-center biomedical studies that cannot share individual data due to privacy and ownership concerns, we develop communication-efficient iterative distributed algorithms for estimation and inference in the high-dimensional sparse Cox proportional hazards model. We demonstrate that our estimator, even with a relatively small number of iterations, achieves the same convergence rate as the ideal full-sample estimator under very mild conditions. To construct confidence intervals for linear combinations of high-dimensional hazard regression coefficients, we introduce a novel debiased method, establish central limit theorems, and provide consistent variance estimators that yield asymptotically valid distributed confidence intervals. In addition, we provide valid and powerful distributed hypothesis tests for any coordinate element based on a decorrelated score test. We allow time-dependent covariates as well as censored survival times. Extensive numerical experiments on both simulated and real data lend further support to our theory and demonstrate that our communication-efficient distributed estimators, confidence intervals, and hypothesis tests improve upon alternative methods.

LISA is an upcoming ESA mission that will detect gravitational waves in space by interferometrically measuring the separation between free-falling test masses at picometer precision. To reach the desired performance, LISA will employ the noise reduction technique time-delay interferometry (TDI), in which multiple raw interferometric readouts are time shifted and combined into the final scientific observables. Evaluating the performance in terms of these TDI variables requires careful tracking of how different noise sources propagate through TDI, as noise correlations might affect the performance in unexpected ways. One example of such potentially correlated noise is the relative intensity noise (RIN) of the six lasers aboard the three LISA satellites, which will couple into the interferometric phase measurements. In this article, we calculate the expected RIN levels based on the current mission architecture and the envisaged mitigation strategies. We find that strict requirements on the technical design reduce the effect from approximately 8.7 pm/rtHz per inter-spacecraft interferometer to that of a much lower sub-1 pm/rtHz noise, with typical characteristics of an uncorrelated readout noise after TDI. Our investigations underline the importance of sufficient balanced detection of the interferometric measurements., Comment: 15 pages, 10 figures, 2 tables

LISA is a space-based mHz gravitational-wave observatory, with a planned launch in 2034. It is expected to be the first detector of its kind, and will present unique challenges in instrumentation and data analysis. An accurate preflight simulation of LISA data is a vital part of the development of both the instrument and the analysis methods. The simulation must include a detailed model of the full measurement and analysis chain, capturing the main features that affect the instrument performance and processing algorithms. Here, we propose a new model that includes, for the first time, proper relativistic treatment of reference frames with realistic orbits; a model for onboard clocks and clock synchronization measurements; proper modeling of total laser frequencies, including laser locking, frequency planning and Doppler shifts; better treatment of onboard processing and updated noise models. We then introduce two implementations of this model, LISANode and LISA Instrument. We demonstrate that TDI processing successfully recovers gravitational-wave signals from the significantly more realistic and complex simulated data. LISANode and LISA Instrument are already widely used by the LISA community and, for example, currently provide the mock data for the LISA Data Challenges., Comment: 29 pages, 16 figures, 4 tables

The LISA mission is the future space-based gravitational wave (GW) observatory of the European Space Agency. It is formed by 3 spacecraft exchanging laser beams in order to form multiple real and virtual interferometers. The data streams to be used in order to extract the large number and variety of GW sources are Time-Delay Interferometry (TDI) data. One important processing to produce these data is the TDI on-ground processing which recombines multiple interferometric on-board measurements to remove certain noise sources from the data such as laser frequency noise or spacecraft jitter. The LISA noise budget is therefore expressed at the TDI level in order to account for the different TDI transfer functions applied for each noise source and thus estimate their real weight on mission performance. In order to derive a usable form of these transfer functions, a model of the beams, the measurements, and TDI have been developed, and several approximation have been made. A methodology for such a derivation has been established, as well as verification procedures. It results in a set of transfer functions, which are now used by the LISA project, in particular in its performance model. Using these transfer functions, realistic noise curves for various instrumental configurations are provided to data analysis algorithms and used for instrument design., Comment: Version with corrections, 15 pages, 7 figures

A prevalent feature of high-dimensional data is the dependence among covariates, and model selection is known to be challenging when covariates are highly correlated. To perform model selection for the high-dimensional Cox proportional hazards model in presence of correlated covariates with factor structure, we propose a new model, Factor-Augmented Regularized Model for Hazard Regression (FarmHazard), which builds upon latent factors that drive covariate dependence and extends Cox's model. This new model generates procedures that operate in two steps by learning factors and idiosyncratic components from high-dimensional covariate vectors and then using them as new predictors. Cox's model is a widely used semi-parametric model for survival analysis, where censored data and time-dependent covariates bring additional technical challenges. We prove model selection consistency and estimation consistency under mild conditions. We also develop a factor-augmented variable screening procedure to deal with strong correlations in ultra-high dimensional problems. Extensive simulations and real data experiments demonstrate that our procedures enjoy good performance and achieve better results on model selection, out-of-sample C-index and screening than alternative methods.

The analysis of gravitational wave (GW) datasets is based on the comparison of measured time series with theoretical templates of the detector's response to a variety of source parameters. For LISA, the main scientific observables will be the so-called time-delay interferometry (TDI) combinations, which suppress the otherwise overwhelming laser noise. Computing the TDI response to GW involves projecting the GW polarizations onto the LISA constellation arms, and then combining projections delayed by a multiple of the light propagation time along the arms. Both computations are difficult to perform efficiently for generic LISA orbits and GW signals. Various approximations are currently used in practice, e.g., assuming constant and equal armlengths, which yields analytical TDI expressions. In this article, we present 'fastlisaresponse', a new efficient GPU-accelerated code that implements the generic TDI response to GWs in the time domain. We use it to characterize the parameter-estimation bias incurred by analyzing loud Galactic-binary signals using the equal-armlength approximation. We conclude that equal-armlength parameter-estimation codes should be upgraded to the generic response if they are to achieve optimal accuracy for high (but reasonable) SNR sources within the actual LISA data., Comment: 15 pages, 7 figures, 2 tables

Time-delay interferometry (TDI) is a data processing technique for LISA designed to suppress the otherwise overwhelming laser noise by several orders of magnitude. It is widely believed that TDI can only be applied once all phase or frequency measurements from each spacecraft have been synchronized to a common time frame. We demonstrate analytically, using as an example the commonly-used Michelson combination X, that TDI can be computed using the raw, unsynchronized data, thereby avoiding the need for an initial synchronization processing step and significantly simplifying the initial noise reduction pipeline. Furthermore, the raw data is free of any potential artifacts introduced by clock synchronization and reference frame transformation algorithms, which allows to operate directly on the MHz beatnotes. As a consequence, in-band clock noise is directly suppressed as part of TDI, in contrast to the approach previously proposed in the literature (in which large trends in the beatnotes are removed before the main laser-noise reduction step, and clock noise is suppressed in an extra processing step). We validate our algorithm with full-scale numerical simulations that use LISA Instrument and PyTDI and show that we reach the same performance levels as the previously proposed methods, ultimately limited by the clock sideband stability., Comment: 14 pages, 6 figures

The LISA Pathfinder (LPF) mission succeeded outstandingly in demonstrating key technological aspects of future space-borne gravitational-wave detectors, such as the Laser Interferometer Space Antenna (LISA). Specifically, LPF demonstrated with unprecedented sensitivity the measurement of the relative acceleration of two free-falling cubic test masses. Although most disruptive non-gravitational forces have been identified and their effects mitigated through a series of calibration processes, some faint transient signals of yet unexplained origin remain in the measurements. If they appear in the LISA data, these perturbations (also called glitches) could skew the characterization of gravitational-wave sources or even be confused with gravitational-wave bursts. For the first time, we provide a comprehensive census of LPF transient events. Our analysis is based on a phenomenological shapelet model allowing us to derive simple statistics about the physical features of the glitch population. We then implement a generator of synthetic glitches designed to be used for subsequent LISA studies, and perform a preliminary evaluation of the effect of the glitches on future LISA data analyses., Comment: 14 pages, 12 figures, to be submitted to PRD

Gravitational wave science is a dynamical, fast-expanding research field founded on results, tools and methodologies drawn from different research areas and communities. Early career scientists entering this field must learn and combine knowledge and techniques from a range of disciplines. The Workshop on Gravitational-Wave Astrophysics for Early Career Scientists (GWAECS), held virtually in May 2021, planted the seeds of an interdisciplinary, well-connected and all-inclusive community of early career scientists working on gravitational waves, able to exchange relevant information and ideas, build a healthy professional and international environment, share and learn valuable skills, and ensure that ongoing research efforts are perpetuated and expanded in order to attain the main scientific goals envisioned by the whole community. GWAECS was the first event unifying early career scientists belonging to different communities, historically associated with different large-scale gravitational wave experiments. It provided a broad perspective on the future of gravitational waves, offered training on soft and transferable skills and allowed ample time for informal discussions between early career scientists and well-known research experts. The essence of those activities is summarised and collected in the present document, which presents a recap of each session of the workshop and aims to provide all early career scientists with a long-lasting, useful reference which constitutes the legacy of all the ideas that circulated at GWAECS., Comment: 60 pages, 12 figures; Legacy document of the Workshop on Gravitational Wave Astrophysics for Early Career Scientists (GW{\AE}CS), Lorentz Center Workshop, online, 3-7 May 2021

This work develops central limit theorems for cross-validation and consistent estimators of its asymptotic variance under weak stability conditions on the learning algorithm. Together, these results provide practical, asymptotically-exact confidence intervals for $k$-fold test error and valid, powerful hypothesis tests of whether one learning algorithm has smaller $k$-fold test error than another. These results are also the first of their kind for the popular choice of leave-one-out cross-validation. In our real-data experiments with diverse learning algorithms, the resulting intervals and tests outperform the most popular alternative methods from the literature., Comment: 34th Conference on Neural Information Processing Systems (NeurIPS 2020); 40 pages, 15 figures