Search for: All records

Dynamical interactions in dense star clusters could significantly influence the properties of black holes, leaving imprints on their gravitational-wave signatures. While previous studies have mostly focused on repeated black hole mergers for spin and mass growth, this work examines the impact of physical collisions and close encounters between black holes and (noncompact) stars. Using Monte CarloN-body models of dense star clusters, we find that a large fraction of black holes retained upon formation undergo collisions with stars. Within our explored cluster models, the proportion of binary black hole mergers affected by stellar collisions ranges from 10%–60%. If all stellar-mass black holes are initially nonspinning, we find that up to 40% of merging binary black holes may have components with dimensionless spin parameterχ ≳ 0.2 because of prior stellar collisions, while typically about 10% have spins nearχ = 0.7 from prior black hole mergers. We demonstrate that young star clusters are especially important environments, as they can produce collisions of black holes with very massive stars, allowing for significant spin-up of the black holes through accretion. Our predictions for black hole spin distributions from these stellar collisions highlight their sensitivity to accretion efficiency, underscoring the need for detailed hydrodynamic calculations to better understand the accretion physics following these interactions.

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 (∼10⁴ M_⊙ pc⁻²), 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 10⁴¹ − 10⁴⁴ erg/s, depending on the GMC properties and specific feedback model assumed. We find that the maximum possible mass growth of seed BHs (ΔM^max_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 ∼10⁴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.

The existence of compact stellar remnants in the mass range 2–5M_⊙has long been debated. This so-called lower-mass gap (LMG) was initially suggested by the lack of low-mass X-ray binary observations with accretors about 2–5M_⊙, but it has recently been called into question following newer observations, including an LMG candidate with a millisecond pulsar (MSP) companion in the dense globular cluster NGC 1851. Here, we model NGC 1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte CarloN-body code Cluster Monte Carlo, and we specifically study the formation of LMG black holes (BHs). We demonstrate that both massive star evolution and dynamical interactions can contribute to forming LMG BHs. In general, the collapse of massive remnants formed through mergers of neutron stars (NSs) or massive white dwarfs produces the largest number of LMG BHs among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC 1851-like models can reproduce MSP—LMG BH binaries similar to the observed system. Additionally, the LMG BHs can also become components of dynamically assembled binaries, and some will be in merging BH–NS systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably ≲1 Gpc⁻³yr⁻¹.

The origins and mergers of supermassive black holes (SMBHs) remain a mystery. We describe a scenario from a novel multiphysics simulation featuring rapid (≲1 Myr) hyper-Eddington gas capture by a ∼1000M_⊙“seed” black hole (BH) up to supermassive (≳10⁶M_⊙) masses in a massive, dense molecular cloud complex typical of high-redshift starbursts. Due to the high cloud density, stellar feedback is inefficient, and most of the gas turns into stars in star clusters that rapidly merge hierarchically, creating deep potential wells. Relatively low-mass BH seeds at random positions can be “captured” by merging subclusters and migrate to the center in ∼1 freefall time (vastly faster than dynamical friction). This also efficiently produces a paired BH binary with ∼0.1 pc separation. The centrally concentrated stellar density profile (akin to a “protobulge”) allows the cluster as a whole to capture and retain gas and build up a large (parsec-scale) circumbinary accretion disk with gas coherently funneled to the central BH (even when the BH radius of influence is small). The disk is “hypermagnetized” and “flux-frozen”: dominated by a toroidal magnetic field with plasmaβ∼ 10⁻³, with the fields amplified by flux-freezing. This drives hyper-Eddington inflow rates ≳1M_⊙yr⁻¹, which also drive the two BHs to nearly equal masses. The late-stage system appears remarkably similar to recently observed high-redshift “little red dots.” This scenario can provide an explanation for rapid SMBH formation, growth, and mergers in high-redshift galaxies.

Gravitational-wave observations provide the unique opportunity of studying black hole formation channels and histories—but only if we can identify their origin. One such formation mechanism is the dynamical synthesis of black hole binaries in dense stellar systems. Given the expected isotropic distribution of component spins of binary black holes in gas-free dynamical environments, the presence of antialigned or in-plane spins with respect to the orbital angular momentum is considered a tell-tale sign of a merger’s dynamical origin. Even in the scenario where birth spins of black holes are low, hierarchical mergers attain large component spins due to the orbital angular momentum of the prior merger. However, measuring such spin configurations is difficult. Here, we quantify the efficacy of the spin parameters encoding aligned-spin (χ_eff) and in-plane spin (χ_p) at classifying such hierarchical systems. Using Monte Carlo cluster simulations to generate a realistic distribution of hierarchical merger parameters from globular clusters, we can infer mergers’χ_effandχ_p. The cluster populations are simulated using Advanced LIGO-Virgo sensitivity during the detector network’s third observing period and projections for design sensitivity. Using a “likelihood-ratio”-based statistic, we find that ∼2% of the recovered population by the current gravitational-wave detector network has a statistically significantχ_pmeasurement, whereas noχ_effmeasurement was capable of confidently determining a system to be antialigned with the orbital angular momentum at current detector sensitivities. These results indicate that measuring spin-precession throughχ_pis a more detectable signature of hierarchical mergers and dynamical formation than antialigned spins.

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 1000M_⊙. 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.

We apply for the first time orbit-averaged Monte Carlo star cluster simulations to study tidal tail and stellar stream formation from globular clusters (GCs), assuming a circular orbit in a time-independent spherical Galactic potential. Treating energetically unbound bodies—potential escapers (PEs)—as collisionless enables this fast but spherically symmetric method to capture asymmetric extratidal phenomena with exquisite detail. Reproducing stream features such as epicyclic overdensities, we show howreturning tidal tailscan form after the stream fully circumnavigates the Galaxy, enhancing the stream's velocity dispersion by several kilometers per second in our ideal case. While a truly clumpy, asymmetric, and evolving Galactic potential would greatly diffuse such tails, they warrant scrutiny as potentially excellent constraints on the Galaxy’s history and substructure. Reexamining the escape timescale Δtof PEs, we find new behavior related to chaotic scattering in the three-body problem; the Δtdistribution features sharp plateaus corresponding to distinct locally smooth patches of the chaotic saddle separating the phase-space basins of escape. We study for the first time Δtin an evolving cluster, finding that$\mathrm{\Delta }t\sim (E_{\mathrm{J}}^{-0.1},E_{\mathrm{J}}^{-0.4})$for PEs with (low, high) Jacobi energyE_J, flatter than for a static cluster ($E_{\mathrm{J}}^{-2}$). Accounting for cluster mass loss and internal evolution lowers the median Δtfrom ∼10 Gyr to ≲100 Myr. We finally outline potential improvements to escape in the Monte Carlo method intended to enable the first large grids of tidal tail/stellar stream models from full GC simulations and detailed comparison to stream observations.

Recently, we demonstrated self-consistent formation of strongly-magnetized quasar accretion disks (QADs) from cosmological radiation-magnetohydrodynamic-thermochemical galaxy-star formation simulations, including the full STARFORGE physics shown previously to produce a reasonable IMF under typical ISM conditions. Here we study star formation and the stellar IMF in QADs, on scales from 100 au to 10 pc from the SMBH. We show it is critical to include physics often previously neglected, including magnetic fields, radiation, and (proto)stellar feedback. Closer to the SMBH, star formation is suppressed, but the (rare) stars that do form exhibit top-heavy IMFs. Stars can form only in special locations (e.g. magnetic field switches) in the outer QAD. Protostars accrete their natal cores rapidly but then dynamically decouple from the gas and ‘wander,’ ceasing accretion on timescales ~100 yr. Their jets control initial core accretion, but the ejecta are ‘swept up’ into the larger-scale QAD flow without much dynamical effect. The strong tidal environment strongly suppresses common-core multiplicity. The IMF shape depends sensitively on un-resolved dynamics of protostellar disks (PSDs), as the global dynamical times can become incredibly short (< yr) and tidal fields are incredibly strong, so whether PSDs can efficiently transport angular momentum or fragment catastrophically at <10 au scales requires novel PSD simulations to properly address. Most analytic IMF models and analogies with planet formation in PSDs fail qualitatively to explain the simulation IMFs, though we discuss a couple of viable models.

Globular clusters (GCs) are particularly efficient at forming millisecond pulsars. Among these pulsars, about half lack a companion star, a significantly higher fraction than in the Galactic field. This fraction increases further in some of the densest GCs, especially those that have undergone core collapse, suggesting that dynamical interaction processes play a key role. For the first time, we createN-body models that reproduce the ratio of single-to-binary pulsars in Milky Way–like GCs. We focus especially on NGC 6752, a typical core-collapsed cluster with many observed millisecond pulsars. Previous studies suggested that an increased rate of neutron star binary disruption in the densest clusters could explain the overabundance of single pulsars in these systems. Here, we demonstrate that binary disruption is ineffective and instead we propose that two additional dynamical processes play dominant roles: (1) tidal disruption of main-sequence stars by neutron stars and (2) gravitational collapse of heavy white dwarf binary merger remnants. Neutron stars formed through these processes may also be associated with fast radio bursts similar to those observed recently in an extragalactic GC.