Stellar Black Holes Research Articles

Mono- and oligochromatic extreme-mass-ratio inspirals

The gravitational capture of a stellar-mass object by a supermassive black hole represents a unique probe of warped spacetime. The small object, typically a stellar-mass black hole, describes a very large number of cycles before crossing the event horizon. Because of the mass difference, we call these captures extreme-mass-ratio inspirals (EMRIs). Their merger event rate at the Galactic Center is negligible, but the amount of time spent in the early inspiral is not. Early EMRIs (E-EMRIs) spend hundreds of thousands of years in band during this phase. At very early stages, the peak of the frequency will not change during an observational time. At later stages in the evolution, it will change a bit and finally the EMRI explores a wide range of them when it is close to merger. We distinguish between “monocromatic” E-EMRIs, which do not change their (peak) frequency, oligochromatic E-EMRIs, which explore a short range and polychromatic ones, and the EMRIs which have been discussed so far in the literature. We derive the number of E-EMRIs at the Galactic Center, and we also calculate their signal-to-noise ratios (SNR) and perform a study of parameter extraction. We show that parameters such as the spin and the mass can be extracted with an error which can be as small as 10−11 and 10−5M⊙. There are between hundreds and thousands of E-EMRIs in their monochromatic stage at the Galactic Centre (GC), and tens in their oligochromatic phase. The SNR ranges from a minimum of 10 (larger likelihood) to a maximum of 106 (smaller likelihood). Moreover, we derive the contribution signal corresponding to the incoherent sum of continuous sources of oligochromatic E-EMRIs with two representatives masses, 10M⊙ and 40M⊙, and show that their curves will cover a significant part of Laser Interferometer Space Antennas (LISAs) sensitivity curve. Depending on their level of circularization, they might be detected as individual sources or form a foreground population. Published by the American Physical Society 2024

Galactic Stellar Black Hole Binaries: Spin Effects on Jet Emissions of High-Energy Gamma-Rays

In the last few decades, galactic stellar black hole X-ray binary systems (BHXRBs) have aroused intense observational and theoretical research efforts specifically focusing on their multi-messenger emissions (radio waves, X-rays, γ-rays, neutrinos, etc.). In this work, we investigate jet emissions of high-energy neutrinos and gamma-rays created through several hadronic and leptonic processes taking place within the jets. We pay special attention to the effect of the black hole’s spin (Kerr black holes) on the differential fluxes of photons originating from synchrotron emission and inverse Compton scattering and specifically on their absorption due to the accretion disk’s black-body radiation. The black hole’s spin (dimensionless spin parameter a*) enters into the calculations through the radius of the innermost circular orbit around the black hole, the RISCO parameter, assumed to be the inner radius of the accretion disk, which determines its optical depth τdisk. In our results, the differential photon fluxes after the absorption effect are depicted as a function of the photon energy in the range 1GeV ≤E≤103GeV. It is worth noting that when the black holes’ spin (α*) increases, the differential photon flux becomes significantly lower.

Studying Binary Formation under Dynamical Friction Using Hill’s Problem

Using the equations of motion from Hill’s problem, with added accelerations for different forms of dynamical friction, we provide the (to-date) broadest scale-free study of friction-driven binary formation in gaseous disks and stellar clusters. We focus mainly on binary formation between stellar-mass black holes in active galactic nuclei (AGNs), considering both gas dynamical friction (GDF) from AGN disks and stellar dynamical friction (SDF) from the nuclear star cluster. We first find simple, dimensionless friction coefficients that approximate the effects of standard models for GDF and SDF. We perform extensive simulations of Hill’s problem under such friction, and we present a picture of binary formation through encounters between single stars on nearby orbits, as a function of friction parameter, eccentricity, and inclination. Notably, we find that the local binary formation rate is a linear function of the friction coefficient so long as the friction is weak. Due to the dimensionless nature of our model problem, our findings are generalizable to binary formation at all scales (e.g., intermediate-mass black holes in a star cluster, planetesimals in a gaseous disk).

Close Encounters of Wide Binaries Induced by the Galactic Tide: Implications for Stellar Mergers and Gravitational-wave Sources

A substantial fraction of stars can be found in wide binaries with projected separations between ∼102 and 105 au. In the standard lore of binary physics, these would evolve as effectively single stars that remotely orbit one another on stationary Keplerian ellipses. However, embedded in their Galactic environment, the low binding energy of wide binaries makes them exceptionally prone to perturbations from the gravitational potential of the Milky Way and encounters with passing stars. Employing a fully relativistic N-body integration scheme, we study the impact of these perturbations on the orbital evolution of wide binaries along their trajectory through the Milky Way. Our analysis reveals that the torques exerted by the Galaxy can cause large-amplitude oscillations of the binary eccentricity to 1 − e ≲ 10−8. As a consequence, the wide binary members pass close to each other at periapsis, which, depending on the type of binary, potentially leads to a mass transfer or collision of stars or to an inspiral and subsequent merger of compact remnants due to gravitational-wave radiation. Based on a simulation of 105 wide binaries across the Galactic field, we find that this mechanism could significantly contribute to the rate of stellar collisions and binary black hole mergers as inferred from observations of luminous red novae and gravitational-wave events by LIGO/Virgo/Kagra. We conclude that the dynamics of wide binaries, despite their large mean separation, can give rise to extreme interactions between stars and compact remnants.

Accretion and dynamical evolutions of stellar mass black holes in active galactic nucleus disks

Accretion and dynamical evolutions of stellar mass black holes in active galactic nucleus disks

Influence of Black Hole Kick Velocity on Microlensing Distributions

The natal kick velocity distribution for black holes (BHs) is unknown regardless of its importance for understanding the BH formation process. Gravitational microlensing is a unique tool for studying the distribution of BHs in our Galaxy, and the first isolated stellar-mass BH event, OGLE-2011-BLG-0462/MOA-2011-BLG-191 (OB110462), was recently identified by astrometric microlensing. This study investigates how the natal kick velocity for Galactic BHs affects the microlensing event rate distribution. We consider a Maxwell distribution with various average kick velocities, as well as the consequent variation of the spatial distribution of BHs. We find that the event rate for the BH lenses toward the Galactic bulge decreases as v avg increases, mainly due to the scale height inflation. We focus on the unique microlensing parameters measured for OB110462, with microlens parallax π E larger than 0.06 for its long timescale of t E > 200 days. We calculate the expected number of BH events occurring with parameters similar to OB110462 during the OGLE-IV survey by Mróz et al. and compare it with the actual number that occurred, at least one. Our fiducial model predicts 0.52, 0.38, 0.18, 0.042, and 4.0 × 10−3 events occurring for v avg = 25, 50, 100, 200, and 400 km s−1, respectively, which suggests that the average kick velocity is likely to be v avg ≲ 100 km s−1. The expected number smaller than unity even at maximum might indicate our luck in finding OB110462, which can be tested with future surveys by, e.g., the Roman Space Telescope.

Rapid identification of time-frequency domain gravitational wave signals from binary black holes using deep learning* *Supported by the National SKA Program of China (2022SKA0110200, 2022SKA0110203), the National Natural Science Foundation of China (12473001, 11975072, 11875102, 11835009), and the National 111 Project (B16009)

Recent developments in deep learning techniques have provided alternative and complementary approaches to the traditional matched-filtering methods for identifying gravitational wave (GW) signals. The rapid and accurate identification of GW signals is crucial to the advancement of GW physics and multi-messenger astronomy, particularly considering the upcoming fourth and fifth observing runs of LIGO-Virgo-KAGRA. In this study, we used the 2D U-Net algorithm to identify time-frequency domain GW signals from stellar-mass binary black hole (BBH) mergers. We simulated BBH mergers with component masses ranging from 7 to 50 and accounted for the LIGO detector noise. We found that the GW events in the first and second observation runs could all be clearly and rapidly identified. For the third observing run, approximately 80% of the GW events could be identified. In contrast to traditional convolutional neural networks, the U-Net algorithm can output time-frequency domain signal images corresponding to probabilities, providing a more intuitive analysis. In conclusion, the U-Net algorithm can rapidly identify the time-frequency domain GW signals from BBH mergers.

Apparently Ultralong Period Radio Signals from Self-lensed Pulsar–Black Hole Binaries

Pulsar–black hole (BH) close binary systems, which have not been found yet, are unique laboratories for testing theories of gravity and understanding the formation channels of gravitational-wave sources. We study the self-gravitational lensing effect in a pulsar–BH system on the pulsar’s emission. Because this effect occurs once per orbital period for almost edge-on binaries, we find that it could generate apparently ultralong period (minutes to hours) radio signals when the intrinsic pulsar signal is too weak to detect. Each of such lensed signals, or “pulse,” is composed of a number of amplified intrinsic pulsar pulses. We estimate that a radio telescope with a sensitivity of 10 mJy could detect a few systems that emit such signals in our Galaxy. The model is applied to three recently found puzzling long-period radio sources: GLEAM-X J1627, PSR J0901-4046, and GPM J1839-10. To explain their observed signal durations and periods, the masses of their lensing components are estimated as ∼104 M ⊙, ∼4 M ⊙, and 103−6 M ⊙, respectively, with their binary coalescence times ranging from a few tens to thousands of years. However, the implied merger rates (as high as ∼103−4 Myr−1 per galaxy) and the large period decay rates (>10−8 s s−1) tend to disfavor this self-lensing scenario for these three sources. Despite this, our work still provides observational characteristics for self-lensed pulsar–BH binaries, which could help the detection of related sources in the future. Finally, for a binary containing a millisecond pulsar and a stellar-mass BH, the Shapiro delay effect would cause a ≥10% variation of the profile width for the subpulses in such lensed signals.

Radiative plasma simulations of black hole accretion flow coronae in the hard and soft states

Stellar-mass black holes in x-ray binary systems are powered by mass transfer from a companion star. The accreted gas forms an accretion disk around the black hole and emits x-ray radiation in two distinct modes: hard and soft state. The origin of the states is unknown. We perform radiative plasma simulations of the electron-positron-photon corona around the inner accretion flow. Our simulations extend previous efforts by self-consistently including all the prevalent quantum electrodynamic processes. We demonstrate that when the plasma is turbulent, it naturally generates the observed hard-state emission. In addition, we show that when soft x-ray photons irradiate the system—mimicking radiation from an accretion disk—the turbulent plasma transitions into a new equilibrium state that generates the observed soft-state emission. Our findings demonstrate that turbulent motions of magnetized plasma can power black-hole accretion flow coronae and that quantum electrodynamic processes control the underlying state of the plasma.

A Pilot Search for Gravitational Self-Lensing Binaries with the Zwicky Transient Facility

Binary systems containing a compact object may exhibit periodic brightening episodes due to gravitational lensing as the compact object transits the companion star. Such “self-lensing” signatures have been detected before for white dwarf binaries. We attempt to use these signatures to identify detached stellar-mass neutron star and black hole binaries using data from the Zwicky Transient Facility (ZTF). We present a systematic search for self-lensing signals in Galactic binaries from a subset of high-cadence ZTF data taken in 2018. We identify 12 plausible candidates from the search, although because each candidate is observed to only brighten once, other origins such as stellar flares are more likely. We discuss prospects for more comprehensive future searches of the ZTF data.

Disc novae: thermodynamics of gas-assisted binary black hole formation in AGN discs

ABSTRACT We investigate the thermodynamics of close encounters between stellar mass black holes (BHs) in the gaseous discs of active galactic nuclei (AGNs), during which binary black holes (BBHs) may form. We consider a suite of 2D viscous hydrodynamical simulations within a shearing box prescription using the Eulerian grid code athena++. We study formation scenarios where the fluid is either an isothermal gas or an adiabatic mixture of gas and radiation in local thermal equilibrium. We include the effects of viscous and shock heating, as well as optically thick cooling. We co-evolve the embedded BHs with the gas, keeping track of the energetic dissipation and torquing of the BBH by gas and inertial forces. We find that compared to the isothermal case, the minidiscs formed around each BH are significantly hotter and more diffuse, though BBH formation is still efficient. We observe massive blast waves arising from collisions between the radiative minidiscs during both the initial close encounter and subsequent periapsis periods for successfully bound BBHs. These ‘disc novae’ have a profound effect, depleting the BBH Hill sphere of gas and injecting energy into the surrounding medium. In analysing the thermal emission from these events, we observe periodic peaks in local luminosity associated with close encounters/periapses, with emission peaking in the optical/near-infrared (IR). In the AGN outskirts, these outbursts can reach 4 per cent of the AGN luminosity in the IR band, with flares rising over 0.5–1 yr. Collisions in different disc regions, or when treated in 3D with magnetism, may produce more prominent flares.

Deci-Hz gravitational waves from the self-interacting axion cloud around a rotating stellar mass black hole

Deci-Hz gravitational waves from the self-interacting axion cloud around a rotating stellar mass black hole

Concurrent Accretion and Migration of Giant Planets in Their Natal Disks with Consistent Accretion Torque

Migration commonly occurs during the epoch of planet formation. For emerging gas giant planets, it proceeds concurrently with their growth through the accretion of gas from their natal protoplanetary disks. A similar migration process should also be applied to the stellar-mass black holes embedded in active galactic nucleus disks. In this work, we perform high-resolution 3D and 2D numerical hydrodynamical simulations to study the migration dynamics for accreting embedded objects over the disk viscous timescales in a self-consistent manner. We find that an accreting planet embedded in a predominantly viscous disk has a tendency to migrate outward, in contrast to the inward orbital decay of nonaccreting planets. 3D and 2D simulations find the consistent outward migration results for the accreting planets. Under this circumstance, the accreting planet’s outward migration is mainly due to the asymmetric spiral arms feeding from the global disk into the Hill radius. This is analogous to the unsaturated corotation torque although the imbalance is due to material accretion within the libration timescale rather than diffusion onto the inner disk. In a disk with a relatively small viscosity, the accreting planets clear deep gaps near their orbits. The tendency of inward migration is recovered, albeit with suppressed rates. By performing a parameter survey with a range of disks’ viscosity, we find that the transition from outward to inward migration occurs with the effective viscous efficiency factor α ∼ 0.003 for Jupiter-mass planets.

Detecting Gravitational-wave Bursts from Black Hole Binaries in the Galactic Center with LISA

Stellar-mass black hole binaries (BHBs) in galactic nuclei are gravitationally perturbed by the central supermassive black hole (SMBH) of the host galaxy, potentially inducing strong eccentricity oscillations through the eccentric Kozai–Lidov mechanism. These highly eccentric binaries emit a train of gravitational-wave (GW) bursts detectable by the Laser Interferometer Space Antenna (LISA)—a planned space-based GW detector—with signal-to-noise ratios up to ∼100 per burst. In this work, we study the GW signature of BHBs orbiting our galaxy’s SMBH, Sgr A*, which are consequently driven to very high eccentricities. We demonstrate that an unmodeled approach using a wavelet decomposition of the data effectively yields the time-frequency properties of each burst, provided that the GW frequency peaks between 10−3 and 10−1 Hz. The wavelet parameters may be used to infer the eccentricity of the binary, measuring log10(1−e) within an error of 20%. Our proposed search method can thus constrain the parameter space to be sampled by complementary Bayesian inference methods, which use waveform templates or orthogonal wavelets to reconstruct and subtract the signal from LISA data.

Rethinking Thorne–Żytkow Object Formation: Assembly via Common Envelope in Field Binaries

Thorne–Żytkow objects (TŻOs), hypothetical merger products in which a neutron star is embedded in a stellar core, are traditionally considered steady-state configurations. Their assembly, especially through dynamical channels, is not well understood. The predominant focus in the literature has been on the observational signatures related to the evolution and long-term fate of TŻOs, with their initial formation often treated as a given. However, the foundational calculations supporting the existence of TŻOs assume nonrotating spherically symmetric initial conditions that we find to be inconsistent with a binary merger scenario. In this work, we explore the implications of postmerger dynamics in TŻO formation scenarios with field binary progenitors, specifically the role that angular momentum transport during the common envelope phase plays in constraining the possible merger products, using the tools of stellar evolution and three-dimensional hydrodynamics. We also propose an alternative steady-state outcome for these mergers: the thin-envelope TŻO, an equilibrium solution consisting of a low-mass spherical envelope supported by the accretion disk luminosity of a central stellar-mass black hole. These configurations may be of interest to upcoming time-domain surveys as potential X-ray sources that may be preceded by a series of bright transient events.

Whispering in the dark: Faint X-ray emission from black holes with OB star companions

Recently, astrometric and spectroscopic surveys of OB stars revealed a few stellar-mass black holes (BHs) with orbital periods of as low as 10 days. Contrary to wind-fed BH high-mass X-ray binaries (HMXBs), no X-ray counterpart was detected, probably because of the absence of a radiatively efficient accretion disc around the BH. Nevertheless, dissipative processes in the hot, dilute, and strongly magnetised plasma around the BH (so-called BH corona) can still lead to non-thermal X-ray emission (e.g. synchrotron). We determine the X-ray luminosity distribution from BH+OB star binaries up to orbital periods of a few thousand days. odot $ and orbital periods of 1-3000\,d. We computed the X-ray luminosity for a broad range of radiative efficiencies that depend on the mass accretion rate and flow geometry. odot $, which aligns well with our predictions and may be interesting sources for follow-up observations. The predicted luminosities of the OB companions to these X-ray-emitting BHs are 10$^ odot These findings advocate for prolonged X-ray observations of the stellar-mass black hole candidates identified in the vicinity of OB stars. Such long exposures could reveal the underlying population of X-ray-faint BHs and provide constraints for the evolution from single to double degenerate binaries and identify the progenitors of gravitational wave mergers.

How black hole activity may influence exoplanetary evolution in our Galaxy

ABSTRACT An increasing number of exoplanets have been discovered in the Milky Way galaxy, which is also known to harbour a super-massive black hole (Sagittarius A*) at its centre. Here, we investigate how the central black hole (BH) activity may affect the evolution of exoplanets in our Galaxy. Accreting BHs emit high-energy radiation – extreme ultraviolet and X-rays – which can lead to XUV photoevaporation of the planetary atmospheres. We evaluate the atmospheric mass-loss using both theoretical estimates of the BH radiative output and observational constraints on the past activity history of Sgr A*. The resulting mass-loss is analysed as a function of the galactocentric distance. For the first time, we compute the exoplanet atmospheric evolution under BH irradiation by explicitly including the temporal evolution of the central luminosity output (i.e. the BH activity history). We obtain that Sgr A* could have a major impact on exoplanets located in the inner region of the Galaxy (e.g. Galactic bulge); a significant fraction of the atmospheric mass can be removed by BH irradiation; and in extreme cases, the initial atmosphere may be completely stripped away. Such mass-loss can have important consequences on the atmospheric chemistry and potential biological evolution. We discuss the physical implications for planetary habitability, and we also briefly consider the case of stellar-mass BHs. Overall, accreting black holes may play a significant role in the evolution of exoplanets in our Galaxy across cosmic time.

High-energy Neutrinos from Outflows Powered by the Kicked Remnants of Binary Black Hole Mergers in Active Galactic Nucleus Accretion Disks

Merging of stellar-mass binary black holes (BBHs) could take place within the accretion disks of active galactic nuclei (AGN). The resulting BH remnant is likely to accrete the disk gas at a super-Eddington rate, launching a fast, quasi-spherical outflow (wind). Particles will be accelerated by shocks driven by the wind, subsequently interacting with the shocked disk gas or radiation field through hadronic processes and resulting in the production of high-energy neutrinos and potential electromagnetic (EM) emission. This study delves into the intricate evolution of the shock driven by a merged BH wind within an AGN disk. Subsequently, we calculated the production of neutrinos and the expected detection numbers for a single event, along with their contributions to the overall diffuse neutrino background. Our analysis, which considers various scenarios, reveals considerable neutrino production and possible detection by IceCube for nearby events. The contribution of merged BH winds on the diffuse neutrino background is minor due to the low event rate density, but it can be improved to some extent for some optimistic parameters. We also propose that there could be two neutrino/EM bursts, one originating from the premerger BBH wind and the other from the merged BH wind, with the latter typically having a delay to the gravitational wave (GW) event of around tens of days. When combined with the anticipated GWs emitted during the BBH merger, such a system emerges as a promising candidate for joint observations involving neutrinos, GWs, and EM signals.

Radiation hydrodynamical simulations of super-Eddington mass transfer and black hole growth in close binaries

ABSTRACT Radiation-driven outflows play a crucial role in extracting mass and angular momentum from binary systems undergoing rapid mass transfer at super-Eddington rates. To study the mass transfer process from a massive donor star to a stellar-mass black hole (BH), we perform multidimensional radiation-hydrodynamical simulations that follow accretion flows from the first Lagrange point down to about a hundred times the Schwarzschild radius of the accreting BH. Our simulations reveal that rapid mass transfer occurring at over a thousand times the Eddington rate leads to significant mass-loss from the accretion disc via radiation-driven outflows. Consequently, the inflow rates at the innermost radius are regulated by two orders of magnitude smaller than the transfer rates. We find that convective motions within the accretion disc drive outward energy and momentum transport, enhancing the radiation pressure in the outskirts of the disc and ultimately generating large-scale outflows with sufficient energy to leave the binary. Furthermore, we observe strong anisotropy in the outflows, which occur preferentially toward both the closest and furthest points from the donor star. However, when averaged over all directions, the specific angular momentum of the outflows is nearly comparable to the value predicted in the isotropic emission case. Based on our simulation results, we propose a formula that quantifies the mass growth rates on BHs and the mass-loss rates from binaries due to radiation-driven outflows. This formula provides important implications for the binary evolution and the formation of merging binary BHs.

The maximum black hole mass at solar metallicity

We analyse the current knowledge and uncertainties in detailed stellar evolution and wind modelling to evaluate the mass of the most massive stellar black hole (BH) at solar metallicity. Contrary to common expectations that it is the most massive stars that produce the most massive BHs, we find that the maximum $M_ BH Max is found in the canonical intermediate range between ZAMS $ simeq 30 and 50 instead. The prime reason for this seemingly counter-intuitive finding is that very massive stars (VMS) have increasingly high mass-loss rates that lead to substantial mass evaporation before they expire as stars and end as lighter BHs than their canonical O-star counterparts.