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 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.

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.

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 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.

For several decades now, open clusters have been used to study the structure and chemical evolution of the disk of our Galaxy. Due to the fact that their ages and metallicities can be determined with relatively good precision, and since they can be observed even at great distances, they are excellent tracers of the variations in the abundance of heavy chemical elements with age and position in the Galactic disk. In the present work we analyze the star formation history and the chemical evolution of the disk of the Galaxy using numerical simulations of the dynamical evolution of the system of open clusters in the Milky Way. Starting from hypotheses on the history of cluster formation and the chemical enrichment of the disk, we model the present properties of the Galactic open cluster system. The comparison of these models with the observations allows us to examine the validity of the assumed hypotheses and to improve our knowledge about the initial conditions of the chemical evolution of the Galactic disk.