No Arabic abstract
We present results from textsc{GigaEris}, a cosmological, $N$-body hydrodynamical ``zoom-in simulation of the formation of a Milky Way-sized galaxy with unprecedented resolution, encompassing of order a billion particles within the refined region. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion. We focus on the early assembly of the galaxy, down to redshift $z=4.4$. The simulated galaxy has properties consistent with extrapolations of the main sequence of star-forming galaxies to higher redshifts and levels off to a star formation rate of $sim$60$, M_{odot}$yr$^{-1}$ at $z=4.4$. A compact, thin rotating stellar disk with properties analogous to those of low-redshift systems arises already at $z sim 8$-9. The galaxy rapidly develops a multi-component structure, and the disk, at least at these early stages, does not grow upside-down as often reported in the literature. Rather, at any given time, newly born stars contribute to sustain a thin disk, while the thick disk grows from stars that are primarily added through accretion and mergers. The kinematics reflect the early, ubiquitous presence of a thin disk, as a stellar disk component with $v_phi/sigma_R$ larger than unity is already present at $z sim 9$-10. Our results suggest that high-resolution spectro-photometric observations of very high-redshift galaxies should find thin rotating disks, consistent with the recent discovery of cold rotating gas disks by ALMA. Finally, we present synthetic images for the JWST NIRCam camera, showing how the early disk would be easily detectable already at $z sim 7$.
We compare the star-formation history and dynamics of the Milky Way (MW) with the properties of distant disk galaxies. During the first ~4 Gyr of its evolution, the MW formed stars with a high star-formation intensity (SFI), Sigma_SFR~0.6 Msun/yr/kpc2 and as a result, generated outflows and high turbulence in its interstellar medium. This intense phase of star formation corresponds to the formation of the thick disk. The formation of the thick disk is a crucial phase which enables the MW to have formed approximately half of its total stellar mass by z~1 which is similar to MW progenitor galaxies selected by abundance matching. This agreement suggests that the formation of the thick disk may be a generic evolutionary phase in disk galaxies. Using a simple energy injection-kinetic energy relationship between the 1-D velocity dispersion and SFI, we can reproduce the average perpendicular dispersion in stellar velocities of the MW with age. This relationship, its inferred evolution, and required efficiency are consistent with observations of galaxies from z~0-3. The high turbulence generated by intense star formation naturally resulted in a thick disk, a chemically well-mixed ISM, and is the mechanism that links the evolution of MW to the observed characteristics of distant disk galaxies.
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.
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.
By means of idealized, dissipationless N-body simulations which follow the formation and subsequent buckling of a stellar bar, we study the characteristics of boxy/peanut-shaped bulges and compare them with the properties of the stellar populations in the Milky Way bulge. The main results of our modeling, valid for the general family of boxy/peanut shaped bulges, are the following: (i) because of the spatial redistribution in the disk initiated at the epoch of bar formation, stars from the innermost regions to the outer Lindblad resonance of the stellar bar are mapped into a boxy bulge; (ii) the contribution of stars to the local bulge density depends on their birth radius: stars born in the innermost disk tend to dominate the innermost regions of the boxy bulge, while stars originating closer to the OLR are preferentially found in the outer regions of the boxy/peanut structure; (iii) stellar birth radii are imprinted in the bulge kinematics, the larger the birth radii of stars ending up in the bulge, the greater their rotational support and the higher their line-of- sight velocity dispersions (but note that this last trend depends on the bar viewing angle); (iv) the higher the classical bulge-over-disk ratio, the larger its fractional contribution of stars at large vertical distance from the galaxy mid-plane. (ABRIDGED) On the basis of their chemical and kinematic characteristics, the results of our modeling suggests that the populations A, B and C, as defined by the ARGOS survey, can be associated, respectively, with the inner thin disk, to the young thick and to the old thick disk, following the nomenclature recently suggested for stars in the solar neighborhood by Haywood et al. (2013).
Bulges are commonly believed to form in the dynamical violence of galaxy collisions and mergers. Here we model the stellar kinematics of the Bulge Radial Velocity Assay (BRAVA), and find no sign that the Milky Way contains a classical bulge formed by scrambling pre-existing disks of stars in major mergers. Rather, the bulge appears to be a bar, seen somewhat end-on, as hinted from its asymmetric boxy shape. We construct a simple but realistic N-body model of the Galaxy that self-consistently develops a bar. The bar immediately buckles and thickens in the vertical direction. As seen from the Sun, the result resembles the boxy bulge of our Galaxy. The model fits the BRAVA stellar kinematic data covering the whole bulge strikingly well with no need for a merger-made classical bulge. The bar in our best fit model has a half-length of ~ 4kpc and extends 20 degrees from the Sun-Galactic Center line. We use the new kinematic constraints to show that any classical bulge contribution cannot be larger than ~ 8% of the disk mass. Thus the Galactic bulge is a part of the disk and not a separate component made in a prior merger. Giant, pure-disk galaxies like our own present a major challenge to the standard picture in which galaxy formation is dominated by hierarchical clustering and galaxy mergers.