Search

It has long been discovered that environment plays a major role in galaxy evolution, with galaxies in dense regions of the cosmos more likely to have their star formation suppressed or \textquote{quenched} relative to their more isolated counterparts. Despite this remarkable discovery, our current understanding of the physical processes responsible for the suppression (or \textquote{quenching}) of star formation remains woefully incomplete. This is especially true for galaxies that are members of massive galaxy groups and clusters – i.e. \textquote{satellite} galaxies – given that our best cosmological models struggle to reproduce key observations, such as the satellite quenched fraction. This dissertation aims to advance our knowledge in the areas of environmental quenching that are poorly understood, focusing on the quenching of low-mass dwarf satellite galaxies and satellite populations at intermediate and high redshifts. Observational studies of low-mass satellite quenching have been limited by their inability to robustly characterize the local environment and star-formation activity of faint systems. This work overcomes these limitations by combining supervised machine learning and statistical background subtraction techniques to constrain the satellite quenched fraction of group populations ($\mhalo = 10^{13-14}~\msun$) down to a satellite stellar mass limit of ${\sim}10^7~\msun$. This approach successfully reproduces the established quenched fraction trends at high-masses while finding that ram-pressure stripping -- the rapid removal of cold gas from the interstellar medium of satellites as they move through the hot and dense intra-group or intracluster medium permeating the host halo -- is the likely dominant quenching mechanism responsible for shutting down star formation in the low-mass regime. The importance of investigating intermediate and high redshift is that the vast majority of environmental quenching studies are limited to the very local Universe ($z \lesssim 0.1$). This work addresses this shortcoming by combining multi-wavelength spectroscopic observations of massive clusters ($\mhalo = 10^{14-15}~\msun$) at $z\gtrsim1$ with the infall histories of analogs from $N$-body simulations to model satellite quenching at this poorly studied epoch. The quenching timescales derived from this study indicate that \textquote{starvation} -- the slow depletion of cold gas in the absence of cosmological accretion after a satellite becomes a member of a more massive host -- is the dominant driver of environmental quenching at this epoch. However, this study lends support to the scenario where ram-pressure stripping potentially acts as a secondary quenching mechanism in massive clusters at this epoch. Overall, this work provides valuable insights into the quenching timescales and underlying mechanisms of satellite galaxies in massive clusters at $z\gtrsim1$, while introducing an impactful methodology that combines supervised machine learning and statistical background subtraction to study the quenching of low-mass satellites beyond the Local Group.

As the closest $L^{\ast}$ galaxy to our own Milky Way, the Andromeda galaxy (M31) is an ideal laboratory for studies of galaxy evolution. The AMIGA project has recently provided observations of the cool ($T\sim10^4$ K) phase of the circumgalactic medium (CGM) of M31, using HST/COS absorption spectra along $\sim40$ background QSO sightlines, located up to and beyond the galaxy virial radius. Based on these data, and by the means of semi-analytic models and Bayesian inference, we provide here a physical description of the origin and dynamics of the cool CGM of M31. We investigate two competing scenarios, in which (i) the cool gas is mostly produced by supernova(SN)-driven galactic outflows or (ii) it mostly originates from infall of gas from the intergalactic medium. In both cases, we take into account the effect of gravity and hydrodynamical interactions with a hot corona, which has a cosmologically motivated angular momentum. We compare the outputs of our models to the observed covering factor, silicon column density and velocity distribution of the AMIGA absorbers. We find that, to explain the observations, the outflow scenario requires an unphysically large (> 100\%) efficiency for SN feedback. Our infall models, on the other hand, can consistently account for the AMIGA observations and the predicted accretion rate, angular momentum and metallicity are consistent with a cosmological infall from the intergalactic medium., Accepted for publication in MNRAS

In this paper, we present data from the Blind Ultra-Deep H i Environmental Survey (BUDHiES), which is a blind 21-cm H i spectral line imaging survey undertaken with the Westerbork Synthesis Radio Telescope. Two volumes were surveyed, each with a single pointing and covering a redshift range of 0.164 < z < 0.224. Within these two volumes, this survey targeted the clusters Abell 963 and Abell 2192, which are dynamically different and offer unique environments to study the process of galaxy evolution within clusters. With an integration time of 117 × 12 h on Abell 963 and 72 × 12 h on Abell 2192, a total of 166 galaxies were detected and imaged in H i. While the clusters themselves occupy only 4 per cent of the 73 400 Mpc3 surveyed by BUDHiES, most of the volume consists of large-scale structures in which the clusters are embedded, including foreground and background overdensities and voids. We present the data processing and source detection techniques and counterpart identification based on a wide-field optical imaging survey using the Isaac Newton Telescope and deep ultraviolet (UV) Galaxy Evolution Explorer (GALEX) imaging. Finally, we present H i and optical catalogues of the detected sources as well as atlases of their global H i properties, which include integrated column density maps, position–velocity diagrams, global H i profiles, and optical and UV images of the H i sources.

A crucial question in galaxy formation is what role new accretion has in star formation. Theoretical models have predicted a wide range of correlation strengths between halo accretion and galaxy star formation. Previously, we presented a technique to observationally constrain this correlation strength for isolated Milky Way-mass galaxies at $z\sim 0.12$, based on the correlation between halo accretion and the density profile of neighbouring galaxies. By applying this technique to both observational data from the Sloan Digital Sky Survey and simulation data from the UniverseMachine, where we can test different correlation strengths, we ruled out positive correlations between dark matter accretion and recent star formation activity. In this work, we expand our analysis by (1) applying our technique separately to red and blue neighbouring galaxies, which trace different infall populations, (2) correlating dark matter accretion rates with $D_{n}4000$ measurements as a longer-term quiescence indicator than instantaneous star-formation rates, and (3) analyzing higher-mass isolated central galaxies with $10^{11.0} < M_*/M_\odot < 10^{11.5}$ out to $z\sim 0.18$. In all cases, our results are consistent with non-positive correlation strengths with $\gtrsim 85$ per cent confidence, suggesting that processes such as gas recycling dominate star formation in massive $z=0$ galaxies., 16 pages; accepted for publication in MNRAS

Dual-Active Galactic Nuclei (AGN) are a natural consequence of the hierarchical structure formation scenario, and can provide an important test of various models for black hole growth. However, due to their rarity and difficulty to find at high redshift, very few confirmed dual-AGN are known at the epoch where galaxy formation peaks. Here we report the discovery of a gravitationally lensed dual-AGN system at redshift 2.37 comprising two optical/IR quasars separated by 6.5+/-0.6 kpc, and a third compact (R_eff = 0.45+/-0.02 kpc) red galaxy that is offset from one of the quasars by 1.7+/-0.1 kpc. From Very Large Array imaging at 3 GHz, we detect 600 and 340 pc-scale radio emission that is associated with both quasars. The 1.4 GHz luminosity densities of the radio sources are about 10^24.35 W / Hz, which is consistent with weak jets. However, the low brightness temperature of the emission is also consistent with star-formation at the level of 850 to 1150 M_sun / yr. Although this supports the scenario where the AGN and/or star-formation is being triggered through an ongoing triple-merger, a post-merger scenario where two black holes are recoiling is also possible, given that neither has a detected host galaxy., Comment: 5 pages, 3 figures, accepted for publication in MNRAS Letters

We now have to set the scene for our attack on the major problems of astrophysical cosmology — how do we account for the origin and evolution of the galaxies and the large-scale structure of our Universe? This is one of the most exciting areas of modern cosmology. The first thing we have to do is to work out the thermal history of the matter and radiation content of the standard Big Bang picture. In this and the succeeding Chap. 10, we will develop many concepts which are essential for understanding the problems of galaxy formation. Let us first consider the dynamics of radiation-dominated universes.

To make quantitative comparisons between the theories of galaxy formation described in Chap. 13 and the observed distribution of galaxies, we need to make assumptions about the spectrum of the fluctuations from which galaxies and larger scale structures formed. The natural way of describing the distribution of galaxies on the large scale is in terms of correlation functions and their associated power spectra. The objective is to relate these properties of the distribution of galaxies at the present day to the spectrum of initial fluctuations which must have been present in the very early Universe. Let us first review the properties of the correlation function of galaxies at the present day, extending our introduction of Sect. 2.2.

Let us recall the reasons for taking dark matter, and in particular non-baryonic dark matter, very seriously in the context of galaxy formation.

The demographics of dark matter substructure depend sensitively on the nature of dark matter. Optimally leveraging this probe requires accurate theoretical predictions regarding the abundance of subhaloes. These predictions are hampered by artificial disruption in numerical simulations, by large halo-to-halo variance, and by the fact that the results depend on the baryonic physics of galaxy formation. In particular, numerical simulations have shown that the formation of a central disc can drastically reduce the abundance of substructure compared to a dark matter-only simulation, which has been attributed to enhanced destruction of substructure due to disc shocking. We examine the impact of discs on substructure using the semi-analytical subhalo model SatGen, which accurately models the tidal evolution of substructure free of the numerical disruption that still hampers $N$-body simulations. Using a sample of 10,000 merger trees of Milky Way-like haloes, we study the demographics of subhaloes that are evolved under a range of composite halo-disc potentials with unprecedented statistical power. We find that the overall subhalo abundance is relatively insensitive to properties of the disc aside from its total mass. For a disc that contains $5\%$ of $M_\mathrm{vir}$, the mean subhalo abundance within $r_\mathrm{vir}$ is suppressed by ${\lesssim}10\%$ relative to the no-disc case, a difference that is dwarfed by halo-to-halo variance. For the same disc mass, the abundance of subhaloes within 50 kpc is reduced by ${\sim}30\%$. We argue that the disc mainly drives excess mass loss for subhaloes with small pericentric radii and that the impact of disc shocking is negligible., 13 pages, 8 figures. Accepted to MNRAS

The ratio of the mass of interstellar dust to the total mass of metals (the dust-to-metals/DTM ratio) tends to increase with metallicity. This can be explained by the increasing efficiency of grain growth in the interstellar medium (ISM) at higher metallicities, with a corollary being that the low DTM ratios seen at low metallicities are due to inefficient stellar dust production. This interpretation assumes that the efficiency of dust destruction in the ISM is constant, whereas it might be expected to increase at low metallicity; the decreased cooling efficiency of low-metallicity gas should result in more post-shock dust destruction via thermal sputtering. We show that incorporating a sufficiently strong metallicity dependence into models of galaxy evolution removes the need for low stellar dust yields. The contribution of stellar sources to the overall dust budget may be significantly underestimated, and that of grain growth overestimated, by models assuming a constant destruction efficiency., Comment: 5 pages, 2 figures. MNRAS accepted

We study the population of active galaxies in void environment in the Sloan Digital Sky Survey. We use optical spectroscopic information to analyse characteristics of the emission lines of galaxies, accomplished by WHAN and BPT diagrams. Also, we study Wide-field Infrared Survey Explorer(WISE) mid-IR colours to assess active galactic nucleus (AGN) activity. We investigate these different AGN classification schemes, both optical and mid-IR, and their dependence on the spatial location with respect to the void centres. To this end, we define three regions: void, the spherical region defined by voidcentric distance relative to void radius (distance/rvoid) smaller than 0.8, comprising overdensities lesser than −0.9, an intermediate/transition shell region (namely, void-wall) 0.8 < distance/rvoid < 1.2, and a region sufficiently distant from voids, the field: distance/rvoid > 2. We find statistical evidence for a larger fraction of AGN and star-forming galaxies in the void region, regardless of the classification scheme addressed (either BPT, WHAN, or WISE). Moreover, we obtain a significantly stronger nuclear activity in voids compared to the field. We find an unusually large fraction of the most massive black holes undergoing strong accretion when their host galaxies reside in voids. Our results suggest a strong influence of the void environment on AGN mechanisms associated with galaxy evolution.

In this paper, we investigate the properties of binary neutron stars (BNSs) and their mergers by combining population synthesis models for binary stellar evolution (BSE) with cosmological galaxy formation and evolution models. We obtain constraints on BSE model parameters by using the observed Galactic BNSs and local BNS merger rate density ($R_0$) inferred from Gravitational Wave (GW) observations, and consequently estimate the host galaxy distributions of BNS mergers. We find that the Galactic BNS observations imply efficient energy depletion in the common envelope (CE) phase, a bimodal kick velocity distribution, and low mass ejection during the secondary supernova explosion. However, the inferred $R_0$ does not necessarily require an extremely high CE ejection efficiency and low kick velocities, different from the previous claims, mainly because the latest inferred $R_0$ is narrowed to a lower value ($320_{-240}^{+490}\,{\rm Gpc^{-3}\,yr^{-1}}$). The BNS merger rate density resulting from the preferred model can be described by $R(z)\sim R_0(1+z)^{\zeta}$ at low redshift ($z\lesssim0.5$), with $R_0\sim316$-$784\,{\rm Gpc^{-3}\,yr^{-1}}$ and $\zeta\sim1.34$-$2.03$, respectively. Our results also show that $R_{0}$ and $\zeta$ depend on settings of BSE model parameters, and thus accurate estimates of these parameters by future GW detections will put strong constraints on BSE models. We further estimate that the fractions of BNS mergers hosted in spiral and elliptical galaxies at $z\sim0$ are $\sim81$%-$84$% and $\sim16$%-$19$%, respectively. The BNS merger rate per galaxy can be well determined by the host galaxy stellar mass, star formation rate, and metallicity, which provides a guidance in search for most probable candidates of BNS host galaxies., Comment: 32 pages, 17 figures, accepted for publication in MNRAS

The locations of Ly-$\alpha$ emitting galaxies (LAEs) at the end of the Epoch of Reionisation (EoR) are expected to correlate with regions of ionised hydrogen, traced by the redshifted 21~cm hyperfine line. Mapping the neutral hydrogen around regions with detected and localised LAEs offers an avenue to constrain the brightness temperature of the Universe within the EoR by providing an expectation for the spatial distribution of the gas, thereby providing prior information unavailable to power spectrum measurements. We use a test set of 12 hours of observations from the Murchison Widefield Array (MWA) in extended array configuration, to constrain the neutral hydrogen signature of 58 LAEs, detected with the Subaru Hypersuprime Cam in the \textit{Silverrush} survey, centred on $z$=6.58. We assume that detectable emitters reside in the centre of ionised HII bubbles during the end of reionization, and predict the redshifted neutral hydrogen signal corresponding to the remaining neutral regions using a set of different ionised bubble radii. A prewhitening matched filter detector is introduced to assess detectability. We demonstrate the ability to detect, or place limits upon, the amplitude of brightness temperature fluctuations, and the characteristic HII bubble size. With our limited data, we constrain the brightness temperature of neutral hydrogen to $\Delta{\rm T}_B, Comment: 9 pages, 6 figures, accepted for publication in MNRAS

The population of satellite galaxies in a host galaxy is a combination of the cumulative accretion of subhaloes and their associated star formation efficiencies, therefore, the luminosity distribution of satellites provides valuable information of both dark matter properties and star formation physics. Recently, the luminosity function of satellites in nearby Milky Way-mass galaxies has been well measured to satellites as faint as Leo I with $M_{V} \sim -8$. In addition to the finding of the diversity in the satellite luminosity functions, it has been noticed that there is a big gap among the magnitude of satellites in some host galaxies, such as M101, where the gap is around 5 in magnitude, noticeably larger than the prediction from the halo abundance matching method. The reason of this gap is still unknown. In this paper, we use a semi-analytical model of galaxy formation, combined with high-resolution N-body simulation, to investigate the probability and origin of such big gap in M101-alike galaxies. We found that, although M101 analogues are very rare with probability of \sim 0.1%-0.2% in the local universe, their formation is a natural outcome of the CDM model. The gap in magnitude is mainly due to the mass of the accreted subhaloes, not from the stochastic star formation in them. We also found that the gap is correlated with the total satellite mass and host halo mass. By tracing the formation history of M101 type galaxies, we find that they likely formed after $z \sim 1$ due to the newly accreted bright satellites. The gap is not in a stable state, and it will disappear in ~7 Gyr due to mergers of bright satellites with the central galaxy., Comment: 8 pages, 11 figures, accepted by MNRAS

We study the evolution of nitrogen in the Galactic halo, thick disc, thin disc and bulge by comparing detailed chemical evolution models with recent observations. The models used in this work have already been constrained to explain the abundance patterns of $\alpha$-elements and the metallicity distribution functions of halo, disc and bulge stars; here, we adopt them to investigate the origin and evolution of N in the different Galactic components. First, we consider different sets of yields and study the importance of the various channels proposed for N production. Secondly, we apply the reference models to study the evolution of both the Galactic discs and bulge. We conclude that: i) primary N produced by rotating massive stars is required to reproduce the plateau in log(N/O) and [N/Fe] ratios at low metallicity, as well as the secondary and primary production from low- and intermediate-mass stars to reproduce the data of the thin disc. ii) The parallel model can provide a good explanation of the evolution of N abundance in the thick and thin discs; we confirm that the thick disc has evolved much faster than the thin disc, in agreement with the results from the abundance patterns of other chemical elements. iii) Finally, we present new model predictions for N evolution in the Galactic bulge, and we show that the observations in bulge stars can be explained if massive stars rotate fast during the earliest phases of Galactic evolution, in agreement with findings from the abundance pattern of carbon., Comment: 10 pages, 4 figures. Accepted for publication in MNRAS

Upcoming galaxy surveys and 21-cm experiments targeting high redshifts $z\gtrsim 6$ are highly complementary probes of galaxy formation and reionization. However, in order to expedite the large volume simulations relevant for 21-cm observations, many models of galaxies within reionization codes are entirely subgrid and/or rely on halo abundances only. In this work, we explore the extent to which resolving and modeling individual galaxy formation histories affects predictions both for the galaxy populations detectable by upcoming surveys and the signatures of reionization accessible to upcoming 21-cm experiments. We find that a common approach, in which galaxy luminosity is assumed to be a function of halo mass only, is biased with respect to models in which galaxy properties are evolved through time via semi-analytic modeling and thus reflective of the diversity of assembly histories that naturally arise in $N$-body simulations. The diversity of galaxy formation histories also results in scenarios in which the brightest galaxies do \textit{not} always reside in the centers of large ionized regions, as there are often relatively low-mass halos undergoing dramatic, but short-term, growth. This has clear implications for attempts to detect or validate the 21-cm background via cross correlation. Finally, we show that a hybrid approach -- in which only halos hosting galaxies bright enough to be detected in surveys are modeled in detail, with the rest modeled as an unresolved field of halos with abundance related to large-scale overdensity -- is a viable way to generate large-volume `simulations` well-suited to wide-area surveys and current-generation 21-cm experiments targeting relatively large $k \lesssim 1 \ h \ \text{Mpc}^{-1}$ scales., Comment: 17 pages, 14 figures, accepted for publication in MNRAS

Understanding the impact of halo properties beyond halo mass on the clustering of galaxies (namely galaxy assembly bias) remains a challenge for contemporary models of galaxy clustering. We explore the use of machine learning to predict the halo occupations and recover galaxy clustering and assembly bias in a semi-analytic galaxy formation model. For stellar-mass selected samples, we train a Random Forest algorithm on the number of central and satellite galaxies in each dark matter halo. With the predicted occupations, we create mock galaxy catalogues and measure the clustering and assembly bias. Using a range of halo and environment properties, we find that the machine learning predictions of the occupancy variations with secondary properties, galaxy clustering and assembly bias are all in excellent agreement with those of our target galaxy formation model. Internal halo properties are most important for the central galaxies prediction, while environment plays a critical role for the satellites. Our machine learning models are all provided in a usable format. We demonstrate that machine learning is a powerful tool for modelling the galaxy-halo connection, and can be used to create realistic mock galaxy catalogues which accurately recover the expected occupancy variations, galaxy clustering and galaxy assembly bias, imperative for cosmological analyses of upcoming surveys., 21 pages, 17 figures, submitted to MNRAS

Next generation telescopes such as the James Webb Space Telescope (JWST) and the Nancy Grace Roman Space Telescope (NGRST) will enable us to study the first billion years of our Universe in unprecedented detail. In this work we use the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework, that couples galaxy formation and reionization (for a wide range of reionization feedback models), to estimate the cosmic variance expected in the UV Luminosity Function (UV LF) and the Stellar Mass Function (SMF) in JWST surveys. We find that different reionization scenarios play a minor role in the cosmic variance. Most of the cosmic variance is completely driven by the underlying density field and increases above $100\%$ for ${\rm M}_{\rm UV} \sim -17.5 (-20)$ at $z =12 (6)$ for the JADES-deep survey (the deep JWST Advanced Extragalactic Survey with an area of 46 arcmin$^2$); the cosmic variance decreases with an increasing survey area roughly independently of redshift. We find that the faint-end (${\rm M}_{\rm UV} > -17$) slope of the Lyman Break Galaxies (LBGs) UV LF becomes increasingly shallower with increasing reionization feedback and show how JWST observations will be able to distinguish between different models of reionization feedback at $z>9$, even accounting for cosmic variance. We also show the environments (in terms of density and ionization fields)of Lyman Break Galaxies during the EoR, finding that the underlying over-density and ionization fraction scale positively with the UV luminosity. Finally, we also provide a public software tool to allow interested readers to compute cosmic variance for different redshifts and survey areas., 16 pages, 13 figures. Accepted for publication in MNRAS

In this work, we use the {\sc astraeus} (seminumerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework which couples galaxy formation and reionization in the first billion years. Exploring a number of models for reionization feedback and the escape fraction of ionizing radiation from the galactic environment ($f_\mathrm{esc}$), we quantify how the contribution of star-forming galaxies {(with halo masses $M_h>10^{8.2}$M$_\odot$)} to reionization depends on the radiative feedback model, $f_\mathrm{esc}$, and the environmental over-density. Our key findings are: (i) for constant $f_\mathrm{esc}$ models, intermediate-mass galaxies (with halo masses of $M_h\simeq10^{9-11}$M$_\odot$ and absolute UV magnitudes of $M_{UV} \sim -15$ to $-20$) in intermediate-density regions drive reionization; (ii) scenarios where $f_\mathrm{esc}$ increases with decreasing halo mass shift the galaxy population driving reionization to lower-mass galaxies ($M_h\lesssim10^{9.5}$M$_\odot$) with lower luminosities ($M_{UV} \gtrsim-16$) and over-densities; (iii) reionization imprints its topology on the ionizing emissivity of low-mass galaxies ($M_h\lesssim10^{9}$M$_\odot$) through radiative feedback. Low-mass galaxies experience a stronger suppression of star formation by radiative feedback and show lower ionizing emissivities in over-dense regions; (iv) a change in $f_\mathrm{esc}$ with galaxy properties has the largest impact on the sources of reionization and their detectability, with the radiative feedback strength and environmental over-density playing a sub-dominant role; (v) JWST-surveys (with a limiting magnitude of $M_{UV} = -16$) will be able to detect the galaxies providing $\sim 60-70\%$ ($\sim 10\%$) of reionization photons at $z=7$ for constant $f_\mathrm{esc}$ models (scenarios where $f_\mathrm{esc}$ increases with decreasing halo mass)., 14 pages, 13 figures, accepted for publication in MNRAS

The diffuse extragalactic background light (EBL) is formed by ultraviolet (UV), optical, and infrared (IR) photons mainly produced by star formation processes over the history of the Universe, and contains essential information about galaxy evolution and cosmology. Here, we present a new determination of the evolving EBL spectral energy distribution using a novel approach purely based on galaxy data aiming to reduce current uncertainties on the higher redshifts and IR intensities. Our calculations use multiwavelength observations from the UV to the far-IR of a sample of approximately 150,000 galaxies detected up to $z\sim 6$ in the five fields of the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) from the Hubble Space Telescope. This is one of the most comprehensive and deepest multi-wavelength galaxy datasets ever obtained. These unprecedented resources allow us to derive the overall EBL evolution up to $z\sim 6$ and its uncertainties. Our results agree with cosmic observables estimated from galaxy surveys and $\gamma$-ray attenuation such as monochromatic luminosity densities, including those in the far-IR, and star formation rate densities, also at the highest redshits. Optical depths from our EBL approximation, which will be robust at high redshifts and for $\gamma$ rays up to tens of TeV, will be reported in a companion paper., Comment: 17 pages, 13 figures, 4 tables; accepted for publication in MNRAS on August 6, 2021. Online material available at https://www.ucm.es/blazars/ebl

We develop a model to establish the interconnection between galaxies and their dark matter halos. We use Principal Component Analysis (PCA) to reduce the dimensionality of both the mass assembly histories of halos/subhalos and the star formation histories of galaxies, and Gradient Boosted Decision Trees (GBDT) to transform halo/subhalo properties into galaxy properties. We use two sets of hydrodynamic simulations to motivate our model architecture and to train the transformation. We then apply the two sets of trained models to dark matter only (DMO) simulations to show that the transformation is reliable and statistically accurate. The model trained by a high-resolution hydrodynamic simulation, or by a set of such simulations implementing the same physics of galaxy formation, can thus be applied to large DMO simulations to make "mock" copies of the hydrodynamic simulation. The model is both flexible and interpretable, which paves the way for future applications in which we will constrain the model using observations at different redshifts simultaneously and explore how galaxies form and evolve in dark matter halos empirically., Comment: 21 pages, 16 figures, 2 tables, MNRAS accepted

Sparse regression algorithms have been proposed as the appropriate framework to model the governing equations of a system from data, without needing prior knowledge of the underlying physics. In this work, we use sparse regression to build an accurate and explainable model of the stellar mass of central galaxies given properties of their host dark matter (DM) halo. Our data set comprises 9,521 central galaxies from the EAGLE hydrodynamic simulation. By matching the host halos to a DM-only simulation, we collect the halo mass and specific angular momentum at present time and for their main progenitors in 10 redshift bins from $z=0$ to $z=4$. The principal component of our governing equation is a third-order polynomial of the host halo mass, which models the stellar-mass halo-mass relation. The scatter about this relation is driven by the halo mass evolution and is captured by second and third-order correlations of the halo mass evolution with the present halo mass. An advantage of sparse regression approaches is that unnecessary terms are removed. Although we include information on halo specific angular momentum, these parameters are discarded by our methodology. This suggests that halo angular momentum has little connection to galaxy formation efficiency. Our model has a root mean square error (RMSE) of $0.167 \log_{10}(M^*/M_\odot)$, and accurately reproduces both the stellar mass function and central galaxy correlation function of EAGLE. The methodology appears to be an encouraging approach for populating the halos of DM-only simulations with galaxies, and we discuss the next steps that are required., Comment: 20 pages, 15 figures

In order to answer the open questions of modern cosmology and galaxy evolution theory, robust algorithms for calculating photometric redshifts (photo-z) for very large samples of galaxies are needed. Correct estimation of the various photo-z algorithms' performance requires attention to both the performance metrics and the data used for the estimation. In this work, we use the supervised machine learning algorithm MLPQNA to calculate photometric redshifts for the galaxies in the COSMOS2015 catalogue and the unsupervised Self-Organizing Maps (SOM) to determine the reliability of the resulting estimates. We find that for spec-z, 26 pages, 15 figures, accepted for publication in MNRAS

Galaxy clusters provide an excellent probe in various research fields in astrophysics and cosmology. However, the number of galaxy clusters detected so far in the $AKARI$ North Ecliptic Pole (NEP) field is limited. In this work, we provide galaxy cluster candidates in the $AKARI$ NEP field with the minimum requisites based only on coordinates and photometric redshift (photo-$z$) of galaxies. We used galaxies detected in 5 optical bands ($g$, $r$, $i$, $z$, and $Y$) by the Subaru Hyper Suprime-Cam (HSC), assisted with $u$-band from Canada-France-Hawaii Telescope (CFHT) MegaPrime/MegaCam, and IRAC1 and IRAC2 bands from the $Spitzer$ space telescope for photo-$z$ estimation. We calculated the local density around every galaxy using the 10$^{th}$-nearest neighbourhood. Cluster candidates were determined by applying the friends-of-friends algorithm to over-densities. 88 cluster candidates containing 4390 member galaxies below redshift 1.1 in 5.4 deg$^2$ have been detected. The reliability of our method was examined through false detection tests, redshift uncertainty tests, and applications on the COSMOS data, giving false detection rates of 0.01 to 0.05 and recovery rate of 0.9 at high richness. 3 X-ray clusters previously observed by $ROSAT$ and $Chandra$ were recovered. The cluster galaxies show higher stellar mass and lower star formation rate (SFR) compared to the field galaxies in two-sample Z-tests. These cluster candidates are useful for environmental studies of galaxy evolution and future astronomical surveys in the NEP, where $AKARI$ has performed unique 9-band mid-infrared photometry for tens of thousands of galaxies., Comment: 19 pages, 16 figures, has been accepted for publication in MNRAS

It has long been argued that the radial distribution of globular clusters (GCs) in the Fornax dwarf galaxy requires its dark matter halo to have a core of size $\sim 1$ kpc. We revisit this argument by investigating analogues of Fornax formed in E-MOSAICS, a cosmological hydrodynamical simulation that self-consistently follows the formation and evolution of GCs in the EAGLE galaxy formation model. In EAGLE, Fornax-mass haloes are cuspy and well described by the Navarro-Frenk-White profile. We post-process the E-MOSAICS to account for GC orbital decay by dynamical friction, which is not included in the original model. Dynamical friction causes 33 per cent of GCs with masses $M_{\rm GC}\geq4\times10^4 {~\rm M_\odot}$ to sink to the centre of their host where they are tidally disrupted. Fornax has a total of five GCs, an exceptionally large number compared to other galaxies of similar stellar mass. In the simulations, we find that only 3 per cent of the Fornax analogues have five or more GCs, while 30 per cent have only one and 35 per cent have none. We find that GC systems in satellites are more centrally concentrated than in field dwarfs, and that those formed in situ (45 per cent) are more concentrated than those that were accreted. The survival probability of a GC increases rapidly with the radial distance at which it formed ($r_{\rm init}$): it is 37 per cent for GCs with $r_{\rm init} \leq 1$ kpc and 92 per cent for GCs with $r_{\rm init} \geq 1$ kpc. The present-day radial distribution of GCs in E-MOSAICS turns out to be indistinguishable from that in Fornax, demonstrating that, contrary to claims in the literature, the presence of five GCs in the central kiloparsec of Fornax does not exclude a cuspy DM halo., 15 pages, 13 figures. submitted to MNRAS, comments welcome!

This document describes the general process of setting up, running, and analysing disc galaxy simulations using the freely available program Phantom of RAMSES (PoR). This implements Milgromian Dynamics (MOND) with a patch to the RAMSES grid-based $N$-body and hydrodynamical code that uses adaptive mesh refinement. We discuss the procedure of setting up isolated and interacting disc galaxy initial conditions for PoR, running the simulations, and analysing the results. This manual also concisely documents all previously developed MOND simulation codes and the results obtained with them., 10 pages. The article has been accepted for publication in the Canadian Journal of Physics. Files mentioned in the article can be found here, https://bitbucket.org/SrikanthTN/bonnpor/src/master/

We study the distribution of accretion rates onto supermassive black holes (SMBH) in AGNs of the local Universe ($z0.01$) in the local Universe also depends weakly on $M_\ast$ and is 0.2-1%. These results confirm the trends found in previous studies for the earlier Universe and refine the parameters of the dependence $f(\lambda_{Edd}|M_\ast)$ at $z, Comment: 16 pages, 12 figures. Published in Astronomy Letters

The $([\alpha/{\rm Fe}],[{\rm Fe/H}])$ distribution of Milky Way stars shows at least two distinct sequences, which have traditionally been associated with the thin and thick disc components. The abundance distribution varies systematically with location $R$ and $|z|$ across the Galaxy. Using an analytical chemodynamical model that includes the effects of radial migration and kinematic heating, we show that it is possible to reproduce the observed abundance distribution at different locations. Unlike some earlier models, our scheme does not require a distinct thick disc component emerging from a separate evolutionary path. The proposed model has a continuous star formation history and a continuous age velocity dispersion relation. Moreover, [$\alpha$/Fe] is constant for stellar ages less than 8 Gyr, but increases sharply for older stars over a time scale of 1.5 Gyr. The gap between the two sequences is due to this sharp transition. We show that the high-[$\alpha$/Fe] sequence at the low metallicity end is simply a pile-up of old stars. At the high metallicity end, we find a sequence of stars having different ages with a similar birth radius that originated in the inner disc. Our model successfully explains the uniformity of the locus of the high-[$\alpha$/Fe] sequence across different locations. The low-[$\alpha$/Fe] sequence contains stars with different birth radii that owes its existence to radial migration. For the low-[$\alpha$/Fe] sequence, angular momentum is anti-correlated with [Fe/H], and the model can reproduce this trend at different Galactic locations. If radial migration is not included, the model fails to generate the double sequence and instead shows only a single sequence. Our simple scheme has major advantages over earlier chemodynamical models, as we show., Comment: 22 pages, 12 figures

Post-starburst galaxies (PSBs) belong to a short-lived transition population between star-forming (SF) and quiescent galaxies. Deciphering their heavily discussed evolutionary pathways is paramount to understanding galaxy evolution. We aim to determine the dominant mechanisms governing PSB evolution in both the field and in galaxy clusters. Using the cosmological hydrodynamical simulation suite Magneticum Pathfinder, we identify 647 PSBs with $z \sim 0$ stellar mass $M_* \geq 5 \cdot 10^{10} \, \mathrm{M_{\odot}}$. We track their galactic evolution, merger history, and black hole activity over a time-span of 3.6Gyr. Additionally, we study cluster PSBs identified at different redshifts and cluster masses. Independent of environment and redshift, we find that PSBs, like SF galaxies, have frequent mergers. At z=0, 89% of PSBs have experienced mergers and 65% had at least one major merger within the last 2.5Gyr, leading to strong star formation episodes. In fact, 23% of z=0 PSBs were rejuvenated during their starburst. Following the mergers, field PSBs are generally shutdown via a strong increase in AGN feedback (power output $P_{AGN,PSB} \geq 10^{56}\,$erg/Myr). We find agreement with observations for both stellar mass functions and z = 0.9 line-of-sight phase space distributions of PSBs in galaxy clusters. Finally, we find that $z \lesssim 0.5$ cluster PSBs are predominantly infalling, especially in high mass clusters and show no signs of enhanced AGN activity. Thus, we conclude that the majority of cluster PSBs are shutdown via an environmental quenching mechanism such as ram-pressure stripping, while field PSBs are mainly quenched by AGN feedback., 28 pages, 15 figures, accepted by MNRAS

The change of the total mass density slope, $\gamma$, of early-type galaxies through cosmic time is a probe of evolutionary pathways. Hydrodynamical cosmological simulations show that at high redshifts density profiles of early-type galaxies were on average steep ($\gamma \sim -3$). As redshift approaches zero, gas-poor mergers progressively cause the total mass density slope to approach the `isothermal' slope of $\gamma \sim -2$. Simulations therefore predict steep density slopes at high redshifts, with little to no evolution in density slopes below $z \sim 1$. Gravitational lensing results in the same redshift range find the opposite, namely a significant trend of shallow density slopes at high redshifts, becoming steeper as redshift approaches zero. Gravitational lensing results indicate a different evolutionary mechanism for early-type galaxies than dry merging, such as continued gas accretion or off-axis mergers. At redshift zero, isothermal solutions are obtained by both simulations and dynamical modelling. This work applies the Jeans dynamical modelling technique to observations of galaxies at intermediate redshifts ($0.29 < z < 0.55 $) in order to derive density slopes to address the tension between observations and simulations. We combine two-dimensional kinematic fields from MUSE data with Hubble Space Telescope photometry. The density slopes of 90 early-type galaxies from the Frontier Fields project are presented. The total sample has a median of ${\gamma = -2.11 \pm 0.03}$ (standard error), in agreement with dynamical modelling studies at redshift zero. The lack of evolution in total density slopes in the past 4-6 Gyrs supports a dry merging model for early-type galaxy evolution., Comment: 27 pages, 10 figures, accepted for publication in MNRAS

In the present work, the galaxy luminosity function (LF) has been studied for a sample of seven clusters in the redshift range (0.0 ≲ z ≲ 0.1), within Abell radius (1.5 h−1 Mpc), in the five SDSS passbands ugriz. In each case, the absolute magnitude distribution is found and then fitted with a Schechter function. The fitting is done, using the χ2 – minimization method to find the best values of Schechter parameters Ф* (normalization constant), M* (characteristic absolute magnitude), and α (faint-end slope). No remarkable changes are found in the values of M* and α, for any cluster, in any passband. Furthermore, the LF does not seem to vary with such cluster parameters as richness, velocity dispersion, and Bautz–Morgan morphology. Finally, it is found that M* becomes brighter toward redder bands, whereas almost no variation is seen in the value of α with passband, being around (−1.00).

Atomic hydrogen is a fundamental constituent of galaxies, and is the primary fuel for star formation. Understanding how the atomic hydrogen content of galaxies evolves with time is thus critical to understanding galaxy evolution. In this article, we discuss how one measures the atomic hydrogen content of distant galaxies, and what these measurements tell us about the evolution of our Universe.

For the first time, a direct NLTE analysis of carbon and nitrogen lines in the spectra of nine RR Lyrae stars was carried out. We have determined the abundances of these elements together with oxygen, and have shown that the nitrogen content is increased in metallicity deficient program stars. We conclude that this is a sign of the first dredge up, which occurred at the previous stage of the red giant branch, and brought material processed in an incomplete CNO cycle to the surface of the star. This effect is significantly enhanced by thermohaline (extra-) mixing, which is more effective for metal-poor RR Lyrae stars. This is clearly seen in the plot showing that C/N ration in our sample of stars gradually decreases as metallicity decreases from about --0.2 to --2. Oxygen abundance depends on metallicity in a similar way to what we see in the Population II stars., Comment: 12 pages, 8 figures Accepted by Astronomische Nachrichten / Astronomical Notes

We present the radial gas-phase, mass-weighted metallicity profiles and gradients of the TNG50 star-forming galaxy population measured at redshifts $z=$ 0--3. We investigate the redshift evolution of gradients and examine relations between gradient steepness and galaxy properties. We find that TNG50 gradients are predominantly negative at all redshifts, although we observe significant diversity among these negative gradients. We determine that the gradient steepness of all galaxies increases approximately monotonically with redshift at a roughly constant rate. This rate does not vary significantly with galaxy mass. We observe a weak negative correlation between gradient steepness and galaxy stellar mass at redshifts $z\leq2$. However, when we normalize gradients by a characteristic radius defined by the galactic star formation distribution, we find that these normalized gradients remain invariant with both stellar mass and redshift. We place our results in the context of previous simulations and show that TNG50 high-redshift gradients are steeper than those of models featuring burstier feedback, which may further highlight high-redshift gradients as important discriminators of galaxy formation models. We also find that redshift $z=0$ and $z=0.5$ TNG50 gradients are consistent with the gradients observed in galaxies at these redshifts, although the preference for flat gradients observed in redshift $z\gtrsim1$ galaxies is not present in TNG50. If future JWST and ELT observations validate these flat gradients, it may indicate a need for simulation models to implement more powerful radial gas mixing within the ISM, possibly via turbulence and/or stronger winds, 24 pages, 14 figures, submitted to MNRAS

The flux ratios of high-ionization lines are commonly assumed to indicate the metallicity of the broad emission line region in luminous quasars. When accounting for the variation in their kinematic profiles, we show that the NV/CIV, (SiIV+OIV])/CIV and NV/Lya line ratios do not vary as a function of the quasar continuum luminosity, black hole mass, or accretion rate. Using photoionization models from CLOUDY , we further show that the observed changes in these line ratios can be explained by emission from gas with solar abundances, if the physical conditions of the emitting gas are allowed to vary over a broad range of densities and ionizing fluxes. The diversity of broad line emission in quasar spectra can be explained by a model with emission from two kinematically distinct regions, where the line ratios suggest that these regions have either very different metallicity or density. Both simplicity and current galaxy evolution models suggest that near-solar abundances, with parts of the spectrum forming in high-density clouds, are more likely. Within this paradigm, objects with stronger outflow signatures show stronger emission from gas which is denser and located closer to the ionizing source, at radii consistent with simulations of line-driven disc-winds. Studies using broad-line ratios to infer chemical enrichment histories should consider changes in density and ionizing flux before estimating metallicities., 12 pages, 7 figures, accepted for publication in MNRAS. The interactive version of Figs. 3 and 4 will be made available as supplementary online material with the journal, and will be available at https://people.ast.cam.ac.uk/~alr53/CIV_yvar/all_plots.html until August 2021

This paper outlines the science case for line-intensity mapping with a space-borne instrument targeting the sub-millimeter (microwaves) to the far-infrared (FIR) wavelength range. Our goal is to observe and characterize the large-scale structure in the Universe from present times to the high redshift Epoch of Reionization. This is essential to constrain the cosmology of our Universe and form a better understanding of various mechanisms that drive galaxy formation and evolution. We argue that the proposed frequency range would make it possible to probe important metal cooling lines such as [CII] up to very high redshift as well as a large number of rotational lines of the CO molecule. These can be used to trace molecular gas and dust evolution and constrain the buildup in both the cosmic star formation rate density and the cosmic infrared background (CIB). Moreover, surveys at the highest frequencies will detect FIR lines which are used as diagnostics of galaxies and AGN. Tomography of these lines over a wide redshift range will enable invaluable measurements of the cosmic expansion history at epochs inaccessible to other methods, competitive constraints on the parameters of the standard model of cosmology, and numerous tests of dark matter, dark energy, modified gravity and inflation. To reach these goals, large-scale structure must be mapped over a wide range in frequency to trace its time evolution over a reasonable fraction of the volume of the observable Universe. In addition, the surveyed area needs to be very large to beat cosmic variance and to probe the largest scales where its easier to separate the astrophysical and cosmological contributions to the observed signal. Only, a space-borne mission can properly meet these requirements., 25 pages, 6 figures, ESA Voyage-2050 White Paper

We propose a dedicated, ultra-deep spectroscopic survey in the near infrared (NIR), that will target a mass-limited sample of galaxies during two of the most fundamental epochs of cosmic evolution: the formation of the first galaxies (at z $\gtrsim $ ≳ 6; cosmic dawn), and at the peak of galaxy formation activity (at redshift z$\sim $ ∼ 1–3; cosmic noon). By way of NIR observations (λ= 0.8–2μ m), it is possible to study the UV Lyman-α region in the former, and the optical rest-frame in the latter, allowing us to extract fundamental observables such as gas and stellar kinematics, chemical abundances, and ages, providing a unique legacy database covering these two crucial stages of cosmic evolution. The need to work in the NIR at extremely low flux levels makes a ground-based approach unfeasible due to atmospheric emission and absorption. Only with the largest facilities of the future (e.g. ELT) will be possible to observe a reduced set of targets, comprising at most of order thousands of galaxies. Likewise, from space, the small field of view of JWST and its use as a general purpose facility will yield a rather small set of high quality NIR spectra of distant galaxies (in the thousands, at best). Our project (codename ) aims to produce $\sim $ ∼ 1 million high quality spectra, with a high S/N in the continuum, where information about the underlying stellar populations is encoded. The main science drivers are: The connection between the star formation history and the mass assembly history. The role of AGN and supernova feedback in shaping the formation histories of galaxies, with a quantitative estimate of quenching timescales. The formation of the first galaxies. The source of reionization. Evolution of the metallicity-mass relation, including [α/Fe] and individual abundances. Precision cosmology through detailed studies of the “baryon physics” of galaxy formation, probing the power spectrum over scales k$\sim $ ∼ 1 Mpc− 1.

The shapes of galaxy N-point correlation functions can be used as standard rulers to constrain the distance-redshift relationship and thence the expansion rate of the Universe. The cosmological density fields traced by late-time galaxy formation are initially nearly Gaussian, and hence all the cosmological information can be extracted from their 2-Point Correlation Function (2PCF) or its Fourier-space analog the power spectrum. Subsequent nonlinear evolution under gravity, as well as halo and then galaxy formation, generate higher-order correlation functions. Since the mapping of the initial to the final density field is, on large scales, invertible, it is often claimed that the information content of the initial field's power spectrum is equal to that of all the higher-order functions of the final, nonlinear field. This claim implies that reconstruction of the initial density field from the nonlinear field renders analysis of higher-order correlation functions of the latter superfluous. We here show that this claim is false when the N-point functions are used as standard rulers. Constraints available from joint analysis of the galaxy power spectrum and bispectrum (Fourier-space analog of the 3-Point Correlation Function) can, in some cases, exceed those offered by the initial power spectrum even when the reconstruction is perfect. We provide a mathematical justification for this claim and also demonstrate it using a large suite of N-body simulations. In particular, we show that for the z = 0 real-space matter field in the limit of vanishing shot noise, taking modes up to k_max = 0.2 h/Mpc, using the bispectrum alone offers a factor of two reduction in the variance on the cosmic distance scale relative to that available from the power spectrum., 14 pages, 5 figures. Summitted to the Monthly Notices of the Royal Astronomical Society

Mergers play an important role in galaxy evolution. In particular, major mergers are able to have a transformative effect on galaxy morphology. In this paper, we investigate the role of magnetic fields in gas-rich major mergers. To this end, we run a series of high-resolution magnetohydrodynamic (MHD) zoom-in simulations with the moving-mesh code Arepo and compare the outcome with hydrodynamic simulations run from the same initial conditions. This is the first time that the effect of magnetic fields in major mergers has been investigated in a cosmologically-consistent manner. In contrast to previous non-cosmological simulations, we find that the inclusion of magnetic fields has a substantial impact on the production of the merger remnant. Whilst magnetic fields do not strongly affect global properties, such as the star formation history, they are able to significantly influence structural properties. Indeed, MHD simulations consistently form remnants with extended discs and well-developed spiral structure, whilst hydrodynamic simulations form more compact remnants that display distinctive ring morphology. We support this work with a resolution study and show that whilst global properties are broadly converged across resolution and physics models, morphological differences only develop given sufficient resolution. We argue that this is due to the more efficient excitement of a small-scale dynamo in higher resolution simulations, resulting in a more strongly amplified field that is better able to influence gas dynamics., Comment: 27 pages, 18 figures. Published in MNRAS on 19 May 2021

Voids may affect galaxy formation via weakening mass infall or increasing disk sizes, which could potentially play a role in the formation of giant low surface brightness galaxies (LSBGs). If a dark matter halo forms at the potential hill corresponding to a void of the cosmic web, which we denote the 'elaphrocentre' in contrast to a barycentre, then the elaphrocentre should weaken the infall rate to the halo when compared to infall rates towards barycentres. We investigate this hypothesis numerically. We present a complete software pipeline to simulate galaxy formation, starting from a power spectrum of initial perturbations and an N-body simulation through to merger-history-tree based mass infall histories. The pipeline is built from well-established, free-licensed cosmological software packages, and aims at highly portable long-term reproducibility. We find that the elaphrocentric accelerations tending to oppose mass infall are modest. We do not find evidence of location in a void or elaphrocentric position weakening mass infall towards a galaxy. However, we find indirect evidence of voids influencing galaxy formation: while void galaxies are of lower mass compared to galaxies in high density environments, their spin parameters are typically higher. For a fixed mass, the implied disk scale length would be greater. Tangential accelerations in voids are found to be high and might significantly contribute to the higher spin parameters. We find significantly later formation epochs for void galaxies; this should give lower matter densities and may imply lower surface densities of disk galaxies. Thus, void galaxies have higher spin parameters and later formation epochs; both are factors that may increase the probability of forming LSBGs in voids., 16 pages, 16 figures, 4 tables, reproducibility: zenodo.4699702 at https://zenodo.org/record/4699702, git: https://codeberg.org/boud/elaphrocentre ; v3: matches accepted version

We present the main sequence for all galaxies and star-forming galaxies for a sample of 28,469 massive ($M_\star \ge 10^{11}$M$_\odot$) galaxies at cosmic noon ($1.5 < z < 3.0$), uniformly selected from a 17.5 deg$^2$ area (0.33 Gpc$^3$ comoving volume at these redshifts). Our large sample allows for a novel approach to investigating the galaxy main sequence that has not been accessible to previous studies. We measure the main sequence in small mass bins in the SFR-M$_{\star}$ plane without assuming a functional form for the main sequence. With a large sample of galaxies in each mass bin, we isolate star-forming galaxies by locating the transition between the star-forming and green valley populations in the SFR-M$_{\star}$ plane. This approach eliminates the need for arbitrarily defined fixed cutoffs when isolating the star-forming galaxy population, which often biases measurements of the scatter around the star-forming galaxy main sequence. We find that the main sequence for all galaxies becomes increasingly flat towards present day at the high-mass end, while the star-forming galaxy main sequence does not. We attribute this difference to the increasing fraction of the collective green valley and quiescent galaxy population from $z=3.0$ to $z=1.5$. Additionally, we measure the total scatter around the star-forming galaxy main sequence and find that it is $\sim0.5-1.0$ dex with little evolution as a function of mass or redshift. We discuss the implications that these results have for pinpointing the physical processes driving massive galaxy evolution., 16 pages, 10 figures, Accepted for publication in MNRAS

We present a comparison of the physical properties of the ionized gas in the circumgalactic medium and intergalactic medium (IGM) at z ∼ 0 between observations and four cosmological hydrodynamical simulations: Illustris, TNG300 of the IllustrisTNG project, EAGLE, and one of the Magneticum simulations. For the observational data, we use the gas properties that are inferred from cross-correlating the Sunyaev–Zel’dovich effect (SZE) from the Planck CMB maps with haloes and large-scale structure. Both the observational and simulation results indicate that the integrated gas pressure in haloes deviates from the self-similar case, showing that feedback impacts haloes with $M_{500}\sim 10^{12\!-\!13}\, {\rm M_\odot }$. The simulations predict that more than half the baryons are displaced from haloes, while the gas fraction inferred from our observational data roughly equals the cosmic baryon fraction throughout the $M_{500}\sim 10^{12\!-\!14.5}\, {\rm M_\odot }$ halo mass range. All simulations tested here predict that the mean gas temperature in haloes is about the virial temperature, while that inferred from the SZE is up to one order of magnitude lower than that from the simulations (and also from X-ray observations). While a remarkable agreement is found for the average properties of the IGM between the observation and some simulations, we show that their dependence on the large-scale tidal field can break the degeneracy between models that show similar predictions otherwise. Finally, we show that the gas pressure and the electron density profiles from simulations are not well described by a generalized NFW profile. Instead, we present a new model with a mass-dependent shape that fits the profiles accurately.

We study the radial motions of cold, star-forming gas in the secular evolution phase of a set of 14 magnetohydrodynamical cosmological zoom-in simulations of Milky Way-mass galaxies. We study the radial transport of material within the disc plane in a series of concentric rings. For the gas in each ring at a given time we compute two quantities as a function of time and radius: 1) the radial bulk flow of the gas; and 2) the radial spread of the gas relative to the bulk flow. Averaging the data from all the halos, we find that the radial spread increases with radius in the form of a power law with strong secondary dependencies on the fraction of accreted material and the local radial velocity dispersion of the gas. We find that the bulk motion of gas is well described in the inner disc regions by a radially-independent mean inward flow speed of $-$2.4 km s$^{-1}$. The spread around this value relates to the change in angular momentum of the gas and also the amount of accreted material. These scalings from fully cosmological, MHD simulations of galaxy formation can then be used in semi-analytic models to better parameterise the radial flow of gas in discs., 17 pages, 13 figures, Accepted for publication in MNRAS