Massive black hole binaries (MBHBs) of $10^5 \, \rm M_\odot - 3 \times 10^7 \, \rm M_\odot $ merging in low redshift galaxies ($z\le4$) are sufficiently loud to be detected weeks before coalescence with the Laser Interferometer Space Antenna (LISA). This allows us to perform the parameter estimation $on$ $the$ $fly$, i.e. as a function of the time to coalescence during the inspiral phase, relevant for early warning of the planned LISA protected periods and for searches of electromagnetic signals. In this work, we study the evolution of the sky position, luminosity distance, chirp mass and mass ratio uncertainties as function of time left before merger. Overall, light systems with total intrinsic mass $\rm M_{\rm tot} = 3 \times 10^5 \, \rm M_\odot$ are characterized by smaller uncertainties than heavy ones ($\rm M_{\rm tot} = 10^7 \, \rm M_\odot$) during the inspiral. Luminosity distance, chirp mass and mass ratio are well constrained at the end of the inspiral. Concerning sky position, at $z=1$, MBHBs with $\rm M_{\rm tot} = 3 \times 10^5 \, \rm M_\odot$ can be localized with a median precision of $\simeq 10^2 \, \rm deg^2 (\simeq 1 \, \rm deg^2)$ at 1 month (1 hour) from merger, while the sky position of heavy MBHBs can be determined to $10 \, \rm deg^2$ only 1 hour before merger. However the uncertainty around the median values broadens with time, ranging in between 0.04 -- 20 $\rm deg^2$ (0.3 -- 3 $\times 10^3 \, \rm deg^2$) for light (heavy) systems at 1 hour before merger. At merger the sky localization improves down to $\simeq 10^{-1} \, \rm deg^2$ for all masses. For the benefit of the observer community, we provide the full set of data from our simulations and simple and ready-to-use analytical fits to describe the time evolution of uncertainties in the aforementioned parameters, valid for systems with total mass between $10^5$--$10^7 \, \rm M_\odot$ and redshift $0.3$--$3$.

One of the central goals of LISA is the detection of gravitational waves from the merger of supermassive black holes. Contrary to stellar-mass black hole mergers, such events are expected to be rich X-ray sources due to the accretion of material from the circumbinary disks onto the black holes. The orbital dynamics before merger is also expected to modulate the resulting X-ray emission via Doppler boosting in a quasi-periodic way, and in a simple phase relation with the gravitational wave from the inspiral of the black holes. Detecting the X-ray source would enable a precise and early localization of the binary, thus allowing many telescopes to observe the very moment of the merger. Although identifying the correct X-ray source in the relatively large LISA sky localization will be challenging due to the presence of many confounding point sources, the quasi-periodic modulation may aid in the identification. We explore the practical feasibility of such idea. We simulate populations of merging supermassive black holes, their detection with LISA and their X-ray lightcurves using a simple model. Taking the parameters of the X-ray Telescope on the proposed NASA Transient Astrophysics Probe, we then design and simulate an observation campaign that searches for the modulated X-ray source while LISA is still observing the inspiral of the black holes. Assuming a fiducial LISA detection rate of $10$ mergers per year at redshift closer than $3.5$, we expect a few detections of modulated X-ray counterparts over the nominal duration of the LISA mission.

In current stellar evolutionary models, the occurrence of pair instability supernovae implies the lack of stellar black holes (BHs) with masses between about $[60, \, 120] \, \rm{M}_\odot$, resulting in the presence of an upper mass gap in the BH mass distribution. In this Letter, we develop a fiducial model describing BHs beyond the pair instability gap, by convolving the initial mass function and star formation rate with the metallicity evolution across cosmic time. Under the ansatz that the underlying physics of binary formation does not change beyond the gap, we then construct the cosmic population of merging BH binaries. The detection rate of BH binaries with both mass components above the gap, is found to range between $\simeq [0.4,\,7] \, \rm{yr}^{-1}$ for LIGO/Virgo at design sensitivity and $[10, \, 460] \, \rm{yr}^{-1}$ for third-generation ground based detectors, considering the most pessimistic and optimistic scenarios. LISA can detect individually these binaries up to thousands of years from coalescence. The number of events merging in less than four years, which enable multi-band observation in sequence, is expected in the range $[1, \, 20]$. While ET will detect all these events, LIGO/Virgo is expected to detect $\lesssim 50\%$ of them. Finally, we estimate that the gravitational-wave background from unresolved sources in the LISA band may be in principle detected with a signal-to-noise ratio between $ \simeq 2.5$ and $\simeq 80$.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

Next-generation observations will revolutionize our understanding of binary black holes and will detect new sources, such as intermediate-mass black holes. Primary science goals include: Discover binary black holes throughout the observable Universe; Reveal the fundamental properties of black holes; Uncover the seeds of supermassive black holes.

The Laser Interferometer Space Antenna (LISA) will observe black hole binaries of stellar origin during their gravitational wave inspiral, months to years before coalescence. Due to the long duration of the signal in the LISA band, a faithful waveform is necessary in order to keep track of the binary phase. This is crucial to extract the signal from the data and to perform an unbiased estimation of the source parameters. We consider Post-Newtonian (PN) waveforms, and analyze the PN order needed to keep the bias caused by the PN approximation negligible relative to the statistical parameter estimation error, as a function of the source parameters. By considering realistic population models, we conclude that for $\sim 90\%$ of the stellar black hole binaries detectable by LISA, waveforms at low Post-Newtonian (PN) order (PN $\le 2$) are sufficiently accurate for an unbiased recovery of the source parameters. Our results provide a first estimate of the trade-off between waveform accuracy and information recovery for this class of LISA sources.