We present the Extreme-Horizon (EH) cosmological simulation: EH models galaxy formation with stellar and AGN feedback and uses a very high resolution in the intergalactic and circumgalactic medium. The high resolution in low-density regions results in smaller-size massive galaxies at redshift $z=2$, in better agreement with observations compared to other simulations. This results from the improved modeling of cold gas flows accreting onto galaxies. Besides, the EH simulation forms a population of particularly compact galaxies with stellar masses of $10^{10-11}$\,M$_\sun$ that are reminiscent of observed ultracompact galaxies at $z\simeq2$. These objects form mainly through repeated major mergers of low-mass progenitors, independently of baryonic feedback mechanisms. This formation process can be missed in simulations using a too low resolution in low-density intergalactic regions.

Accretion of stars on massive black holes (MBHs) can feed MBHs and generate tidal disruption events (TDEs). We introduce a new physically motivated model to self-consistently treat TDEs in cosmological simulations, and apply it to the assembly of a galaxy with final mass $3\times 10^{10}\mathrm{\, M}_\odot$ at $z=6$. This galaxy exhibits a TDE rate of $\sim 10^{-5}\mathrm{\, yr}^{-1}$, consistent with local observations but already in place when the Universe was one billion year old. A fraction of the disrupted stars participate in the growth of MBHs, dominating it until the MBH reaches mass $\sim 5 \times 10^5 \mathrm{\, M}_\odot$, but their contribution then becomes negligible compared to gas. TDEs could be a viable mechanism to grow light MBH seeds, but fewer TDEs are expected when the MBH becomes sufficiently massive to reach the luminosity of, and be detected as, an active galactic nucleus. Galaxy mergers bring multiple MBHs in the galaxy, resulting in an enhancement of the global TDE rate in the galaxy by $\sim 1$ order of magnitude during $100\mathrm{\, Myr}$ around mergers. This enhancement is not on the central MBH, but caused by the presence of MBHs in the infalling galaxies. This is the first self consistent study of TDEs in a cosmological environment and highlights that accretion of stars and TDEs are a natural process occurring in a Milky~Way-mass galaxy at early cosmic times.

Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large volume and low resolution which tracks the high-mass (> 1e7 Msun) MBH population, and NewHorizon with a smaller volume but higher resolution that traces the low-mass (< 1e7 Msun) MBH population. While Horizon-AGN can be used to estimate the rate of inspirals for Pulsar Timing Arrays, NewHorizon can investigate MBH mergers in a statistical sample of dwarf galaxies for LISA, which is sensitive to low-mass MBHs. After the numerical MBH merger at the resolution limit, we post-process MBH dynamics to account for time delays mostly determined by dynamical friction and stellar hardening. In both simulations, MBHs typically merge long after the galaxies do, so that the galaxy morphology at the time of the MBH merger is no longer determined by the galaxy merger from which the MBH merger originated. These time delays cause a loss of high-z MBH coalescences, shifting the peak of the MBH merger rate to z~1-2. This study shows how tracking MBH mergers in low-mass galaxies is crucial to probing the MBH merger rate for LISA and investigate the properties of the host galaxies.

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first 1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass halos. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant and scaling with stellar mass) from both star formation ($\langle f_{\rm esc}^{\rm sf} \rangle$) and AGN ($f_{\rm esc}^{\rm bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_*<10^9 {\rm M_\odot}$ ) galaxies at $z>6$ with the AGN contribution (driven by $M_{bh}>10^6 {\rm M_\odot}$ black holes in $M_* > 10^9 {\rm M_\odot}$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25\%$ to the cumulative ionizing emissivity by $z=4$ for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of $\langle f_{\rm esc}^{\rm sf} \rangle$ is shallower than $f_{\rm esc}^{\rm bh}$, at $z<7$ a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at $z \sim 6.8$, they take-over at stellar masses that are a tenth of the knee mass by $z=4$.

Recent deep surveys have unravelled a population of faint active galactic nuclei (AGN) in the high redshift Universe, leading to various discussions on their nature and their role during the Epoch of Reionization. We use cosmological radiation-hydrodynamics simulations of a bright galaxy at z = 6 ($M_\star \gtrsim 10^{10} M_\odot$) hosting an actively growing super-massive black hole to study the properties of these objects. In particular, we study how the black hole and the galaxy co-evolve and what is the relative contribution of the AGN and of the stellar populations to the luminosity budget of the system. We find that the feedback from the AGN has no strong effect on the properties of the galaxy, and does not increase the total ionizing luminosity of the host. The average escape fraction of our galaxy is around $f_{\rm esc} \sim 5\%$. While our galaxy would be selected as an AGN in deep X-ray surveys, most of the UV luminosity is originating from stellar populations. This confirms that there is a transition in the galaxy population from star forming galaxies to quasar hosts, with bright Lyman-Break Galaxies (LBGs) with $M_{\rm UV}$ around -22 falling in the overlap region. Our results also suggest that faint AGN do not contribute significantly to reionizing the Universe.

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$.

Resonant lines are powerful probes of the interstellar and circumgalactic medium of galaxies. Their transfer in gas being a complex process, the interpretation of their observational signatures, either in absorption or in emission, is often not straightforward. Numerical radiative transfer simulations are needed to accurately describe the travel of resonant line photons in real and in frequency space, and to produce realistic mock observations. This paper introduces RASCAS, a new public 3D radiative transfer code developed to perform the propagation of any resonant line in numerical simulations of astrophysical objects. RASCAS was designed to be easily customisable and to process simulations of arbitrarily large sizes on large supercomputers. RASCAS performs radiative transfer on an adaptive mesh with an octree structure using the Monte Carlo technique. RASCAS features full MPI parallelisation, domain decomposition, adaptive load-balancing, and a standard peeling algorithm to construct mock observations. The radiative transport of resonant line photons through different mixes of species (e.g. \ion{H}{i}, \ion{Si}{ii}, \ion{Mg}{ii}, \ion{Fe}{ii}), including their interaction with dust, is implemented in a modular fashion to allow new transitions to be easily added to the code. RASCAS is very accurate and efficient. It shows perfect scaling up to a minimum of a thousand cores. It has been fully tested against radiative transfer problems with analytic solutions and against various test cases proposed in the literature. Although it was designed to describe accurately the many scatterings of line photons, RASCAS may also be used to propagate photons at any wavelength (e.g. stellar continuum or fluorescent lines), or to cast millions of rays to integrate the optical depths of ionising photons, making it highly versatile.

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.

Dense gas in the centre of galaxies feeds massive black holes, but can also become a source of obscuration and limit our ability to find faint Active Galactic Nuclei (AGN). We use a high resolution cosmological radiative hydrodynamics simulation to connect the properties of the gas in the central region (a few tens of parsecs) of a high redshift galaxy to the growth of a massive black hole during the first billion years of the Universe. We find that the feedback from the AGN efficiently controls the growth of the black hole and limits the duration of the high accretion episodes by emptying the gas reservoir. As the galaxy grows in mass, the production of metals results in the presence of dust-enriched gas in the galaxy centre that can obscure highly accreting black holes enough to strongly reduce their UV/optical visibility. We also find that the gas outside the very centre of the galaxy can contribute to the total column density and obscuration at a level at least comparable to the gas in the nuclear region. We suggest that this explains the different duty cycles required to explain the masses of high redshift quasars and the observed UV/optical luminosity functions: in our case, the AGN would be observed with an X-ray luminosity above $L_X = 10^{42}$ erg/s around 30% of the time, but with UV magnitude brighter than $M_{1450}$ = -23 only 4% of the time.

While low mass, star forming galaxies are often considered as the primary driver of reionization, their actual contribution to the cosmic ultraviolet background is still uncertain, mostly because the escape fraction of ionizing photons is only poorly constrained. Theoretical studies have shown that efficient supernova feedback is a necessary condition to create paths through which ionizing radiation can escape into the intergalactic medium. We investigate the possibility that accreting supermassive black holes in early dwarf galaxies may provide additional feedback and enhance the leakage of ionizing radiation. We use a series of high resolution cosmological radiation hydrodynamics simulations where we isolate the different sources of feedback. We find that supernova feedback prevents the growth of the black hole, thus quenching its associated feedback. Even in cases where the black hole can grow, the structure of the interstellar medium is strongly dominated by supernova feedback. We conclude that, in the dwarf galaxy regime, supermassive black holes do not appear to play a significant role in enhancing the escape fraction and in contributing to the early UV background.

The first billion years of the Universe is a pivotal time: stars, black holes (BHs) and galaxies form and assemble, sowing the seeds of galaxies as we know them today. Detecting, identifying and understand- ing the first galaxies and BHs is one of the current observational and theoretical challenges in galaxy formation. In this paper we present a population synthesis model aimed at galaxies, BHs and Active Galactic Nuclei (AGNs) at high redshift. The model builds a population based on empirical relations. Galaxies are characterized by a spectral energy distribution determined by age and metallicity, and AGNs by a spectral energy distribution determined by BH mass and accretion rate. We validate the model against observational constraints, and then predict properties of galaxies and AGN in other wavelength and/or luminosity ranges, estimating the contamination of stellar populations (normal stars and high-mass X-ray binaries) for AGN searches from the infrared to X-rays, and vice-versa for galaxy searches. For high-redshift galaxies, with stellar ages < 1 Gyr, we find that disentangling stellar and AGN emission is challenging at restframe UV/optical wavelengths, while high-mass X-ray binaries become more important sources of confusion in X-rays. We propose a color-color selection in JWST bands to separate AGN vs star-dominated galaxies in photometric observations. We also esti- mate the AGN contribution, with respect to massive, hot, metal-poor stars, at driving high ionization lines, such as C IV and He II. Finally, we test the influence of the minimum BH mass and occupa- tion fraction of BHs in low mass galaxies on the restframe UV/near-IR and X-ray AGN luminosity function.

Low mass galaxies are thought to provide the bulk of the ionizing radiation necessary to reionize the Universe. The amount of photons escaping the galaxies is poorly constrained theoretically, and difficult to measure observationally. Yet it is an essential parameter of reionization models. We study in detail how ionizing radiation can leak from high redshift galaxies. For this purpose, we use a series of high resolution radiation hydrodynamics simulations, zooming on three dwarf galaxies in a cosmological context. We find that the energy and momentum input from the supernova explosions has a pivotal role in regulating the escape fraction, by disrupting dense star forming clumps, and clearing sight lines in the halo. In the absence of supernovae, photons are absorbed very locally, within the birth clouds of massive stars. We follow the time evolution of the escape fraction, and find that it can vary by more than six orders of magnitude. This explains the large scatter in the value of the escape fraction found by previous studies. This fast variability also impacts the observability of the sources of reionization: a survey even as deep as $M_{\rm UV} = -14$ would miss about half of the underlying population of Lyman-continuum emitters.

Lyman-$\alpha$ nebulae are usually found in massive environments at high redshift ($z > 2$). The origin of their Lyman-$\alpha$ (Lya) emission remains debated. Recent polarimetric observations showed that at least some Lya sources are polarized. This is often interpreted as a proof that the photons are centrally produced, and opposed to the scenario in which the Lya emission is the cooling radiation emitted by gas heated during the accretion onto the halo. We suggest that this scenario is not incompatible with the polarimetric observations. In order to test this idea, we post-process a radiative hydrodynamics simulation of a blob with the MCLya Monte Carlo transfer code. We compute radial profiles for the surface brightness and the degree of polarization and compare them to existing observations. We find that both are consistent with a significant contribution of the extragalactic gas to the Lya emission. Most of the photons are centrally emitted and scattered inside the filament afterwards, producing the observed high level of polarization. We argue that the contribution of the extragalactic gas to the Lya emission does not prevent polarization to arise. On the contrary, we find that pure galactic emission causes the polarization profile to be too steep to be consistent with observations.

Small galaxies are thought to be the main contributors to the ionising budget of the Universe before reionisation was complete. There have been a number of numerical studies trying to quantify their ionising efficiency through the escape fraction $f_{esc}$. While there is a clear trend that $f_{esc}$ is higher for smaller haloes, there is a large scatter in the distribution of $f_{esc}$ for a single halo mass. We propose that this is due to the intrinsic burstiness of star formation in low mass galaxies. We performed high resolution radiative hydrodynamics simulations with Ramses-RT to model the evolution of three galaxies and their ionising efficiency. We found that the variability of $f_{esc}$ follows that of the star formation rate. We then discuss the consequences of this variability on the observability of such galaxies by JWST.

A generic prediction of string theory is the existence of many axion fields. It has recently been argued that many of these fields should be light and, like the well known QCD axion, lead to observable cosmological consequences. In this paper we study in detail the effect of the so-called string axiverse on large scale structure, focusing on the morphology and evolution of density perturbations, anisotropies in the cosmic microwave background and weak gravitational lensing of distant galaxies. We quantify specific effects that will arise from the presence of the axionic fields and highlight possible degeneracies that may arise in the presence of massive neutrinos. We take particular care understanding the different physical effects and scales that come into play. We then forecast how the string axiverse may be constrained and show that with a combination of different observations, it should be possible to detect a fraction of ultralight axions to dark matter of a few percent.