The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter density slope, $\alpha$, on the ratio between stellar-to-halo mass. We show that this relation holds approximately at all redshifts (with an intrinsic scatter of ~0.18 in $\alpha$). This implies that in practically all haloes the shape of their inner density profile changes quite substantially over cosmic time, as they grow in stellar and total mass. Thus, depending on their final stellar-to-halo mass ratio, haloes can either form and keep a substantial density core (size~1 kpc), or form and then destroy the core and re-contract the halo, going back to a cuspy profile, which is even steeper than CDM predictions for massive galaxies (~1e12 Msun). We show that results from the NIHAO suite are in good agreement with recent observational measurements of $\alpha$ in dwarf galaxies. Overall our results suggest that the notion of a universal density profile for dark matter haloes is no longer valid in the presence of galaxy formation.

Observational and theoretical arguments increasingly suggest that the initial mass function (IMF) of stars may depend systematically on environment, yet most galaxy formation models to date assume a universal IMF. Here we investigate simulations of the formation of Milky Way analogues run with an empirically derived metallicity-dependent IMF and the moving-mesh code AREPO in order to characterize the associated uncertainties. In particular, we compare a constant Chabrier and a varying metallicity-dependent IMF in cosmological, magneto-hydrodynamical zoom-in simulations of Milky Way-sized halos. We find that the non-linear effects due to IMF variations typically have a limited impact on the morphology and the star formation histories of the formed galaxies. Our results support the view that constraints on stellar-to-halo mass ratios, feedback strength, metallicity evolution and metallicity distributions are in part degenerate with the effects of a non-universal, metallicity-dependent IMF. Interestingly, the empirical relation we use between metallicity and the high mass slope of the IMF does not aid in the quenching process. It actually produces up to a factor of 2-3 more stellar mass if feedback is kept constant. Additionally, the enrichment history and the z = 0 metallicity distribution are significantly affected. In particular, the alpha enhancement pattern shows a steeper dependence on iron abundance in the metallicity-dependent model, in better agreement with observational constraints.

To fully harvest the rich library of stellar elemental abundance data available, we require reliable models that facilitate our interpretation of them. Galactic chemical evolution (GCE) models are one such set, and a key part of which are the selection of chemical yields from different nucleosynthetic enrichment channels, predominantly asymptotic giant branch (AGB) stars, Type Ia supernovae (SNe Ia), and core-collapse supernovae (CC-SNe). Here, we present a scoring system for yield tables based on their ability to reproduce proto-solar abundances within a simple parametrisation of the GCE modelling software Chempy, which marginalises over galactic parameters describing simple stellar populations (SSPs) and interstellar medium physics. Two statistical scoring methods are presented, based on Bayesian evidence and leave-one-out cross-validation and are applied to five CC-SN tables, (a) for all mutually available elements and (b) for a subset of the 9 most abundant elements. We find that the yields of Prantzos et al. (P18, including stellar rotation) and Chieffi & Limongi (C04) best reproduce proto-solar abundances for the two cases, respectively. The inferred best-fit SSP parameters for (b) are $\alpha_\mathrm{IMF}=-2.45^{+0.15}_{-0.11}$ for the initial mass function high-mass slope and $\mathrm{N}_\mathrm{Ia}=1.29^{+0.45}_{-0.31}\times10^{-3}$ M$_\odot^{-1}$ for the SN Ia normalisation, which are broadly consistent across tested yield tables. Additionally, we demonstrate how Chempy can be used to dramatically improve elemental abundance predictions of hydrodynamical simulations by plugging tailored best-fit SSP parameters into a Milky Way analogue from Gutcke & Springel. Our code, including a comprehensive tutorial, is freely available and can additionally provide SSP enrichment tables for any combination of parameters and yield tables.

We use a sample of 83 high-resolution cosmological zoom-in simulations and a semi-analytic model (SAM) to study the stochasticity of galaxy formation in haloes ranging from dwarf to Milky Way masses. Our simulated galaxies reproduce the observed inefficiency of galaxy formation as expressed through the stellar, gas and baryonic Tully-Fisher relations. For HI velocities in the range (70 < V < 220 km/s), the scatter is just 0.08 to 0.14 dex, consistent with the observed intrinsic scatter at these scales. At low velocities (20 < V < 70 km/s), the simulated scatter is 0.2-0.25 dex, which could be tested with future observations. The scatter in the stellar mass versus dark halo velocity relation is constant for 30 < V < 180 km/s, and smaller (~0.17 dex) when using the maximum circular velocity of the dark matter only simulation, Vmax, compared to the virial velocity ($V_{200}$ or $V_{200}^{DMO}$). The scatter in stellar mass is correlated with halo concentration, and is minimized when using a circular velocity at a fixed fraction of the virial radius $\simeq 0.4 R_{200}$ or with $V_{\alpha}=V_{200}^{DMO} (V_{max}^{DMO} / V_{200}^{DMO})^\alpha$ with \alpha ~ 0.7, consistent with constraints from halo clustering. Using the SAM we show the correlation between halo formation time and concentration is essential in order to reproduce this result. This uniformity in galaxy formation efficiency we see in our hydrodynamical simulations and a semi-analytic model proves the simplicity and self-regulating nature of galaxy formation in a \Lambda Cold Dark Matter universe.

We use the NIHAO galaxy formation simulations to make predictions for the baryonic budget in present day galaxies ranging from dwarf to Milky Way masses. The sample is made of 88 independent high resolution cosmological zoom-in simulations. NIHAO galaxies reproduce key properties of observed galaxies, such as the stellar mass vs halo mass and cold gas vs stellar mass relations. Thus they make plausible predictions for the baryon budget. We present the mass fractions of stars, cold gas ($T<10^4$K), cool gas ($10^4 < T < 10^5$K), warm-hot gas ($10^5 < T < 5\times10^6$K), and hot gas (T$> 5\times10^6$K) inside the virial radius, $R_{200}$. Compared to the predicted baryon mass, using the dark halo mass and the universal baryon fraction, $f_{\rm b}\equiv \Omega_{\rm b}/\Omega_{\rm m}=0.15$, we find that all of our haloes are missing baryons. The missing mass has been relocated past 2 virial radii, and cool gas dominates the corona at low mass (M$_{200} < 3 \times 10^{11} M_\odot$) while the warm-hot gas dominates at high mass (M$_{200} > 3 \times 10^{11} M_\odot$). Haloes of mass $\sim 10^{10} M_\odot$ are missing $\sim 90\%$ of their baryons. More massive haloes ($\sim 10^{12} M_\odot$) retain a higher fraction of their baryons, with $\sim 30\%$ missing, consistent with recent observational estimates. Moreover, these more massive haloes reproduce the observed fraction of cold, warm-hot and hot gas. The fraction of cool gas we predict ($0.11\pm0.06$) is significantly lower than the observation from COS-HALOs (0.3-0.47), but agrees with the alternative analysis of Stern et al. 2016.

The decrease in star formation (SF) and the morphological change necessary to produce the $z=0$ elliptical galaxy population are commonly ascribed to a sudden quenching event, which is able to rid the central galaxy of its cold gas reservoir in a short time. Following this event, the galaxy is able to prevent further SF and stay quiescent via a maintenance mode. We test whether such a quenching event is truly necessary using a simple model of quiescence. In this model, hot gas (all gas above a temperature threshold) in a $\sim10^{12} M_{\odot}$ halo mass galaxy at redshift $z\sim3$ is prevented from cooling. The cool gas continues to form stars at a decreasing rate and the galaxy stellar mass, morphology, velocity dispersion and position on the color magnitude diagram (CMD) proceed to evolve. By $z=0$, the halo mass has grown to $10^{13} M_{\odot}$ and the galaxy has attained characteristics typical of an observed $z=0$ elliptical galaxy. Our model is run in the framework of a cosmological, smooth particle hydrodynamic code which includes SF, early stellar feedback, supernova feedback, metal cooling and metal diffusion. Additionally, we post-process our simulations with a radiative transfer code to create a mock CMD. In contrast to previous assumptions that a pure "fade away" model evolves too slowly to account for the sparsity of galaxies in the "green valley", we demonstrate crossing times of $<$1 Gyr. We conclude that no sudden quenching event is necessary to produce such rapid colour transitions.

We study the hot and cold circum-galactic medium (CGM) of 86 galaxies of the cosmological, hydrodynamical simulation suite NIHAO. NIHAO allows a study of how the $z=0$ CGM varies across 5 orders of magnitude of stellar mass using OVI and HI as proxies for hot and cold gas. The cool HI covering fraction and column density profiles match observations well, particularly in the inner CGM. OVI shows increasing column densities with mass, a trend seemingly echoed in the observations. As in multiple previous simulations, the OVI column densities in simulations are lower than observed and optically thick HI does not extend as far out as in observations. We take a look at the collisional ionisation fraction of OVI as a function of halo mass. We make observable predictions of the bipolarity of outflows and their effect on the general shape of the CGM. Bipolar outflows can be seen out to around 40 kpc in intermediate and low mass halos ($M_{\mathrm{halo}}<10^{11}M_{\mathrm{sun}}$), but outside that radius, the CGM is too well mixed to detect an elongated shape. Larger halos have extended gas discs beyond the stellar disc that dominate the shape of the inner CGM. The simulated CGM is remarkably spherical even in low mass simulations. The chemical enrichment of both halo and disc gas follow expected increasing trends as a function of halo mass that are well fit with power laws. These relations can be used in non-hydrodynamic models, such as semi-analytic models.

The Waltz Spectrograph is a fiber-fed high-resolution \'echelle spectrograph for the 72 cm Waltz Telescope at the Landessternwarte, Heidelberg. It uses a 31.6 lines/mm 63.5$^{\circ}$ blaze angle \'echelle grating in white-pupil configuration, providing a spectral resolving power of $R\sim$65,000 covering the spectral range between 450$-$800\,nm in one CCD exposure. A prism is used for cross-dispersion of \'echelle orders. The spectrum is focused by a commercial apochromat onto a 2k$\times$2k CCD detector with 13.5$\mu$m per pixel. An exposure meter will be used to obtain precise photon-weighted midpoints of observations, which will be used in the computation of the barycentric corrections of measured radial velocities. A stabilized, newly designed iodine cell is employed for measuring radial velocities with high precision. Our goal is to reach a radial velocity precision of better than 5 m/s, providing an instrument with sufficient precision and sensitivity for the discovery of giant exoplanets. Here we describe the design of the Waltz spectrograph and early on-sky results.

We use cosmological hydrodynamical zoom-in simulations with the SPH code gasoline of four haloes of mass M_{200} \sim 10^{13}\Msun to study the response of the dark matter to elliptical galaxy formation. Our simulations include metallicity dependent gas cooling, star formation, and feedback from massive stars and supernovae, but not active galactic nuclei (AGN). At z=2 the progenitor galaxies have stellar to halo mass ratios consistent with halo abundance matching, assuming a Salpeter initial mass function. However by z=0 the standard runs suffer from the well known overcooling problem, overpredicting the stellar masses by a factor of > 4. To mimic a suppressive halo quenching scenario, in our forced quenching (FQ) simulations, cooling and star formation are switched off at z=2. The resulting z=0 galaxies have stellar masses, sizes and circular velocities close to what is observed. Relative to the control simulations, the dark matter haloes in the FQ simulations have contracted, with central dark matter density slopes d\log\rho/d\log r \sim -1.5, showing that dry merging alone is unable to fully reverse the contraction that occurs at z>2. Simulations in the literature with AGN feedback however, have found expansion or no net change in the dark matter halo. Thus the response of the dark matter halo to galaxy formation may provide a new test to distinguish between ejective and suppressive quenching mechanisms.

In a Universe where AGN feedback regulates star formation in massive galaxies, a strong correlation between these two quantities is expected. If the gas causing star formation is also responsible for feeding the central black hole, then a positive correlation is expected. If powerful AGNs are responsible for the star formation quenching, then a negative correlation is expected. Observations so far have mainly found a mild correlation or no correlation at all (i.e. a flat relation between star formation rate (SFR) and AGN luminosity), raising questions about the whole paradigm of "AGN feedback". In this paper, we report the predictions of the GALFORM semi-analytical model, which has a very strong coupling between AGN activity and quenching of star formation. The predicted SFR-AGN luminosity correlation appears negative in the low AGN luminosity regime, where AGN feedback acts, but becomes strongly positive in the regime of the brightest AGN. Our predictions reproduce reasonably well recent observations by Rosario et al., yet there is some discrepancy in the normalisation of the correlation at low luminosities and high redshifts. Though this regime could be strongly influenced by observational biases, we argue that the disagreement could be ascribed to the fact that GALFORM neglects AGN variability effects. Interestingly, the galaxies that dominate the regime where the observations imply a weak correlation are massive early-type galaxies that are subject to AGN feedback. Nevertheless, these galaxies retain high enough molecular hydrogen contents to maintain relatively high star formation rates and strong infrared emission.