We use the Simba cosmological hydrodynamic simulation suite to explore the impact of feedback on the circumgalactic medium (CGM) and intergalactic medium (IGM) around $2 \leq z \leq 3$ quasars. We identify quasars in Simba as the most rapidly-accreting black holes, and show that they are well-matched in bolometric luminosity and correlation strength to real quasars. We extract Lyman-alpha (Ly-a) absorption in spectra passing at different transverse distances (10 kpc $\lesssim b \lesssim$ 10 Mpc) around those quasars, and compare to observations of the mean Ly-a absorption profile. The observations are well reproduced, except within 100 kpc from the foreground quasar, where Simba overproduces absorption; this could potentially be mitigated by including ionisation from the quasar itself. By comparing runs with different feedback modules activated, we find that (mechanical) AGN feedback has little impact on the surrounding CGM even around these most highly luminous black holes, while stellar feedback has a significant impact. By further investigating thermodynamic and kinematic properties of CGM gas, we find that stellar feedback, and not AGN feedback, is the primary physical driver in determining the average properties of the CGM around $z\sim 2-3$ quasars. We also compare our results with previous works, and find that Simba predicts much more absorption within 100 kpc than the Nyx and Illustris simulations, showing that the Ly-a absorption profile can be a powerful constraint on simulations. Instruments such as VLT-MUSE and upcoming surveys (e.g., WEAVE and DESI) promise to further improve such constraints.

We examine the impact of black hole jet feedback on the properties of the low-redshift intergalactic medium (IGM) in the Simba simulation, with a focus on the Ly$\alpha$ forest mean flux decrement $D_A$ and the inferred HI photo-ionisation rate $\Gamma_{HI}$. Without jet feedback, we confirm the Photon Underproduction Crisis (PUC; Kollmeier et al. 2014) in which $\Gamma_{HI}$ at $z=0$ must be increased by $\times 6$ over the Haardt & Madau (2012) value in order to match the observed $D_A$. Turning on jet feedback lowers this discrepancy to $\sim\times 2.5$, and additionally using the recent Faucher-Gigu\`ere 2019 background results in even better agreement, nearly solving the PUC. The PUC becomes apparent at late epochs ($z \sim < 1$) where the jet and no-jet simulations diverge; at higher redshifts Simba reproduces the observed $D_A$ with no adjustment, with or without jets. The main impact of jet feedback is to lower the cosmic baryon fraction in the diffuse IGM from 39% to 16% at $z=0$, while increasing the warm-hot intergalactic medium (WHIM) baryon fraction from 30% to 70%; the lowering of the diffuse IGM content directly translates into a lowering of $D_A$ by a similar factor. Comparing to the older Mufasa simulation that employs different quenching feedback but is otherwise similar to Simba, Mufasa matches $D_A$ less well than Simba, suggesting that low-redshift measurements of $D_A$ and $\Gamma_{HI}$ could provide constraints on feedback mechanisms. Our results suggest that widespread IGM heating at late times is a plausible solution to the PUC, and that Simba's jet AGN feedback model, constrained to reproduce quenched massive galaxies, approximately yields this required heating.

In this work we present a new hybrid method to simulate the thermal effects of the reionization in cosmological hydrodynamical simulations. The method improves upon the standard approach used in simulations of the intergalactic medium (IGM) and galaxy formation without a significant increase of the computational cost allowing for efficient exploration of the parameter space. The method uses a small set of phenomenological input parameters and combines a semi-numerical reionization model to solve for the topology of reionization and an approximate model of how reionization heats the IGM, with the massively parallel \texttt{Nyx} hydrodynamics code, specifically designed to solve for the structure of diffuse IGM gas. We have produced several large-scale high resolution cosmological hydrodynamical simulations ($2048^3$, $L_{\rm box} = 40$ Mpc/h) with different instantaneous and inhomogeneous HI reionization models that use this new methodology. We study the IGM thermal properties of these models and find that large scale temperature fluctuations extend well beyond the end of reionization. Analyzing the 1D flux power spectrum of these models, we find up to $\sim 50\%$ differences in the large scale properties (low modes, $k\lesssim0.01$ s/km) of the post-reionization power spectrum due to the thermal fluctuations. We show that these differences could allow one to distinguish between different reionization scenarios already with existing Ly$\alpha$ forest measurements. Finally, we explore the differences in the small-scale cutoff of the power spectrum and we find that, for the same heat input, models show very good agreement provided that the reionization redshift of the instantaneous reionization model happens at the midpoint of the inhomogeneous model.

Galaxy formation depends critically on the physical state of gas in the circumgalactic medium (CGM) and its interface with the intergalactic medium (IGM), determined by the complex interplay between inflows from the IGM and outflows from supernovae or AGN feedback. The average Lyman-alpha (Ly-a) absorption profile around galactic halos represents a powerful tool to probe their gaseous environments. We compare predictions from Illustris and Nyx hydrodynamical simulations with the observed absorption around foreground quasars, damped Ly-a systems, and Lyman-break galaxies. We show how large-scale BOSS and small-scale quasar pair measurements can be combined to precisely constrain the absorption profile over three decades in transverse distance 20kpc$\lesssim b\lesssim$20Mpc. Far from galaxies $\gtrsim2$Mpc, the simulations converge to the same profile and provide a reasonable match to the observations. This asymptotic agreement arises because the $\Lambda$CDM model successfully describes the ambient IGM, and represents a critical advantage of studying the mean absorption profile. However, significant differences between the simulations, and between simulations and observations are present on scales 20kpc$\lesssim b\lesssim$2Mpc, illustrating the challenges of accurately modeling and resolving galaxy formation physics. It is noteworthy that these differences are observed as far out as $\sim2$Mpc, indicating that the `sphere-of-influence' of galaxies could extend to approximately $\sim20$ times the halo virial radius ($\sim100$kpc). Current observations are very precise on these scales and can thus strongly discriminate between different galaxy formation models. We demonstrate that the Ly-a absorption profile is primarily sensitive to the underlying temperature-density relationship of diffuse gas around galaxies, and argue that it thus provides a fundamental test of galaxy formation models.

Measuring the clustering of galaxies from surveys allows us to estimate the power spectrum of matter density fluctuations, thus constraining cosmological models. This requires careful modelling of observational effects to avoid misinterpretation of data. In particular, signals coming from different distances encode information from different epochs. This is known as "light-cone effect" and is going to have a higher impact as upcoming galaxy surveys probe larger redshift ranges. Generalising the method by Feldman et al. (1994), I define a minimum-variance estimator of the linear power spectrum at a fixed time, properly taking into account the light-cone effect. An analytic expression for the estimator is provided, and that is consistent with the findings of previous works in the literature. I test the method within the context of the halo model, assuming Planck 2014 cosmological parameters. I show that the estimator presented recovers the fiducial linear power spectrum at present time within 5% accuracy up to $k \sim 0.80\;h\,\mathrm{Mpc}^{-1}$ and within 10% up to $k \sim 0.94\;h\,\mathrm{Mpc}^{-1}$, well into the non-linear regime of the growth of density perturbations. As such, the method could be useful in the analysis of the data from future large-scale surveys, like Euclid.

Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we present Iteratively Matched Statistics (IMS), a novel method to accurately model the Lyman-alpha forest with collisionless N-body simulations, where the relevant density fluctuations are unresolved. We use a small-box, high-resolution hydrodynamic simulation to obtain the probability distribution function (PDF) and the power spectrum of the real-space Lyman-alpha forest flux. These two statistics are iteratively mapped onto a pseudo-flux field of an N-body simulation, which we construct from the matter density. We demonstrate that our method can perfectly reproduce line-of-sight observables, such as the PDF and power spectrum, and accurately reproduce the 3D flux power spectrum (5-20%). We quantify the performance of the commonly used Gaussian smoothing technique and show that it has significantly lower accuracy (20-80%), especially for N-body simulations with achievable mean inter-particle separations in large-volume simulations. In addition, we show that IMS produces reasonable and smooth spectra, making it a powerful tool for modeling the IGM in large cosmological volumes and for producing realistic "mock" skies for Lyman-alpha forest surveys.