No Arabic abstract
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 Ways, 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.
We study the structure, age and metallicity gradients, and dynamical evolution using a cosmological zoom-in simulation of a Milky Way-mass galaxy from the Feedback in Realistic Environments project. In the simulation, stars older than 6 Gyr were formed in a chaotic, bursty mode and have the largest vertical scale heights (1.5-2.5 kpc) by z=0, while stars younger than 6 Gyr were formed in a relatively calm, stable disk. The vertical scale height increases with stellar age at all radii, because (1) stars that formed earlier were thicker at birth, and (2) stars were kinematically heated to an even thicker distribution after formation. Stars of the same age are thicker in the outer disk than in the inner disk (flaring). These lead to positive vertical age gradients and negative radial age gradients. The radial metallicity gradient is neg- ative at the mid-plane, flattens at larger disk height |Z|, and turns positive above |Z|~1.5kpc. The vertical metallicity gradient is negative at all radii, but is steeper at smaller radii. These trends broadly agree with observations in the Milky Way and can be naturally understood from the age gradients. The vertical stellar density profile can be well-described by two components, with scale heights 200-500 pc and 1-1.5 kpc, respectively. The thick component is a mix of stars older than 4 Gyr which formed through a combination of several mechanisms. Our results also demonstrate that it is possible to form a thin disk in cosmological simulations even with strong stellar feedback.
We analyze 494 main sequence turnoff and subgiant stars from the AMBRE:HARPS survey. These stars have accurate astrometric information from textit{Gaia}/DR1, providing reliable age estimates with relative uncertainties of $pm1-2$ Gyr and allowing precise orbital determinations. The sample is split based on chemistry into a low-[Mg/Fe] sequence, which are often identified as thin disk stellar populations, and a high-[Mg/Fe] sequence, which are often associated with the thick disk. We find that the high-[Mg/Fe] chemical sequence has extended star formation for several Gyr and is coeval with the oldest stars of the low-[Mg/Fe] chemical sequence: both the low- and high-[Mg/Fe] sequences were forming stars at the same time. The high-[Mg/Fe] stellar populations are only vertically extended for the oldest, most-metal poor and highest [Mg/Fe] stars. When comparing vertical velocity dispersion for both sequences, the high-[Mg/Fe] sequence has lower velocity dispersion than the low-[Mg/Fe] sequence for stars of similar age. Identifying either group as thin or thick disk based on chemistry is misleading. The stars belonging to the high-[Mg/Fe] sequence have perigalacticons that originate in the inner disk, while the perigalacticons of stars on the low-[Mg/Fe] sequence are generally around the solar neighborhood. From the orbital properties of the stars, the high-and low-[Mg/Fe] sequences are most likely a reflection of the chemical enrichment history of the inner and outer disk populations; radial mixing causes both populations to be observed in situ at the solar position. Based on these results, we emphasize that it is important to be clear in defining what populations are being referenced when using the terms thin and thick disk, and that ideally the term thick disk should be reserved for purely geometric definitions to avoid confusion and be consistent with definitions in external galaxies.
The Milky Way disk consists of two prominent components - a thick, alpha-rich, low-metallicity component and a thin, metal-rich, low-alpha component. External galaxies have been shown to contain thin and thick disk components, but whether distinct components in the [$alpha$/Fe]-[Z/H] plane exist in other Milky Way-like galaxies is not yet known. We present VLT-MUSE observations of UGC 10738, a nearby, edge-on Milky Way-like galaxy. We demonstrate through stellar population synthesis model fitting that UGC 10738 contains alpha-rich and alpha-poor stellar populations with similar spatial distributions to the same components in the Milky Way. We discuss how the finding that external galaxies also contain chemically distinct disk components may act as a significant constraint on the formation of the Milky Ways own thin and thick disk.
We map the stellar structure of the Galactic thick disk and halo by applying color-magnitude diagram (CMD) fitting to photometric data from the SEGUE survey, allowing, for the first time, a comprehensive analysis of their structure at both high and low latitudes using uniform SDSS photometry. Incorporating photometry of all relevant stars simultaneously, CMD fitting bypasses the need to choose single tracer populations. Using old stellar populations of differing metallicities as templates we obtain a sparse 3D map of the stellar mass distribution at |Z|>1 kpc. Fitting a smooth Milky Way model comprising exponential thin and thick disks and an axisymmetric power-law halo allows us to constrain the structural parameters of the thick disk and halo. The thick-disk scale height and length are well constrained at 0.75+-0.07 kpc and 4.1+-0.4 kpc, respectively. We find a stellar halo flattening within ~25 kpc of c/a=0.88+-0.03 and a power-law index of 2.75+-0.07 (for 7<R_{GC}<~30 kpc). The model fits yield thick-disk and stellar halo densities at the solar location of rho_{thick,sun}=10^{-2.3+-0.1} M_sun pc^{-3} and rho_{halo,sun}=10^{-4.20+-0.05} M_sun pc^{-3}, averaging over any substructures. Our analysis provides the first clear in situ evidence for a radial metallicity gradient in the Milky Ways stellar halo: within R<~15 kpc the stellar halo has a mean metallicity of [Fe/H]=-1.6, which shifts to [Fe/H]=-2.2 at larger radii. Subtraction of the best-fit smooth and symmetric model from the overall density maps reveals a wealth of substructures at all latitudes, some attributable to known streams and overdensities, and some new. A simple warp cannot account for the low latitude substructure, as overdensities occur simultaneously above and below the Galactic plane. (abridged)
We present MWFitting, a method to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) from TRILEGAL models to real data. We apply MWFitting to photometric data from the first three years of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud and the Sagittarius Stream, our main sample spans a total area of $sim$2,300 deg$^2$ distributed across the DES footprint. We further explore a smaller subset ($sim$ 1,300 deg$^2$) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3%. Based on the best-fit models, we create simulated stellar catalogues covering the whole DES footprint down to $g = 24$ magnitude. Comparisons of data and simulations provide evidence for a break in the power law index describing the stellar density of the Milky Way (MW) halo. Several previously discovered stellar over-densities are recovered in the residual stellar density map, showing the reliability of MWFitting in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict MW star counts for surveys such as LSST and Euclid.