Black Hole Binaries Research Articles

Investigating the Cosmological Rate of Compact Object Mergers from Isolated Massive Binary Stars

Gravitational-wave (GW) detectors are observing compact object mergers from increasingly far distances, revealing the redshift evolution of the binary black hole (BBH)—and soon the black hole–neutron star (BHNS) and binary neutron star (BNS)—merger rate. To help interpret these observations, we investigate the expected redshift evolution of the compact object merger rate from the isolated binary evolution channel. We present a publicly available catalog of compact object mergers and their accompanying cosmological merger rates from population synthesis simulations conducted with the COMPAS software. To explore the impact of uncertainties in stellar and binary evolution, our simulations use two-parameter grids of binary evolution models that vary the common-envelope efficiency with mass transfer accretion efficiency and supernova (SN) remnant mass prescription with SN natal kick velocity, respectively. We quantify the redshift evolution of our simulated merger rates using the local (z ∼ 0) rate, the redshift at which the merger rate peaks, and the normalized differential rates (as a proxy for slope). We find that although the local rates span a range of ∼103 across our model variations, their redshift evolutions are remarkably similar for BBHs, BHNSs, and BNSs, with differentials typically within a factor 3 and peaks of z ≈ 1.2–2.4 across models. Furthermore, several trends in our simulated rates are correlated with the model parameters we explore. We conclude that future observations of the redshift evolution of the compact object merger rate can help constrain binary models for stellar evolution and GW formation channels.

Searching for Gravitational Wave Optical Counterparts with the Zwicky Transient Facility: Summary of O4a

During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network, the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star–black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (False Alarm Rate less than 1 yr−1) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing, Fritz, an instance of SkyPortal, instrumental in coordinating and managing our telescope scheduling, candidate vetting, and follow-up observations through a user-friendly interface. ZTF covered a total of 2841 deg2 within the skymaps of the high-significance GW events, reaching a median depth of g ≈ 20.2 mag. We circulated 15 candidates, but found no viable KN counterpart to any of the GW events. Based on the ZTF non-detections of the high-significance events in O4a, we used a Bayesian approach, nimbus, to quantify the posterior probability of KN model parameters that are consistent with our non-detections. Our analysis favors KNe with initial absolute magnitude fainter than −16 mag. The joint posterior probability of a GW170817-like KN associated with all our O4a follow-ups was 64%. Additionally, we use a survey simulation software, simsurvey, to determine that our combined filtered efficiency to detect a GW170817-like KN is 36%, when considering the 5 confirmed astrophysical events in O3 (1 BNS and 4 NSBH events), along with our O4a follow-ups. Following Kasliwal et al., we derived joint constraints on the underlying KN luminosity function based on our O3 and O4a follow-ups, determining that no more than 76% of KNe fading at 1 mag day−1 can peak at a magnitude brighter than −17.5 mag.

The Potential for Long-lived Intermediate-mass Black Hole Binaries in the Lowest Density Dwarf Galaxies

Intermediate-mass black hole (IMBH) mergers with masses 104–106 M ⊙ are expected to produce gravitational waves detectable by the Laser Interferometer Space Antenna (LISA) with high signal-to-noise ratios from the present day to cosmic dawn. IMBH mergers are expected to take place within dwarf galaxies; however, the dynamics, timescales, and effect on their hosts are largely unexplored. In a previous study, we examined how IMBHs would pair and merge within nucleated dwarf galaxies. IMBHs in nucleated hosts evolve very efficiently, forming a binary system and coalescing within a few hundred million years. Although the fraction of dwarf galaxies (107 M ⊙ ≤ M ⋆ ≤ 1010 M ⊙) hosting nuclear star clusters is between 60% and 100%, this fraction drops to 20%–70% for lower-mass dwarfs (M ⋆ ≈ 107 M ⊙), with the largest drop in low-density environments. Here, we extend our previous study by performing direct N-body simulations to explore the dynamics and evolution of IMBHs within nonnucleated dwarf galaxies, under the assumption that IMBHs exist within these dwarfs. To our surprise, none of the IMBHs in our simulation suite merge within a Hubble time, despite many attaining high eccentricities e ∼ 0.7–0.95. We conclude that extremely low stellar density environments in the centers of nonnucleated dwarfs do not provide an ample supply of stars to interact with an IMBH binary, resulting in its stalling, in spite of triaxiality and high eccentricity, common means to drive a binary to coalescence. Our findings underline the importance of considering all detailed host properties to predict IMBH merger rates for LISA.

Dual Active Galactic Nuclei: Precursors of Binary Supermassive Black Hole Formation and Mergers

The presence of dual active galactic nuclei (AGNs) on scales of a few tens of kiloparsecs can be used to study merger-induced accretion on supermassive black holes (SMBHs) and offer insights about SMBH mergers, using dual AGNs as merger precursors. This study uses the Romulus25 cosmological simulation to investigate the properties and evolution of dual AGNs. We first analyze the properties of AGNs (L bol > 1043 erg s−1) and their neighboring SMBHs (any SMBHs closer than 30 pkpc to an AGN) at z ≤ 2. This is our underlying population. We then applied the luminosity threshold of L bol > 1043 erg s−1 to the neighboring SMBHs thereby identifying dual and multiple AGNs. Our findings indicate an increase in the number of both single and dual AGNs from lower to higher redshifts. We also find that the number of dual AGNs with separations of 0.5–4 kpc is twice the number of duals with separations of 4–30 kpc. All dual AGNs in our sample resulted from major mergers. Compared to single AGNs, duals have a lower black hole-to-halo mass ratio. We found that the properties of dual AGN host halos, including halo mass, stellar mass, star formation rate, and gas mass, are generally consistent with those of single AGN halos, albeit tending toward the higher end of their respective property ranges. Our analysis uncovered a diverse array of evolutionary patterns among dual AGNs, including rapidly evolving systems, slower ones, and instances where SMBH mergers are ineffective.

Rates and Beaming Angles of Gamma-Ray Bursts Associated with Compact Binary Coalescences

Some, if not all, binary neutron star (BNS) coalescences, and a fraction of neutron star–black hole (NSBH) mergers, are thought to produce sufficient mass ejection to power gamma-ray bursts (GRBs). However, this fraction, as well as the distribution of beaming angles of BNS-associated GRBs, is poorly constrained from observation. Recent work applied machine learning tools to analyze GRB light curves observed by Fermi/Gamma-Ray Burst Monitor (GBM) and Swift/Burst Alert Telescope (BAT). GRBs were segregated into multiple distinct clusters, with the tantalizing possibility that one of them (BNS cluster) could be associated with BNSs and another (NSBH cluster) with NSBHs. As a proof of principle, assuming that all GRBs detected by Fermi/GBM and Swift/BAT associated with BNSs (NSBHs) lie in the BNS (NSBH) cluster, we estimate their rates (Gpc−3 yr−1). We compare these rates with corresponding BNS and NSBH rates estimated by the LIGO–Virgo–KAGRA (LVK) collaboration from the first three observing runs (O1, O2, O3). We find that the BNS rates are consistent with LVK’s rate estimates, assuming a uniform distribution of beaming fractions (f b ∈ [0.01, 0.1]). Conversely, using the LVK’s BNS rate estimates, assuming all BNS mergers produce GRBs, we are able to constrain the beaming angle distribution to θ j ∈ [0.°8, 33.°5] at 90% confidence. We similarly place limits on the fraction of GRB-bright NSBHs as f B ∈ [1.3%, 63%] (f B ∈ [0.4%, 15%]) with Fermi/GBM (Swift/BAT) data.

Black holes - from relativity to Hawking radiation

Abstract. Black Hole (BH for short) is a kind of compact body in the universe predicted by the general theory of relativity. The gravitational pull of a black hole is extremely strong. A black hole is a celestial body whose curvature of space-time is so great that light cannot escape from its event horizon. There are four main ideas about how black holes form: collapsing stars, exploding supernovae, merging multiple black holes, and primordial black holes. In this article, we will start with Einstein's theory of relativity to analyze black holes, and focus on Schwarzschild black holes.

Joint Observations of Late Universe Probes: Cosmological Parameter Constraints from Gravitational Wave and Type Ia Supernova Data

The growing tensions between the early Universe and the late Universe increasingly highlight the importance of developing precise probes for late cosmology. As significant late-Universe probes, Type Ia supernovae (SNe Ia) and gravitational waves (GWs) can provide measurements of relative and absolute distances, respectively. Their complementary nature is likely to break the degeneracies among cosmological parameters, thereby yielding more precise constraints. In this study, we use 43 gravitational-wave sources from the Third LIGO–Virgo–KAGRA Gravitational-Wave Transient Catalog (GWTC-3) and 1590 light curves from Pantheon+ compilation to constrain the dark energy models, as an attempt to achieve precise late-Universe cosmological constraints. For the dark siren GW event, we estimate the corresponding redshift using the binary black hole redshift distribution model. The combination of GW and SNe Ia data could provide the precision on the Hubble constant H 0 and the present matter density Ω m of approximately 16% and 8% for the ΛCDM model. If we consider the equation of state of dark energy w, the combination sample constrains the precision of w to approximately 24%. Although the combination of GWs and SNe Ia observations effectively breaks degeneracies among various cosmological parameters, yielding more stringent constraints, the precision of these constraints still does not meet the stringent standards required by precision cosmology. However, it is reasonable to anticipate that, in the near future, the joint observations of GWs and SNe Ia will become a powerful tool, particularly in the late Universe, for the precise measurement of cosmological parameters.

The Detection of Possible Quasiperiodic Oscillations in the BL Lac 4FGL J2139.4−4235

We present periodicity search analyses on the long-term γ-ray light curve of the BL Lacertae object 4FGL J2139.4−4235 observed by the Fermi Large Area Telescope, over a period of more than 15 yr, from 2008 August 4 to 2023 December 10. To determine the quasiperiodic oscillation (QPO) behavior of 4FGL J2139.4−4235 in the 0.3–300 GeV energy range, we used four methods, namely the Lomb–Scargle periodogram, the weighted wavelet z-transform, the phase dispersion minimization, and the autoregressive integrated moving average model. A Monte Carlo simulation technique is used to evaluate the significance level of the QPO signal. Significant levels above 3.5σ were detected in the γ-ray light curve at about 650 days QPO, which is presented throughout the observation period. Interestingly, there was some correlation between the three bands in the discrete correlation function method calculations, which may be an indication that the variability trends between the three bands are similar. We explore the possible physical models and show that a supermassive binary black hole system or a jet helical motion model seem to be reasonable explanations for the potential QPO behavior.

The gravitational redshift effect of quantum matter waves passing a binary black hole

The gravitational redshift effect of quantum matter waves passing a binary black hole

Physical Models for the Astrophysical Population of Black Holes: Application to the Bump in the Mass Distribution of Gravitational-wave Sources

Abstract Gravitational-wave observations of binary black holes have revealed unexpected structure in the black hole mass distribution. Previous studies employ physically motivated phenomenological models and infer the parameters that control the features of the mass distribution that are allowed in their model, associating the constraints on those parameters with their physical motivations a posteriori. In this work, we take an alternative approach in which we introduce a model parameterizing the underlying stellar and core-collapse physics and obtaining the remnant black hole distribution as a derived by-product. In doing so, we constrain the stellar physics necessary to explain the astrophysical distribution of black hole properties under a given model. We apply this to the mapping between initial mass and remnant black hole mass, accounting for mass-dependent mass loss using a simple parameterized description. Allowing the parameters of the initial mass–remnant mass relationship to evolve with redshift permits correlated and physically reasonable changes to features in the mass function. We find that the current data are consistent with no redshift evolution in the core–remnant mass relationship, but place only weak constraints on the change of these parameters. This procedure can be applied to modeling any physical process underlying the astrophysical distribution. We illustrate this by applying our model to the pulsational pair instability supernova (PPISN) process, previously proposed as an explanation for the observed excess of black holes at ∼35 M ⊙. Placing constraints on the reaction rates necessary to explain the PPISN parameters, we concur with previous results in the literature that the peak observed at ∼35 M ⊙ is unlikely to be a signature from the PPISN process as presently understood.

Efficient Stochastic Template Bank Using Inner Product Inequalities

Gravitational wave searches are crucial for studying compact sources such as neutron stars and black holes. Many sensitive modeled searches use matched filtering to compare gravitational strain data to a set of waveform models known as template banks. We introduce a new stochastic placement method for constructing template banks, offering efficiency and flexibility to handle arbitrary parameter spaces, including orbital eccentricity, tidal deformability, and other extrinsic parameters. This method can be computationally limited by the ability to compare proposal templates with the accepted templates in the bank. To alleviate this computational load, we introduce the use of inner product inequalities to reduce the number of required comparisons. We also introduce a novel application of Gaussian Kernel Density Estimation to enhance waveform coverage in sparser regions. Our approach has been employed to search for eccentric binary neutron stars, low-mass neutron stars, primordial black holes, and supermassive black hole binaries. We demonstrate that our method produces self-consistent banks that recover the required minimum fraction of signals. For common parameter spaces, our method shows comparable computational performance and similar template bank sizes to geometric placement methods and stochastic methods, while easily extending to higher-dimensional problems. The time to run a search exceeds the time to generate the bank by a factor of O(105) for dedicated template banks, such as geometric, mass-only stochastic, and aligned spin cases, O(104) for eccentric and O(103) for the tidal deformable bank. With the advent of efficient template bank generation, the primary area for improvement is developing more efficient search methodologies.

Occurrence of Gravitational Collapse in the Accreting Neutron Stars of Binary-driven Hypernovae

The binary-driven hypernova (BdHN) model proposes long gamma-ray bursts (GRBs) originate in binaries composed of a carbon–oxygen (CO) star and a neutron star (NS) companion. The CO core collapse generates a newborn NS and a supernova that triggers the GRB by accreting onto the NSs, rapidly transferring mass and angular momentum to them. This article aims to determine the conditions under which a black hole (BH) forms from NS collapse induced by the accretion and the impact on the GRB’s observational properties and taxonomy. We perform three-dimensional, smoothed particle hydrodynamics simulations of BdHNe using up-to-date NS nuclear equations of state, with and without hyperons, and calculate the structure evolution in full general relativity. We assess the binary parameters leading either NS in the binary to the critical mass for gravitational collapse into a BH and its occurrence time, t col. We include a nonzero angular momentum of the NSs and find that t col ranges from a few tens of seconds to hours for decreasing NS initial angular momentum values. BdHNe I are the most compact (about 5 minute orbital period), promptly form a BH, and release ≳1052 erg of energy. They form NS–BH binaries with tens of kiloyears merger timescales by gravitational-wave emission. BdHNe II and III do not form BHs, and release ∼1050–1052 erg and ≲1050 erg of energy, respectively. They form NS–NS binaries with a range of merger timescales larger than for NS–BH binaries. In some compact BdHNe II, either NS can become supramassive, i.e., above the critical mass of a nonrotating NS. Magnetic braking by a 1013 G field can delay BH formation, leading to BH–BH or NS–BH with tens of kiloyears merger timescales.

Quasiperiodic γ-Ray Modulations in the Blazars PKS 2155-83 and PKS 2255-282

While there has been an increase in interest in the possibility of quasiperiodic oscillations (QPOs) in blazars, the search has hitherto been restricted to sources with well-sampled light curves. Objects with light curves that include gaps have been, to our knowledge, overlooked. Here, we study two such curves, which have the interesting feature of pertaining to relatively high-redshift blazars—FSRQs, PKS 2155-83, and PKS 2255-282—observed by the Fermi Large Area Telescope. Their redshifts border the “cosmic noon” era of galaxy formation and merging, and their light curves exhibit a distinctive pattern of repetitive high and low (gap dominant) states for 15.6 yr. To accommodate for the gaps in the curves, data are integrated over extended time intervals of 1 month and 2 months. The resulting curves were also examined using methods suitable for sparsely sampled data. This investigation of PKS 2155-83 and PKS 2255-282 suggests QPOs with periods of 4.69 ± 0.79 yr (3σ) and 6.82 ± 2.25 yr (2.8σ), respectively. The probability density functions of the blazars’ fluxes, along with the correlation between their flux and spectral index, were also analyzed. Given the epochs in which the objects are observed, the plausibility of a binary black hole scenario as an origin of the apparent periodicity was examined. We estimated the prospective parameters of such a system using a simple geometric model. The total masses were estimated and found to be consistent, in principle, with independent (dynamical) measurements of the central black hole masses in the two host galaxies.

Detecting Detached Black Hole Binaries through Photometric Variability

Understanding the connection between the properties of black holes (BHs) and their progenitors is interesting in many branches of astrophysics. Discovering BHs in detached orbits with luminous companions (LCs) promises to help establish this connection since the LC and BH progenitor are expected to have the same metallicity and formation time. We explore the possibility of detecting BH–LC binaries in detached orbits using photometric variations of the LC flux, induced by tidal ellipsoidal variation, relativistic beaming, and self-lensing. We create realistic present-day populations of detached BH–LC binaries in the Milky Way (MW) using binary population synthesis where we adopt observationally motivated initial stellar and binary properties, star formation history, and the present-day distribution of these sources in the MW based on detailed cosmological simulations. We test detectability of these sources via photometric variability by Gaia and TESS missions by incorporating their respective detailed detection biases as well as interstellar extinction. We find that Gaia is expected to resolve ∼300–1000 (∼700–1500) detached BH–LC binaries with signal-to-noise ratio (SNR) ≥ 10 (1) depending on the photometric precision and details of supernova physics. Similarly, the number of resolved BH–LC binaries with TESS is ∼50–200 (∼140–350). We find that 136−15+15 BH–LC binaries would be common between Gaia and TESS. Moreover, ∼60–70 (∼50–200) BH–LC binaries identifiable using photometry with SNR ≥ 10 may also be resolved using Gaia’s radial velocity (astrometry).

Stripped Helium Star and Compact Object Binaries in Coeval Populations: Predictions Based on Detailed Binary Evolution Models

Massive stars mainly form in close binaries, where their mutual interactions can profoundly alter their evolutionary paths. Evolved binaries consisting of a massive OB-type main-sequence star with a stripped helium star or a compact companion represent a crucial stage in the evolution toward double compact objects, whose mergers are (potentially) detectable via gravitational waves. The recent detection of X-ray-quiet OB+black hole binaries and OB+stripped helium star binaries has set the stage for discovering more of these systems in the near future. In this work, based on 3670 detailed binary-evolution models and using empirical distributions of initial binary parameters, we compute the expected population of such evolved massive binaries in coeval stellar populations, including stars in star clusters and in galaxies with starburst activities, for ages up to 100 Myr. Our results are vividly illustrated in an animation that shows the evolution of these binaries in the color–magnitude diagram over time. We find that the number of OB+black hole binaries peaks around 10 Myr, and OB+neutron star binaries are most abundant at approximately 20 Myr. Both black holes and neutron stars can potentially be found in populations with ages up to 90 Myr. Additionally, we analyze the properties of such binaries at specific ages. We find that OB+helium stars and OB+black hole binaries are likely to be identifiable as single-lined spectroscopic binaries. Our research serves as a guide for future observational efforts to discover such binaries in young star clusters and starburst environments.

Hydrodynamic 3D Simulation of Roche Lobe Overflow in High-mass X-Ray Binaries

Abstract While binary merger events have been an active area of study in both simulations and observational work, the formation channels by which a high-mass star extends from Roche lobe overflow (RLO) in a decaying orbit of a black-hole (BH) companion to a binary black-hole (BBH) system merits further investigation. Variable length-scales must be employed to accurately represent the dynamical fluid transfer and morphological development of the primary star as it conforms to a diminishing Roche lobe under the runaway influence of the proximal BH. We have simulated and evolved binary mass flow under these conditions to better identify the key transitional processes from RLO to BBHs. We demonstrate a new methodology to model RLO systems to unprecedented resolution simultaneously across the envelope, donor wind, tidal stream, and accretion disk regimes without reliance upon previously universal symmetry, mass flux, and angular momentum flux assumptions. We have applied this method to the semidetached high-mass X-ray binary M33 X-7 in order to provide a direct comparison to recent observations of an RLO candidate system at two overflow states of overfilling factors f = 1.01 and f = 1.1. We found extreme overflow (f = 1.1) to be entirely conservative in both mass and angular momentum transport, forming a conical L1 tidal stream of density and deflected angle comparable to existing predictions. This case lies within the unstable mass transfer (MT) regime as recently proposed of M33 X-7. The f = 1.01 case differed in stream geometry, accretion disk size, and efficiency, demonstrating nonconservative stable MT through a ballistic uniform-width stream. The nonconservative and stable nature of the f = 1.01 case MT also suggests that existing assumptions of semidetached binaries undergoing RLO may mischaracterize their role and distribution as progenitors of BBHs and common envelopes.

Circumbinary Disk Spectra Irradiated by Two Central Accretion Disks in a Binary Black Hole System

We study the effect of irradiation from two accretion disks (minidisks) around respective black holes of stellar-to-intermediate masses in a circular binary on the spectrum of a circumbinary disk (CBD) surrounding them. We assume the CBD to be a standard disk and adopt the orbit-averaged irradiation flux because the viscous timescale is much longer than the orbital period. We then solve the energy equation both analytically and numerically to compute the CBD temperature distribution and the corresponding disk spectrum. We find that the analytically calculated spectra are in good agreement with the numerical ones. The CBD spectrum is almost independent of the binary mass ratio. We also find that the combined spectra of two minidisks and the CBD have double peaks, one peak in the soft X-ray band and the other in the infrared (IR) band. The former peak comes from the two minidisks, while the latter peak from the CBD. The observed flux density increases with frequency as ν 1/3 toward the soft X-ray peak, while it decreases with frequency away from the IR peak as ν −5/3. The latter feature is testable with near-IR observations with Subaru and JWST.

Electromagnetic signatures of black hole clusters in the center of super-Eddington galaxies

Supermassive black holes (SMBHs) at the centers of active galaxies are fed by accretion disks that radiate from the infrared or optical to the X-ray bands. Several types of objects can orbit SMBHs, including massive stars, neutron stars, clouds from the broad- and narrow-line regions, and X-ray binaries. Isolated black holes with a stellar origin (BHs of ∼10 M⊙) should also be present in large numbers within the central parsec of the galaxies. These BHs are expected to form a cluster around the SMBH as a result of the enhanced star formation rate in the inner galactic region and the BH migration caused by gravitational dynamical friction. However, except for occasional microlensing effects on background stars or gravitational waves from binary BH mergers, the presence of a BH population is hard to verify. In this paper, we explore the possibility of detecting electromagnetic signatures of a central cluster of BHs when the accretion rate onto the central SMBH is greater than the Eddington rate. In these supercritical systems, the accretion disk launches powerful winds that interact with the objects orbiting the SMBH. Isolated BHs can capture matter from this dense wind, leading to the formation of small accretion disks around them. If jets are produced in these single microquasars, they could be sites of particle acceleration to relativistic energies. These particles in turn are expected to cool by various radiative processes. Therefore, the wind of the SMBH might illuminate the BHs through the production of both thermal and nonthermal radiation.

A dark standard siren measurement of the Hubble constant following LIGO/Virgo/KAGRA O4a and previous runs

Abstract We present a new constraint on the Hubble constant (H0) from the standard dark siren method using a sample of 5 well-covered gravitational waves (GW) alerts reported during the first part of the fourth LIGO/Virgo/KAGRA observing run and with 3 updated standard dark sirens from third observation run in combination the previous constraints from the first three runs. Our methodology relies on the galaxy catalogue method alone. We use a Deep Learning method to derive the full probability density estimation of photometric redshifts using the Legacy Survey catalogues. We add the constraints from well localized Binary Black Hole mergers to the sample of standard dark sirens analysed in our previous work. We combine the H0 posterior for 5 new standard sirens with other 10 previous events (using the most recent available data for the 5 novel events and updated 3 previous posteriors from O3), finding $H_0 = 70.4^{+13.6}_{-11.7}~{\rm km~s^{-1}~Mpc^{-1}}$ (68% Confidence interval) with the catalogue method only. This result represents an improvement of $\sim 23~{{\%}}$ comparing the new 15 dark siren constrain with the previous 10 dark siren constraint and a reduction in uncertainty of $\sim 40~{{\%}}$ from the combination of 15 dark and bright sirens compared with the GW170817 bright siren alone. The combination of dark and bright siren GW170817 with recent jet constraints yields H0 of $68.0^{+4.4}_{-3.8}~{\rm km~s^{-1}~Mpc^{-1}}$, a $\sim 6~{{\%}}$ precision from Standard Sirens, reducing the previous constraint uncertainty by $\sim 10~{{\%}}$ .

Multi-wavelength observations of the luminous fast blue optical transient AT2023fhn. Up to ∼200 days post-explosion

Luminous fast blue optical transients (LFBOTs) are a class of extragalactic transients notable for their rapid rise and fade times, blue colour, and accompanying luminous X-ray and radio emission. Only a handful have been studied in detail since the prototypical example AT\,2018cow. Their origins are currently unknown, but ongoing observations of previous and new events are placing ever stronger constraints on their progenitors. We aim to put further constraints on the LFBOT AT\,2023fhn, and LFBOTs as a class, using information from the multi-wavelength transient light curve, its host galaxy, and local environment. Our primary results were obtained by fitting galaxy models to the spectral energy distribution of AT\,2023fhn's host and local environment, and by modelling the radio light curve of AT\,2023fhn as due to synchrotron self-absorbed emission from an expanding blast wave in the circumstellar medium. We find that neither the host galaxy nor circumstellar environment of AT\,2023fhn are unusual compared with previous LFBOTs, but that AT\,2023fhn has a much lower X-ray to ultraviolet luminosity ratio than previous events. We argue that the variety in ultraviolet-optical to X-ray luminosity ratios among LFBOTs is likely due to viewing angle differences, and that the diffuse, yet young local environment of AT\,2023fhn - combined with a similar circumstellar medium to previous events - favours a progenitor system containing a massive star with strong winds. Plausible progenitor models in this interpretation therefore include the mergers of black holes and Wolf-Rayet stars or failed supernovae.