We trace the formation and advection of several elements within a cosmological adaptive mesh refinement simulation of an L* galaxy. We use nine realisations of the same initial conditions with different stellar Initial Mass Functions (IMFs), mass limits for type-II and type-Ia supernovae (SNII, SNIa) and stellar lifetimes to constrain these sub-grid phenomena. Our code includes self-gravity, hydrodynamics, star formation, radiative cooling and feedback from multiple sources within a cosmological framework. Under our assumptions of nucleosynthesis we find that SNII with progenitor masses of up to 100 Msun are required to match low metallicity gas oxygen abundances. Tardy SNIa are necessary to reproduce the classical chemical evolution knee in [O/Fe]-[Fe/H]: more prompt SNIa delayed time distributions do not reproduce this feature. Within our framework of hydrodynamical mixing of metals and galaxy mergers we find that chemical evolution is sensitive to the shape of the IMF and that there exists a degeneracy with the mass range of SNII. We look at the abundance plane and present the properties of different regions of the plot, noting the distinct chemical properties of satellites and a series of nested discs that have greater velocity dispersions, are more alpha-rich and metal poor with age.

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 present a new chemodynamical code - Ramses-CH - for use in simulating the self-consistent evolution of chemical and hydrodynamical properties of galaxies within a fully cosmological framework. We build upon the adaptive mesh refinement code Ramses, which includes a treatment of self-gravity, hydrodynamics, star formation, radiative cooling, and supernovae feedback, to trace the dominant isotopes of C, N, O, Ne, Mg, Si, and Fe. We include the contribution of Type Ia and II supernovae, in addition to low- and intermediate-mass asymptotic giant branch stars, relaxing the instantaneous recycling approximation. The new chemical evolution modules are highly flexible and portable, lending themselves to ready exploration of variations in the underpining stellar and nuclear physics. We apply Ramses-CH to the cosmological simulation of a typical L* galaxy, demonstrating the successful recovery of the basic empirical constraints regarding, [{\alpha}/Fe]-[Fe/H] and Type Ia/II supernovae rates.

We highlight two research strands related to our ongoing chemodynamical Galactic Archaeology efforts: (i) the spatio-temporal infall rate of gas onto the disk, drawing analogies with the infall behaviour imposed by classical galactic chemical evolution models of inside-out disk growth; (ii) the radial age gradient predicted by spectrophometric models of disk galaxies. In relation to (i), at low-redshift, we find that half of the infall onto the disk is gas associated with the corona, while half can be associated with cooler gas streams; we also find that gas enters the disk preferentially orthogonal to the system, rather than in-plane. In relation to (ii), we recover age gradient troughs/inflections consistent with those observed in nature, without recourse to radial migrations.

We present a new chemodynamical code based on the adaptive mesh refinement code RAMSES. The new code uses Eulerian hydrodynamics and N-body dynamics in a cosmological framework to trace the production and advection of several chemical species. It is the first such code to follow the self-consistent evolution of chemical elements in cosmological volumes while maintaining sub-kiloparsec resolution. The code will be used to simulate disk galaxies and explore the influence of chemical evolution models and star formation on galactic abundance ratios.

We present for the first time predictions for UV line emission of intergalactic and circumgalactic gas from Adaptive Mesh Resolution (AMR) Large Scale Structure (LSS) simulations at redshifts 0.3<z<1.2, with specific emphasis on observability with current and near-future UV instrumentation. In three transitions of interest (Lya, OVI and CIV) there is a clear bimodality in the type of objects : the overwhelming majority of flux stems from discrete, compact sources, while a much larger volume fraction is filled by more tenuous gas. We characterise both object types with regard to number density, physical size and shape, brightness, luminosity, velocity structure, mass, temperature, ionisation state, and metal content. Degrading AMR grids to characteristic resolutions of available (such as FIREBall) or foreseeable instrumentation, allows to assess which inferences can be drawn from currently possible observations, and set foundations to prepare observing strategies for future missions. In general, the faint emission of the IGM and filamentary structure remains beyond capabilities of instruments with only short duration exposure potential (stratospheric balloons), even for optimistic assumptions for Lya, while the yet fainter metal line transitions for these structures will remain challenging for long duration exposures (space-based telescopes), mostly due to low metallicity pushing them more than three orders of magnitudes in brightness below Lya radiation. For the circum-galactic medium (CGM) the situation is more promising, and it is foreseeable that in the near future we will not only just dectect such sources, but the combination of all three lines in addition to velocity information will yield valuable insight into the physical processes at hand, illuminating important mechanisms during the formation of galaxies and their backreaction onto the IGM from whence they formed. (abrigded)

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 analyze the physical properties and infall rates of the circum-galactic gas around disks obtained in multi-resolved, cosmological, AMR simulations. At intermediate and low redshifts, disks are embedded into an extended, hot, tenuous corona that contributes largely in fueling the disk with non-enriched gas whereas the accretion of enriched gas from tidal streams occurs throughout episodic events. We derive an infall rate close to the disk of the same value as the one of the star formation rate in the disk and its temporal evolution as a function of galacto-centric radius nicely shows that the growth of galactic disks proceeds according to an inside-out formation scenario.

The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMA

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.

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.

We explore galaxy properties in general and properties of host galaxies of gamma-ray bursts (GRBs) in particular, using N-body/Eulerian hydrodynamic simulations and the stellar population synthesis model, Starburst99, to infer observable properties. We identify simulated galaxies that have optical star formation rate (SFR) and SFR-to-luminosity ratio similar to those observed in a well-defined sample of ten host galaxies. Each of the numerical counterparts are found in catalogs at the same redshifts as the observed hosts. The counterparts are found to be low-mass galaxies, with low mass-to-light ratio, recent epoch of formation, and high ratio between the SFR and the average of the SFR. When compared to the overall galaxy population, they have colors much bluer than the high-mass star-forming galaxy population. Although their SFRs span a range of values, the specific rates of the numerical counterparts are equal to or higher than the median values estimated at the different redshifts. We also emphasize the strong relationships between the specific star formation rate (SFR) and quantities known to reflect the star formation history of galaxies, i.e. color and mass-to-light ratio: At intermediate redshift, the faintest and bluest galaxies are also the objects with the highest specific rates. These results suggest that GRB host galaxies are likely to be drawn from the high specific SFR sub-population of galaxies, rather than the high SFR galaxy population. Finally, as indicated by our catalogs, in an extended sample, the majority of GRB host galaxies is expected to have specific SFRs higher than found in the magnitude-limited sample studied here.

We review recent developments in the field of chemodynamical simulations of elliptical galaxies, highlighting (in an admittedly biased fashion) the work conducted with our cosmological N-body/SPH code GCD+. We have demonstrated previously the recovery of several primary integrated early-type system scaling relations (e.g. colour-magnitude relation, L_X-T_X-[Fe/H]_X) when employing a phenomenological AGN heating scheme in conjunction with a self-consistent treatment of star formation, supernovae feedback, radiative cooling, chemical enrichment, and stellar/X-ray population synthesis. Here we emphasise characteristics derived from the full spatial information contained within the simulated dataset, including stellar and coronal morphologies, metallicity distribution functions, and abundance gradients.

The problem of galaxy formation and its dependence on thermodynamic properties is addressed by using Eulerian hydrodynamic numerical simulations of large scale structure formation. Global galaxy properties are explored in simulations including gravitation, shock heating and cooling processes, and following self-consistently the chemical evolution of a primordial composition hydrogen-helium plasma without assuming collisional ionization equilibrium. The galaxy formation model is mainly based on the identification of converging dense cold gas regions. We show that the evolution at low redshift of the observed cosmic star formation rate density is reproduced, and that the galaxy-like object mass function is dominated by low-mass objects. The galaxy mass functions are well described by a two power-law Schechter function whose parameters are in good agreement with observational fits of the galaxy luminosity function. The high-mass end of the galaxy mass function includes objects formed at early epochs and residing in high-mass dark matter halos whereas the low-mass end includes galaxies formed at later epochs and active in their ``stellar'' mass formation. Finally, the influence of two other physical processes, photoionization and non-equipartition processes between electrons, ions and neutrals of the cosmological plasma is discussed and the modifications on galaxy formation are examined.

We use numerical simulations of large scale structure formation to explore the cosmological properties of Gamma-Ray Burst (GRB) host galaxies. Among the different sub-populations found in the simulations, we identify the host galaxies as the most efficient star-forming objects, i.e. galaxies with high specific star formation rates. We find that the host candidates are low-mass, young galaxies with low to moderate star formation rate. These properties are consistent with those observed in GRB hosts, most of which are sub-luminous, blue galaxies. Assuming that host candidates are galaxies with high star formation rates would have given conclusions inconsistent with the observations. The specific star formation rate, given a galaxy mass, is shown to increase as the redshift increases. The low mass of the putative hosts makes them difficult to detect with present day telescopes and the probability density function of the specific star formation rate is predicted to change depending on whether or not these galaxies are observed.

Using N-body/hydrodynamic simulations, the influence of non-equipartition processes on the thermal and dynamical properties of cosmological baryonic gas is investigated. We focus on a possible departure from equilibrium between electrons, ions and neutral atoms in low temperature (10^4-10^6 K) and weakly ionized regions of the intergalactic medium. The simulations compute the energy exchanges between ions, neutrals and electrons, without assuming thermal equilibrium. They include gravitation, shock heating and cooling processes, and follow self-consistently the chemical evolution of a primordial composition hydrogen-helium plasma without assuming collisional ionization equilibrium. At high redshift, a significant fraction of the intergalactic medium is found to be warmer and weakly ionized in simulations with non-equipartition processes than in simulations in which the cosmological plasma is considered to be in thermodynamic equilibrium. With a semi-analytical study of the out of equilibrium regions we show that, during the formation of cosmic structures, departure from equilibrium in accreted plasma results from the competition between the atomic cooling processes and the elastic processes between heavy particles and electrons. Our numerical results are in agreement with this semi-analytical model. Therefore, since baryonic matter with temperatures around 10^4 K is a reservoir for galaxy formation, non-equipartition processes are expected to modify the properties of the objects formed.