A three-phase amplification of the cosmic magnetic field in galaxies
Monthly Notices of the Royal Astronomical Society Oxford University Press 479:3 (2018) 3343-3365
Abstract:
Arguably the main challenge of galactic magnetism studies is to explain how the interstellar medium of galaxies reaches energetic equipartition despite the extremely weak cosmic primordial magnetic fields that are originally predicted to thread the inter-galactic medium. Previous numerical studies of isolated galaxies suggest that a fast dynamo amplification might suffice to bridge the gap spanning many orders of magnitude in strength between the weak early Universe magnetic fields and the ones observed in high redshift galaxies. To better understand their evolution in the cosmological context of hierarchical galaxy growth, we probe the amplification process undergone by the cosmic magnetic field within a spiral galaxy to unprecedented accuracy by means of a suite of constrained transport magnetohydrodynamical adaptive mesh refinement cosmological zoom simulations with different stellar feedback prescriptions. A galactic turbulent dynamo is found to be naturally excited in this cosmological environment, being responsible for most of the amplification of the magnetic energy. Indeed, we find that the magnetic energy spectra of simulated galaxies display telltale inverse cascades. Overall, the amplification process can be divided in three main phases, which are related to different physical mechanisms driving galaxy evolution: an initial collapse phase, an accretion-driven phase, and a feedback-driven phase. While different feedback models affect the magnetic field amplification differently, all tested models prove to be subdominant at early epochs, before the feedback-driven phase is reached. Thus the three-phase evolution paradigm is found to be quite robust vis-a-vis feedback prescriptions.Bondi or not Bondi: The impact of resolution on accretion and drag force modelling for Supermassive Black Holes
Monthly Notices of the Royal Astronomical Society Oxford University Press 478:1 (2018) 995-1016
Abstract:
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.Gas flows in the circumgalactic medium around simulated high-redshift galaxies
Monthly Notices of the Royal Astronomical Society Oxford University Press 474:4 (2017) 4279-4301
Abstract:
We analyse the properties of circumgalactic gas around simulated galaxies in the redshift range z ≥ 3, utilizing a new sample of cosmological zoom simulations. These simulations are intended to be representative of the observed samples of Lyman α (Ly α) emitters recently obtained with the multi unit spectroscopic explorer (MUSE) instrument (halo masses ~ 10 10 - 10 11 M⊙). We show that supernova feedback has a significant impact on both the inflowing and outflowing circumgalactic medium (CGM) by driving outflows, reducing diffuse inflow rates, and by increasing the neutral fraction of inflowing gas. By temporally stacking simulation outputs, we find that significant net mass exchange occurs between inflowing and outflowing phases: none of the phases are mass-conserving. In particular, we find that the mass in neutral outflowing hydrogen declines exponentially with radius as gas flows outwards from the halo centre. This is likely caused by a combination of both fountain-like cycling processes and gradual photoionization/collisional ionization of outflowing gas. Our simulations do not predict the presence of fast-moving neutral outflows in the CGM. Neutral outflows instead move with modest radial velocities (~ 50 km s -1 ), and the majority of the kinetic energy is associated with tangential rather than radial motion.Cosmic evolution of stellar quenching by AGN feedback: clues from the Horizon-AGN simulation
Monthly Notices of the Royal Astronomical Society Oxford University Press 472:1 (2017) 949-965
Abstract:
The observed massive end of the local galaxy stellar mass function is steeper than its predicted dark matter (DM) halo counterpart in the standard $\Lambda $CDM paradigm. We investigate how active galactic nuclei (AGN) feedback can account for such a reduction in the stellar content of massive galaxies, through an influence on the gas content of their interstellar (ISM) and circum-galactic medium (CGM). We isolate the impact of AGNs by comparing two simulations from the HORIZON suite, which are identical except that one includes super massive black holes (SMBH) and related feedback. This allows us to cross-identify individual galaxies between these simulations and quantify the effect of AGN feedback on their properties, such as stellar mass and gas outflows. We find that the most massive galaxies ($ \rm M_{*} \geq 3 \times 10^{11} M_\odot $) are quenched to the extent that their stellar masses decrease by about 80% at $z=0$. More generally, SMBHs affect their host halo through a combination of outflows that reduce their baryonic mass, particularly for galaxies in the mass range $ \rm 10^9 M_\odot \leq M_{*} \leq 10^{11} M_\odot $, and a disruption of central gas inflows, which limits in-situ star formation, particularly massive galaxies with $ \rm M_{*} \approx10^{11} M_\odot $. As a result of these processes, net gas inflows onto massive galaxies drop by up to 70%. Finally, we measure a redshift evolution in the stellar mass ratio of twin galaxies with and without AGN feedback, with galaxies of a given stellar mass showing stronger signs of quenching earlier on. This evolution is driven by a progressive flattening of the $\rm M_{SMBH}-M_* $ relation for galaxies with $\rm M_{*} \leq 10^{10} M_\odot $ as redshift decreases, which translates into smaller SBMHs being harboured by galaxies of any fixed stellar mass, and indicates stronger AGN feedback at higher redshift.Fluctuating feedback-regulated escape fraction of ionizing radiation in low-mass, high-redshift galaxies
Monthly Notices of the Royal Astronomical Society Oxford University Press 470:1 (2017) 224-239