Search

Supermassive black holes are prevalent at the centers of massive galaxies, and their masses scale with galaxy properties, increasing evidence suggesting that these trends continue to low stellar masses. Seeds are needed for supermassive black holes, especially at the highest redshifts explored by the James Webb Space Telescope. We study the hierarchical merging of galaxies via cosmological merger trees and argue that the seeds of supermassive black holes formed in nuclear star clusters via stellar black hole mergers at early epochs. Observable tracers include intermediate-mass black holes, nuclear star clusters, and early gas accretion in host dwarf galaxies, along with a potentially detectable stochastic gravitational wave background, ejection of intermediate and supermassive black holes, and consequences of a significant population of tidal disruption events and extreme-mass ratio inspirals., Comment: 20 pages, 12 figures

Observations of supermassive black holes at high redshift challenge our understanding of the evolution of the first generation of black holes (BHs) in proto-galactic environments. One possibility is that they grow much more rapidly than current estimates of feedback and accretion efficiency permit. Following our previous analysis of super-Eddington accretion onto stellar-mass black holes in mini-haloes under no-feedback conditions, we now investigate whether this can be sustained when thermal feedback is included. We use four sets of cosmological simulations at sub-pc resolution with initial black hole masses varying from $1 \times 10^3 - 6 \times 10^4 M_\odot$, exploring a range of feedback efficiencies. We also vary the feedback injection radius to probe the threshold of numerical overcooling. We find that super-Eddington growth sustained on the order of $\sim$$100 \, \rm kyr$ is possible with very weak thermal feedback efficiency in all environments and moderate efficiency for the $6 \times 10^4 M_\odot$ BH. Trans-Eddington growth is possible for a $3 \times 10^3 - 6 \times 10^6 M_\odot$ BH at moderate feedback efficiencies. We discuss the effectiveness of thermal feedback in heating the gas, suppressing accretion, and driving outflows at these parameter configurations. Our results suggest that super-Eddington growth may be possible in the presence of thermal feedback for black holes formed from the first stars., Comment: Submitted to MNRAS. 17 pages, 9 figures

After a merger of two massive black holes (MBHs), the remnant receives a gravitational wave (GW) recoil kick that can have a strong effect on its future evolution. The magnitude of the kick ($v_\mathrm{recoil}$) depends on the mass ratio and the alignment of the spins and orbital angular momenta, therefore on the previous evolution of the MBHs. We investigate the cosmic effect of GW recoil by running for the first time a high-resolution cosmological simulation including GW recoil that depends on the MBH spins (evolved through accretion and mergers), masses and dynamics computed self-consistently. We also run a twin simulation without GW recoil. The simulations are run down to $z=4.4$. We find that GW recoil reduces the growth of merger remnants, and can have a significant effect on the MBH-galaxy correlations and the merger rate. We find large recoil kicks across all galaxy masses in the simulation, up to a few $10^{11}\,\rm M_\odot$. The effect of recoil can be significant even if the MBHs are embedded in a rotationally supported gaseous structure. We investigate the dynamics of recoiling MBHs and find that MBHs remain in the centre of the host galaxy for low $v_\mathrm{recoil}/v_\mathrm{esc}$ and escape rapidly for high $v_\mathrm{recoil}/v_\mathrm{esc}$. Only if $v_\mathrm{recoil}$ is comparable to $v_\mathrm{esc}$ the MBHs escape the central region of the galaxy but might remain as wandering MBHs until the end of the simulation. Recoiling MBHs are a significant fraction of the wandering MBH population. Although the dynamics of recoiling MBHs may be complex, some retain their initial radial orbits but are difficult to discern from other wandering MBHs on radial orbits. Others scatter with the halo substructure or circularise in the asymmetric potential. Our work highlights the importance of including GW recoil in cosmological simulation models., Comment: submitted to A&A

The detection of gravitational waves (GWs) from massive black hole binary (MBHB) coalescence motivates the development of a sub-grid model. We present RAMCOAL, integrated into the RAMSES code, which simulates the orbital evolution of MBHBs, accounting for stellar and gaseous dynamical friction (DF), stellar scattering, circumbinary disk interactions, and GW emission at scales below the simulation resolution. Unlike post-processing approaches, RAMCOAL tracks the real-time evolution of MBHBs within hydrodynamical simulations of galaxies using local quantities to model dynamics and accretion. This enables more accurate predictions of both GW signals and the properties of merging black holes. We validate RAMCOAL across isolated and merging galaxy setups at resolutions of 10, 50, and 100 pc, with and without black hole accretion and feedback. In addition, we test the model in seven galaxy merger scenarios at 100 pc resolution. These tests demonstrate that RAMCOAL is largely resolution-independent and successfully captures the effects of DF from stars, dark matter, and gas, loss-cone scattering, viscous drag from circumbinary disks, and GW emission -- all within a realistic galactic environment, even at low resolutions. With RAMCOAL, we can better estimate MBHB coalescence rates and the GW background, while providing insights into the electromagnetic counterparts of GW sources. This approach bridges the gap between electromagnetic observations and GW detection, offering a more comprehensive understanding of MBHB evolution in cosmological simulations., Comment: Submitted to Astronomy & Astrophysics, 20 pages, 11 figures

Astrophysical black holes (BHs) have two fundamental properties: mass and spin. While the mass-evolution of BHs has been extensively studied, much less work has been done on predicting the distribution of BH spins. In this paper we present the spin evolution for a sample of intermediate-mass and massive BHs from the newHorizon simulation, which evolved BH spin across cosmic time in a full cosmological context through gas accretion, BH-BH mergers and BH feedback including jet spindown. As BHs grow, their spin evolution alternates between being dominated by gas accretion and BH mergers. Massive BHs are generally highly spinning. Accounting for the spin energy extracted through the Blandford-Znajek mechanism increases the scatter in BH spins, especially in the mass range $10^{5-7} \rm \ M_\odot$, where BHs had previously been predicted to be almost universally maximally spinning. We find no evidence for spin-down through efficient chaotic accretion. As a result of their high spin values, massive BHs have an average radiative efficiency of $<\varepsilon_{\rm r}^{\rm thin}> \approx 0.19$. As BHs spend much of their time at low redshift with a radiatively inefficient thick disc, BHs in our sample remain hard to observe. Different observational methods probe different sub-populations of BHs, significantly influencing the observed distribution of spins. Generally, X-ray-based methods and higher luminosity cuts increase the average observed BH spin. When taking BH spin evolution into account, BHs inject on average between 3 times (in quasar mode) and 8 times (in radio mode) as much feedback energy into their host galaxy as previously assumed., Comment: 21 pages, accepted for publication in MNRAS

We introduce a model for dust evolution in the RAMSES code for simulations of galaxies with a resolved multiphase interstellar medium. Dust is modelled as a fluid transported with the gas component, and is decomposed into two sizes, 5 nm and 0.1 $\mu\rm m$, and two chemical compositions for carbonaceous and silicate grains. Using a suite of isolated disc simulations with different masses and metallicities, the simulations can explore the role of these processes in shaping the key properties of dust in galaxies. The simulated Milky Way analogue reproduces the dust-to-metal mass ratio (DTM), depletion factors, size distribution and extinction curves of the Milky Way. Galaxies with lower metallicities reproduce the observed decrease in the DTM with metallicity at around a few 0.1 $\rm Z_\odot$. This break in the DTM corresponds to a galactic gas metallicity threshold that marks the transition from an ejecta-dominated to an accretion-dominated grain growth, and that is different for silicate and carbonaceous grains, with $\simeq$ 0.1 $\rm Z_\odot$ and $\simeq$ 0.5 $\rm Z_\odot$ respectively. This leads to more Magellanic Cloud-like extinction curves, i.e. with steeper slopes in the ultraviolet and a weaker bump feature at 217.5 nm, in galaxies with lower masses and lower metallicities. Steeper slopes in these galaxies are caused by the combination of the higher efficiency of gas accretion by silicate relative to carbonaceous grains and by the low rates of coagulation that preserves the amount of small silicate grains. Weak bumps are due to the overall inefficient accretion growth of carbonaceous dust at low metallicity, whose growth is mostly supported by the release of large grains in SN ejecta. We also show that the formation of CO molecules is a key component to limit the ability of carbonaceous dust to grow, in particular in low-metallicity gas-rich galaxies., Comment: A&A accepted

(Reduced) Using the recent cosmological high-resolution zoom-in simulations, NewHorizon and Galactica, in which the evolution of black hole spin is followed on the fly, we have tracked the cosmic history of a hundred of black holes (BHs) with a mass greater than 2x10^4 Ms. For each of them, we have studied the variations of the three dimensional angle (Psi) subtended between the BH spins and the angular momentum vectors of their host galaxies. The analysis of the individual evolution of the most massive BHs suggests that they are generally passing by three different regimes. First, for a short period after their birth, low mass BHs (<3x10^4 Ms) are rapidly spun up by gas accretion and their spin tends to be aligned with their host galaxy spin. Then follows a second phase in which the accretion of gas onto low mass BHs (<10^5 Ms) is quite chaotic and inefficient, reflecting the complex and disturbed morphologies of forming proto-galaxies at high redshifts. The variations of Psi are rather erratic during this phase and are mainly driven by the rapid changes of the direction of the galaxy angular momentum. Then, in a third and long phase, BHs are generally well settled in the center of galaxies around which the gas accretion becomes much more coherent (>10^5 Ms). In this case, the BH spins tend to be well aligned with the angular momentum of their host galaxy and this configuration is generally stable even though BH merger episodes can temporally induce misalignment. We have also derived the distributions of cos(Psi) at different redshifts and found that BHs and galaxy spins are generally aligned. Finally, based on a Monte Carlo method, we also predict statistics for the 2-d projected spin-orbit angles lambda. In particular, the distribution of lambda traces well the alignment tendency in the 3-d analysis. Such predictions provide an interesting background for future observational analyses., Comment: Accepted for publication in A&A (12/03/2024)

In galaxy clusters, the hot intracluster medium (ICM) can develop a striking multi-phase structure around the brightest cluster galaxy. Much work has been done on understanding the origin of this central nebula, but less work has studied its eventual fate after the originally filamentary structure is broken into individual cold clumps. In this paper we perform a suite of 30 (magneto-)hydrodynamical simulations of kpc-scale cold clouds with typical parameters as found by galaxy cluster simulations, to understand whether clouds are mixed back into the hot ICM or can persist. We investigate the effects of radiative cooling, small-scale heating, magnetic fields, and (anisotropic) thermal conduction on the long-term evolution of clouds. We find that filament fragments cool on timescales shorter than the crushing timescale, fall out of pressure equilibrium with the hot medium, and shatter, forming smaller clumplets. These act as nucleation sites for further condensation, and mixing via Kelvin-Helmholtz instability, causing cold gas mass to double within 75 Myr. Cloud growth depends on density, as well as on local heating processes, which determine whether clouds undergo ablation- or shattering-driven evolution. Magnetic fields slow down but don't prevent cloud growth, with the evolution of both cold and warm phase sensitive to the field topology. Counter-intuitively, anisotropic thermal conduction increases the cold gas growth rate compared to non-conductive clouds, leading to larger amounts of warm phase as well. We conclude that dense clumps on scales of $500$ pc or more cannot be ignored when studying the long-term cooling flow evolution of galaxy clusters., Comment: 21 pages, 26 figures, Accepted for Publication in MNRAS

The orbital decay of massive black holes in galaxies in the aftermath of mergers is at the heart of whether massive black holes successfully pair and merge, leading to emission of low-frequency gravitational waves. The role of dynamical friction sourced from the gas distribution has been uncertain because many analytical and numerical studies have either focussed on a homogeneous medium or have not reached resolutions below the scales relevant to the problem, namely the Bondi-Hoyle-Lyttleton radius. We performed numerical simulations of a massive black hole moving in a turbulent medium in order to study dynamical friction from turbulent gas. We find that the black hole slows down to the sound speed, rather than the turbulent speed, and that the orbital decay is well captured if the Bondi-Hoyle-Lyttleton radius is resolved with at least five resolution elements. We find that the larger the turbulent eddies, the larger the scatter in dynamical friction magnitude, because of the stochastic nature of the problem, and also because of the larger over- and under-densities encountered by the black hole along its trajectory. Compared to the classic solution in a homogeneous medium, the magnitude of the force depends more weakly on the Mach number, and dynamical friction is overall more efficient for high Mach numbers, but less efficient towards and at the transonic regime., Comment: 14 pages, 16 figures

By performing three-dimensional hydrodynamical simulations of a galaxy in an isolated dark matter halo, we follow the evolution of the spin parameter $a$ of a black hole (BH) undergoing super-Eddington phases throughout its growth. This regime, suspected to be accompanied by powerful jet outflows, is expected to decrease the BH spin magnitude. We combine super-Eddington accretion with sub-Eddington phases (quasar and radio modes) and follow the BH spin evolution. Due to the low frequency of super-Eddington episodes, relativistic jets in this regime are not able to decrease the magnitude of the spin effectively, as thin disc accretion in the quasar mode inevitably increases the BH spin. The combination of super- and sub-Eddington accretion does not lead to a simple explicit expression for the spin evolution because of feedback from super-Eddington events. An analytical expression can be used to calculate the evolution for $a\lesssim0.3$, assuming the super-Eddington feedback is consistently weak. Finally, BHs starting with low spin magnitude are able to grow to the highest mass, and if they initially start misaligned with the galactic disc, they get a small boost of accretion through retrograde accretion., Comment: Updated the paper to the published A&A version

Centres of galaxy clusters must be efficiently reheated to avoid a cooling catastrophe. One potential reheating mechanism is anisotropic thermal conduction, which could transport thermal energy from intermediate radii to the cluster center. However, if fields are not re-randomised, anisotropic thermal conduction drives the heat buoyancy instability (HBI) which reorients magnetic field lines and shuts off radial heat fluxes. We revisit the efficiency of thermal conduction under the influence of spin-driven AGN jets in idealised magneto-hydrodynamical simulations with anisotropic thermal conduction. Despite the black hole spin's ability to regularly re-orientate the jet so that the jet-induced turbulence is driven in a quasi-isotropic fashion, the HBI remains efficient outside the central 50 kpc of the cluster, where the reservoir of heat is the largest. As a result, conduction plays no significant role in regulating the cooling of the intra-cluster medium if central active galactic nuclei are the sole source of turbulence. Whistler-wave driven saturation of thermal conduction reduces the magnitude of the HBI but does not prevent it., Comment: Resubmitted to A&A. Replaced the whistler wave based simulation after correcting a code bug in the original simulation, and updated results accordingly

One of the key physical processes that helps prevent strong cooling flows in galaxy clusters is the continued energy input from the central active galactic nucleus (AGN) of the cluster. However, it remains unclear how this energy is thermalised so that it can effectively prevent global thermal instability. One possible option is that a fraction of the AGN energy is converted into cosmic rays (CRs), which provide non-thermal pressure support, and can retain energy even as thermal energy is radiated away. By means of magneto-hydrodynamical simulations, we investigate how CR injected by the AGN jet influence cooling flows of a massive galaxy cluster. We conclude that converting a fraction of the AGN luminosity as low as 10\% into CR energy prevents cooling flows on timescales of billion years, without significant changes in the structure of the multi-phase intra-cluster medium. CR-dominated jets, by contrast, lead to the formation of an extended, warm central nebula that is supported by CR pressure. We report that the presence of CRs is not able to suppress the onset of thermal instability in massive galaxy clusters, but CR-dominated jets do significantly change the continued evolution of gas as it continues to cool from isobaric to isochoric. The CR redistribution in the cluster is dominated by advection rather than diffusion or streaming, but the heating by CR streaming helps maintain gas in the hot and warm phase. Observationally, self-regulating, CR-dominated jets produce a \gammaray~ flux in excess of current observational limits, but low CR fractions in the jet are not ruled out., Comment: 24 pages, accepted in A&A. Corrected typo in coauthor name

Super-Eddington accretion is one scenario that may explain the rapid assembly of $\sim 10^9\rm\, M_\odot$ supermassive black holes (BHs) within the first billion year of the Universe. This critical regime is associated with radiatively inefficient accretion and accompanied by powerful outflows in the form of winds and jets. By means of hydrodynamical simulations of BH evolution in an isolated galaxy and its host halo with 12 pc resolution, we investigate how super-Eddington feedback affects the mass growth of the BH. It is shown that super-Eddington feedback efficiently prevents BH growth within a few Myr. The super-Eddington accretion events remain relatively mild with typical rates of about 2-3 times the Eddington limit, because of the efficient regulation by jets in that regime. We find that these jets are powerful enough to eject gas from the centre of the host galaxy all the way up to galactic scales at a few kpc, but do not significantly impact gas inflows at those large scales. By varying the jet feedback efficiency, we find that weaker super-Eddington jets allow for more significant BH growth through more frequent episodes of super-Eddington accretion. We conclude that effective super-Eddington growth is possible, as we find that simulations with weak jet feedback efficiencies provide a slightly larger BH mass evolution over long periods of time ($\sim 80\,\rm Myr$) than that for a BH accreting at the Eddington limit., Comment: Updated the paper to the published A&A version

The occurrence of dual active galactic nuclei (AGN) on scales of a few tens of kpc can be used to study merger-induced accretion on massive black holes (MBHs) and to derive clues on MBH mergers, using dual AGN as a parent population of precursors. We investigate the properties of dual AGN in the cosmological simulation Horizon-AGN. We create catalogs of dual AGN selected with distance and luminosity criteria, plus sub-catalogs where further mass cuts are applied. We divide the sample into dual AGN hosted in different galaxies, on the way to a merger, and into those hosted in one galaxy, after the galaxy merger has happened. We find that the relation between MBH and galaxy mass is similar to that of general AGN population and we compare the properties of dual AGN also with a control sample, discussing differences and similarities in masses and Eddington ratios. The typical mass ratio of galaxy mergers associated to dual AGN is 0.2, with mass loss in the smaller galaxy decreasing the mass ratio as the merger progresses. Between 30 and 80 per cent of dual AGN with separations between 4 and 30~kpc can be matched to an ensuing MBH merger. The dual AGN fraction increases with redshift and with separation threshold, although above 50~kpc the increase of multiple AGN limits that of duals. Multiple AGN are generally associated with massive halos, and mass loss of satellites shapes the galaxy-halo relation., Comment: Published in MNRAS

Hydrodynamical cosmological simulations are increasing their level of realism by considering more physical processes and having greater resolution or larger statistics. However, usually either the statistical power of such simulations or the resolution reached within galaxies are sacrificed. Here, we introduce the NewHorizon project in which we simulate at high resolution a zoom-in region of $\sim(16\,\rm Mpc)^3$ that is larger than a standard zoom-in region around a single halo and is embedded in a larger box. A resolution of up to 34 pc is reached within galaxies; this allows the simulation to capture the multi-phase nature of the interstellar medium and the clumpy nature of the star formation process in galaxies. In this introductory paper, we present several key fundamental properties of galaxies and their black holes, including the galaxy mass function, cosmic star formation rate, galactic metallicities, the Kennicutt-Schmidt relation, the stellar-to-halo mass relation, galaxy sizes, stellar kinematics and morphology, gas content within galaxies and its kinematics, and the black hole mass and spin properties over time. The various scaling relations are broadly reproduced by NewHorizon with some differences with the standard observables. Owing to its exquisite spatial resolution, NewHorizon captures the inefficient process of star formation in galaxies, which evolve over time from being more turbulent, gas rich, and star bursting at high redshift. These high-redshift galaxies are also more compact, and they are more elliptical and clumpier until the level of internal gas turbulence decays enough to allow for the formation of discs. The NewHorizon simulation gives access to a broad range of galaxy formation and evolution physics at low-to-intermediate stellar masses, which is a regime that will become accessible in the near future through surveys such as the LSST., Comment: 29 pages, 28 figures, A&A accepted

We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from the big bang, down to $z = 3.5$. The simulation resolves haloes down to the atomic cooling limit, and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The Obelisk project is designed specifically to study the coevolution of high redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to high redshift protoclusters. We find that hydrogen reionization happens inside-out and is completed by $z \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$., Comment: 24+2 pages, 18+1 figures, published in Astronomy & Astrophysics

The Lyman-$\alpha$ forest is a powerful probe for cosmology, but it is also strongly impacted by galaxy evolution and baryonic processes such as Active Galactic Nuclei (AGN) feedback, which can redistribute mass and energy on large scales. We constrain the signatures of AGN feedback on the 1D power spectrum of the Lyman-$\alpha$ forest using a series of eight hydro-cosmological simulations performed with the Adaptative Mesh Refinement code RAMSES. This series starts from the Horizon-AGN simulation and varies the sub-grid parameters for AGN feeding, feedback and stochasticity. These simulations cover the whole plausible range of feedback and feeding parameters according to the resulting galaxy properties. AGNs globally suppress the Lyman-$\alpha$ power at all scales. On large scales, the energy injection and ionization dominate over the supply of gas mass from AGN-driven galactic winds, thus suppressing power. On small scales, faster cooling of denser gas mitigates the suppression. This effect increases with decreasing redshift. We provide lower and upper limits of this signature at nine redshifts between $z=4.25$ and $z=2.0$, making it possible to account for it at post-processing stage in future work given that running simulations without AGN feedback can save considerable amounts of computing resources. Ignoring AGN feedback in cosmological inference analyses leads to strong biases with 2\% shift on $\sigma_8$ and 1\% shift on $n_s$, which represents twice the standards deviation of the current constraints on $n_s$., Comment: Accepted for publication in MNRAS

The intrinsic correlations of galaxy shapes and orientations across the large-scale structure of the Universe are a known contaminant to weak gravitational lensing. They are known to be dependent on galaxy properties, such as their mass and morphologies. The complex interplay between alignments and the physical processes that drive galaxy evolution remains vastly unexplored. We assess the sensitivity of intrinsic alignments (shapes and angular momenta) to Active Galactic Nuclei -AGN- feedback by comparing galaxy alignment in twin runs of the cosmological hydrodynamical Horizon simulation, which do and do not include AGN feedback respectively. We measure intrinsic alignments in three dimensions and in projection at z=0 and z=1. We find that the projected alignment signal of all galaxies with resolved shapes with respect to the density field in the simulation is robust to AGN feedback, thus giving similar predictions for contamination to weak lensing. The relative alignment of galaxy shapes around galaxy positions is however significantly impacted, especially when considering high-mass ellipsoids. Using a sample of galaxy "twins" across simulations, we determine that AGN changes both the galaxy selection and their actual alignments. Finally, we measure the alignments of angular momenta of galaxies with their nearest filament. Overall, these are more significant in the presence of AGN as a result of the higher abundance of massive pressure-supported galaxies., Comment: 17 pages, 13 figures, accepted in MNRAS

We describe a physical model of the outflows produced as a result of gas accretion onto a black hole, and the resultant changes to star formation rates and efficiencies in galaxies, using the Radio-SAGE semi-analytic galaxy formation model. We show that the ratio of outflow rate to SFR of galaxies is mainly driven by black hole mass and virial halo mass, and show that the SFR is higher than the outflow rate at low black hole masses. The model consistently reproduces the observed evolution of star formation rate density from z = 6 to z = 0, as well as the trend of the stellar mass - halo mass relations. We show the characteristic growth of massive galaxies influenced by AGN feedback at different redshifts. We find feedback to be prevalent in the most massive galaxy halos, inhibiting the cooling catastrophe., Comment: 15 pages, 8 figures. Accepted for publication in MNRAS

Super-Eddington accretion onto massive black hole seeds may be commonplace in the early Universe, where the conditions exist for rapid accretion. Direct collapse black holes are often invoked as a possible solution to the observation of super massive black holes (SMBHs) in the pre-reionisation Universe. We investigate here how feedback, mainly in the form of bipolar jets, from super-Eddington accreting seed black holes will affect their subsequent growth. We find that, nearly independent of the mass loading of the bipolar jets, the violent outflows generated by the jets evacuate a region of approximately 0.1 pc surrounding the black hole seed. However, the jet outflows are unable to break free of the halo and their impact is limited to the immediate vicinity of the black hole. The outflows suppress any accretion for approximately a dynamical time. The gas then cools, recombines and falls back to the centre where high accretion rates are again observed. The overall effect is to create an effective accretion rate with values of between 0.1 and 0.5 times the Eddington rate. If this episodic accretion rate is maintained for order 500 million years then the black hole will increase in mass by a factor of between 3 and 300 but far short of the factor of $10^4$ required for the seeds to become the SMBHs observed at $z>6$. Therefore, direct collapse black holes born into atomic cooling haloes and which experience strong negative mechanical feedback will require external influences (e.g. rapid major mergers with other haloes) to promote efficient accretion and reach SMBH masses within a few hundred million years., Comment: 16 pages. MNRAS accepted

Born in rapidly evolving mini-halos during the first billion years of the Universe, super- massive black holes (SMBH) feed from gas flows spanning many orders of magnitude, from the cosmic web in which they are embedded to their event horizon. As such, accretion onto SMBHs constitutes a formidable challenge to tackle numerically, and currently requires the use of sub-grid models to handle the flow on small, unresolved scales. In this paper, we study the impact of resolution on the accretion pattern of SMBHs initially inserted at the heart of dense galactic gas clouds, using a custom super-Lagrangian refinement scheme to resolve the black hole (BH) gravitational zone of influence. We find that once the self-gravitating gas cloud host is sufficiently well re- solved, accretion onto the BH is driven by the cloud internal structure, independently of the BH seed mass, provided dynamical friction is present during the early stages of cloud collapse. For a pristine gas mix of hydrogen and helium, a slim disc develops around the BH on sub-parsec scales, turning the otherwise chaotic BH accretion duty cycle into an episodic one, with potentially important consequences for BH feedback. In the presence of such a nuclear disc, BH mass growth predominantly occurs when infalling dense clumps trigger disc instabilities, fuelling intense albeit short-lived gas accretion episodes., Comment: Resubmitted to mnras after reviewer comments, 24 pages

Whilst in galaxy-size simulations, supermassive black holes (SMBH) are entirely handled by sub-grid algorithms, computational power now allows the accretion radius of such objects to be resolved in smaller scale simulations. In this paper, we investigate the impact of resolution on two commonly used SMBH sub-grid algorithms; the Bondi-Hoyle-Lyttleton (BHL) formula for accretion onto a point mass, and the related estimate of the drag force exerted onto a point mass by a gaseous medium. We find that when the accretion region around the black hole scales with resolution, and the BHL formula is evaluated using local mass-averaged quantities, the accretion algorithm smoothly transitions from the analytic BHL formula (at low resolution) to a supply limited accretion (SLA) scheme (at high resolution). However, when a similar procedure is employed to estimate the drag force it can lead to significant errors in its magnitude, and/or apply this force in the wrong direction in highly resolved simulations. At high Mach numbers and for small accretors, we also find evidence of the advective-acoustic instability operating in the adiabatic case, and of an instability developing around the wake's stagnation point in the quasi-isothermal case. Moreover, at very high resolution, and Mach numbers above $\mathcal{M}_\infty \geq 3$, the flow behind the accretion bow shock becomes entirely dominated by these instabilities. As a result, accretion rates onto the black hole drop by about an order of magnitude in the adiabatic case, compared to the analytic BHL formula., Comment: 24 pages, accepted for publication in MNRAS

Accurate cosmology from upcoming weak lensing surveys relies on knowledge of the total matter power spectrum at percent level at scales $k < 10$ $h$/Mpc, for which modelling the impact of baryonic physics is crucial. We compare measurements of the total matter power spectrum from the Horizon cosmological hydrodynamical simulations: a dark matter-only run, one with full baryonic physics, and another lacking Active Galactic Nuclei (AGN) feedback. Baryons cause a suppression of power at $k\simeq 10$ $h/$Mpc of $<15\%$ at $z=0$, and an enhancement of a factor of a few at smaller scales due to the more efficient cooling and star formation. The results are sensitive to the presence of the highest mass haloes in the simulation and the distribution of dark matter is also impacted up to a few percent. The redshift evolution of the effect is non-monotonic throughout $z=0-5$ due to an interplay between AGN feedback and gas pressure, and the growth of structure. We investigate the effectiveness of an analytic `baryonic correction model' in describing our results. We require a different redshift evolution and propose an alternative fitting function with $4$ free parameters that reproduces our results within $5\%$. Compared to other simulations, we find the impact of baryonic processes on the total matter power spectrum to be smaller at $z=0$. Correspondingly, our results suggest that AGN feedback is not strong enough in the simulation. Total matter power spectra from the Horizon simulations are made publicly available at https://www.horizon-simulation.org/catalogues.html, Comment: 19 pages, 14 figures, matches published version

Intrinsic alignments of galaxies are a significant astrophysical systematic affecting cosmological constraints from weak gravitational lensing. Obtaining numerical predictions from hydrodynamical simulations of expected survey volumes is expensive, and a cheaper alternative relies on populating large dark matter-only simulations with accurate models of alignments calibrated on smaller hydrodynamical runs. This requires connecting the shapes and orientations of galaxies to those of dark matter haloes and to the large-scale structure. In this paper, we characterise galaxy-halo alignments in the Horizon-AGN cosmological hydrodynamical simulation. We compare the shapes and orientations of galaxies in the redshift range $0