Search for: All records

The upcoming Laser Interferometer Space Antenna (LISA) is expected to detect gravitational waves (GWs) from massive black hole binaries (MBHB). Finding the electromagnetic (EM) counterparts for these GW events will be crucial for understanding how and where MBHBs merge, measuring their redshifts, constraining the Hubble constant and the graviton mass, and for other novel science applications. However, due to poor GW sky localization, multiwavelength, time-dependent EM models are needed to identify the right host galaxy. We studied merging MBHBs embedded in a circumbinary disc (CBD) using high-resolution two-dimensional simulations, with a Γ-law equation of state, incorporating viscous heating, shock heating, and radiative cooling. We simulate the binary from large separation until after merger, allowing us to model the decoupling of the binary from the CBD. We compute the EM signatures and identify distinct features before, during, and after the merger. Our main result is a multiband EM signature: we find that the MBHB produces strong thermal X-ray emission until 1–2 d prior to the merger. However, as the binary decouples from the CBD, the X-ray-bright minidiscs rapidly shrink in size, become disrupted, and the accretion rate drops precipitously. As a result, the thermal X-ray luminosity drops by orders of magnitude, and the source remains X-ray dark for several days, regardless of any post-merger effects such as GW recoil or mass-loss. Looking for the abrupt spectral change where the thermal X-ray disappears is a tell-tale EM signature of LISA mergers that does not require extensive pre-merger monitoring.

 

We explore the role of galactic feedback on the low-redshift Lyα(Lyα) forest (z≲ 2) statistics and its potential to alter the thermal state of the intergalactic medium. Using the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) suite, we explore variations of the AGN and stellar feedback models in the IllustrisTNG and Simba subgrid models. We find that both AGN and stellar feedback in Simba play a role in setting the Lyαforest column density distribution function (CDD) and the Doppler width (b-value) distribution. The Simba AGN jet feedback mode is able to efficiently transport energy out to the diffuse IGM, causing changes in the shape and normalization of the CDD and a broadening of theb-value distribution. We find that stellar feedback plays a prominent role in regulating supermassive black hole growth and feedback, highlighting the importance of constraining stellar and AGN feedback simultaneously. In IllustrisTNG, the AGN feedback variations explored in CAMELS do not affect the Lyαforest, but varying the stellar feedback model does produce subtle changes. Our results imply that the low-zLyαforest can be sensitive to changes in the ultraviolet background, stellar and black hole feedback, and that AGN jet feedback in particular can have a strong effect on the thermal state of the IGM.

 

The shrinking of a binary orbit driven by the interaction with a gaseous circumbinary disc, initially advocated as a potential way to catalyse the binary merger, has recently been debated in the case of geometrically thick (i.e. with H/R ≳ 0.1) discs. However, a clear consensus is still missing mainly owing to numerical limitations, such as fixed orbit binaries or lack of resolution inside the cavity carved by the binary in its circumbinary disc. In this work, we assess the importance of evolving the binary orbit by means of hydrodynamic simulations performed with the code gizmo in meshless finite mass mode. In order to model the interaction between equal mass circular binaries and their locally isothermal circumbinary discs, we enforce hyper-Lagrangian resolution inside the cavity. We find that fixing the binary orbit ultimately leads to an overestimate of the gravitational torque that the gas exerts on the binary and an underestimate of the torque due to the accretion of material on to the binary components. Furthermore, we find that the modulation of the accretion rate on the binary orbital period is strongly suppressed in the fixed orbit simulation, while it is clearly present in the live binary simulations. This has potential implications for the prediction of the observable periodicities in massive black hole binary candidates.

 

Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated haloes. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution function only require resolutions sufficient to keep the jet-inflated cavities stable. Comparing different model variations, as well as an independent jet model using a different hydrodynamics code, we show that the dominant uncertainties are the choices of jet properties within a given model. Independent of implementation, we find that light, thermal jets with low momentum flux tend to delay the onset of a cooling flow more efficiently on a 50 Myr time-scale than heavy, kinetic jets. The delay of the cooling flow originates from a displacement and boost in entropy of the central gas. If the jet kinetic luminosity depends on accretion rate, collimated, light, hydrodynamic jets are able to reduce cooling flows in haloes, without a need for jet precession or wide opening angles. Comparing the jet feedback with a ‘kinetic wind’ implementation shows that equal amounts of star formation rate reduction can be achieved by different interactions with the halo gas: the jet has a larger effect on the hot halo gas while leaving the denser, star-forming phase in place, while the wind acts more locally on the star-forming phase, which manifests itself in different time-variability properties.

 

The early growth of black holes (BHs) in high-redshift galaxies is likely feedback regulated. While radiative feedback has been extensively studied, the role of mechanical feedback has received less scrutiny to date. Here, we use high-resolution parsec-scale hydrodynamical simulations to study jet propagation and its effect on 100 M⊙ BH accretion in the dense, low-metallicity gas expected in early protogalaxies. As the jet propagates, it shocks the surrounding gas forming a jet cocoon. The cocoon consists of a rapidly cooling cold phase at the interface with the background gas and an overpressured subsonic phase of reverse shock-heated gas filling the interior. We vary the background gas density and temperature, BH feedback efficiency, and the jet model. We found that the width of the jet cocoon roughly follows a scaling derived by assuming momentum conservation in the jet-propagation direction and energy conservation in the lateral directions. Depending on the assumed gas and jet properties, the cocoon either stays elongated to large radii or isotropizes before reaching the Bondi radius, forming a nearly spherical bubble. Lower jet velocities and higher background gas densities result in self-regulation to higher momentum fluxes and elongated cocoons. In all cases, the outward cocoon momentum flux balances the inward inflowing gas momentum flux near the Bondi radius, which ultimately regulates BH accretion. The time-averaged accretion rate always remains below the Bondi rate, and exceeds the Eddington rate only if the ambient medium is dense and cold, and/or the jet is weak (low velocity and mass loading).

 

Fuzzy dark matter (FDM) is a proposed modification for the standard cold dark matter (CDM) model motivated by small-scale discrepancies in low-mass galaxies. Composed of ultralight (mass ∼ 10²²eV) axions with kiloparsec-scale de Broglie wavelengths, this is one of a class of candidates that predicts that the first collapsed objects form in relatively massive dark matter halos. This implies that the formation history of the first stars and galaxies would be very different, potentially placing strong constraints on such models. Here we numerically simulate the formation of the first stars in an FDM cosmology, following the collapse in a representative volume all the way down to primordial protostar formation including a primordial nonequilibrium chemical network and cooling for the first time. We find two novel results: first, the large-scale collapse results in a very thin and flat gas “pancake”; second, despite the very different cosmology, this pancake fragments until it forms protostellar objects indistinguishable from those in CDM. Combined, these results indicate that the first generation of stars in this model are also likely to be massive and, because of the sheet morphology, do not self-regulate, resulting in a massive Population III starburst. We estimate the total number of first stars forming in this extended structure to be 10⁴over 20 Myr using a simple model to account for the ionizing feedback from the stars, and should be observable with the James Webb Space Telescope. These predictions provide a potential smoking gun signature of FDM and similar dark matter candidates.

 

The detection of starlight from the host galaxies of quasars during the reionization epoch (z > 6) has been elusive, even with deep HST observations1,2. The current highest redshift quasar host detected3, at z = 4.5, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars4,5,6 from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP)7 mitigate the challenge of detecting their underlying, previously-undetected host galaxies. Here we report rest-frame optical images and spectroscopy of two HSC-SSP quasars at z > 6 with JWST. Using NIRCam imaging at 3.6μm and 1.5μm and subtracting the light from the unresolved quasars, we find that the host galaxies are massive (stellar masses of 13 × and 3.4 × 1010 M⊙, respectively), compact, and disk-like. NIRSpec medium-resolution spectroscopy shows stellar absorption lines in the more massive quasar, confirming the detection of the host. Velocity-broadened gas in the vicinity of these quasars enables measurements of their black hole masses (1.4 × 109 and 2.0 × 108 M⊙, respectively). Their location in the black hole mass - stellar mass plane is consistent with the distribution at low redshift, suggesting that the relation between black holes and their host galaxies was already in place less than a billion years after the Big Bang. 

Abstract Linear analysis of gas flows around orbiting binaries suggests that a centrifugal barrier ought to clear a low-density cavity around the binary and inhibit mass transfer onto it. Modern hydrodynamics simulations have confirmed the low-density cavity, but show that any mass flowing from large scales into the circumbinary disk is eventually transferred onto the binary components. Even though many numerical studies confirm this picture, it is still not understood precisely how gas parcels overcome the centrifugal barrier and ultimately accrete. We present a detailed analysis of the binary accretion process, using an accurate prescription for evolving grid-based hydrodynamics with Lagrangian tracer particles that track the trajectories of individual gas parcels. We find that binary accretion can be described in four phases: (1) gas is viscously transported through the circumbinary disk up to the centrifugal barrier at the cavity wall, (2) the cavity wall is tidally distorted into accretion streams consisting of near-ballistic gas parcels on eccentric orbits, (3) the portion of each stream moving inwards of an accretion horizon radius r ¯ ≃ a —the radius beyond which no material is returned to the cavity wall—becomes bound to a minidisk orbiting an individual binary component, and (4) the minidisk gas accretes onto the binary component through the combined effect of viscous and tidal stresses. 

Supermassive stars (SMSs) with masses of 𝑀∗ ≃ 104–105 M⊙ are invoked as possible seeds of high-redshift supermassive black holes, but it remains under debate whether their protostar indeed acquires sufficient mass via gas accretion overcoming radiative feedback. We investigate protostellar growth in dynamically heated atomic-cooling haloes (ACHs) found in recent cosmological simulations, performing three-dimensional radiation hydrodynamical (RHD) simulations that consider stellar evolution under variable mass accretion. We find that one of the ACHs feeds the central protostar at rates exceeding a critical value, above which the star evolves in a cool bloating phase and hardly produces ionizing photons. Consequently, the stellar mass reaches 𝑀∗ 􏰁 104 M⊙ unimpeded by radiative feedback. In the other ACH, where the mass supply rate is lower, the star spends most of its life as a hot main-sequence star, emitting intense ionizing radiation. Then, the stellar mass growth is terminated around 500 M⊙ by photoevaporation of the circumstellar disk. A series of our RHD simulations provide a formula of the final stellar mass determined either by stellar feedback or their lifetime as a function of the mass supply rate from the parent cloud in the absence of stellar radiation. Combining the results with the statistical properties of SMS-forming clouds in high-redshift quasar progenitor haloes, we construct a top-heavy mass distribution of primordial stars over 𝑀∗ ≃ 100–105 M⊙, approximately following a power-law spectrum of ∝ 𝑀−1.3 with a steeper decline at 𝑀 􏰁 2 × 104 M . Their massive BH remnants would be ∗∗⊙ further fed via the dense debris disk, powering “milli-quasars" with a bolometric luminosity of 𝐿bol 􏰁 1043 erg s−1. 

The existence of 109 M⊙ supermassive black holes (SMBHs) within the first billion years of the universe remains a puzzle in our conventional understanding of black hole formation and growth. The so-called direct-collapse scenario suggests that the formation of supermassive stars (SMSs) can yield the massive seeds of early SMBHs. This scenario leads to an overly massive BH galaxy (OMBG), whose nuclear black hole’s mass is comparable to or even greater than the surrounding stellar mass: a 104 − 106 M⊙ seed black hole is born in a dark matter halo with a mass as low as 107 − 108 M⊙. The black hole to stellar mass ratio is 𝑀bh/𝑀∗ ≫ 10−3, well in excess of the typical values at lower redshift. We investigate how long these newborn BHs remain outliers in the 𝑀bh − 𝑀∗ relation, by exploring the subsequent evolution of two OMBGs previously identified in the Renaissance simulations. We find that both OMBGs have𝑀bh/𝑀∗>1 during their entire life, from their birth at 𝑧≈15 until they merge with much more massive haloes at 𝑧 ≈ 8. We find that the OMBGs are spatially resolvable from their more massive, 1011 M⊙, neighboring haloes until their mergers are complete at 𝑧 ≈ 8. This affords a window for future observations with JWST and sensitive X-ray telescopes to diagnose the direct-collapse scenario, by detecting similar OMBGs and establishing their uniquely high black hole-to-stellar mass ratio.