Black Hole Mass Research Articles

Global analysis of LISA data with Galactic binaries and massive black hole binaries

Global analysis of LISA data with Galactic binaries and massive black hole binaries

Black Hole–Halo Mass Relation from UNIONS Weak Lensing

This Letter presents, for the first time, direct constraints on the black hole–halo mass relation using weak gravitational-lensing measurements. We construct type I and type II active galactic nucleus (AGN) samples from the Sloan Digital Sky Survey, with a mean redshift of 0.4 (0.1) for type I (type II) AGNs. This sample is cross correlated with weak-lensing shear from the Ultraviolet Near Infrared Optical Northern Survey. We compute the excess surface mass density of the halos associated with 36,181 AGNs from 94,308,561 lensed galaxies and fit the halo mass in bins of black hole mass. We find that more massive AGNs reside in more massive halos. The relation between halo mass and black hole mass is well described by a power law of slope 0.6 for both type I and type II samples, in agreement with models that link black hole growth to baryon feedback. We see no dependence on AGN type or redshift in the black hole–halo mass relation below a black hole mass of 108.5 M ⊙. Above that mass, we find more massive halos for the low-z type II sample compared to the high-z type I sample, but this difference may be interpreted as systematic error in the black hole mass measurements. Our results are consistent with previous measurements for non-AGN galaxies. At a fixed black hole mass, our weak-lensing halo masses are consistent with galaxy rotation curves but significantly lower than galaxy-clustering measurements. Finally, our results are broadly consistent with state-of-the-art hydrodynamical cosmological simulations, providing a new constraint for black hole masses in simulations.

Discovery of a 100 kpc Narrow Curved Twin Jet in the S-shaped Giant Radio Galaxy J0644+1043

We report the discovery of an S-shaped morphology of the radio galaxy J0644+1043 imaged with a 30 μJy beam−1 sensitive 525 MHz broadband (bands 3 + 4) uGMRT map. Dedicated spectroscopic observations of the host galaxy carried out with the 2 m Rozhen telescope yielded a redshift of 0.0488, giving a projected linear size of the peculiar radio structure of over 0.7 Mpc. This giant radio galaxy is powered by a black hole of mass 4.1−2.87+9.39×108 M ⊙, from the vicinity of which emanate well-collimated and knotty jets, each ∼100 kpc long. The entire radio structure, presumably due to the effective jet precession, is less than 50 Myr old, has a power of ∼6 × 1024 W Hz−1 at 1.4 GHz, and has observed morphological characteristics that do not strictly conform to the traditional Fanaroff–Riley (FR) FR I or FR II categories.

Forming off-center massive black hole binaries in dwarf galaxies through Jacobi capture

It is well established that black holes reside in the central regions of virtually all types of known galaxies. Recent observational and numerical studies however challenge this picture, suggesting that intermediate-mass black holes in dwarf galaxies may be found on orbits far from the center. In particular, constant-density cores minimize orbital energy losses due to dynamical friction, and allow black holes to settle on stable off-center orbits. Using controlled simulations, we study the dynamics of off-center black holes in dwarf galaxies with such cores. We propose a new scenario to describe off-center mergers of massive black holes, starting with a Jacobi capture. We focus on initially circular and co-planar black hole orbits and explore a large parameter space of black hole masses and orbital parameters. We find that Jacobi captures are a complex and chaotic phenomenon that occurs in about 13% of cases in this simplified setup, and we quantify how the likelihood of capture depends on the simulation parameters. We note that this percentage is likely an upper limit of the general case. Nevertheless, we show that Jacobi captures in cored dwarf galaxies can facilitate the formation of off-center black hole binaries, and that this process is sufficiently common to have a substantial effect on black hole growth. While our setup only allows for temporary captures, we expect dissipative forces from baryons and post-Newtonian corrections to maintain the captures over time and to lead to the formation of stable binary systems. This motivates future studies of the effectiveness of such dissipative forces, within stripped nuclei or globular clusters, in forming stable bound systems.

X-Ray Quasi-periodic Eruptions and Tidal Disruption Events Prefer Similar Host Galaxies

In the past 5 yr, six X-ray quasi-periodic eruption (QPE) sources have been discovered in the nuclei of nearby galaxies. Their origin remains an open question. We present Multi Unit Spectroscopic Explorer integral field spectroscopy of five QPE host galaxies to characterize their properties. We find that 3/5 galaxies host extended emission-line regions (EELRs) up to 10 kpc in size. The EELRs are photoionized by a nonstellar continuum, but the current nuclear luminosity is insufficient to power the observed emission lines. The EELRs are decoupled from the stars both kinematically and in projected sky position, and the low velocities and velocity dispersions (<100 km s−1 and ≲75 km s−1, respectively) are inconsistent with being driven by active galactic nuclei (AGNs) or shocks. The origin of the EELRs is likely a previous phase of nuclear activity. QPE host galaxies share several similarities with tidal disruption event (TDE) hosts, including an overrepresentation of galaxies with strong Balmer absorption and little ongoing star formation, as well as a preference for a short-lived (the typical EELR lifetime is ∼15,000 yr), gas-rich phase where the nucleus has recently faded significantly. This suggests that QPEs and TDEs may share a common formation channel, disfavoring AGN accretion disk instabilities as the origin of QPEs. If QPEs are related to extreme mass ratio inspiral systems (EMRIs), e.g., stellar-mass objects on bound orbits about massive black holes, the high incidence of EELRs and recently faded nuclei could be used to localize the hosts of EMRIs discovered by low-frequency gravitational-wave observatories.

A New Timescale–Mass Scaling for the Optical Variation of Active Galactic Nuclei across the Intermediate-mass to Supermassive Scales

The variability of active galactic nuclei (AGNs) has long been servicing as an essential avenue of exploring the accretion physics of a black hole (BH). There are two commonly used methods for analyzing AGN variability. First, the AGN variability, characterized by the structure function (SF) of a single band, can be well described by a damped random walk (DRW) process on timescales longer than ∼weeks, shorter than which departures have been reported. Second, the color variation (CV) between two bands behaves timescale dependent, raising challenges to the widely accepted reprocessing scenario. However, both the departure from the DRW process and the timescale-dependent CV are hitherto limited to AGNs, mainly quasars, at the supermassive scale. Here, utilizing the high-cadence multiwavelength monitoring on NGC 4395 harboring an intermediate-mass BH, we unveil at the intermediate-mass scale for the first time, prominent departures from the DRW process at timescales shorter than ∼hours in all three nights and bands, and plausible timescale-dependent CVs in the two longest nights of observation. Furthermore, comparing SFs of NGC 4395 to four AGNs at the supermassive scale, we suggest a new scaling relation between the timescale (τ; across nearly 3 orders of magnitude) and the BH mass (M BH): τ∝MBHγ where the exponent γ is likely ≃0.6–0.8. This exponent differs from most previous measurements, but confirms a few, and is consistent with a recent theoretical prediction, suggesting a similar accretion process in AGNs across different mass scales.

Constraining the Low-mass End of the Black Hole Mass Function and the Active Fraction of the Intermediate-mass Black Holes

We investigate the black hole mass function (BHMF) and the Eddington ratio distribution function (ERDF), focusing on the intermediate-mass black holes (IMBHs) with masses down to M • ∼ 104 M ⊙. Based on the active galactic nuclei (AGNs) with a detected broad Hα emission line, we construct a sample of 14,242 AGNs at redshift z < 0.35, including 243 IMBHs with M • < 106 M ⊙. By jointly modeling the BHMF and ERDF via the maximum posterior estimation, we find that the BHMF peaks at ∼106 M ⊙ and exhibits a relatively constant value of 10−4 Mpc−3 dex−1 at the low-mass end. By comparing the derived BHMF of type 1 AGNs with the galaxy mass function based on the updated black hole mass–host galaxy stellar mass relation, we derive the active fraction. We also determine the active fraction for all AGNs using the upper and lower limit of the type 1 fraction. The active fraction decreases from 15%–40% for massive galaxies (M ⋆ > 1010 M ⊙) to lower than ∼2% for dwarf galaxies with M ⋆ ∼ 108 M ⊙. These results suggest that the black hole occupation fraction is expected to be ∼50% for low-mass galaxies (M ⋆ ∼ 108.5–109 M ⊙) if the duty cycle is similar between IMBHs and supermassive black holes.

Ensemble Variability Properties of Active Galactic Nuclei in the SDSS DR17

We present the results from a study of ∼9600 broadline-selected active galactic nuclei (AGNs) with host galaxies detected from the Sloan Digital Sky Survey (SDSS) Data Release 17. We compute ensemble variability statistics based on a comparison of the original SDSS photometric data with spectrophotometric measurements obtained days to decades later in the Sloan g, r, and i bands. Galaxy and AGN templates have been fitted to the SDSS spectra to isolate the AGN component from the host galaxy. The sources have absolute magnitudes in the range −24 < M i < −18 and lie at redshifts less than z ∼ 0.9. A variability analysis reveals that the anticorrelation between luminosity and variability amplitude continues down to log(Lbol [erg s−1]) = 43.5, demonstrating that the relationship extends by 4 orders of magnitude in AGN luminosity. To further explore the connection between AGN luminosity and variability, we determine the black hole mass and the accretion rate through measurement of the Hβ line width and the monochromatic luminosity at rest frame 5100 Å. Our results suggest that the accretion rate is the dominant parameter impacting the amplitude of variability and that the anticorrelation between the accretion rate and amplitude extends to rates as low as 1% Eddington. Moreover, we also identify an anticorrelation between the variability amplitude and black hole mass, with the correlation appearing strongest among the AGNs with low accretion rates.

Big Bang Nucleosynthesis with f(R) Gravity Scalarons and Astrophysical Consequences

f(R) gravity is one of the serious alternatives of general relativity with a large range of astronomical consequences. In this work, we study Big Bang nucleosynthesis (BBN) in f(R) gravity theory. We consider a modification to gravity due to the existence of primordial black holes (PBHs) in the radiation era that introduce additional degrees of freedom known as scalarons. We calculate the light element abundances by using the BBN code PArthENoPE. It is found that for a range of scalaron mass (2.2 − 3.5) × 104 eV, the abundance of lithium is lowered by 3−4 times the value predicted by general relativistic BBN, which is a level desired to address the cosmological lithium problem. For the above scalaron mass range, the helium abundance is within the observed bound. However, the deuterium abundance is found to be increased by 3−6 times the observed primordial abundance. It calls for a high efficiency of stellar formation and evolution processes for the destruction of primordial deuterium, which is suggested as possible in scalaron gravity. A novel relation between scalaron mass and black hole mass has been used to show that the above scalaron mass range corresponds to PBHs of subplanetary mass (∼1019 g) serving as one of the potential candidates of nonbaryonic dark matter. We infer Big Bang equivalence of power-law f(R) gravity with PBHs that are detectable with upcoming gravitational wave detectors.

Comprehensive view of a z ∼ 6.5 radio-loud quasi-stellar object: From the radio to the optical/NIR to the X-ray band

We present a multi-wavelength analysis, from the radio to the X-ray band, of the redshift z = 6.44 VIK J2318−31 radio-loud quasi-stellar object, one of the most distant currently known of this class. The work is based on newly obtained observations (uGMRT, ATCA, and Chandra) as well as dedicated archival observations that have not yet been published (GNIRS and X-shooter). Based on the observed X-ray and radio emission, its relativistic jets are likely young and misaligned from our line of sight. Moreover, we can confirm, with simultaneous observations, the presence of a turnover in the radio spectrum at νpeak ∼ 650 MHz that is unlikely to be associated with self-synchrotron absorption. From the near-infrared spectrum we derived the mass of the central black hole, $ M_{\mathrm{BH}}=8.1^{+6.8}_{-5.6} \times 10^8\,{M_{\odot}} $, and the Eddington ratio, $ \lambda_{\mathrm{EDD}} = 0.8^{+0.8}_{-0.6} $, using broad emission lines as well as an accretion disc model fit to the continuum emission. Given the high accretion rate, the presence of a ∼8 × 108 M⊙ black hole at z = 6.44 can be explained by a seed black hole (∼104 M⊙) that formed at z ∼ 25, assuming a radiative efficiency ηd ∼ 0.1. However, by assuming ηd ∼ 0.3, as expected for jetted systems, the mass observed would challenge current theoretical models of black hole formation.

Reconstruction of Cosmic Black Hole Growth and Mass Distribution from Quasar Luminosity Functions at z > 4: Implications for Faint and Low-mass Populations in JWST

The evolution of the quasar luminosity function (QLF) is fundamental to understanding the cosmic evolution of black holes (BHs) through their accretion phases. In the era of the James Webb Space Telescope (JWST), Euclid, and Nancy Grace Roman Space Telescope, their unprecedented detection sensitivity and wide survey area can unveil the low-luminosity quasar and low-mass BH population, and provide new insights into quasar host galaxies. We present a theoretical model describing BH growth from initial seeding at z ≳ 20 to ∼ 4, incorporating the duration of accretion episodes, the distribution of Eddington ratios, and the mass dependency of BH accretion rates. By constraining the model parameters with the observed QLFs at 4 ≤ z ≤ 6 across a wide UV luminosity range, we find that the high-redshift BH population grows rapidly at z ≳ 6, and decelerates the pace in subsequent epochs. Toward lower redshifts (z < 6), mass-dependent accretion inhibits the growth of high-mass BHs with M • > 108 M ⊙, leading to mass saturation at M • ≳ 1010 M ⊙. We predict the BH mass function down to M • ∼ 106 M ⊙ for both unobscured and obscured quasar populations at 4 ≤ z ≤ 11, offering a benchmark for future observational tests. Our model accounts for the presence of both bright and faint quasars at z > 4, including those discovered by JWST. Furthermore, our findings suggest two distinct pathways for the early assembly of the BH–galaxy mass correlation: the population with a BH-to-stellar-mass ratio near the local value of M •/M ⋆ ≃ 5 × 10−3 maintains proximity to the relation via moderate growth, while the population that begins to grow above the local relation becomes as overmassive as M •/M ⋆ ∼ 0.01–0.1 by z ∼ 6 via rapid mass accretion.

X-Ray View of Little Red Dots: Do They Host Supermassive Black Holes?

The discovery of Little Red Dots (LRDs)—a population of compact, high-redshift, dust-reddened galaxies—is one of the most surprising results from JWST. However, the nature of LRDs is still debated: does the near-infrared emission originate from accreting supermassive black holes (SMBHs), or intense star formation? In this work, we utilize ultra-deep Chandra observations and study LRDs residing behind the lensing galaxy cluster, A2744. We probe the X-ray emission from individual galaxies but find that they remain undetected and provide SMBH mass upper limits of ≲(1.5–16) × 106 M ⊙ assuming Eddington limited accretion. To increase the signal-to-noise ratios, we conduct a stacking analysis of the full sample with a total lensed exposure time of ≈87 Ms. We also bin the galaxies based on their stellar mass, lensing magnification, and detected broad-line Hα emission. For the LRDs exhibiting broad-line Hα emission, there is a hint of a stacked signal (∼2.6σ), corresponding to an SMBH mass of ∼3.2 × 106 M ⊙. Assuming unobscured, Eddington-limited accretion, this black hole (BH) mass is at least 1.5 orders of magnitude lower than that inferred from virial mass estimates using JWST spectra. Given galaxy-dominated stellar mass estimates, our results imply that LRDs do not host overmassive SMBHs and/or accrete at a few percent of their Eddington limit. However, alternative stellar mass estimates may still support that LRDs host overmassive BHs. The significant discrepancy between the JWST and Chandra data hints that the scaling relations used to infer the SMBH mass from the Hα line and virial relations may not be applicable for high-redshift LRDs.

Viscous torque in turbulent magnetized active galactic nucleus accretion disks and its effects on the gravitational waves of extreme mass ratio inspirals

The merger of supermassive black holes produces millihertz gravitational waves (GWs), which are potentially detectable by the future Laser Interferometer Space Antenna (LISA). Such binary systems are usually embedded in an accretion disk environment at the center of an active galactic nucleus (AGN). Recent studies suggest the plasma environment imposes measurable imprints on the GW signal if the mass ratio of the binary is around q ∼ 10−4 − 10−3. The effect of the gaseous environment on the GW signal is strongly dependent on the disk’s parameters; therefore, it is believed that future low-frequency GW detections will provide us with precious information about the physics of AGN accretion disks. We investigated this effect by measuring the viscous torque via modeling of the evolution of magnetized tori around the primary massive black hole. Using the general relativistic magnetohydrodynamic HARM-COOL code, we performed 2D and 3D simulations of weakly magnetized, thin accretion disks, with a possible truncation and transition to advection-dominated accretion flow. We studied the angular momentum transport and turbulence generated by the magnetorotational instability. We quantified the disk’s effective alpha viscosity and its evolution over time. We applied our numerical results to quantify the relativistic viscous torque on a hypothetical low-mass secondary black hole via a 1D analytical approach, and we estimated the GW phase shift due to the gas environment.

Correspondence of Concept of Circular Orbits in Schwarzschild Metric and Kerr Metric to those in Tidally Evolving Binaries in Weak Gravitational Regime and its Application in Determination of Spin of Supermassive Black Hole

This paper rederives diameters of the Innermost Stable Circular Orbits (ISCO) in Schwarzschild Metric (RDISCO = 6RSch for non-rotating BH with a =0) and in Kerr metric with particular reference to recently studied Super Massive Black Hole in NGC1365 (RDISCO = 9RSch for retrograde orbits and RDISCO = 1RSch for prograde orbits around maximally spinning BH of a=1). Next the same derivations are made in tidally evolving binaries based on kinematic model. Kinematic Model derivation corresponds to the derivations made in Kerr Metric for Circular Orbits in all its details at weak gravitation regime. At strong gravitation regime in Kerr metric, the two orbits converge to Innermost Stable Circular Orbit due to the Space-Time Curvature and Frame-Dragging Effect in the vicinity of SMBH. Kinematic model fails to arrive at ISCO since it has not included the relativistic mechanics.

Tidal Disruption Event AT2020ocn: Early Time X-Ray Flares Caused by a Possible Disk Alignment Process

A tidal disruption event (TDE) may occur when a star is torn apart by the tidal force of a black hole (BH). Eventually, an accretion disk is thought to form out of stellar debris falling back toward the BH. If the star’s orbital angular momentum vector prior to disruption is not aligned with the BH spin angular momentum vector, the disk will be tilted with respect to the BH equatorial plane. The disk will eventually be drawn into the BH equatorial plane due to a combination of the Bardeen–Petterson effect and internal torques. Here, we analyze the X-ray and UV observations of the TDE AT2020ocn obtained by Swift, XMM-Newton, and Neutron star Interior Composition ExploreR. The X-ray light curve shows strong flares during the first ≈100 days, while, over the same period, the UV emission decays gradually. We find that the X-ray flares can be explained by a model that also explains the spectral evolution. This model includes a slim disk viewed under a variable inclination plus an inverse-Comptonization component processing the slim disk emission. A scenario where the ongoing Lense–Thirring precession during the disk alignment process is responsible for the observed inclination variations is consistent with the data. In later observations, we find that the X-ray spectrum of AT2020ocn becomes harder, while the mass accretion rate remains at super-Eddington levels, suggesting the formation of a corona in line with accretion onto other compact objects. We constrain the BH mass to be (7−3+13)×105 M ⊙ at the 1σ (68%) confidence level.

Spectral distortions from acoustic dissipation with non-Gaussian (or not) perturbations

A well-known route to form primordial black holes in the early universe relies on the existence of unusually large primordial curvature fluctuations, confined to a narrow range of wavelengths that would be too small to be constrained by Cosmic Microwave Background (CMB) anisotropies. This scenario would however boost the generation of μ-type spectral distortions in the CMB due to an enhanced dissipation of acoustic waves. Previous studies of μ-distortion bounds on the primordial spectrum were based on the assumptions of Gaussian primordial fluctuations. In this work, we push the calculation of μ-distortions to one higher order in photon anisotropies. We discuss how to derive bounds on primordial spectrum peaks obeying non-Gaussian statistics under the assumption of local (perturbative or not) non-Gaussianity. We find that, depending on the value of the peak scale, the bounds may either remain stable or get tighter by several orders of magnitude, but only when the departure from Gaussian statistics is very strong. Our results are translated in terms of bounds on primordial supermassive black hole mass in a companion paper.

Inferring Small Neutron Star Spins with Neutron Star–Black Hole Mergers

The precise measurement of neutron star (NS) spins can provide important insight into the formation and evolution of compact binaries containing NSs. While traditional methods of NS spin measurement rely on pulsar observations, gravitational-wave detections offer a complementary avenue. However, determining component spins with gravitational waves is hindered by the small dimensionless spins of the NSs and the degeneracy in the mass and spin parameters. This degeneracy can be addressed by the inclusion of higher-order modes in the waveform, which are important for systems with unequal masses. This study shows the suitability of NS–black hole mergers, which are naturally mass-asymmetric, for precise NS spin measurements. We explore the effects of the black hole masses and spins, higher-mode content, inclination angles, and detector sensitivity on the measurement of NS spin. We find that networks with next-generation observatories like the Cosmic Explorer and the Einstein Telescope can distinguish NS dimensionless spin of 0.04 (0.1) from zero at 1σ confidence for events within ∼350 (∼1000) Mpc. Networks with A+ and A♯ detectors achieve similar distinction within ∼30 (∼70) Mpc and ∼50 (∼110) Mpc, respectively.

Low mass naked singularities from dark core collapse

Near-solar mass black holes (BHs) could have been involved in the two recent gravitational wave events, GW190425 and GW190814. Since such a low mass BH cannot be formed via stellar evolution, a model has been proposed based on the core collapse of a neutron star initiated by a certain number of dark matter (DM) particles. In this process, the accumulated DM particles collapse to form a tiny BH inside the neutron star, and the entire neutron star is transmuted into a BH after a certain time due to the accretion of matter by the endoparasitic BH from its host. Here, we argue that, depending on the initial conditions,a dark core collapse could give rise to either a BH or a naked singularity.For example, if the accumulated cloud of DM particles in the core of a neutron star can be modeled as an anisotropic fluid and it fulfils the criterion for collapse, an endoparasitic naked singularity could form instead of an endoparasitic BH.Immediately after its formation, the naked singularity should begin accreting matter from the host neutron star,thus eventually transmuting the entire host into a near-solar mass, relatively slowly-spinning naked singularity.We also propose a general technique to constrain the DM particle-neutron scattering cross section using the lack of pulsars near the Galactic centre and assuming that these missing pulsars have already been transmuted into BHs and/or naked singularities.Thus, the missing pulsars also indicate the existence of many such singularities near the Galactic center.

Simulations predict intermediate-mass black hole formation in globular clusters.

Intermediate-mass black holes (IMBHs) are those between 100 and 105 solar masses (M⨀); their formation process is debated. One potential origin is the growth of less massive black holes, by merging with stars and compact objects within globular clusters (GCs). However, previous simulations have indicated that this process only produces IMBHs <500 M⨀, before gravitational wave recoil ejects them from the GC. We perform star-by-star simulations of GC formation, finding that high-density star formation in a GC's parent giant molecular cloud can produce sufficient mergers of massive stars to overcome that mass threshold. We conclude that GCs can form with IMBHs ≳103M⨀, which is sufficiently massive to be retained within the GC even with the expected gravitational wave recoil.

Intermediate-mass Black Hole Progenitors from Stellar Collisions in Dense Star Clusters

Very massive stars (VMSs) formed via a sequence of stellar collisions in dense star clusters have been proposed as the progenitors of massive black hole seeds. VMSs could indeed collapse to form intermediate-mass black holes, which would then grow by accretion to become the supermassive black holes observed at the centers of galaxies and powering high-redshift quasars. Previous studies have investigated how different cluster initial conditions affect the formation of a VMS, including mass segregation, stellar collisions, and binaries, among others. In this study, we investigate the growth of VMSs with a new grid of Cluster Monte Carlo star cluster simulations—the most expansive to date. The simulations span a wide range of initial conditions, varying the number of stars, cluster density, stellar initial mass function (IMF), and primordial binary fraction. We find a gradual shift in the mass of the most massive collision product across the parameter space; in particular, denser clusters born with top-heavy IMFs provide strong collisional regimes that form VMSs with masses easily exceeding 1000 M ⊙. Our results are used to derive a fitting formula that can predict the typical mass of a VMS formed as a function of the star cluster properties. Additionally, we study the stochasticity of this process and derive a statistical distribution for the mass of the VMS formed in one of our models, recomputing the model 50 times with different initial random seeds.