The Milky Way, like other disc galaxies, underwent violent mergers and accretion of smaller satellite galaxies in its early history. Thanks to $\mathit Gaia$-DR2 and spectroscopic surveys, the stellar remnants of such mergers have been identified. The chronological dating of such events is crucial to uncover the formation and evolution of the Galaxy at high redshift, but it has so far been challenging owing to difficulties in obtaining precise ages for these oldest stars. Here we combine asteroseismology -- the study of stellar oscillations -- with kinematics and chemical abundances, to estimate precise stellar ages ($\sim$ 11%) for a sample of stars observed by the $\mathit{Kepler}$ space mission. Crucially, this sample includes not only some of the oldest stars that were formed inside the Galaxy, but also stars formed externally and subsequently accreted onto the Milky Way. Leveraging this unprecedented temporal resolution in age, we provide overwhelming evidence in favour of models in which the Galaxy had already formed a significant population of its stars (which now reside mainly in its thick disc) before the in-fall of the satellite galaxy $\mathit{Gaia}$-$\mathit{Enceladus}$ 10 billions years ago.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We make use of APOGEE and $Gaia$ data to identify stars that are consistent with being born in the same association or star cluster as the Sun. We limit our analysis to stars that match solar abundances within their uncertainties, as they could have formed from the same Giant Molecular Cloud (GMC) as the Sun. We constrain the range of orbital actions that solar siblings can have with a suite of simulations of solar birth clusters evolved in static and time-dependent tidal fields. The static components of each galaxy model are the bulge, disk, and halo, while the various time-dependent components include a bar, spiral arms, and GMCs. In galaxy models without GMCs, simulated solar siblings all have $J_R < 122$ km $\rm s^{-1}$ kpc, $990 < L_z < 1986$ km $\rm s^{-1}$ kpc, and $0.15 < J_z < 0.58$ km $\rm s^{-1}$ kpc. Given the actions of stars in APOGEE and $Gaia$, we find 104 stars that fall within this range. One candidate in particular, Solar Sibling 1, has both chemistry and actions similar enough to the solar values that strong interactions with the bar or spiral arms are not required for it to be dynamically associated with the Sun. Adding GMCs to the potential can eject solar siblings out of the plane of the disk and increase their $J_z$, resulting in a final candidate list of 296 stars. The entire suite of simulations indicate that solar siblings should have $J_R < 122$ km $\rm s^{-1}$ kpc, $353 < L_z < 2110$ km $\rm s^{-1}$ kpc, and $J_z < 0.8$ km $\rm s^{-1}$ kpc. Given these criteria, it is most likely that the association or cluster that the Sun was born in has reached dissolution and is not the commonly cited open cluster M67.

We summarise recent results from analysis of APOGEE/Gaia data for stellar populations in the Galactic halo, disk, and bulge, leading to constraints on the contribution of dwarf galaxies and globular clusters to the stellar content of the Milky Way halo. Interpretation of the extant data in light of cosmological numerical simulations suggests that the Milky Way has been subject to an unusually intense accretion history at z >~ 1.5.

Studies of the kinematics and chemical compositions of Galactic globular clusters (GCs) enable the reconstruction of the history of star formation, chemical evolution, and mass assembly of the Galaxy. Using the latest data release (DR16) of the SDSS/APOGEE survey, we identify 3,090 stars associated with 46 GCs. Using a previously defined kinematic association, we break the sample down into eight separate groups and examine how the kinematics-based classification maps into chemical composition space, considering only $\alpha$ (mostly Si and Mg) elements and Fe. Our results show that: (i) The loci of both in situ and accreted subgroups in chemical space match those of their field counterparts; (ii) GCs from different individual accreted subgroups occupy the same locus in chemical space. This could either mean that they share a similar origin or that they are associated with distinct satellites which underwent similar chemical enrichment histories; (iii) The chemical compositions of the GCs associated with the low orbital energy subgroup defined by Massari and collaborators is broadly consistent with an in situ origin. However, at the low metallicity end, the distinction between accreted and in situ populations is blurred; (iv) Regarding the status of GCs whose origin is ambiguous, we conclude the following: the position in Si-Fe plane suggests an in situ}origin for Liller 1 and a likely accreted origin for NGC 5904 and NGC 6388. The case of NGC 288 is unclear, as its orbital properties suggest an accretion origin, its chemical composition suggests it may have formed in situ.

Over the course of its history, the Milky Way has ingested multiple smaller satellite galaxies. While these accreted stellar populations can be forensically identified as kinematically distinct structures within the Galaxy, it is difficult in general to precisely date the age at which any one merger occurred. Recent results have revealed a population of stars that were accreted via the collision of a dwarf galaxy, called \textit{Gaia}-Enceladus, leading to a substantial pollution of the chemical and dynamical properties of the Milky Way. Here, we identify the very bright, naked-eye star $\nu$\,Indi as a probe of the age of the early in situ population of the Galaxy. We combine asteroseismic, spectroscopic, astrometric, and kinematic observations to show that this metal-poor, alpha-element-rich star was an indigenous member of the halo, and we measure its age to be $11.0 \pm 0.7$ (stat) $\pm 0.8$ (sys)$\,\rm Gyr$. The star bears hallmarks consistent with it having been kinematically heated by the \textit{Gaia}-Enceladus collision. Its age implies that the earliest the merger could have begun was 11.6 and 13.2 Gyr ago at 68 and 95% confidence, respectively. Input from computations based on hierarchical cosmological models tightens (i.e. reduces) slightly the above limits.

The stellar mass in the halo of the Milky Way is notoriously difficult to determine, owing to the paucity of its stars in the solar neighbourhood. With tentative evidence from $\it{Gaia}$ that the nearby stellar halo is dominated by a massive accretion event -- referred to as $\it{Gaia-Enceladus}$ or Sausage -- these constraints are now increasingly urgent. We measure the mass in kinematically selected mono-abundance populations (MAPs) of the stellar halo between $-3 <$ [Fe/H] $< -1$ and $0.0 <$ [Mg/Fe] $< 0.4$ using red giant star counts from APOGEE DR14. We find that MAPs are well fit by single power laws on triaxial ellipsoidal surfaces, and we show that that the power law slope $\alpha$ changes such that high $\mathrm{[Mg/Fe]}$ populations have $\alpha \sim 4$, whereas low [Mg/Fe] MAPs are more extended with shallow slopes, $\alpha \sim 2$. We estimate the total stellar mass to be $M_{*,\mathrm{tot}} = 1.3^{+0.3}_{-0.2}\times10^{9}\ \mathrm{M_{\odot}}$, of which we estimate $\sim 0.9^{+0.2}_{-0.1} \times 10^{9} \mathrm{M_{\odot}}$ to be accreted. We estimate that the mass of accreted stars with $e > 0.7$ is $M_{*,\mathrm{accreted},e>0.7} = 3 \pm 1\ \mathrm{(stat.)}\pm 1\ \mathrm{(syst.)}\times10^{8}\ \mathrm{M_\odot}$, or $\sim 30-50\%$ of the accreted halo mass. If the majority of these stars $\it{are}$ the progeny of a massive accreted dwarf, this places an upper limit on its stellar mass, and implies a halo mass for the progenitor of $\sim 10^{10.2\pm0.2}\ \mathrm{M_\odot}$. This constraint not only shows that the $\it{Gaia-Enceladus}$/Sausage progenitor may not be as massive as originally suggested, but that the majority of the Milky Way stellar halo was accreted. These measurements are an important step towards fully reconstructing the assembly history of the Milky Way.

Studies of the ages, abundances, and motions of individual stars in the Milky Way provide one of the best ways to study the evolution of disk galaxies over cosmic time. The formation of the Milky Way's barred inner region in particular is a crucial piece of the puzzle of disk galaxy evolution. Using data from APOGEE and Gaia, we present maps of the kinematics, elemental abundances, and age of the Milky Way bulge and disk that show the barred structure of the inner Milky Way in unprecedented detail. The kinematic maps allow a direct, purely kinematic determination of the bar's pattern speed of 41+/-3 km/s/kpc and of its shape and radial profile. We find the bar's age, metallicity, and abundance ratios to be the same as those of the oldest stars in the disk that are formed in its turbulent beginnings, while stars in the bulge outside of the bar are younger and more metal-rich. This implies that the bar likely formed ~8 Gyr ago, when the decrease in turbulence in the gas disk allowed a thin disk to form that quickly became bar-unstable. The bar's formation therefore stands as a crucial epoch in the evolution of the Milky Way, a picture that is in line with the evolutionary path that emerges from observations of the gas kinematics in external disk galaxies over the last ~10 Gyr.

We investigate how the stellar and gas-phase He abundances evolve as functions of time within simulated star-forming disc galaxies with different star formation histories. We make use of a cosmological chemodynamical simulation for galaxy formation and evolution, which includes star formation, as well as energy and chemical enrichment feedback from asymptotic giant branch stars, core-collapse supernovae, and Type Ia supernovae. The predicted relations between the He mass fraction, $Y$, and the metallicity, $Z$, in the interstellar medium of our simulated disc galaxies depend on the past galaxy star formation history. In particular, $dY/dZ$ is not constant and evolves as a function of time, depending on the specific chemical element that we choose to trace $Z$; in particular, $dY/dX_{\text{O}}$ and $dY/dX_{\text{C}}$ increase as functions of time, whereas $dY/dX_{\text{N}}$ decreases. In the gas-phase, we find negative radial gradients of $Y$, due to the inside-out growth of our simulated galaxy discs as a function of time; this gives rise to longer chemical enrichment time scales in the outer galaxy regions, where we find lower average values for $Y$ and $Z$. Finally, by means of chemical evolution models, in the galactic bulge and inner disc, we predict steeper $Y$ versus age relations at high $Z$ than in the outer galaxy regions. We conclude that, for calibrating the assumed $Y$-$Z$ relation in stellar models, C, N, and C+N are better proxies for the metallicity than O, because they show steeper and less scattered relations.

Action space synthesizes the orbital information of stars and is well-suited to analyse the rich amount of kinematic disk substructure in the Gaia DR2 radial velocity sample (RVS). In this work, we revisit one of the strongest perturbers in the Milky Way (MW) disk: the m=2 bar. We investigate how its resonances affect the actions of individual test particle stars, i.e., (JR,Lz,Jz) estimated in an axisymmetric MW potential. We confirm that the stars' behaviour is well approximated by scattering and oscillation along a slope Delta JR / Delta Lz = l/m centered on the l:m resonance lines. The Outer Lindblad Resonance (OLR, l=+1,m=2) creates signatures in the stellar action space that can be used to identify the Galactic bar's OLR in the Gaia DR2 RVS data: (a) The JR dependence of the oscillation causes an overdensity ridge (underdensity region) at Lz larger (smaller) than the resonance line in the (Lz,JR) plane. (b) For the first time, we demonstrate that the OLR is expected to cause a gradient in average Jz with Lz across the resonance. (c) We show that the change of predominantly outward to inward motions at the OLR occurs along the resonance line in action space. The latter signature allows us to identify three candidates for the bar's OLR - and therefore its pattern speed Omega_bar - in the Gaia data within 3 kpc from the Sun: 1.85 Omega0, 1.2 Omega0, and 1.6 Omega0 (with ~0.1 Omega0 uncertainty). This demonstrates that (i) the local Gaia action data is consistent with both the short-fast and long-slow bar models in the literature, and that (ii) axisymmetrically estimated actions are a powerful diagnostic even in non-axisymmetric systems.

The kinematics of the Milky Way disc as a function of age are well measured at the solar radius, but have not been studied over a wider range of Galactocentric radii. Here, we measure the kinematics of mono-age, mono-$\mathrm{[Fe/H]}$ populations in the low and high $\mathrm{[\alpha/Fe]}$ discs between $4 \lesssim R \lesssim 13$ kpc and $|z| \lesssim 2$ kpc using 65,719 stars in common between APOGEE DR14 and $\it{Gaia}$ DR2 for which we estimate ages using a Bayesian neural network model trained on asteroseismic ages. We determine the vertical and radial velocity dispersions, finding that the low and high $\mathrm{[\alpha/Fe]}$ discs display markedly different age--velocity-dispersion relations (AVRs) and shapes $\sigma_z/\sigma_R$. The high $\mathrm{[\alpha/Fe]}$ disc has roughly flat AVRs and constant $\sigma_z/\sigma_R = 0.64\pm 0.04$, whereas the low $\mathrm{[\alpha/Fe]}$ disc has large variations in this ratio which positively correlate with the mean orbital radius of the population at fixed age. The high $\mathrm{[\alpha/Fe]}$ disc component's flat AVRs and constant $\sigma_z/\sigma_R$ clearly indicates an entirely different heating history. Outer disc populations also have flatter radial AVRs than those in the inner disc, likely due to the waning effect of spiral arms. Our detailed measurements of AVRs and $\sigma_z/\sigma_R$ across the disc indicate that low $\mathrm{[\alpha/Fe]}$, inner disc ($R \lesssim 10\,\mathrm{kpc}$) stellar populations are likely dynamically heated by both giant molecular clouds and spiral arms, while the observed trends for outer disc populations require a significant contribution from another heating mechanism such as satellite perturbations. We also find that outer disc populations have slightly positive mean vertical and radial velocities, likely because they are part of the warped disc.

Gaia DR2 has provided an unprecedented wealth of information about the kinematics of stars in the Solar neighbourhood, and has highlighted the degree of features in the Galactic disc. We confront the data with a range of bar and spiral models in both action-angle space, and the $R_{\mathrm{G}}-v_{\phi}$ plane. We find that the phase mixing induced by transient spiral structure creates ridges and arches in the local kinematics which are consistent with the Gaia data. We are able to produce a qualitatively good match to the data when combined with a bar with a variety of pattern speeds, and show that it is non trivial to decouple the effects of the bar and the spiral structure.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

Recent work indicates that the nearby Galactic halo is dominated by the debris from a major accretion event. We confirm that result from an analysis of APOGEE-DR14 element abundances and $\textit{Gaia}$-DR2 kinematics of halo stars. We show that $\sim$2/3 of nearby halo stars have high orbital eccentricities ($e \gtrsim 0.8$), and abundance patterns typical of massive Milky Way dwarf galaxy satellites today, characterised by relatively low [Fe/H], [Mg/Fe], [Al/Fe], and [Ni/Fe]. The trend followed by high $e$ stars in the [Mg/Fe]-[Fe/H] plane shows a change of slope at [Fe/H]$\sim-1.3$, which is also typical of stellar populations from relatively massive dwarf galaxies. Low $e$ stars exhibit no such change of slope within the observed [Fe/H] range and show slightly higher abundances of Mg, Al and Ni. Unlike their low $e$ counterparts, high $e$ stars show slightly retrograde motion, make higher vertical excursions and reach larger apocentre radii. By comparing the position in [Mg/Fe]-[Fe/H] space of high $e$ stars with those of accreted galaxies from the EAGLE suite of cosmological simulations we constrain the mass of the accreted satellite to be in the range $10^{8.5}\lesssim M_*\lesssim 10^{9}\mathrm{M_\odot}$. We show that the median orbital eccentricities of debris are largely unchanged since merger time, implying that this accretion event likely happened at $z\lesssim1.5$. The exact nature of the low $e$ population is unclear, but we hypothesise that it is a combination of $\textit{in situ}$ star formation, high $|z|$ disc stars, lower mass accretion events, and contamination by the low $e$ tail of the high $e$ population. Finally, our results imply that the accretion history of the Milky Way was quite unusual.

Orbital parameters, such as eccentricity and maximum vertical excursion, of stars in the Milky Way are an important tool for understanding its dynamics and evolution, but calculation of such parameters usually relies on computationally-expensive numerical orbit integration. We present and test a fast method for estimating these parameters using an application of the St\"ackel fudge, used previously for the estimation of action-angle variables. We show that the method is highly accurate, to a level of $<1\%$ in eccentricity, over a large range of relevant orbits and in different Milky Way-like potentials, and demonstrate its validity by estimating the eccentricity distribution of the RAVE-TGAS data set and comparing it to that from orbit integration. Using the method, the orbital characteristics of the $\sim 7$ million $\textit{Gaia}$ DR2 stars with radial velocity measurements are computed with Monte Carlo sampled errors in $\sim 116$ hours of parallelised cpu time, at a speed that we estimate to be $\sim 3$ to $4$ orders of magnitude faster than using numerical orbit integration. We demonstrate using this catalogue that $\textit{Gaia}$ DR2 samples a large range of orbits in the solar vicinity, down to those with $r_\mathrm{peri} \lesssim 2.5$ kpc, and out to $r_\mathrm{ap} \gtrsim 13$ kpc. We also show that many of the features present in orbital parameter space have a low mean $z_\mathrm{max}$, suggesting that they likely result from disk dynamical effects.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

Spectroscopic surveys of the Galaxy reveal that its disc stars exhibit a spread in $\mathrm{[\alpha/Fe]}$ at fixed $\mathrm{[Fe/H]}$, manifest at some locations as a bimodality. The origin of these diverse, and possibly distinct, stellar populations in the Galactic disc is not well understood. We examine the Fe and $\alpha$-element evolution of 133 Milky Way-like galaxies from the EAGLE simulation, to investigate the origin and diversity of their $\mathrm{[\alpha/Fe]}$-$\mathrm{[Fe/H]}$ distributions. We find that bimodal $\mathrm{[\alpha/Fe]}$ distributions arise in galaxies whose gas accretion histories exhibit episodes of significant infall at both early and late times, with the former fostering more intense star formation than the latter. The shorter characteristic consumption timescale of gas accreted in the earlier episode suppresses its enrichment with iron synthesised by Type Ia SNe, resulting in the formation of a high-$\mathrm{[\alpha/Fe]}$ sequence. We find that bimodality in $\mathrm{[\alpha/Fe]}$ similar to that seen in the Galaxy is rare, appearing in approximately 5 percent of galaxies in our sample. We posit that this is a consequence of an early gas accretion episode requiring the mass accretion history of a galaxy's dark matter halo to exhibit a phase of atypically-rapid growth at early epochs. The scarcity of EAGLE galaxies exhibiting distinct sequences in the $\mathrm{[\alpha/Fe]}$-$\mathrm{[Fe/H]}$ plane may therefore indicate that the Milky Way's elemental abundance patterns, and its accretion history, are not representative of the broader population of $\sim L^\star$ disc galaxies.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

The best way to trace back the history of star formation and mass assembly of the Milky Way disc is by combining chemical compositions, ages and phase-space information for a large number of disc stars. With the advent of large surveys of the stellar populations of the Galaxy, such data have become available and can be used to pose constraints on sophisticated models of galaxy formation. We use SDSS-III/APOGEE data to derive the first detailed 3D map of stellar density in the Galactic disc as a function of age, [Fe/H] and $\mathrm{[\alpha/Fe]}$. We discuss the implications of our results for the formation and evolution of the disc, presenting new constraints on the disc structural parameters, stellar radial migration and its connection with disc flaring. We also discuss how our results constrain the inside out formation of the disc, and determine the surface-mass density contributions at the solar radius for mono-age, mono-[Fe/H] populations.

The measurement of the structure of stellar populations in the Milky Way disk places fundamental constraints on models of galaxy formation and evolution. Previously, the disk's structure has been studied in terms of populations defined geometrically and/or chemically, but a decomposition based on stellar ages provides a more direct connection to the history of the disk, and stronger constraint on theory. Here, we use positions, abundances and ages for 31,244 red giant branch stars from the SDSS-APOGEE survey, spanning $3 < R_{\mathrm{gc}} < 15$ kpc, to dissect the disk into mono-age and mono-[Fe/H] populations at low and high [$\alpha$/Fe]. For each population, with $\Delta \mathrm{age} < 2$ Gyr and $\Delta \mathrm{[Fe/H]} < 0.1$ dex, we measure the structure and surface-mass density contribution. We find that low [$\alpha$/Fe] mono-age populations are fit well by a broken exponential, which increases to a peak radius and decreases thereafter. We show that this profile becomes broader with age, interpreted here as a new signal of disk heating and radial migration. High [$\alpha$/Fe] populations are well fit as single exponentials within the radial range considered, with an average scale length of $1.9\pm 0.1$ kpc. We find that the relative contribution of high to low [$\alpha$/Fe] populations at $R_0$ is $f_\Sigma = 18\% \pm 5\%$; high [$\alpha$/Fe] contributes most of the mass at old ages, and low [$\alpha$/Fe] at young ages. The low and high [$\alpha$/Fe] populations overlap in age at intermediate [Fe/H], although both contribute mass at $R_{0}$ across the full range of [Fe/H]. The mass weighted scale height $h_Z$ distribution is a smoothly declining exponential function. High [$\alpha$/Fe] populations are thicker than low [$\alpha$/Fe], and the average $h_Z$ increases steadily with age, between 200 and 600 pc.

Formation of globular clusters (GCs), the Galactic bulge, or galaxy bulges in general, are important unsolved problems in Galactic astronomy. Homogeneous infrared observations of large samples of stars belonging to GCs and the Galactic bulge field are one of the best ways to study these problems. We report the discovery by APOGEE of a population of field stars in the inner Galaxy with abundances of N, C, and Al that are typically found in GC stars. The newly discovered stars have high [N/Fe], which is correlated with [Al/Fe] and anti-correlated with [C/Fe]. They are homogeneously distributed across, and kinematically indistinguishable from, other field stars in the same volume. Their metallicity distribution is seemingly unimodal, peaking at [Fe/H]~-1, thus being in disagreement with that of the Galactic GC system. Our results can be understood in terms of different scenarios. N-rich stars could be former members of dissolved GCs, in which case the mass in destroyed GCs exceeds that of the surviving GC system by a factor of ~8. In that scenario, the total mass contained in so-called "first-generation" stars cannot be larger than that in "second-generation" stars by more than a factor of ~9 and was certainly smaller. Conversely, our results may imply the absence of a mandatory genetic link between "second generation" stars and GCs. Last, but not least, N-rich stars could be the oldest stars in the Galaxy, the by-products of chemical enrichment by the first stellar generations formed in the heart of the Galaxy.