Understanding the galaxy-halo connection is fundamental for contemporary models of galaxy clustering. The extent to which the haloes' assembly history and environment impact galaxy clustering (a.k.a. galaxy assembly bias; GAB), remains a complex and challenging problem. Using a semi-analytic galaxy formation model, we study the individual contributions of different secondary halo properties to the GAB signal. These are obtained by comparing the clustering of stellar-mass selected samples to that of shuffled samples where the galaxies are randomly reassigned to haloes of fixed mass and a specified secondary halo property. We explore a large range of internal halo properties and environmental measures. We find that commonly-used properties like the halo age or concentration amount to only 20-30% of the signal, while the smoothed matter density or the tidal anisotropy can explain the full level of GAB (though care should be given to the specific definition). For the ``successful'' measures, we examine the occupancy variations and the associated changes in the halo occupation function parameters. These are used to create mock catalogues which reproduce the full level of GAB. Finally, we propose a practical modification of the standard halo occupation distribution model, which can be tuned to any level of assembly bias. Fitting the parameters to our semi-analytic model, we demonstrate that the corresponding mock catalogue recovers the target level of GAB as well as the occupancy variations. Our results enable producing realistic mock catalogues and directly inform theoretical modelling of assembly bias and attempts to detect it in the Universe.

We use the {\sc Illustris TNG300} magneto-hydrodynamic simulation, the {\sc SAGE} semi-analytical model, and the subhalo abundance matching technique (SHAM) to examine the diversity in predictions for galaxy assembly bias. We consider samples of galaxies selected according to their stellar mass or star formation rate at various redshifts. We find that all models predict an assembly bias signal of different magnitude, redshift evolution, and dependence with selection criteria and number density. To model these non-trivial dependences, we propose an extension to the standard SHAM technique so it can include arbitrary amounts of assembly bias. We do this by preferentially selecting subhaloes with the same internal property but different {\it individual} large-scale bias. We find that with this model, we can successfully reproduce the galaxy assembly bias signal in either {\sc SAGE} or the {\sc TNG}, for all redshifts and galaxy number densities. We anticipate that this model can be used to constrain the level of assembly bias in observations and aid in the creation of more realistic mock galaxy catalogues.

We present the BACCO project, a simulation framework specially designed to provide highly-accurate predictions for the distribution of mass, galaxies, and gas as a function of cosmological parameters. In this paper, we describe our main suite of simulations (L $\sim2$ Gpc and $4320^3$ particles) and present various validation tests. Using a cosmology-rescaling technique, we predict the nonlinear mass power spectrum over the redshift range 0 < z < 1.5 and over scales $10^{-2} < k/(h Mpc^{-1} ) < 5$ for 800 points in an 8-dimensional cosmological parameter space. For an efficient interpolation of the results, we build a Gaussian emulator which we test against other widely-used methods. Over the whole range of scales considered, we expect our predictions to be accurate at the 2% level for parameters in the minimal ${\Lambda}$CDM model and to 3% when extended to dynamical dark energy and massive neutrinos. We make our emulator publicly available under http://www.dipc.org/bacco

We use the EAGLE simulations to study the effects of the intra-cluster medium (ICM) on the spatially resolved star-formation activity in galaxies. We study three cases of galaxy asymmetry dividing each galaxy in two halves using the plane (i) perpendicular to the $\rm \texttt{velocity}$ direction, differentiating the galaxy part approaching to the cluster center, hereafter dubbed as the "leading half", and the opposite one "trailing half", (ii) perpendicular to the $\rm \texttt{radial}$ position of the satellite to the centre of the cluster, (iii) that maximizes the star-formation rate ($\rm SFR$) difference between the two halves. For (i), we find an enhancement of the $\rm SFR$, star formation efficiency ($\rm SFE$), and interstellar medium pressure in the leading half with respect to the trailing one and normal star-forming galaxies in the EAGLE simulation, and a clear overabundance of gas particles in their trailing. These results suggest that ram pressure (RP) is boosting the star formation by gas compression in the leading half, and transporting the gas to the trailing half. This effect is more pronounced in satellites of intermediate stellar masses $\rm 10^{9.5-10.5} M_{\odot}$, with gas masses above $\rm 10^{9} M_{\odot}$, and located within one virial radius or in the most massive clusters. In (iii) we find an alignment between the velocity and the vector perpendicular to the plane that maximizes the $\rm SFR$ difference between the two halves. It suggests that finding this plane in real galaxies can provide insights into the velocity direction.

The interpretation of cosmological observables requires the use of increasingly sophisticated theoretical models. Since these models are becoming computationally very expensive and display non-trivial uncertainties, the use of standard Bayesian algorithms for cosmological inferences, such as MCMC, might become inadequate. Here, we propose a new approach to parameter estimation based on an iterative Gaussian emulation of the target likelihood function. This requires a minimal number of likelihood evaluations and naturally accommodates for stochasticity in theoretical models. We apply the algorithm to estimate 9 parameters from the monopole and quadrupole of a mock power spectrum in redshift space. We obtain accurate posterior distribution functions with approximately 100 times fewer likelihood evaluations than an affine invariant MCMC, roughly independently from the dimensionality of the problem. We anticipate that our parameter estimation algorithm will accelerate the adoption of more accurate theoretical models in data analysis, enabling more comprehensive exploitation of cosmological observables.

Current and future cosmological surveys are targeting star-forming galaxies at z~1 with nebular emission lines. In this work we use a state-of-the-art model of galaxy formation and evolution to explore the large scale environment of star-forming emission line galaxies (ELGs). Model galaxies are selected with cuts such that the samples can be directly compared with the DEEP2, VVDS, eBOSS-SGC and DESI surveys. Their large scale environment have been classified using a velocity-shear-tensor and a tidal-tensor algorithms. Half of model ELGs live in filaments and about a third in sheets. Model ELGs in knots have the largest satellite fractions. We find that the shape of the mean halo occupation distribution of model ELGs varies widely for different large scale environments. To put these results in context, we have compared fixed number density samples of galaxies, generated by either imposing an extra cut in stellar mass, star formation rate (SFR) or L[OII] to the ELGs, or by imposing a single cut in these quantities to the total model galaxy population. For the fixed number density samples, ELGs are close to L[OII] and SFR selected samples for densities above $10^{-4.2}h^{3}{\rm Mpc}^{-3}$. ELGs with an extra cut in stellar mass to fix their number density, present differences in sheets and knots with respect to the other samples. ELGs, SFR and L[OII] selected samples with equal number density, have similar large scale bias but their clustering below separations of $1h^{-1}$Mpc is different.

We present and test a framework that models the three-dimensional distribution of mass in the Universe as a function of cosmological and astrophysical parameters. Our approach combines two different techniques: a rescaling algorithm that modifies the cosmology of gravity-only N-body simulations, and a baryonification algorithm which mimics the effects of astrophysical processes induced by baryons, such as star formation and AGN feedback. We show how this approach can accurately reproduce the effects of baryons on the matter power spectrum of various state-of-the-art hydro-dynamical simulations (EAGLE, Illustris, Illustris-TNG, Horizon-AGN, and OWLS,Cosmo-OWLS and BAHAMAS), to percent level from very large down to small, highly nonlinear scales, k= 5 h/Mpc, and from z=0 up to z=2. We highlight that, thanks to the heavy optimisation of the algorithms, we can obtain these predictions for arbitrary baryonic models and cosmology (including massive neutrinos and dynamical dark energy models) with an almost negligible CPU cost. Therefore, this approach is efficient enough for cosmological data analyses. With these tools in hand we explore the degeneracies between cosmological and astrophysical parameters in the nonlinear mass power spectrum. Our findings suggest that after marginalising over baryonic physics, cosmological constraints inferred from weak gravitational lensing should be moderately degraded.

Predicting the spatial distribution of objects as a function of cosmology is an essential ingredient for the exploitation of future galaxy surveys. In this paper we show that a specially-designed suite of gravity-only simulations together with cosmology-rescaling algorithms can provide the clustering of dark matter, haloes, and subhaloes with high precision. Specifically, with only 3 $N$-body simulations we obtain the power spectrum of dark matter at $z=0$ and $z=1$ to better than 3\% precision for essentially all currently viable values of 8 cosmological parameters, including massive neutrinos and dynamical dark energy, and over the whole range of scales explored, 0.03 < $k/h^{-1}Mpc$ < 5. This precision holds at the same level for mass-selected haloes and for subhaloes selected according to their peak maximum circular velocity. As an initial application of these predictions, we successfully constrain $\Omega_{\rm m}$, $\sigma_8$, and the scatter in subhalo-abundance-matching employing the projected correlation function of mock SDSS galaxies.

Providing accurate predictions for the spatial distribution of matter and luminous tracers in the presence of massive neutrinos is an important task, given the imminent arrival of highly accurate large-scale structure observations. In this work, we address this challenge by extending cosmology-rescaling algorithms to massive neutrino cosmologies. In this way, a $\Lambda$CDM simulation can be modified to provide nonlinear structure formation predictions in the presence a hot component of arbitrary mass, and, if desired, to include non-gravitational modifications to the clustering of matter on large scales. We test the accuracy of the method by comparing its predictions to a suite of simulations carried out explicitly including a neutrino component in its evolution equations. We find that, for neutrino masses in the range $M_\nu \in [0.06, 0.3] ~ \mathrm{eV}$ the matter power spectrum is recovered to better than $1\%$ on all scales $k<2~h~\mathrm{Mpc}^{-1}$. Similarly, the halo mass function is predicted at a few percent level over the range $M_{\rm halo} \in [10^{12}, 10^{15}] ~ h^{-1} ~ \mathrm{M}_{\odot}$, and so do also the multipoles of the galaxy 2-point correlation function in redshift space over $r \in [0.1, 200] ~ h^{-1} ~ \mathrm{Mpc}$. We provide parametric forms for the necessary transformations, as a function of $\Omega_{\rm m}$ and $\Omega_{\nu}$ for various target redshifts.

We investigate a conceptual modification of the halo occupation distribution approach, using the halos' present-day maximal circular velocity, $\vmax$, as an alternative to halo mass. In particular, using a semi-analytic galaxy formation model applied to the Millennium WMAP7 simulation, we explore the extent that switching to $\vmax$ as the primary halo property incorporates the effects of assembly bias into the formalism. We consider fixed number density galaxy samples ranked by stellar mass and examine the variations in the halo occupation functions with either halo concentration or formation time. We find that using $\vmax$ results in a significant reduction in the occupancy variation of the central galaxies, particularly for concentration. The satellites occupancy variation on the other hand increases in all cases. We find effectively no change in the halo clustering dependence on concentration, for fixed bins of $\vmax$ compared to fixed halo mass. Most crucially, we calculate the impact of assembly bias on galaxy clustering by comparing the amplitude of clustering to that of a shuffled galaxy sample, finding that the level of galaxy assembly bias remains largely unchanged. Our results suggest that while using $\vmax$ as a proxy for halo mass diminishes some of occupancy variations exhibited in the galaxy-halo relation, it is not able to encapsulate the effects of assembly bias potentially present in galaxy clustering. The use of other more complex halo properties, such as $\vpeak$, the peak value of $\vmax$ over the assembly history, provides some improvement and warrants further investigation.

We test different implementations of the halo occupation distribution (HOD) model to reconstruct the spatial distribution of galaxies as predicted by a publicly available semi-analytical model (SAM). We compare the measured two-point correlation functions of the HOD mock catalogues and the SAM samples to quantify the fidelity of the reconstruction. We use fixed number density galaxy samples selected according to stellar mass or star formation rate (SFR). We develop three different schemes to populate haloes with galaxies with increasing complexity, considering the scatter of the satellite HOD as an additional parameter in the modelling. We first modify the SAM output, removing assembly bias and using a standard Navarro-Frenk-White density profile for the satellite galaxies as the target to reproduce with our HOD mocks. We find that all models give similar reproductions of the two-halo contribution to the clustering signal, but there are differences in the one-halo term. In particular, the HOD mock reproductions work equally well using either the HOD of central and satellites separately or using a model that also accounts for whether or not the haloes contain a central galaxy. We find that the HOD scatter does not have an important impact on the clustering predictions for stellar mass selected samples. For SFR selections, we obtain the most accurate results assuming a negative binomial distribution for the number of satellites in a halo. The scatter in the satellites HOD is a key consideration for HOD mock catalogues that mimic ELG or SFR selected samples in future galaxy surveys.

We examine the evolution of assembly bias using a semi-analytical model of galaxy formation implemented in the Millennium-WMAP7 N-body simulation. We consider fixed number density galaxy samples ranked by stellar mass or star formation rate. We investigate how the clustering of haloes and their galaxy content depend on halo formation time and concentration, and how these relationships evolve with redshift. At $z=0$ the dependences of halo clustering on halo concentration and formation time are similar. However, at higher redshift, halo assembly bias weakens for haloes selected by age, and reverses and increases for haloes selected by concentration. The variation of the halo occupation with concentration and formation time is also similar at $z=0$ and changes at higher redshifts. In this case, the occupancy variation with halo age stays mostly constant with redshift but decreases for concentration. Finally, we look at the evolution of assembly bias reflected in the galaxy distribution by examining the galaxy correlation functions relative to those of shuffled galaxy samples which remove the occupancy variation. This correlation functions ratio monotonically decreases with larger redshift and for lower number density samples, going below unity in some cases, leading to reduced galaxy clustering. While the halo occupation functions themselves vary, the assembly bias trends are similar whether selecting galaxies by stellar mass or star formation rate. Our results provide further insight into the origin and evolution of assembly bias. Our extensive occupation function measurements and fits are publicly available and can be used to create realistic mock catalogues.

We study potential systematic effects of assembly bias on cosmological parameter constraints from redshift space distortion measurements. We use a semi-analytic galaxy formation model applied to the Millennium N-body WMAP-7 simulation to study the effects of halo assembly bias on the redshift space distortions of the galaxy correlation function. We look at the pairwise velocities of galaxies living in haloes with concentrations and ages in the upper and lower quintiles, and find that the velocity differences between these are consistent with those reported for real-space clustering analyses, i.e. samples with higher clustering also exhibit stronger infall pairwise motions. This can also be seen in the monopole and quadrupole of the redshift-space correlation function. We find that regardless of the method of measurement, the changes in the $\beta$ parameter due to different secondary halo parameters fully tracks the change in the bias Parameter. Hence, assembly bias does not introduce detectable systematics in the inferred logarithmic growth factor.

We investigate the variations in galaxy occupancy of the dark matter haloes with the large-scale environment and halo formation time, using two state-of-the-art hydrodynamical cosmological simulations, EAGLE and Illustris. For both simulations, we use three galaxy samples with a fixed number density ranked by stellar mass. For these samples we find that low-mass haloes in the most dense environments are more likely to host a central galaxy than those in the least dense environments. When splitting the halo population by formation time, these relations are stronger. Hence, at a fixed low halo mass, early-formed haloes are more likely to host a central galaxy than late-formed haloes since they have had more time to assemble. The satellite occupation shows a reverse trend where early-formed haloes host fewer satellites due to having more time to merge with the central galaxy. We also analyse the stellar mass -- halo mass relation for central galaxies in terms of the large-scale environment and formation time of the haloes. We find that low mass haloes in the most dense environment host relatively more massive central galaxies. This trend is also found when splitting the halo population by age, with early-formed haloes hosting more massive galaxies. Our results are in agreement with previous findings from semi-analytical models, providing robust predictions for the occupancy variation signature in the halo occupation distribution of galaxy formation models.

We study the correlation between the specific star formation rate of central galaxies and neighbour galaxies, also known as 'galactic conformity', out to 20 Mpc/h using three semi-analytic models (SAMs, one from L-GALAXIES and other two from GALFORM). The aim is to establish whether SAMs are able to show galactic conformity using different models and selection criteria. In all the models, when the selection of primary galaxies is based on an isolation criterion in real space, the mean fraction of quenched galaxies around quenched primary galaxies is higher than that around star-forming primary galaxies of the same stellar mass. The overall signal of conformity decreases when we remove satellites selected as primary galaxies, but the effect is much stronger in GALFORM models compared with the L-GALAXIES model. We find this difference is partially explained by the fact that in GALFORM once a galaxy becomes a satellite remains as such, whereas satellites can become centrals at a later time in L-GALAXIES. The signal of conformity decreases down to 60% in the L-GALAXIES model after removing central galaxies that were ejected from their host halo in the past. Galactic conformity is also influenced by primary galaxies at fixed stellar mass that reside in dark matter haloes of different masses. Finally, we explore a proxy of conformity between distinct haloes. In this case the conformity is weak beyond ~ 3 Mpc/h (<3% in L-GALAXIES, <1-2% in GALFORM models). Therefore, it seems difficult that conformity is directly related with a long-range effect.

Emission line galaxies (ELGs) are used in several ongoing and upcoming surveys (SDSS-IV/eBOSS, DESI) as tracers of the dark matter distribution. Using a new galaxy formation model, we explore the characteristics of [OII] emitters, which dominate optical ELG selections at $z\simeq 1$. Model [OII] emitters at $0.5<z<1.5$ are selected to mimic the DEEP2, VVDS, eBOSS and DESI surveys. The luminosity functions of model [OII] emitters are in reasonable agreement with observations. The selected [OII] emitters are hosted by haloes with $M_{\rm halo}\geq 10^{10.3}h^{-1}{\rm M}_{\odot}$, with ~90% of them being central star-forming galaxies. The predicted mean halo occupation distributions of [OII] emitters has a shape typical of that inferred for star-forming galaxies, with the contribution from central galaxies, $\langle N \rangle_{\left[OII\right]\, cen}$, being far from the canonical step function. The $\langle N \rangle_{\left[OII\right]\, cen}$ can be described as the sum of an asymmetric Gaussian for disks and a step function for spheroids, which plateaus below unity. The model [OII] emitters have a clustering bias close to unity, which is below the expectations for eBOSS and DESI ELGs. At $z\sim 1$, a comparison with observed g-band selected galaxy, which are expected to be dominated by [OII] emitters, indicates that our model produces too few [OII] emitters that are satellite galaxies. This suggests the need to revise our modelling of hot gas stripping in satellite galaxies.

We use EAGLE to quantify the effect galaxy mergers have on the stellar specific angular momentum of galaxies, $j_{\rm stars}$. We split mergers into: dry (gas-poor)/wet (gas-rich), major/minor, and different spin alignments and orbital parameters. Wet (dry) mergers have an average neutral gas-to-stellar mass ratio of $1.1$ ($0.02$), while major (minor) mergers are those with stellar mass ratios $\ge 0.3$ ($0.1-0.3$). We correlate the positions of galaxies in the $j_{\rm stars}$-stellar mass plane at $z=0$ with their merger history, and find that galaxies of low spins suffered dry mergers, while galaxies of normal/high spins suffered predominantly wet mergers, if any. The radial $j_{\rm stars}$ profiles of galaxies that went through dry mergers are deficient by $\approx 0.3$~dex at $r\lesssim 10\,r_{50}$ (with $r_{50}$ being the half-stellar mass radius), compared to galaxies that went through wet mergers. Studying the merger remnants reveals that dry mergers reduce $j_{\rm stars}$ by $\approx 30$\%, while wet mergers increase it by $\approx 10$\%, on average. The latter is connected to the build-up of the bulge by newly formed stars of high rotational speed. Moving from minor to major mergers accentuates these effects. When the spin vectors of the galaxies prior to the dry merger are misaligned, $j_{\rm stars}$ decreases to a greater magnitude, while in wet mergers co-rotation and high orbital angular momentum efficiently spun-up galaxies. We predict what would be the observational signatures in the $j_{\rm stars}$ profiles driven by dry mergers: (i) shallow radial profiles and (ii) profiles that rise beyond $\approx 10\,r_{50}$, both of which are significantly different from spiral galaxies.

We use three different cosmological dark matter simulations to study how the orientation of the angular momentum vector (AM) in dark matter haloes evolve with time. We find that haloes in this kind of simulations are constantly affected by a spurious change of mass, which translates into an artificial change in the orientation of the AM. After removing the haloes affected by artificial mass change, we found that the change in the orientation of the AM vector is correlated with time. The change in its angle and direction (i.e. the angle subtended by the AM vector in two consecutive timesteps) that affect the AM vector has a dependence on the change of mass that affects a halo, the time elapsed in which the change of mass occurs and the halo mass. We create a Monte-Carlo simulation that reproduces the change of angle and direction of the AM vector. We reproduce the angular separation of the AM vector since a look back time of 8.5 Gyrs to today ( $\rm \alpha$) with an accuracy of approximately 0.05 in $\rm cos(\alpha)$. We are releasing this Monte-Carlo simulation together with this publication. We also create a Monte Carlo simulation that reproduces the change of the AM modulus. We find that haloes in denser environments display the most dramatic evolution in their AM direction, as well as haloes with a lower specific AM modulus. These relations could be used to improve the way we follow the AM vector in low-resolution simulations.

We study the dependence of the galaxy content of dark matter halos on large-scale environment and halo formation time using semi-analytic galaxy models applied to the Millennium simulation. We analyze subsamples of halos at the extremes of these distributions and measure the occupation functions for the galaxies they host. We find distinct differences in these occupation functions. The main effect with environment is that central galaxies (and in one model also the satellites) in denser regions start populating lower-mass halos. A similar, but significantly stronger, trend exists with halo age, where early-forming halos are more likely to host central galaxies at lower halo mass. We discuss the origin of these trends and the connection to the stellar mass -- halo mass relation. We find that, at fixed halo mass, older halos and to some extent also halos in dense environments tend to host more massive galaxies. Additionally, we see a reverse trend for the satellite galaxies occupation where early-forming halos have fewer satellites, likely due to having more time for them to merge with the central galaxy. We describe these occupancy variations also in terms of the changes in the occupation function parameters, which can aid in constructing realistic mock galaxy catalogs. Finally, we study the corresponding galaxy auto- and cross-correlation functions of the different samples and elucidate the impact of assembly bias on galaxy clustering. Our results can inform theoretical models of assembly bias and attempts to detect it in the real universe.

We use the hydrodynamic, cosmological EAGLE simulations to investigate how hot gas in haloes condenses to form and grow galaxies. We select haloes from the simulations that are actively cooling and study the temperature, distribution, and metallicity of their hot, cold, and transitioning `cooling' gas, placing these in context of semi-analytic models. Our selection criteria lead us to focus on Milky Way-like haloes. We find the hot-gas density profiles of the haloes form a progressively stronger core over time, the nature of which can be captured by a beta profile that has a simple dependence on redshift. In contrast, the hot gas that will cool over a time-step is broadly consistent with a singular isothermal sphere. We find that cooling gas carries a few times the specific angular momentum of the halo and is offset in spin direction from the rest of the hot gas. The gas loses ~60% of its specific angular momentum during the cooling process, generally remaining greater than that of the halo, and it precesses to become aligned with the cold gas already in the disc. We find tentative evidence that angular-momentum losses are slightly larger when gas cools onto dispersion-supported galaxies. We show that an exponential surface density profile for gas arriving on a disc remains a reasonable approximation, but a cusp containing ~20% of the mass is always present, and disc scale radii are larger than predicted by a vanilla Fall & Efstathiou model. These scale radii are still closely correlated with the halo spin parameter, for which we suggest an updated prescription for galaxy formation models.

We use the halo occupation distribution (HOD) framework to characterise the predictions from two independent galaxy formation models for the galactic content of dark matter haloes and its evolution with redshift. Our galaxy samples correspond to a range of fixed number densities defined by stellar mass and span $0 \le z \le 3$. We find remarkable similarities between the model predictions. Differences arise at low galaxy number densities which are sensitive to the treatment of heating of the hot halo by active galactic nuclei. The evolution of the form of the HOD can be described in a relatively simple way, and we model each HOD parameter using its value at $z=0$ and an additional evolutionary parameter. In particular, we find that the ratio between the characteristic halo masses for hosting central and satellite galaxies can serve as a sensitive diagnostic for galaxy evolution models. Our results can be used to test and develop empirical studies of galaxy evolution and can facilitate the construction of mock galaxy catalogues for future surveys.

We use the EAGLE cosmological hydrodynamic simulation suite to study the specific angular momentum of galaxies, $j$, with the aims of (i) investigating the physical causes behind the wide range of $j$ at fixed mass and (ii) examining whether simple, theoretical models can explain the seemingly complex and non-linear nature of the evolution of $j$. We find that $j$ of the stars, $j_{\rm stars}$, and baryons, $j_{\rm bar}$, are strongly correlated with stellar and baryon mass, respectively, with the scatter being highly correlated with morphological proxies such as gas fraction, stellar concentration, (u-r) intrinsic colour, stellar age and the ratio of circular velocity to velocity dispersion. We compare with available observations at $z=0$ and find excellent agreement. We find that $j_{\rm bar}$ follows the theoretical expectation of an isothermal collapsing halo under conservation of specific angular momentum to within $\approx 50$%, while the subsample of rotation-supported galaxies are equally well described by a simple model in which the disk angular momentum is just enough to maintain marginally stable disks. We extracted evolutionary tracks of the stellar spin parameter of EAGLE galaxies and found that the fate of their $j_{\rm stars}$ at $z=0$ depends sensitively on their star formation and merger histories. From these tracks, we identified two distinct physical channels behind low $j_{\rm stars}$ galaxies at $z=0$: (i) galaxy mergers, and (ii) early star formation quenching. The latter can produce galaxies with low $j_{\rm stars}$ and early-type morphologies even in the absence of mergers.

We demonstrate how the properties of a galaxy depend on the mass of its host dark matter subhalo, using two independent models of galaxy formation. For the cases of stellar mass and black hole mass, the median property value displays a monotonic dependence on subhalo mass. The slope of the relation changes for subhalo masses for which heating by active galactic nuclei becomes important. The median property values are predicted to be remarkably similar for central and satellite galaxies. The two models predict considerable scatter around the median property value, though the size of the scatter is model dependent. There is only modest evolution with redshift in the median galaxy property at a fixed subhalo mass. Properties such as cold gas mass and star formation rate, however, are predicted to have a complex dependence on subhalo mass. In these cases subhalo mass is not a good indicator of the value of the galaxy property. We illustrate how the predictions in the galaxy property - subhalo mass plane differ from the assumptions made in empirical models of galaxy clustering by reconstructing the model output using a subhalo abundance matching scheme. In its simplest form, abundance matching generally does not reproduce the clustering predicted by the models, typically resulting in an overprediction of the clustering signal. We show how the basic abundance matching scheme can be extended to reproduce the model predictions more faithfully, which has implications for the analysis of galaxy clustering, particularly for low abundance samples.

We study the origin of the wide distribution of angles between the angular momenta of the stellar and gas components, $\alpha_{\rm G,S}$, in early-type galaxies (ETGs). We use the GALFORM model of galaxy formation, set in the $\Lambda$ cold dark matter framework, and coupled it with a Monte-Carlo simulation to follow the angular momenta flips driven by matter accretion onto haloes and galaxies. We consider a gas disk to be misaligned with respect to the stellar body if $\alpha_{\rm G,S}>30$~degrees. By assuming that the only sources of misaligments in galaxies are galaxy mergers, we place a lower limit of $2-5$ per cent on the fraction of ETGs with misaligned gas/stellar components. These low fractions are inconsistent with the observed value of $\approx 42\pm 6$ per cent in ATLAS$^{\rm 3D}$. In the more general case, in which smooth gas accretion in addition to galaxy mergers can drive misalignments, our calculation predicts that $\approx 46$ per cent of ETGs have $\alpha_{\rm G,S}>30$~degrees. In this calculation, we find correlations between $\alpha_{\rm G,S}$ and stellar mass, cold gas fraction and star formation rate, such that ETGs with high masses, low cold gas fractions and low star formation rates are more likely to display aligned cold gas and stellar components. We confirm these trends observationally for the first time using ATLAS$^{\rm 3D}$ data. We argue that the high fraction of misaligned gas discs observed in ETGs is mostly due to smooth gas accretion (e.g. cooling from the hot halo of galaxies) which takes place after most of the stellar mass of the galaxy is in place and comes misaligned with respect to the stellar component. Galaxies that have accreted most of their cold gas content prior to the time where most of the stellar mass was in place show aligned components.

The angular momentum of galactic discs in semi-analytic models of galaxy formation is usually updated in time as material is accreted to the disc by adopting a constant dimensionless spin parameter and little attention is paid to the effects of accretion with misaligned angular momenta. These effects are the subject of this paper, where we adopt a Monte-Carlo simulation for the changes in the direction of the angular momentum of a galaxy disc as it accretes matter based on accurate measurements from dark-matter haloes in the Millennium II simulation. In our semi-analytic model implementation, the flips seen the dark matter haloes are assumed to be the same for the cold baryons; however, we also assume that in the latter the flip also entails a difficulty for the disc to increase its angular momentum which causes the disc to become smaller relative to a no-flip case. This makes star formation to occur faster, specially in low mass galaxies at all redshifts allowing galaxies to reach higher stellar masses faster. We adopt a new condition for the triggering of starbursts during mergers. As these produce the largest flips it is natural to adopt the disc instability criterion to evaluate the triggering of bursts in mergers instead of one based on mass ratios as in the original model. The new implementation reduces the average lifetimes of discs by a factor of 2, while still allowing old ages for the present-day discs of large spiral galaxies. It also provides a faster decline of star formation in massive galaxies and a better fit to the bright end of the luminosity function at z = 0.

We study the alignment of dark matter haloes with the cosmic web characterized by the tidal and velocity shear fields. We focus on the alignment of their shape, angular momentum and peculiar velocities. We use a cosmological N-body simulation that allows to study dark matter halos spanning almost five orders of magnitude in mass ($10^{9}$-$10^{14}$) $h^{-1}$$M_{\odot}$ and spatial scales of $(0.5$-$1.0)$ $h^{-1}$ Mpc to define the cosmic web. We find that the halo shape presents the strongest alignment along the smallest tidal eigenvector, e.g. along filaments and walls, with a signal that gets stronger as the halo mass increases. In the case of the velocity shear field only massive halos $>10^{12}$ $h^{-1}$$M_{\odot}$ tend to have their shapes aligned along the largest tidal eigenvector; that is, perpendicular to filaments and walls. For the angular momentum we find alignment signals only for halos more massive than $10^{12}$ $h^{-1}$$M_{\odot}$ both in the tidal and velocity shear webs where the preferences are for it to be parallel to the middle eigenvector; perpendicular to filaments and parallel to walls. Finally, the peculiar velocities show a strong alignment along the smallest tidal eigenvector for all halo masses; halos move along filaments and walls. In the velocity shear the same alignment is present but weaker and only for haloes less massive than $10^{12}$ $h^{-1}$$M_{\odot}$. Our results clearly show that the two different algorithms we used to define the cosmic web describe different physical aspects of non-linear collapse and should be used in a complementary way to understand the effect of the cosmic web on galaxy evolution.

We use the Millennium Simulation database to compare how different versions of the Durham and Munich semi-analytical galaxy formation models populate dark matter haloes with galaxies. The models follow the same physical processes but differ in how these are implemented. All of the models we consider use the Millennium N-body Simulation; however, the Durham and Munich groups use independent algorithms to construct halo merger histories. We compare the predicted halo occupation distributions (HODs) and correlation functions for galaxy samples defined by stellar mass, cold gas mass and star formation rate. The predictions for the HOD are remarkably similar for samples ranked by stellar mass. The predicted bias averaged over 5-25 Mpc/h is consistent between models to within 10%. At small pair separations there is a difference in the predicted clustering. This arises because the Durham models allow some satellite galaxies to merge with the central galaxy in a halo when they are still associated with resolved subhaloes. The agreement between the models is less good for samples defined by cold gas mass or star formation rate, with the small scale clustering differing by an order of magnitude, reflecting the uncertainty in the modelling of star formation. The model predictions in these cases are qualitatively similar, with a markedly shallower slope for the correlation function than is found for stellar mass selected samples and with the HOD displaying an asymmetric peak for central galaxies. We provide illustrative parametric fits to the HODs predicted by the models. Our results reveal the current limitations on how well we can predict galaxy bias in a fixed cosmology, which has implications for the interpretation of constraints on the physics of galaxy formation from galaxy clustering measurements and the ability of future galaxy surveys to measure dark energy.