No Arabic abstract
A short overview is presented of several topics concerning the evolution of the Milky Way (MW) in a cosmological context. In particular, the metallicity distribution of the MW halo is derived analytically and the halo metallicity and abundance patterns are compared to those of Local Group galaxies. The inside-out formation of the MW disk is supported by the observed evolution of the abundance gradients, while their magnitude suggests that the role of the Galactic bar has been negligible. Finally, the empirical foundations (age-metallicity relation and metallicity distribution) of the evolution of the solar neighborhood, which is the best studied galactic sub-system, have been seriously questioned recently.
We combine the Santa-Cruz Semi-Analytic Model (SAM) for galaxy formation and evolution with the circumgalactic medium (CGM) model presented in Faerman et al. (2020) to explore the CGM properties of $L^{*}$ galaxies. We use the SAM to generate a sample of galaxies with halo masses similar to the Milky Way (MW) halo, $M_{rm vir} approx 10^{12}~{rm M_{sun}}$, and find that the CGM mass and mean metallicity in the sample are correlated. We use the CGM masses and metallicities of the SAM galaxies as inputs for the FSM20 model, and vary the amount of non-thermal support. The density profiles in our models can be approximated by power-law functions with slopes in the range of $0.75 < a_n < 1.25$, with higher non-thermal pressure resulting in flatter distributions. We explore how the gas pressure, dispersion measure, OVI-OVIII column densities, and cooling rates behave with the gas distribution and total mass. We show that for CGM masses below $sim 3 times 10^{10}~{rm M_{sun}}$, photoionization has a significant effect on the column densities of OVI and OVIII. The combination of different MW CGM observations favors models with similar fractions in thermal pressure, magnetic fields/cosmic rays, and turbulent support, and with $M_{rm gas} sim 3-10 times 10^{10}~{rm M_{sun}}$. The MW OVI column requires $t_{rm cool}/t_{rm dyn} sim 4$, independent of the gas distribution. The AGN jet-driven heating rates in the SAM are enough to offset the CGM cooling, although exact balance is not required in star-forming galaxies. We provide predictions for the columns densities of additional metal ions - NV, NeVIII, and MgX.
We investigate the impact of cosmic rays (CR) and different modes of CR transport on the properties of Milky Way-mass galaxies in cosmological magneto-hydrodynamical simulations in the context of the AURIGA project. We systematically study how advection, anisotropic diffusion and additional Alfven-wave cooling affect the galactic disc and the circum-galactic medium (CGM). Global properties such as stellar mass and star formation rate vary little between simulations with and without various CR transport physics, whereas structural properties such as disc sizes, CGM densities or temperatures can be strongly affected. In our simulations, CRs affect the accretion of gas onto galaxies by modifying the CGM flow structure. This alters the angular momentum distribution which manifests itself as a difference in stellar and gaseous disc size. The strength of this effect depends on the CR transport model: CR advection results in the most compact discs while the Alfven-wave model resembles more the AURIGA model. The advection and diffusion models exhibit large ($rsim50$ kpc) CR pressure-dominated gas haloes causing a smoother and partly cooler CGM. The additional CR pressure smoothes small-scale density peaks and compensates for the missing thermal pressure support at lower CGM temperatures. In contrast, the Alfven-wave model is only CR pressure dominated at the disc-halo interface and only in this model the gamma-ray emission from hadronic interactions agrees with observations. In contrast to previous findings, we conclude that details of CR transport are critical for accurately predicting the impact of CR feedback on galaxy formation.
We study the effect of different Type Ia SN nucleosynthesis prescriptions on the Milky Way chemical evolution. To this aim, we run detailed one-infall and two-infall chemical evolution models, adopting a large compilation of yield sets corresponding to different white dwarf progenitors (near-Chandrasekar and sub-Chandrasekar) taken from the literature. We adopt a fixed delay time distribution function for Type Ia SNe , in order to avoid degeneracies in the analysis of the different nucleosynthesis channels. We also combine yields for different Type Ia SN progenitors in order to test the contribution to chemical evolution of different Type Ia SN channels. The results of the models are compared with recent LTE and NLTE observational data. We find that classical W7 and WDD2 models produce Fe masses and [$alpha$/Fe] abundance patterns similar to more recent and physical near-Chandrasekar and sub- Chandrasekar models. For Fe-peak elements, we find that the results strongly depend either on the white dwarf explosion mechanism (deflagration-to-detonation, pure deflagration, double detonation) or on the initial white dwarf conditions (central density, explosion pattern). The comparison of chemical evolution model results with observations suggests that a combination of near-Chandrasekar and sub-Chandrasekar yields is necessary to reproduce the data of V, Cr, Mn and Ni, with different fractions depending on the adopted massive stars stellar yields. This comparison also suggests that NLTE and singly ionised abundances should be definitely preferred when dealing with most of Fe-peak elements at low metallicity.
We apply a semi-analytic galaxy formation model to two high resolution cosmological N-body simulations to investigate analogues of the Milky Way system. We select these according to observed properties of the Milky Way rather than by halo mass as in most previous work. For disk-dominated central galaxies with stellar mass (5--7) x 10d10Msun, the median host halo mass is 1.4 x 10d12Msun, with 1 sigma dispersion in the range [0.86, 3.1] x 10d12Msun, consistent with dynamical measurements of the Milky Way halo mass. For any given halo mass, the probability of hosting a Milky Way system is low, with a maximum of ~20% in haloes of mass ~10d12Msun. The model reproduces the V-band luminosity function and radial profile of the bright (MV < -9) Milky Way satellites. Galaxy formation in low mass haloes is found to be highly stochastic, resulting in an extremely large scatter in the relation between MV (or stellar mass) for satellites and the depth of the subhalo potential well in which they live, as measured by the maximum of the rotation curve, Vmax. We conclude that the too big to fail problem is an artifact of selecting satellites in N-body simulations according to subhalo properties: in 10% of cases we find that three or fewer of the brightest (or most massive) satellites have Vmax > 30 km/s. Our model predicts that around half of the dark matter subhaloes with Vmax > 20 km/s host satellites fainter than MV = -9 and so may be missing from existing surveys.
We study the evolution of Milky Way thick and thin discs in the light of the most recent observational data. In particular, we analyze abundance gradients of O, N, Fe and Mg along the thin disc as well as the [Mg/Fe] vs. [Fe/H] relations and the metallicity distribution functions at different Galactocentric distances. We run several models starting from the two-infall paradigm, assuming that the thick and thin discs formed by means of two different infall episodes, and we explore several physical parameters, such as radial gas flows, variable efficiency of star formation, different times for the maximum infall onto the disc, different distributions of the total surface mass density of the thick disc and enriched gas infall. Our best model suggests that radial gas flows and variable efficiency of star formation should be acting together with the inside-out mechanism for the thin disc formation. The timescale for maximum infall onto the thin disc, which determines the gap between the formation of the two discs, should be $t_{max}simeq 3.25$ Gyr. The thick disc should have an exponential, small scale length density profile and gas infall on the inner thin disc should be enriched. We compute also the evolution of Gaia-Enceladus system and study the effects of possible interactions with the thick and thin discs. We conclude that the gas lost by Enceladus or even part of it could have been responsible for the formation of the thick disc but not the thin disc.