Zooming in on supermassive black holes: how resolving their gas cloud host renders their accretion episodic

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 page

Building Merger Trees from Cosmological N-body Simulations

Although a fair amount of work has been devoted to growing Monte-Carlo merger trees which resemble those built from an N-body simulation, comparatively little effort has been invested in quantifying the caveats one necessarily encounters when one extracts trees directly from such a simulation. To somewhat revert the tide, this paper seeks to provide its reader with a comprehensive study of the problems one faces when following this route. The first step to building merger histories of dark matter haloes and their subhaloes is to identify these structures in each of the time outputs (snapshots) produced by the simulation. Even though we discuss a particular implementation of such an algorithm (called AdaptaHOP) in this paper, we believe that our results do not depend on the exact details of the implementation but extend to most if not all (sub)structure finders. We then highlight different ways to build merger histories from AdaptaHOP haloes and subhaloes, contrasting their various advantages and drawbacks. We find that the best approach to (sub)halo merging histories is through an analysis that goes back and forth between identification and tree building rather than one which conducts a straightforward sequential treatment of these two steps. This is rooted in the complexity of the merging trees which have to depict an inherently dynamical process from the partial temporal information contained in the collection of instantaneous snapshots available from the N-body simulation.Comment: 19 pages, 28 figure

Probing Cosmic Dawn: Modelling the Assembly History, SEDs, and Dust Content of Selected z∼9z\sim9 Galaxies

The presence of spectroscopically confirmed Balmer breaks in galaxy spectral energy distributions (SEDs) at $z>9$ provides one of the best probes of the assembly history of the first generations of stars in our Universe. Recent observations of the gravitationally lensed source, MACS 1149_JD1 (JD1), indicate that significant amounts of star formation likely occurred at redshifts as high as $z\simeq15$. The inferred stellar mass, dust mass, and assembly history of JD1, or any other galaxy at these redshifts that exhibits a strong Balmer break, can provide a strong test of our best theoretical models from high-resolution cosmological simulations. In this work, we present the results from a cosmological radiation-hydrodynamics simulation of the region surrounding a massive Lyman-break galaxy. For two of our most massive systems, we show that dust preferentially resides in the vicinity of the young stars thereby increasing the strength of the measured Balmer break such that the simulated SEDs are consistent with the photometry of JD1 and two other $z>9$ systems (GN-z10-3 and GN-z9-1) that have proposed Balmer breaks at high redshift. We find strong variations in the shape and luminosity of the SEDs of galaxies with nearly identical stellar and halo masses, indicating the importance of morphology, assembly history, and dust distribution in making inferences on the properties of individual galaxies at high redshifts. Our results stress the importance that dust may play in modulating the observable properties of galaxies, even at the extreme redshifts of $z>9$.Comment: 16 pages, 13 Figures, Accepted to MNRA

Reionization history constraints from neural network based predictions of high-redshift quasar continua

Observations of the early Universe suggest that reionization was complete by $z\sim6$, however, the exact history of this process is still unknown. One method for measuring the evolution of the neutral fraction throughout this epoch is via observing the Ly$\alpha$ damping wings of high-redshift quasars. In order to constrain the neutral fraction from quasar observations, one needs an accurate model of the quasar spectrum around Ly$\alpha$, after the spectrum has been processed by its host galaxy but before it is altered by absorption and damping in the intervening IGM. In this paper, we present a novel machine learning approach, using artificial neural networks, to reconstruct quasar continua around Ly$\alpha$. Our QSANNdRA algorithm improves the error in this reconstruction compared to the state-of-the-art PCA-based model in the literature by 14.2% on average, and provides an improvement of 6.1% on average when compared to an extension thereof. In comparison with the extended PCA model, QSANNdRA further achieves an improvement of 22.1% and 16.8% when evaluated on low-redshift quasars most similar to the two high-redshift quasars under consideration, ULAS J1120+0641 at $z=7.0851$ and ULAS J1342+0928 at $z=7.5413$, respectively. Using our more accurate reconstructions of these two $z>7$ quasars, we estimate the neutral fraction of the IGM using a homogeneous reionization model and find $\bar{x}_\mathrm{HI} = 0.25^{+0.05}_{-0.05}$ at $z=7.0851$ and $\bar{x}_\mathrm{HI} = 0.60^{+0.11}_{-0.11}$ at $z=7.5413$. Our results are consistent with the literature and favour a rapid end to reionization

Massive spheroids can form in single minor mergers

Accepted for publication in MNRAS, 12 pages, 6 figuresUnderstanding how rotationally supported discs transform into dispersion-dominated spheroids is central to our comprehension of galaxy evolution. Morphological transformation is largely merger-driven. While major mergers can efficiently create spheroids, recent work has highlighted the significant role of other processes, like minor mergers, in driving morphological change. Given their rich merger histories, spheroids typically exhibit large fractions of ‘ex situ’ stellar mass, i.e. mass that is accreted, via mergers, from external objects. This is particularly true for the most massive galaxies, whose stellar masses typically cannot be attained without a large number of mergers. Here, we explore an unusual population of extremely massive (M ∗ > 10 11M) spheroids, in the Horizon-AGN simulation, which exhibit anomalously low ex situ mass fractions, indicating that they form without recourse to significant merging. These systems form in a single minor-merger event (with typical merger mass ratios of 0.11–0.33), with a specific orbital configuration, where the satellite orbit is virtually co-planar with the disc of the massive galaxy. The merger triggers a catastrophic change in morphology, over only a few hundred Myr, coupled with strong in situ star formation. While this channel produces a minority (∼5 per cent) of such galaxies, our study demonstrates that the formation of at least some of the most massive spheroids need not involve major mergers – or any significant merging at all – contrary to what is classically believed.Peer reviewedFinal Accepted Versio

Simulations of Dust in Interacting Galaxies I: Dust Attenuation

A new Monte-Carlo radiative-transfer code, Sunrise, is used in conjunction with hydrodynamic simulations of major galaxy mergers to calculate the effects of dust in such systems. The simulations are in good agreement with observations of dust absorption in starburst galaxies, and the dust has a profound effect on their appearance. The dust attenuation increases with luminosity such that at peak luminosities ~90% of the bolometric luminosity is absorbed by dust. In general, the detailed appearance of the merging event depends on the stage of the merger and the geometry of the encounter. The fraction of bolometric energy absorbed by the dust, however, is a robust quantity that can be predicted from the intrinsic properties bolometric luminosity, baryonic mass, star-formation rate, and metallicity of the system. This paper presents fitting formulae, valid over a wide range of masses and metallicities, from which the absorbed fraction of luminosity (and consequently also the infrared dust luminosity) can be predicted. The attenuation of the luminosity at specific wavelengths can also be predicted, albeit with a larger scatter due to the variation with viewing angle. These formulae for dust attenuation appear to be valid for both isolated and interacting galaxies, are consistent with earlier studies, and would be suitable for inclusion in theoretical models, e.g. semi-analytic models of galaxy formation.Comment: 12 pages, 10 figures, submitted to Ap

The elliptical galaxy colour-magnitude relation as a discriminant between the monolithic and merger paradigms

The colour-magnitude relation (CMR) of cluster elliptical galaxies has been widely used to constrain their star formation histories (SFHs) and to discriminate between the monolithic collapse and merger paradigms of elliptical galaxy formation. We use a Λ cold dark matter hierarchical merger model of galaxy formation to investigate the existence and redshift evolution of the elliptical galaxy CMR in the merger paradigm. We show that the SFH of cluster ellipticals predicted by the model is quasi-monolithic, with only ∼10 per cent of the total stellar mass forming after z∼ 1. The quasi-monolithic SFH results in a predicted CMR that agrees well with its observed counterpart in the redshift range 0 < z < 1.27. We use our analysis to argue that the elliptical-only CMR can be used to constrain the SFHs of present-day cluster ellipticals only if we believe a priori in the monolithic collapse model. It is not a meaningful tool for constraining the SFH in the merger paradigm, since a progressively larger fraction of the progenitor set of present-day cluster ellipticals is contained in late-type star-forming systems at higher redshift, which cannot be ignored when deriving the SFHs. Hence, the elliptical-only CMR is not a useful discriminant between the two competing theories of elliptical galaxy evolutio