Laser Interferometer Space Antenna Research Articles

A stochastic gravitational wave background in LISA from unresolved white dwarf binaries in the large magellanic cloud

Abstract The Laser Interferometer Space Antenna (LISA) is expected to detect a wide variety of gravitational wave sources in the mHz band. Some of these signals will elude individual detection, instead contributing as confusion noise to one of several stochastic gravitational-wave backgrounds (SGWBs) – notably including the ‘Galactic foreground’, a loud signal resulting from the superposition of millions of unresolved double white dwarf binaries (DWDs) in the Milky Way. It is possible that similar, weaker SGWBs will be detectable from other DWD populations in the local universe, including the Large Magellanic Cloud (LMC). We use the Bayesian LISA Inference Package (${\tt BLIP}$) to investigate the possibility of an anisotropic SGWB generated by unresolved DWDs in the LMC. To do so, we compute the LMC SGWB from a realistic DWD population generated via binary population synthesis, simulate four years of time-domain data with BLIP comprised of stochastic contributions from the LMC SGWB and the LISA detector noise, and analyze this data with BLIP’s spherical harmonic anisotropic SGWB search. We also consider the case of spectral separation from the Galactic foreground. We present the results of these analyses and show, for the first time, that the unresolved DWDs in the LMC will comprise a significant SGWB for LISA.

Distinguishability of binary extreme-mass-ratio inspirals in low frequency band

The inspiral of compact stellar objects into massive black holes are one of the main astrophysical sources for the Laser Interferometer Space Antenna (LISA) and Taiji. These extreme-mass-ratio inspirals (EMRIs) have great potential for cosmology and fundamental physics. A binary extreme-mass-ratio inspiral (b-EMRI) describes the case where binary black holes (BBHs) are captured by a supermassive black hole. The b-EMRIs serve as multi-band gravitational wave sources and provide insights into the dynamics of nuclei and tests of general relativity. However, if the b-EMRIs can be distinguished from the normal EMRIs or not is still not clear. In this work, with a few of assumptions, and using the Teukolsky equation, we calculate the approximate gravitational waves of b-EMRIs and assess their detectability by space-based detectors. We also decouple the secondary object information from the Teukolsky equation, enabling us to calculate the energy fluxes and generate the waveforms more conveniently. Variations in the quadrupole of the binary result in small but non-negligible changes in energy fluxes and waveforms, making it possible to distinguish b-EMRI signals with data analysis. This opens up the potential of using b-EMRIs to test gravity theories and for further astrophysical studies.

Measuring source properties and quasinormal mode frequencies of heavy massive black-hole binaries with LISA

The Laser Interferometer Space Antenna (LISA) will be launched in the mid-2030s. It promises to observe the coalescence of massive black-hole (BH) binaries with signal-to-noise ratios (SNRs) reaching thousands. Crucially, it will detect some of these binaries with high SNR in both the inspiral and the merger-ringdown stages. Such signals are ideal for tests of general relativity (GR) using information from the whole waveform. Here, we consider astrophysically motivated binary systems at the high-mass end of the population observable by LISA and simulate their signals using the newly developed multipolar effective-one-body model: pSEOBNRv5HM. The merger-ringdown signal in this model depends on the binary properties (masses and spins) and also on parameters that describe fractional deviations from the GR quasinormal mode (complex) frequencies of the remnant BH. Performing full Bayesian analyses, we assess to which accuracy LISA will be able to constrain deviations from GR in the ringdown signal when using information from the whole signal. We find that these deviations can typically be constrained to within 10% and in the best cases to within 1%. We also show that with this model we can measure the binary masses and spins with great accuracy even for very massive BH systems with low SNR in the inspiral. In particular, individual source-frame masses can typically be constrained to within 10% and as precisely as 1%, and individual spins can typically be constrained to within 0.1 and, in the best cases, to within 0.001. We also probe the accuracy of the SEOBNRv5HM waveform family by performing synthetic injections of GR numerical-relativity waveforms. Using a novel method that we develop here to quantify the impact of systematic errors, we show that, already for sources with SNR O(100), we would measure erroneous deviations from GR due to waveform model inaccuracies. One of the main sources of error is the mismodeling of the relative alignment between harmonics. These results confirm the need for improving waveform models to perform tests of GR with binary BHs observed at high SNR by LISA. Published by the American Physical Society 2024

Gravitational wave peeps from EMRIs and their implication for LISA signal confusion noise

Scattering events around the center of massive galaxies will occasionally toss a stellar-mass compact object into an orbit around the massive black hole (MBH) at the center, beginning an extreme mass ratio inspiral (EMRI). The early stages of such a highly eccentric orbit are not likely to produce detectable gravitational waves (GWs), as the source will only be in a suitable frequency band briefly when it is close to periapsis during each long-period orbit. This repeated burst of emission, firmly in the millihertz band, is the GW peep. While a single peep is not likely to be detectable, if we consider an ensemble of such subthreshold sources, spread across the Universe, together they may produce an unresolvable background noise that could obscure sources otherwise detectable by the Laser Interferometer Space Antenna. Previous studies of the extreme mass ratio signal confusion background focused either on parabolic orbits near the MBH or events closer to merger. We seek to improve this characterization by implementing numerical kludge waveforms that can calculate highly eccentric orbits with relativistic effects. Our focus is on orbits at the point of capture that are farther away from the MBH. Here we present the waveforms and spectra of peeps generated from recent calculations of EMRIs/extreme mass ratio bursts capture parameters and discuss how these can be used to estimate the signal confusion noise generated by such events. We demonstrate the effects of changing the orbital parameters on the resulting spectra as well as showing direct comparisons to parabolic orbits and why the GW ‘peep’ needs to be studied further. The results of this study will be expanded upon in a further paper that aims to provide an update on the EMRI signal confusion noise problem.

Viscous torque in turbulent magnetized active galactic nucleus accretion disks and its effects on the gravitational waves of extreme mass ratio inspirals

The merger of supermassive black holes produces millihertz gravitational waves (GWs), which are potentially detectable by the future Laser Interferometer Space Antenna (LISA). Such binary systems are usually embedded in an accretion disk environment at the center of an active galactic nucleus (AGN). Recent studies suggest the plasma environment imposes measurable imprints on the GW signal if the mass ratio of the binary is around q $ $. The effect of the gaseous environment on the GW signal is strongly dependent on the disk’s parameters; therefore, it is believed that future low-frequency GW detections will provide us with precious information about the physics of AGN accretion disks. We investigated this effect by measuring the viscous torque via modeling of the evolution of magnetized tori around the primary massive black hole. Using the general relativistic magnetohydrodynamic HARM-COOL code, we performed 2D and 3D simulations of weakly magnetized, thin accretion disks, with a possible truncation and transition to advection-dominated accretion flow. We studied the angular momentum transport and turbulence generated by the magnetorotational instability. We quantified the disk’s effective alpha viscosity and its evolution over time. We applied our numerical results to quantify the relativistic viscous torque on a hypothetical low-mass secondary black hole via a 1D analytical approach, and we estimated the GW phase shift due to the gas environment.

The evaluation for plasma noise in arbitrary time-delay interferometry combinations

The laser interferometer space antenna (LISA) uses laser interferometry to measure gravitational wave-induced distance changes between freely falling test masses on separate spacecraft. In practice, the space-borne gravitational wave detector operates in a plasma medium, and subsequently, the variations in electron density affect the refractive index and add displacement noise to measurements. Geometric time-delay interferometry (TDI) is the TDI combinations searched by geometric method, which is employed to mitigate overwhelming laser phase noise. This study explores all forty-five second-generation geometric TDI combinations up to sixteen links, analyzing plasma noise effects via power spectral density and cross-spectral density calculations for different optical links. Our findings reveal that plasma noise can exceed the noise floor, which comprises optical metrology system and test mass acceleration noise, in certain low-frequency regions due to different noise transfer behaviors. Further, we establish electron density requirements for various TDI combinations, showing that densities below 100cm−3/Hz @ 10mHz satisfy LISA’s criteria across all forty-five combinations. This allows for optimized TDI selection based on specific plasma densities.

Study of quantum electron diffraction for the LISA test-mass charging

The future Laser Interferometer Space Antenna (LISA) for gravitational wave detection in space and the precursor LISA Pathfinder missions rely on metal test masses (TMs) playing the role of free-falling geodesic-reference and interferometer end mirrors. Charge deposited on the TMs by high-energy particles of galactic and solar origin couples with stray electric fields thus perturbing the free-falling geodesic motion. In previous works it has been shown that sub-keV electrons play a relevant role for TM charging. In particular, electrons with energies below 100 eV show a wave-like behavior when impacting on the gold coatings of the TMs and electrode-housing (EH), the capacitive position sensor/actuator that surrounds the TMs. We aim at studying the effects of quantum diffraction of electrons on LISA Pathfinder and LISA TMs and EH gold lattice. We calculate the quantum diffraction probability of electrons with energies < 100 eV impinging on gold lattice from different directions by using the methods of quantum mechanics. The effects of electron quantum diffraction on LISA TM charging are reported here for the first time. This work is also of interest for the second-generation LISA-like space interferometers and for all space missions whose performance relies on free-falling metal test masses.

LISA and LISA-like mission test-mass charging for gamma-ray burst detection

Cubic gold-platinum test-masses (TMs) play the role of free-falling geodesic reference and interferometer end mirrors of the future Laser Interferometer Space Antenna (LISA) observatory for low-frequency gravitational wave detection in space. A similar arrangement has been proposed for the Chinese missions Taiji and TianQin and for the LISA follow-on missions such as ALIA and BBO. The TMs are charged by high-energy particles and photons. The deposited charge couples with stray electric fields surrounding the TMs thus inducing spurious forces that limit the sensitivity of the mission mainly at low frequencies. The TM charging was measured in space in 2016-2017 with LISA Pathfinder (LPF), meant for the testing of the LISA instrumentation. Unfortunately, during the time LPF remained in orbit, no solar energetic particle events or major astrophysical phenomena were observed. We aim at estimating the LISA TM charging attributable to long, short gamma-ray bursts (GRBs) and magnetar flares in comparison to that of charged particles of galactic and solar origin. The contribution of these major astrophysical phenomena to the LISA TM charging is discussed here for the first time. The results found here can be extended to LISA-like missions. The response of the radiation monitors hosted on the three LISA S/C is also evaluated. We show that long, intense extragalactic GRBs and galactic magnetar flares at kpc distances can be detected and monitored through a sudden change from positive to negative charging of the TMs and an increase of the TM charging noise. This is a unique signature since both galactic and solar particles charge positively the LISA TMs. This peculiar behavior of the TM charging would allow monitoring the whole dynamics of GRBs. The suggestion reported in the literature, about the detection of long GRBs and gravitational waves from the same sources, in principle, may apply to LISA and other LISA-like missions since the increase of the TM charging noise during these extreme transient phenomena is estimated to remain below the mission sensitivity while particle detectors are expected to saturate.

A Superconducting Tensor Detector for Mid-Frequency Gravitational Waves: Its Multichannel Nature and Main Astrophysical Targets

Abstract Mid-frequency band gravitational-wave detectors will be complementary to the existing Earth-based detectors (sensitive above 10 Hz or so) and the future space-based detectors such as the Laser Interferometer Space Antenna (LISA), which will be sensitive below around 10 mHz. A ground-based superconducting omnidirectional gravitational radiation observatory (SOGRO) has recently been proposed along with several design variations for the frequency band of 0.1–10 Hz. For two conceptual designs of SOGRO (i.e. SOGRO and advanced SOGRO [aSOGRO]), we examine their multichannel natures, sensitivities, and science cases. One of the key characteristics of the SOGRO concept is its six detection channels. The response functions of each channel are calculated for all possible gravitational wave (GW) polarizations including scalar and vector modes. Combining these response functions, we also confirm the omnidirectional nature of SOGRO. Hence, even a single SOGRO detector will be able to determine the position of a source and polarizations of GWs, if detected. Taking into account SOGRO’s sensitivity and technical requirements, two main targets are most plausible: GWs from compact binaries and stochastic backgrounds. Based on assumptions we consider in this work, detection rates for intermediate-mass binary black holes (in the mass range of hundreds up to $10^{5}\, M_\odot$) are expected to be 0.0065–8.1 yr−1. In order to detect the stochastic GW background, multiple detectors are required. Two aSOGRO detector networks may be able to put limits on the stochastic background beyond the indirect limit from cosmological observations.

Bridging machine learning and diagnostics of the ESA LISA space mission with equation discovery via explainable artificial intelligence

The Laser Interferometer Space Antenna (LISA) of the European Space Agency will be the first interferometer for low-frequency (10-4–10-1 Hz) gravitational wave detection in space. LISA will detect spurious accelerations of the test masses, the end mirrors of the interferometer, of the order of femto-g. Amongst spurious signals, Coulomb forces due to stray electric fields coupling with charge deposited on the test masses by galactic and solar particles and the noise associated to the charging process must be monitored during the mission lifetime. Precious clues on spurious forces acting on the test masses have been studied with the LISA Pathfinder mission, meant for the testing of the instrumentation that will be placed on board LISA, in 2016–2017. In this work we present the design and implementation of a workflow leading to equation discovery for the relationship between solar wind speed, galactic cosmic-ray variations and interplanetary magnetic field intensity observations for the diagnostics of LISA. The workflow exploits explainable artificial intelligence tools to build a bridge between the opaque predictions obtained with machine learning (ML) models and data analysis performed by humans with space mission observations. The core step of the workflow is the implementation, tuning and optimisation of an opaque ML regressor explicitly designed for the future LISA mission, based on input observations of the galactic cosmic-ray flux and of the interplanetary magnetic field intensity. The regressor is aimed at reconstructing the solar wind speed, a parameter of fundamental importance, but not available, for the mission environment monitoring. The workflow ends with the application of explainable clustering techniques to the opaque model in order to (i) obtain human-interpretable outputs instead of opaque ones without any noticeable loss in the predictive performance and (ii) discover an equation describing the quantitative relationship between the involved variables, currently missing in the literature. The correlation amongst the transit of interplanetary structures, galactic cosmic-ray variations and LISA test-mass charging is illustrated here.

Probing scalar induced gravitational waves with PTA and LISA: the importance of third order correction

We revisit the calculation of third order scalar induced gravitational waves (SIGWs) and extend it from a monochromatic primordial power spectrum to a more general log-normal one. We investigate the impact of third order SIGWs on signal-to-noise ratio (SNR) of Laser Interferometer Space Antenna (LISA) and pulsar timing array (PTA) observations, and find that third order SIGWs significantly contribute to the total energy density spectrum of gravitational waves (GW) in high-frequency region. For a primordial power spectrum amplitude of Aζ = 10-2 ∼ 10-1, the effects of third order SIGWs lead to a 40% to 400% increase in the SNR for LISA. Additionally, our PTA data analysis reveals that third order SIGWs diminish both the amplitude Aζ and the peak frequency f * of the primordial power spectrum.

Prospects of detecting deviations to Kerr geometry with radiation reaction effects in EMRIs

Direct detection of gravitational waves and binary black hole mergers have proven to be remarkable investigations of general relativity. In order to have a definitive answer as to whether the black hole spacetime under test is the Kerr or non-Kerr, one requires accurate mapping of the metric. Since EMRIs are perfect candidates for space-based detectors, Laser Interferometer Space Antenna (LISA) observations will serve a crucial purpose in mapping the spacetime metric. In this article, we consider such a study with the Johannsen spacetime that captures the deviations from the Kerr black hole and further discuss their detection prospects. We analytically derive the leading order post-Newtonian corrections in the average loss of energy and angular momentum fluxes generated by a stellar-mass object exhibiting eccentric equatorial motion in the Johannsen background. We further study the orbital evolution of the inspiralling object within the adiabatic approximation. We lastly provide the possible detectability of deviations from the Kerr black hole by estimating gravitational wave dephasing and highlight the crucial role of LISA observations.

Probing the vector charge of Sagittarius A* with pulsar timing

Timing a pulsar orbiting around Sagittarius A* (Sgr A*) can provide uswith a unique opportunity of testing gravity theories. We investigate thedetectability of a vector charge carried by the Sgr A* black hole (BH) in thebumblebee gravity model with simulated future pulsar timing observations. Thespacetime of a bumblebee BH introduces characteristic changes to the orbitaldynamics of the pulsar and the light propagation of radio signals. Assuming atiming precision of 1 ms, our simulation shows that a 5-yr observation of apulsar with an orbital period Pb ∼ 0.5 yr and an orbitaleccentricity e ∼ 0.8 can probe a vector charge-to-mass ratio as small asQ/M ∼ 10-3, which is much more stringent than the current constraint fromthe Event Horizon Telescope (EHT) observations, and comparable to theprospective constraint from extreme mass-ratio inspirals with the LaserInterferometer Space Antenna (LISA).

Neutron star – white dwarf binaries: probing formation pathways and natal kicks with LISA

ABSTRACT Neutron star–white dwarf (NS + WD) binaries offer a unique opportunity for studying NS-specific phenomena with gravitational waves. In this paper, we employ the binary population synthesis technique to study the Galactic population of NS + WD binaries with the future Laser Interferometer Space Antenna (LISA). We anticipate approximately $\mathcal {O}(10^2)$ detectable NS + WD binaries by LISA, encompassing both circular and eccentric ones formed via different pathways. Despite the challenge of distinguishing these binaries from more prevalent double white dwarfs (especially at frequencies below 2 mHz), we show that their eccentricity and chirp mass distributions may provide avenues to explore the NS natal kicks and common envelope evolution. Additionally, we investigate the spatial distribution of detectable NS + WD binaries relative to the Galactic plane and discuss prospects for identifying electromagnetic counterparts at radio wavelengths. Our results emphasise LISA’s capability to detect and characterize NS + WD binaries and to offer insights into the properties of the underlying population. Our conclusions carry significant implications for shaping LISA data analysis strategies and future data interpretation.

Impacts of gravitational-wave background from supermassive black hole binaries on the detection of compact binaries by LISA* *Supported by the NSFC (12250010, 11991052) and Key Research Program of Frontier Sciences, CAS, (ZDBS-LY-7009)

In the frequency band of the Laser Interferometer Space Antenna (LISA), extensive research has been conducted on the impact of foreground confusion noise generated by galactic binaries within the Milky Way Galaxy. Additionally, recent evidence of a stochastic signal, announced by the NANOGrav, EPTA, PPTA, CPTA, and InPTA, indicates that the stochastic gravitational-wave background (SGWB) generated by supermassive black hole binaries (SMBHBs) can contribute strong background noise within the LISA band. Given the presence of such strong noise, it is expected to have significant impacts on LISA's scientific missions. In this study, we investigate the impacts of the SGWB generated by SMBHBs on the detection of individual massive black hole binaries, verified galactic binaries, and extreme mass ratio inspirals in the context of LISA. We find it essential to resolve and eliminate the excess noise from the SGWB to guarantee the success of LISA's missions.

MEMS miniaturized low-noise magnetic field sensor for the observation of sub-millihertz magnetic fluctuations in space exploration

The objective of this paper is to show that, using magnetic field modulation with a MEMS resonator, it is possible to reduce the noise floor of a commercial Tunneling Magnetic Resistance by one order of magnitude, in the ultra-low frequency range. Low noise at this frequency range is a strong requirement in planetary or space exploration missions like LISA (Laser Interferometer Space Antenna), where the observed gravitational waves will be in the 0.1 mHz to 0.1 Hz range. It will be shown that the noise spectral amplitude after demodulation from 0.1 mHz to 100 mHz is well below the 100 nT/Hz requirement for the future space-borne gravitational wave detector LISA – that we take as reference – being at 20 nT/Hz at 0.1 mHz, and reaching an average 1.5nT/Hz at higher frequencies, from 100 mHz to 1 Hz.

Bilby in space: Bayesian inference for transient gravitational-wave signals observed with LISA

ABSTRACT The Laser Interferometer Space Antenna (LISA) is scheduled to launch in the mid-2030s, and is expected to observe gravitational-wave candidates from massive black hole binary mergers, extreme mass ratio inspirals, and more. Accurately inferring the source properties from the observed gravitational-wave signals is crucial to maximize the scientific return of the LISA mission. bilby, the user-friendly Bayesian inference library, is regularly used for performing gravitational-wave inference on data from existing ground-based gravitational-wave detectors. Given that Bayesian inference with LISA includes additional subtitles and complexities beyond its ground-based counterpart, in this work we introduce bilby_lisa , a python package that extends bilby to perform parameter estimation with LISA. We show that full nested sampling can be performed to accurately infer the properties of LISA sources from transient gravitational-wave signals in (a) zero noise and (b) idealized instrumental noise. By focusing on massive black hole binary mergers, we demonstrate that higher order multipole waveform models can be used to analyse a year’s worth of simulated LISA data, and discuss the computational cost and performance of full nested sampling compared with techniques for optimizing likelihood calculations, such as the heterodyned likelihood.

Constraining mass transfer models with galactic neutron star–white dwarf binaries as gravitational wave sources

ABSTRACT Neutron star–white dwarf (NSWD) binaries are one of the most abundant sources of gravitational waves (GWs) in the Milky Way. These GW sources are the evolutionary products of primordial binaries that experienced many processes of binary interaction. We employ a binary population synthesis method to investigate the properties of Galactic NSWD binaries detectable by the Laser Interferometer Space Antenna (LISA). In this paper, only the NSWD systems with a COWD or ONeWD component are included. We consider various models related to mass-transfer efficiencies during primordial binary evolution, supernova explosion mechanisms at NS formation, common envelope ejection efficiencies, and critical WD masses that determine the stability of mass transfer between WDs and NSs. Based on our calculations, we estimate that tens to hundreds of LISA NSWD binaries exist in the Milky Way. We find that the detection of LISA NSWD binaries is able to provide profound insights into mass-transfer efficiencies during the evolution of primordial binaries and critical WD masses during mass transfer from a WD to an NS.

LISA Galactic Binaries with Astrometry from Gaia DR3

Galactic compact binaries with orbital periods shorter than a few hours emit detectable gravitational waves (GWs) at low frequencies. Their GW signals can be detected with the future Laser Interferometer Space Antenna (LISA). Crucially, they may be useful in the early months of the mission operation in helping to validate LISA's performance in comparison to prelaunch expectations. We present an updated list of 55 candidate LISA-detectable binaries with measured properties, for which we derive distances based on Gaia Data Release 3 astrometry. Based on the known properties from electromagnetic observations, we predict the LISA detectability after 1, 3, 6, and 48 months using Bayesian analysis methods. We distinguish between verification and detectable binaries as being detectable after 3 and 48 months, respectively. We find 18 verification binaries and 22 detectable sources, which triples the number of known LISA binaries over the last few years. These include detached double white dwarfs, AM CVn binaries, one ultracompact X-ray binary, and two hot subdwarf binaries. We find that across this sample the GW amplitude is expected to be measured to ≈10% on average, while the inclination is expected to be determined with ≈15° precision. For detectable binaries, these average errors increase to ≈50% and ≈40°, respectively.

Low Frequency Hierarchical Cooperative Impulse Control for Gravitational Wave Detector Formation Keeping

Configuration stability is crucial for the detection accuracy of space-based gravitational wave (GW) detectors. This paper proposes a low-frequency hierarchical cooperative impulse control strategy and a high-fidelity orbit optimization method to achieve a formation of GW detectors with high detection performance. The control strategy employs a hierarchical structure to divide the configuration control into position and shape layers. In the upper layer, the position layer, a state transition tensor method is employed to analyze the sensitivity and feasible range of the formation center position on the configuration stability. In the lower layer, the shape layer, a novel consistency protocol cooperatively uses the states of all individuals in the formation to determine the control expectations. The constraint on formation position is relaxed to reduce the number of impulses, achieving a low-frequency impulse control. The optimization method, characterized by a careful selection of the optimization variables and a design of adaptive model continuation from a low-fidelity model to a high-fidelity model, is proposed to reduce the problem-solving difficulty. The proposed strategy and method are applied to the Laser Interferometer Space Antenna (LISA) mission, and the results demonstrate that the configuration stability is effectively enhanced with several low-frequency impulses.