In recent years a new window on galaxy evolution opened, thanks to the increasing discovery of galaxies with a low surface brightness, such as Ultra Diffuse Galaxies (UDGs). The formation mechanism of these systems is still a much debated question, and so are their kinematical properties. In this work, we address this topic by analyzing the stellar kinematics of isolated UDGs formed in the hydrodynamical simulation suite NIHAO. We construct projected line-of-sight velocity and velocity dispersion maps to compute the projected specific angular momentum, $\lambda_{\rm R}$, to characterize the kinematical support of the stars in these galaxies. We found that UDGs cover a broad distribution, ranging from dispersion to rotation supported galaxies, with similar abundances in both regimes. The degree of rotation support of simulated UDGs correlates with several properties such as galaxy morphology, higher HI fractions and larger effective radii with respect to the dispersion supported group, while the dark matter halo spin and mass accretion history are similar amongst the two populations. We demonstrate that the alignment of the infalling baryons into the protogalaxy at early $z$ is the principal driver of the $z$=0 stellar kinematic state: pressure supported isolated UDGs form via mis-aligned gas accretion while rotation supported ones build-up their baryons in an ordered manner. Accounting for random inclination effects, we predict that a comprehensive survey will find nearly half of field UDGs to have rotationally supported stellar disks.

The Milky Way underwent its last significant merger ten billion years ago, when the Gaia-Enceladus-Sausage (GES) was accreted. Accreted GES stars and progenitor stars born prior to the merger make up the bulk of the inner halo. Even though these two main populations of halo stars have similar $durations$ of star formation prior to their merger, they differ in [$\alpha$/Fe]-[Fe/H] space, with the GES population bending to lower [$\alpha$/Fe] at a relatively low value of [Fe/H]. We use cosmological simulations of a 'Milky Way' to argue that the different tracks of the halo stars through the [$\alpha$/Fe]-[Fe/H] plane are due to a difference in their star formation history and efficiency, with the lower mass GES having its low and constant star formation regulated by feedback whilst the higher mass main progenitor has a higher star formation rate prior to the merger. The lower star formation efficiency of GES leads to lower gas pollution levels, pushing [$\alpha$/Fe]-[Fe/H] tracks to the left. In addition, the increasing star formation rate maintains a higher relative contribution of Type~II SNe to Type~Ia SNe for the main progenitor population that formed during the same time period, thus maintaining a relatively high [$\alpha$/Fe]. Thus the different positions of the downturns in the [$\alpha$/Fe]-[Fe/H] plane for the GES stars are not reflective of different star formation durations, but instead reflect different star formation efficiencies. We argue that cosmological simulations match a wide range of independent observations, breaking degeneracies that exist in simpler models.

We use Gaia DR2 systemic proper motions of 45 satellite galaxies to constrain the mass of the Milky Way using the scale free mass estimator of Watkins et al. (2010). We first determine the anisotropy parameter $\beta$, and the tracer satellites' radial density index $\gamma$ to be $\beta$=$-0.67^{+0.45}_{-0.62}$ and $\gamma=2.11\pm0.23$. When we exclude possible former satellites of the Large Magellanic Cloud, the anisotropy changes to $\beta$=$-0.21^{+0.37}_{-0.51}$. We find that the index of the Milky Way's gravitational potential $\alpha$, which is dependent on the mass itself, is the parameter with the largest impact on the mass determination. Via comparison with cosmological simulations of Milky Way-like galaxies, we carried out a detailed analysis of the estimation of the observational uncertainties and their impact on the mass estimator. We found that the mass estimator is biased when applied naively to the satellites of simulated Milky Way halos. Correcting for this bias, we obtain for our Galaxy a mass of $0.58^{+0.15}_{-0.14}\times10^{12}$M$_\odot$ within 64 kpc, as computed from the inner half of our observational sample, and $1.43^{+0.35}_{-0.32}\times10^{12}$M$_\odot$ within 273 kpc, from the full sample; this latter value extrapolates to a virial mass of $M_\mathrm{vir\,\Delta=97}$=$1.51^{+0.45}_{-0.40} \times 10^{12}M_{\odot}$ corresponding to a virial radius of R$_\mathrm{vir}$=$308\pm29$ kpc. This value of the Milky Way mass lies in-between other mass estimates reported in the literature, from various different methods.

Knowledge of ages for stars formed over a galaxy's lifetime is fundamental to understand its formation and evolution. However, stellar ages are difficult to obtain since they cannot be measured from observations, being comparison with stellar models (Soderblom 2010) required. Alternatively, age distributions can be derived applying the robust technique of colour-magnitude diagram fitting (Gallart et al. 2005), till now mainly employed to study nearby galaxies. The new distances to individual Milky Way stars from the Gaia mission (Brown et al. 2018) have allowed us to use this technique to derive ages from a thick disk colour-magnitude diagram, and from the enigmatic, two-sequenced colour-magnitude diagram of the kinematically hot local halo (Babusiaux et al. 2018), which blue-sequence has been linked to a major accretion event (Haywood et al. 2018, Helmi et al. 2018). Because accurate ages were lacking, the time of the merger and its role on our Galaxy's early evolution remained unclear. We show that the stars in both halo sequences share identical age distributions, and are older than the bulk of thick disc stars. The sharp halo age cut 10 Gyr ago can be identified with the accretion of Gaia-Enceladus. Along with state-of-the-art cosmological simulations of galaxy formation (Brook et al. 2012), these robust ages allow us to order the early sequence of events that shaped our Galaxy, identifying the red-sequence as the first stars formed within the Milky Way progenitor which, because of their kinematics, can be described as its long sought in-situ halo.

We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter density slope, $\alpha$, on the ratio between stellar-to-halo mass. We show that this relation holds approximately at all redshifts (with an intrinsic scatter of ~0.18 in $\alpha$). This implies that in practically all haloes the shape of their inner density profile changes quite substantially over cosmic time, as they grow in stellar and total mass. Thus, depending on their final stellar-to-halo mass ratio, haloes can either form and keep a substantial density core (size~1 kpc), or form and then destroy the core and re-contract the halo, going back to a cuspy profile, which is even steeper than CDM predictions for massive galaxies (~1e12 Msun). We show that results from the NIHAO suite are in good agreement with recent observational measurements of $\alpha$ in dwarf galaxies. Overall our results suggest that the notion of a universal density profile for dark matter haloes is no longer valid in the presence of galaxy formation.

The existence of galaxies with a surface brightness $\mu$ lower than the night sky has been known since three decades. Yet, their formation mechanism and emergence within a $\rm\Lambda CDM$ universe has remained largely undetermined. For the first time, we investigated the origin of Low Surface Brightness (LSB) galaxies with M$_{\star}$$\sim$10$^{9.5-10}$M$_{\odot}$, which we are able to reproduce within hydrodynamical cosmological simulations from the NIHAO suite. The simulated and observed LSBs share similar properties, having large HI reservoir, extended star formation histories and effective radii, low S\'{e}rsic index and slowly rising rotation curves. The formation mechanism of these objects is explored: simulated LSBs form as a result of co-planar co-rotating mergers and aligned accretion of gas at early times, while perpendicular mergers and mis-aligned gas accretion result in higher $\mu$ galaxies by $z$=0. The larger the merger, the stronger the correlation between merger orbital configuration and final $\mu$. While the halo spin parameter is consistently high in simulated LSB galaxies, the impact of halo concentration, feedback-driven gas outflows and merger time only plays a minor-to-no role in determining $\mu$. Interestingly, the formation scenario of such `classical' LSBs differs from the one of less massive, M$_{\star}$$\sim$10$^{7-9}$M$_{\odot}$, Ultra-Diffuse Galaxies, the latter resulting from the effects of SNae driven gas outflows: a M$_{\star}$ of $\sim$10$^9$M$_{\odot}$ thus represents the transition regime between a feedback-dominated to an angular momentum-dominated formation scenario in the LSB realm. Observational predictions are offered regarding spatially resolved star formation rates through LSB discs: these, together with upcoming surveys, can be used to verify the proposed emergence scenario of LSB galaxies.

A proper understanding of the Milky Way (MW) dwarf galaxies in a cosmological context requires knowledge of their 3D velocities and orbits. However, proper motion (PM) measurements have generally been of limited accuracy and available only for more massive dwarfs. We therefore present a new study of the kinematics of the MW dwarf galaxies. We use the Gaia DR2 for those dwarfs that have been spectroscopically observed in the literature. We derive systemic PMs for 39 galaxies and galaxy candidates out to 420 kpc, and generally find good consistency for the subset with measurements available from other studies. We derive the implied Galactocentric velocities, and calculate orbits in canonical MW halo potentials of "low" ($0.8 \times 10^{12} M_\odot$) and "high" mass ($1.6 \times 10^{12} M_\odot$). Comparison of the distributions of orbital apocenters and 3D velocities to the halo virial radius and escape velocity, respectively, suggests that the satellite kinematics are best explained in the high-mass halo. Tuc III, Crater II, and additional candidates have orbital pericenters small enough to imply significant tidal influences. Relevant to the missing satellite problem, the fact that fewer galaxies are observed to be near apocenter than near pericenter implies that there must be a population of distant dwarf galaxies yet to be discovered. Of the 39 dwarfs: 12 have orbital poles that do not align with the MW plane of satellites (given reasonable assumptions about its intrinsic thickness); 10 have insufficient PM accuracy to establish whether they align; and 17 satellites align, of which 11 are co-orbiting and (somewhat surprisingly, in view of prior knowledge) 6 are counter-orbiting. Group infall might have contributed to this, but no definitive association is found for the members of the Crater-Leo group.

HST/ACS observations along the major axis of M33 show that the mean age of its stars decreases with increasing distance from the galaxy center. Such a behavior is consistent with an inside-out growth of the disc. However, in the outermost observed field, at r$\simeq$11.6 kpc, a reversal of this gradient is detected, with old stars found in high percentages beyond this radius. In this work we investigate the origin of such a reversal in stellar age gradient, by using a simulated M33 analogue from the Constrained Local UniversE Simulations (CLUES). The simulated M33 is similar to the observed one in terms of mass, rotation velocity, surface brightness and, similar to what has been reported in observations, shows a stellar age turnaround at large radii. We demonstrate that this reversal is mostly a result of stellar accretion from old satellite galaxies and, to a lesser extent, of stellar migration of in-situ stars. The old accreted stars, with formation times tf < 4 Gyrs, are kinematically hot and can be differentiated from the in-situ stars by their high velocity dispersion and the fact that they do not have rotationally-supported orbits. In the future, obtaining kinematic information of the stars in the outskirt of M33 will help to verify this scenario.

In a $\Lambda$CDM cosmology, galaxy formation is a globally inefficient process: it is often the case that far fewer baryons are observed in galaxy disks than expected from the cosmic baryon fraction. The location of these "missing baryons" is unclear. By fitting halo profiles to the rotation curves of galaxies in the SPARC data set, we measure the "missing baryon" mass for individual late-type systems. Assuming that haloes initially accrete the cosmological baryon fraction, we show that the maximum energy available from supernovae is typically not enough to completely eject these "missing baryons" from a halo, but it is often sufficient to heat them to the virial temperature. The energy available from supernovae has the same scaling with galaxy mass as the energy needed to heat or eject the "missing baryons", indicating that the coupling efficiency of the feedback to the ISM may be constant with galaxy virial mass. We further find that the energy available from supernova feedback is always enough to convert a primordial cusp into a core and has magnitude consistent with what is required to heat the "missing baryons" to the virial temperature. Taking a census of the baryon content of galaxies with ${\rm 10^9<M_{vir}/M_{\odot}<10^{12}}$ reveals that $\sim86\%$ of baryons are likely to be in a hot phase surrounding the galaxies and possibly observable in the X-ray, $\sim7\%$ are in the form of cold gas, and $\sim7\%$ are in stars.

Understanding the peculiar properties of Ultra Diffuse Galaxies (UDGs) via spectroscopic analysis is a challenging task requiring very deep observations and exquisite data reduction. In this work we perform one of the most complete characterisations of the stellar component of UDGs to date using deep optical spectroscopic data from OSIRIS at GTC. We measure radial and rotation velocities, star formation histories (SFH) and mean population parameters, such as ages and metallicities, for a sample of five UDG candidates in the Coma cluster. From the radial velocities, we confirm the Coma membership of these galaxies. We find that their rotation properties, if detected at all, are compatible with dwarf-like galaxies. The SFHs of the UDG are dominated by old (~ 7 Gyr), metal-poor ([M/H] ~ -1.1) and alpha-enhanced ([Mg/Fe] ~ 0.4) populations followed by a smooth or episodic decline which halted ~ 2 Gyr ago, possibly a sign of cluster-induced quenching. We find no obvious correlation between individual SFH shapes and any UDG morphological properties. The recovered stellar properties for UDGs are similar to those found for DDO44, a local UDG analogue resolved into stars. We conclude that the UDGs in our sample are extended dwarfs whose properties are likely the outcome of both internal processes, such as bursty SFHs and/or high-spin haloes, as well as environmental effects within the Coma cluster.

According to star formation histories (SFHs), Local Group dwarf galaxies can be broadly classified in two types: those forming most of their stars before $z=2$ (${\it fast}$) and those with more extended SFHs (${\it slow}$). The most precise SFHs are usually derived from deep but not very spatially extended photometric data; this might alter the ratio of old to young stars when age gradients are present. Here we correct for this effect and derive the mass formed in stars by $z=2$ for a sample of 16 Local Group dwarf galaxies. We explore early differences between ${\it fast}$ and ${\it slow}$ dwarfs, and evaluate the impact of internal feedback by supernovae (SN) on the baryonic and dark matter (DM) component of the dwarfs. ${\it Fast}$ dwarfs assembled more stellar mass at early times and have larger amounts of DM within the half-light radius than ${\it slow}$ dwarfs. By imposing that ${\it slow}$ dwarfs cannot have lost their gas by $z=2$, we constrain the maximum coupling efficiency of SN feedback to the gas and to the DM to be $\sim$10%. We find that internal feedback alone appears insufficient to quench the SFH of ${\it fast}$ dwarfs by gas deprivation, in particular for the fainter systems. Nonetheless, SN feedback can core the DM halo density profiles relatively easily, producing cores of the sizes of the half-light radius in ${\it fast}$ dwarfs by $z=2$ with very low efficiencies. Amongst the "classical" Milky Way satellites, we predict that the smallest cores should be found in Draco and Ursa Minor, while Sculptor and Fornax should host the largest ones.

The significant diversity of rotation curve (RC) shapes in dwarf galaxies has recently emerged as a challenge to LCDM: in dark matter (DM) only simulations, DM halos have a universal cuspy density profile that results in self-similar RC shapes. We compare RC shapes of simulated galaxies from the NIHAO project with observed galaxies from the homogeneous SPARC dataset. The DM halos of the NIHAO galaxies can expand to form cores, with the degree of expansion depending on their stellar-to-halo mass ratio. By means of the $V_{\rm 2kpc}-V_{\rm Rlast}$ relation (where $V_{\rm Rlast}$ is the outermost measured rotation velocity), we show that both the average trend and the scatter in RC shapes of NIHAO galaxies are in reasonable agreement with SPARC: this represents a significant improvement compared to simulations that do not result in DM core formation, suggesting that halo expansion is a key process in matching the diversity of dwarf galaxy RCs. Note that NIHAO galaxies can reproduce even the extremely slowly rising RCs of IC 2574 and UGC 5750. Revealingly, the range where observed galaxies show the highest diversity corresponds to the range where core formation is most efficient in NIHAO simulations, 50$<V_{\rm Rlast}$\km s$^{-1}<$100. A few observed galaxies in this range cannot be matched by any NIHAO RC nor by simulations that predict a universal halo profile. Interestingly, the majority of these are starbursts or emission-line galaxies, with steep RCs and small effective radii. Such galaxies represent an interesting observational target providing new clues to the process/viability of cusp-core transformation, the relationship between starburst and inner potential well, and the nature of DM.

We have developed a novel Markov Chain Mote Carlo (MCMC) chemical "painting" technique to explore possible radial and vertical metallicity gradients for the thick disc progenitor. In our analysis we match an N-body simulation to the data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We assume that the thick disc has a constant scale-height and has completed its formation at an early epoch, after which time radial mixing of its stars has taken place. Under these assumptions, we find that the initial radial metallicity gradient of the thick disc progenitor should not be negative, but either flat or even positive, to explain the current negative vertical metallicity gradient of the thick disc. Our study suggests that the thick disc was built-up in an inside-out and upside-down fashion, and older, smaller and thicker populations are more metal poor. In this case, star forming discs at different epochs of the thick disc formation are allowed to have different radial metallicity gradients, including a negative one, which helps to explain a variety of slopes observed in high redshift disc galaxies. This scenario helps to explain the positive slope of the metallicity-rotation velocity relation observed for the Galactic thick disc. On the other hand, radial mixing flattens the slope of an existing gradient.

Recent determinations of the radial distributions of mono-metallicity populations (MMPs, i.e., stars in narrow bins in [Fe/H] within wider [$\alpha$/Fe] ranges) by the SDSS-III/APOGEE DR12 survey cast doubts on the classical thin - thick disk dichotomy. The analysis of these observations lead to the non-$[\alpha$/Fe] enhanced populations splitting into MMPs with different surface densities according to their [Fe/H]. By contrast, $[\alpha$/Fe] enhanced (i.e., old) populations show an homogeneous behaviour. We analyze these results in the wider context of disk formation within non-isolated halos embedded in the Cosmic Web, resulting in a two-phase mass assembly. By performing hydrodynamical simulations in the context of the $\rm \Lambda CDM$ model, we have found that the two phases of halo mass assembly (an early, fast phase, followed by a slow one, with low mass assembly rates) are very relevant to determine the radial structure of MMP distributions, while radial mixing has only a secondary role, depending on the coeval dynamical and/or destabilizing events. Indeed, while the frequent dynamical violent events occuring at high redshift remove metallicity gradients, and imply efficient stellar mixing, the relatively quiescent dynamics after the transition keeps [Fe/H] gaseous gradients and prevents newly formed stars to suffer from strong radial mixing. By linking the two-component disk concept with the two-phase halo mass assembly scenario, our results set halo virialization (the event marking the transition from the fast to the slow phases) as the separating event marking periods characterized by different physical conditions under which thick and thin disk stars were born.

The spectral energy distributions (SEDs) of low-mass low-metallicity (dwarf) galaxies are a challenging piece of the puzzle of galaxy formation in the near Universe. These SEDs show some particular features in the submillimeter to far-infrared wavelength range compared to normal, larger, and metal-richer galaxies that cannot be explained by the current models. These can be summarized as: a broadening of the IR peak, which implies a warmer dust component; an excess of emission in the submm ($\sim$500 $\mu$m), that causes a flattening of the submm/FIR slope; and a very low intensity of PAH emission features. With the aim of explaining these emission patterns, the SEDs of a sample of 27 simulated dwarf galaxies were calculated using the GRASIL-3D radiation transfer code. This code has the particularity that it separately treats the radiative transfer through dust grains within molecular clouds and within the cirrus, the dense and diffuse components of the gas phase, respectively. The simulated galaxies have stellar masses ranging from 10$^6$-10$^9$ M$_\odot$, and were obtained from a single simulation run within a Local Group environment with initial conditions from the CLUES project. We report a study of the IRAS, Spitzer and Herschel bands luminosities, and of the SFRs, dust, and gas (HI and H$_2$) mass contents. We find a satisfactory agreement with observational data, with GRASIL-3D naturally reproducing the specific spectral features mentioned above. We conclude that the GRASIL-3D two-component dust model gives a physical interpretation of the emission of dwarf galaxies: molecular clouds and cirrus represent the warm and cold dust components, respectively, needed to reproduce observational data.

Cosmological $N$-body simulations predict dark matter (DM) haloes with steep central cusps (e.g. NFW, Navarro et al. 1996). This contradicts observations of gas kinematics in low-mass galaxies that imply the existence of shallow DM cores. Baryonic processes such as adiabatic contraction and gas outflows can, in principle, alter the initial DM density profile, yet their relative contributions to the halo transformation remain uncertain. Recent high resolution, cosmological hydrodynamic simulations (Di Cintio et al. 2014, DC14) predict that inner density profiles depend systematically on the ratio of stellar to DM mass (M$_*$/M$_{\text{halo}}$). Using a Markov Chain Monte Carlo approach, we test the NFW and the M$_*$/M$_{\text{halo}}$-dependent DC14 halo models against a sample of 147 galaxy rotation curves from the new {\it Spitzer} Photometry and Accurate Rotation Curves (SPARC) data set. These galaxies all have extended H{\small I} rotation curves from radio interferometry as well as accurate stellar mass density profiles from near-infrared photometry. The DC14 halo profile provides markedly better fits to the data compared to the NFW profile. Unlike NFW, the DC14 halo parameters found in our rotation curve fits naturally fall within two standard deviations of the mass-concentration relation predicted by $\Lambda$CDM and the stellar mass-halo mass relation inferred from abundance matching with few outliers. Halo profiles modified by baryonic processes are therefore more consistent with expectations from $\Lambda$ cold dark matter ($\Lambda$CDM) cosmology and provide better fits to galaxy rotation curves across a wide range of galaxy properties than do halo models that neglect baryonic physics. Our results offer a solution to the decade long cusp-core discrepancy.

We address the origin of Ultra-Diffuse Galaxies (UDGs), which have stellar masses typical of dwarf galaxies but effective radii of Milky Way-sized objects. Their formation mechanism, and whether they are failed $\rm L_{\star}$ galaxies or diffuse dwarfs, are challenging issues. Using zoom-in cosmological simulations from the NIHAO project, we show that UDG analogues form naturally in medium-mass haloes due to episodes of gas outflows associated with star formation. The simulated UDGs live in isolated haloes of masses $10^{10-11}\rm M_{\odot}$, have stellar masses of $10^{7-8.5}\rm M_{\odot}$, effective radii larger than 1 kpc and dark matter cores. They show a broad range of colors, an average S\'ersic index of 0.83, a typical distribution of halo spin and concentration, and a non-negligible HI gas mass of $10^{7-9}\rm M_{\odot}$, which correlates with the extent of the galaxy. Gas availability is crucial to the internal processes that form UDGs: feedback driven gas outflows, and subsequent dark matter and stellar expansion, are the key to reproduce faint, yet unusually extended, galaxies. This scenario implies that UDGs represent a dwarf population of low surface brightness galaxies and should exist in the field. The largest isolated UDGs should contain more HI gas than less extended dwarfs of similar $\rm M_{\star}$.

Using idealised N-body simulations of a Milky Way-sized disc galaxy, we qualitatively study how the metallicity distributions of the thin disc star particles are modified by the formation of the bar and spiral arm structures. The thin disc in our numerical experiments initially has a tight negative radial metallicity gradient and a constant vertical scale-height. We show that the radial mixing of stars drives a positive vertical metallicity gradient in the thin disc. On the other hand, if the initial thin disc is flared, with vertical scale-height increasing with galactocentric radius, the metal poor stars originally in the outer disc become dominant in regions above the disc plane at every radii. This process can drive a negative vertical metallicity gradient, which is consistent with the current observed trend. This model mimics a scenario where the star-forming thin disc was flared in the outer region at earlier epochs. Our numerical experiment with an initial flared disc predicts that the negative vertical metallicity gradient of the mono-age relatively young thin disc population should be steeper in the inner disc, and the radial metallicity gradient of the mono-age population should be shallower at greater heights above the disc plane. We also predict that the metallicity distribution function of mono-age young thin disc populations above the disc plane would be more positively skewed in the inner disc compared to the outer disc.

We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ~10 changes in the enclosed dark matter mass at one per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: e_SF=(M_{star}/M_{200})/(\Omega_b/\Omega_m). In addition we report a new correlation with the compactness of the stellar system: e_R=r_{1/2}/R_{200}. We provide an analytic formula depending on e_SF and e_R for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes.

We compare the Baryonic Tully-Fisher relation (BTFR) of simulations and observations of galaxies ranging from dwarfs to spirals, using various measures of rotational velocity Vrot. We explore the BTFR when measuring Vrot at the flat part of the rotation curve, Vflat, at the extent of HI gas, Vlast, and using 20% (W20) and 50% (W50) of the width of HI line profiles. We also compare with the maximum circular velocity of the parent halo, Vmax, within dark matter only simulations. The different BTFRs increasingly diverge as galaxy mass decreases. Using Vlast one obtains a power law over four orders of magnitude in baryonic mass, with slope similar to the observed BTFR. Measuring Vflat gives similar results as Vlast when galaxies with rising rotation curves are excluded. However, higher rotation velocities would be found for low mass galaxies if the cold gas extended far enough for Vrot to reach a maximum. W20 gives a similar slope as Vlast but with slightly lower values of Vrot for low mass galaxies, although this may depend on the extent of the gas in your galaxy sample. W50 bends away from these other relations toward low velocities at low masses. By contrast, Vmax bends toward high velocities for low mass galaxies, as cold gas does not extend out to the radius at which halos reach Vmax. Our study highlights the need for careful comparisons between observations and models: one needs to be consistent about the particular method of measuring Vrot, and precise about the radius at which velocities are measured.

Using cosmological galaxy formation simulations from the MaGICC project, spanning more than three magnitudes in stellar mass (~10^7-3x10^{10} Msun), we trace the baryonic cycle of infalling gas from the virial radius through to its participation in the star formation process. An emphasis is placed upon the temporal history of chemical enrichment during its passage through the corona and CGM. We derive the distributions of time between gas crossing the virial radius and being accreted to the star forming region (which allows mixing within the corona), as well as the time between gas being accreted to the star forming region and then forming stars (which allows mixing within the disc). Significant numbers of stars are formed from gas that cycles back through the hot halo after first accreting to the star forming region. Gas entering high mass galaxies is pre-enriched in low mass proto-galaxies prior to entering the virial radius of the central progenitor, with only small amounts of primordial gas accreted, even at high redshift (z~5). After entering the virial radius, significant further enrichment occurs prior to the accretion of the gas to the star forming region, with gas that is feeding the star forming region surpassing 0.1Z by z=0. Mixing with halo gas, itself enriched via galactic fountains, is thus crucial in determining the metallicity at which gas is accreted to the disc. The lowest mass simulation (Mvir~2x10^{10}Msun, with M*~10^7Msun), by contrast, accretes primordial gas through the virial radius and onto the disc at all times. Much like the classical analytical solutions to the `G-dwarf problem', overproduction of low-metallicity stars is ameliorated by the inefficiency of star formation. Finally, gas outflow/metal removal rates from star forming regions as a function of galactic mass are presented.

Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulated spiral galaxies with masses comparable to the Milky Way. The simulations span a range of star formation and energy feedback strengths and prescriptions, particle- and grid-based hydrodynamical implementations, as well as initial conditions/assembly history. Results. Consistently, we find that the steeper, negative, radial metallicity gradients seen in the mid-plane flatten with increasing height away from the plane. In simulations with stronger (and/or more spatially-extended) feedback, the negative radial gradients invert, becoming positive for heights in excess of 1 kpc. Such behaviour is consistent with that inferred from recent observations. Our measurements of the vertical metallicity gradients show no clear correlation with galactocentric radius, and are in good agreement with those observed in the Milky Way's thick disc (locally). Conclusions. Simulations employing stronger/more extended feedback prescriptions possess radial and vertical metallicity and kinematic gradients more in line with recent observations. The inverted, positive, radial metallicity gradients seen in the simulated thick stellar discs originate from a population of younger, more metal-rich, stars formed in-situ, superimposed upon a background population of older migrators from the inner disc; the contrast provided by the former increases radially, due to the inside-out growth of the disc. A similar behaviour may be responsible for the same flattening seen in the radial gradients with scaleheight in the Milky Way.

Rotation curves of galaxies show a wide range of shapes, which can be paramaterized as scatter in Vrot(1kpc)/Vmax i.e.the ratio of the rotation velocity measured at 1kpc and the maximum measured rotation velocity. We examine whether the observed scatter can be accounted for by combining scatters in disc scale-lengths, the concentration-halo mass relation, and the M*-Mhalo relation. We use these scatters to create model galaxy populations; when housed within dark matter halos that have universal, NFW density profiles, the model does not match the lowest observed values of Vrot(1kpc)/Vmax and has too little scatter in Vrot(1kpc)/Vmax compared to observations. By contrast, a model using a mass dependent dark matter profile, where the inner slope is determined by the ratio of M*/Mhalo, produces galaxies with low values of Vrot(1kpc)/Vmax and a much larger scatter, both in agreement with observation. We conclude that the large observed scatter in Vrot(1kpc)/Vmax favours density profiles that are significantly affected by baryonic processes. Alternative dark matter core formation models such as SIDM may also account for the observed variation in rotation curve shapes, but these observations may provide important constraints in terms of core sizes, and whether they vary with halo mass and/or merger history.

Using 22 hydrodynamical simulated galaxies in a LCDM cosmological context we recover not only the observed baryonic Tully-Fisher relation, but also the observed "mass discrepancy--acceleration" relation, which reflects the distribution of the main components of the galaxies throughout their disks. This implies that the simulations, which span the range 52 < V$_{\rm flat}$ < 222 km/s where V$_{\rm flat}$ is the circular velocity at the flat part of the rotation curve, and match galaxy scaling relations, are able to recover the observed relations between the distributions of stars, gas and dark matter over the radial range for which we have observational rotation curve data. Furthermore, we explicitly match the observed baryonic to halo mass relation for the first time with simulated galaxies. We discuss our results in the context of the baryon cycle that is inherent in these simulations, and with regards to the effect of baryonic processes on the distribution of dark matter.

Dark matter cores within galaxy haloes can be formed by energy feedback from star forming regions: an energy balance suggests that the maximum core formation efficiency arises in galaxies with M$_{\star}\sim10^{8.5}$M$_{\odot}$. We show that a model population of galaxies, in which the density profile has been modified by such baryonic feedback, is able to explain the observed galaxy velocity function and Tully-Fisher relations significantly better than a model in which a universal cuspy density profile is assumed. Alternative models, namely warm or self-interacting dark matter, also provide a better match to these observed relations than a universal profile model does, but make different predictions for how halo density profiles vary with mass compared to the baryonic feedback case. We propose that different core formation mechanisms may be distinguished based on the imprint they leave on galaxy populations over a wide range of mass. Within the current observational data we find evidence of the expected signatures of the mass dependence of core formation generated by baryonic feedback.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

The velocity function derived from large scale surveys can be compared with the predictions of LCDM cosmology, by matching the measured rotation velocities Vrot of galaxies to the maximum circular velocity of dark matter (DM) halos Vmax. For Vrot<50km/s, a major discrepancy arises between the observed and LCDM velocity functions. However, the manner in which different observational measures of Vrot are associated with Vmax is not straight forward in dwarf galaxies. We instead relate galaxies to DM halos using the empirical baryon- mass to halo-mass relation, and show that different observational measures of Vrot result in very different velocity functions. We show how the W50 velocity function, i.e. using the HI profile line width at 50% of peak HI flux to measure Vrot, can be reconciled with a LCDM cosmology. Our semi-empirical methodology allows us to determine the region of rotation curves that are probed by HI measurements (RHI), and shows that the Vrot of dwarfs are generally measured at a fraction of Rmax, explaining their tendency to have rising rotation curves. We provide fitting formulae for relating RHI and Reff (the effective radius) to the virial radius of DM halos. To continue to use velocity functions as a probe of LCDM cosmology, it is necessary to be precise about how the different measures of rotation velocity are probing the mass of the DM halos, dropping the assumption that any measure of rotational velocity can be equally used as a proxy for Vmax.

Observed kinematical data of 40 Local Group (LG) members are used to derive the dark matter halo mass of such galaxies. Haloes are selected from the theoretically expected LG mass function and two different density profiles are assumed, a standard universal cuspy model and a mass dependent profile which accounts for the effects of baryons in modifying the dark matter distribution within galaxies. The resulting relations between stellar and halo mass are compared with expectations from abundance matching. Using a universal cuspy profile, the ensemble of LG galaxies is fit in relatively low mass haloes, leaving "dark" many massive haloes of \mhalo$\gtrsim$10$^{10}$\msun: this reflects the "too big to fail" problem and results in a \mstar-\mhalo\ relation that differs from abundance matching predictions. Moreover, the star formation efficiency of isolated LG galaxies increases with decreasing halo mass when adopting a cuspy model. By contrast, using the mass dependent density profile, dwarf galaxies with \mstar$\gtrsim$10$^{6}$\msun are assigned to more massive haloes, which have a central cored distribution of dark matter: the "too big to fail" problem is alleviated, the resultant \mstar-\mhalo\ relation follows abundance matching predictions down to the completeness limit of current surveys, and the star formation efficiency of isolated members decreases with decreasing halo mass, in agreement with theoretical expectations. Finally, the cusp/core space of LG galaxies is presented, providing a framework to understand the non-universality of their density profiles.

We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range of disk scaling relationships. We find that the best fit parameters of a generic double power-law density profile vary in a systematic manner that depends on the stellar-to-halo mass ratio of each galaxy. Thus, the quantity Mstar/Mhalo constrains the inner ($\gamma$) and outer ($\beta$) slopes of dark matter density, and the sharpness of transition between the slopes($\alpha$), reducing the number of free parameters of the model to two. Due to the tight relation between stellar mass and halo mass, either of these quantities is sufficient to describe the dark matter halo profile including the effects of baryons. The concentration of the haloes in the hydrodynamical simulations is consistent with N-body expectations up to Milky Way mass galaxies, at which mass the haloes become twice as concentrated as compared with pure dark matter runs. This mass dependent density profile can be directly applied to rotation curve data of observed galaxies and to semi analytic galaxy formation models as a significant improvement over the commonly used NFW profile.

Using cosmological galaxy simulations from the MaGICC project, we study the evolution of the stellar masses, star formation rates and gas phase abundances of star forming galaxies. We derive the stellar masses and star formation rates using observational relations based on spectral energy distributions by applying the new radiative transfer code GRASIL-3D to our simulated galaxies. The simulations match well the evolution of the stellar mass-halo mass relation, have a star forming main sequence that maintains a constant slope out to redshift z $\sim$ 2, and populate projections of the stellar mass - star formation - metallicity plane, similar to observed star forming disc galaxies. We discuss small differences between these projections in observational data and in simulations, and the possible causes for the discrepancies. The light-weighted stellar masses are in good agreement with the simulation values, the differences between the two varying between 0.06 dex and 0.20 dex. We also find a good agreement between the star formation rate tracer and the true (time-averaged) simulation star formation rates. Regardless if we use mass- or light-weighted quantities, our simulations indicate that bursty star formation cycles can account for the scatter in the star forming main sequence.

We present an analysis of the role of feedback in shaping the neutral hydrogen (HI) content of simulated disc galaxies. For our analysis, we have used two realisations of two separate Milky Way-like (~L*) discs - one employing a conservative feedback scheme (MUGS), the other significantly more energetic (MaGICC). To quantify the impact of these schemes, we generate zeroth moment (surface density) maps of the inferred HI distribution; construct power spectra associated with the underlying structure of the simulated cold ISM, in addition to their radial surface density and velocity dispersion profiles. Our results are compared with a parallel, self-consistent, analysis of empirical data from THINGS (The HI Nearby Galaxy Survey). Single power-law fits (P~k^gamma) to the power spectra of the stronger-feedback (MaGICC) runs (over spatial scales corresponding to 0.5 kpc to 20 kpc) result in slopes consistent with those seen in the THINGS sample (gamma = -2.5). The weaker-feedback (MUGS) runs exhibit shallower power law slopes (gamma = -1.2). The power spectra of the MaGICC simulations are more consistent though with a two-component fit, with a flatter distribution of power on larger scales (i.e., gamma = -1.4 for scales in excess of 2 kpc) and a steeper slope on scales below 1 kpc (gamma = -5), qualitatively consistent with empirical claims, as well as our earlier work on dwarf discs. The radial HI surface density profiles of the MaGICC discs show a clear exponential behaviour, while those of the MUGS suite are essentially flat; both behaviours are encountered in nature, although the THINGS sample is more consistent with our stronger (MaGICC) feedback runs.

We introduce a new model for the spectral energy distribution of galaxies, GRASIL-3D, which includes a careful modelling of the dust component of the interstellar medium. GRASIL-3D is an entirely new model based on the formalism of an existing and widely applied spectrophotometric model, GRASIL, but specifically designed to be interfaced with galaxies with any arbitrarily given geometry, such as galaxies calculated by theoretical hydrodynamical galaxy formation codes. GRASIL-3D is designed to separately treat radiative transfer in molecular clouds and in the diffuse cirrus component. The code has a general applicability to the outputs of simulated galaxies, either from Lagrangian or Eulerian hydrodynamic codes. As an application, the new model has been interfaced to the P-DEVA and GASOLINE smoothed-particle hydrodynamic codes, and has been used to calculate the spectral energy distribution for a variety of simulated galaxies from UV to sub-millimeter wavelengths, whose comparison with observational data gives encouraging results. In addition, GRASIL-3D allows 2D images of such galaxies to be obtained, at several angles and in different bands.

We contend that a single power law halo mass distribution is appropriate for direct matching to the stellar masses of observed Local Group dwarf galaxies, allowing the determination of the slope of the stellar mass-halo mass relation for low mass galaxies. Errors in halo masses are well defined as the Poisson noise of simulated local group realisations, which we determine using constrained local universe simulations (CLUES). For the stellar mass range 10$^7$<M*<10$^8$M$_\odot$, for which we likely have a complete census of observed galaxies, we find that the stellar mass-halo mass relation follows a power law with slope of 3.1, significantly steeper than most values in the literature. The steep relation between stellar and halo masses indicates that Local Group dwarf galaxies are hosted by dark matter halos with a small range of mass. Our methodology is robust down to the stellar mass to which the census of observed Local Group galaxies is complete, but the significant uncertainty in the currently measured slope of the stellar-to halo mass relation will decrease dramatically if the Local Group completeness limit was $10^{6.5}$M$\odot$ or below, highlighting the importance of pushing such limit to lower masses and larger volumes.

We study the effect of warm dark matter (WDM) on hydrodynamic simulations of galaxy formation as part of the Making Galaxies in a Cosmological Context (MaGICC) project. We simulate three different galaxies using three WDM candidates of 1, 2 and 5 keV and compare results with pure cold dark matter simulations. WDM slightly reduces star formation and produces less centrally concentrated stellar profiles. These effects are most evident for the 1 keV candidate but almost disappear for $m_{\mathrm{WDM}}>2$ keV. All simulations form similar stellar discs independent of WDM particle mass. In particular, the disc scale length does not change when WDM is considered. The reduced amount of star formation in the case of 1 keV particles is due to the effects of WDM on merging satellites which are on average less concentrated and less gas rich. The altered satellites cause a reduced starburst during mergers because they trigger weaker disc instabilities in the main galaxy. Nevertheless we show that disc galaxy evolution is much more sensitive to stellar feedback than it is to WDM candidate mass. Overall we find that WDM, especially when restricted to current observational constraints ($m_{\mathrm{WDM}}>2$ keV), has a minor impact on disc galaxy formation.

We present a cosmological hydrodynamical simulation of a representative volume of the Universe, as part of the Making Galaxies in a Cosmological Context (MaGICC) project. MaGICC uses a thermal implementation for supernova and early stellar feedback. This work tests the feedback model at lower resolution across a range of galaxy masses, morphologies and merger histories. The simulated sample compares well with observations of high redshift galaxies ($z \ge 2$) including the stellar mass - halo mass ($M_\star - M_h$ ) relation, the Galaxy Stellar Mass Function (GSMF) at low masses ($M_\star \lt 5 \times 10^{10} M_\odot$ ) and the number density evolution of low mass galaxies. The poor match of $M_\star - M_h$ and the GSMF at high masses ($M_\star \ge 5 \times 10^{10} M_\odot$ ) indicates supernova feedback is insufficient to limit star formation in these haloes. At $z = 0$, our model produces too many stars in massive galaxies and slightly underpredicts the stellar mass around $L_\star$ mass galaxy. Altogether our results suggest that early stellar feedback, in conjunction with supernovae feedback, plays a major role in regulating the properties of low mass galaxies at high redshift.

We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9<Mhalo/Msun<7.8e11, hosting galaxies with stellar masses: 5.0e5<M*/Msun<8.3e10, i.e. from dwarf to L*. The galaxies are simulated with several baryonic prescriptions, including a range of stellar feedbacks. The main result is a clear dependence of the inner slope of the DM density profile, \alpha\ in \rho r^\alpha, on the ratio between stellar-to-halo mass (M*/Mhalo). This relation is independent of the stellar feedback scheme, allowing a prediction for cusp vs core formation. When M*/Mhalo is low, ~0.01%, energy from stellar feedback is insufficient to significantly alter the inner DM density and the galaxy retains a cuspy profile. At higher M*/Mhalo, feedback drives the expansion of the DM and generates cored profiles. The flattest profiles form where M*/Mhalo~0.5%. Above this ratio, stars formed in the central regions deepen the gravitational potential enough to oppose this supernova-driven expansion process, resulting in smaller cores and cuspier profiles. Combining the dependence of \alpha\ on M*/Mhalo with the abundance matching relation between M* and Mhalo provides a prediction for how \alpha\ varies with M*. Further, using the Tully-Fisher relation allows a prediction for the dependence of the DM inner slope on the observed rotation velocity of galaxies. The most cored galaxies are expected to have Vrot~50km/s, with \alpha\ decreasing for more massive disc galaxies: spirals with Vrot~150km/s have central slopes \alpha<-0.8, approaching the NFW profile. This novel prediction for the dependence of \alpha\ on disc galaxy mass can be tested using current observational data sets, and can be applied to theoretical modeling of mass profiles and populations of disc galaxies.

We analyse the structure and chemical enrichment of a Milky Way-like galaxy with a stellar mass of 2 10^{10} M_sun, formed in a cosmological hydrodynamical simulation. It is disk-dominated with a flat rotation curve, and has a disk scale length similar to the Milky Way's, but a velocity dispersion that is ~50% higher. Examining stars in narrow [Fe/H] and [\alpha/Fe] abundance ranges, we find remarkable qualitative agreement between this simulation and observations: a) The old stars lie in a thickened distribution with a short scale length, while the young stars form a thinner disk, with scale lengths decreasing, as [Fe/H] increases. b) Consequently, there is a distinct outward metallicity gradient. c) Mono-abundance populations exist with a continuous distribution of scale heights (from thin to thick). However, the simulated galaxy has a distinct and substantive very thick disk (h_z~1.5 kpc), not seen in the Milky Way. The broad agreement between simulations and observations allows us to test the validity of observational proxies used in the literature: we find in the simulation that mono-abundance populations are good proxies for single age populations (<1 Gyr) for most abundances.

Using a cosmological hydrodynamical simulation of a galaxy of similar mass to the Large Magellanic Cloud (LMC), we examine the predicted characteristics of its lowest metallicity populations. In particular, we emphasise the spatial distributions of first (Pop III) and second (polluted by only immediate Pop III ancestors) generation stars. We find that primordial composition stars form not only in the central galaxy's progenitor, but also in locally collapsed sub-halos during the early phases of galaxy formation. The lowest metallicity stars in these sub-halos end up in a relatively extended distribution around the host, with these accreted stars possessing present-day galactocentric distances as great as ~40kpc. By contrast, the earliest stars formed within the central galaxy remain in the inner region, where the vast majority of star formation occurs, for the entirety of the simulation. Consequently, the fraction of stars that are from the earliest generation increases strongly with radius.

We examine the role of energy feedback in shaping the distribution of metals within cosmological hydrodynamical simulations of L* disc galaxies. While negative abundance gradients today provide a boundary condition for galaxy evolution models, in support of inside-out disc growth, empirical evidence as to whether abundance gradients steepen or flatten with time remains highly contradictory. We made use of a suite of L* discs, realised with and without `enhanced' feedback. All the simulations were produced using the smoothed particle hydrodynamics code Gasoline, and their in situ gas-phase metallicity gradients traced from redshift z~2 to the present-day. Present-day age-metallicity relations and metallicity distribution functions were derived for each system. The `enhanced' feedback models, which have been shown to be in agreement with a broad range of empirical scaling relations, distribute energy and re-cycled ISM material over large scales and predict the existence of relatively `flat' and temporally invariant abundance gradients. Enhanced feedback schemes reduce significantly the scatter in the local stellar age-metallicity relation and, especially, the [O/Fe]-[Fe/H] relation. The local [O/Fe] distribution functions for our L* discs show clear bimodality, with peaks at [O/Fe]=-0.05 and +0.05 (for stars with [Fe/H]>-1), consistent with our earlier work on dwarf discs. Our results with `enhanced' feedback are inconsistent with our earlier generation of simulations realised with `conservative' feedback. We conclude that spatially-resolved metallicity distributions, particularly at high-redshift, offer a unique and under-utilised constraint on the uncertain nature of stellar feedback processes.

We use dark matter only and full hydrodynamical Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group to study the density profile of subhaloes of the simulated Milky Way and Andromeda galaxies. We show that the Einasto model provides the best description of the subhaloes' density profile, as opposed to the more commonly used NFW profile or any generalisation of it. We further find that the Einasto shape parameter \nEin\ is strongly correlated with the total subhalo mass, pointing towards the notion of a non-universality of the subhaloes' density profile. We observe that the effect of mass loss due to tidal stripping, in both the dark matter only and the hydrodynamical run, is the reduction of the shape parameter \nEin\ between the infall and the present time. Assuming now that the dSphs of our Galaxy follow the Einasto profile and using the maximum and minimum values of \nEin\ from our hydrodynamical simulation as a gauge, we can improve the observational constraints on the \Rmax-\Vmax\ pairs obtained for the brightest satellite galaxies of the Milky Way. When considering only the subhaloes with $-13.2\lesssim M_V\lesssim-8.8$, i.e. the range of luminosity of the classical dwarfs, we find that all our simulated objects are consistent with the observed dSphs if their haloes follow the Einasto model with $1.6\lesssim n_{\rm E} \lesssim5.3$. The numerically motivated Einasto profile for the observed dSphs will alleviate the recently presented "massive failures" problem.

We analyze and compare the bulges of a sample of L* spiral galaxies in hydrodynamical simulations in a cosmological context, using two different codes, P-DEVA and GASOLINE. The codes regulate star formation in very different ways, with P-DEVA simulations inputing low star formation efficiency under the assumption that feedback occurs on subgrid scales, while the GASOLINE simulations have feedback which drives large scale outflows. In all cases, the marked knee-shape in mass aggregation tracks, corresponding to the transition from an early phase of rapid mass assembly to a later slower one, separates the properties of two populations within the simulated bulges. The bulges analyzed show an important early starburst resulting from the collapse-like fast phase of mass assembly, followed by a second phase with lower star formation, driven by a variety of processes such as disk instabilities and/or mergers. Classifying bulge stellar particles identified at z=0 into old and young according to these two phases, we found bulge stellar sub-populations with distinct kinematics, shapes, stellar ages and metal contents. The young components are more oblate, generally smaller, more rotationally supported, with higher metallicity and less alpha-element enhanced than the old ones. These results are consistent with the current observational status of bulges, and provide an explanation for some apparently paradoxical observations, such as bulge rejuvenation and metal-content gradients observed. Our results suggest that bulges of L* galaxies will generically have two bulge populations which can be likened to classical and pseudo-bulges, with differences being in the relative proportions of the two, which may vary due to galaxy mass and specific mass accretion and merger histories.

We introduce the Making Galaxies in a Cosmological Context (MaGICC) program of smoothed particle hydrodynamics (SPH) simulations. We describe a parameter study of galaxy formation simulations of an L* galaxy that uses early stellar feedback combined with supernova feedback to match the stellar mass--halo mass relationship. While supernova feedback alone can reduce star formation enough to match the stellar mass--halo mass relationship, the galaxy forms too many stars before z=2 to match the evolution seen using abundance matching. Our early stellar feedback is purely thermal and thus operates like a UV ionization source as well as providing some additional pressure from the radiation of massive, young stars. The early feedback heats gas to >10^6 K before cooling to 10^4 K. The pressure from this hot gas creates a more extended disk and prevents more star formation prior to z=1 than supernovae feedback alone. The resulting disk galaxy has a flat rotation curve, an exponential surface brightness profile, and matches a wide range of disk scaling relationships. The disk forms from the inside-out with an increasing exponential scale length as the galaxy evolves. Overall, early stellar feedback helps to simulate galaxies that match observational results at low and high redshifts.

We test the validity of comparing simulated field disk galaxies with the empirical properties of systems situated within environments more comparable to loose groups, including the Milky Way's Local Group. Cosmological simulations of Milky Way-mass galaxies have been realised in two different environment samples: in the field and in environments with similar properties to the Local Group. Apart from the environments of the galaxies, the samples are kept as homogeneous as possible with equivalent ranges in last major merger time, halo mass and halo spin. Comparison of these two samples allow for systematic differences in the simulations to be identified. Metallicity gradients, disk scale lengths, colours, magnitudes and age-velocity dispersion relations are studied for each galaxy in the suite and the strength of the link between these and environment of the galaxies is studied. The bulge-to-disk ratio of the galaxies show that these galaxies are less spheroid dominated than many other simulated galaxies in literature with the majority of both samples being disk dominated. We find that secular evolution and mergers dominate the spread of morphologies and metallicity gradients with no visible differences between the two environment samples. In contrast with this consistency in the two samples there is tentative evidence for a systematic difference in the velocity dispersion-age relations of galaxies in the different environments. Loose group galaxies appear to have more discrete steps in their velocity dispersion-age relations. We conclude that at the current resolution of cosmological galaxy simulations field environment galaxies are sufficiently similar to those in loose groups to be acceptable proxies for comparison with the Milky Way provided that a similar assembly history is considered.

We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ~ 0.1 and 1 L* generated using the smooth particle hydrodynamic code GASOLINE. These simulations are part of the Making Galaxies In a Cosmological Context (MAGICC) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a 'lower feedback' (LF) model which has strength comparable to other implementations in the literature. The 'MAGICC feedback' (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O VI absorption. The simulations do not yet address the question of what happens to the O VI when the galaxies stop forming stars. Our models also predict a reservoir of cool H I clouds that show strong Ly\alpha absorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O VI and H I gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H I gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.

By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born "in situ" relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.

We examine the properties and evolution of a simulated polar disc galaxy. This galaxy is comprised of two orthogonal discs, one of which contains old stars (old stellar disc), and the other, containing both younger stars and the cold gas (polar disc) of the galaxy. By exploring the shape of the inner region of the dark matter halo, we are able to confirm that the halo shape is a oblate ellipsoid flattened in the direction of the polar disc. We also note that there is a twist in the shape profile, where the innermost 3 kpc of the halo flattens in the direction perpendicular to the old disc, and then aligns with the polar disc out until the virial radius. This result is then compared to the halo shape inferred from the circular velocities of the two discs. We also use the temporal information of the simulation to track the system's evolution, and identify the processes which give rise to this unusual galaxy type. We confirm the proposal that the polar disc galaxy is the result of the last major merger, where the angular moment of the interaction is orthogonal to the angle of the infalling gas. This merger is followed by the resumption of coherent gas infall. We emphasise that the disc is rapidly restored after the major merger and that after this event the galaxy begins to tilt. A significant proportion of the infalling gas comes from filaments. This infalling gas from the filament gives the gas its angular momentum, and, in the case of the polar disc galaxy, the direction of the gas filament does not change before or after the last major merger.

Within a cosmological hydrodynamical simulation, we form a disc galaxy with sub- components which can be assigned to a thin stellar disc, thick disk, and a low mass stellar halo via a chemical decomposition. The thin and thick disc populations so selected are distinct in their ages, kinematics, and metallicities. Thin disc stars are young (<6.6 Gyr), possess low velocity dispersion ({\sigma}U,V,W = 41, 31, 25 km/s), high [Fe/H], and low [O/Fe]. The thick disc stars are old (6.6<age<9.8 Gyrs), lag the thin disc by \sim21 km/s, possess higher velocity dispersion ({\sigma}U,V,W = 49, 44, 35 km/s), relatively low [Fe/H] and high [O/Fe]. The halo component comprises less than 4% of stars in the "solar annulus" of the simulation, has low metallicity, a velocity ellipsoid defined by ({\sigma}U,V,W = 62, 46, 45 km/s) and is formed primarily in-situ during an early merger epoch. Gas-rich mergers during this epoch play a major role in fuelling the formation of the old disc stars (the thick disc). This is consistent with studies which show that cold accretion is the main source of a disc galaxy's baryons. Our simulation initially forms a relatively short (scalelength \sim1.7 kpc at z=1) and kinematically hot disc, primarily from gas accreted during the galaxy's merger epoch. Far from being a competing formation scenario, migration is crucial for reconciling the short, hot, discs which form at high redshift in {\Lambda}CDM, with the properties of the thick disc at z=0. The thick disc, as defined by its abundances maintains its relatively short scale-length at z = 0 (2.31 kpc) compared with the total disc scale-length of 2.73 kpc. The inside-out nature of disc growth is imprinted the evolution of abundances such that the metal poor {\alpha}-young population has a larger scale-length (4.07 kpc) than the more chemically evolved metal rich {\alpha}-young population (2.74 kpc).

We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass -- halo mass relationship. Each realisation employed identical initial conditions and assembly histories, but differed from one another in their underlying baryonic physics prescriptions, including (a) the efficiency with which each supernova energy couples to the ISM, (b) the impact of feedback associated with massive star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling of Type II SNe remnants, (d) the treatment of metal diffusion, and (e) varying the IMF. Our analysis focusses on the resulting stellar metallicity distribution functions (MDFs) in each simulated (analogous) `solar neighbourhood' and central `bulge' region. We compare the simulated MDFs' skewness, kurtosis, and dispersion (inter-quartile, inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empirical solar neighbourhood MDF and Local Group dwarfs. We find that the MDFs of the simulated discs are more negatively skewed, with higher kurtosis, than those observed locally. We can trace this difference to the simulations' tight and correlated age-metallicity relations (compared with that of the Milky Way), suggesting that these relations within `dwarf' discs might be steeper than in L* discs and/or the degree of stellar orbital re-distribution and migration inferred locally has not been captured in their entirety, at the resolution of our simulations. The important role of metal diffusion in ameliorating the over-production of extremely metal-poor stars is highlighted.

We use the same physical model to simulate four galaxies that match the relation between stellar and total mass, over a mass range that includes the vast majority of disc galaxies. The resultant galaxies, part of the Making Galaxies in a Cosmological Context (MaGICC) program, also match observed relations between luminosity, rotation velocity, size, colour, star formation rate, HI mass, baryonic mass, and metallicity. Radiation from massive stars and supernova energy regulate star formation and drive outflows, balancing the complex interplay between cooling gas, star formation, large scale outflows, and recycling of gas in a manner which correctly scales with the mass of the galaxy. Outflows also play a key role in simulating galaxies with exponential surface brightness profiles, flat rotation curves and dark matter cores. Our study implies that large scale outflows are the primary driver of the dependence of disc galaxy properties on mass. We show that the amount of outflows invoked in our model is required to meet the constraints provided by observations of OVI absorption lines in the circum-galactic-media of local galaxies.

We examine radial and vertical metallicity gradients using a suite of disk galaxy simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient and reconcile differences between extant models and observations within the `inside-out' disk growth paradigm. A sample of 25 disks is used, consisting of 19 from our RaDES (Ramses Disk Environment Study) sample, realised with the adaptive mesh refinement code RAMSES. Four disks are selected from the MUGS (McMaster Unbiased Galaxy Simulations) sample, generated with the smoothed particle hydrodynamics (SPH) code GASOLINE, alongside disks from Rahimi et al. (GCD+) and Kobayashi & Nakasato (GRAPE-SPH). Two chemical evolution models of inside-out disk growth were employed to contrast the temporal evolution of their radial gradients with those of the simulations. We find that systematic differences exist between the predicted evolution of radial abundance gradients in the RaDES and chemical evolution models, compared with the MUGS sample; specifically, the MUGS simulations are systematically steeper at high-redshift, and present much more rapid evolution in their gradients. We find that the majority of the models predict radial gradients today which are consistent with those observed in late-type disks, but they evolve to this self-similarity in different fashions, despite each adhering to classical `inside-out' growth. We find that radial dependence of the efficiency with which stars form as a function of time drives the differences seen in the gradients; systematic differences in the sub-grid physics between the various codes are responsible for setting these gradients. Recent, albeit limited, data at redshift z=1.5 are consistent with the steeper gradients seen in our SPH sample, suggesting a modest revision of the classical chemical evolution models may be required.

A clear prediction of the Cold Dark Matter model is the existence of cuspy dark matter halo density profiles on all mass scales. This is not in agreement with the observed rotation curves of spiral galaxies, challenging on small scales the otherwise successful CDM paradigm. In this work we employ high resolution cosmological hydro-dynamical simulations to study the effects of dissipative processes on the inner distribution of dark matter in Milky-Way like objects (M~1e12 Msun). Our simulations include supernova feedback, and the effects of the radiation pressure of massive stars before they explode as supernovae. The increased stellar feedback results in the expansion of the dark matter halo instead of contraction with respect to N-body simulations. Baryons are able to erase the dark matter cuspy distribution creating a flat, cored, dark matter density profile in the central several kpc of a massive Milky-Way like halo. The profile is well fit by a Burkert profile, with fitting parameters consistent with the observations. In addition, we obtain flat rotation curves as well as extended, exponential stellar disk profiles. While the stellar disk we obtain is still partially too thick to resemble the MW thin disk, this pilot study shows that there is enough energy available in the baryonic component to alter the dark matter distribution even in massive disc galaxies, providing a possible solution to the long standing problem of cusps vs. cores.

Encoded within the morphological structure of galaxies are clues related to their formation and evolutionary history. Recent advances pertaining to the statistics of galaxy morphology include sophisticated measures of concentration (C), asymmetry (A), and clumpiness (S). In this study, these three parameters (CAS) have been applied to a suite of simulated galaxies and compared with observational results inferred from a sample of nearby galaxies. The simulations span a range of late-type systems, with masses between ~1e10 Msun and ~1e12 Msun, and employ star formation density thresholds between 0.1 cm^-3 and 100 cm^-3. We have found that the simulated galaxies possess comparable concentrations to their real counterparts. However, the results of the CAS analysis revealed that the simulated galaxies are generally more asymmetric, and that the range of clumpiness values extends beyond the range of those observed. Strong correlations were obtained between the three CAS parameters and colour (B-V), consistent with observed galaxies. Furthermore, the simulated galaxies possess strong links between their CAS parameters and Hubble type, mostly in-line with their real counterparts.

Within a fully cosmological hydrodynamical simulation, we form a galaxy which rotates at 140 km/s, and is characterised by two loose spiral arms and a bar, indicative of a Hubble Type SBc/d galaxy. We show that our simulated galaxy has no classical bulge, with a pure disc profile at z=1, well after the major merging activity has ended. A long-lived bar subsequently forms, resulting in the formation of a secularly-formed "pseudo" bulge, with the final bulge-to-total light ratio B/T=0.21. We show that the majority of gas which loses angular momentum and falls to the central region of the galaxy during the merging epoch is blown back into the hot halo, with much of it returning later to form stars in the disc. We propose that this mechanism of redistribution of angular momentum via a galactic fountain, when coupled with the results from our previous study which showed why gas outflows are biased to have low angular momentum, can solve the angular momentum/bulgeless disc problem of the cold dark matter paradigm.

Using a suite of simulations (Governato et al. 2010) which successfully produce bulgeless (dwarf) disk galaxies, we provide an analysis of their associated cold interstellar media (ISM) and stellar chemical abundance patterns. A preliminary comparison with observations is undertaken, in order to assess whether the properties of the cold gas and chemistry of the stellar components are recovered successfully. To this end, we have extracted the radial and vertical gas density profiles, neutral hydrogen velocity dispersion, and the power spectrum of structure within the ISM. We complement this analysis of the cold gas with a brief examination of the simulations' metallicity distribution functions and the distribution of alpha-elements-to-iron.

We present an analysis of the neutral hydrogen (HI) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al. (2010), the high resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the HI distribution and kinematics of this simulated bulgeless disk with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g., velocity ellipsoid, turbulence), and the power spectrum of structure within the cold interstellar medium from the simulations. The highest resolution dwarf, when using a high density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the SMC. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within pc-scale molecular clouds.

We make a direct comparison of the derived dark matter (DM) distributions between hydrodynamical simulations of dwarf galaxies assuming a LCDM cosmology and the observed dwarf galaxies sample from the THINGS survey in terms of (1) the rotation curve shape and (2) the logarithmic inner density slope alpha of mass density profiles. The simulations, which include the effect of baryonic feedback processes, such as gas cooling, star formation, cosmic UV background heating and most importantly physically motivated gas outflows driven by supernovae (SNe), form bulgeless galaxies with DM cores. We show that the stellar and baryonic mass is similar to that inferred from photometric and kinematic methods for galaxies of similar circular velocity. Analyzing the simulations in exactly the same way as the observational sample allows us to address directly the so-called "cusp/core" problem in the LCDM model. We show that the rotation curves of the simulated dwarf galaxies rise less steeply than CDM rotation curves and are consistent with those of the THINGS dwarf galaxies. The mean value of the logarithmic inner density slopes alpha of the simulated galaxies' dark matter density profiles is ~ -0.4 +- 0.1, which shows good agreement with \alpha = -0.29 +- 0.07 of the THINGS dwarf galaxies. The effect of non-circular motions is not significant enough to affect the results. This confirms that the baryonic feedback processes included in the simulations are efficiently able to make the initial cusps with \alpha ~ -1.0 to -1.5 predicted by dark-matter-only simulations shallower, and induce DM halos with a central mass distribution similar to that observed in nearby dwarf galaxies.

Semi-analytic models (SAMs) are currently one of the primary tools with which to model statistically significant ensembles of galaxies. The underlying physical prescriptions inherent to each SAM are, in many cases, different from one another. Several SAMs have been applied to the dark matter merger trees extracted from the Millennium Run, including those associated with the well-known Munich and Durham lineages. We compare the predicted luminosity distributions of galaxy groups using four publicly available SAMs (De Lucia et al. 2006; Bower et al. 2006; Bertone et al. 2007; Font et al. 2008), in order to explore a galactic environment in which the models have not been explored to the same degree as they have in the field or in rich clusters. We identify a characteristic "wiggle" in the group galaxy luminosity function generated using the De Lucia et al. (2006) SAM, that is not present in the Durham-based models, consistent to some degree with observations. However, a comparison between conditional luminosity functions of groups between the models and observations of Yang et al. (2007) suggest that neither model is a particularly good match. The luminosity function wiggle is interpreted as the result of the two-mode AGN feedback implementation used in the Munich models, which itself results in flattened magnitude gap distribution. An associated analysis of the magnitude gap distribution between first- and second-ranked group galaxies shows that while the Durham models yield distributions with approximately equal luminosity first- and second-ranked galaxies, in agreement with observations, the De Lucia et al. models favours the scenario in which the second-ranked galaxy is approximately one magnitude fainter than the primary,especially when the dynamic range of the mock data is limited to 3 magnitudes.

We present the analysis of a suite of simulations run with different particle-and grid-based cosmological hydrodynamical codes and compare them with observational data of the Milky Way. This is the first study to make comparisons of properties of galaxies simulated with particle and grid-based codes. Our analysis indicates that there is broad agreement between these different modelling techniques. We study the velocity dispersion - age relation for disc stars at z=0 and find that four of the simulations are more consistent with observations by Holmberg et al. (2008) in which the stellar disc appears to undergo continual/secular heating. Two other simulations are in better agreement with the Quillen & Garnett (2001) observations that suggest a "saturation" in the heating profile for young stars in the disc. None of the simulations have thin discs as old as that of the Milky Way. We also analyse the kinematics of disc stars at the time of their birth for different epochs in the galaxies' evolution and find that in some simulations old stars are born cold within the disc and are subsequently heated, while other simulations possess old stellar populations which are born relatively hot. The models which are in better agreement with observations of the Milky Way's stellar disc undergo significantly lower minor-merger/assembly activity after the last major merger - i.e. once the disc has formed. All of the simulations are significantly "hotter" than the Milky Way disc; on top of the effects of mergers, we find a "floor" in the dispersion that is related to the underlying treatment of the heating and cooling of the interstellar medium, and the low density threshold which such codes use for star formation. This finding has important implications for all studies of disc heating that use hydrodynamical codes.

We analyse for the first time the radial abundance gradients of the disc stars of a disc galaxy simulated with our three dimensional, fully cosmological chemodynamical galaxy evolution code GCD+. We study how [Fe/H], [N/O], [O/Fe], [Mg/Fe] and [Si/Fe] vary with galactocentric radius. For the young stars of the disc, we found a negative slope for [Fe/H] and [N/O] but a positive [O/Fe], [Mg/Fe] and [Si/Fe] slope with radius. By analysing the star formation rate (SFR) at different radii, we found that the simulated disc contains a greater fraction of young stars in the outer regions, while the old stars tend to be concentrated in the inner parts of the disc. This can explain the positive [alpha/Fe] gradient as well as the negative [N/O] gradient with radius. This radial trend is a natural outcome of an inside-out formation of the disc, regardless of its size and can thus explain the recently observed positive [alpha/Fe] gradients in the Milky Way disc open clusters.

Using high resolution, fully cosmological smoothed particle hydro-dynamical simulations of dwarf galaxies in a Lambda cold dark matter Universe, we show how baryons attain a final angular momentum distribution which allows pure disc galaxies to form. Blowing out substantial amounts of gas through supernovae and stellar winds, which is well supported observationally, is a key ingredient in forming bulgeless discs. We outline why galactic outflows preferentially remove low angular momentum material, and show that this is a natural result when structure forms in a cold dark matter cosmology. The driving factors are a) the mean angular momentum of accreted material increases with time, b) lower potentials at early times, c) the existence of an extended reservoir of high angular momentum gas which is not within star forming regions, meaning that only gas from the inner region (low angular momentum gas) is expelled and d) the tendency for outflows to follow the path of least resistance which is perpendicular to the disc. We also show that outflows are enhanced during mergers, thus expelling much of the gas which has lost its angular momentum during these events, and preventing the formation of "classical", merger driven bulges in low mass systems. Stars formed prior to such mergers form a diffuse, extended stellar halo component.

We performed numerical simulations of mergers between gas-rich disc galaxies, which result in the formation of late-type galaxies. Stars formed during the merger end up in a thick disc that is partially supported by velocity dispersion and has high [alpha/Fe] ratios at all metallicities. Stars formed later end up in a thin, rotationally supported disc which has lower [alpha/Fe] ratios. While the structural and kinematical properties of the merger remnants depend strongly upon the orbital parameters of the mergers, we find a clear chemical signature of gas-rich mergers.

A sample of very high resolution cosmological disk galaxy simulations is used to investigate the evolution of galaxy disk sizes back to redshift 1 within the Lambda CDM cosmology. Artificial images in the rest frame B band are generated, allowing for a measurement of disk scale lengths using surface brightness profiles as observations would, and avoiding any assumption that light must follow mass as previous models have assumed. We demonstrate that these simulated disks are an excellent match to the observed magnitude - size relation for both local disks, and for disks at z=1 in the magnitude/mass range of overlap. We disentangle the evolution seen in the population as a whole from the evolution of individual disk galaxies. In agreement with observations, our simulated disks undergo roughly 1.5 magnitudes/arcsec^2 of surface brightness dimming since z=1. We find evidence that evolution in the magnitude - size plane varies by mass, such that galaxies with M* > 10^9 M_sun undergo more evolution in size than luminosity, while dwarf galaxies tend to evolve potentially more in luminosity. The disks grow in such a way as to stay on roughly the same stellar mass - size relation with time. Finally, due to an evolving stellar mass - SFR relation, a galaxy at a given stellar mass (or size) at z=1 will reside in a more massive halo and have a higher SFR, and thus a higher luminosity, than a counterpart of the same stellar mass at z=0.

We present a detailed analysis of a disc galaxy forming in a high-resolution fully cosmological simulation to investigate the nature of the outer regions of discs and their relevance for the disc formation process. Specifically, we focus on the phenomenon of misaligned disc components and find that the outer disc warp is a consequence of the misalignment between the inner disc and the surrounding hot gaseous halo. As the infalling cold gas sinks toward the centre of the galaxy, it is strongly torqued by the hot gas halo. By the time the fresh gas reaches the central disc-forming region its angular momentum is completely aligned with the spin of the hot gas halo. If the spin of the hot gas halo, in turn, is not aligned with that of the inner disc, a misaligned outer disc forms comprised of newly accreted material. The inner and outer components are misaligned with each other because they respond differently to infalling substructure and accretion. The warped disc feeds the main gas disc due to viscous angular momentum losses, but small amounts of star formation in the warp itself form a low-metallicity thick disc. We show that observations of resolved stellar populations in warped galaxies in the local universe could provide evidence for the presence of these processes and therefore indirectly reveal ongoing gas accretion and the existence of hot gas halos.

We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$\times10^{11}$ M$_\odot$ to 2$\times10^{12}$ M$_\odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.

For almost two decades the properties of "dwarf" galaxies have challenged the Cold Dark Matter (CDM) paradigm of galaxy formation. Most observed dwarf galaxies consists of a rotating stellar disc embedded in a massive DM halo with a near constant-density core. Yet, models based on the CDM scenario invariably form galaxies with dense spheroidal stellar "bulges" and steep central DM profiles, as low angular momentum baryons and DM sink to the center of galaxies through accretion and repeated mergers. Processes that decrease the central density of CDM halos have been identified, but have not yet reconciled theory with observations of present day dwarfs. This failure is potentially catastrophic for the CDM model, possibly requiring a different DM particle candidate. This Letter presents new hydrodynamical simulations in a Lambda$CDM framework where analogues of dwarf galaxies, bulgeless and with a shallow central DM profile, are formed. This is achieved by resolving the inhomogeneous interstellar medium, resulting in strong outflows from supernovae explosions which remove low angular momentum gas. This inhibits the formation of bulges and decreases the dark-matter density to less than half within the central kiloparsec. Realistic dwarf galaxies are thus shown to be a natural outcome of galaxy formation in the CDM scenario.

We used an N-body smoothed particle hydrodynamics algorithm, with a detailed treatment of star formation, supernovae feedback, and chemical enrichment, to perform eight simulations of mergers between gas-rich disc galaxies. We vary the mass ratio of the progenitors, their rotation axes, and their orbital parameters and analyze the kinematic, structural, and chemical properties of the remnants. Six of these simulations result in the formation of a merger remnant with a disc morphology as a result of the large gas-fraction of the remnants. We show that stars formed during the merger (a sudden starburst occur in our simulation and last for 0.2-0.3 Gyr) and those formed after the merger have different kinematical and chemical properties. The first ones are located in thick disc or the halo. They are partially supported by velocity dispersion and have high [alpha/Fe] ratios even at metallicities as high as [Fe/H]=-0.5. The former ones -- the young component -- are located in a thin disc rotationally supported and have lower [alpha/Fe] ratios. The difference in the rotational support of both components results in the rotation of the thick disc lagging that of the thin disc by as much as a factor of two, as recently observed.We find that, while the kinematic and structural properties of the merger remnant depends strongly upon the orbital parameters of the mergers, there is a remarkable uniformity in the chemical properties of the mergers. This suggests that general conclusions about the chemical signature of gas-rich mergers can be drawn.

We analyse the kinematics and chemistry of the bulge stars of two simulated disc galaxies using our chemodynamical galaxy evolution code GCD+. First we compare stars that are born inside the galaxy with those that are born outside the galaxy and are accreted into the centre of the galaxy. Stars that originate outside of the bulge are accreted into it early in its formation within 3 Gyrs so that these stars have high [alpha/Fe] as well as having a high total energy reflecting their accretion to the centre of the galaxy. Therefore, higher total energy is a good indicator for finding accreted stars. The bulges of the simulated galaxies formed through multiple mergers separated by about a Gyr. Since [alpha/Fe] is sensitive to the first few Gyrs of star formation history, stars that formed during mergers at different epochs show different [alpha/Fe]. We show that the [Mg/Fe] against star formation time relation can be very useful to identify a multiple merger bulge formation scenario, provided there is sufficiently good age information available. Our simulations also show that stars formed during one of the merger events retain a systematically prograde rotation at the final time. This demonstrates that the orbit of the ancient merger that helped to form the bulge could still remain in the kinematics of bulge stars.

We study the orbital properties of stars in four (published) simulations of thick disks formed by: i) accretion from disrupted satellites, ii) heating of a pre-existing thin disk by a minor merger, iii) radial migration and iv) gas rich mergers. We find that the distribution of orbital eccentricities are predicted to be different for each model: a prominent peak at low eccentricity is expected for the heating, migration and gas-rich merging scenarios, while the eccentricity distribution is broader and shifted towards higher values for the accretion model. These differences can be traced back to whether the bulk of the stars in each case is formed 'in-situ' or is 'accreted', and are robust to the peculiarities of each model. A simple test based on the eccentricity distribution of nearby thick disk stars may thus help elucidate the dominant formation mechanism of the Galactic thick disk.

Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major merger and the present day properties of the final remnant are typical of early type spiral galaxies, with an i band B/D ~0.65, a disk scale length of 7.2 kpc, g-r = 0.5 mag, an HI line width (W_{20}/2) of 238 km/sec and total magnitude i = -22.4. The key ingredients for the formation of a dominant stellar disk component after a major merger are: i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, ii) supernova feedback that is able to partially suppress star formation during mergers and iii) relative fading of the spheroidal component. The gas fraction of the progenitors' disks does not exceed 25% at z<3, emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of disks and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disk galaxies by the high likelihood of interactions and mergers for galaxy sized halos at relatively low z.

We investigate the formation of the stellar halos of four simulated disk galaxies using high resolution, cosmological SPH + N-Body simulations. These simulations include a self-consistent treatment of all the major physical processes involved in galaxy formation. The simulated galaxies presented here each have a total mass of ~10^12 M_sun, but span a range of merger histories. These simulations allow us to study the competing importance of in-situ star formation (stars formed in the primary galaxy) and accretion of stars from subhalos in the building of stellar halos in a LambdaCDM universe. All four simulated galaxies are surrounded by a stellar halo, whose inner regions (r < 20 kpc) contain both accreted stars, and an in-situ stellar population. The outer regions of the galaxies' halos were assembled through pure accretion and disruption of satellites. Most of the in-situ halo stars formed at high redshift out of smoothly accreted cold gas in the inner 1 kpc of the galaxies' potential wells, possibly as part of their primordial disks. These stars were displaced from their central locations into the halos through a succession of major mergers. We find that the two galaxies with recently quiescent merger histories have a higher fraction of in-situ stars (~20-50%) in their inner halos than the two galaxies with many recent mergers (~5-10% in-situ fraction). Observational studies concentrating on stellar populations in the inner halo of the Milky Way will be the most affected by the presence of in-situ stars with halo kinematics, as we find that their existence in the inner few tens of kpc is a generic feature of galaxy formation.

We analyse a high-resolution, fully cosmological, hydrodynamical disc galaxy simulation, to study the source of the double-exponential light profiles seen in many stellar discs, and the effects of stellar radial migration upon the spatio-temporal evolution of both the disc age and metallicity distributions. We find a "break" in the pure exponential stellar surface brightness profile, and trace its origin to a sharp decrease in the star formation per unit surface area, itself produced by a decrease in the gas volume density due to a warping of the gas disc. Star formation in the disc continues well beyond the break. We find that the break is more pronounced in bluer wavebands. By contrast, we find little or no break in the mass density profile. This is, in part, due to the net radial migration of stars towards the external parts of the disc. Beyond the break radius, we find that ~60% of the resident stars migrated from the inner disc, while ~25% formed in situ. Our simulated galaxy also has a minimum in the age profile at the break radius but, in disagreement with some previous studies, migration is not the main mechanism producing this shape. In our simulation, the disc metallicity gradient flattens with time, consistent with an "inside-out" formation scenario. We do not find any difference in the intensity or the position of the break with inclination, suggesting that perhaps the differences found in empirical studies are driven by dust extinction.

We use high resolution cosmological hydrodynamical simulations to demonstrate that cold flow gas accretion, particularly along filaments, modifies the standard picture of gas accretion and cooling onto galaxy disks. In the standard picture, all gas is initially heated to the virial temperature of the galaxy as it enters the virial radius. Low mass galaxies are instead dominated by accretion of gas that stays well below the virial temperature, and even when a hot halo is able to develop in more massive galaxies there exist dense filaments that penetrate inside of the virial radius and deliver cold gas to the central galaxy. For galaxies up to ~L*, this cold accretion gas is responsible for the star formation in the disk at all times to the present. Even for galaxies at higher masses, cold flows dominate the growth of the disk at early times. Within this modified picture, galaxies are able to accrete a large mass of cold gas, with lower initial gas temperatures leading to shorter cooling times to reach the disk. Although star formation in the disk is mitigated by supernovae feedback, the short cooling times allow for the growth of stellar disks at higher redshifts than predicted by the standard model.

To date, fully cosmological hydrodynamic disk simulations to redshift zero have only been undertaken with particle-based codes, such as GADGET, Gasoline, or GCD+. In light of the (supposed) limitations of traditional implementations of smoothed particle hydrodynamics (SPH), or at the very least, their respective idiosyncrasies, it is important to explore complementary approaches to the SPH paradigm to galaxy formation. We present the first high-resolution cosmological disk simulations to redshift zero using an adaptive mesh refinement (AMR)-based hydrodynamical code, in this case, RAMSES. We analyse the temporal and spatial evolution of the simulated stellar disks' vertical heating, velocity ellipsoids, stellar populations, vertical and radial abundance gradients (gas and stars), assembly/infall histories, warps/lopsideness, disk edges/truncations (gas and stars), ISM physics implementations, and compare and contrast these properties with our sample of cosmological SPH disks, generated with GCD+. These preliminary results are the first in our long-term Galactic Archaeology Simulation program.

Polar Ring Galaxies, such as NGC4650A, are a class of galaxy which have two kinematically distinct components that are inclined by almost 90 degrees to each other. These striking galaxies challenge our understanding of how galaxies form; the origin of their distinct components has remained uncertain, and the subject of much debate. We use high-resolution cosmological simulations of galaxy formation to show that Polar Ring Galaxies are simply an extreme example of the angular moment misalignment that occurs during the hierarchical structure formation characteristic of Cold Dark Matter cosmology. In our model, Polar Ring Galaxies form through the continuous accretion of gas whose angular momentum is misaligned with the central galaxy.

We present a new set of multi-million particle SPH simulations of the formation of disk dominated galaxies in a cosmological context. Some of these galaxies are higher resolution versions of the models already described in Governato et al (2007). To correctly compare simulations with observations we create artificial images of our simulations and from them measure photometric Bulge to Disk (B/D) ratios and disk scale lengths. We show how feedback and high force and mass resolution are necessary ingredients to form galaxies that have flatter rotation curves, larger I band disk scale lengths and smaller B/D ratios. A new simulated disk galaxy has an I-band disk scale length of 9.2 kpc and a B/D flux ratio of 0.64 (face on, dust reddened).

The difference in density profiles of the contributions from different density peaks to dark matter halos results in certain expectations about the Milky Way's stellar halo. We cut our simulated halo stars into two populations: those forming before/during the last major merger, and those accreted after the last major merger. The former population are more centrally located at z=0, while stars forming in low mass late forming proto-galaxies are spread through the halo. A difference in observed binding energy distinguishes these two populations. We look at possible chemical abundance signatures of the two populations. We also show that galaxies forming in isolated low sigma peaks will form from primordial material. Thus, even though the oldest stars are centrally concentrated as they originated in the early collapsing, densest regions, primordial stars would be found distributed throughout the halo. Thus, the lack of observed metal free stars can be taken as directly constraining the Population III IMF, and the lowest metallicity observed stars can be interpreted as holding clues to the chemical yields of Pop III stars.

We use cosmological, chemo-dynamical, smoothed particle hydrodynamical simulations of Milky-Way-analogue galaxies to find the expected present-day distributions of both metal-free stars that formed from primordial gas and the oldest star populations. We find that metal-free stars continue to form until z~4 in halos that are chemically isolated and located far away from the biggest progenitor of the final system. As a result, if the Population III initial mass function allows stars with low enough mass to survive until z=0 (< 0.8 Msol), they would be distributed throughout the Galactic halo. On the other hand, the oldest stars form in halos that collapsed close to the highest density peak of the final system, and at z=0 they are located preferentially in the central region of the Galaxy, i.e., in the bulge. According to our models, these trends are not sensitive to the merger histories of the disk galaxies or the implementation of supernova feedback. Furthermore, these full hydrodynamics results are consistent with our N-body results in Paper I, and lend further weight to the conclusion that surveys of low-metallicity stars in the Galactic halo can be used to directly constrain the properties of primordial stars. In particular, they suggest that the current lack of detections of metal-free stars implies that their lifetimes were shorter than a Hubble time, placing constraints on the metal-free initial mass function.

We study the spatial distribution of Galactic metal-free stars by combining an extremely high-resolution (7.8 X 10^5 solar masses per particle) Cold Dark Matter N-body simulation of the Milky-Way with a semi-analytic model of metal enrichment. This approach allows us to resolve halos with virial temperatures down to the 10^4K atomic cooling limit, and it is sufficiently flexible to make a number of robust conclusions, despite the extremely uncertain properties of the first stars. Galactic metal-free stars are formed over a large redshift range, which peaks at z~10, but continues down to z~5, contributing stars at wide range of Galactocentric radii. Stars containing only metals generated by primordial stars are similarly widespread. Neither changing the efficiency of metal dispersal by two orders of magnitude, nor drastically changing the approximations in our semi-analytical model can affect these result. Thus, if they have sufficiently long lifetimes, a significant number of stars formed in initially primordial star clusters should be found in the nearby Galactic halo regardless of the specifics of metal-free star formation. Observations of metal abundances in Galactic halo stars should be taken as directly constraining the properties of primordial stars, and the lack of metal-free halo stars today should be taken as strongly suggesting a 0.8 solar mass lower limit on the primordial initial mass function.

We employ N-body, smoothed particle hydrodynamical simulations, including detailed treatment of chemical enrichment, to follow a gas-rich merger which results in a galaxy with disk morphology. We trace the kinematic, structural and chemical properties of stars formed before, during, and after the merger. We show that such a merger produces two exponential disk components, with the older, hotter component having a scale-length 20% larger than the later-forming, cold disk. Rapid star formation during the merger quickly enriches the protogalactic gas reservoir, resulting in high metallicities of the forming stars. These stars form from gas largely polluted by Type II supernovae, which form rapidly in the merger-induced starburst. After the merger, a thin disk forms from gas which has had time to be polluted by Type Ia supernovae. Abundance trends are plotted, and we examine the proposal that increased star formation during gas-rich mergers may explain the high alpha-to-iron abundance ratios which exist in the relatively high-metallicity thick disk component of the Milky Way.

Increasingly large populations of disk galaxies are now being observed at increasingly high redshifts, providing new constraints on our knowledge of how such galaxies evolve. Are these observations consistent with a cosmology in which structures form hierarchically? To probe this question, we employ SPH/N-body galaxy scale simulations of late-type galaxies. We examine the evolution of these simulated disk galaxies from redshift 1 to 0, looking at the mass-size and luminosity-size relations, and the thickness parameter, defined as the ratio of scale-height to scale-length. The structural parameters of our simulated disks settle down quickly, and after redshift z=1 the galaxies evolve to become only slightly flatter. Our present day simulated galaxies are larger, more massive, less bright, and redder than at z=1. The inside-out nature of the growth of our simulated galaxies reduces, and perhaps eliminates, expectations of evolution in the size-mass relation.

In order to explain disk galaxy formation within the hierarchical structure formation, it seems that gas rich mergers must play an important role. We review here our previous studies which have shown the importance of mergers at high redshift being gas rich, in the formation of both the stellar halo and thick disk components of disk galaxies. Regulation of star formation in the building blocks of our galaxy is required to form a low mass low metallicity stellar halo. This regulation results in high redshift, gas rich mergers during which the thick disk forms. In these proceedings, we categorise stars from our simulated disk galaxy into thin and thick disk components by using the Toomre diagram.

Observations suggest that the structural parameters of disk galaxies have not changed greatly since redshift 1. We examine whether these observations are consistent with a cosmology in which structures form hierarchically. We use SPH/N-body galaxy-scale simulations to simulate the formation and evolution of Milky-Way-like disk galaxies by fragmentation, followed by hierarchical merging. The simulated galaxies have a thick disk, that forms in a period of chaotic merging at high redshift, during which a large amount of alpha-elements are produced, and a thin disk, that forms later and has a higher metallicity. Our simulated disks settle down quickly and do not evolve much since redshift z~1, mostly because no major mergers take place between z=1 and z=0. During this period, the disk radius increases (inside-out growth) while its thickness remains constant. These results are consistent with observations of disk galaxies at low and high redshift.

The stellar halos of spiral galaxies bear important chemo-dynamical signatures of galaxy formation. We present here the analysis of 89 semi-cosmological spiral galaxy simulations, spanning ~ 4 magnitudes in total galactic luminosity. These simulations sample a wide variety of merging histories and show significant dispersion in halo metallicity at a given total luminosity - more than a factor of ten in metallicity. Our preliminary analysis suggests that galaxies with a more extended merging history possess halos which have younger and more metal rich stellar populations than the stellar halos associated with galaxies with a more abbreviated assembly. A correlation between halo metallicity and its surface brightness has also been found, reflecting the correlation between halo metallicity and its stellar mass. Our simulations are compared with recent Hubble Space Telescope observations of resolved stellar halos in nearby spirals.

Seven cosmological hydrodynamic simulations of disk galaxy formation are analyzed to determine the alignment of the disk within the dark matter halo and the internal structure of the halo. We find that the orientation of the outer halo, beyond ~0.1 r_vir, is unaffected by the presence of the disk. In contrast, the inner halo is aligned such that the halo minor axis aligns with the disk axis. The relative orientation of these two regions of the halo are uncorrelated. The alignment of the disk and inner halo appears to take place simultaneously through their joint evolution. The disconnect between these two regions of the halo should be taken into account when modelling tidal streams in the halos of disk galaxies and when calculating intrinsic alignments of disk galaxies based on the properties of dark matter halos.

Despite having entered the era of "precision cosmology," the formation of galaxies within the favored CDM cosmological paradigm remains problematic. By relating our N-body/smooth particle hydro-dynamical simulation to an extensive range of Galactic and extragalactic observations, we shed light on the formation and evolution of the Milky Way and other late type galaxies. In light of recent observations of the stellar populations of extragalactic thick disks, we examine the proposal that the thick disk formed in a high redshift period characterized by gas rich merging. We show such a scenario to be consistent with color observations. We then follow the evolution of structural parameters of the subsequently formed disk galaxy, from redshift ~1 to the present. Consistent with observation, little evolution in the ratio of scale-height to scale-length is found in our simulated galaxy, despite its somewhat chaotic origins. The simulated galaxy in this report forms part of a larger study of a suite of galaxies, with which these issues are to be studied in detail.

The recently emerging conviction that thick disks are prevalent in disk galaxies, and their seemingly ubiquitous old ages, means that the formation of the thick disk, perhaps more than any other component, holds the key to unravelling the evolution of the Milky Way, and indeed all disk galaxies. In Paper I, we proposed that the thick disk was formed in an epoch of gas rich mergers, at high redshift. This hypothesis was based on comparing N-body/SPH simulations to a variety of Galactic and extragalactic observations, including stellar kinematics, ages and chemical properties.Here examine our thick disk formation scenario in light of the most recent observations of extragalactic thick disks. In agreement, our simulted thick disks are old and relatively metal rich, with V-I colors that do not vary significantly with distance from the plane. Further, we show that our proposal results in an enhancement of alpha-elements in thick disk stars as compared with thin disk stars, consistent with observations of the relevant populations of the Milky Way. We also find that our scenario naturally leads to the formation of an old metal weak stellar halo population with high alpha-element abundances.

The disk galaxy simulated using our chemo-dynamical galaxy formation code, GCD+, is shown to have a thick disk component. This is evidenced by the velocity dispersion versus age relation for solar neighbourhood stars, which clearly shows an abrupt increase in velocity dispersion at lookback time of approximately 8 Gyrs, and is in excellent agreement with observation. These thick disk stars are formed from gas which is accreted to the galaxy during a chaotic period of hierarchical clustering at high redshift. This formation scenario is shown to be consistent with observations of both the Milky Way and extragalactic thick disks.

Using stellar population synthesis techniques, we explore the photometric signatures of white dwarf progenitor dominated galactic halos, in order to constrain the fraction of halo mass that may be locked-up in white dwarf stellar remnants. We first construct a 10^9 M_sun stellar halo using the canonical Salpeter initial stellar mass distribution, and then allow for an additional component of low- and intermediate-mass stars, which ultimately give rise to white dwarf remnants. Microlensing observations towards the Large Magellanic Cloud, coupled with several ground-based proper motion surveys, have led to claims that in excess of 20% of the dynamical mass of the halo (10^12 M_sun) might be found in white dwarfs. Our results indicate that (1) even if only 1% of the dynamical mass of the dark halo today could be attributed to white dwarfs, their main sequence progenitors at high redshift (z ~ 3) would have resulted in halos more than 100 times more luminous than those expected from conventional initial mass functions alone, and (2) any putative halo white dwarf progenitor dominated initial mass function component, regardless of its dynamical importance, would be virtually impossible to detect at the present-day, due to its extremely faint surface brightness.

How do late type spiral galaxies form within the context of a CDM cosmology? We contrast N-body, smoothed particle hydrodynamical simulations of galaxy formation which employ two different supernova feedback mechanisms. Observed mass and metallicity distributions of the stellar halos of the Milky Way and M31 provide constraints on these models. A strong feedback model, incorporating an adiabatic phase in star burst regions, better reproduces these observational properties than our comparative thermal feedback model. This demonstrates the importance of energy feedback in regulating star formation in small systems, which collapse at early epochs, in the evolution of late type disk galaxies.

Observational evidence has suggested the possibility of a Galactic halo which is dominated by white dwarfs (WDs). While debate continues concerning the interpretation of this evidence, it is clear that an initial mass function (IMF) biased heavily toward WD precursors (1 < m/Msol < 8), at least in the early Universe, would be necessary in generating such a halo. Within the framework of homogeneous, closed-box models of Galaxy formation, such biased IMFs lead to an unavoidable overproduction of carbon and nitrogen relative to oxygen (as measured against the abundance patterns in the oldest stars of the Milky Way). Using a three-dimensional Tree N-body smoothed particle hydrodynamics code, we study the dynamics and chemical evolution of a galaxy with different IMFs. Both invariant and metallicity-dependent IMFs are considered. Our variable IMF model invokes a WD-precursor-dominated IMF for metallicities less than 5% solar (primarily the Galactic halo), and the canonical Salpeter IMF otherwise (primarily the disk). Halo WD density distributions and C,N/O abundance patterns are presented. While Galactic haloes comprised of ~5% (by mass) of WDs are not supported by our simulations, mass fractions of ~1-2% cannot be ruled out. This conclusion is consistent with the present-day observational constraints.

The present day chemical and dynamical properties of the Milky Way bear the imprint of the Galaxy's formation and evolutionary history. One of the most enduring and critical debates surrounding Galactic evolution is that regarding the competition between ``satellite accretion'' and ``monolithic collapse''; the apparent strong correlation between orbital eccentricity and metallicity of halo stars was originally used as supporting evidence for the latter. While modern-day unbiased samples no longer support the claims for a significant correlation, recent evidence has been presented by Chiba & Beers (2000,AJ,119,2843) for the existence of a minor population of high-eccentricity metal-deficient halo stars. It has been suggested that these stars represent the signature of a rapid (if minor) collapse phase in the Galaxy's history. Employing velocity- and integrals of motion-phase space projections of these stars, coupled with a series of N-body/Smoothed Particle Hydrodynamic (SPH) chemodynamical simulations, we suggest an alternative mechanism for creating such stars may be the recent accretion of a polar orbit dwarf galaxy.

The present-day chemical and dynamical properties of the Milky Way are signatures of the Galaxy's formation and evolution. Using a self consistent chemodynamical evolution code we examine these properties within the currently favoured paradigm for galaxy formation - hierarchical clustering within a CDM cosmology. Our Tree N-body/Smoothed Particle Hydrodynamics code includes a self-consistent treatment of gravity, hydrodynamics, radiative cooling, star formation, supernova feedback and chemical enrichment. Two models are described which explore the role of small-scale density perturbations in driving the evolution of structure within the Milky Way. The relationship between metallicity and kinematics of halo stars are quantified and the implications for galaxy formation discussed. While high-eccentricity halo stars have previously been considered a signature of ``rapid collapse'', we suggest that many such stars may have come from recently accreted satellites.

Halo initial mass functions (IMFs) heavily biased toward white dwarf (WD) precursors ~1-8 Msun have been suggested as a suitable mechanism for explaining MACHO events. However, by simple chemical evolution argument, Gibson & Mould (1997; GM97) pointed out that such WD-heavy IMFs (wdIMF) cause the overproduction of carbon and nitrogen. We re-examine this problem using numerical simulations.

The methodology involved in deriving the Hubble Constant via the calibration of the corrected peak luminosities of Type Ia supernovae (SNe) is reviewed. We first present a re-analysis of the Calan-Tololo (C-T) and Center for Astrophysics (CfA) Type Ia SN surveys. Bivariate linear least squares and quadratic boot-strapped fits in peak apparent magnitude and light curve shape are employed to correct this heterogeneous sample of peak apparent magnitudes, resulting in an homogeneous (and excellent) secondary distance indicator: the so-called corrected peak luminosity. We next provide an empirical calibration for this corrected luminosity, using Cepheid-based distances for seven nearby spiral galaxies host to Type Ia SNe. Included in this sample is the spectroscopically peculiar SN 1991T (in NGC 4527), whose corrected peak luminosity is shown to be indistinguishable from that of so-called ``normal'' SNe. A robust value of the Hubble Constant is derived and shown to be H0=73+/-2(r)+/-7(s) km/s/Mpc.