We study the effect of the gas accretion rate ($\dot M_{\rm accr}$) on the radial gas metallicity profile (RMP) of galaxies using the EAGLE cosmological hydrodynamic simulations, focusing on central galaxies of stellar mass $M_\star \gtrsim 10^9 \, {\rm M_\odot}$ at $z \le 1$. We find clear relations between $\dot M_{\rm accr}$ and the slope of the RMP (measured within an effective radius), where higher $\dot M_{\rm accr}$ are associated with more negative slopes. The slope of the RMPs depends more strongly on $\dot M_{\rm accr}$ than on stellar mass, star formation rate or gas fraction, suggesting $\dot M_{\rm accr}$ to be a more fundamental driver of the RMP slope of galaxies. We find that eliminating the dependence on stellar mass is essential for pinning down the properties that shape the slope of the RMP. Although $\dot M_{\rm accr}$ is the main property modulating the slope of the RMP, we find that it causes other correlations that are more easily testable observationally: at fixed stellar mass, galaxies with more negative RMP slopes tend to have higher gas fractions and SFRs, while galaxies with lower gas fractions and SFRs tend to have flatter metallicity profiles within an effective radius.

We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6<z<1.2$. We compare model predictions with different observational data sets, including DEEP2--Firefly galaxies with absolute magnitudes. We estimate the [OII] luminosity, L[OII], using simple relations derived both from the models and observations and also using a public code. This code ideally uses as input instantaneous star formation rates (SFRs), which are only provided by one of the SAMs under consideration. We use this SAM to study the feasibility of inferring galaxies' L[OII] for models that only provide average SFRs. We find that the post-processing computation of L[OII] from average SFRs is accurate for model galaxies with dust attenuated L[OII]$\lesssim10^{42.2}$erg s$^{-1}$ ($<5\%$ discrepancy). We also explore how to derive the [OII] luminosity from simple relations using global properties usually output by SAMs. Besides the SFR, the model L[OII] is best correlated with the observed-frame $u$ and $g$ broad-band magnitudes. These correlations have coefficients (r-values) above 0.64 and a dispersion that varies with L[OII]. We use these correlations and an observational one based on SFR and metallicity to derive L[OII]. These relations result in [OII] luminosity functions and halo occupation distributions with shapes that vary depending on both the model and the method used. Nevertheless, for all the considered models, the amplitude of the clustering at scales above 1$h^{-1}$Mpc remains unchanged independently of the method used to derive L[OII].

One possible and natural derivation from the collisionless cold dark matter (CDM) standard cosmological framework is the assumption of the existence of interactions between dark matter (DM) and photons or neutrinos. Such possible interacting dark matter (IDM) model would imply a suppression of small-scale structures due to a large collisional damping effect, even though the weakly interacting massive particle (WIMP) can still be the DM candidate. Because of this, IDM models can help alleviate alleged tensions between standard CDM predictions and observations at small mass scales. In this work, we investigate the properties of DM halo substructure or subhalos formed in a high-resolution cosmological N-body simulation specifically run within these alternative models. We also run its CDM counterpart, which allowed us to compare subhalo properties in both cosmologies. We show that, in the lower mass range covered by our simulation runs, both subhalo concentrations and abundances are systematically lower in IDM compared to the CDM scenario. Yet, as in CDM, we find that median IDM subhalo concentration values increase towards the innermost regions of their hosts for same mass subhalos. Also similarly to CDM, we find IDM subhalos to be more concentrated than field halos of the same mass. Our work has a direct application on studies aimed at the indirect detection of DM where subhalos are expected to boost the DM signal of their host halos significantly. From our results, we conclude that the role of halo substructure in DM searches will be less important in interacting scenarios than in CDM, but is nevertheless far from being negligible.

We use the semi-analytic model of galaxy formation SAG to study the relevance of mass and environmental quenching on satellite galaxies. We find that environmental processes dominate the star formation (SF) quenching of low-mass satellites ($M_{\star} \lesssim 10^{10.5}\, {\rm M}_{\odot}$), whereas high-mass galaxies typically quench as centrals. High-mass galaxies that remain actively forming stars while being accreted are found to be mainly affected by mass quenching after their first infall. For a given stellar mass, our model predicts SF quenching to be less efficient in low-mass haloes both before and after infall, in contradiction with common interpretations of observational data. Our model supports a two-stage scenario to explain the SF quenching. Initially, the SF of satellites resembles that of centrals until the gas cooling rate is reduced to approximately half its value at infall. Then, the SF fades through secular processes that exhaust the cold gas reservoir. This reservoir is not replenished efficiently due to the action of either ram-pressure stripping (RPS) of the hot gas in low-mass satellites, or feedback from the active galactic nucleus (AGN) in high-mass satellites. The delay times for the onset of SF quenching are found to range from $\approx 3\,{\rm Gyr}$ to $\approx 1\,{\rm Gyr}$ for low-mass ($M_{\star} \approx 10^{10}\, {\rm M}_{\odot}$) and high-mass ($M_{\star} \approx 10^{11}\, {\rm M}_{\odot}$) satellites, respectively. SF fades in $\approx 1.5\,{\rm Gyr}$, largely independent of stellar mass. We find that the SF quenching of low-mass satellites supports the so-called delay-then-rapid quenching scenario. However, the SF history of $z=0$ passive satellites of any stellar mass is better described by a delay-then-fade quenching scenario.

We present results from the semi-analytic model of galaxy formation SAG applied on the MultiDark simulation MDPL2. SAG features an updated supernova (SN) feedback scheme and a robust modelling of the environmental effects on satellite galaxies. This incorporates a gradual starvation of the hot gas halo driven by the action of ram pressure stripping (RPS), that can affect the cold gas disc, and tidal stripping (TS), which can act on all baryonic components. Galaxy orbits of orphan satellites are integrated providing adequate positions and velocities for the estimation of RPS and TS. The star formation history and stellar mass assembly of galaxies are sensitive to the redshift dependence implemented in the SN feedback model. We discuss a variant of our model that allows to reconcile the predicted star formation rate density at $z \gtrsim 3$ with the observed one, at the expense of an excess in the faint end of the stellar mass function at $z=2$. The fractions of passive galaxies as a function of stellar mass, halo mass and the halo-centric distances are consistent with observational measurements. The model also reproduces the evolution of the main sequence of star forming central and satellite galaxies. The similarity between them is a result of the gradual starvation of the hot gas halo suffered by satellites, in which RPS plays a dominant role. RPS of the cold gas does not affect the fraction of quenched satellites but it contributes to reach the right atomic hydrogen gas content for more massive satellites ($M_{\star}\gtrsim 10^{10}\,{\rm M}_{\odot}$).

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.

We use the semi-analytic model (SAM) of galaxy formation and evolution SAG coupled with the MULTIDARK simulation MDPL2 to study the evolution of the stellar mass-gas metallicity relation of galaxies (MZR). We test several implementations of the dependence of the mass loading due to supernovae (SNe). We find that no evolution in the normalization of the MZR is obtained unless we introduce an explicit scaling of the reheated and ejected mass with redshift as $(1+z)^\beta$. The latter is in agreement with results from the FIRE simulations, and it should encompass small scale properties of the interstellar medium varying over time, which are not captured in SAMs, as well as other energy sources in addition to SNe. Increasing $\beta$ leads to stronger evolution of the MZR normalization; $\beta = 1.9$ reproduces the observed MZR in the range $0 < z < 3.5$. A stronger redshift dependence of outflows reduces the levels of star formation at earlier epochs with the consequent decrease of metal production. This leads to a slower increase of the gas metallicity compared to the stellar mass build-up. The cold gas can be contaminated either by receiving a direct injection of the material recycled by stellar winds and SNe or by gas cooling. The relative role of each process for a given stellar mass depends on the criterion adopted to regulate the fate of the recycled material. However, modifying the metal loading of the outflows has mild impact on the zero-point evolution and does not affect our conclusions.

We present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into 'passive' and 'star-forming' classes and find that several of the models produce too many low-mass star-forming galaxies at high redshift compared to observations, in some cases by nearly a factor of 10 in the redshift range $2.5 < z < 3.0$. We also find important differences in the implied mass of the dark matter haloes the galaxies inhabit, by comparing with halo masses inferred from observations. Galaxies at high redshift in the models are in lower mass haloes than suggested by observations, and the star formation efficiency in low-mass haloes is higher than observed. We conclude that many of the models require a physical prescription that acts to dissociate the growth of low-mass galaxies from the growth of their dark matter haloes at high redshift.

The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs.

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available.

We study the stellar populations of bulges of Milky Way-like (MW-like) galaxies with the aim of identifying the physical processes involved in the formation of the bulge of our Galaxy. We use the semi-analytic model of galaxy formation and evolution SAG adapted to this aim; this kind of models can trace the properties of galaxies and their components like stellar discs, bulges and halos, but resolution limits prevent them from reaching the scale of stellar populations (SPs). Properties of groups of stars formed during single star formation events are stored and tracked in the model and results are compared with observations of stars in the galactic bulge. MW-like galaxies are selected using two different criteria. One of them considers intrinsic photo-metric properties and the second is focused on the cosmological context of the local group of galaxies (LG). We compare our model results with spectroscopic and photometric stellar metallicity distributions. We find that 87% of stars in bulges of MWtype galaxies in our model are accreted and formed in starbursts during disc instability events. Mergers contribute to 13% of the mass budget of the bulge and are responsible for the low metallicity tail of the distribution. Abundance ratios of {\alpha} elements with respect to iron, [{\alpha}/Fe], are measured in SPs of model galaxies. The patterns found in the model for SPs with different origins help to explain the lack of a gradient of [{\alpha}/Fe] ratios in observed stars along the minor axis of the bulge.

We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.

The Large Synoptic Survey Telescope (LSST) will survey the southern sky from 2022--2032 with unprecedented detail. Since the observing strategy can lead to artifacts in the data, we investigate the effects of telescope-pointing offsets (called dithers) on the $r$-band coadded 5$\sigma$ depth yielded after the 10-year survey. We analyze this survey depth for several geometric patterns of dithers (e.g., random, hexagonal lattice, spiral) with amplitude as large as the radius of the LSST field-of-view, implemented on different timescales (per season, per night, per visit). Our results illustrate that per night and per visit dither assignments are more effective than per season. Also, we find that some dither geometries (e.g., hexagonal lattice) are particularly sensitive to the timescale on which the dithers are implemented, while others like random dithers perform well on all timescales. We then model the propagation of depth variations to artificial fluctuations in galaxy counts, which are a systematic for large-scale structure studies. We calculate the bias in galaxy counts caused by the observing strategy, accounting for photometric calibration uncertainties, dust extinction, and magnitude cuts; uncertainties in this bias limit our ability to account for structure induced by the observing strategy. We find that after 10 years of the LSST survey, the best dither strategies lead to uncertainties in this bias smaller than the minimum statistical floor for a galaxy catalog as deep as $r$$<$27.5. A few of these strategies bring the uncertainties close to the statistical floor for $r$$<$25.7 after only one year of survey.

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

Galaxy surveys targeting emission lines are characterising the evolution of star-forming galaxies, but there is still little theoretical progress in modelling their physical properties. We predict nebular emission from star-forming galaxies within a cosmological galaxy formation model. Emission lines are computed by combining the semi-analytical model \sag\ with the photoionisation code \mapp. We characterise the interstellar medium (ISM) of galaxies by relating the ionisation parameter of gas in galaxies to their cold gas metallicity, obtaining a reasonable agreement with the observed \ha, \oii, \oiii\ luminosity functions, and the the BPT diagram for local star-forming galaxies. The average ionisation parameter is found to increase towards low star-formation rates and high redshifts, consistent with recent observational results. The predicted link between different emission lines and their associated star-formation rates is studied by presenting scaling relations to relate them. Our model predicts that emission line galaxies have modest clustering bias, and thus reside in dark matter haloes of masses below $M_{\rm halo} \lesssim 10^{12} {[\rm h^{-1} M_{\odot}]}$. Finally, we exploit our modelling technique to predict galaxy number counts up to $z\sim 10$ by targeting far-infrared (FIR) emission lines detectable with submillimetre facilities

We present a fast and accurate method to select an optimal set of parameters in semi-analytic models of galaxy formation and evolution (SAMs). Our approach compares the results of a model against a set of observables applying a stochastic technique called Particle Swarm Optimization (PSO), a self-learning algorithm for localizing regions of maximum likelihood in multidimensional spaces that outperforms traditional sampling methods in terms of computational cost. We apply the PSO technique to the SAG semi-analytic model combined with merger trees extracted from a standard $\Lambda$CDM N-body simulation. The calibration is performed using a combination of observed galaxy properties as constraints, including the local stellar mass function and the black hole to bulge mass relation. We test the ability of the PSO algorithm to find the best set of free parameters of the model by comparing the results with those obtained using a MCMC exploration. Both methods find the same maximum likelihood region, however the PSO method requires one order of magnitude less evaluations. This new approach allows a fast estimation of the best-fitting parameter set in multidimensional spaces, providing a practical tool to test the consequences of including other astrophysical processes in SAMs.

We present a new technique for modeling submillimeter galaxies (SMGs): the "Count Matching" approach. Using lightcones drawn from a semi-analytic model of galaxy formation, we choose physical galaxy properties given by the model as proxies for their submillimeter luminosities, assuming a monotonic relationship. As recent interferometric observations of the Extended Chandra Deep Field South show that the brightest sources detected by single-dish telescopes are comprised by emission from multiple fainter sources, we assign the submillimeter fluxes so that the combined LABOCA plus bright-end ALMA observed number counts for this field are reproduced. After turning the model catalogs given by the proxies into submillimeter maps, we perform a source extraction to include the effects of the observational process on the recovered counts and galaxy properties. We find that for all proxies, there are lines of sight giving counts consistent with those derived from LABOCA observations, even for input sources with randomized positions in the simulated map. Comparing the recovered redshift, stellar mass and host halo mass distributions for model SMGs with observational data, we find that the best among the proposed proxies is that in which the submillimeter luminosity increases monotonically with the product between dust mass and SFR. This proxy naturally reproduces a positive trend between SFR and bolometric IR luminosity. The majority of components of blended sources are spatially unassociated.

We make use of a semi-analytical model of galaxy formation to investigate the origin of the observed correlation between [a/Fe] abundance ratios and stellar mass in elliptical galaxies. We implement a new galaxy-wide stellar initial mass function (Top Heavy Integrated Galaxy Initial Mass Function, TH-IGIMF) in the semi-analytic model SAG and evaluate its impact on the chemical evolution of galaxies. The SFR-dependence of the slope of the TH-IGIMF is found to be key to reproducing the correct [a/Fe]-stellar mass relation. Massive galaxies reach higher [a/Fe] abundance ratios because they are characterized by more top-heavy IMFs as a result of their higher SFR. As a consequence of our analysis, the value of the minimum embedded star cluster mass and of the slope of the embedded cluster mass function, which are free parameters involved in the TH-IGIMF theory, are found to be as low as 5 solar masses and 2, respectively. A mild downsizing trend is present for galaxies generated assuming either a universal IMF or a variable TH-IGIMF. We find that, regardless of galaxy mass, older galaxies (with formation redshifts > 2) are formed in shorter time-scales (< 2 Gyr), thus achieving larger [a/Fe] values. Hence, the time-scale of galaxy formation alone cannot explain the slope of the [a/Fe]-galaxy mass relation, but is responsible for the big dispersion of [a/Fe] abundance ratios at fixed stellar mass.We further test the hyphothesis of a TH-IGIMF in elliptical galaxies by looking into mass-to-light ratios, and luminosity functions. Models with a TH-IGIMF are also favoured by these constraints. In particular, mass-to-light ratios agree with observed values for massive galaxies while being overpredicted for less massive ones; this overprediction is present regardless of the IMF considered.

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 test and present the application of the full rescaling method by Angulo & White (2010) to change the cosmology of halo catalogues in numerical simulations for cosmological parameter search using semi-analytic galaxy properties. We show that a reduced form of the method can be applied in small simulations with box side of ~50/h Mpc. We perform statistical tests on the accuracy of the properties of rescaled individual haloes, and also on the rescaled population as a whole. We find that individual positions and velocities are recovered with almost no detectable biases. The dispersion in the recovered halo mass does not seem to depend on the resolution of the simulation. Regardless of the halo mass, the individual accretion histories, spin parameter evolution and fraction of mass in substructures are well recovered. The mass of rescaled haloes can be underestimated (overestimated) for negative (positive) variations of either sigma_8 or Omega_m, in a way that does not depend on the halo mass. Statistics of abundances and correlation functions of haloes show also small biases of <10 percent when moving away from the base simulation by up to 2 times the uncertainty in the WMAP7 cosmological parameters. The merger tree properties related to the final galaxy population in haloes also show small biases; the time since the last major merger, the assembly time-scale, and a time-scale related to the stellar ages show correlated biases which indicate that the spectral shapes of galaxies would only be affected by global age changes of ~150 Myr. We show some of these biases for different separations in the cosmological parameters with respect to the desired cosmology so that these can be used to estimate the expected accuracy of the resulting halo population. We also present a way to construct grids of simulations to provide stable accuracy across the Omega_m vs sigma_8 parameter space.

We investigate the origin of the color-magnitude relation (CMR) observed in cluster galaxies by using a combination of a cosmological N-body simulation of a cluster of galaxies and a semi-analytic model of galaxy formation. The departure of galaxies in the bright end of the CMR with respect to the trend defined by less luminous galaxies could be explained by the influence of minor mergers.

We investigate the development of the red sequence (RS) of cluster galaxies by using a semi-analytic model of galaxy formation. Results show good agreement between the general trend of the simulated RS and the observed relation in early-type galaxies. However, the most luminous galaxies ($M_V \lesssim -20$) depart from the linear fit to observed data, displaying almost constant colours. We analyze the dependence with redshift of the fraction of stellar mass contributed to each galaxy by different processes (i.e., quiescent star formation, disc instability and mergers), finding that the evolution of the bright end, since $z\approx 2$, is mainly driven by minor and major dry mergers. Since the most luminous galaxies have a narrow spread in ages ($1.0\times 10^{10}$ yr $<t<1.2\times 10^{10}$ yr), their metallicities are the main factor that affects their colours. Galaxies in the bright end reach an upper limit in metallicity as a result of the competition of the mass of stars and metals provided by the star formation within the galaxies and by the accretion of merging satellites. Star formation activity in massive galaxies (M_\star \gtrsim 10^{10} M_{\odot}$) contribute with stellar components of high metallicity, but this fraction of stellar mass is negligible. Mergers contribute with a larger fraction of stellar mass ($\approx 10-20$ per cent), but the metallicity of the accreted satellites is lower by $\approx 0.2$ dex than the mean metallicity of galaxies they merge with. The effect of dry mergers is to increase the mass of galaxies in the bright end, without significantly altering their metallicities, and hence,their colours, giving rise to the break in the RS. These results are found for clusters with different virial masses, supporting the idea of the universality of the CMR in agreement with observational results.

We investigate the development of the colour-magnitude re- lation (CMR) of cluster galaxies. This study is carried out using a semi- analytic model of galaxy formation and evolution coupled to a sample of simulated galaxy clusters of different masses, reinforcing the conclusions reached by Jim\'enez et al. (2009). We compare both simulated and obeserved CMRs in different colour-magnitude planes, finding a very good agreement in all cases. This indicates that model parameters are correctly tuned, giving accurate values of the main properties of galaxies for further use in our study. In the present work, we perform a statistical analysis of the relative contribution to the stellar mass and metallicity of galaxies along the CMR by the different processes involved in their formation and evolution (i.e. quiescent star formation, disc instability events and galaxy mergers). Our results show that a mix of minor and major dry mergers at low redshifts is relevant in the evolution of the most luminous galaxies in the CMR. These processes contribute with low metallicity stars to the remnant galaxies, thus increasing the galaxy masses without significantly altering their colours. These results are found for all simulated clusters, supporting the idea of the universality of the CMR in agreement with observational results.

We investigate the origin of the color-magnitude relation (CMR) observed in cluster galaxies by using a combination of cosmological N-body/SPH simulations of galaxy clusters, and a semi-analaytic model of galaxy formation (Lagos, Cora & Padilla 2008). Simulated results are compared with the photometric properties of early-type galaxies in the Antlia cluster (Smith Castelli et al. 2008). The good agreement obtained between observations and simulations allows us to use the information provided by the model for unveiling the physical processes that yield the tigh observed CMR.

We investigate the origin of the colour-magnitude relation (CMR) observed in cluster galaxies by using a combination of a cosmological N-body simulation of a cluster of galaxies and a semi-analytic model of galaxy formation. The departure of galaxies in the bright end of the CMR with respect to the trend denoted by less luminous galaxies could be explained by the influence of minor mergers

Using a hybrid method which combines non-radiative hydrodynamical simulations with a semi-analytic model of galaxy formation, we determine the ram pressure as a function of halocentric distance experienced by galaxies in haloes with virial masses 12.5 <= log (M_200 h/M_Sun) < 15.35, for redshifts 0 <= z <= 3. The ram pressure is calculated with a self-consistent method which uses the simulation gas particles to obtain the properties of the intergalactic medium. The ram pressure profiles obtained can be well described by beta profile models, with parameters that depend on redshift and halo virial mass in a simple fashion. The fitting formulae provided here will prove useful to include ram pressure effects into semi-analytic models based on methods which lack gas physics, such as dark matter-only simulations or the Press-Schechter formalism.

We investigate the origin of the colour-magnitude relation (CMR) followed by early-type cluster galaxies by using a combination of cosmological N-body simulations of cluster of galaxies and a semi-analytic model of galaxy formation (Lagos, Cora & Padilla 2008). Results show good agreement between the general trend of the simulated and observed CMR. However, in many clusters, the most luminous galaxies depart from the linear fit to observed data displaying almost constant colours. With the aim of understanding this behaviour, we analyze the dependence with redshift of the stellar mass contributed to each galaxy by different processes, i.e., quiescent star formation, and starburst during major/minor and wet/dry mergers, and disk instability events. The evolution of the metallicity of the stellar component, contributed by each of these processes, is also investigated. We find that the major contribution of stellar mass at low redshift is due to minor dry merger events, being the metallicity of the stellar mass accreted during this process quite low. Thus, minor dry merger events seem to increase the mass of the more luminous galaxies without changing their colours.

We determine the intrinsic shapes and orientations of 27,450 type I and II active galactic nucleus (AGN) galaxies in the spectroscopic sample of the SDSS DR7, by studying the distribution of projected axis ratios of AGN hosts. Our aim is to study possible alignments between the AGN and host galaxy systems (e.g. the accretion disc and the galaxy angular momentum) and the effect of dust obscuration geometry on the AGN type. We define control samples of non-AGN galaxies that mimic the morphology, colour, luminosity and concentration distributions of the AGN population, taking into account the effects of dust extinction and reddening. Assuming that AGN galaxies have the same underlying three-dimensional shape distribution as their corresponding control samples, we find that the spiral and elliptical type I AGN populations are strongly skewed toward face-on galaxies, while ellipticals and spirals type II AGN are skewed toward edge-on orientations. These findings rule out random orientations for AGN hosts at high confidence for type I spirals (delta chi^2~230) and type II ellipticals (delta chi^2~15), while the signal for type I ellipticals and type II spirals is weaker (delta chi^2~3 and delta chi^2~6, respectively). We obtain a much stronger tendency for the type II spirals to be edge-on when just high [OIII] equivalent width (EW) AGN are considered, suggesting that >20% of low [OIII] EW edge-on type II AGN may be missing from the optical sample. Galactic dust absorption of the broad-line region alone cannot explain the observed inclination angle and projected axis ratio distributions of type I and II Seyfert types, implying that obscuration by a small-scale circumnuclear torus is necessary. These results favour a scenario in which the angular momentum of the material which feeds the black hole retains a memory of its original gas source at least to some small, non-negligible degree.

We develop a new numerical approach to describe the action of ram pressure stripping (RPS) within a semi-analytic model of galaxy formation and evolution which works in combination with non-radiative hydrodynamical simulations of galaxy clusters. The new feature in our method is the use of the gas particles to obtain the kinematical and thermodynamical properties of the intragroup and intracluster medium (ICM). This allows a self-consistent estimation of the RPS experienced by satellite galaxies. We find that the ram pressure (RP) in the central regions of clusters increases approximately one order of magnitude between z = 1 and 0, consistent with the increase in the density of the ICM. The mean RP experienced by galaxies within the virial radius increases with decreasing redshift. In clusters with virial masses M_vir ~10^15 h^-1 M_Sun, over 50 per cent of satellite galaxies have experienced RP ~10^(-11) h^2 dyn cm^-2 or higher for z <= 0.5. In smaller clusters (M_vir ~10^14 h^-1 M_Sun) the mean RP are approximately one order of magnitude lower at all redshifts. RPS has a strong effect on the cold gas content of galaxies for all cluster masses. At z = 0, over 70 per cent of satellite galaxies within the virial radius are completely depleted of cold gas. For the more massive clusters the fraction of depleted galaxies is already established at z ~ 1, whereas for the smaller clusters this fraction increases appreciably between z = 1 and 0. This indicates that the rate at which the cold gas is stripped depends on the virial mass of the host cluster. Compared to our new approach, the use of an analytic profile to describe the ICM results in an overestimation of the RP larger than 50 per cent for z >= 0.5.

We study the properties of building blocks (BBs, i.e. accreted satellites) and surviving satellites of present-day galaxies using the SAG semi-analytic model of galaxy formation in the context of a concordance Lambda Cold Dark Matter (LCDM) cosmology. We consider large numbers of DM halo merger trees spanning a wide range of masses (~1x10^10 - 2.14x10^15 Msun). We find higher metallicities for BBs with respect to surviving satellites, an effect produced by the same processes behind the build-up of the mass-metallicity relation. We prove that these metallicity differences arise from the higher peak height in the density fluctuation field occupied by BBs and central galaxies which have collapsed into a single object earlier than surviving satellites. BBs start to form stars earlier, during the peak of the merger activity in LCDM, and build-up half of their final stellar mass (measured at the moment of disruption) up to four times faster than surviving satellites. Surviving satellites keep increasing their stellar masses rather quiescently down to z~1. The difference between the metallicities of satellites, BBs and central galaxies depends on the host DM halo mass, in a way that can be used as a further test for the concordance cosmology.

We use a combination of a cosmological N-body simulation of the concordance Lambda cold dark matter (LCDM) paradigm and a semi-analytic model of galaxy formation to investigate the spin development of central supermassive black holes (BHs) and its relation to the BH host galaxy properties. In order to compute BH spins, we use the alpha-model of Shakura & Sunyaev and consider the King et al. warped disc alignment criterion. The orientation of the accretion disc is inferred from the angular momentum of the source of accreted material, which bears a close relationship to the large-scale structure in the simulation. We find that the final BH spin depends almost exclusively on the accretion history and only weakly on the warped disc alignment. The main mechanisms of BH spin-up are found to be gas cooling processes and disc instabilities, a result that is only partially compatible with Monte-Carlo models where the main spin-up mechanisms are major mergers and disc instabilities; the latter results are reproduced when implementing randomly oriented accretion discs in our model. Regarding the BH population, we find that more massive BHs, which are hosted by massive ellipticals, have higher spin values than less-massive BHs, hosted by spiral galaxies. We analyse whether gas accretion rates and BH spins can be used as tracers of the radio loudness of active galactic nuclei (AGN). We find that the current observational indications of an increasing trend of radio-loud AGN fractions with stellar and BH mass can be easily obtained when placing lower limits on the BH spin, with a minimum influence from limits on the accretion rates; a model with random accretion disc orientations is unable to reproduce this trend. (ABRIDGED)

We model the high-energy emission that results from the interaction of relativistic particles with photons and matter in the cluster of galaxies Abell 3376. The presence of relativistic particles is inferred from the recently found radio relics in this cluster, being one of the most prominent examples of double opposite, giant ringlike radio structures. Assuming that diffusive shock acceleration takes place in the cluster regions where radio relics are observed, we calculate the spectral energy distribution resulting from the most relevant non-thermal processes, which are synchrotron radiation, inverse Compton scattering, relativistic Bremsstrahlung, and inelastic proton-proton collisions. In the context of our model, the major radiative component at high energies is inverse Compton scattering, which could reach luminosities L ~ 9x10^{41} erg/s in the energy range between ~ 1 MeV and 10 TeV. Hadronic interactions would yield a minor contribution to the overall non-thermal emission, but would dominate at ultra-high energies. The cluster Abell 3376 might be detectable at gamma-rays by HESS, GLAST satellite and future planned Cherenkov arrays.

We study the effects of Active Galactic Nuclei (AGN) feedback on the formation and evolution of galaxies in a semi-analytic model of galaxy formation. This model is an improved version of the one described by Cora (2006), which now considers the growth of black holes (BHs) as driven by (i) gas accretion during merger-driven starbursts and mergers with other BHs, (ii) accretion during starbursts triggered by disc instabilities, and (iii) accretion of gas cooled from quasi-hydrostatic hot gas haloes. It is assumed that feedback from AGN operates in the later case. The model has been calibrated in order to reproduce observational correlations between BH mass and mass, velocity dispersion, and absolute magnitudes of the galaxy bulge. AGN feedback has a strong impact on reducing or even suppressing gas cooling, an effect that becomes important at lower redshifts. This phenomenon helps to reproduce the observed galaxy luminosity function (LF) in the optical and near IR bands at z=0, and the cosmic star formation rate and stellar mass functions over a wide redshift range (0<z<5). It also allows to have a population of massive galaxies already in place at z>1, which are mostly early-type and have older and redder stellar populations than lower mass galaxies, reproducing the observed bimodality in the galaxy colour distribution, and the morphological fractions. The evolution of the optical QSO LF is also reproduced, provided that the presence of a significant fraction of obscured QSOs is assumed. We explore the effects of AGN feedback during starbursts and new recent prescriptions for dynamical friction time-scales. (ABRIDGED)

We present a study on the origin of the metallicity evolution of the intra-cluster medium (ICM) by applying a semi-analytic model of galaxy formation to N-body/SPH (smoothed particle hydrodynamic) non-radiative numerical simulations of clusters of galaxies. The semi-analytic model includes gas cooling, star formation, supernovae feedback and metal enrichment, and is linked to the diffuse gas of the underlying simulations so that the chemical properties of gas particles are dynamically and consistently generated from stars in the galaxies. This hybrid model let us have information on the spatial distribution of metals in the ICM. The results obtained for a set of clusters with virial masses of ~1.5*10^15 h^{-1} M_sun contribute to the theoretical interpretation of recent observational X-ray data, which indicate a decrease of the average iron content of the intra-cluster gas with increasing redshift. We find that this evolution arises mainly as a result of a progressive increase of the iron abundance within ~0.15 R_vir. The clusters have been considerably enriched by z~1 with very low contribution from recent star formation. Low entropy gas that has been enriched at high redshift sinks to the cluster centre contributing to the evolution of the metallicity profiles.

We develop a model for the metal enrichment of the intracluster medium (ICM) that combines a cosmological non-radiative hydrodynamical N-Body/SPH simulation of a cluster of galaxies, and a semi-analytic model of galaxy formation. The novel feature of our hybrid model is that the chemical properties of the diffuse gas in the underlying simulation are dynamically and consistently generated from stars in the galaxies. We follow the production of several chemical elements, provided by low- and intermediate-mass stars, core collapse and type Ia supernovae. We analyse the spatial distribution of metals in the ICM, investigate the way in which the chemical enrichment proceeds, and use iron emissivity as a tracer of gas motions. Fe and O radial abundance profiles are enhanced in the inner 100 h^-1 kpc in the last Gyr because of the convergence of enriched gas clumps to the cluster centre. Our results support a scenario in which part of the central intracluster gas comes from gas clumps that, in the redshift range of z~0.2 to ~0.5, have been enriched to solar values and are at large distances from the cluster centre (from ~1 to ~6 h^-1 Mpc) moving at very high velocities (from ~1300 to ~2500 km s^-1). The turbulent gas motions within the cluster, originated in the inhomogeneous gas infall during the cluster assembly, are manifested in emission-weighted velocity maps as gradients that can be as large as ~1000 km s^-1 over distances of a few hundred kpc. Gradients of this magnitude are also seen in velocity distributions along sightlines through the cluster centre. Doppler shifting and broadening suffered by the Fe K 6.7 keV emission line along such sightlines could be used to probe these gas large-scale motions when they are produced within an area characterised by high iron line emissivity.

We carried out a comprehensive analysis of the chemical properties of the interstellar medium (ISM) and the stellar population (SP) of current normal galaxies and their progenitors in a hierarchical clustering scenario. We compared the results with observations of Damped Lyman-$\alpha$ systems (DLAs) under the hypothesis that, at least, part of the observed DLAs could originate in the building blocks of today normal galaxies. We used a hydrodynamical cosmological code which includes star formation and chemical enrichment. Galaxy-like objects are identified at $z=0$ and then followed back in time. Random line-of-sights (LOS) are drawn through these structures in order to mimic Damped Lyman $\alpha$ systems. We then analysed the chemical properties of the ISM and SP along the LOS. We found that the progenitors of current galaxies in the field with mean $L <0.5 L^* $ and virial circular velocity of $100-250 {\rm km/sec}$ could be the associated DLA galaxies. For these systems we detected a trend for $<L/L^*>$ to increase with redshift.(Abridged)

The combined enrichment of Supernovae II and I in a hierarchical clustering scenario could produce regions with low N content respect to $\alpha$-elements consistent with observed values measured in Damped Ly-$\alpha$ (DLAs). We have studied the formation of DLAs in a hierarchical clustering scenario under the hypothesis that the building blocks of current field galaxies could be part of the structures mapped by DLAs. In our models the effects of the non-linear evolution of the structure (which produces bursty star formation histories, gas infall, etc.) and the contributions of SNIa and SNII are found to be responsible of producing these N regions with respect to the $\alpha$-elements. Although SNIa are not main production sites for Si or O, because of the particular timing Consistently, we found the simulated low nitrogen DLAs to have sub-solar [Fe/H]. We show that low nitrogen DLAs have experienced important star formation activity in the past with higher efficiency than normal DLAs. Our chemical model suggests that SNIa play a relevant role in the determination of the abundance pattern of DLA and, that the observed low nitrogen DLA frequency could be explained taking into account the time-delay of $\approx $ 0.5 Gyr introduced by these supernova to release metals.

We use cosmological hydrodynamical simulations including star formation and metal enrichment to study the evolution of the chemical properties of galaxy-like objects at high redshift in the range $0.25<z< 2.35$ in a hierarchical clustering scenario. As the galactic objects are assembled we find that their gaseous components exhibit neutral Hydrogen column densities with abundances and scatter comparable to those observed in damped Lyman-$\alpha$ systems (DLAs).The unweighted mean of abundance ratios and least square linear regressions through the simulated DLAs yield intrinsic metallicity evolution for the [Zn/H] and [Fe/H], consistent with results obtained from similar analysis of available observations. Our model statistically reproduces the mild evolution detected in the metallicity of the neutral hydrogen content of the Universe, given by mass-weighted means,if observational constraints are considered (as suggested by Boiss\'ee et al. 1998). For the $\alpha$-elements in the simulated DLAs, we find neither enhancement nor dependence on metallicity. Our results support the hypotheses that DLAs trace a variety of galactic objects with different formation histories and that both SNI and SNII are contributing to the chemical enrichment of the gas component at least since $z \approx 2$. This study indicates that DLAs could be understood as the building blocks that merged to form today normal galaxies within a hierarchical clustering scenario.

We develop a model to implement metal enrichment in a cosmological context based on the hydrodynamical AP3MSPH code described by Tissera, Lambas and Abadi (1997).The star formation model is based on the Schmidt law and has been modified in order to describe the transformation of gas into stars in more detail. The enrichment of the interstellar medium due to supernovae I and II explosions is taken into account by assuming a Salpeter Initial Mass Function and different nucleosynthesis models.The different chemical elements are mixed within the gaseous medium according to the Smooth Particle Hydrodynamics technique. We present tests of the code that assess the effects of resolution and model parameters on the results. scenario and we present results of the analysis of the star formation and chemical properties of the interstellar medium and stellar population of the simulated galactic objects. We show that these systems reproduce abundance ratios for primary and secondary elements of the interstellar medium, and the correlation between the (O/H) abundance and the gas fraction of galaxies. The numerical simulations performed provide a detailed description of the chemical properties of galactic objects formed in hierarchical clustering scenarios and proved to be useful tools to deepen our understanding of galaxy formation and evolution.

We present preliminary results on the effects of mergers on the chemical properties of galactic objects in hierarchical clustering scenarios. We adopt a hydrodynamical chemical code that allows to describe the coupled evolution of dark matter and baryons within a cosmological context. We found that disk-like and spheroid-like objects have distinctive metallicity patterns that may be the result of different evolution.

We report first results of an implementation of a chemical model in a cosmological code, based on the Smoothed Particle Hydrodynamics (SPH) technique. We show that chemical SPH simulations are a promising tool to provide clues for the understanding of the chemical properties of galaxies in relation to their formation and evolution in a cosmological framework.