Stellar Black Holes Research Articles

Nuclear burning in an accretion flow around a stellar-mass black hole embedded within an AGN disk

Abstract A stellar-mass black hole, embedded within the accretion disk of an active galactic nuclei (AGN), has the potential to accrete gas at a rate that can reach approximately ∼109 times the Eddington limit. This study explores the potential for nuclear burning in the rapidly accreting flow towards this black hole and studies how nucleosynthesis affects metal production. Using numerical methods, we have obtained the disk structure while considering nuclear burning and assessed the stability of the disk. In contrast to gas accretion onto the surface of a neutron star or white dwarf, the disk remains stable against the thermal and secular instabilities because advection cooling offsets the nuclear heating effects. The absence of a solid surface for a black hole prevents excessive mass accumulation in the inner disk region. Notably, nuclear fusion predominantly takes place in the inner disk region, resulting in substantial burning of $\rm ^{12}C$ and $\rm ^{3}He$, particularly for black holes around M = 10 M⊙ with accretion rates exceeding approximately ∼107 times the Eddington rate. The ejection of carbon-depleted gas through outflows can lead to an increase in the mass ratio of oxygen or nitrogen to carbon, which may be reflected in observed line ratios such as $\rm N\, V/C\, IV$ and $\rm O\, IV/C\, IV$. Consequently, these elevated spectral line ratios could be interpreted as indications of super-solar metallicity in the broad line region.

Investigating 39 Galactic Wolf-Rayet stars with VLTI/GRAVITY. Uncovering a long-period binary desert

Wolf-Rayet stars (WRs) represent one of the final evolutionary stages of massive stars and are thought to be the immediate progenitors of stellar-mass black holes. Their multiplicity characteristics form an important anchor point in single and binary population models for predicting gravitational-wave progenitors. Recent spectroscopic campaigns have suggested incompatible multiplicity fractions and period distributions for N- and C-rich Galactic WRs (WNs and WCs) at both short and long orbital periods, in contradiction with evolutionary model predictions. In this work, we employed long-baseline infrared interferometry to investigate the multiplicity of WRs at long periods and explored the nature of their companions. We present a magnitude-limited ($K<9$; $V<14$) survey of 39 Galactic WRs, including 11 WN, 15 WC, and 13 H-rich WN (WNh) stars. We used the $K$-band instrument GRAVITY at the Very Large Telescope Interferometer (VLTI) in Chile. The sensitivity of GRAVITY at spatial scales of sim 1--200 milliarcseconds and flux contrast of $1<!PCT!>$ allowed an exploration of periods in the range $10^ $ d and companions down to sim 5$\,M_ We carried out a companion search for all our targets, with the aim of either finding wide companions or calculating detection limits. We also explored the rich GRAVITY dataset beyond a multiplicity search to look for other interesting properties of the WR sample. We detected wide companions with VLTI/GRAVITY for only four stars in our sample: WR 48, WR 89, WR 93, and WR 115. Combining our results with spectroscopic studies, we arrived at observed multiplicity fractions of $f^ WN obs WC obs and WNh obs The multiplicity fractions and period distributions of WNs and WCs are consistent in our sample. For single WRs, we placed upper limits on the mass of potential companions down to sim 5$\,M_ for WNs and WCs, and sim 7$\,M_ for WNh stars. In addition, we also found other features in the GRAVITY dataset, such as (i) a diffuse extended component contributing significantly to the $K$-band flux in over half the WR sample; (ii) five known spectroscopic binaries resolved in differential phase data, which constitutes an alternative detection method for close binaries; and (iii) spatially resolved winds in four stars: WR 16, WR 31a, WR 78, and WR 110. Our survey reveals a lack of intermediate- (a few hundred days) and long- (a few years to decades) period WR systems. The 200d peak in the period distributions of WR+OB and BH+OB binaries predicted by Case B mass-transfer binary evolution models is not seen in our data. The rich companionship of their O-type progenitors in this separation range suggests that the WR progenitor stars expand and interact with their companions, most likely through unstable mass transfer, resulting in either a short-period system or a merger.

Probing Blackbody Components in Gamma-Ray Bursts from Black Hole Neutrino-dominated Accretion Flows

A stellar-mass black hole (BH) surrounded by a neutrino-dominated accretion flow (NDAF) is generally considered to be the central engine of gamma-ray bursts (GRBs). Neutrinos escaping from the disk will annihilate outside the disk to produce the fireball that could power GRBs with blackbody (BB) components. The initial GRB jet power and fireball launch radius are related to the annihilation luminosity and annihilation height of the NDAFs, respectively. In this paper, we collect seven GRBs with known redshifts and identified BB components to test whether the NDAF model works. We find that, in most cases, the values of the accretion rates and the central BH properties are all in the reasonable range, suggesting that these BB components indeed originate from the neutrino annihilation process.

Astroparticles from X-Ray Binary Coronae

The recent observation of high-energy neutrinos from the Galactic plane implies an abundant population of hadronic cosmic-ray sources in the Milky Way. We explore the role of the coronae of accreting stellar-mass black holes as such astroparticle emitters. We show that the particle acceleration and interaction timescales in the coronal region are tied to the compactness of the X-ray source. Thus, neutrino emission processes may similarly happen in the cores of active galactic nuclei and black hole X-ray binaries (XRBs), despite their drastically different masses and physical sizes. We apply the model to the well-measured XRB Cygnus X-1 and find that the cascaded gamma rays accompanying the neutrino emission naturally explain the GeV emission that only presents during the source’s hard state, while the state-averaged gamma-ray emission explains the LHAASO observation above 20 TeV. We show that XRB coronae could contribute significantly to the Galactic cosmic-ray and Galactic plane neutrino fluxes. Our model predicts variable high-energy neutrino emission from bright Galactic XRBs that may be observed by IceCube and future neutrino observatories.

Migration of Accreting Planets and Black Holes in Disks

Nascent planets are thought to lose angular momentum (AM) to the gaseous protoplanetary disk via gravitational interactions, leading to inward migration. A similar migration process also applies to stellar-mass black holes (BHs) embedded in the disks of active galactic nuclei. However, AM exchange via accretion onto the planet/BH may strongly influence the migration torque. In this study, we perform 2D global hydrodynamic simulations of an accreting planet/BH embedded in a disk, where AM exchange between the planet/BH and disk via gravity, accretion, pressure, and viscosity are considered. When accretion is turned off, we recover the linear estimate for Type I migration torque. However, for the two planet masses we investigated with our accreting simulations, we find outward migration due to the positive AM deposited onto the accreting body by the disk gas. Our simulations achieve the global steady state for the transport of mass and AM: The mass and AM fluxes are constant across the disk except for jumps ( ΔṀ and ΔJ̇ ) at the planet’s location, and the jumps match the accretion rate and torque on the planet. Our findings suggest that caution is needed when applying the standard results of disk migration to accreting planets and BHs.

An X-Ray Shell Reveals the Supernova Explosion for Galactic Microquasar SS 433

How black holes are formed remains an open and fundamental question in astrophysics. Despite theoretical predictions, it lacks observations to understand whether the black hole formation experiences a supernova explosion. Here we report the discovery of an X-ray shell north of the Galactic microquasar SS 433 harboring a stellar-mass black hole spatially associated with radio continuum and polarization emissions and an H i cloud. Its spectrum can be reproduced by a 1 keV underionized plasma, from which the shell is inferred to have been created by a supernova explosion 20–30 kyr ago, and its properties constitute evidence for canonical supernova explosions to create some black holes. Our analysis precludes other possible origins including heated by jets or blown by disk winds. According to the lower mass limit of the compact object in SS 433, we roughly deduced that the progenitor should be more massive than 25 M ⊙. The existence of such a young remnant in SS 433 can also lead to new insights into the supercritical accretion in young microquasars and the γ-ray emission of this system. The fallback ejecta may provide accretion materials within tens of thousands of years, while the shock of the supernova remnant may play a crucial role in the cosmic-ray (re)acceleration.

Probing Populations of Dark Stellar Remnants in the Globular Clusters 47 Tuc and Terzan 5 Using Pulsar Timing

We present a new method to combine multimass equilibrium dynamical models and pulsar timing data to constrain the mass distribution and remnant populations of Milky Way globular clusters (GCs). We first apply this method to 47 Tuc, a cluster for which there exists an abundance of stellar kinematic data and which is also host to a large population of millisecond pulsars. We demonstrate that the pulsar timing data allow us to place strong constraints on the overall mass distribution and remnant populations even without fitting on stellar kinematics. Our models favor a small population of stellar-mass black holes (BHs) in this cluster (with a total mass of 446−72+75M⊙ ), arguing against the need for a large (>2000 M ⊙) central intermediate-mass BH. We then apply the method to Terzan 5, a heavily obscured bulge cluster that hosts the largest population of millisecond pulsars of any Milky Way GC and for which the collection of conventional stellar kinematic data is very limited. We improve existing constraints on the mass distribution and structural parameters of this cluster and place stringent constraints on its black hole content, finding an upper limit on the mass in BHs of ∼4000 M ⊙. This method allows us to probe the central dynamics of GCs even in the absence of stellar kinematic data and can be easily applied to other GCs with pulsar timing data, for which data sets will continue to grow with the next generation of radio telescopes.

Feedback-regulated seed black hole growth in star-forming molecular clouds and galactic nuclei

Context. The detection of supermassive black holes (SMBHs) in high-redshift luminous quasars may require a phase of rapid accretion, and as a precondition, substantial gas influx toward seed black holes (BHs) from kiloparsec or parsec scales. Our previous research demonstrated the plausibility of such gas supply for BH seeds within star-forming giant molecular clouds (GMCs) with high surface density (∼104 M⊙ pc−2), facilitating “hyper-Eddington” accretion via efficient feeding by dense clumps, which are driven by turbulence and stellar feedback. Aims. This article presents an investigation of the impacts of feedback from accreting BHs on this process, including radiation, mechanical jets, and highly relativistic cosmic rays. Methods. We ran a suite of numerical simulations to explore diverse parameter spaces of BH feedback, including the subgrid accretion model, feedback energy efficiency, mass loading factor, and initial metallicity. Results. Using radiative feedback models inferred from the slim disk, we find that hyper-Eddington accretion is still achievable, yielding BH bolometric luminosities of as high as 1041 − 1044 erg/s, depending on the GMC properties and specific feedback model assumed. We find that the maximum possible mass growth of seed BHs (ΔMmax BH) is regulated by the momentum-deposition rate from BH feedback, ṗfeedback/(ṀBHc), which leads to an analytic scaling that agrees well with simulations. This scenario predicts the rapid formation of ∼104 M⊙ intermediate-massive BHs (IMBHs) from stellar-mass BHs within ∼1 Myr. Furthermore, we examine the impacts of subgrid accretion models and how BH feedback may influence star formation within these cloud complexes.

Massive scalar clouds and black hole spacetimes in Gauss-Bonnet gravity

We study static black holes in scalar-Gauss-Bonnet (sGB) gravity with a massive scalar field as an example of higher curvature gravity. The scalar mass introduces an additional scale and leads to a strong suppression of the scalar field beyond its Compton wavelength. We numerically compute sGB black hole spacetimes and scalar configurations and also compare with perturbative results for small couplings, where we focus on a dilatonic coupling function. We analyze the constraints on the parameters from requiring the curvature singularity to be located inside the black hole horizon r_hrh and the relation to the regularity condition for the scalar field. For scalar field masses m r_h \gtrsim 10^{-1}mrh≳10−1, this leads to a new and currently most stringent bound on sGB coupling constant \alphaα of \alpha/r_h^2\sim 10^{-1}α/rh2∼10−1 in the context of stellar mass black holes. Lastly, we look at several properties of the black hole configurations relevant for further work on observational consequences, including the scalar monopole charge, Arnowitt–Deser–Misner mass, curvature invariants and the frequencies of the innermost stable circular orbit and light ring.

Dark lens candidates from Gaia Data Release 3

Gravitational microlensing is a phenomenon that allows us to observe the dark remnants of stellar evolution, even if these bodies are no longer emitting electromagnetic radiation. In particular, it can be useful to observe solitary neutron stars or stellar-mass black holes, providing a unique window through which to understand stellar evolution. Obtaining direct mass measurements with this technique requires precise observations of both the change in brightness and the position of the microlensed star. The European Space Agency's satellite can provide both. Using publicly available data from different surveys, we analysed events published in the Data Release 3 ( ) microlensing catalogue. Here, we describe our selection of candidate dark lenses, where we suspect the lens is a white dwarf (WD), a neutron star (NS), a black hole (BH), or a mass-gap object, with a mass in the range between the heaviest NS and the least massive BH. We estimated the mass of the lenses using information obtained from the best-fitting microlensing models, source star, Galactic model, and the expected parameter distributions. We found eleven candidates for dark remnants: one WDs, three NSs, three mass-gap objects, and four BHs.

Constraints on metric-Palatini gravity from QPO data

AbstractIn this work, we study metric-Palatini gravity extended by the antisymmetric part of the affine curvature. This gravity theory leads to general relativity plus a geometric Proca field. Using our previous construction of its static spherically-symmetric AdS solution (Eur Phys J. C 83(4):318, 2023), we perform a detailed analysis in this work using the observational quasiperiodic oscillations (QPOs) data. To this end, we use the latest data from stellar-mass black hole GRO J1655-40, intermediate-mass black hole in M82-X1, and the super-massive black hole in SgA* (our Milky Way) and perform a Monte-Carlo-Markov-Chain (MCMC) analysis to determine or bound the model parameters. Our results shed light on the allowed ranges of the Proca mass and other parameters. The results imply that our solutions can cover all three astrophysical black holes. Our analysis can also be extended to more general metric-affine gravity theories.

Area of Ignorance in Stellar Physics: Stellar Mass Black Hole Distribution; Lowest Initial Progenitor Mass Limit for Black Hole Evolution

Area of Ignorance in Stellar Physics: Stellar Mass Black Hole Distribution; Lowest Initial Progenitor Mass Limit for Black Hole Evolution

Systematic collapse of the accretion disc across the supermassive black hole population

ABSTRACT The structure of the accretion flow on to supermassive black holes is not well understood. Standard disc models match to zeroth-order in predicting substantial energy dissipation within optically thick material producing a characteristic strong blue/ultraviolet continuum. However, they fail at reproducing more detailed comparisons to the observed spectral shapes along with their observed variability. Based on stellar mass black holes within our Galaxy, accretion discs should undergo a transition into an X-ray hot, radiatively inefficient flow, below a (mass scaled) luminosity of $\sim 0.02\, L_{\rm {Edd}}$. While this has been seen in limited samples of nearby low-luminosity active galactic nuclei (AGN) and few rare changing-look AGN, it is not at all clear whether this transition is present in the wider AGN population across cosmic time. A key issue is the difficulty in disentangling a change in spectral state from increased dust obscuration and/or host galaxy contamination, effectively drowning out the AGN emission. Here, we use the new eROSITA eFEDS Survey to identify unobscured AGN from their X-ray emission, matched to excellent optical imaging from Subaru’s Hyper Suprime-Cam; allowing the subtraction of the host galaxy contamination. The resulting, uncontaminated, AGN spectra reveal a smooth transition from a strongly disc-dominated state in bright AGN, to the collapse of the disc into an inefficient X-ray plasma in the low-luminosity AGN, with the transition occurring at $\sim 0.02\, L_{\rm {Edd}}$; revealing fundamental aspects of accretion physics in AGN.

Revised spin for the black hole in GRS 1716-249 given a new distance determination

GRS 1716--249 is a stellar-mass black hole in a low-mass X-ray binary that underwent a giant outburst in 2016-17. In this paper, we use simultaneous observations from the Hard X-ray Modulation Telescope (Insight-HXMT) and the Nuclear Spectroscopic Telescope Array ( NuSTAR ) to determine its basic parameters. The observations were performed during the softest part of the outburst and the spectra show clear thermal disk emission and reflection features. We fit the X-ray energy spectra using the joint fitting method of the continuum and reflection components with the kerrbb2 relxill model. Since there is a possibility that the distance to this source was previously underestimated, we used the latest distance parameter of 6.9\,kpc in our study, in contrast to previous works, where the distance was set at 2.4\,kpc. Through a spectral fitting of the black hole mass at odot $, we observe a strong dependence of the derived spin on the distance: $a_ $ at an assumed distance of 2.4\,kpc and $ at an assumed distance of 6.9\,kpc, at a confidence level of 90<!PCT!>. When considering the uncertainties in the distance and black hole mass, there will be a wider range of spin with $ < 0.78. The fitting results with the new distance indicate that GRS 1716--249 harbors a moderate spin black hole with an inclined ($i circ $) accretion disk around it. Additionally, we have also found that solely using the method of reflection component fitting, while ignoring the constraints on the spin from the accretion disk component will result in an extremely high spin.

The formation and evolution of dark star clusters – II. The impact of primordial mass segregation

ABSTRACT We investigate the impact of primordial mass segregation on the formation and evolution of dark star clusters (DSCs). Considering a wide range of initial conditions, we conducted N-body simulations of globular clusters (GCs) around the Milky Way. In particular, we assume a canonical initial mass function for all GCs without natal kicks for supernova remnants, namely neutron stars or black holes. Our results demonstrate that clusters with larger degrees of primordial mass segregation reach their DSC phase earlier and spend a larger fraction of their dissolution time in such a phase, compared to clusters without mass segregation. In primordially segregated clusters, the maximum Galactocentric distance that the clusters can have to enter the DSC phase is almost twice that of the clusters without primordial mass segregation. Primordially segregated clusters evolve with a higher number of stellar mass black holes, accelerating energy creation in their central regions and consequently increasing evaporation rates and cluster sizes during dark phases. The simulations reveal that aggregating heavy components at the centre doubles the time spent in the dark phase. Additionally, the study identifies potential links between simulated dark clusters and initial conditions of Milky Way GCs, suggesting some may transition to dark phases before dissolution. Higher primordial mass segregation coefficients amplify the average binary black hole formation rate by 2.5 times, raising higher expectations for gravitational-wave emissions.

A Comparison of the X-Ray Polarimetric Properties of Stellar and Supermassive Black Holes

X-ray polarization provides a new way to probe accretion geometry in black hole systems. If the accretion geometry of black holes is similar regardless of mass, we should expect the same to be true of their polarization properties. We compare the polarimetric properties of all nonblazar black holes observed with the Imaging X-ray Polarimetry Explorer. We find that their polarization properties are very similar, particularly in the hard state, where the corona dominates. This tentatively supports the idea that stellar and supermassive black holes share a common coronal geometry.

Inference of black-hole mass fraction in Galactic globular clusters

Context. Globular clusters (GCs) are suggested to host many stellar-mass black holes (BHs) at their centers, thus resulting in ideal testbeds for BH formation and retention theories. BHs are expected to play a major role in GC structural and dynamical evolution and their study has attracted a lot of attention. In recent years, several works attempted to constrain the BH mass fraction in GCs typically by comparing a single observable (for example, mass segregation proxies) with scaling relations obtained from numerical simulations. Aims. We aim to uncover the possible intrinsic degeneracies in determining the BH mass fraction from single dynamical parameters and identify the possible parameter combinations that are able to break these degeneracies. Methods. We used a set of 101 Monte Carlo simulations sampling a large grid of initial conditions. In particular, we explored the impact of different BH natal kick prescriptions on widely adopted scaling relations. We then compared the results of our simulations with observations obtained using state-of-the-art HST photometric and astrometric catalogs for a sample of 30 Galactic GCs. Results. We find that using a single observable to infer the present-day BH mass fraction in GCs is degenerate, as similar values could be attained by simulations including different BH mass fractions. We argue that the combination of mass-segregation indicators with GC velocity dispersion ratios could help us to break this degeneracy efficiently. We show that such a combination of parameters can be derived with currently available data. However, the limited sample of stars with accurate kinematic measures and its impact on the overall errors do not allow us to discern fully different scenarios yet.

Neutrino Oscillation Effects on the Luminosity of Neutrino-dominated Accretion Flows around Black Holes

A stellar-mass black hole surrounded by the neutrino-dominated accretion flow (NDAF) has been proposed as the central engine of gamma-ray bursts. In this work, we investigate the neutrino/antineutrino luminosity and annihilation luminosity from the NDAF and pair annihilation model, taking into account the neutrino oscillation above the accretion disk. The disk's hydrodynamical properties are modeled using an empirical solution previously derived with boundary conditions, including the effect of electron degeneracy and neutrino trapping. Our key parameters are the mass accretion rate and the black hole spin given in the ranges of Ṁ=0.1 –10 M ⊙ s−1 and 0 ≤ a < 1, respectively. Without neutrino oscillation, the obtained neutrino/antineutrino luminosity is ∼1051–1053 erg s−1, while the neutrino annihilation luminosity is found to be ∼1046–1051 erg s−1. In the presence of neutrino oscillation in the vacuum limit, the electron neutrino annihilation luminosity decreases by ≲22% through the flavor transformation, while the muon and tau neutrino luminosity can increase by up to ∼45% and 60%, respectively. As a result, the total annihilation luminosity can be reduced by up to ∼19% due to the oscillation process above the disk. Finally, we also investigate the case whereby the charge-conjugation-parity-symmetry-violating (CP-violating) phase is changed from δCP = 0° to 245°. However, our results reveal that the CP-violating phase has minimal impact on the neutrino annihilation luminosity.

The effect of stellar rotation on black hole mass and spin

Abstract The gravitational wave signature of a binary black hole (BBH) merger is dependent on its component mass and spin. If such black holes originate from rapidly rotating progenitors, the large angular momentum reserve in the star could drive a collapsar-like supernova explosion, hence substantially impacting these characteristics of the black holes in the binary. To examine the effect of stellar rotation on the resulting black hole mass and spin, we conduct a 1D general relativistic study of the end phase of the stellar collapse. We find that the resulting black hole mass at times differs significantly from the previously assumed values. We quantify the dependence of the black hole spin magnitude on the hydrodynamics of the accretion flow, providing analytical relations for calculating the mass and spin based on the progenitor’s pre-collapse properties. Depending on the nature of the accretion flow, our findings have implications for the black hole upper mass gap resulting from pair-instability supernovae, the maximum mass of a maximally rotating stellar black hole, and the maximum effective spin of a BBH formed in a tidally locked helium star - black hole binary.

Eccentric signatures of stellar-mass binary black holes with circumbinary disks in LISA

Abstract Stellar-mass binary black holes may have circumbinary disks if formed through common-envelope evolution or within gaseous environments. Disks can drive binaries into wider and more eccentric orbits, while gravitational waves harden and circularise them. We combine cutting-edge evolution prescriptions for disk-driven binaries with well-known equations for gravitational-wave-driven evolution, and study the evolution of stellar-mass binary black holes. We find that binaries are driven by their disk to an equilibrium eccentricity, 0.2 ≲ eeq ≲ 0.5, that dominates their evolution. Once they transition to the GW-dominated regime their eccentricity decreases rapidly; we find that stellar-mass binary black holes with long-lived disks will likely be observed in LISA with detectable eccentricities ∼10−2 at 0.01 Hz, with the precise value closely correlating with the binary’s initial mass ratio. This may lead stellar-mass binary black holes with CBDs observed in LISA to be confused with dynamically-formed binary black holes.