Large stellar surveys are revealing the chemodynamical structure of the Galaxy across a vast spatial extent. However, the many millions of low-resolution spectra observed to date are yet to be fully exploited. We employ The Cannon, a data-driven approach to estimating abundances, to obtain detailed abundances from low-resolution (R = 1800) LAMOST spectra, using the GALAH survey as our reference. We deliver five (for dwarfs) or six (for giants) estimated abundances representing five different nucleosynthetic channels, for 3.9 million stars, to a precision of 0.05 - 0.23 dex. Using wide binary pairs, we demonstrate that our abundance estimates provide chemical discriminating power beyond metallicity alone. We show the coverage of our catalogue with radial, azimuthal and dynamical abundance maps, and examine the neutron capture abundances across the disk and halo, which indicate different origins for the in-situ and accreted halo populations. LAMOST has near-complete Gaia coverage and provides an unprecedented perspective on chemistry across the Milky Way.

The age and chemical characteristics of the Galactic bulge link to the formation and evolutionary history of the Galaxy. Data-driven methods and large surveys enable stellar ages and precision chemical abundances to be determined for vast regions of the Milky Way, including the bulge. Here, we use the data-driven approach of The Cannon, to infer the ages and abundances for 125,367 stars in the Milky Way, using spectra from APOGEE DR14. We examine the ages and metallicities of 1654 bulge stars within $R_{\text{GAL}}<3.5$ kpc. We focus on fields with $b<12^\circ$, and out to longitudes of $l<15^\circ$. We see that stars in the bulge are about twice as old ($\tau=8$ Gyrs), on average, compared to those in the solar neighborhood ($\tau=4$ Gyrs), with a larger dispersion in [Fe/H] ($\approx0.38$ compared to 0.23 dex). This age gradient comes primarily from the low-$\alpha$ stars. Looking along the Galactic plane, the very central field in the bulge shows by far the largest dispersion in [Fe/H] ($\sigma_{[Fe/H]}\approx0.4$ dex) and line of sight velocity ($\sigma_{vr}\approx90$ km/s), and simultaneously the smallest dispersion in age. Moving out in longitude, the stars become kinematically colder and less dispersed in [Fe/H], but show a much broader range of ages. We see a signature of the X-shape within the bulge at a latitude of $b=8^\circ$, but not at $b=12^\circ$. Future APOGEE and other survey data, with larger sampling, affords the opportunity to extend our approach and study in more detail, to place stronger constraints on models of the Milky Way.

Stellar variability is driven by a multitude of internal physical processes that depend on fundamental stellar properties. These properties are our bridge to reconciling stellar observations with stellar physics, and for understanding the distribution of stellar populations within the context of galaxy formation. Numerous ongoing and upcoming missions are charting brightness fluctuations of stars over time, which encode information about physical processes such as rotation period, evolutionary state (such as effective temperature and surface gravity), and mass (via asteroseismic parameters). Here, we explore how well we can predict these stellar properties, across different evolutionary states, using only photometric time series data. To do this, we implement a convolutional neural network, and with data-driven modeling we predict stellar properties from light curves of various baselines and cadences. Based on a single quarter of \textit{Kepler} data, we recover stellar properties, including surface gravity for red giant stars (with an uncertainty of $\lesssim$ 0.06 dex), and rotation period for main sequence stars (with an uncertainty of $\lesssim$ 5.2 days, and unbiased from $\approx$5 to 40 days). Shortening the \textit{Kepler} data to a 27-day TESS-like baseline, we recover stellar properties with a small decrease in precision, $\sim$0.07 dex for log $g$ and $\sim$5.5 days for $P_{\rm rot}$, unbiased from $\approx$5 to 35 days. Our flexible data-driven approach leverages the full information content of the data, requires minimal feature engineering, and can be generalized to other surveys and datasets. This has the potential to provide stellar property estimates for many millions of stars in current and future surveys.

We explore the chemodynamical properties of a sample of barred galaxies in the Auriga magneto-hydrodynamical cosmological zoom-in simulations, which form boxy/peanut (b/p) bulges, and compare these to the Milky Way (MW). We show that the Auriga galaxies which best reproduce the chemodynamical properties of stellar populations in the MW bulge have quiescent merger histories since redshift $z\sim3.5$: their last major merger occurs at $t_{\rm lookback}>12\,\rm Gyrs$, while subsequent mergers have a stellar mass ratio of $\leq$1:20, suggesting an upper limit of a few percent for the mass ratio of the recently proposed Gaia Sausage/Enceladus merger. These Auriga MW-analogues have a negligible fraction of ex-situ stars in the b/p region ($<1\%$), with flattened, thick disc-like metal-poor stellar populations. The average fraction of ex-situ stars in the central regions of all Auriga galaxies with b/p's is 3% -- significantly lower than in those which do not host a b/p or a bar. While the central regions of these barred galaxies contain the oldest populations, they also have stars younger than 5Gyrs (>30%) and exhibit X-shaped age and abundance distributions. Examining the discs in our sample, we find that in some cases a star-forming ring forms around the bar, which alters the metallicity of the inner regions of the galaxy. Further out in the disc, bar-induced resonances lead to metal-rich ridges in the $V_{\phi}-r$ plane -- the longest of which is due to the Outer Lindblad Resonance. Our results suggest the Milky Way has an uncommonly quiet merger history, which leads to an essentially in-situ bulge, and highlight the significant effects the bar can have on the surrounding disc.

A major goal in the field of galaxy formation is to understand the formation of the Milky Way's disk. The first step toward doing this is to empirically describe its present state. We use the new high-dimensional dataset of 19 abundances from 27,135 red clump APOGEE stars to examine the distribution of clusters defined using abundances. We explore different dimensionality reduction techniques and implement a non-parametric agglomerate hierarchical clustering method. We see that groups defined using abundances are spatially separated, as a function of age. Furthermore, the abundance groups represent different distributions in the [Fe/H]-age plane. Ordering our clusters by age reveals patterns suggestive of the sequence of chemical enrichment in the disk over time. Our results indicate that a promising avenue to trace the details of the disk's assembly is via a full interpretation of the empirical connections we report.

Planetary engulfment events have long been proposed as a lithium (Li) enrichment mechanism contributing to the population of Li-rich giants $(A(\mathrm{Li}) \geq 1.5$ dex). Using $GALAH$ survey data and MESA stellar models, we calculate the strength and duration of the Li enrichment signature produced in the convective envelope of a host star that has engulfed a hot Jupiter (HJ) companion. We consider solar-metallicity stars in the mass range of $1.0-2.0~\mathrm{M_{\odot}}$ and the Li supplied by a HJ of $1.0~\mathrm{M_{J}}$. We explore engulfment events that occur near the main sequence turn-off (MSTO) and out to orbital separations of $R_{\star}{\sim}~0.1~\mathrm{AU}= 22~\mathrm{R_{\odot}}$. We map our results onto the Hertzsprung-Russell (H-R) Diagram, revealing the parameter space where planetary engulfment events produce significant Li enrichment signatures. We also map the associated survival times of these signatures, which range across 9 orders of magnitude. Our calculations indicate that if the HJ engulfment event occurs near the MSTO or on the subgiant branch, Li enrichment can be measured at a $5\sigma$ confidence level and with strengths that exceed meteoritic abundance measurements. Moreover, for stars of $1.4~\mathrm{M_{\odot}}$, these signatures are predicted to survive for up to 1 Gyr. We determine that Li enrichment beyond the subgiant branch must be produced by other mechanisms, such as the Cameron-Fowler process or accretion of material from an AGB companion.

Upcoming million-star spectroscopic surveys have the potential to revolutionize our view of the formation and chemical evolution of the Milky Way. Realizing this potential requires automated approaches to optimize estimates of stellar properties, such as chemical element abundances, from the spectra. The volume and quality of the observations strongly motivate that these approaches should be data-driven. With this in mind, we introduce SSSpaNG: a data-driven Gaussian Process model of stellar spectra. We demonstrate the capabilities of SSSpaNG using a sample of APOGEE red clump stars, whose model parameters we infer via Gibbs sampling. Pooling information between stars to infer their covariance, we permit clear identification of the correlations between spectral pixels. Harnessing these correlations, we infer the true spectrum of each star, inpainting missing regions and denoising by a factor of at least 2-3 for low-signal-to-noise stars. As we marginalize over the covariance matrix of the spectra, the effective prior on these true spectra is non-Gaussian, non-convex and sparsifying, promoting typically small but occasionally large excursions from the mean. The high-fidelity true spectra produced will enable improved elemental abundance measurements for individual stars. Our model also allows us to quantify the information gained by observing portions of a star's spectrum, and thereby define the most mutually informative spectral regions. Using 25 windows centred on elemental absorption lines, we demonstrate that the iron-peak and alpha-process elements are particularly mutually informative for these spectra, and that the majority of information about a target window is contained in the 10-or-so most informative windows. Our information-gain metric has the potential to inform models of nucleosynthetic yields and optimize the design of future observations.

Here we present the discovery of 895 s-process-rich candidates from 454,180 giant stars observed by the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) using a data-driven approach. This sample constitutes the largest number of s-process enhanced stars ever discovered. Our sample includes 187 s-process-rich candidates that are enhanced in both barium and strontium, 49 stars with significant barium enhancement only and 659 stars that show only a strontium enhancement. Most of the stars in our sample are in the range of effective temperature and log g typical of red giant branch (RGB) populations, which is consistent with our observational selection bias towards finding RGB stars. We estimate that only a small fraction (0.5 per cent) of binary configurations are favourable for s-process enriched stars. The majority of our s-process-rich candidates (95 per cent) show strong carbon enhancements, whereas only five candidates (less than 3 percent) show evidence of sodium enhancement. Our kinematic analysis reveals that 97 percent of our sample are disc stars, with the other 3 percent showing velocities consistent with the Galactic halo. The scaleheight of the disc is estimated to be zh=0.634kpc, comparable with values in the literature. A comparison with the yields from asymptotic giant branch (AGB) models suggests that the main neutron source responsible for the Ba and Sr enhancements is the 13C(alpha,n)16O reaction. We conclude that the s-process-rich candidates may have received their overabundances via mass transfer from a previous ABG companion with an initial mass in the range 1-3Msun.

Theoretical models of stellar evolution predict that most of the lithium inside a star is destroyed as the star becomes a red giant. However, observations reveal that about 1% of red giants are peculiarly rich in lithium, often exceeding the amount in the interstellar medium or predicted from the Big Bang. With only about 150 lithium-rich giants discovered in the past four decades, and no distinguishing properties other than lithium enhancement, the origin of lithium-rich giant stars is one of the oldest problems in stellar astrophysics. Here we report the discovery of 2,330 low-mass (1 to 3$\,M_\odot$) lithium-rich giant stars, which we argue are consistent with internal lithium production that is driven by tidal spin-up by a binary companion. Our sample reveals that most lithium-rich giants have helium-burning cores ($80^{+7}_{-6}\%$), and that the frequency of lithium-rich giants rises with increasing stellar metallicity. We find that while planet accretion may explain some lithium-rich giants, it cannot account for the majority that have helium-burning cores. We rule out most other proposed explanations as the primary mechanism for lithium-rich giants, including all stages related to single star evolution. Our analysis shows that giants remain lithium-rich for only about two million years. A prediction from this lithium depletion timescale is that most lithium-rich giants with a helium-burning core have a binary companion.

The Orion complex is arguably the most studied star-forming region in the Galaxy. While stars are still being born in the Orion nebula, the oldest part was believed to be no more than 13 Myr old. In order to study the full hierarchy of star formation across the Orion complex, we perform a clustering analysis of the Ori OB1a region using new stellar surveys and derive robust ages for each identified stellar aggregate. We use Gaia DR2 parameters supplemented with radial velocities from the GALAH and APOGEE surveys to perform clustering of the Ori OB1a association. Five overdensities are resolved in a six-dimensional parameter space (positions, distance, proper motions, and radial velocity). Most correspond to previously known structures (ASCC 16, 25 Orionis, ASCC 20, ASCC 21). We use Gaia DR2, Pan-STARRS1 and 2MASS photometry to fit isochrones to the colour-magnitude diagrams of the identified clusters. The ages of the clusters can thus be measured with ~10% precision. While four of the clusters have ages between 11 and 13 Myr, the ASCC 20 cluster stands out at an age of 21 $\pm$ 3 Myr. This is significantly greater than the age of any previously known component of the Orion complex. To some degree, all clusters overlap in at least one of the six phase-space dimensions. We argue that the formation history of the Orion complex, and its relation to the Gould belt, must be reconsidered. A significant challenge in reconstructing the history of the Ori OB1a association is to understand the impact of the newly discovered 21 Myr old population on the younger parts of the complex, including their formation.

It is well established that the chemical structure of the Milky Way exhibits a bimodality with respect to the $\alpha$-enhancement of stars at a given [Fe/H]. This has been studied largely based on a bulk $\alpha$ abundance, computed as a summary of several individual $\alpha$-elements. Inspired by the expected subtle differences in their nucleosynthetic origins, here we probe the higher level of granularity encoded in the inter-family [Mg/Si] abundance ratio. Using a large sample of stars with APOGEE abundance measurements, we first demonstrate that there is additional information in this ratio beyond what is already apparent in [$\alpha$/Fe] and [Fe/H] alone. We then consider Gaia astrometry and stellar age estimates to empirically characterize the relationships between [Mg/Si] and various stellar properties. We find small but significant trends between this ratio and $\alpha$-enhancement, age, [Fe/H], location in the Galaxy, and orbital actions. To connect these observed [Mg/Si] variations to a physical origin, we attempt to predict the Mg and Si abundances of stars with the galactic chemical evolution model Chempy. We find that we are unable to reproduce abundances for the stars that we fit, which highlights tensions between the yield tables, the chemical evolution model, and the data. We conclude that a more data-driven approach to nucleosynthetic yield tables and chemical evolution modeling is necessary to maximize insights from large spectroscopic surveys.

We quantify the inside-out growth of the Milky Way's low-alpha stellar disk, modelling the ages, metallicities and Galactocentric radii of APOGEE red clump stars with 6 < R < 13 kpc. The current stellar distribution differs significantly from that expected from the star formation history due to the redistribution of stars through radial orbit mixing. We propose and fit a global model for the Milky Way disk, specified by an inside-out star formation history, radial orbit mixing, and an empirical, parametric model for its chemical evolution. We account for the spatially complex survey selection function, and find that the model fits all data well. We find distinct inside-out growth of the Milky Way disk; the best fit model implies that the half-mass radius of the Milky Way disk has grown by 43\% over the last 7 Gyr. Yet, such inside-out growth still results in present-day age gradient weaker than 0.1 Gyr/kpc. Our model predicts the half-mass and half-light sizes of the Galactic disk at earlier epochs, which can be compared to the observed redshift -size relations of disk galaxies. We show that radial orbit migration can reconcile the distinct disk-size evolution with redshift, also expected from cosmological simulations, with the modest present-day age gradients seen in the Milky Way and other galaxies.

Orbital properties of stars, computed from their six-dimensional phase space measurements and an assumed Galactic potential, are used to understand the structure and evolution of the Galaxy. Stellar actions, computed from orbits, have the attractive quality of being invariant under certain assumptions and are therefore used as quantitative labels of a star's orbit. We report a subtle but important systematic error that is induced in the actions as a consequence of local midplane variations expected for the Milky Way. This error is difficult to model because it is non-Gaussian and bimodal, with neither mode peaking on the null value. An offset in the vertical position of the Galactic midplane of $\sim15\,\text{pc}$ for a thin disk-like orbit or $\sim 120\,\text{pc}$ for a thick disk-like orbit induces a $25\%$ systematic error in the vertical action $J_z$. In FIRE simulations of Milky Way-mass galaxies, these variations are on the order of $\sim100\,\text{pc}$ at the solar circle. From observations of the mean vertical velocity variation of $\sim5\text{--}10\,\text{km}\,\text{s}^{-1}$ with radius, we estimate that the Milky Way midplane variations are $\sim60\text{--}170\,\text{pc}$, consistent with three-dimensional dust maps. Action calculations and orbit integrations, which assume the global and local midplanes are identical, are likely to include this induced error, depending on the volume considered. Variation in the local standard of rest or distance to the Galactic center causes similar issues. The variation of the midplane must be taken into account when performing dynamical analysis across the large regions of the disk accessible to Gaia and future missions.

We explore the structure of the element abundance--age--orbit distribution of the stars in the Milky Way's low-$\alpha$ disk, by (re-)deriving precise [Fe/H], [X/Fe] and ages, along with orbits, for red clump stars from the APOGEE survey. There has been a long-standing theoretical expectation and observational evidence that metallicity ([Fe/H]) and age are informative about a star's orbit, e.g. about its angular momentum and the corresponding mean Galactocentric distance or its vertical motion. Indeed, our analysis of the APOGEE data confirms that [Fe/H] or age alone can predict the stars' orbits far less well than the combination of the two. Remarkably, we find and show explicitly, that for known [Fe/H] and age, the other abundances [X/Fe] of Galactic disk stars can be predicted well (on average to 0.02 dex) across a wide range of Galactocentric radii, and therefore provide little additional information, e.g. for predicting their orbit. While the age-abundance space for metal poor stars and potentially for stars near the Galactic center is rich or complex, for the bulk of the Galaxy's low-$\alpha$ disk it is simple: [Fe/H] and age contain most information, unless [X/Fe] can be measured to 0.02, or better. Consequently, we do not have the precision with current (and likely near-future) data to assign stars to their individual (coeval) birth clusters, from which the disk is presumably formed. We can, however, place strong constraints on future models of galactic evolution, chemical enrichment and mixing.

The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.

Chemistry and kinematic studies can determine the origins of stellar population across the Milky Way. The metallicity distribution function of the bulge indicates that it comprises multiple populations, the more metal-poor end of which is particularly poorly understood. It is currently unknown if metal-poor bulge stars ([Fe/H] $<$ -1 dex) are part of the stellar halo in the inner most region, or a distinct bulge population or a combination of these. Cosmological simulations also indicate that the metal-poor bulge stars may be the oldest stars in the Galaxy. In this study, we successfully target metal-poor bulge stars selected using SkyMapper photometry. We determine the stellar parameters of 26 stars and their elemental abundances for 22 elements using R$\sim$ 47,000 VLT/UVES spectra and contrast their elemental properties with that of other Galactic stellar populations. We find that the elemental abundances we derive for our metal-poor bulge stars have much lower overall scatter than typically found in the halo. This indicates that these stars may be a distinct population confined to the bulge. If these stars are, alternatively, part of the inner-most distribution of the halo, this indicates that the halo is more chemically homogeneous at small Galactic radii than at large radii. We also find two stars whose chemistry is consistent with second-generation globular cluster stars. This paper is the first part of the Chemical Origins of Metal-poor Bulge Stars (COMBS) survey that will chemo-dynamically characterize the metal-poor bulge population.

This work explores the detailed chemistry of the Milky Way bulge using the HERMES spectrograph on the Anglo-Australian Telescope. Here we present the abundance ratios of 13 elements for 832 red giant branch and clump stars along the minor bulge axis at latitudes $b=-10^{\circ}, -7.5$ and $-5^{\circ}$. Our results show that none of the abundance ratios vary significantly with latitude. We also observe {\color{red}disk-like} [Na/Fe] abundance ratios, which indicates the bulge does not contain helium-enhanced populations as observed in some globular clusters. Helium enhancement is therefore not the likely explanation for the double red-clump observed in the bulge. We confirm that bulge stars mostly follow abundance trends observed in the disk. However, this similarity is not confirmed across for all elements and metallicity regimes. The more metal-poor bulge population at [Fe/H] $\lesssim -0.8$ is enhanced in the elements associated with core collapse supernovae (SNeII). In addition, the [La/Eu] abundance ratio suggests higher $r$-process contribution, and likely higher star formation in the bulge compared to the disk. This highlights the complex evolution in the bulge, which should be investigated further, both in terms of modelling; and with additional observations of the inner Galaxy.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

To better understand the origin and evolution of the Milky Way bulge, we have conducted a survey of bulge red giant branch and clump stars using the HERMES spectrograph on the Anglo-Australian Telescope. We targeted ARGOS survey stars with pre-determined bulge memberships, covering the full metallicity distribution function. The spectra have signal-to-noise ratios comparable to, and were analysed using the same methods as the GALAH survey. In this work, we present the survey design, stellar parameters, distribution of metallicity and alpha-element abundances along the minor bulge axis at latitudes $b$ = $-10^{\circ}, -7.5^{\circ}$ and $-5^{\circ}$. Our analysis of ARGOS stars indicates that the centroids of ARGOS metallicity components should be located $\approx$0.09 dex closer together. The vertical distribution of $\alpha$-element abundances is consistent with the varying contributions of the different metallicity components. Closer to the plane, alpha abundance ratios are lower as the metal-rich population dominates. At higher latitudes, the alpha abundance ratios increase as the number of metal-poor stars increases. However, we find that the trend of alpha-enrichment with respect to metallicity is independent of latitude. Comparison of our results with those of GALAH DR2 revealed that for [Fe/H] $\approx -0.8$, the bulge shares the same abundance trend as the high-$\alpha$ disk population. However, the metal-poor bulge population ([Fe/H] $\lesssim -0.8$) show enhanced alpha abundance ratios compared to the disk/halo. These observations point to fairly rapid chemical evolution in the bulge, and that the metal-poor bulge population does not share the same similarity with the disk as the more metal-rich populations.

Asteroseismology is the only observational tool in astronomy that can probe the interiors of stars, and is a benchmark method for deriving fundamental properties of stars and exoplanets. Over the coming decade, space-based and ground-based observations will provide a several order of magnitude increase of solar-like oscillators, as well as a dramatic increase in the number and quality of classical pulsator observations, providing unprecedented possibilities to study stellar physics and galactic stellar populations. In this white paper, we describe key science questions and necessary facilities to continue the asteroseismology revolution into the 2020's.

The stellar populations in the inner kiloparsecs of the Milky Way (MW) show complex kinematical and chemical structures. The origin and evolution of these structures is still under debate. Here we study the central region of a fully cosmological hydrodynamical simulation of a disk galaxy that reproduces key properties of the inner kiloparsecs of the MW: it has a boxy morphology and shows an overall rotation and dispersion profile in agreement with observations. We use a clustering algorithm on stellar kinematics to identify a number of discrete kinematic components: a high- and low-spin disk, a stellar halo and two bulge components; one fast rotating and one slow-rotating. We focus on the two bulge components and show that the slow rotating one is spherically symmetric while the fast rotating component shows a boxy/peanut morphology. Although the two bulge components are kinematically discrete populations at present-day, they are both mostly formed over similar time scales, from disk material. We find that stellar particles with lower initial birth angular momentum (most likely thick disc stars) end up in the slow-rotating low-spin bulge, while stars with higher birth angular momentum (most likely thin disc stars) are found in the high-spin bulge. This has the important consequence that a bulge population with a spheroidal morphology does not necessarily indicate a merger origin. In fact, we do find that only $\sim2.3$\% of the stars in the bulge components are ex-situ stars brought in by accreted dwarf galaxies early on. We identify these ex-situ stars as the oldest and most metal-poor stars on highly radial orbits with large vertical excursions from the disk.

There is a dichotomy in the Milky Way in the $[\alpha/$Fe]-[Fe/H] plane, in which stars fall into high-$\alpha$, and low-$\alpha$ sequences. The high-$\alpha$ sequence comprises mostly old stars, and the low-$\alpha$ sequence comprises primarily young stars. The origin of this dichotomy is uncertain. To better understand how the high- and low-$\alpha$ stars are affiliated, we examine if the high- and low-$\alpha$ sequences have distinct orbits at all ages, or if age sets the orbital properties of stars irrespective of their $\alpha$-enhancement. Orbital actions $J_R$, $J_z$, and $J_\phi$ (or $L_z$) are our labels of stellar dynamics. We use ages for 58,278 LAMOST stars (measured to a precision of 40\%) within $\leq$2kpc of the Sun and we calculate orbital actions from proper motions and parallaxes given by Gaia's DR2. We find that \emph{at all ages}, the high- and low-$\alpha$ sequences are dynamically distinct. This implies separate formation and evolutionary histories for the two sequences; a star's membership in the high- or low-$\alpha$ sequence indicates its dynamical properties at a given time. We use action space to make an efficient selection of halo stars and subsequently report a group of old, low-$\alpha$ stars in the halo, which may be a discrete population from an infall event.

We use the second data releases of the ESA Gaia astrometric survey and the high-resolution GALAH spectroscopic survey to analyse the structure of our Galaxy's disc components. With GALAH, we separate the alpha-rich and alpha-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars in the zV_z phase plane, for both V_phi and V_R, and recover the remarkable "phase spiral" discovered by Gaia. We identify the anticipated quadrupole signature in zV_z of a tilted velocity ellipsoid for stars above and below the Galactic plane. By connecting our work with earlier studies, we show that the phase spiral is likely to extend well beyond the narrow solar neighbourhood cylinder in which it was found. The phase spiral is a signature of corrugated waves that propagate through the disc, and the associated non-equilibrium phase mixing. The radially asymmetric distribution of stars involved in the phase spiral reveals that the corrugation, which is mostly confined to the alpha-poor disc, grows in z-amplitude with increasing radius. We present new simulations of tidal disturbance of the Galactic disc by the Sagittarius (Sgr) dwarf. The effect on the zV_z phase plane lasts >2 Gyr but a subsequent disc crossing wipes out the coherent structure. We find that the phase spiral was excited < 0.5 Gyr ago by an object like Sgr with total mass 3 x 10^10 Msun (stripped down from 5 x 10^10 Msun when it first entered the halo) passing through the plane.

The disc structure of the Milky Way is marked by a chemical dichotomy, with high-alpha and low-alpha abundance sequences, traditionally identified with the geometric thick and thin discs. This identification is aided by the old ages of the high-alpha stars, and lower average ages of the low-alpha ones. Recent large scale surveys such as APOGEE have provided a wealth of data on this chemical structure, including showing that an identification of chemical and geometric thick discs is not exact, but the origin of the chemical dichotomy has remained unclear. Here we demonstrate that a dichotomy arises naturally if the early gas-rich disc fragments, leading to some fraction of the star formation occurring in clumps of the type observed in high-redshift galaxies. These clumps have high star formation rate density. They, therefore, enrich rapidly, moving from the low-alpha to the high-alpha sequence, while more distributed star formation produces the low-alpha sequence. We demonstrate that this model produces a chemically-defined thick disc that has many of the properties of the Milky Way's thick disc. Because clump formation is common in high redshift galaxies, we predict that chemical bimodalities are common in massive galaxies.

The orbital properties of stars in the disk are signatures of their formation, but they are also expected to change over time due to the dynamical evolution of the Galaxy. Stellar orbits can be quantified by three dynamical actions, J_r, L_z, and J_z, which provide measures of the orbital eccentricity, guiding radius, and non-planarity, respectively. Changes in these dynamical actions over time reflect the strength and efficiency of the evolutionary processes that drive stellar redistributions. We examine how dynamical actions of stars are correlated with their age using two samples of stars with well-determined ages: 78 solar twin stars (with ages to ~5%) and 4376 stars from the APOKASC2 sample (~20%). We compute actions using spectroscopic radial velocities from previous surveys and parallax and proper motion measurements from Gaia DR2. We find weak gradients in all actions with stellar age, of (7.51 +/- 0.52, -29.0 +/- 1.83, 1.54 +/- 0.18) kpc km/s/Gyr for J_r, L_z, and J_z, respectively. There is, however, significant scatter in the action-age relation. We caution that our results will be affected by the restricted spatial extent of our sample, particularly in the case of J_z. Nevertheless, these action-age gradients and their associated variances provide strong constraints on the efficiency of the mechanisms that drive the redistribution of stellar orbits over time and demonstrate that actions are informative as to stellar age. The shallow action-age gradients combined with the large dispersion in each action at a given age, however, renders the prospect of age inference from orbits of individual stars bleak. Using the precision measurements of [Fe/H] and [$\alpha$/Fe] for our stars we investigate the abundance-action relationship and find weak correlations. Similar to our stellar age results, dynamical actions afford little discriminating power between low- and high-$\alpha$ stars.

We measure the circular velocity curve $v_{\rm c}(R)$ of the Milky Way with the highest precision to date across Galactocentric distances of $5\leq R \leq 25$ kpc. Our analysis draws on the $6$-dimensional phase-space coordinates of $\gtrsim 23,000$ luminous red-giant stars, for which we previously determined precise parallaxes using a data-driven model that combines spectral data from APOGEE with photometric information from WISE, 2MASS, and Gaia. We derive the circular velocity curve with the Jeans equation assuming an axisymmetric gravitational potential. At the location of the Sun we determine the circular velocity with its formal uncertainty to be $v_{\rm c}(R_{\odot}) = (229.0\pm0.2)\rm\,km\,s^{-1}$ with systematic uncertainties at the $\sim 2-5\%$ level. We find that the velocity curve is gently but significantly declining at $(-1.7\pm0.1)\rm\,km\,s^{-1}\,kpc^{-1}$, with a systematic uncertainty of $0.46\rm\,km\,s^{-1}\,kpc^{-1}$, beyond the inner $5$ kpc. We exclude the inner $5$ kpc from our analysis due to the presence of the Galactic bar, which strongly influences the kinematic structure and requires modeling in a non-axisymmetric potential. Combining our results with external measurements of the mass distribution for the baryonic components of the Milky Way from other studies, we estimate the Galaxy's dark halo mass within the virial radius to be $M_{\rm vir} = (7.25\pm0.26)\cdot 10^{11}M_{\odot}$ and a local dark matter density of $\rho_{\rm dm}(R_{\odot}) = 0.30\pm0.03\,\rm GeV\,cm^{-3}$.

The overlap between the spectroscopic Galactic Archaeology with HERMES (GALAH) survey & $Gaia$ provides a high-dimensional chemodynamical space of unprecedented size. We present a first analysis of a subset of this overlap, of 7066 dwarf, turn-off, & sub-giant stars. [...] We investigate correlations between chemical compositions, ages, & kinematics for this sample. Stellar parameters & elemental abundances are derived from the GALAH spectra with the spectral synthesis code SME. [...] We report Li, C, O, Na, Mg, Al, Si, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, as well as Ba & we note that we employ non-LTE calculations for Li, O, Al, & Fe. We show that the use of astrometric & photometric data improves the accuracy of the derived spectroscopic parameters, especially $\log g$. [...] we recover the result that stars of the high-$\alpha$ sequence are typically older than stars in the low-$\alpha$ sequence, the latter spanning $-0.7<$[Fe/H]$<+0.5$. While these two sequences become indistinguishable in [$\alpha$/Fe] vs. [Fe/H] at the metal-rich regime, we find that age can be used to separate stars from the extended high-$\alpha$ & the low-$\alpha$ sequence even in this regime. [...] we find that the old stars ($>8$ Gyr have lower angular momenta $L_z$ than the Sun, which implies that they are on eccentric orbits & originate from the inner disk. Contrary to some previous smaller scale studies we find a continuous evolution in the high-$\alpha$-sequence up to super-solar [Fe/H] rather than a gap, which has been interpreted as a separate "high-$\alpha$ metal-rich" population. Stars in our sample that are younger than 10 Gyr, are mainly found on the low $\alpha$-sequence & show a gradient in $L_z$ from low [Fe/H] ($L_z>L_{z,\odot}$) towards higher [Fe/H] ($L_z<L_{z,\odot}$), which implies that the stars at the ends of this sequence are likely not originating from the close solar vicinity.

GALAH is a large-scale magnitude-limited southern stellar spectroscopic survey. Its second data release (GALAH DR2) provides values of stellar parameters and abundances of 23 elements for 342,682 stars (Buder et al.). Here we add a description of the public release of radial velocities with a typical accuracy of 0.1 km/s for 336,215 of these stars, achievable due to the large wavelength coverage, high resolving power and good signal to noise ratio of the observed spectra, but also because convective motions in stellar atmosphere and gravitational redshift from the star to the observer are taken into account. In the process we derive medians of observed spectra which are nearly noiseless, as they are obtained from between 100 and 1116 observed spectra belonging to the same bin with a width of 50 K in temperature, 0.2 dex in gravity, and 0.1 dex in metallicity. Publicly released 1181 median spectra have a resolving power of 28,000 and trace the well-populated stellar types with metallicities between -0.6 and +0.3. Note that radial velocities from GALAH are an excellent match to the accuracy of velocity components along the sky plane derived by Gaia for the same stars. The level of accuracy achieved here is adequate for studies of dynamics within stellar clusters, associations and streams in the Galaxy. So it may be relevant for studies of the distribution of dark matter.

Stars with unusual elemental abundances offer clues about rare astrophysical events or nucleosynthetic pathways. Stars with significantly depleted magnesium and enhanced potassium ([Mg/Fe] < -0.5; [K/Fe] > 1) have to date only been found in the massive globular cluster NGC 2419 and, to a lesser extent, NGC 2808. The origin of this abundance signature remains unknown, as does the reason for its apparent exclusivity to these two globular clusters. Here we present 112 field stars, identified from 454,180 LAMOST giants, that show significantly enhanced [K/Fe] and possibly depleted [Mg/Fe] abundance ratios. Our sample spans a wide range of metallicities (-1.5 < [Fe/H] < 0.3), yet none show abundance ratios of [K/Fe] or [Mg/Fe] that are as extreme as those observed in NGC 2419. If confirmed, the identified sample of stars represents evidence that the nucleosynthetic process producing the anomalous abundances ratios of [K/Fe] and [Mg/Fe] probably occurs at a wide range of metallicities. This would suggest that pollution scenarios that are limited to early epochs (such as Population III supernovae) are an unlikely explanation, although they cannot be ruled out entirely. This sample is expected to help guide modelling attempts to explain the origin of the Mg-K abundance signature.

We present a study using the second data release of the GALAH survey of stellar parameters and elemental abundances of 15 pairs of stars identified by Oh et al 2017. They identified these pairs as potentially co-moving pairs using proper motions and parallaxes from Gaia DR1. We find that 11 very wide (>1.7 pc) pairs of stars do in fact have similar Galactic orbits, while a further four claimed co-moving pairs are not truly co-orbiting. Eight of the 11 co-orbiting pairs have reliable stellar parameters and abundances, and we find that three of those are quite similar in their abundance patterns, while five have significant [Fe/H] differences. For the latter, this indicates that they could be co-orbiting because of the general dynamical coldness of the thin disc, or perhaps resonances induced by the Galaxy, rather than a shared formation site. Stars such as these, wide binaries, debris of past star formation episodes, and coincidental co-orbiters, are crucial for exploring the limits of chemical tagging in the Milky Way.

We show for the first time, that a fully cosmological hydrodynamical simulation can reproduce key properties of the innermost region of the Milky Way. Our high resolution simulation matches the profile and kinematics of the Milky Way's boxy/peanut-shaped bulge, and hence we can use it to reconstruct and understand the bulge assembly. In particular, the age dependence of the X-shape morphology of the simulated bulge parallels the observed metallicity dependent split in the red clump stars of the inner Galaxy. We use this feature to derive an observational metric that allows us to quantify when the bulge formed from the disk. The metric we propose can be employed with upcoming survey data to constrain the age of the Milky Way bar. From the split in stellar counts we estimate the formation of the 4~kpc scale bar in the simulation to have happened $t^{\rm bar}_{\rm form}\sim8^{+2}_{-2}$ Gyr ago, in good agreement with conventional methods to measure bar formation in simulations. We test the prospects for observationally differentiating the stars that belong to the bulge/bar compared to the surrounding disk, and find that the inner disk and bulge are practically indistinguishable in both chemistry and ages.

We develop and apply a model to quantify the global efficiency of radial orbit migration among stars in the Milky Way disk. This model parameterizes the possible star formation and enrichment histories, radial birth profiles, and combines them with a migration model that relates present-day orbital radii to birth radii through a Gaussian probability, broadening with age $\tau$ as $\sigma_\mathrm{RM8}~\sqrt{\tau/8{~\mathrm{Gyr}}}$. Guided by observations, we assume that stars are born with an initially tight age--metallicity relation at given radius, which becomes subsequently scrambled by radial orbit migration, thereby providing a direct observational constraint on radial orbit migration strength $\sigma_\mathrm{RM8}$. We fit this model with MCMC to the observed age--metallicity distribution of low-$\alpha$ red clump stars with Galactocentric radii between 5 and 14 kpc from APOGEE DR12, sidestepping the complex spatial selection function and accounting for the considerable age uncertainties. This simple model reproduces well the observed data, and we find a global (in radius and time) radial orbit migration efficiency in the Milky Way of $\sigma_\mathrm{RM8}=3.6\pm 0.1$ kpc when marginalizing over all other model aspects. This shows that radial orbit migration in the Milky Way's main disk is indeed rather strong, in line with theoretical expectations: stars migrate by about a half-mass radius over the age of the disk. The model finds the Sun's birth radius at $\sim 5.2$ kpc. If such strong radial orbit migration is typical, this mechanism plays indeed an important role in setting the structural regularity of disk galaxies.

It has been demonstrated that the time variability of a star's brightness at different frequencies can be used to infer its surface gravity, radius, mass, and age. With large samples of light curves now available from Kepler and K2, and upcoming surveys like TESS, we wish to quantify the overall information content of this data and identify where the information resides. As a first look into this question we ask which stellar parameters we can predict from the brightness variations in red-giant stars data and to what precision, using a data-driven model. We demonstrate that the long-cadence (30-minute) Kepler light curves for 2000 red-giant stars can be used to predict their $T_{\rm eff}$ and $\log g$. Our inference makes use of a data-driven model of a part of the autocorrelation function (ACF) of the light curve, where we posit a polynomial relationship between stellar parameters and the ACF pixel values. We find that this model, trained using 1000 stars, can be used to recover the temperature $T_{\rm eff}$ to $<$100 K, the surface gravity to $<$ 0.1 dex, and the asteroseismic power-spectrum parameters $\rm \nu_{max}$ and $\rm \Delta{\nu}$ to $<11$ $\mu$Hz and $<0.9$ $\mu$Hz ($\lesssim$ 15\%). We recover $T_{\rm eff}$ from range of time-lags 0.045 $<$ $T_{\rm lag}$ $<$ 370 days and the $\log g$, $\rm \nu_{max}$ and $\rm \Delta{\nu}$ from the range 0.045 $<$ $T_{\rm lag}$ $<$ 35 days. We do not discover any information about stellar metallicity. The information content of the data about each parameter is empirically quantified using this method, enabling comparisons to theoretical expectations about convective granulation.

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.

The Galactic Archaeology with HERMES (GALAH) survey is a large-scale stellar spectroscopic survey of the Milky Way and designed to deliver chemical information complementary to a large number of stars covered by the $Gaia$ mission. We present the GALAH second public data release (GALAH DR2) containing 342,682 stars. For these stars, the GALAH collaboration provides stellar parameters and abundances for up to 23 elements to the community. Here we present the target selection, observation, data reduction and detailed explanation of how the spectra were analysed to estimate stellar parameters and element abundances. For the stellar analysis, we have used a multi-step approach. We use the physics-driven spectrum synthesis of Spectroscopy Made Easy (SME) to derive stellar labels ($T_\mathrm{eff}$, $\log g$, $\mathrm{[Fe/H]}$, $\mathrm{[X/Fe]}$, $v_\mathrm{mic}$, $v \sin i$, $A_{K_S}$) for a representative training set of stars. This information is then propagated to the whole survey with the data-driven method of $The~Cannon$. Special care has been exercised in the spectral synthesis to only consider spectral lines that have reliable atomic input data and are little affected by blending lines. Departures from local thermodynamic equilibrium (LTE) are considered for several key elements, including Li, O, Na, Mg, Al, Si, and Fe, using 1D MARCS stellar atmosphere models. Validation tests including repeat observations, Gaia benchmark stars, open and globular clusters, and K2 asteroseismic targets lend confidence in our methods and results. Combining the GALAH DR2 catalogue with the kinematic information from $Gaia$ will enable a wide range of Galactic Archaeology studies, with unprecedented detail, dimensionality, and scope.

Using data from the GALAH pilot survey, we determine properties of the Galactic thin and thick disks near the solar neighbourhood. The data cover a small range of Galactocentric radius ($7.9 \leq R_\mathrm{GC} \leq 9.5$ kpc), but extend up to 4 kpc in height from the Galactic plane, and several kpc in the direction of Galactic anti-rotation (at longitude $260 ^\circ \leq \ell \leq 280^\circ$). This allows us to reliably measure the vertical density and abundance profiles of the chemically and kinematically defined `thick' and `thin' disks of the Galaxy. The thin disk (low-$\alpha$ population) exhibits a steep negative vertical metallicity gradient, at d[M/H]/d$z=-0.18 \pm 0.01$ dex kpc$^{-1}$, which is broadly consistent with previous studies. In contrast, its vertical $\alpha$-abundance profile is almost flat, with a gradient of d[$\alpha$/M]/d$z$ = $0.008 \pm 0.002$ dex kpc$^{-1}$. The steep vertical metallicity gradient of the low-$\alpha$ population is in agreement with models where radial migration has a major role in the evolution of the thin disk. The thick disk (high-$\alpha$ population) has a weaker vertical metallicity gradient d[M/H]/d$z = -0.058 \pm 0.003$ dex kpc$^{-1}$. The $\alpha$-abundance of the thick disk is nearly constant with height, d[$\alpha$/M]/d$z$ = $0.007 \pm 0.002$ dex kpc$^{-1}$. The negative gradient in metallicity and the small gradient in [$\alpha$/M] indicate that the high-$\alpha$ population experienced a settling phase, but also formed prior to the onset of major SNIa enrichment. We explore the implications of the distinct $\alpha$-enrichments and narrow [$\alpha$/M] range of the sub-populations in the context of thick disk formation.

Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present first results from the K2-HERMES project, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra of up to 360 stars in one exposure. This ongoing project aims to derive self-consistent spectroscopic parameters for about half of K2 target stars. We present complete stellar parameters and isochrone-derived masses and radii for 46 stars hosting 57 K2 candidate planets in Campaigns 1-3. Our revised host-star radii cast severe doubt on three candidate planets: EPIC\,201407812.01, EPIC\,203070421.01, and EPIC\,202843107.01, all of which now have inferred radii well in excess of the largest known inflated Jovian planets.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) provides spectroscopic information of regions of the inner Milky Way inaccessible to optical surveys. We present the first large study of the metallicity distribution of the innermost Galactic regions based on homogeneous measurements from the SDSS Data Release 12 for 7545 red giant stars within 4.5 kpc of the Galactic center, with the goal to shed light on the structure and origin of the Galactic bulge. Stellar metallicities are found, through multiple-Gaussian decompositions, to be distributed in several components indicative of the presence of various stellar populations such as the bar, or the thin and the thick disk. A super-solar ([Fe/H]=+0.32) and a solar ([Fe/H]=+0.00) metallicity components, tentatively associated with the thin disk and the Galactic bar, respectively, seem to be the major contributors near the midplane. The solar-metallicity component extends outwards in the midplane but is not observed in the innermost regions. The central regions (within 3 kpc of the Galactic center) reveal, on the other hand, the presence of a significant metal-poor population ([Fe/H] =-0.46), tentatively associated with the thick disk, and which becomes the dominant component far from the midplane ($|Z| > +0.75 kpc). Varying contributions from these different components produce a transition region at +0.5 kpc < |Z| < +1.0 kpc characterized by a significant vertical metallicity gradient.

SDSS-V will be an all-sky, multi-epoch spectroscopic survey of over six million objects. It is designed to decode the history of the Milky Way, trace the emergence of the chemical elements, reveal the inner workings of stars, and investigate the origin of planets. It will also create an integral-field spectroscopic map of the gas in the Galaxy and the Local Group that is 1,000x larger than the current state of the art and at high enough spatial resolution to reveal the self-regulation mechanisms of galactic ecosystems. SDSS-V will pioneer systematic, spectroscopic monitoring across the whole sky, revealing changes on timescales from 20 minutes to 20 years. The survey will thus track the flickers, flares, and radical transformations of the most luminous persistent objects in the universe: massive black holes growing at the centers of galaxies. The scope and flexibility of SDSS-V will be unique among extant and future spectroscopic surveys: it is all-sky, with matched survey infrastructures in both hemispheres; it provides near-IR and optical multi-object fiber spectroscopy that is rapidly reconfigurable to serve high target densities, targets of opportunity, and time-domain monitoring; and it provides optical, ultra-wide-field integral field spectroscopy. SDSS-V, with its programs anticipated to start in 2020, will be well-timed to multiply the scientific output from major space missions (e.g., TESS, Gaia, eROSITA) and ground-based projects. SDSS-V builds on the 25-year heritage of SDSS's advances in data analysis, collaboration infrastructure, and product deliverables. The project is now refining its science scope, optimizing the survey strategies, and developing new hardware that builds on the SDSS-IV infrastructure. We present here an overview of the current state of these developments as we seek to build our worldwide consortium of institutional and individual members.

In Debattista et al. (2015), we proposed that a kiloparsec-scale nuclear disc is responsible for the high-velocity secondary peak in the stellar line-of-sight velocity distributions (LOSVDs) seen at positive longitudes in the bulge by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). Here, we make further qualitative but distinctive predictions of the kinematic properties of a nuclear disc, including for the LOSVDs at negative longitudes (which APOGEE-2 will observe) and examine the proper motions throughout the disc. Since a nuclear ring is also able to produce similar high-velocity LOSVD peaks, we present predictions for the proper motion signatures which distinguish between a nuclear disc and a nuclear ring. We also demonstrate that the stars in a nuclear disc, which would be on x2 orbits perpendicular to the bar, can remain on these orbits for a long time and can therefore be old. We show that such (old) nuclear discs of comparable size exist in external galaxies.

The Transiting Exoplanet Survey Satellite (TESS) will provide high precision time-series photometry for millions of stars with at least a half-hour cadence. Of particular interest are the circular regions of 12-degree radius centered around the ecliptic poles that will be observed continuously for a full year. Spectroscopic stellar parameters are desirable to characterize and select suitable targets for TESS, whether they are focused on exploring exoplanets, stellar astrophysics, or Galactic archaeology. Here, we present spectroscopic stellar parameters ($T_{\rm eff}$, $\log g$, [Fe/H], $v \sin i$, $v_{\rm micro}$) for about 16,000 dwarf and subgiant stars in TESS' southern continuous viewing zone. For almost all the stars, we also present Bayesian estimates of stellar properties including distance, extinction, mass, radius, and age using theoretical isochrones. Stellar surface gravity and radius are made available for an additional set of roughly 8,500 red giants. All our target stars are in the range $10<V<13.1$. Among them, we identify and list 227 stars belonging to the Large Magellanic Cloud. The data were taken using the the High Efficiency and Resolution Multi-Element Spectrograph (HERMES, R $\sim 28,000$) at the Anglo-Australian Telescope as part of the TESS-HERMES survey. Comparing our results with the TESS Input Catalog (TIC) shows that the TIC is generally efficient in separating dwarfs and giants, but it has flagged more than hundred cool dwarfs ($T_{\rm eff}< 4800$ K) as giants, which ought to be high-priority targets for the exoplanet search. The catalog can be accessed via http://www.physics.usyd.edu.au/tess-hermes/ , or at MAST via https://archive.stsci.edu/prepds/tess-hermes/ .

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.

APOGEE-2 is a high-resolution, near-infrared spectroscopic survey observing roughly 300,000 stars across the entire sky. It is the successor to APOGEE and is part of the Sloan Digital Sky Survey IV (SDSS-IV). APOGEE-2 is expanding upon APOGEE's goals of addressing critical questions of stellar astrophysics, stellar populations, and Galactic chemodynamical evolution using (1) an enhanced set of target types and (2) a second spectrograph at Las Campanas Observatory in Chile. APOGEE-2 is targeting red giant branch (RGB) and red clump (RC) stars, RR Lyrae, low-mass dwarf stars, young stellar objects, and numerous other Milky Way and Local Group sources across the entire sky from both hemispheres. In this paper, we describe the APOGEE-2 observational design, target selection catalogs and algorithms, and the targeting-related documentation included in the SDSS data releases.

The technique of chemical tagging uses the elemental abundances of stellar atmospheres to `reconstruct' chemically homogeneous star clusters that have long since dispersed. The GALAH spectroscopic survey --which aims to observe one million stars using the Anglo-Australian Telescope -- allows us to measure up to 30 elements or dimensions in the stellar chemical abundance space, many of which are not independent. How to find clustering reliably in a noisy high-dimensional space is a difficult problem that remains largely unsolved. Here we explore t-distributed stochastic neighbour embedding (t-SNE) -- which identifies an optimal mapping of a high-dimensional space into fewer dimensions -- whilst conserving the original clustering information. Typically, the projection is made to a 2D space to aid recognition of clusters by eye. We show that this method is a reliable tool for chemical tagging because it can: (i) resolve clustering in chemical space alone, (ii) recover known open and globular clusters with high efficiency and low contamination, and (iii) relate field stars to known clusters. t-SNE also provides a useful visualization of a high-dimensional space. We demonstrate the method on a dataset of 13 abundances measured in the spectra of 187,000 stars by the GALAH survey. We recover 7 of the 9 observed clusters (6 globular and 3 open clusters) in chemical space with minimal contamination from field stars and low numbers of outliers. With chemical tagging, we also identify two Pleiades supercluster members (which we confirm kinematically), one as far as 6$^\circ$ -- one tidal radius away from the cluster centre.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We present the first self-consistent chemodynamical model fitted to reproduce data for the galactic bulge, bar and inner disk. We extend the Made-to-Measure method to an augmented phase-space including the metallicity of stars, and show its first application to the bar region of the Milky Way. Using data from the ARGOS and APOGEE (DR12) surveys, we adapt the recent dynamical model from Portail et al. to reproduce the observed spatial and kinematic variations as a function of metallicity, thus allowing the detailed study of the 3D density distributions, kinematics and orbital structure of stars in different metallicity bins. We find that metal-rich stars with [Fe/H] > -0.5 are strongly barred and have dynamical properties that are consistent with a common disk origin. Metal-poor stars with [Fe/H] < -0.5 show strong kinematic variations with metallicity, indicating varying contributions from the underlying stellar populations. Outside the central kpc, metal-poor stars are found to have the density and kinematics of a thick disk while in the inner kpc, evidence for an extra concentration of metal-poor stars is found. Finally, the combined orbit distributions of all metallicities in the model naturally reproduce the observed vertex deviations in the bulge. This paper demonstrates the power of Made-to-Measure chemodynamical models, that when extended to other chemical dimensions will be very powerful tools to maximize the information obtained from large spectroscopic surveys such as APOGEE, GALAH and MOONS.

We measure carbon and nitrogen abundances to $\lesssim$ 0.1 dex for 450,000 giant stars from their low-resolution (R$\sim$1800) LAMOST DR2 survey spectra. We use these [C/M] and [N/M] measurements, together with empirical relations based on the APOKASC sample, to infer stellar masses and implied ages for 230,000 of these objects to 0.08 dex and 0.2 dex respectively. We use The Cannon, a data-driven approach to spectral modeling, to construct a predictive model for LAMOST spectra. Our reference set comprises 8125 stars observed in common between the APOGEE and LAMOST surveys, taking seven APOGEE DR12 labels (parameters) as ground truth: Teff, logg, [M/H], [$\alpha$/M], [C/M], [N/M], and Ak. We add seven colors to the Cannon model, based on the g, r, i, J, H, K, W1, and W2 magnitudes from APASS, 2MASS & WISE, which improves our constraints on Teff and logg by up to 20% and on Ak by up to 70%. Cross-validation of the model demonstrates that, for high-SNR objects, our inferred labels agree with the APOGEE values to within 50 K in temperature, 0.04 magnitudes in Ak, and < 0.1 dex in logg, [M/H], [C/M], [N/M], and [$\alpha$/M]. We apply the model to 450,000 giants in LAMOST DR2 that have not been observed by APOGEE. This demonstrates that precise individual abundances can be measured from low-resolution spectra, and represents the largest catalog of [C/M], [N/M], masses and ages to date. As as result, we greatly increase the number and sky coverage of stars with mass and age estimates.

We present a novel interpretation of the previously puzzling different behaviours of stellar populations of the Milky Way's bulge. We first show, by means of pure N-body simulations, that initially co-spatial stellar populations with different in-plane random motions separate when a bar forms. The radially cooler populations form a strong bar, and are vertically thin and peanut-shaped, while the hotter populations form a weaker bar and become a vertically thicker box. We demonstrate that it is the radial, not the vertical, velocity dispersion that dominates this evolution. Assuming that early stellar discs heat rapidly as they form, then both the in-plane and vertical random motions correlate with stellar age and chemistry, leading to different density distributions for metal-rich and metal-poor stars. We then use a high-resolution simulation, in which all stars form out of gas, to demonstrate that this is what happens. When we apply these results to the Milky Way we show that a very broad range of observed trends for ages, densities, kinematics and chemistries, that have been presented as evidence for contradictory paths to the formation of the bulge, are in fact consistent with a bulge which formed from a continuum of disc stellar populations which were kinematically separated by the bar. For the first time we are able to account for the bulge's main trends via a model in which the bulge formed largely in situ. Since the model is generic, we also predict the general appearance of stellar population maps of external edge-on galaxies.

We explore to which extent stars within Galactic disk open clusters resemble each other in the high-dimensional space of their photospheric element abundances, and contrast this with pairs of field stars. Our analysis is based on abundances for 20 elements, homogeneously derived from APOGEE spectra (with carefully quantified uncertainties, with a median value of $\sim 0.03$ dex). We consider 90 red giant stars in seven open clusters and find that most stars within a cluster have abundances in most elements that are indistinguishable (in a $\chi^2$-sense) from those of the other members, as expected for stellar birth siblings. An analogous analysis among pairs of $>1000$ field stars shows that highly significant abundance differences in the 20-dimensional space can be established for the vast majority of these pairs, and that the APOGEE-based abundance measurements have high discriminating power. However, pairs of field stars whose abundances are indistinguishable even at 0.03~dex precision exist: $\sim 0.3$ percent of all field star pairs, and $\sim 1.0$ percent of field star pairs at the same (solar) metallicity [Fe/H] = $0 \pm 0.02$. Most of these pairs are presumably not birth siblings from the same cluster, but rather doppelganger. Our analysis implies that 'chemical tagging' in the strict sense, identifying birth siblings for typical disk stars through their abundance similarity alone, will not work with such data. However, our approach shows that abundances have extremely valuable information for probabilistic chemo-orbital modeling and combined with velocities, we have identified new cluster members from the field.

In this era of large-scale stellar spectroscopic surveys, measurements of stellar attributes ("labels," i.e. parameters and abundances) must be made precise and consistent across surveys. Here, we demonstrate that this can be achieved by a data-driven approach to spectral modeling. With The Cannon, we transfer information from the APOGEE survey to determine precise Teff, log g, [Fe/H], and [$\alpha$/M] from the spectra of 450,000 LAMOST giants. The Cannon fits a predictive model for LAMOST spectra using 9952 stars observed in common between the two surveys, taking five labels from APOGEE DR12 as ground truth: Teff, log g, [Fe/H], [\alpha/M], and K-band extinction $A_k$. The model is then used to infer Teff, log g, [Fe/H], and [$\alpha$/M] for 454,180 giants, 20% of the LAMOST DR2 stellar sample. These are the first [$\alpha$/M] values for the full set of LAMOST giants, and the largest catalog of [$\alpha$/M] for giant stars to date. Furthermore, these labels are by construction on the APOGEE label scale; for spectra with S/N > 50, cross-validation of the model yields typical uncertainties of 70K in Teff, 0.1 in log g, 0.1 in [Fe/H], and 0.04 in [$\alpha$/M], values comparable to the broadly stated, conservative APOGEE DR12 uncertainties. Thus, by using "label transfer" to tie low-resolution (LAMOST R $\sim$ 1800) spectra to the label scale of a much higher-resolution (APOGEE R $\sim$ 22,500) survey, we substantially reduce the inconsistencies between labels measured by the individual survey pipelines. This demonstrates that label transfer with The Cannon can successfully bring different surveys onto the same physical scale.

Bulges of edge-on galaxies are often boxy/peanut-shaped (B/PS), and unsharp masks reveal the presence of an X shape. Simulations show that these shapes can be produced by dynamical processes driven by a bar which vertically thickens the centre. In the Milky Way, which contains such a B/PS bulge, the X-shaped structure is traced by the metal-rich stars but not by the metal-poor ones. Recently Debattista et al. (2016) interpreted this property as a result of the varying effect of the bar on stellar populations with different starting kinematics. This kinematic fractionation model predicts that cooler populations at the time of bar formation go on to trace the X shape, whereas hotter populations are more uniformly distributed. As this prediction is not specific to the Milky Way, we test it with MUSE observations of the B/PS bulge in the nearby galaxy NGC 4710. We show that the metallicity map is more peanut-shaped than the density distribution itself, in good agreement with the prediction. This result indicates that the X-shaped structure in B/PS bulges is formed of relatively metal-rich stars that have been vertically redistributed by the bar, whereas the metal-poor stars have a more uniform, box-shaped distribution.

We construct a large set of dynamical models of the galactic bulge, bar and inner disk using the Made-to-Measure method. Our models are constrained to match the red clump giant density from a combination of the VVV, UKIDSS and 2MASS infrared surveys together with stellar kinematics in the bulge from the BRAVA and OGLE surveys, and in the entire bar region from the ARGOS survey. We are able to recover the bar pattern speed and the stellar and dark matter mass distributions in the bar region, thus recovering the entire galactic effective potential. We find a bar pattern speed of $39.0 \pm 3.5 \,\rm{km\,s^{-1}\,kpc^{-1}}$, placing the bar corotation radius at $6.1 \pm 0.5 \rm{kpc}$ and making the Milky Way bar a typical fast rotator. We evaluate the stellar mass of the long bar and bulge structure to be $M_{\rm{bar/bulge}} = 1.88 \pm 0.12 \times 10^{10} \, \rm{M}_{\odot}$, larger than the mass of disk in the bar region, $M_{\rm{inner\ disk}} = 1.29\pm0.12 \times 10^{10} \, \rm{M}_{\odot}$. The total dynamical mass in the bulge volume is $1.85\pm0.05\times 10^{10} \, \rm{M}_{\odot}$. Thanks to more extended kinematic data sets and recent measurement of the bulge IMF our models have a low dark matter fraction in the bulge of $17\%\pm2\%$. We find a dark matter density profile which flattens to a shallow cusp or core in the bulge region. Finally, we find dynamical evidence for an extra central mass of $\sim0.2\times10^{10} \,\rm{M}_{\odot}$, probably in a nuclear disk or disky pseudobulge.

Much of the inner Milky Way's (MW) global rotation and velocity dispersion patterns can be reproduced by models of secularly-evolved, bar-dominated bulges. More sophisticated constraints, including the higher moments of the line-of-sight velocity distributions (LOSVDs) and limits on the chemodynamical substructure, are critical for interpreting observations of the unresolved inner regions of extragalactic systems and for placing the MW in context with other galaxies. Here, we use SDSS-APOGEE data to develop these constraints, by presenting the first maps of the LOSVD skewness and kurtosis of metal-rich and metal-poor inner MW stars (divided at [Fe/H] = -0.4), and comparing the observed patterns to those that are seen both in N-body models and in extragalactic bars. Despite closely matching the mean velocity and dispersion, the models do not reproduce the observed LOSVD skewness patterns in different ways, which demonstrates that our understanding of the detailed orbital structure of the inner MW remains an important regime for improvement. We find evidence in the MW of the skewness-velocity correlation that is used as a diagnostic of extragalactic bar/bulges. This correlation appears in metal-rich stars only, providing further evidence for different evolutionary histories of chemically differentiated populations. We connect these skewness measurements to previous work on high-velocity "peaks" in the inner Galaxy, confirming the presence of that phenomenon, and we quantify the cylindrical rotation of the inner Galaxy, finding that the latitude-independent rotation vanishes outside of lon ~ 7 deg. Finally, we evaluate the MW data in light of select extragalactic bar diagnostics and discuss progress and challenges of using the MW as a resolved analog of unresolved stellar populations.

The orbits, atmospheric parameters, chemical abundances, and ages of individual stars in the Milky Way provide the most comprehensive illustration of galaxy formation available. The Tycho-Gaia Astrometric Solution (TGAS) will deliver astrometric parameters for the largest ever sample of Milky Way stars, though its full potential cannot be realized without the addition of complementary spectroscopy. Among existing spectroscopic surveys, the RAdial Velocity Experiment (RAVE) has the largest overlap with TGAS ($\gtrsim$200,000 stars). We present a data-driven re-analysis of 520,781 RAVE spectra using The Cannon. For red giants, we build our model using high-fidelity APOGEE stellar parameters and abundances for stars that overlap with RAVE. For main-sequence and sub-giant stars, our model uses stellar parameters from the K2/EPIC. We derive and validate effective temperature $T_{\rm eff}$, surface gravity $\log{g}$, and chemical abundances of up to seven elements (O, Mg, Al, Si, Ca, Fe, Ni). We report a total of 1,685,851 elemental abundances with a typical precision of 0.07 dex, a substantial improvement over previous RAVE data releases. The synthesis of RAVE-on and TGAS is the most powerful data set for chemo-dynamic analyses of the Milky Way ever produced.

The Galactic Archaeology with HERMES (GALAH) Survey is a massive observational project to trace the Milky Way's history of star formation, chemical enrichment, stellar migration and minor mergers. Using high-resolution (R$\simeq$28,000) spectra taken with the High Efficiency and Resolution Multi-Element Spectrograph (HERMES) instrument at the Anglo-Australian Telescope (AAT), GALAH will determine stellar parameters and abundances of up to 29 elements for up to one million stars. Selecting targets from a colour-unbiased catalogue built from 2MASS, APASS and UCAC4 data, we expect to observe dwarfs at 0.3 to 3 kpc and giants at 1 to 10 kpc. This enables a thorough local chemical inventory of the Galactic thin and thick disks, and also captures smaller samples of the bulge and halo. In this paper we present the plan, process and progress as of early 2016 for GALAH survey observations. In our first two years of survey observing we have accumulated the largest high-quality spectroscopic data set at this resolution, over 200,000 stars. We also present the first public GALAH data catalogue: stellar parameters (Teff, log(g), [Fe/H], [alpha/Fe]), radial velocity, distance modulus and reddening for 10680 observations of 9860 Tycho-2 stars that may be included in the first Gaia data release.

In the Milky Way, the thick disk can be defined using individual stellar abundances, kinematics, or age; or geometrically, as stars high above the mid-plane. In nearby galaxies, where only a geometric definition can be used, thick disks appear to have large radial scale-lengths, and their red colors suggest that they are uniformly old. The Milky Way's geometrically thick disk is also radially extended, but it is far from chemically uniform: alpha-enhanced stars are confined within the inner Galaxy. In simulated galaxies, where old stars are centrally concentrated, geometrically thick disks are radially extended, too. Younger stellar populations flare in the simulated disks' outer regions, bringing those stars high above the mid-plane. The resulting geometrically thick disks therefore show a radial age gradient, from old in their central regions to younger in their outskirts. Based on our age estimates for a large sample of giant stars in the APOGEE survey, we can now test this scenario for the Milky Way. We find that the geometrically-defined thick disk in the Milky Way has indeed a strong radial age gradient: the median age for red clump stars goes from ~9 Gyr in the inner disk to 5 Gyr in the outer disk. We propose that at least some nearby galaxies could also have thick disks that are not uniformly old, and that geometrically thick disks might be complex structures resulting from different formation mechanisms in their inner and outer parts.

Chemical tagging promises to use detailed abundance measurements to identify spatially separated stars that were in fact born together (in the same molecular cloud), long ago. This idea has not yielded much practical success, presumably because of the noise and incompleteness in chemical-abundance measurements. We have succeeded in substantially improving spectroscopic measurements with The Cannon, which has now delivered 15 individual abundances for ~100,000 stars observed as part of the APOGEE spectroscopic survey, with precisions around 0.04 dex. We test the chemical-tagging hypothesis by looking at clusters in abundance space and confirming that they are clustered in phase space. We identify (by the k-means algorithm) overdensities of stars in the 15-dimensional chemical-abundance space delivered by The Cannon, and plot the associated stars in phase space. We use only abundance-space information (no positional information) to identify stellar groups. We find that clusters in abundance space are indeed clusters in phase space. We recover some known phase-space clusters and find other interesting structures. This is the first-ever project to identify phase-space structures at survey-scale by blind search purely in abundance space; it verifies the precision of the abundance measurements delivered by The Cannon; the prospects for future data sets appear very good.

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.

Context. Efforts to look for signatures of the first stars have concentrated on metal-poor halo objects. However, the low end of the bulge metallicity distribution has been shown to host some of the oldest objects in the Milky Way and hence this Galactic component potentially offers interesting targets to look at imprints of the first stellar generations. As a pilot project, we selected bulge field stars already identified in the ARGOS survey as having [Fe/H] ~ -1 and oversolar [alpha/Fe] ratios, and we used FLAMES-UVES to obtain detailed abundances of key elements that are believed to reveal imprints of the first stellar generations. Aims. The main purpose of this study is to analyse selected ARGOS stars using new high-resolution (R~45,000) and high-signal-to-noise (S/N >100) spectra. We aim to derive their stellar parameters and elemental ratios, in particular the abundances of C, N, the alpha-elements O, Mg, Si, Ca, and Ti, the odd-Z elements Na and Al, the neutron-capture s-process dominated elements Y, Zr, La, and Ba, and the r-element Eu. Methods. High-resolution spectra of five field giant stars were obtained at the 8m VLT UT2-Kueyen telescope with the UVES spectrograph in FLAMES-UVES configuration. Spectroscopic parameters were derived based on the excitation and ionization equilibrium of Fe I and Fe II. The abundance analysis was performed with a MARCS LTE spherical model atmosphere grid and the Turbospectrum spectrum synthesis code.

We have shown that data-driven models are effective for inferring physical attributes of stars (labels; Teff, logg, [M/H]) from spectra, even when the signal-to-noise ratio is low. Here we explore whether this is possible when the dimensionality of the label space is large (Teff, logg, and 15 abundances: C, N, O, Na, Mg, Al, Si, S, K, Ca, Ti, V, Mn, Fe, Ni) and the model is non-linear in its response to abundance and parameter changes. We adopt ideas from compressed sensing to limit overall model complexity while retaining model freedom. The model is trained with a set of 12,681 red-giant stars with high signal-to-noise spectroscopic observations and stellar parameters and abundances taken from the APOGEE Survey. We find that we can successfully train and use a model with 17 stellar labels. Validation shows that the model does a good job of inferring all 17 labels (typical abundance precision is 0.04 dex), even when we degrade the signal-to-noise by discarding ~50% of the observing time. The model dependencies make sense: the spectral derivatives with respect to abundances correlate with known atomic lines, and we identify elements belonging to atomic lines that were previously unknown. We recover (anti-)correlations in abundance labels for globular cluster stars, consistent with the literature. However we find the intrinsic spread in globular cluster abundances is 3--4 times smaller than previously reported. We deliver 17 labels with associated errors for 87,563 red giant stars, as well as open-source code to extend this work to other spectroscopic surveys.

We used the new ESO VLT instrument MUSE to obtain spectral and imaging coverage of NGC 4710. The wide area and excellent sampling of the MUSE integral field spectrograph allows us to investigate the dynamical properties of the X-shaped bulge of NGC 4710 and compare it with the properties of the Milky Way's own X-shaped bulge. We measured the radial velocities, velocity dispersion, and stellar populations using a penalized pixel full spectral fitting technique adopting simple stellar populations models, on a 1' x 1' area centred on the bulge of NGC 4710. We have constructed the velocity maps of the bulge of NGC 4710 and we investigated the presence of vertical metallicity gradients. These properties were compared to those of the Milky Way bulge and as well as to a simulated galaxy with boxy/peanut bulge. We find the line-of-sight velocity maps and 1D rotation curves of the bulge of NGC 4710 to be remarkably similar to those of the Milky Way bulge. Some specific differences that were identified are in good agreement with the expectations from variations in the bar orientation angle. The bulge of NGC 4710 has a boxy-peanut morphology with a pronounced X-shape, showing no indication of any additional spheroidally distributed bulge population, in which we measure a vertical metallicity gradient of 0.35 dex/kpc. The general properties of NGC 4710 are very similar to those observed in the MW bulge. However, it has been suggested that the MW bulge has an additional component that is comprised of the oldest, most metal-poor stars and which is not part of the boxy-peanut bulge structure. Such a population is not observed in NGC 4710, but could be hidden in the integrated light we observed.

The Milky Way bulge has a boxy/peanut morphology and an X-shaped structure. This X-shape has been revealed by the `split in the red clump' from star counts along the line of sight toward the bulge, measured from photometric surveys. This boxy, X-shaped bulge morphology is not unique to the Milky Way and such bulges are observed in other barred spiral galaxies. N-body simulations show that boxy and X-shaped bulges are formed from the disk via dynamical instabilities. It has also been proposed that the Milky Way bulge is not X-shaped, but rather, the apparent split in the red clump stars is a consequence of different stellar populations, in an old classical spheroidal bulge. We present a WISE image of the Milky Way bulge, produced by downsampling the publicly available "unWISE" coadds. The WISE image of the Milky Way bulge shows that the X-shaped nature of the Milky Way bulge is self-evident and irrefutable. The X-shape morphology of the bulge in itself and the fraction of bulge stars that comprise orbits within this structure has important implications for the formation history of the Milky Way, and, given the ubiquity of boxy X-shaped bulges, spiral galaxies in general.

We present the stellar kinematics across the Galactic bulge and into the disk at positive longitudes from the SDSS-III APOGEE spectroscopic survey of the Milky Way. APOGEE includes extensive coverage of the stellar populations of the bulge along the mid-plane and near-plane regions. From these data, we have produced kinematic maps of 10,000 stars across longitudes 0 deg < l < 65 deg, and primarily across latitudes of |b| < 5 deg in the bulge region. The APOGEE data reveal that the bulge is cylindrically rotating across all latitudes and is kinematically hottest at the very centre of the bulge, with the smallest gradients in both kinematic and chemical space inside the inner-most region (l,|b|) < (5,5) deg. The results from APOGEE show good agreement with data from other surveys at higher latitudes and a remarkable similarity to the rotation and dispersion maps of barred galaxies viewed edge on. The thin bar that is reported to be present in the inner disk within a narrow latitude range of |b| < 2 deg appears to have a corresponding signature in [Fe/H] and [alpha/Fe]. Stars with [Fe/H] > -0.5 have dispersion and rotation profiles that are similar to that of N-body models of boxy/peanut bulges. There is a smooth kinematic transition from the thin bar and boxy bulge (l,|b|) < (15,12) deg out into the disk for stars with [Fe/H] > -1.0, and the chemodynamics across (l,b) suggests the stars in the inner Galaxy with [Fe/H] > -1.0 have an origin in the disk.

The mass of a star is arguably its most fundamental parameter. For red giant stars, tracers luminous enough to be observed across the Galaxy, mass implies a stellar evolution age. It has proven to be extremely difficult to infer ages and masses directly from red giant spectra using existing methods. From the KEPLER and APOGEE surveys, samples of several thousand stars exist with high-quality spectra and asteroseismic masses. Here we show that from these data we can build a data-driven spectral model using The Cannon, which can determine stellar masses to $\sim$ 0.07 dex from APOGEE DR12 spectra of red giants; these imply age estimates accurate to $\sim$ 0.2 dex (40 percent). We show that The Cannon constrains these ages foremost from spectral regions with CN absorption lines, elements whose surface abundances reflect mass-dependent dredge-up. We deliver an unprecedented catalog of 80,000 giants (including 20,000 red-clump stars) with mass and age estimates, spanning the entire disk (from the Galactic center to R $\sim$ 20 kpc). We show that the age information in the spectra is not simply a corollary of the birth-material abundances [Fe/H] and [$\alpha$/Fe], and that even within a mono-abundance population of stars, there are age variations that vary sensibly with Galactic position. Such stellar age constraints across the Milky Way open up new avenues in Galactic archeology.

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 1,475 red giants with asteroseismic mass estimates from Kepler that also have spectroscopic labels from APOGEE DR12. We assess the accuracy of our model, and find that it predicts stellar masses with fractional r.m.s. errors of about 14% (typically 0.2 Msun). From these masses, we derive ages with r.m.s errors of 40%. 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 disk, and that the relation of age with [alpha/M] and metallicity is broadly consistent with established expectations based on detailed studies of the solar neighbourhood.

The Galactic bulge of the Milky Way is made up of stars with a broad range of metallicity, -3.0 < [Fe/H] < 1 dex. The mean of the Metallicity Distribution Function (MDF) decreases as a function of height z from the plane and, more weakly, with galactic radius. The most metal rich stars in the inner Galaxy are concentrated to the plane and the more metal poor stars are found predominantly further from the plane, with an overall vertical gradient in the mean of the MDF of about -0.45 dex/kpc. This vertical gradient is believed to reflect the changing contribution with height of different populations in the inner-most region of the Galaxy. The more metal rich stars of the bulge are part of the boxy/peanut structure and comprise stars in orbits which trace out the underlying X-shape. There is still a lack of consensus on the origin of the metal poor stars ([Fe/H] < -0.5) in the region of the bulge. Some studies attribute the more metal poor stars of the bulge to the thick disk and stellar halo that are present in the inner region, and other studies propose that the metal poor stars are a distinct "old spheroid" bulge population. Understanding the origin of the populations that make up the MDF of the bulge, and identifying if there is a unique bulge population which has formed separately from the disk and halo, has important consequences for identifying the relevant processes in the the formation and evolution of the Milky Way.

The first stars are predicted to have formed within 200 million years after the Big Bang, initiating the cosmic dawn. A true first star has not yet been discovered, although stars with tiny amounts of elements heavier than helium ('metals') have been found in the outer regions ('halo') of the Milky Way. The first stars and their immediate successors should, however, preferentially be found today in the central regions ('bulges') of galaxies, because they formed in the largest over-densities that grew gravitationally with time. The Milky Way bulge underwent a rapid chemical enrichment during the first 1-2 billion years, leading to a dearth of early, metal-poor stars. Here we report observations of extremely metal-poor stars in the Milky Way bulge, including one star with an iron abundance about 10,000 times lower than the solar value without noticeable carbon enhancement. We confirm that the most metal-poor bulge stars are on tight orbits around the Galactic Centre, rather than being halo stars passing through the bulge, as expected for stars formed at redshifts greater than 15. Their chemical compositions are in general similar to typical halo stars of the same metallicity although intriguing differences exist, including lower abundances of carbon.

The Apache Point Observatory Galactic Evolution Experiment has measured the stellar velocities of red giant stars in the inner Milky Way. We confirm that the line of sight velocity distributions (LOSVDs) in the mid-plane exhibit a second peak at high velocities, whereas those at |b| = 2degrees do not. We use a high resolution simulation of a barred galaxy, which crucially includes gas and star formation, to guide our interpretation of the LOSVDs. We show that the data are fully consistent with the presence of a thin, rapidly rotating, nuclear disk extending to ~1 kpc. This nuclear disk is orientated perpendicular to the bar and is likely to be composed of stars on x2 orbits. The gas in the simulation is able to fall onto such orbits, leading to stars populating an orthogonal disk.

New spectroscopic surveys offer the promise of consistent stellar parameters and abundances ('stellar labels') for hundreds of thousands of stars in the Milky Way: this poses a formidable spectral modeling challenge. In many cases, there is a sub-set of reference objects for which the stellar labels are known with high(er) fidelity. We take advantage of this with The Cannon, a new data-driven approach for determining stellar labels from spectroscopic data. The Cannon learns from the 'known' labels of reference stars how the continuum-normalized spectra depend on these labels by fitting a flexible model at each wavelength; then, The Cannon uses this model to derive labels for the remaining survey stars. We illustrate The Cannon by training the model on only 542 stars in 19 clusters as reference objects, with Teff, log g and [Fe/H] as the labels, and then applying it to the spectra of 56,000 stars from APOGEE DR10. The Cannon is very accurate. Its stellar labels compare well to the stars for which APOGEE pipeline (ASPCAP) labels are provided in DR10, with rms differences that are basically identical to the stated ASPCAP uncertainties. Beyond the reference labels, The Cannon makes no use of stellar models nor any line-list, but needs a set of reference objects that span label-space. The Cannon performs well at lower signal-to-noise, as it delivers comparably good labels even at one ninth the APOGEE observing time. We discuss the limitations of The Cannon and its future potential, particularly, to bring different spectroscopic surveys onto a consistent scale of stellar labels.

The Ophiuchus stream is a recently discovered stellar tidal stream in the Milky Way. We present high-quality spectroscopic data for 14 stream member stars obtained using the Keck and MMT telescopes. We confirm the stream as a fast moving ($v_{los}\sim290$ km s$^{-1}$), kinematically cold group ($\sigma_{v_{los}}\lesssim1$ km s$^{-1}$) of $\alpha$-enhanced and metal-poor stars (${\rm [\alpha/Fe]\sim0.4}$ dex, ${\rm [Fe/H]\sim-2.0}$ dex). Using a probabilistic technique, we model the stream simultaneously in line-of-sight velocity, color-magnitude, coordinate, and proper motion space, and so determine its distribution in 6D phase-space. We find that that the stream extends in distance from 7.5 to 9 kpc from the Sun; it is 50 times longer than wide, merely appearing highly foreshortened in projection. The analysis of the stellar population contained in the stream suggests that it is $\sim12$ Gyr old, and that its initial stellar mass was $\sim2\times10^4$ $M_{\odot}$ (or at least $\gtrsim7\times10^3$ $M_{\odot}$). Assuming a fiducial Milky Way potential, we fit an orbit to the stream which matches the observed phase-space distribution, except for some tension in the proper motions: the stream has an orbital period of $\sim350$ Myr, and is on a fairly eccentric orbit ($e\sim0.66$) with a pericenter of $\sim3.5$ kpc and an apocenter of $\sim17$ kpc. The phase-space structure and stellar population of the stream show that its progenitor must have been a globular cluster that was disrupted only $\sim240$ Myr ago. We do not detect any significant overdensity of stars along the stream that would indicate the presence of a progenitor, and conclude that the stream is all that is left of the progenitor.

We model the split red clump of the Galactic bulge in OGLE-III photometry, and compare the results to predictions from two N-body models. Our analysis yields precise maps of the brightness of the two red clumps, the fraction of stars in the more distant peak, and their combined surface density. We compare the observations to predictions from two N-body models previously used in the literature. Both models correctly predict several features as long as one assumes an angle $\alpha_{\rm{Bar}} \approx 30^{\circ}$ between the Galactic bar's major axis and the line of sight to the Galactic centre. In particular that the fraction of stars in the faint red clump should decrease with increasing longitude. The biggest discrepancies between models and data are in the rate of decline of the combined surface density of red clump stars toward negative longitudes and of the brightness difference between the two red clumps toward positive longitudes, with neither discrepancy exceeding $\sim$25% in amplitude. Our analysis of the red giant luminosity function also yields an estimate of the red giant branch bump parameters toward these high-latitude fields, and evidence for a high rate ($\sim$25%) of disk contamination in the bulge at the colour and magnitude of the red clump, with the disk contamination rate increasing toward sightlines further distant from the plane.

The details of bulge formation via collapse, mergers, secular processes or their interplay remain unresolved. To start answering this question and quantify the importance of distinct mechanisms, we mapped a sample of three galactic bulges using data from the integral field spectrograph WiFeS on the ANU 2.3m telescope in Siding Spring Observatory. Its high resolution gratings (R=7000) allow us to present a detailed kinematic and stellar population analysis of their inner structures with classical and novel techniques. The comparison of those techniques calls for the necessity of inversion algorithms in order to understand complex substructures and separate populations. We use line-strength indices to derive SSP-equivalent ages and metallicities. Additionally, we use full spectral fitting methods, here the code STECKMAP, to extract their star formation histories. The high quality of our data allows us to study the 2D distribution of different stellar populations (i.e. young, intermediate, and old). We can identify their dominant populations based on these age-discriminated 2D light and mass contribution. In all galactic bulges studied, at least 50% of the stellar mass already existed 10 Gyrs ago, more than currently predicted by simulations. A younger component (age between 1 to 8 Gyrs) is also prominent and its present day distribution seems to be affected much more strongly by morphological structures, especially bars, than the older one. This in-depth analysis of the three bulges supports the notion of increasing complexity in their evolution, likely to be found in numerous bulge structures if studied at this level of detail, which cannot be achieved by mergers alone and require a non-negligible contribution of secular evolution.

By analyzing a N-body simulation of a bulge formed simply via a bar instability mechanism operating on a kinematically cold stellar disk, and by comparing the results of this analysis with the structural and kinematic properties of the main stellar populations of the Milky Way bulge, we conclude that the bulge of our Galaxy is not a pure stellar bar formed from a pre-existing thin stellar disk, as some studies have recently suggested. On the basis of several arguments emphasized in this paper, we propose that the bulge population which, in the Milky Way, is observed not to be part of the peanut structure corresponds to the old galactic thick disk, thus implying that the Milky Way is a pure thin+thick disk galaxy, with only a possible limited contribution of a classical bulge.

We present an abundance analysis of eight potential member stars of the old Galactic bulge globular cluster NGC6522. The same stars have previously been studied by Chiappini et al. (2011), who found very high abundances of the slow neutron capture elements compared with other clusters and field stars of similar metallicity, which they interpreted as reflecting nucleosynthesis in rapidly rotating, massive Population III stars. In contrast to their analysis, we do not find any unusual enhancements of the neutron capture elements Sr, Y, Ba and Eu and conclude that previous claims result mainly from not properly accounting for blending lines. Instead we find NGC6522 to be an unremarkable globular cluster with comparable abundance trends to other Galactic globular clusters at the same metallicity ([Fe/H] = -1.15 +/- 0.16). The stars are also chemically similar to halo and bulge field stars at the same metallicity, spanning a small range in [Y/Ba] and with normal {\alpha}-element abundances. We thus find no observational evidence for any chemical signatures of rapidly rotating Population III stars in NGC 6522.

The center of our disk galaxy, the Milky Way, is dominated by a boxy/peanut-shaped bulge. Numerous studies of the bulge based on stellar photometry have concluded that the bulge stars are exclusively old. The perceived lack of young stars in the bulge strongly constrains its likely formation scenarios, providing evidence that the bulge is a unique population that formed early and separately from the disk. However, recent studies of individual bulge stars using the microlensing technique have reported that they span a range of ages, emphasizing that the bulge may not be a monolithic structure. In this letter we demonstrate that the presence of young stars that are located predominantly near the plane is expected for a bulge that has formed from the disk via dynamical instabilities. Using an N-body+SPH simulation of a disk galaxy forming out of gas cooling inside a dark matter halo and forming stars, we find a qualitative agreement between our model and the observations of young metal-rich stars in the bulge. We are also able to partially resolve the apparent contradiction in the literature between results that argue for a purely old bulge population and those which show a population comprised of a range in ages; the key is where to look.

We present the kinematic results from our ARGOS spectroscopic survey of the Galactic bulge of the Milky Way. Our aim is to understand the formation of the Galactic bulge. We examine the kinematics of about 17,400 stars in the bulge located within 3.5 kpc of the Galactic centre, identified from the 28,000 star ARGOS survey. We aim to determine if the formation of the bulge has been internally driven from disk instabilities as suggested by its boxy shape, or if mergers have played a significant role as expected from Lambda CDM simulations. From our velocity measurements across latitudes b = -5 deg, -7.5 deg and -10 deg we find the bulge to be a cylindrically rotating system that transitions smoothly out into the disk. Within the bulge, we find a kinematically distinct metal-poor population ([Fe/H] < -1.0) that is not rotating cylindrically. The 5% of our stars with [Fe/H] < -1.0 are a slowly rotating spheroidal population, which we believe are stars of the metal weak thick disk and halo which presently lie in the inner Galaxy. The kinematics of the two bulge components that we identified in ARGOS paper III (mean [Fe/H] = -0.25 and [Fe/H] = +0.15, respectively) demonstrate that they are likely to share a common formation origin and are distinct from the more metal poor populations of the thick disk and halo which are colocated inside the bulge. We do not exclude an underlying merger generated bulge component but our results favour bulge formation from instabilities in the early thin disk.

We describe the motivation, field locations and stellar selection for the ARGOS spectroscopic survey of 28,000 stars in the bulge and inner disk of the Milky Way galaxy across latitudes of b = -5 deg to -10 deg. The primary goal of this survey is to constrain the formation processes of the bulge and establish whether it is predominantly a merger or instability remnant. From the spectra (R = 11,000), we have measured radial velocities and determined stellar parameters, including metallicities and [alpha/Fe] ratios. Distances were estimated from the derived stellar parameters and about 14,000 stars are red giants within 3.5 kpc of the Galactic centre. In this paper, we present the observations and analysis methods. Subsequent papers (III and IV) will discuss the stellar metallicity distribution and kinematics of the Galactic bulge and inner disk, and the implications for the formation of the bulge.

We present the metallicity results from the ARGOS spectroscopic survey of the Galactic bulge. Our aim is to understand the formation of the Galactic bulge: did it form via mergers, as expected from Lambda CDM theory, or from disk instabilities, as suggested by its boxy/peanut shape, or both? We have obtained spectra for 28,000 stars at a spectral resolution of R = 11,000. From these spectra, we have determined stellar parameters and distances to an accuracy of < 1.5 kpc. The stars in the inner Galaxy span a large range in [Fe/H], -2.8 < [Fe/H] < +0.6. From the spatial distribution of the red clump stars as a function of [Fe/H] (Ness et al. 2012a), we propose that the stars with [Fe/H] > -0.5 are part of the boxy/peanut bar/bulge. We associate the lower metallicity stars ([Fe/H] < -0.5) with the thick disk, which may be puffed up in the inner region, and with the inner regions of the metal-weak thick disk and inner halo. For the bulge stars with [Fe/H] > -0.5, we find two discrete populations; (i) stars with [Fe/H] ~ -0.25 which provide a roughly constant fraction of the stars in the latitude interval b = -5 deg to -10 deg, and (ii) a kinematically colder, more metal-rich population with mean [Fe/H] ~ +0.15 which is more prominent closer to the plane. The changing ratio of these components with latitude appears as a vertical abundance gradient of the bulge. We attribute both of these bulge components to instability-driven bar/bulge formation from the thin disk. We do not exclude a weak underlying classical merger-generated bulge component, but see no obvious kinematic association of any of our bulge stars with such a classical bulge component. [abridged]

Near the minor axis of the Galactic bulge, at latitudes b < -5 degrees, the red giant clump stars are split into two components along the line of sight. We investigate this split using the three fields from the ARGOS survey that lie on the minor axis at (l,b) = (0,-5), (0,-7.5), (0,-10) degrees. The separation is evident for stars with [Fe/H] > -0.5 in the two higher-latitude fields, but not in the field at b = -5 degrees. Stars with [Fe/H] < -0.5 do not show the split. We compare the spatial distribution and kinematics of the clump stars with predictions from an evolutionary N-body model of a bulge that grew from a disk via bar-related instabilities. The density distribution of the peanut-shaped model is depressed near its minor axis. This produces a bimodal distribution of stars along the line of sight through the bulge near its minor axis, very much as seen in our observations. The observed and modelled kinematics of the two groups of stars are also similar. We conclude that the split red clump of the bulge is probably a generic feature of boxy/peanut bulges that grew from disks, and that the disk from which the bulge grew had relatively few stars with [Fe/H] < -0.5