-
Identifying the electromagnetic counterparts of LISA massive black hole binaries in archival LSST data
Authors:
Chengcheng Xin,
Zoltan Haiman
Abstract:
The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will catalogue the light-curves of up to 100 million quasars. Among these there can be up to approximately 100 ultra-compact massive black hole (MBH) binaries, which 5-15 years later can be detected in gravitational waves (GWs) by the Laser Interferometer Space Antenna (LISA). Here we assume that GWs from a MBH binary have been…
▽ More
The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will catalogue the light-curves of up to 100 million quasars. Among these there can be up to approximately 100 ultra-compact massive black hole (MBH) binaries, which 5-15 years later can be detected in gravitational waves (GWs) by the Laser Interferometer Space Antenna (LISA). Here we assume that GWs from a MBH binary have been detected by LISA, and we assess whether or not its electromagnetic (EM) counterpart can be uniquely identified in archival LSST data as a periodic quasar. We use the binary's properties derived from the LISA waveform, such as the past evolution of its orbital frequency, its total mass, distance and sky localization, to predict the redshift, magnitude and historical periodicity of the quasar expected in the archival LSST data. We then use Monte Carlo simulations to compute the false alarm probability, i.e. the number of quasars in the LSST catalogue matching these properties by chance, based on the (extrapolated) quasar luminosity function, the sampling cadence of LSST, and intrinsic ``damped random walk (DRW)" quasar variability. We perform our analysis on four fiducial LISA binaries, with total masses and redshifts of $(M_{\rm bin}/{\rm M_{\odot}},z) = (3\times10^5,0.3)$, $(3\times10^6,0.3)$, $(10^7,0.3)$ and $(10^7,1)$. While DRW noise and aliasing due to LSST's cadence can produce false periodicities by chance, we find that the frequency chirp of the LISA source during the LSST observations washes out these noise peaks and allows the genuine source to stand out in Lomb-Scargle periodograms. We find that all four fiducial binaries yield excellent chances to be uniquely identified, with false alarm probabilities below $10^{-5}$, a week or more before their merger. This then enables deep follow-up EM observations targeting the individual EM counterparts during their inspiral stage.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
Massive Black Hole Binaries as LISA Precursors in the Roman High Latitude Time Domain Survey
Authors:
Zoltán Haiman,
Chengcheng Xin,
Tamara Bogdanović,
Pau Amaro Seoane,
Matteo Bonetti,
J. Andrew Casey-Clyde,
Maria Charisi,
Monica Colpi,
Jordy Davelaar,
Alessandra De Rosa,
Daniel J. D'Orazio,
Kate Futrowsky,
Poshak Gandhi,
Alister W. Graham,
Jenny E. Greene,
Melanie Habouzit,
Daryl Haggard,
Kelly Holley-Bockelmann,
Xin Liu,
Alberto Mangiagli,
Alessandra Mastrobuono-Battisti,
Sean McGee,
Chiara M. F. Mingarelli,
Rodrigo Nemmen,
Antonella Palmese
, et al. (5 additional authors not shown)
Abstract:
With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12…
▽ More
With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12 day orbital periods, which emit copious gravitational radiation and will inevitably merge on timescales of $10^{3-5}$ years. During the last few months of their merger, such binaries are observable with the Laser Interferometer Space Antenna (LISA), a joint ESA/NASA gravitational wave mission set to launch in the mid-2030s. Roman can thus find LISA precursors, provide uniquely robust constraints on the LISA source population, help identify the host galaxies of LISA mergers, and unlock the potential of multi-messenger astrophysics with massive black hole binaries.
△ Less
Submitted 26 June, 2023;
originally announced June 2023.
-
Tidal Peeling Events: low-eccentricity tidal disruption of a star by a stellar-mass black hole
Authors:
Chengcheng Xin,
Zoltan Haiman,
Rosalba Perna,
Yihan Wang,
Taeho Ryu
Abstract:
Close encounters between stellar-mass black holes (BHs) and stars occur frequently in dense star clusters and in the disks of active galactic nuclei (AGNs). Recent studies have shown that in highly eccentric close encounters, the star can be tidally disrupted by the BH (micro-tidal disruption event, or micro-TDE), resulting in rapid mass accretion and possibly bright electromagnetic signatures. He…
▽ More
Close encounters between stellar-mass black holes (BHs) and stars occur frequently in dense star clusters and in the disks of active galactic nuclei (AGNs). Recent studies have shown that in highly eccentric close encounters, the star can be tidally disrupted by the BH (micro-tidal disruption event, or micro-TDE), resulting in rapid mass accretion and possibly bright electromagnetic signatures. Here we consider a scenario in which the star might approach the stellar-mass BH in a gradual, nearly circular inspiral, under the influence of dynamical friction on a circum-binary gas disk or three-body interactions in a star cluster. We perform hydro-dynamical simulations of this scenario using the smoothed particle hydrodynamics code PHANTOM. We find that the mass of the star is slowly stripped away by the BH. We call this gradual tidal disruption a "tidal-peeling event", or a TPE. Depending on the initial distance and eccentricity of the encounter, TPEs might exhibit significant accretion rates and orbital evolution distinct from those of a typical (eccentric) micro-TDE.
△ Less
Submitted 22 March, 2023;
originally announced March 2023.
-
Dissecting the microphysics behind the metallicity-dependence of massive stars radii
Authors:
Chengcheng Xin,
Mathieu Renzo,
Brian D. Metzger
Abstract:
Understanding the radii of massive stars throughout their evolution is important to answering numerous questions about stellar physics, from binary interactions on the main sequence to the pre-supernova radii. One important factor determining a star's radius is the fraction of its mass in elements heavier than Helium (metallicity, $Z$). However, the metallicity enters stellar evolution through sev…
▽ More
Understanding the radii of massive stars throughout their evolution is important to answering numerous questions about stellar physics, from binary interactions on the main sequence to the pre-supernova radii. One important factor determining a star's radius is the fraction of its mass in elements heavier than Helium (metallicity, $Z$). However, the metallicity enters stellar evolution through several distinct microphysical processes, and which dominates can change throughout stellar evolution and with the overall magnitude of $Z$. We perform a series of numerical experiments with 15M$_{\odot}$ MESA models computed doubling separately the metallicity entering the radiative opacity, the equation of state, and the nuclear reaction network to isolate the impact of each on stellar radii. We explore separately models centered around two metallicity values: one near solar $Z=0.02$ and another sub-solar $Z\sim10^{-3}$, and consider several key epochs from the end of the main sequence to core carbon depletion. We find that the metallicity entering the opacity dominates at most epochs for the solar metallicity models, contributing to on average $\sim$60 - 90% of the total change in stellar radius. Nuclear reactions have a larger impact ($\sim$50 - 70%) during most epochs in the subsolar $Z$ models. The methodology introduced here can be employed more generally to propagate known microphysics errors into uncertainties on macrophysical observables including stellar radii.
△ Less
Submitted 24 June, 2022; v1 submitted 22 June, 2022;
originally announced June 2022.
-
Ultra-Short-Period Massive Black Hole Binary Candidates in LSST as LISA "Verification Binaries"
Authors:
Chengcheng Xin,
Zoltan Haiman
Abstract:
The Legacy Survey of Space and Time (LSST) by the Vera C. Rubin Observatory is expected to discover tens of millions of quasars. A significant fraction of these could be powered by coalescing massive black hole (MBH) binaries, since many quasars are believed to be triggered by mergers. We show that under plausible assumptions about the luminosity functions, lifetimes, and binary fractions of quasa…
▽ More
The Legacy Survey of Space and Time (LSST) by the Vera C. Rubin Observatory is expected to discover tens of millions of quasars. A significant fraction of these could be powered by coalescing massive black hole (MBH) binaries, since many quasars are believed to be triggered by mergers. We show that under plausible assumptions about the luminosity functions, lifetimes, and binary fractions of quasars, we expect the full LSST quasar catalogue to contain between 20-100 million compact MBH binaries with masses $M=10^{5-9}M_{\odot}$, redshifts $z=0-6$, and orbital periods $P=1-70$ days. Their light-curves are expected to be distinctly periodic, which can be confidently distinguished from stochastic red-noise variability, because LSST will cover dozens, or even hundreds of cycles. A very small subset of 10-150 ultra-compact ($P\lesssim1$ day) binary quasars among these will, over $\sim$5-15 years, evolve into the mHz gravitational-wave (GW) frequency band and can be detected by $\textit{LISA}$. They can therefore be regarded as "$\textit{LISA}$ verification binaries", analogous to short-period Galactic compact-object binaries. The practical question is how to find these handful of "needles in the haystack" among the large number of quasars: this will likely require a tailored co-adding analysis optimised for this purpose.
△ Less
Submitted 30 April, 2021;
originally announced May 2021.
-
Multimessenger pulsar timing array constraints on supermassive black hole binaries traced by periodic light curves
Authors:
Chengcheng Xin,
Chiara M. F. Mingarelli,
Jeffrey S. Hazboun
Abstract:
Supermassive black hole binary systems (SMBHBs) emitting gravitational waves may be traced by periodic light curves. We assembled a catalog of 149 such periodic light curves, and using their masses, distances, and periods, predicted the gravitational-wave strain and detectability of each binary candidate using all-sky detection maps. We found that the International Pulsar Timing Array (IPTA) provi…
▽ More
Supermassive black hole binary systems (SMBHBs) emitting gravitational waves may be traced by periodic light curves. We assembled a catalog of 149 such periodic light curves, and using their masses, distances, and periods, predicted the gravitational-wave strain and detectability of each binary candidate using all-sky detection maps. We found that the International Pulsar Timing Array (IPTA) provides almost uniform sky coverage -- a unique ability of the IPTA -- and by 2025 will improve NANOGrav's current minimum detectable strain by a factor of 6, and its volume by a factor of 216. Moreover, IPTA will reach detection sensitivities for three candidates by 2025, and 13 by the end of the decade, enabling us to constrain the underlying empirical relations used to estimate SMBH masses. We find that we can in fact already constrain the mass of a binary in Mrk 504 to $M<3.3\times 10^9~M_\odot$. We also identify 24 high-mass high-redshift galaxies which, according to our models, should not be able to host SMBHBs. Importantly the GW detection of even one of these candidates would be an essentially eternal multimessenger system, and identifying common false positive signals from non-detections will be useful to filter the data from future large-scale surveys such as LSST.
△ Less
Submitted 19 May, 2021; v1 submitted 24 September, 2020;
originally announced September 2020.
-
Correlation between Optical and UV Variability of Quasars
Authors:
Chengcheng Xin,
Maria Charisi,
Zoltán Haiman,
David Schiminovich
Abstract:
The variability of quasars across multiple wavelengths is a useful probe of physical conditions in active galactic nuclei. In particular, variable accretion rates, instabilities, and reverberation effects in the accretion disk of a supermassive black hole (SMBH) are expected to produce correlated flux variations in UV and optical bands. Recent work has further argued that binary quasars should exh…
▽ More
The variability of quasars across multiple wavelengths is a useful probe of physical conditions in active galactic nuclei. In particular, variable accretion rates, instabilities, and reverberation effects in the accretion disk of a supermassive black hole (SMBH) are expected to produce correlated flux variations in UV and optical bands. Recent work has further argued that binary quasars should exhibit strongly correlated UV and optical periodicities. Strong UV-optical correlations have indeed been established in small samples of up to approximately 30 quasars with well-sampled light curves, and have extended the "bluer-when-brighter" trend previously found within the optical bands. Here we further test the nature of quasar variability by examining the observed-frame UV-optical correlations in a large sample of 1,315 bright quasars with overlapping UV and optical light curves for the Galaxy Evolution Explorer (GALEX) and the Catalina Real-time Transient Survey (CRTS), respectively. We find that strong correlations exist in this much larger sample, but we rule out, at approximately 95% confidence, the simple hypothesis that the intrinsic UV and optical variations of all quasars are fully correlated. Our results therefore imply the existence of physical mechanism(s) that can generate uncorrelated optical and UV flux variations.
△ Less
Submitted 9 January, 2020;
originally announced January 2020.
-
Testing the relativistic Doppler boost hypothesis for the binary candidate quasar PG1302-102 with multi-band Swift data
Authors:
Chengcheng Xin,
Maria Charisi,
Zoltan Haiman,
David Schiminovich,
Matthew J. Graham,
Daniel Stern,
Daniel J. D'Orazio
Abstract:
The bright quasar PG1302-102 has been identified as a candidate supermassive black hole binary from its near-sinusoidal optical variability. While the significance of its optical periodicity has been debated due to the stochastic variability of quasars, its multi-wavelength variability in the ultraviolet (UV) and optical bands is consistent with relativistic Doppler boost caused by the orbital mot…
▽ More
The bright quasar PG1302-102 has been identified as a candidate supermassive black hole binary from its near-sinusoidal optical variability. While the significance of its optical periodicity has been debated due to the stochastic variability of quasars, its multi-wavelength variability in the ultraviolet (UV) and optical bands is consistent with relativistic Doppler boost caused by the orbital motion in a binary. However, this conclusion was based previously on sparse UV data which were not taken simultaneously with the optical data. Here we report simultaneous follow-up observations of PG1302-102 with the Ultraviolet Optical Telescope on the Neil Gehrels Swift Observatory in six optical + UV bands. The additional nine Swift observations produce light curves roughly consistent with the trend under the Doppler boost hypothesis, which predicts that UV variability should track the optical, but with a ~2.2 times higher amplitude. We perform a statistical analysis to quantitatively test this hypothesis. We find that the data are consistent with the Doppler boost hypothesis when we compare the the amplitudes in optical B-band and UV light curves. However, the ratio of UV to V-band variability is larger than expected and is consistent with the Doppler model, only if either the UV/optical spectral slopes vary, the stochastic variability makes a large contribution in the UV, or the sparse new optical data underestimate the true optical variability. We have evidence for the latter from comparison with the optical light curve from ASAS-SN. Additionally, the simultaneous analysis of all four bands strongly disfavors the Doppler boost model whenever Swift V-band is involved. Additional, simultaneous optical + UV observations tracing out another cycle of the 5.2-year proposed periodicity should lead to a definitive conclusion.
△ Less
Submitted 22 June, 2020; v1 submitted 25 July, 2019;
originally announced July 2019.