Beecroft Institute for Particle Astrophysics and Cosmology

Abstract:

Supermassive black holes (BH) accrete gas from their surroundings and coalesce with companions during galaxy mergers, and both processes change the BH mass and spin. By means of high-resolution hydrodynamical simulations of galaxies, either idealised or embedded within the cosmic web, we explore the effects of interstellar gas dynamics and external perturbations on BH spin evolution. All these physical quantities were evolved on-the-fly in a self-consistent manner. We use a 'maximal' model to describe the turbulence induced by stellar feedback to highlight its impact on the angular momentum of the gas accreted by the BH. Periods of intense star formation are followed by phases where stellar feedback drives large-scale outflows and hot bubbles. We find that BH accretion is synchronised with star formation, as only when gas is cold and dense do both processes take place. During such periods, gas motion is dominated by consistent rotation. On the other hand, when stellar feedback becomes substantial, turbulent motion randomises gas angular momentum. However BH accretion is strongly suppressed in that case, as cold and dense gas is lacking. In our cosmological simulation, at very early times (z>6), the galactic disc has not yet settled and no preferred direction exists for the angular momentum of the accreted gas, so the BH spin remains low. As the gas settles into a disc (6>z>3), the BH spin then rapidly reaches its maximal value. At lower redshifts (z<3), even when galaxy mergers flip the direction of the angular momentum of the accreted gas, causing it to counter-rotate, the BH spin magnitude only decreases modestly and temporarily. Should this be a typical evolution scenario for BH, it potentially has dramatic consequences regarding their origin and assembly, as accretion on maximally spinning BH embedded in thin Shakura-Sunyaev disc is significantly reduced.

Abstract:

Supermassive black holes (BH) accrete gas from their surroundings and coalesce with companions during galaxy mergers, and both processes change the BH mass and spin. By means of high-resolution hydrodynamical simulations of galaxies, either idealised or embedded within the cosmic web, we explore the effects of interstellar gas dynamics and external perturbations on BH spin evolution. All these physical quantities were evolved on-the-fly in a self-consistent manner. We use a 'maximal' model to describe the turbulence induced by stellar feedback to highlight its impact on the angular momentum of the gas accreted by the BH. Periods of intense star formation are followed by phases where stellar feedback drives large-scale outflows and hot bubbles. We find that BH accretion is synchronised with star formation, as only when gas is cold and dense do both processes take place. During such periods, gas motion is dominated by consistent rotation. On the other hand, when stellar feedback becomes substantial, turbulent motion randomises gas angular momentum. However BH accretion is strongly suppressed in that case, as cold and dense gas is lacking. In our cosmological simulation, at very early times (z>6), the galactic disc has not yet settled and no preferred direction exists for the angular momentum of the accreted gas, so the BH spin remains low. As the gas settles into a disc (6>z>3), the BH spin then rapidly reaches its maximal value. At lower redshifts (z<3), even when galaxy mergers flip the direction of the angular momentum of the accreted gas, causing it to counter-rotate, the BH spin magnitude only decreases modestly and temporarily. Should this be a typical evolution scenario for BH, it potentially has dramatic consequences regarding their origin and assembly, as accretion on maximally spinning BH embedded in thin Shakura-Sunyaev disc is significantly reduced.

Abstract:

Cosmic ray (CR) feedback has been proposed as a powerful mechanism for driving warm gas outflows in galaxies. We use cosmological magnetohydrodynamic simulations to investigate the impact of CR feedback on neutral hydrogen (H i) in a 1011 M⊙ dark matter halo at 2 < z < 4. To this end, we postprocess the simulations with ionizing radiative transfer and perform Monte Carlo Lyman-α (Lyα) transfer calculations. CR feedback reduces H i column densities around young stars, thereby allowing more Lyα photons to escape and consequently offering a better match to the Lyα luminosities of observed Lyα emitters. Although galaxies with CR-driven outflows have more extended H i in the circumgalactic medium, two Lyα line properties sensitive to optical depth and gas kinematic—the location of the red peak relative to the Lyα line center in velocity space (vred) and relative strength of the blue-to-red peaks (B/R)—cannot distinguish between the CR-driven and non-CR simulations. This is because Lyα photons propagate preferentially along low H i density channels created by the ionizing radiation, thereby limiting the scattering with volume-filling H i. In contrast, the observed low flux ratios between the valley and peak and the surface brightness profiles are better reproduced in the model with CR-driven outflows because the Lyα photons interact more before escaping, rather than being destroyed by dust as is the case in the non-CR simulation. We discuss the potential cause of the paucity of sightlines in simulations that exhibit prominent red peaks and large vred, which may require the presence of more volume-filling H i.

Abstract:

Standard cosmological weak lensing analyses using cosmic shear are inevitably sensitive to small-scale, non-linear clustering from low-redshift structures. The need to adequately model the clustering of matter on this non-linear regime, accounting for both gravitational and baryonic effects, adds significant uncertainty to weak lensing studies, particularly in the context of near-future Stage-IV datasets. In this paper, inspired by previous work on so-called “nulling” techniques, we present a general method that selects the linear combinations of a given tomographic cosmic shear dataset that are least sensitive to small-scale non-linearities, by essentially suppressing the contribution from low-redshift structures. We apply this method to the latest public cosmic shear data from the Dark Energy Survey, DES-Y3, that corresponds to 3 years of observation, and show: a) that a large fraction of the signal is dominated by the single mode that is most affected by non-linear scales, and b) that removing this mode leads to a ∼ 1σ upwards shift in the preferred value of S 8 ≡ σ 8√(ΩM/0.3), alleviating the tension with current CMB data. However, the removal of the most contaminated mode also results in a significant increase in the statistical uncertainties. Taking this into account, we find this shift to be compatible with a random fluctuation caused by removing this most-contaminated mode at the ∼ 1.4σ level. We also show that this technique may be used by future Stage-IV surveys to mitigate the sensitivity of the final constraints to baryonic effects, trading precision for robustness.