No Arabic abstract
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.
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 examine the distribution of young stars associated with the spiral arms of a simulated L* cosmological disk galaxy. We find age patterns orthogonal to the arms which are not inconsistent with the predictions of classical density wave theory, a view further supported by recent observations of face-on Grand Design spirals such as M51. The distribution of metals within a simulated ~0.1L* disk is presented, reinforcing the link between star formation, the age-metallicity relation, and the metallicity distribution function.
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 use cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies to study the relative importance of the main stellar components, i.e., discs, bulges, and bars, at redshift zero. The main aim of this work is to understand if estimates of the structural parameters of these components determined from kinematics (as is usually done in simulations) agree well with those obtained using a photometric bulge/disc/bar decomposition (as done in observations). To perform such a comparison, we have produced synthetic observations of the simulation outputs with the Monte-Carlo radiative transfer code SUNRISE and used the BUDDA code to make 2D photometric decompositions of the resulting images (in the i and g bands). We find that the kinematic disc-to-total ratio (D/T) estimates are systematically and significantly lower than the photometric ones. While the maximum D/T ratios obtained with the former method are of the order of 0.2, they are typically >0.4, and can be as high as 0.7, according to the latter. The photometric decomposition shows that many of the simulated galaxies have bars, with Bar/T ratios in the range 0.2-0.4, and that bulges have in all cases low Sersic indices, resembling observed pseudo-bulges instead of classical ones. Simulated discs, bulges and bars generally have similar (g-i) colours, which are in the blue tail of the distribution of observed colours. This is not due to the presence of young stars, but rather to low metallicities and poor gas content in the simulated galaxies, which makes dust extinction low. Photometric decompositions thus match the component ratios usually quoted for spiral galaxies better than kinematic decompositions, but the shift is insufficient to make the simulations consistent with observed late-type systems.
We have examined the resolved stellar populations at large galactocentric distances along the minor axis (from 10 kpc up to between 40 and 75 kpc), with limited major axis coverage, of six nearby highly-inclined Milky Way-mass disc galaxies using HST data from the GHOSTS survey. We select red giant branch stars to derive stellar halo density profiles. The projected minor axis density profiles can be approximated by power laws with projected slopes of between $-2$ and $-3.7$ and a diversity of stellar halo masses of $1-6times 10^{9}M_{odot}$, or $2-14%$ of the total galaxy stellar masses. The typical intrinsic scatter around a smooth power law fit is $0.05-0.1$ dex owing to substructure. By comparing the minor and major axis profiles, we infer projected axis ratios $c/a$ at $sim 25$ kpc between $0.4-0.75$. The GHOSTS stellar haloes are diverse, lying between the extremes charted out by the (rather atypical) haloes of the Milky Way and M31. We find a strong correlation between the stellar halo metallicities and the stellar halo masses. We compare our results with cosmological models, finding good agreement between our observations and accretion-only models where the stellar haloes are formed by the disruption of dwarf satellites. In particular, the strong observed correlation between stellar halo metallicity and mass is naturally reproduced. Low-resolution hydrodynamical models have unrealistically high stellar halo masses. Current high-resolution hydrodynamical models appear to predict stellar halo masses somewhat higher than observed but with reasonable metallicities, metallicity gradients and density profiles.