Simpson's paradox, or Yule-Simpson effect, arises when a trend appears in different subsets of data but disappears or reverses when these subsets are combined. We describe here seven cases of this phenomenon for chemo-kinematical relations believed to constrain the Milky Way disk formation and evolution. We show that interpreting trends in relations, such as the radial and vertical chemical abundance gradients, the age-metallicity relation, and the metallicity-rotational velocity relation (MVR), can lead to conflicting conclusions about the Galaxy past if analyses marginalize over stellar age and/or birth radius. It is demonstrated that the MVR in RAVE giants is consistent with being always strongly negative, when narrow bins of [Mg/Fe] are considered. This is directly related to the negative radial metallicity gradients of stars grouped by common age (mono-age populations) due to the inside out disk formation. The effect of the asymmetric drift can then give rise to a positive MVR trend in high-[alpha/Fe] stars, with a slope dependent on a given survey's selection function and observational uncertainties. We also study the variation of lithium abundance, A(Li), with [Fe/H] of AMBRE:HARPS dwarfs. A strong reversal in the positive A(Li)-[Fe/H] trend of the total sample is found for mono-age populations, flattening for younger groups of stars. Dissecting by birth radius shows strengthening in the positive A(Li)-[Fe/H] trend, shifting to higher [Fe/H] with decreasing birth radius; these observational results suggest new constraints on chemical evolution models. This work highlights the necessity for precise age estimates for large stellar samples covering wide spatial regions.

One of the fundamental assumptions of chemical evolution models (CEMs) of the Milky Way (MW) and other spirals is that higher gas accretion rates are expected in the past, and in the inner regions of the Galaxy. This leads to the so-called `inside-out disc formation scenario'. Yet, these are probably the most unconstrained inputs of such models. In the present paper, we aim at investigating these main assumptions by studying how gas is accreted in four simulated MW-like galaxies assembled within the $\Lambda$CDM scenario. The galaxies were obtained using two different simulation techniques, cosmological setups and initial conditions. Two of them are MW candidates corresponding to the chemodynamical model of Minchev et al. (2013, 2014) (known as MCM) and the Local Group cosmological simulation of Nuza et al. (2014). We investigate vertical and radial gas accretion on to galaxy discs as a function of cosmic time and disc radius. We find that accretion in the MW-like galaxies seem to happen in two distinct phases, namely: an early, more violent period; followed by a subsequent, slowly declining phase. Our simulations seem to give support to the assumption that the amount of gas incorporated into the MW disc exponentially decreases with time, leading to current net accretion rates of $0.6-1\,$M$_\odot\,$yr$^{-1}$. In particular, accretion timescales on to the simulated thin-disc-like structures are within $\sim5-7\,$Gyr, consistent with expectations from CEMs. Moreover, our simulated MW discs are assembled from the inside-out with gas in the inner disc regions accreted in shorter timescales than in external ones, in qualitative agreement with CEMs of the Galaxy. However, this type of growth is not general to all galaxies and it is intimately linked to their particular merger and gas accretion history.

We present an update to the chemical enrichment component of the smoothed-particle hydrodynamics model for galaxy formation presented in Scannapieco et al. (2005) in order to address the needs of modelling galactic chemical evolution in realistic cosmological environments. Attribution of the galaxy-scale abundance patterns to individual enrichment mechanisms such as the winds from asymptotic giant branch (AGB) stars or the presence of a prompt fraction of Type Ia supernovae is complicated by the interaction between them and gas cooling, subsequent star formation and gas ejection. In this work we address the resulting degeneracies by extending our implementation to a suite of mechanisms that encompasses different IMFs, models for yields from the aforementioned stars, models for the prompt component of the delay-time-distribution (DTDs) for Type Ia SNe and metallicity-dependent gas cooling rates, and then applying these to both isolated initial conditions and cosmological hydrodynamical zoom simulations. We find DTDs with a large prompt fraction (such as the bimodal and power-law models) have, at z=0, similar abundance patterns compared to the low-prompt component time distributions (uniform or wide Gaussian models). However, some differences appear, such as the former having systematically higher [X/Fe] ratios and narrower [O/Fe] distributions compared to the latter, and a distinct evolution of the [Fe/H] abundance.

In this work we present a set of synthetic observations that mimic the properties of the Integral Field Spectroscopy (IFS) survey CALIFA, generated using radiative transfer techniques applied to hydrodynamical simulations of galaxies in a cosmological context. The simulated spatially-resolved spectra include stellar and nebular emission, kinematic broadening of the lines, and dust extinction and scattering. The results of the radiative transfer simulations have been post-processed to reproduce the main properties of the CALIFA V500 and V1200 observational setups. The data has been further formatted to mimic the CALIFA survey in terms of field of view size, spectral range and sampling. We have included the effect of the spatial and spectral Point Spread Functions affecting CALIFA observations, and added detector noise after characterizing it on a sample of 20 galaxies. The simulated datacubes are suited to be analyzed by the same algorithms used on real IFS data. In order to provide a benchmark to compare the results obtained applying IFS observational techniques to our synthetic datacubes, and test the calibration and accuracy of the analysis tools, we have computed the spatially-resolved properties of the simulations. Hence, we provide maps derived directly from the hydrodynamical snapshots or the noiseless spectra, in a way that is consistent with the values recovered by the observational analysis algorithms. Both the synthetic observations and the product datacubes are public and can be found in the collaboration website http://astro.ft.uam.es/selgifs/data_challenge/.

In chemodynamical evolution models it is usually assumed that the Milky Way galaxy forms from the inside-out implying that gas inflows onto the disk decrease with galactocentric distance. Similarly, to reproduce differences between chemical abundances of the thick disk and bulge with respect to those of the thin disk, higher accretion fluxes at early times are postulated. By using a suite of Milky Way-like galaxies extracted from cosmological simulations, we investigate the accretion of gas on the simulated stellar disks during their whole evolution. In general, we find that the picture outlined above holds, although the detailed behavior depends on the assembly history of the Galaxy and the complexities inherent to the physics of galaxy formation.

We study the effects of applying observational techniques to derive the properties of simulated galaxies, with the aim of making an unbiased comparison between observations and simulations. For our study, we used fifteen galaxies simulated in a cosmological context using three different feedback and chemical enrichment models, and compared their z=0 properties with data from the Sloan Digital Sky Survey (SDSS). We show that the physical properties obtained directly from the simulations without post-processing can be very different to those obtained mimicking observational techniques. In order to provide simulators a way to reliably compare their galaxies with SDSS data, for each physical property that we studied - colours, magnitudes, gas and stellar metallicities, mean stellar ages and star formation rates - we give scaling relations that can be easily applied to the values extracted from the simulations. These scalings have in general a high correlation, except for the galaxy mean stellar ages and gas oxygen metallicities. Our simulated galaxies are photometrically similar to galaxies in the blue cloud/green valley, but in general they appear older, passive and with lower metal content compared to most of the spirals in SDSS. As a careful assessment of the agreement/disagreement with observations is the primary test of the baryonic physics implemented in hydrodynamical codes, our study shows that considering the observational biases in the derivation of the galaxies' properties is of fundamental importance to decide on the failure/success of a galaxy formation model.

In this work, we investigate the contribution of dust scattering to the diffuse H-alpha emission observed in nearby galaxies. As initial conditions for the spatial distribution of HII regions, gas, and dust, we take three Milky Way-like galaxies from state-of-the-art cosmological hydrodynamical simulations that implement different prescriptions for star formation, feedback, and chemical enrichment. Radiative transfer has been solved a posteriori, using the publicly-available Monte Carlo code Sunrise to take into account dust absorption and scattering of the H-alpha photons, originating exclusively from the HII regions. No contribution from recombinations in the diffuse ionized gas (DIG) component is explicitly or implicitly included in our model. Our main result is that the flux arising from scattered light is of the order of 1-2 per cent of the H-alpha flux coming directly from the HII regions. Building upon previous studies, we conclude that the DIG contributes lass than 50 per cent of the total H-alpha emission.

We study the sources of biases and systematics in the derivation of galaxy properties of observational studies, focusing on stellar masses, star formation rates, gas/stellar metallicities, stellar ages and magnitudes/colors. We use hydrodynamical cosmological simulations of galaxy formation, for which the real quantities are known, and apply observational techniques to derive the observables. We also make an analysis of biases that are relevant for a proper comparison between simulations and observations. For our study, we post-process the simulation outputs to calculate the galaxies' spectral energy distributions (SEDs) using Stellar Population Synthesis models and also generating the fully-consistent far UV-submillimeter wavelength SEDs with the radiative transfer code SUNRISE. We compared the direct results of simulations with the observationally-derived quantities obtained in various ways, and found that systematic differences in all studied galaxy properties appear, which are caused by: (1) purely observational biases, (2) the use of mass-weighted/luminosity-weighted quantities, with preferential sampling of more massive/luminous regions, (3) the different ways to construct the template of models when a fit to the spectra is performed, and (4) variations due to the use of different calibrations, most notably in the cases of the gas metallicities and star formation rates. Our results show that large differences can appear depending on the technique used to derive galaxy properties. Understanding these differences is of primary importance both for simulators, to allow a better judgement on similarities/differences with observations, and for observers, to allow a proper interpretation of the data.

Recent spectroscopic observations in the Milky Way suggest that the chemically defined thick disk (stars with high [alpha/Fe] ratios and thus old) has a significantly smaller scale-length than the thin disk. This is in apparent contradiction with observations of external edge-on galaxies, where the thin and thick components have comparable scale-lengths. Moreover, while observed disks do not flare (scale-height does not increase with radius), numerical simulations suggest that disk flaring is unavoidable, resulting from both environmental effects and secular evolution. Here we address these problems by studying two different suites of simulated galactic disks formed in the cosmological context. We show that the scale-heights of coeval populations always increase with radius. However, the total population can be decomposed morphologically into thin and thick disks, which do not flare. We relate this to the disk inside-out formation, where younger populations have increasingly larger scale-lengths and flare at progressively larger radii. In this new picture, thick disks are composed of the imbedded flares of mono-age stellar populations. Assuming that disks form inside out, we predict that morphologically defined thick disks must show a decrease in age (or [alpha/Fe] ratios) with radius and that coeval populations should always flare. This also explains the observed inversion in the metallicity and [alpha/Fe] gradients for stars away from the disk midplane in the Milky Way. The results of this work are directly linked to, and can be seen as evidence of, inside-out disk growth.

We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way (MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. We use a state-of-the-art hydrodynamical code which includes star formation, feedback and chemical enrichment to carry out our study. We run two simulations, where we include or neglect the effects of radiation pressure from stars, to investigate the impact of this process on the morphologies and star formation rates of the simulated galaxies. We find that the simulated M31 and MW have different formation histories, even though both inhabit, at z=0, the same environment. These differences directly translate into and explain variations in their star formation rates, in-situ fractions and final morphologies. The M31 candidate has an active merger history, as a result of which its stellar disc is unable to survive unaffected until the present time. In contrast, the MW candidate has a smoother history with no major mergers at late times, and forms a disc that grows steadily; at z=0 the simulated MW has an extended, rotationally-supported disc which is dominant over the bulge. Our two feedback implementations predict similar evolution of the galaxies and their discs, although some variations are detected, the most important of which is the formation time of the discs: in the model with weaker/stronger feedback the discs form earlier/later. In summary, by comparing the formation histories of the two galaxies, we conclude that the particular merger/accretion history of a galaxy rather than its environment at the LG-scales is the main driver of the formation and subsequent growth or destruction of galaxy discs.

In this letter we present, for the first time, a study of star formation rate, gas fraction and galaxy morphology of a constrained simulation of the Milky Way (MW) and Andromeda (M31) galaxies, compared to other MW-mass galaxies. By combining with unconstrained simulations we cover a sufficient volume to compare these galaxies environmental densities ranging from the field to that of the Local Group (LG). This is particularly relevant as it has been shown that, quite generally, galaxy properties depend intimately upon their environment, most prominently when galaxies in clusters are compared to those in the field. For galaxies in loose groups such as the LG, however, environmental effects have been less clear. We consider the galaxy's environmental density in spheres of 1200 kpc (comoving) and find that whilst environment does not appear to directly affect morphology, there is a positive trend with star formation rates. This enhancement in star formation occurs systematically for galaxies in higher density environments, regardless whether they are part of the LG or in filaments. Our simulations suggest that the richer environment at Mpc-scales may help replenish the star-forming gas, allowing higher specific star formation rates in galaxies such as the MW.

We analyse the density profiles of the stellar halo populations in eight Milky-Way mass galaxies, simulated within the $\Lambda$-Cold Dark Matter scenario. We find that accreted stars can be well-fitted by an Einasto profile, as well as any subsample defined according to metallicity. We detect a clear correlation between the Einasto fitting parameters of the low-metallicity stellar populations and those of the dark matter haloes. The correlations for stars with [Fe/H]$<-3$ allow us to predict the shape of the dark matter profiles within residuals of $\sim 10 $ per cent, in case the contribution from in situ stars remains small. Using Einasto parameters estimated for the stellar halo of the Milky Way and assuming the later formed with significant contributions from accreted low-mass satellite, our simulations predict $\alpha \sim 0.15 $ and $r_2 \sim 15$ kpc for its dark matter profile. These values, combined with observed estimations of the local dark matter density, yield an enclosed dark matter mass at $\sim 8$ kpc in the range $3.9 - 6.7 \times 10^{10}$ M$_{\odot}$, in agreement with recent observational results. These findings suggest that low-metallicity stellar haloes could store relevant information on the DM haloes. Forthcoming observations would help us to further constrain our models and predictions.

We study the gas distribution in the Milky Way and Andromeda using a constrained cosmological simulation of the Local Group (LG) within the context of the CLUES (Constrained Local UniversE Simulations) project. We analyse the properties of gas in the simulated galaxies at $z=0$ for three different phases: `cold', `hot' and HI, and compare our results with observations. The amount of material in the hot halo ($M_{hot}\approx 4-5\times10^{10}\,$M$_{\odot}$), and the cold ($M_{cold}(r\lesssim10\,$kpc$)\approx10^{8}\,$M$_{\odot}$) and HI ($M_{HI}(r\lesssim50\,$kpc$)\approx 3-4\times10^8\,$M$_{\odot}$) components display a reasonable agreement with observations. We also compute the accretion/ejection rates together with the HI (radial and all-sky) covering fractions. The integrated HI accretion rate within $r=50\,$kpc gives $\sim$$0.2-0.3\,$M$_{\odot}\,$yr$^{-1}$, i.e. close to that obtained from high-velocity clouds in the Milky Way. We find that the global accretion rate is dominated by hot material, although ionized gas with $T\lesssim10^5\,$K can contribute significantly too. The $net$ accretion rates of $all$ material at the virial radii are $6-8\,$M$_{\odot}\,$yr$^{-1}$. At $z=0$, we find a significant gas excess between the two galaxies, as compared to any other direction, resulting from the overlap of their gaseous haloes. In our simulation, the gas excess first occurs at $z\sim1$, as a consequence of the kinematical evolution of the LG.

We present a comprehensive study of the chemical properties of the stellar haloes of Milky-Way mass galaxies, analysing the transition between the inner to the outer haloes. We find the transition radius between the relative dominance of the inner-halo and outer-halo stellar populations to be ~15-20 kpc for most of our haloes, similar to that inferred for the Milky Way from recent observations. While the number density of stars in the simulated inner-halo populations decreases rapidly with distance, the outer-halo populations contribute about 20-40 per cent in the fiducial solar neighborhood, in particular at the lowest metallicities. We have determined [Fe/H] profiles for our simulated haloes; they exhibit flat or mild gradients, in the range [-0.002, -0.01 ] dex/kpc. The metallicity distribution functions exhibit different features, reflecting the different assembly history of the individual stellar haloes. We find that stellar haloes formed with larger contributions from massive subgalactic systems have steeper metallicity gradients. Very metal-poor stars are mainly contributed to the halo systems by lower-mass satellites. There is a clear trend among the predicted metallicity distribution functions that a higher fraction of low-metallicity stars are found with increasing radius. These properties are consistent with the range of behaviours observed for stellar haloes of nearby galaxies.

The velocity dispersions of stars near the Sun are known to increase with stellar age, but age can be difficult to determine so a proxy like the abundance of alpha elements (e.g., Mg) with respect to iron, [alpha/Fe], is used. Here we report an unexpected behavior found in the velocity dispersion of a sample of giant stars from the RAdial Velocity Experiment (RAVE) survey with high quality chemical and kinematical information, in that it decreases strongly for stars with [Mg/Fe] > 0.4 dex (i.e., those that formed in the first Gyr of the Galaxy's life). These findings can be explained by perturbations from massive mergers in the early Universe, which have affected more strongly the outer parts of the disc, and the subsequent radial migration of stars with cooler kinematics from the inner disc. Similar reversed trends in velocity dispersion are also found for different metallicity subpopulations. Our results suggest that the Milky Way disc merger history can be recovered by relating the observed chemo-kinematic relations to the properties of past merger events.

We construct new estimates on the Galactic escape speed at various Galactocentric radii using the latest data release of the Radial Velocity Experiment (RAVE DR4). Compared to previous studies we have a database larger by a factor of 10 as well as reliable distance estimates for almost all stars. Our analysis is based on the statistical analysis of a rigorously selected sample of 90 high-velocity halo stars from RAVE and a previously published data set. We calibrate and extensively test our method using a suite of cosmological simulations of the formation of Milky Way-sized galaxies. Our best estimate of the local Galactic escape speed, which we define as the minimum speed required to reach three virial radii $R_{340}$, is $533^{+54}_{-41}$ km/s (90% confidence) with an additional 5% systematic uncertainty, where $R_{340}$ is the Galactocentric radius encompassing a mean over-density of 340 times the critical density for closure in the Universe. From the escape speed we further derive estimates of the mass of the Galaxy using a simple mass model with two options for the mass profile of the dark matter halo: an unaltered and an adiabatically contracted Navarro, Frenk & White (NFW) sphere. If we fix the local circular velocity the latter profile yields a significantly higher mass than the un-contracted halo, but if we instead use the statistics on halo concentration parameters in large cosmological simulations as a constraint we find very similar masses for both models. Our best estimate for $M_{340}$, the mass interior to $R_{340}$ (dark matter and baryons), is $1.3^{+0.4}_{-0.3} \times 10^{12}$ M$_\odot$ (corresponding to $M_{200} = 1.6^{+0.5}_{-0.4} \times 10^{12}$ M$_\odot$). This estimate is in good agreement with recently published independent mass estimates based on the kinematics of more distant halo stars and the satellite galaxy Leo I.

We present an update to the multiphase SPH galaxy formation code by Scannapieco et al. We include a more elaborate treatment of the production of metals, cooling rates based on individual element abundances, and a scheme for the turbulent diffusion of metals. Our SN feedback model now transfers energy to the ISM in kinetic and thermal form, and we include a prescription for the effects of radiation pressure from massive young stars on the ISM. We calibrate our new code on the well studied Aquarius haloes and then use it to simulate a sample of 16 galaxies with halo masses between 1x10^11 and 3x10^12 M_sun. In general, the stellar masses of the sample agree well with the stellar mass to halo mass relation inferred from abundance matching techniques for redshifts z=0-4. There is however a tendency to overproduce stars at z>4 and to underproduce them at z<0.5 in the least massive haloes. Overly high SFRs at z<1 for the most massive haloes are likely connected to the lack of AGN feedback in our model. The simulated sample also shows reasonable agreement with observed star formation rates, sizes, gas fractions and gas-phase metallicities at z=0-3. Remaining discrepancies can be connected to deviations from predictions for star formation histories from abundance matching. At z=0, the model galaxies show realistic morphologies, stellar surface density profiles, circular velocity curves and stellar metallicities, but overly flat metallicity gradients. 15 out of 16 of our galaxies contain disk components with kinematic disk fraction ranging between 15 and 65 %. The disk fraction depends on the time of the last destructive merger or misaligned infall event. Considering the remaining shortcomings of our simulations we conclude that even higher kinematic disk fractions may be possible for LambdaCDM haloes with quiet merger histories, such as the Aquarius haloes.

We investigate the chemical and kinematic properties of the diffuse stellar haloes of six simulated Milky Way-like galaxies from the Aquarius Project. Binding energy criteria are adopted to defined two dynamically distinct stellar populations: the diffuse inner and outer haloes, which comprise different stellar sub-populations with particular chemical and kinematic characteristics. Our simulated inner- and outer-halo stellar populations have received contributions from debris stars (formed in sub-galactic systems while they were outside the virial radius of the main progenitor galaxies) and endo-debris stars (those formed in gas-rich sub-galactic systems inside the dark matter haloes). The inner haloes possess an additional contribution from disc-heated stars in the range $\sim 3 - 30 %$, with a mean of $\sim 20% $. Disc-heated stars might exhibit signatures of kinematical support, in particular among the youngest ones. Endo-debris plus disc-heated stars define the so-called \insitu stellar populations. In both the inner- and outer-halo stellar populations, we detect contributions from stars with moderate to low [$\alpha$/Fe] ratios, mainly associated with the endo-debris or disc-heated sub-populations. The observed abundance gradients in the inner-halo regions are influenced by both the level of chemical enrichment and the relative contributions from each stellar sub-population. Steeper abundance gradients in the inner-halo regions are related to contributions from the disc-heated and endo-debris stars, which tend to be found at lower binding energies than debris stars. (Abridged).

We study the properties of two bars formed in fully cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies. In one case, the bar formed in a system with disc, bulge and halo components and is relatively strong and long, as could be expected for a system where the spheroid strongly influences the evolution. The second bar is less strong, shorter, and formed in a galaxy with no significant bulge component. We study the strength and length of the bars, the stellar density profiles along and across the bars and the velocity fields in the bar region. We compare them with the results of dynamical (idealised) simulations and with observations, and find, in general, a good agreement, although we detect some important differences as well. Our results show that more or less realistic bars can form naturally in a $\Lambda$CDM cosmology, and open up the possibility to study the bar formation process in a more consistent way than previously done, since the host galaxies grow, accrete matter and significantly evolve during the formation and evolution of the bar.

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.

We study the chemical properties of the stellar populations in eight simulations of the formation of Milky-Way mass galaxies in a LCDM Universe. Our simulations include metal-dependent cooling and an explicitly multiphase treatment of the effects on the gas of cooling, enrichment and supernova feedback. We search for correlations between formation history and chemical abundance patterns. Differing contributions to spheroids and discs from in situ star formation and from accreted populations are reflected in differing chemical properties. Discs have younger stellar populations, with most stars forming in situ and with low alpha-enhancement from gas which never participated in a galactic outflow. Up to 15 per cent of disc stars can come from accreted satellites. These tend to be alpha-enhanced, older and to have larger velocity dispersions than the in situ population. Inner spheroids have old, metal-rich and alpha-enhanced stars which formed primarily in situ, more than 40 per cent from material recycled through earlier galactic winds. Few accreted stars are found in the inner spheroid unless a major merger occurred recently. Such stars are older, more metal-poor and more alpha-enhanced than the in situ population. Stellar haloes tend to have low metallicity and high alpha-enhancement. The outer haloes are made primarily of accreted stars. Their mean metallicity and alpha-enhancement reflect the masses of the disrupted satellites where they formed: more massive satellites typically have higher [Fe/H] and lower [alpha/Fe]. Surviving satellites have distinctive chemical patterns which reflect their extended, bursty star formation histories. These produce lower alpha-enhancement at given metallicity than in the main galaxy, in agreement with observed trends in the Milky Way.

We present dynamical scaling relations for a homogeneous and representative sample of ~500 massive galaxies, selected only by stellar mass (>10^10 Msun) and redshift (0.025<z<0.05) as part of the ongoing GALEX Arecibo SDSS Survey. We compare baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations for this sample, and investigate how galaxies scatter around the best fits obtained for pruned subsets of disk-dominated and bulge-dominated systems. The BFJ relation is significantly less scattered than the BTF when the relations are applied to their maximum samples, and is not affected by the inclination problems that plague the BTF. Disk-dominated, gas-rich galaxies systematically deviate from the BFJ relation defined by the spheroids. We demonstrate that by applying a simple correction to the stellar velocity dispersions that depends only on the concentration index of the galaxy, we are able to bring disks and spheroids onto the same dynamical relation -- in other words, we obtain a generalized BFJ relation that holds for all the galaxies in our sample, regardless of morphology, inclination or gas content, and has a scatter smaller than 0.1 dex. We find that disks and spheroids are offset in the stellar dispersion-size relation, and that the offset is removed when corrected dispersions are used instead. The generalized BFJ relation represents a fundamental correlation between the global dark matter and baryonic content of galaxies, which is obeyed by all (massive) systems regardless of morphology. [abridged]

We study the stellar discs and spheroids in eight simulations of galaxy formation within Milky Way-mass haloes in a Lambda Cold Dark Matter cosmology. A first paper in this series concentrated on disc properties. Here, we extend this analysis to study how the formation history, structure and dynamics of discs and spheroids relate to the assembly history and structure of their haloes. We find that discs are generally young, with stars spanning a wide range in stellar age: the youngest stars define thin discs and have near-circular orbits, while the oldest stars form thicker discs which rotate ~2 times slower than the thin components, and have 2-3 times larger velocity dispersions. Unlike the discs, spheroids form early and on short time-scales, and are dominated by velocity dispersion. We find great variety in their structure. The inner regions are bar- or bulge-like, while the extended outer haloes are rich in complex non-equilibrium structures such as stellar streams, shells and clumps. Our discs have very high in-situ fractions, i.e. most of their stars formed in the disc itself. Nevertheless, there is a non-negligible contribution (~15 percent) from satellites that are accreted on nearly coplanar orbits. The inner regions of spheroids also have relatively high in-situ fractions, but 65-85 percent of their outer stellar population is accreted. We analyse the circular velocities, rotation velocities and velocity dispersions of our discs and spheroids, both for gas and stars, showing that the dynamical structure is complex as a result of the non-trivial interplay between cooling and SN heating.

We investigate the formation and evolution of dwarf galaxies in a high resolution, hydrodynamical cosmological simulation of a Milky Way sized halo and its environment. Our simulation includes gas cooling, star formation, supernova feedback, metal enrichment and UV heating. In total, 90 satellites and more than 400 isolated dwarf galaxies are formed in the simulation, allowing a systematic study of the internal and environmental processes that determine their evolution. We find that 95% of satellite galaxies are gas-free at z=0, and identify three mechanisms for gas loss: supernova feedback, tidal stripping, and photo-evaporation due to re-ionization. Gas-rich satellite galaxies are only found with total masses above ~ 5x10^9 solar masses. In contrast, for isolated dwarf galaxies, a total mass of ~ 10^9 solar masses constitutes a sharp transition; less massive galaxies are predominantly gas-free at z=0, more massive, isolated dwarf galaxies are often able to retain their gas. In general, we find that the total mass of a dwarf galaxy is the main factor which determines its star formation, metal enrichment, and its gas content, but that stripping may explain the observed difference in gas content between field dwarf galaxies and satellites with total masses close to 10^9 solar masses. We also find that a morphological transformation via tidal stripping of infalling, luminous dwarf galaxies whose dark matter is less concentrated than their stars, cannot explain the high total mass-light ratios of the faint dwarf spheroidal galaxies.

We present cosmological hydrodynamical simulations of the formation of dwarf galaxies in a representative sample of haloes extracted from the Millennium-II Simulation. Our six haloes have a z = 0 mass of ~10^10 solar masses and show different mass assembly histories which are reflected in different star formation histories. We find final stellar masses in the range 5 x 10^7 - 10^8 solar masses, consistent with other published simulations of galaxy formation in similar mass haloes. Our final objects have structures and stellar populations consistent with dwarf elliptical and dwarf irregular galaxies. However, in a Lambda CDM universe, 10^10 solar mass haloes must typically contain galaxies with much lower stellar mass than our simulated objects if they are to match observed galaxy abundances. The dwarf galaxies formed in our own and all other current hydrodynamical simulations are more than an order of magnitude more luminous than expected for haloes of this mass. We discuss the significance and possible implications of this result.

We use cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies to study the relative importance of the main stellar components, i.e., discs, bulges, and bars, at redshift zero. The main aim of this work is to understand if estimates of the structural parameters of these components determined from kinematics (as is usually done in simulations) agree well with those obtained using a photometric bulge/disc/bar decomposition (as done in observations). To perform such a comparison, we have produced synthetic observations of the simulation outputs with the Monte-Carlo radiative transfer code SUNRISE and used the BUDDA code to make 2D photometric decompositions of the resulting images (in the i and g bands). We find that the kinematic disc-to-total ratio (D/T) estimates are systematically and significantly lower than the photometric ones. While the maximum D/T ratios obtained with the former method are of the order of 0.2, they are typically >0.4, and can be as high as 0.7, according to the latter. The photometric decomposition shows that many of the simulated galaxies have bars, with Bar/T ratios in the range 0.2-0.4, and that bulges have in all cases low Sersic indices, resembling observed pseudo-bulges instead of classical ones. Simulated discs, bulges and bars generally have similar (g-i) colours, which are in the blue tail of the distribution of observed colours. This is not due to the presence of young stars, but rather to low metallicities and poor gas content in the simulated galaxies, which makes dust extinction low. Photometric decompositions thus match the component ratios usually quoted for spiral galaxies better than kinematic decompositions, but the shift is insufficient to make the simulations consistent with observed late-type systems.

We have resimulated the six galaxy-sized haloes of the Aquarius Project including metal-dependent cooling, star formation and supernova feedback. This allows us to study not only how dark matter haloes respond to galaxy formation, but also how this response is affected by details of halo assembly history. In agreement with previous work, we find baryon condensation to lead to increased dark matter concentration. Dark matter density profiles differ substantially in shape from halo to halo when baryons are included, but in all cases the velocity dispersion decreases monotonically with radius. Some haloes show an approximately constant dark matter velocity anisotropy with $ \beta \approx 0.1-02$, while others retain the anisotropy structure of their baryon-free versions. Most of our haloes become approximately oblate in their inner regions, although a few retain the shape of their dissipationless counterparts. Pseudo-phase-space densities are described by a power law in radius of altered slope when baryons are included. The shape and concentration of the dark matter density profiles are not well reproduced by published adiabatic contraction models. The significant spread we find in the density and kinematic structure of our haloes appears related to differences in their formation histories. Such differences already affect the final structure in baryon-free simulations, but they are reinforced by the inclusion of baryons, and new features are produced. The details of galaxy formation need to be better understood before the inner dark matter structure of galaxies can be used to constrain cosmological models or the nature of dark matter.

We have studied the dark matter (DM) distribution in a approx 10^12 h^-1 M_sun mass halo extracted from a simulation consistent with the concordance cosmology, where the physics regulating the transformation of gas into stars was allowed to change producing galaxies with different morphologies. The presence of baryons produces the concentration of the DM halo with respect to its corresponding dissipationless run, but we found that this response does not only depend on the amount of baryons gathered in the central region but also on the way they have been assembled. DM and baryons affect each other in a complex way so the formation history of a galaxy plays an important role on its final total mass distribution. Supernova (SN) feedback regulates the star formation and triggers galactic outflows not only in the central galaxy but also in its satellites. Our results suggest that, as the effects of SN feedback get stronger, satellites get less massive and can even be more easily disrupted by dynamical friction, transferring less angular momentum. We found indications that this angular momentum could be acquired not only by the outer part of the DM halo but also by the inner ones and by the stellar component in the central galaxy. The latter effect produces stellar migration which contributes to change the inner potential well, probably working against further DM contraction. As a consequence of the action of these processes, when the halo hosts a galaxy with an important disc structure formed by smooth gas accretion, it is more concentrated than when it hosts a spheroidal system which experienced more massive mergers and interactions. (abridged)

We present results of high resolution hydrodynamical simulations of the formation and evolution of dwarf galaxies. Our simulations start from cosmological initial conditions at high redshift. They include metal-dependent cooling, star formation, feedback from type II and type Ia supernovae and UV background radiation, with physical recipes identical to those applied in a previous study of Milky Way type galaxies. We find that a combination of feedback and the cosmic UV background results in the formation of galaxies with properties similar to the Local Group dwarf spheroidals, and that their effect is strongly moderated by the depth of the gravitational potential. Taking this into account, our models naturally reproduce the observed luminosities and metallicities. The final objects have halo masses between 2.3x10^8 and 1.1x10^9 solar masses, mean velocity dispersions between 6.5 and 9.7 kms-1, stellar masses ranging from 5x10^5 to 1.2x10^7 solar masses, median metallicities between [Fe/H] = -1.8 and -1.1, and half-light radii of the order of 200 to 300 pc, all comparable with Local Group dwarf spheroidals. Our simulations also indicate that the dwarf spheroidal galaxies observed today lie near a halo mass threshold around 10^9 solar masses, in agreement with stellar kinematic data, where supernova feedback not only suffices to completely expel the interstellar medium and leave the residual gas-free, but where the combination of feedback, UV radiation and self-shielding establishes a dichotomy of age distributions similar to that observed in the Milky Way and M31 satellites.

We study the formation of galaxies in a Lambda-CDM Universe using high resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to investigate how the final properties of the simulated galaxies correlate with the formation histories of their haloes. We find that, at z = 0, none of our galaxies contain a disc with more than 20 per cent of its total stellar mass. Four of the eight galaxies nevertheless have well-formed disc components, three have dominant spheroids and very small discs, and one is a spheroidal galaxy with no disc at all. The z = 0 spheroids are made of old stars, while discs are younger and formed from the inside-out. Neither the existence of a disc at z = 0 nor the final disc-to-total mass ratio seems to depend on the spin parameter of the halo. Discs are formed in haloes with spin parameters as low as 0.01 and as high as 0.05; galaxies with little or no disc component span the same range in spin parameter. Except for one of the simulated galaxies, all have significant discs at z > ~2, regardless of their z = 0 morphologies. Major mergers and instabilities which arise when accreting cold gas is misaligned with the stellar disc trigger a transfer of mass from the discs to the spheroids. In some cases, discs are destroyed, while in others, they survive or reform. This suggests that the survival probability of discs depends on the particular formation history of each galaxy. A realistic Lambda-CDM model will clearly require weaker star formation at high redshift and later disc assembly than occurs in our models.

We analyse the dark matter (DM) distribution in a approx 10^12 M_sun halo extracted from a simulation consistent with the concordance cosmology, where the physics regulating the transformation of gas into stars was allowed to change producing galaxies with different morphologies. Although the DM profiles get more concentrated as baryons are collected at the centre of the haloes compared to a pure dynamical run, the total baryonic mass alone is not enough to fully predict the reaction of the DM profile. We also note that baryons affect the DM distribution even outside the central regions. Those systems where the transformation of gas into stars is regulated by Supernova (SN) feedback, so that significant disc structures are able to form, are found to have more concentrated dark matter profiles than a galaxy which has efficiently transformed most of its baryons into stars at early times. The accretion of satellites is found to be associated with an expansion of the dark matter profiles, triggered by angular momentum transfer from the incoming satellites. As the impact of SN feedback increases, the satellites get less massive and are even strongly disrupted before getting close to the main structure causing less angular momentum transfer. Our findings suggest that the response of the DM halo is driven by the history of assembly of baryons into a galaxy along their merger tree.

We study the effects of Supernova (SN) feedback on the formation of galaxies using hydrodynamical simulations in a Lambda-CDM cosmology. We use an extended version of the code GADGET-2 which includes chemical enrichment and energy feedback by Type II and Type Ia SN, metal-dependent cooling and a multiphase model for the gas component. We focus on the effects of SN feedback on the star formation process, galaxy morphology, evolution of the specific angular momentum and chemical properties. We find that SN feedback plays a fundamental role in galaxy evolution, producing a self-regulated cycle for star formation, preventing the early consumption of gas and allowing disks to form at late times. The SN feedback model is able to reproduce the expected dependence on virial mass, with less massive systems being more strongly affected.

We use cosmological simulations in order to study the effects of supernova (SN) feedback on the formation of a Milky Way-type galaxy of virial mass ~10^12 M_sun/h. We analyse a set of simulations run with the code described by Scannapieco et al. (2005, 2006), where we have tested our star formation and feedback prescription using isolated galaxy models. Here we extend this work by simulating the formation of a galaxy in its proper cosmological framework, focusing on the ability of the model to form a disk-like structure in rotational support. We find that SN feedback plays a fundamental role in the evolution of the simulated galaxy, efficiently regulating the star formation activity, pressurizing the gas and generating mass-loaded galactic winds. These processes affect several galactic properties such as final stellar mass, morphology, angular momentum, chemical properties, and final gas and baryon fractions. In particular, we find that our model is able to reproduce extended disk components with high specific angular momentum and a significant fraction of young stars. The galaxies are also found to have significant spheroids composed almost entirely of stars formed at early times. We find that most combinations of the input parameters yield disk-like components, although with different sizes and thicknesses, indicating that the code can form disks without fine-tuning the implemented physics. We also show how our model scales to smaller systems. By analysing simulations of virial masses 10^9 M_sun/h and 10^10 M_sun/h, we find that the smaller the galaxy, the stronger the SN feedback effects.

We use cosmological simulations to study the effects of supernova (SN) feedback on the dark matter distribution in galaxies. We simulate the formation of a Milky-Way type galaxy using a version of the SPH code GADGET2 which includes chemical enrichment and energy feedback by SN, a multiphase model for the gas component and metal-dependent cooling. We analyse the impact of the main three input SN feedback parameters on the amplitude and shape of the dark matter density profiles, focusing on the inner regions of the halo. In order to test the dependence of the results on the halo mass, we simulated a scale-down version of this system. First results of this ongoing work show that the dark matter distribution is affected by the feedback, through the redistribution of the baryons. Our findings suggest that the response of the dark matter halo could be the result of a combination of several physical parameters such as the amount of stellar mass formed at the centre, its shape, and probably the bursty characteristics of the star formation rate. As expected, we find that the dark matter haloes of small galaxies are more sensitive to SN feedback. Higher resolution simulations are being performed to test for numerical effects.

We developed a Monte Carlo code to generate long-duration gamma ray burst (LGRB) events within cosmological hydrodynamical simulations consistent with the concordance model. As structure is assembled, LGRBs are generated in the substructure that formed galaxies today. We adopted the collapsar model so that LGRBs are produced by single, massive stars at the final stage of their evolution. We found that the observed properties of the LGRB host galaxies (HGs) are reproduced if LGRBs are also required to be generated by low metallicity stars. The low metallicity condition imposed on the progenitor stars of LGRBs selects a sample of HGs with mean gas abundances of 12 + log O/H \~ 8.6. For z<1 the simulated HGs of low metallicity LGRB progenitors tend to be faint, slow rotators with high star formation efficiency, compared with the general galaxy population, in agreement with observations. At higher redshift, our results suggest that larger systems with high star formation activity could also contribute to the generation of LGRBs from low metallicity progenitors since the fraction of low metallicity gas available for star formation increases for all systems with look-back time. Under the hypothesis of our LGRB model, our results support the claim that LGRBs could be unbiased tracers of star formation at high redshifts.

We analyse the evolutionary history of galaxies formed in a hierarchical scenario consistent with the concordance $\Lambda$-CDM model focusing on the study of the relation between their chemical and dynamical properties. Our simulations consistently describe the formation of the structure and its chemical enrichment within a cosmological context. Our results indicate that the luminosity-metallicity (LZR) and the stellar mass-metallicity (MZR) relations are naturally generated in a hierarchical scenario. Both relations are found to evolve with redshift. In the case of the MZR, the estimated evolution is weaker than that deduced from observational works by approximately 0.10 dex. We also determine a characteristic stellar mass, $M_c \approx 3 \times 10^{10} M_{\odot}$, which segregates the simulated galaxy population into two distinctive groups and which remains unchanged since $z\sim 3$, with a very weak evolution of its metallicity content. The value and role played by $M_c$ is consistent with the characteristic mass estimated from the SDSS galaxy survey by Kauffmann et al. (2004). Our findings suggest that systems with stellar masses smaller than $M_c$ are responsible for the evolution of this relation at least from $ z\approx 3$. Larger systems are stellar dominated and have formed more than 50 per cent of their stars at $z \ge 2$, showing very weak evolution since this epoch. We also found bimodal metallicity and age distributions from $z\sim3$, which reflects the existence of two different galaxy populations. Although SN feedback may affect the properties of galaxies and help to shape the MZR, it is unlikely that it will significantly modify $M_c$ since, from $z=3$ this stellar mass is found in systems with circular velocities larger than $100 \kms$.

We have developed a new scheme to treat a multiphase interstellar medium in smoothed particle hydrodynamics simulations of galaxy formation. This scheme can represent a co-spatial mixture of cold and hot ISM components, and is formulated without scale-dependent parameters. It is thus particularly suited to studies of cosmological structure formation where galaxies with a wide range of masses form simultaneously. We also present new algorithms for energy and heavy element injection by supernovae, and show that together these schemes can reproduce several important observed effects in galaxy evolution. Both in collapsing systems and in quiescent galaxies our codes can reproduce the Kennicutt relation between the surface densities of gas and of star formation. Strongly metal-enhanced winds are generated in both cases with ratios of mass-loss to star formation which are similar to those observed. This leads to a self-regulated cycle for star formation activity. The overall impact of feedback depends on galaxy mass. Star formation is suppressed at most by a factor of a few in massive galaxies, but in low-mass systems the effects can be much larger, giving star formation an episodic, bursty character. The larger the energy fraction assumed available in feedback, the more massive the outflows and the lower the final stellar masses. Winds from forming disks are collimated perpendicular to the disk plane, reach velocities up to 1000 km/s, and efficiently transport metals out of the galaxies. The asymptotically unbound baryon fraction drops from >95 per cent to ~30 per cent from the least to the most massive of our idealised galaxies, but the fraction of all metals ejected with this component exceeds 60 per cent regardless of mass. Such winds could plausibly enrich the intergalactic medium to observed levels.

We discuss a model for treating chemical enrichment by SNII and SNIa explosions in simulations of cosmological structure formation. Our model includes metal-dependent radiative cooling and star formation in dense collapsed gas clumps. Metals are returned into the diffuse interstellar medium by star particles using a local SPH smoothing kernel. A variety of chemical abundance patterns in enriched gas arise in our treatment owing to the different yields and lifetimes of SNII and SNIa progenitor stars. In the case of SNII chemical production, we adopt metal-dependent yields. Because of the sensitive dependence of cooling rates on metallicity, enrichment of galactic haloes with metals can in principle significantly alter subsequent gas infall and the build up of the stellar components. Indeed, in simulations of isolated galaxies we find that a consistent treatment of metal-dependent cooling produces 25% more stars outside the central region than simulations with a primordial cooling function. In the highly-enriched central regions, the evolution of baryons is however not affected by metal cooling, because here the gas is always dense enough to cool. A similar situation is found in cosmological simulations because we include no strong feedback processes which could spread metals over large distances and mix them into unenriched diffuse gas. We demonstrate this explicitly with test simulations which adopt super-solar cooling functions leading to large changes both in the stellar mass and in the metal distributions. We also find that the impact of metallicity on the star formation histories of galaxies may depend on their particular evolutionary history. Our results hence emphasise the importance of feedback processes for interpreting the cosmic metal enrichment.

Gamma ray bursts (GRBs) belong to the most energetic events in the Universe. Recently, the extragalactic nature of these sources has been confirmed with the discovery of several host galaxies (HGs) and the measurement of their redshifts. To explain the origin of GRBs various models have been proposed, among which the coalescence of compact objects and the "collapsar" scenarios are the most representative, being the collapsar model one of the most accepted to explain the long duration GRBs. A natural consequence of this model is that the GRBs would trace the star formation rate (SFR) of their HGs. In this contributed paper we present preliminary results of the development of a Montecarlo-based code for collapsar event formation which is coupled to chemical-cosmological simulations aiming at studying the properties of HGs in a hierarchical scenario.

We perform an statistical analysis of galaxies in pairs in a Lambda-CDM scenario by using the chemical GADGET-2 of Scannapieco et al. (2005) in order to study the effects of galaxy interactions on colours and metallicities. We find that galaxy-galaxy interactions can produce a bimodal colour distribution with galaxies with significant recent star formation activity contributing mainly to blue colours. In the simulations, the colours and the fractions of recently formed stars of galaxies in pairs depend on environment more strongly than those of galaxies without a close companion, suggesting that interactions play an important role in galaxy evolution. If the metallicity of the stellar populations is used as the chemical indicator, we find that the simulated galaxies determine luminosity-metallicity and stellar mass-metallicity relations which do not depend on the presence of a close companion. However, in the case of the luminosity-metallicity relation, at a given level of enrichment, we detect a systematic displacement of the relation to brighter magnitudes for active star forming systems. Regardless of relative distance and current level of star formation activity, galaxies in pairs have stellar populations with higher level of enrichment than galaxies without a close companion. In the case of the gas component, this is no longer valid for galaxies in pairs with passive star formation which only show an excess of metals for very close pair members, consequence of an important recent past star formation activity. (Abridged).

In this work, we study the Luminosity-Metallicity relation (LMR) and the Stellar Mass -Metallicity relation (MMR) of galactic systems in a hierarhical clustering scenario. We performed numerical hydrodynamical simulations with the chemical GADGET-2 of Scannapieco et al.(2005) in a LCDM universe. We found that our simulated galactic systems reproduce the observed local LMR and its evolution in zero point and slope. The simulated MMR is also in agreement with recent observational results. From the analysis of the evolution of the MMR, we found a characteristic mass at ~10^(10.2) M_sun / h which separates two galactic populations with different astrophysical properties. More massive systems tend to have their stars formed at z > 2 and show less evolution than smaller systems. Hence, this characteristic mass is determined by the formation of the structure in a hierarchical scenario. Our results also suggest the need for important supernova feedback.

We present a study of the effects of mergers and interactions on the mass distribution of galactic systems in hierarchical clustering scenarios using the disc-bulge structural parameters and their dynamical properties to quantify them. We focus on the analysis of the Fundamental Mass Plane relation, finding that secular evolution phases contribute significantly to the determination of a plane with a slope in agreement to that of the observed luminosity relation. In these simulations, secular phases are responsible for the formation of compact stellar bulges with the correct structural parameter combination. We also test that the relations among these parameters agree with observations. Taking into account these results, the hierarchical growth of the structure predicts a bulge formation scenario for typical field spiral galaxies where secular evolution during dissipate mergers plays a fundamental role. Conversely, subsequent mergers can help to enlarge the bulges but do not seem to strongly modify their fundamental mass relations. Systems get to the local mass relations at different stages of evolution (i.e. different redshifts) so that their formation histories introduce a natural scatter in the relations. Our results suggest that the formation mechanisms of the bulge and disc components, satisfying their corresponding fundamental mass relations, might be coupled and that secular evolution could be the possible connecting process.

We study the Luminosity-Metallicity Relation and the Stellar Mass-Metallicity Relation by performing chemo-dynamical simulations in a cosmological scenario. Our results predict a tight linear correlation between oxygen chemical abundance and luminosity for galactic systems up to z=3. The Luminosity-Metallicity Relation evolves with redshift in such a way that systems at a given oxygen abundance were ~3 dex brighter at z=3 compared to local ones, in good agreement with observations. We determin also a characteristic stellar mass M_c~10^(10.2) M_sun/h above which the Stellar Mass-Metallicity Relation starts to flatten. We encounter that this mass arises as a consequence of the hierarchical aggregation of structure in a LCDM universe and segregates two galactic populations with different astophysical properties.

We present a Supernova (SN) feedback model that succeeds at describing the chemical and energetic effects of SN explosions in galaxy formation simulations. This new SN model has been coupled to GADGET-2 and works within a new multiphase scheme which allows the description of a co-spatial mixture of cold and hot interstellar medium phases. No ad hoc scale-dependent parameters are associated to these SN and multiphase models making them particularly suited to studies of galaxy formation in a cosmological framework. Our SN model succeeds not only in setting a self-regulated star formation activity in galaxies but in triggering collimated chemical-enriched galactic winds. The effects of winds vary with the virial mass of the systems so that the smaller the galaxy, the larger the fraction of swept away gas and the stronger the decrease in its star formation activity. The fact that the fraction of ejected metals exceeds 60 per cent regardless of mass, suggests that SN feedback can be the responsible mechanism of the enrichment of the intergalactic medium to the observed levels.

We study the effects of Supernova (SN) feedback on the formation of disc galaxies. For that purpose we run simulations using the extended version of the code GADGET-2 which includes a treatment of chemical and energy feedback by SN explosions. We found that our model succeeds in setting a self-regulated star formation process since an important fraction of the cold gas from the center of the haloes is efficiently heated up and transported outwards. The impact of SN feedback on galactic systems is also found to depend on virial mass: smaller systems are more strongly affected with star formation histories in which several starbursts can develop. Our implementation of SN feedback is also successful in producing violent outflows of chemical enriched material.

We study the evolution of the Stellar Mass-Metallicity Relation and the Luminosity-Metallicity Relation by performing numerical simulations in a cosmological framework. We find that the slope and the zero point of the Luminosity-Metallicity Relation evolve in such a way that, at a given metallicity, systems were ~3 mag brighter at z=3 compared to galaxies in the local universe, which is consistent with the observational trend. The local Stellar Mass-Metallicity Relation shows also a good agreement with recent observations. We identify a characteristic stellar mass M_c ~ 10^(10.2) M_sun/h at which the slope of the Stellar Mass-Metallicity Relation decreases for larger stellar masses. Our results indicate that M_c arises naturally as a consequence of the hierarchical building up of the structure.

We carry out a statistical analysis of galaxy pairs selected from chemical hydrodynamical simulations with the aim at assessing the capability of hierarchical scenarios to reproduce recent observational results for galaxies in pairs. Particularly, we analyse the effects of mergers and interactions on the star formation (SF) activity, the global mean chemical properties and the colour distribution of interacting galaxies. We also assess the effects of spurious pairs.

Using chemical hydrodynamical simulations consistent with a Lambda-CDM model, we study the role played by mergers and interactions in the regulation of the star formation activity, colours and the chemical properties of galaxies in pairs. A statistical analysis of the orbital parameters in galaxy pairs (r <100 kpc/h) shows that the star formation (SF) activity correlates strongly with the relative separation and weakly with the relative velocity, indicating that close encounters (r <30 kpc/h) can increase the SF activity to levels higher than that exhibit in galaxies without a close companion. Analysing the internal properties of interacting systems, we find that their stability properties also play a role in the regulation the SF activity (Perez et al 2005a). Particularly, we find that the passive star forming galaxies in pairs are statistically more stable with deeper potential wells and less leftover gas than active star forming pairs. In order to compare our results with observations, we also build a projected catalog of galaxy pairs (2D-GP: rp <100 kpc/h and Vr <350 km/s), constructed by projecting the 3D sample in different random directions. In good agreement with observations (Lambas et al 2003), our results indicate that galaxies with rp < 25 kpc/h (close pairs) show an enhancement of the SF activity with respect to galaxies without a close companion. (Abridged.)

We carried out a statistical analysis of galaxy pairs in hydrodynamical Lambda-CDM simulations. We focused on the triggering of star formation by interactions and analysed the enhancement of star formation activity in terms of orbital parameters. By comparing to a suitable sample of simulated galaxies without a nearby companion, we find that close encounters (r<30 kpc/h) may effectively induce star formation. However, our results suggest that the stability properties of systems and the spatial proximity are both relevant factors in the process of triggering star formation by tidal interactions. In order to assess the effects of projection and spurious pairs in observational samples, we also constructed and analysed samples of pairs of galaxies in the simulations obtained in projection. We found a good agreement with observational results with a threshold at rp ~ 25 kpc/h for interactions to effectively enhance star formation activity. For pairs within rp < 100 kpc/h, we estimated a ~27% contamination by spurious pairs, reduced to ~19% for close systems. We also found that spurious pairs affect more strongly high density regions with 17% of spurious pairs detected for low density regions compared to 33% found in high density ones. Also, we analysed the dependence of star formation on environment by defining the usual projected density parameter for both pairs and isolated galaxies in the simulations. We find the expected star formation-local density relation for both galaxies in pairs and without a close companion, with a stronger density dependence for close pairs which suggests a relevant role for interactions in driving this relation.

We investigate the generation of galactic outflows by supernova feedback in the context of SPH cosmological simulations. We use a modified version of the code GADGET-2 which includes chemical enrichment and energy feedback by Supernova. We find that energy feedback plays a fundamental role in the evolution of galaxies, heating up the cold material in the centre of the haloes and triggering outflows which efficiently transport gas from the centre to the outskirts of galaxies. The impact of feedback is found to depend on the virial mass of the system with smaller systems, such as dwarf galaxies, being more strongly affected. The outflows help to establish a self-regulated star formation process, and to transport a significant amount of metals into the haloes and even out of the systems. According to our results, energy feedback by supernovae could be the mechanism responsible for the chemical enrichment of the intergalactic medium.

We investigate the star formation activity in galaxy pairs in chemical hydrodynamical simulations consistent with a Lambda-CDM scenario. A statistical analysis of the effects of galaxy interactions on the star formation activity as a function of orbital parameters shows that close encounters (r < 30 kpc/h) can be effectively correlated with an enhancement of star formation activity with respect to galaxies without a close companion. Our results suggest that the stability properties of systems are also relevant in this process. We found that the passive star forming galaxies pairs tend to have deeper potential wells, older stellar populations, and less leftover gas than active star forming ones. In order to assess the effects that projection and interlopers could introduce in observational samples, we have also constructed and analysed projected simulated catalogs of galaxy pairs. In good agreement with observations, our results show a threshold (rp < 25 kpc/h) for interactions to enhance the star formation activity with respect to galaxies without a close companion. Finally, analysing the environmental effect, we detect the expected SFR-local density relation for both pairs and isolated galaxy samples, although the density dependence is stronger for galaxies in pairs suggesting a relevant role for interactions in driving this relation.

We study the mass-metallicity relation of galactic systems with stellar masses larger than 10^9 Mo in Lambda-CDM scenarios by using chemical hydrodynamical simulations. We find that this relation arises naturally as a consequence of the formation of the structure in a hierarchical scenario. The hierarchical building up of the structure determines a characteristic stellar mass at M_c ~10^10.2 Moh^-1 which exhibits approximately solar metallicities from z ~ 3 to z=0. This characteristic mass separates galactic systems in two groups with massive ones forming most of their stars and metals at high redshift. We find evolution in the zero point and slope of the mass-metallicity relation driven mainly by the low mass systems which exhibit the larger variations in the chemical properties. Although stellar mass and circular velocity are directly related, the correlation between circular velocity and metallicity shows a larger evolution with redshift making this relation more appropriate to confront models and observations. The dispersion found in both relations is a function of the stellar mass and reflects the different dynamical history of evolution of the systems.

We use cosmological SPH simulations to investigate the effects of mergers and interactions on the formation of the bulge and disc components of galactic systems. We find that secular evolution during mergers seems to be a key process in the formation of stable disc-bulge systems with observational counterparts and contributes to establish the fundamental relations observed in galaxies. Our findings suggest that the secular evolution phase couples the formation mechanisms of the bulge and disc components. According to our results, depending on the particular stability properties and merger parameters, violents events could drive a morphological loop in which the outcome could be a disc or a spheroid.

We present first results of an implementation of chemical evolution in a cosmological hydrodynamical code, focusing the analysis on the effects of cooling baryons according to their metallicity. We found that simulations with primordial cooling can underestimate the star formation rate from $z < 3$ and by up to $\approx 20 %$. We constructed simulated spectra by combining the star formation and chemical history of galactic systems with spectra synthesis models and assess the impact of chemical evolution on the energy distribution.

We study the effects of mergers on the structural properties of disc-like systems by using Smooth Particle Hydrodynamical (SPH) numerical simulations in hierarchical clustering scenarios. In order to assess the effects of mergers on the mass distributions we performed a bulge-disc decomposition of the projected surface density of the systems at different stages of the merger process. We assumed an exponential law for the disc component and the S\'ersic law for the bulges. We found that simulated objects at $z=0$ have bulge profiles with shape parameters $n\approx 1$, consistent with observational results of spiral galaxies. The complete sample of simulated objects at $z=0$ and $z>0$ shows that $n$ takes values in the range $n\approx 0.4 - 4$. We found that secular evolution tends to produce exponential bulge profiles, while the fusion of baryonic cores tends to increase the $n$ value and helps to generate the correlation between $B/D$ and $n$. We found no dependence on the relative mass of the colliding objects. Our results suggest that mergers, through secular evolution and fusions, could produce the transformation of galactic objects along the Hubble sequence by driving a morphological loop that might also depend on the properties of the central galactic potential wells, which are also affected by mergers.

We use cosmological SPH simulations to study the effects of mergers in the star formation history of galactic objects in hierarchical clustering scenarios. We find that during some merger events, gaseous discs can experience two starbursts: the first one during the orbital decay phase, due to gas inflows driven as the satellite approaches, and the second one, when the two baryonic clumps collide. A trend for these first induced starbursts to be more efficient at transforming the gas into stars is also found. We detect that systems which do not experience early gas inflows have well-formed stellar bulges and more concentrated potential wells, which seem to be responsible for preventing further gas inward transport triggered by tidal forces. Our results constitute the first proof that bulges can form as the product of collapse, collisions and secular evolution in a cosmological framework, and they are consistent with a rejuvenation of the stellar population in bulges at intermediate z with, at least, 50% of the stars (in SCDM) being formed at high z. (Abridged)