Marie Martig

The Thick Disks of Spiral Galaxies as Relics from Gas-Rich, Turbulent, Clumpy Disks at High Redshift

The formation of thick stellar disks in spiral galaxies is studied. Simulations of gas-rich young galaxies show formation of internal clumps by gravitational instabilities, clump coalescence into a bulge, and disk thickening by strong stellar scattering. The bulge and thick disks of modern galaxies may form this way. Simulations of minor mergers make thick disks too, but there is an important difference. Thick disks made by internal processes have a constant scale height with galactocentric radius, but thick disks made by mergers flare. The difference arises because in the first case, perpendicular forcing and disk-gravity resistance are both proportional to the disk column density, so the resulting scale height is independent of this density. In the case of mergers, perpendicular forcing is independent of the column density and the low-density regions get thicker; the resulting flaring is inconsistent with observations. Late-stage gas accretion and thin-disk growth are shown to preserve the constant scale heights of thick disks formed by internal evolution. These results reinforce the idea that disk galaxies accrete most of their mass smoothly and acquire their structure by internal processes, in particular through turbulent and clumpy phases at high redshift.

Chemodynamical evolution of the Milky Way disk - I. The solar vicinity

[Abridged] In this first paper of this series, we present a new approach for studying the chemo-dynamical evolution in disk galaxies, which consists of fusing disk chemical evolution models with compatible numerical simulations of galactic disks. This method avoids known star formation and chemical enrichment problems encountered in simulations. Here we focus on the Milky Way, by using a detailed thin-disk chemical evolution model and a simulation in the cosmological context, with dynamical properties close to those of our Galaxy. We show that, due to radial migration from mergers at high redshift and the central bar at later times, a sizable fraction of old metal-poor high-[alpha/Fe] stars reaches the solar vicinity. This naturally accounts for a number of recent observations related to both the thin and thick disks, despite the fact that we use thin-disk chemistry only. Although significant radial mixing is present, the slope in the age-metallicity relation is only weakly affected, with a scatter compatible with recent observational work. While we find a smooth density distribution in the [O/Fe]-[Fe/H] plane, we can recover the observed discontinuity by selecting particles according to kinematic criteria used in high-resolution samples to define the thin and thick disks. We outline a new method for estimating the birth place of the Sun and predict that the most likely radius lies in the range 4.4 < r < 7.7 kpc (for a location at r = 8 kpc). A new, unifying model for the Milky Way thick disk is offered, where both mergers and radial migration play a role at different stages of the disk evolution. We show that in the absence of early-on massive mergers the vertical velocity dispersion of the oldest stars is underestimated by a factor of ~2 compared with observations. We can, therefore, argue that the Milky Way thick disk is unlikely to have been formed through a quiescent disk evolution.

On the frequency, intensity and duration of starburst episodes triggered by galaxy interactions and mergers

We investigate the intensity enhancement and the duration of starburst episodes triggered by major galaxy interactions and mergers. We analyze two large statistical datasets of numerical simulations. These have been obtained using two independent and different numerical techniques to model baryonic and dark matter evolution that are extensively compared for the first time. One is a Tree-SPH code, the other one is a grid-based N-body sticky-particles code. We show that, at low redshift, galaxy interactions and mergers in general trigger only moderate star formation enhancements. Strong starbursts where the star formation rate is increased by a factor greater than 5 are rare and found only in about 15% of major galaxy interactions and mergers. Merger-driven starbursts are also rather short-lived, with a typical duration of activity of a few 10 8 yr. These conclusions are found to be robust, independent of the numerical techniques and star formation models. At higher redshifts where galaxies contain more gas, gas inflow-induced starbursts are neither stronger nor longer than their local counterparts. In turn, the formation of massive gas clumps, results of local Jeans instability that can occur spontaneously in gas-rich disks or be indirectly favored by galaxy interactions, could play a more important role in determining the duration and intensity of star formation episodes.

Radial migration does little for Galactic disc thickening

Non-axisymmetric components, such as spirals and central bars, play a major role in shaping galactic discs. An important aspect of the disc secular evolution driven by these perturbers is the radial migration of stars. It has been suggested recently that migration can populate a thick-disc component from inner-disc stars with high vertical energies. Since this has never been demonstrated in simulations, we study in detail the effect of radial migration on the disc velocity dispersion and disc thickness, by separating simulated stars into migrators and non-migrators. We apply this method to three isolated barred Tree-SPH N-body galaxies with strong radial migration. Contrary to expectations, we find that as stellar samples migrate, on the average, their velocity dispersion change (by as much as 50%) in such a way as to approximately match the non-migrating population at the radius at which they arrive. We show that, in fact, migrators suppress heating in parts of the disc. To confirm the validity of our findings, we also apply our technique to three cosmological re-simulations, which use a completely different simulation scheme and, remarkably, find very similar results. We believe the inability of migration to thicken discs is a fundamental property of internal disc evolution, irrespective of the migration mechanism at work. We explain this with the approximate conservation of the (average) vertical and radial actions rather than the energy. This “action mixing” can be used to constrain the migration rate in the Milky Way: estimates of the average vertical action in observations for different populations of stars should reveal flattening with radius for older groups of stars.

Chemodynamical evolution of the Milky Way disk II: Variations with Galactic radius and height above the disk plane

the date of receipt and acceptance should be inserted later Abstract. In the first paper of this series (paper I) we presented a new approach for studying the chemodynamical evolution in disk galaxies, focusing on the Milky Way. While in paper I we studied extensively the Solarvicinity, here we extend these results to different distances from the Galactic center, looking for variations of observables that can be related to on-going and future spectroscopic surveys. By separating the effects of kinematic heating and radial migration, we show that migration is much more important, even for the oldest and hottest stellar population. The distributions of stellar birth guiding radii and final guiding radii (signifying contamination from migration and heating, respectively) widen with increasing distance from the Galactic center. As a result, the slope in the age-metallicity relation flattens significantly at Galactic radii larger than solar. We predict that the metallicity distributions of (unbiased) samples at different distances from the Galactic center peak at approximately the same value, (Fe/H) ≈− 0.15 dex, and have similar metal-poor tails extending to (Fe/H) ≈− 1.3 dex. In contrast, the metal-rich tail decreases with increasing radius, thus giving rise to the expected decline of mean metallicity with radius. Similarly, the (Mg/Fe) distribution always peaks at ≈ 0.15 dex, but its low-end tail is lost as radius increases, while the high-end diminishes at (Mg/Fe) ≈ 0.45 dex. The radial metallicity and (Mg/Fe) gradients in our model show significant variations with height above the plane because of changes in the mixture of stellar ages. An inversion in the radial metallicity gradient is found from negative to weakly positive (at r < 10 kpc), and from positive to negative for the (Mg/Fe) gradient, with increasing distance from the disk plane. We relate this to the combined effect of (i) the predominance of young stars close to the disk plane and old stars away from it, (ii) the more concentrated older stellar component, and (iii) the flaring of mono-age disk populations. We also investigate the effect of recycled gas flows on the mean (Fe/H) and find that in the region 4 < r < 12 kpc the introduced errors are less than 0.05-0.1 dex, related to the fact that inward and outward flows mostly cancel in that radial range. We show that radial migration cannot compete with the inside-out formation of the disk, exposed by the more centrally concentrated older disk populations, and consistent with recent observations.

A DIVERSITY OF PROGENITORS AND HISTORIES FOR ISOLATED SPIRAL GALAXIES

We analyze a suite of 33 cosmological simulations of the evolution of Milky-Way-mass galaxies in low-density environments. Our sample spans a broad range of Hubble types at z = 0, from nearly bulgeless disks to bulge-dominated galaxies. Despite the fact that a large fraction of the bulge is typically in place by z ~ 1, we find no significant correlation between the morphology at z = 1 and at z = 0. The z = 1 progenitors of disk galaxies span a range of morphologies, including smooth disks, unstable disks, interacting galaxies, and bulge-dominated systems. By z ~ 0.5, spiral arms and bars are largely in place and the progenitor morphology is correlated with the final morphology. We next focus on late-type galaxies with a bulge-to-total ratio (B/T) 1. We find that the galaxies with the lowest B/T tend to have a quiet baryon input history, with no major mergers at z < 2, and with a low and constant gas accretion rate that keeps a stable angular-momentum direction. More violent merger or gas accretion histories lead to galaxies with more prominent bulges. Most disk galaxies have a bulge Sersic index n ≤ 2. The galaxies with the highest bulge Sersic index tend to have histories of intense gas accretion and disk instability rather than active mergers.

ON THE FORMATION OF GALACTIC THICK DISKS

Recent spectroscopic observations in the Milky Way suggest that the chemically defined thick disk (stars with high [alpha/Fe] ratios and thus old) has a significantly smaller scale-length than the thin disk. This is in apparent contradiction with observations of external edge-on galaxies, where the thin and thick components have comparable scale-lengths. Moreover, while observed disks do not flare (scale-height does not increase with radius), numerical simulations suggest that disk flaring is unavoidable, resulting from both environmental effects and secular evolution. Here we address these problems by studying two different suites of simulated galactic disks formed in the cosmological context. We show that the scale-heights of coeval populations always increase with radius. However, the total population can be decomposed morphologically into thin and thick disks, which do not flare. We relate this to the disk inside-out formation, where younger populations have increasingly larger scale-lengths and flare at progressively larger radii. In this new picture, thick disks are composed of the imbedded flares of mono-age stellar populations. Assuming that disks form inside out, we predict that morphologically defined thick disks must show a decrease in age (or [alpha/Fe] ratios) with radius and that coeval populations should always flare. This also explains the observed inversion in the metallicity and [alpha/Fe] gradients for stars away from the disk midplane in the Milky Way. The results of this work are directly linked to, and can be seen as evidence of, inside-out disk growth.

The ATLAS3D project – XXII. Low-efficiency star formation in early-type galaxies: hydrodynamic models and observations

We study the global efficiency of star formation in high-resolution hydrodynamical simulations of gas discs embedded in isolated early-type and spiral galaxies. Despite using a universal local law to form stars in the simulations, we find that the early-type galaxies are offset from the spirals on the large-scale Kennicutt relation, and form stars two to five times less efficiently. This offset is in agreement with previous results on morphological quenching: gas discs are more stable against star formation when embedded in early-type galaxies due to the lower disc self-gravity and increased shear. As a result, these gas discs do not fragment into dense clumps and do not reach as high densities as in the spiral galaxies. Even if some molecular gas is present, the fraction of very dense gas (typically above 10(4) cm(-3)) is significantly reduced, which explains the overall lower star formation efficiency. We also analyse a sample of local early-type and spiral galaxies, measuring their CO and H i surface densities and their star formation rates as determined by their non-stellar 8 mu m emission. As predicted by the simulations, we find that the early-type galaxies are offset from the Kennicutt relation compared to the spirals, with a twice lower efficiency. Finally, we validate our approach by performing a direct comparison between models and observations. We run a simulation designed to mimic the stellar and gaseous properties of NGC 524, a local lenticular galaxy, and find a gas disc structure and global star formation rate in good agreement with the observations. Morphological quenching thus seems to be a robust mechanism, and is also consistent with other observations of a reduced star formation efficiency in early-type galaxies in the COLD GASS survey. This lower efficiency of star formation is not enough to explain the formation of the whole red sequence, but can contribute to the reddening of some galaxies.

Red giant masses and ages derived from carbon and nitrogen abundances

We show that the masses of red giant stars can be well predicted from their photospheric carbon and nitrogen abundances, in conjunction with their spectroscopic stellar labels log g, Teff, and [Fe/H]. This is qualitatively expected from mass-dependent post-main-sequence evolution. We here establish an empirical relation between these quantities by drawing on 1475 red giants with asteroseismic mass estimates from Kepler that also have spectroscopic labels from Apache Point Observatory Galactic Evolution Experiment (APOGEE) DR12. We assess the accuracy of our model, and find that it predicts stellar masses with fractional rms errors of about 14 per cent (typically 0.2 M⊙). From these masses, we derive ages with rms errors of 40 per cent. This empirical model allows us for the first time to make age determinations (in the range 1–13 Gyr) for vast numbers of giant stars across the Galaxy. We apply our model to ∼52 000 stars in APOGEE DR12, for which no direct mass and age information was previously available. We find that these estimates highlight the vertical age structure of the Milky Way disc, and that the relation of age with [α/M] and metallicity is broadly consistent with established expectations based on detailed studies of the solar neighbourhood.

Young [α/Fe]-enhanced stars discovered by CoRoT and APOGEE: What is their origin?

We report the discovery of a group of apparently young CoRoT red-giant stars exhibiting enhanced [α/Fe] abundance ratios (as determined from APOGEE spectra) with respect to solar values. Their existence is not explained bystandard chemical evolution models of the Milky Way, and shows that the chemical-enrichment history of the Galactic disc is more complex. We find similar stars in previously published samples for which isochrone-ages could be reliably obtained, although in smaller relative numbers. This might explain why these stars have not previously received attention. The young [α/Fe]-rich stars are much more numerous in the CoRoT-APOGEE (CoRoGEE) inner-field sample than in any other high-resolution sample available at present because only CoRoGEE can explore the inner-disc regions and provide ages for its field stars. The kinematic properties of the young [α/Fe]-rich stars are not clearly thick-disc like, despite their rather large distances from the Galactic mid-plane. Our tentative interpretation of these and previous intriguing observations in the Milky Way is that these stars were formed close to the end of the Galactic bar, near corotation – a region where gas can be kept inert for longer times than in other regions that are more frequently shocked by the passage of spiral arms. Moreover, this is where the mass return from older inner-disc stellar generations is expected to be highest (according to an inside-out disc-formation scenario), which additionally dilutes the in-situ gas. Other possibilities to explain these observations (e.g., a recent gas-accretion event) are also discussed.