We present MINESweeper, a tool to measure stellar parameters by jointly fitting observed spectra and broadband photometry to model isochrones and spectral libraries. This approach enables the measurement of spectrophotometric distances, in addition to stellar parameters such as Teff, log(g), [Fe/H], [a/Fe], and radial velocity. MINESweeper employs a Bayesian framework and can easily incorporate a variety of priors, including Gaia parallaxes. Mock data are fit in order to demonstrate how the precision of derived parameters depends on evolutionary phase and SNR. We then fit a selection of data in order to validate the model outputs. Fits to a variety of benchmark stars including Procyon, Arcturus, and the Sun result in derived stellar parameters that are in good agreement with the literature. We then fit combined spectra and photometry of stars in the open and globular clusters M92, M13, M3, M107, M71, and M67. Derived distances, [Fe/H], [a/Fe], and log(g)-Teff, relations are in overall good agreement with literature values, although there are trends between metallicity and log(g), within clusters that point to systematic uncertainties at the ~0.1 dex level. Finally, we fit a large sample of stars from the H3 Spectroscopic Survey in which high quality Gaia parallaxes are also available. These stars are fit without the Gaia parallaxes so that the geometric parallaxes can serve as an independent test of the spectrophotometric distances. Comparison between the two reveals good agreement within their formal uncertainties after accounting for the Gaia zero point uncertainties.

We discuss how to achieve mapping from large sets of imperfect simulations and observational data with unsupervised domain adaptation. Under the hypothesis that simulated and observed data distributions share a common underlying representation, we show how it is possible to transfer between simulated and observed domains. Driven by an application to interpret stellar spectroscopic sky surveys, we construct the domain transfer pipeline from two adversarial autoencoders on each domains with a disentangling latent space, and a cycle-consistency constraint. We then construct a differentiable pipeline from physical stellar parameters to realistic observed spectra, aided by a supplementary generative surrogate physics emulator network. We further exemplify the potential of the method on the reconstructed spectra quality and to discover new spectral features associated to elemental abundances.

Spectroscopy provides an immense amount of information on stellar objects, and this field continues to grow with recent developments in multi-object data acquisition and rapid data analysis techniques. Current automated methods for analyzing spectra are either (a) data-driven models, which require large amounts of data with prior knowledge of stellar parameters and elemental abundances, or (b) based on theoretical synthetic models that are susceptible to the gap between theory and practice. In this study, we present a hybrid generative domain adaptation method to turn simulated stellar spectra into realistic spectra, learning from the large spectroscopic surveys. We use a neural network to emulate computationally expensive stellar spectra simulations, and then train a separate unsupervised domain-adaptation network that learns to relate the generated synthetic spectra to observational spectra. Consequently, the network essentially produces data-driven models without the need for a labeled training set. As a proof of concept, two case studies are presented. The first of which is the auto-calibration of synthetic models without using any standard stars. To accomplish this, synthetic models are morphed into spectra that resemble observations, thereby reducing the gap between theory and observations. The second case study is the identification of the elemental source of missing spectral lines in the synthetic modelling. These sources are predicted by interpreting the differences between the domain-adapted and original spectral models. To test our ability to identify missing lines, we use a mock dataset and show that, even with noisy observations, absorption lines can be recovered when they are absent in one of the domains. While we focus on spectral analyses in this study, this method can be applied to other fields, which use large data sets and are currently limited by modelling accuracy.

In the $\Lambda$CDM paradigm the Galactic stellar halo is predicted to harbor the accreted debris of smaller systems. To identify these systems, the H3 Spectroscopic Survey, combined with $Gaia$, is gathering 6D phase-space and chemical information in the distant Galaxy. Here we present a comprehensive inventory of structure within 50 kpc from the Galactic center using a sample of 5684 giants at $|b|>40^{\circ}$ and $|Z|>2$ kpc. We identify known structures including the high-$\alpha$ disk, the in-situ halo (disk stars heated to eccentric orbits), Sagittarius (Sgr), $Gaia$-Sausage-Enceladus (GSE), the Helmi Streams, Sequoia, and Thamnos. Additionally, we identify the following new structures: (i) Aleph ([Fe/H]$=-0.5$), a low eccentricity structure that rises a surprising 10 kpc off the plane, (ii, iii) Arjuna ([Fe/H]$=-1.2$) and I'itoi ([Fe/H]$<-2$), which comprise the high-energy retrograde halo along with Sequoia, and (iv) Wukong ([Fe/H]$=-1.6$), a prograde phase-space overdensity chemically distinct from GSE. For each structure we provide [Fe/H], [$\alpha$/Fe], and orbital parameters. Stars born within the Galaxy are a major component at $|Z|\sim$2 kpc ($\approx$60$\%$), but their relative fraction declines sharply to $\lesssim$5$\%$ past 15 kpc. Beyond 15 kpc, $>$80$\%$ of the halo is built by two massive ($M_{\star}\sim10^{8}-10^{9}M_{\odot}$) accreted dwarfs: GSE ([Fe/H]$=-1.2$) within 25 kpc, and Sgr ([Fe/H]$=-1.0$) beyond 25 kpc. This explains the relatively high overall metallicity of the halo ([Fe/H]$\approx-1.2$). We attribute $\gtrsim$95$\%$ of the sample to one of the listed structures, pointing to a halo built entirely from accreted dwarfs and heating of the disk.

Parameter estimation with non-Gaussian stochastic fields is a common challenge in astrophysics and cosmology. In this paper we advocate performing this task using the scattering transform, a statistical tool rooted in the mathematical properties of convolutional neural nets. This estimator can characterize a complex field without explicitly computing higher-order statistics, thus avoiding the high variance and dimensionality problems. It generates a compact set of coefficients which can be used as robust summary statistics for non-Gaussian information. It is especially suited for fields presenting localized structures and hierarchical clustering, such as the cosmological density field. To demonstrate its power, we apply this estimator to the cosmological parameter inference problem in the context of weak lensing. Using simulated convergence maps with realistic noise, the scattering transform outperforms the power spectrum and peak counts, and is on par with the state-of-the-art CNN. It retains the advantages of traditional statistical descriptors (it does not require any training nor tuning), has provable stability properties, allows to check for systematics, and importantly, the scattering coefficients are interpretable. It is a powerful and attractive estimator for observational cosmology and, in general, the study of physically-motivated fields.

Increasingly powerful and multiplexed spectroscopic facilities promise detailed chemical abundance patterns for millions of resolved stars in galaxies beyond the Milky Way (MW). Here, we employ the Cram\'er-Rao Lower Bound (CRLB) to forecast the precision to which stellar abundances for metal-poor, low-mass stars outside the MW can be measured for 41 current (e.g., Keck, MMT, VLT, DESI) and planned (e.g., MSE, JWST, ELTs) spectrograph configurations. We show that moderate resolution ($R\lesssim5000$) spectroscopy at blue-optical wavelengths ($\lambda\lesssim4500$ \AA) (i) enables the recovery of 2-4 times as many elements as red-optical spectroscopy ($5000\lesssim\lambda\lesssim10000$ \AA) at similar or higher resolutions ($R\sim 10000$) and (ii) can constrain the abundances of several neutron capture elements to $\lesssim$0.3 dex. We further show that high-resolution ($R\gtrsim 20000$), low S/N ($\sim$10 pixel$^{-1}$) spectra contain rich abundance information when modeled with full spectral fitting techniques. We demonstrate that JWST/NIRSpec and ELTs can recover (i) $\sim$10 and 30 elements, respectively, for metal-poor red giants throughout the Local Group and (ii) [Fe/H] and [$\alpha$/Fe] for resolved stars in galaxies out to several Mpc with modest integration times. We show that select literature abundances are within a factor of $\sim$2 (or better) of our CRLBs. We suggest that, like ETCs, CRLBs should be used when planning stellar spectroscopic observations. We include an open source python package, \texttt{Chem-I-Calc}, that allows users to compute CRLBs for spectrographs of their choosing.

Lithium depletion and enrichment in the cosmos is not yet well understood. To help tighten constraints on stellar and Galactic evolution models, we present the largest high-resolution analysis of Li abundances A(Li) to date, with results for over 100 000 GALAH field stars spanning effective temperatures $5900\,\mathrm{K} \lesssim \rm{T_{eff}} \lesssim7000\,\mathrm{K}$ and metallicities $-3 \lesssim \rm[Fe/H] \lesssim +0.5$. We separated these stars into two groups, on the warm and cool side of the so-called Li-dip, a localised region of the Kiel diagram wherein lithium is severely depleted. We discovered that stars in these two groups show similar trends in the A(Li)-[Fe/H] plane, but with a roughly constant offset in A(Li) of 0.4 dex, the warm group having higher Li abundances. At $\rm[Fe/H]\gtrsim-0.5$, a significant increasing in Li abundance with increasing metallicity is evident in both groups, signalling the onset of significant Galactic production. At lower metallicity, stars in the cool group sit on the Spite plateau, showing a reduced lithium of around 0.4 dex relative to the primordial value predicted from Big Bang nucleosynthesis (BBN). However, stars in the warm group between [Fe/H] = -1.0 and -0.5, form an elevated plateau that is largely consistent with the BBN prediction. This may indicate that these stars in fact preserve the primordial Li produced in the early Universe.

We present $\gtrsim 15,000$ metal-rich (${\rm [Fe/H]}>-0.2$dex) A and F stars whose surface abundances deviate strongly from Solar abundance ratios and cannot plausibly reflect their birth material composition. These stars are identified by their high [Ba/Fe] abundance ratios (${\rm [Ba/Fe]}>1.0$dex) in the LAMOST DR5 spectra analyzed by Xiang et al. (2019). They are almost exclusively main sequence and subgiant stars with $T_{\rm eff}\gtrsim6300$K. Their distribution in the Kiel diagram ($T_{\rm eff}$--$\log g$) traces a sharp border at low temperatures along a roughly fixed-mass trajectory (around $1.4M_\odot)$ that corresponds to an upper limit in convective envelope mass fraction of around $10^{-4}$. Most of these stars exhibit distinctly enhanced abundances of iron-peak elements (Cr, Mn, Fe, Ni) but depleted abundances of Mg and Ca. Rotational velocity measurements from GALAH DR2 show that the majority of these stars rotate slower than typical stars in an equivalent temperature range. These characteristics suggest that they are related to the so-called Am/Fm stars. Their abundance patterns are qualitatively consistent with the predictions of stellar evolution models that incorporate radiative acceleration, suggesting they are a consequence of stellar internal evolution particularly involving the competition between gravitational settling and radiative acceleration. These peculiar stars constitute 40% of the whole population of stars with mass above 1.5$M_\odot$, affirming that "peculiar" photospheric abundances due to stellar evolution effects are a ubiquitous phenomenon for these intermediate-mass stars. This large sample of Ba-enhanced chemically peculiar A/F stars with individual element abundances provides the statistics to test more stringently the mechanisms that alter the surface abundances in stars with radiative envelopes.

We present a neural network autoencoder structure that is able to extract essential latent spectral features from observed spectra and then reconstruct a spectrum from those features. Because of the training with a set of unpeculiar spectra, the network is able to reproduce a spectrum of high signal-to-noise ratio that does not show any spectral peculiarities even if they are present in an observed spectrum. Spectra generated in this manner were used to identify various emission features among spectra acquired by multiple surveys using the HERMES spectrograph at the Anglo-Australian telescope. Emission features were identified by a direct comparison of the observed and generated spectra. Using the described comparison procedure, we discovered 10,364 candidate spectra with a varying degree of H$\alpha$/H$\beta$ emission component produced by different physical mechanisms. A fraction of those spectra belongs to the repeated observation that shows temporal variability in their emission profile. Among emission spectra, we find objects that feature contributions of a nearby rarefied gas (identified through the emission of [NII] and [SII] lines) that was identified in 4004 spectra, which were not all identified as having H$\alpha$ emission. Positions of identified emission-line objects coincide with multiple known regions that harbour young stars. Similarly, detected nebular emission spectra coincide with visually-prominent nebular clouds observable in the red all-sky photographic composites.

We investigate the properties of 1306 red giant stars with high photospheric abundances of lithium observed by the GALAH, K2-HERMES and TESS-HERMES surveys, and discuss them in the context of proposed mechanisms for lithium enrichment in giant stars. We confirm that Li-rich giants are rare, making up only 1.1 per cent of our giant star sample. We use GALAH+ DR3 stellar parameters and a Bayesian isochrone analysis to divide the sample into first-ascent red giant branch and red clump stars, and confirm these classifications using asteroseismic data from K2. We find that red clump stars are more than three times as likely to be lithium-rich as red giant branch stars, and the occurrence rate of lithium richness with metallicity is quite different between the two populations. The probability for a rapidly rotating giant to be lithium-rich is distinctly higher as compared to a non rapidly rotating giant, independent of its evolutionary state. There is also a clear correspondence on the red clump between rapid rotation and extreme levels of lithium enrichment. There are almost no stars on the secondary red clump with lithium abundances above the primordial level, indicating a mass dependence in the red clump lithium enrichment mechanism(s). The complex distribution of lithium-rich giants across evolutionary phase, metallicity, rotation rate and mass implies multiple independent mechanisms for producing lithium enrichment.

Globular cluster progenitors may have been detected by \textit{HST}, and are predicted to be observable with \textit{JWST} and ground-based extremely-large telescopes with adaptive optics. This has the potential to elucidate the issue of globular cluster formation and the origins of significantly helium-enriched subpopulations, a problem in Galactic astronomy with no satisfactory theoretical solution. Given this context, we use model stellar tracks and isochrones to investigate the predicted observational properties of helium-enriched stellar populations in globular cluster progenitors. We find that, relative to helium-normal populations, helium-enriched (${\Delta}Y=+0.12$) stellar populations similar to those inferred in the most massive globular clusters, are expected, modulo some rapid fluctuations in the first $\sim$30 Myr, to be brighter and redder in the rest frame. At fixed age, stellar mass, and metallicity, a helium-enriched population is predicted to converge to being $\sim$0.40 mag brighter at $\lambda \approx 2.0\, {\mu}m$, and to be 0.30 mag redder in the \textit{JWST}-NIRCam colour $(F070W-F200W)$, and to actually be fainter for $\lambda \lesssim 0.50 \, {\mu}m$. Separately, we find that the time-integrated shift in ionizing radiation is a negligible $\sim 5\%$, though we show that the Lyman-$\alpha$ escape fraction could end up higher for helium-enriched stars.

Binary stellar systems form a large fraction of the Galaxy's stars. They are useful as laboratories for studying the physical processes taking place within stars, and must be correctly taken into account when observations of stars are used to study the structure and evolution of the Galaxy. We present a sample of 12760 well-characterised double-lined spectroscopic binaries that are appropriate for statistical studies of the binary populations. They were detected as SB2s using a t-distributed stochastic neighbour embedding (t-SNE) classification and a cross-correlation analysis of GALAH spectra. This sample consists mostly of dwarfs, with a significant fraction of evolved stars and several dozen members of the giant branch. To compute parameters of the primary and secondary star ($T_{\rm eff[1,2]}$, $\log g_{[1,2]}$, [Fe/H], $V_{r[1,2]}$, $v_{\rm mic[1,2]}$, $v_{\rm broad[1,2]}$, $R_{[1,2]}$, and $E(B-V)$), we used a Bayesian approach that includes a parallax prior from Gaia DR2, spectra from GALAH, and apparent magnitudes from APASS, Gaia DR2, 2MASS, and WISE. The derived stellar properties and their distributions show trends that are expected for a population of close binaries (a $<$ 10 AU) with mass ratios $0.5 \leq q \leq 1$. The derived metallicity of these binary stars is statistically lower than that of single dwarf stars from the same magnitude-limited sample.

Large pristine samples of red clump stars are highly sought after given that they are standard candles and give precise distances even at large distances. However, it is difficult to cleanly select red clumps stars because they can have the same T$_{\mathrm{eff}}$ and log $g$ as red giant branch stars. Recently, it was shown that the asteroseismic parameters, $\rm{\Delta}$P and $\rm{\Delta\nu}$, which are used to accurately select red clump stars, can be derived from spectra using the change in the surface carbon to nitrogen ratio ([C/N]) caused by mixing during the red giant branch. This change in [C/N] can also impact the spectral energy distribution. In this study, we predict the $\rm{\Delta}$P, $\rm{\Delta\nu}$, T$_{\mathrm{eff}}$ and log $g$ using 2MASS, AllWISE, \gaia, and Pan-STARRS data in order to select a clean sample of red clump stars. We achieve a contamination rate of $\sim$20\%, equivalent to what is achieved when selecting from T$_{\mathrm{eff}}$ and log $g$ derived from low resolution spectra. Finally, we present two red clump samples. One sample has a contamination rate of $\sim$ 20\% and $\sim$ 405,000 red clump stars. The other has a contamination of $\sim$ 33\% and $\sim$ 2.6 million red clump stars which includes $\sim$ 75,000 stars at distances $>$ 10 kpc. For |b|>30 degrees we find $\sim$ 15,000 stars with contamination rate of $\sim$ 9\%. The scientific potential of this catalog for studying the structure and formation history of the Galaxy is vast given that it includes millions of precise distances to stars in the inner bulge and distant halo where astrometric distances are imprecise.

The archaeological record of stars in the Milky Way opens a uniquely detailed window into the early formation and assembly of galaxies. Here we use 11,000 main-sequence turn-off stars with well-measured ages, [Fe/H], [$\alpha$/Fe], and orbits from the H3 Survey and Gaia to time the major events in the early Galaxy. Located beyond the Galactic plane, $1\lesssim |Z|/\rm kpc \lesssim4$, this sample contains three chemically distinct groups: a low metallicity population, and low-$\alpha$ and high-$\alpha$ groups at higher metallicity. The age and orbit distributions of these populations show that: 1) the high-$\alpha$ group, which includes both disk stars and the in-situ halo, has a star-formation history independent of eccentricity that abruptly truncated $8.3\pm0.1$ Gyr ago ($z\simeq1$); 2) the low metallicity population, which we identify as the accreted stellar halo, is on eccentric orbits and its star formation truncated $10.2.^{+0.2}_{-0.1}$ Gyr ago ($z\simeq2$); 3) the low-$\alpha$ population is primarily on low eccentricity orbits and the bulk of its stars formed less than 8 Gyr ago. These results suggest a scenario in which the Milky Way accreted a satellite galaxy at $z\approx2$ that merged with the early disk by $z\approx1$. This merger truncated star formation in the early high-$\alpha$ disk and perturbed a fraction of that disk onto halo-like orbits. The merger enabled the formation of a chemically distinct, low-$\alpha$ disk at $z\lesssim1$. The lack of any stars on halo-like orbits at younger ages indicates that this event was the last significant disturbance to the Milky Way disk.

We explore the fundamental relations governing the radial and vertical velocity dispersions of stars in the Milky Way, from combined studies of complementary surveys including GALAH, LAMOST, APOGEE, the NASA $Kepler$ and K2 missions, and $Gaia$ DR2. We find that different stellar samples, even though they target different tracer populations and employ a variety of age estimation techniques, follow the same set of fundamental relations. We provide the clearest evidence to date that, in addition to the well-known dependence on stellar age, the velocity dispersions of stars depend on orbital angular momentum $L_z$, metallicity and height above the plane $|z|$, and are well described by a multiplicatively separable functional form. The dispersions have a power-law dependence on age with exponents of 0.441$\pm 0.007$ and 0.251$\pm 0.006$ for $\sigma_z$ and $\sigma_R$ respectively, and the power law is valid even for the oldest stars. For the solar neighborhood stars, the apparent break in the power law for older stars, as seen in previous studies, is due to the anti-correlation of $L_z$ with age. The dispersions decrease with increasing $L_z$ until we reach the Sun's orbital angular momentum, after which $\sigma_z$ increases (implying flaring in the outer disc) while $\sigma_R$ flattens. The dispersions increase with decreasing metallicity, suggesting that the dispersions increase with birth radius. The dispersions also increase linearly with $|z|$. The same set of relations that work in the solar neighborhood also work for stars between $3<R/{\rm kpc}<20$. Finally, the high-[$\alpha$/Fe] stars follow the same relations as the low-[$\alpha$/Fe] stars.

We present the serendipitous discovery of a low luminosity nova occurring in a symbiotic binary star system in the Milky Way. We lay out the extensive archival data alongside new follow-up observations related to the stellar object V$^*$ CN Cha in the constellation of Chamaeleon. The object had long period ($\sim\! 250\,$day), high amplitude ($\sim\! 3\,$mag) optical variability in its recent past, preceding an increase in optical brightness by $\sim\! 8\,$magnitudes and a persistence at this luminosity for about 3 years, followed by a period of $\sim\! 1.4\,{\rm mag}\,{\rm yr}^{-1}$ dimming. The object's current optical luminosity seems to be dominated by H$\alpha$ emission, which also exhibits blue-shifted absorption (a P-Cygni-like profile). After consideration of a number of theories to explain these myriad observations, we determine that V$^*$ CN Cha is most likely a symbiotic (an evolved star-white dwarf binary) system which has undergone a long-duration, low luminosity, nova. Interpreted in this way, the outburst in V$^*$ CN Cha is among the lowest luminosity novae ever observed.

A star in the Milky Way's disk can now be at a Galactocentric radius quite distant from its birth radius for two reasons: either its orbit has become eccentric through radial heating, which increases its radial action $J_R$ (`blurring'); or merely its angular momentum $L_z$ has changed and thereby its guiding radius (`churning'). We know that radial orbit migration is strong in the Galactic low-$\alpha$ disk and set out to quantify the relative importance of these two effects, by devising and applying a parameterized model for the distribution $p(L_z, J_R, \tau, \mathrm[Fe/H])$ in the stellar disk. This model describes the orbit evolution for stars of age $\tau$ and metallicity [Fe/H], presuming coeval stars were initially born on (near-)circular orbits, and with a unique [Fe/H] at a given birth angular momentum and age. We fit this model to APOGEE red clump stars, accounting for the complex selection function of the survey. The best fit model implies changes of angular momentum of $\sqrt{\langle \Delta L_z \rangle^2} \approx 619\, \mathrm{kpc~km/s~}(\tau/\mathrm{6~Gyr})^{0.5}$, and changes of radial action as $\sqrt{\langle \Delta J_R \rangle^2} \approx 63\, \mathrm{kpc~km/s~} (\tau/\mathrm{6~Gyr})^{0.6}$ at 8 kpc. This suggests that the secular orbit evolution of the disk is dominated by diffusion in angular momentum, with radial heating being an order of magnitude lower.

One of the high-level goals of Galactic archaeology is chemical tagging of stars across the Milky Way to piece together its assembly history. For this to work, stars born together must be uniquely chemically homogeneous. Wide binary systems are an important laboratory to test this underlying assumption. Here we present the detailed chemical abundance patterns of 50 stars across 25 wide binary systems comprised of main-sequence stars of similar spectral type identified in Gaia DR2 with the aim of quantifying their level of chemical homogeneity. Using high-resolution spectra obtained with McDonald Observatory, we derive stellar atmospheric parameters and precise detailed chemical abundances for light/odd-Z (Li, C, Na, Al, Sc, V, Cu), $\alpha$ (Mg, Si, Ca), Fe-peak (Ti, Cr, Mn, Fe, Co, Ni, Zn), and neutron capture (Sr, Y, Zr, Ba, La, Nd, Eu) elements. Results indicate that 80% (20 pairs) of the systems are homogeneous in [Fe/H] at levels below 0.02 dex. These systems are also chemically homogeneous in all elemental abundances studied, with offsets and dispersions consistent with measurement uncertainties. We also find that wide binary systems are far more chemically homogeneous than random pairings of field stars of similar spectral type. These results indicate that wide binary systems tend to be chemically homogeneous but in some cases they can differ in their detailed elemental abundances at a level of [X/H] ~ 0.10 dex, overall implying chemical tagging in broad strokes can work.

We present isochrone ages and initial bulk metallicities ($\rm [Fe/H]_{bulk}$, by accounting for diffusion) of 163,722 stars from the GALAH Data Release 2, mainly composed of main sequence turn-off stars and subgiants ($\rm 7000 K>T_{eff}>4000 K$ and $\rm log g>3$ dex). The local age-metallicity relationship (AMR) is nearly flat but with significant scatter at all ages; the scatter is even higher when considering the observed surface abundances. After correcting for selection effects, the AMR appear to have intrinsic structures indicative of two star formation events, which we speculate are connected to the thin and thick disks in the solar neighborhood. We also present abundance ratio trends for 16 elements as a function of age, across different $\rm [Fe/H]_{bulk}$ bins. In general, we find the trends in terms of [X/Fe] vs age from our far larger sample to be compatible with studies based on small ($\sim$ 100 stars) samples of solar twins but we now extend it to both sub- and super-solar metallicities. The $\alpha$-elements show differing behaviour: the hydrostatic $\alpha$-elements O and Mg show a steady decline with time for all metallicities while the explosive $\alpha$-elements Si, Ca and Ti are nearly constant during the thin disk epoch (ages $\lessapprox $ 12 Gyr). The s-process elements Y and Ba show increasing [X/Fe] with time while the r-process element Eu have the opposite trend, thus favouring a primary production from sources with a short time-delay such as core-collapse supernovae over long-delay events such as neutron star mergers.

Using kinematics from Gaia and the large elemental abundance space of the second data release of the GALAH survey, we identify two new members of the Fimbulthul stellar stream, and chemically tag them to massive, multi-metallic globular cluster $\omega$ Centauri. Recent analysis of the second data release of Gaia had revealed the Fimbulthul stellar stream in the halo of the Milky Way. It had been proposed that the stream is associated with the $\omega$ Centauri, but this proposition relied exclusively upon the kinematics and metallicities of the stars to make the association. In this work, we find our two new members of the stream to be metal-poor stars that are enhanced in sodium and aluminium, typical of second population globular cluster stars, but not otherwise seen in field stars. Furthermore, the stars share the s-process abundance pattern seen in $\omega$ Centauri, which is rare in field stars. Apart from one star within 1.5 deg of $\omega$ Centauri, we find no other stars observed by GALAH spatially near $\omega$ Centauri or the Fimbulthul stream that could be kinematically and chemically linked to the cluster. Chemically tagging stars in the Fimbulthul stream to $\omega$ Centauri confirms the earlier work, and further links this tidal feature in the Milky Way halo to $\omega$ Centauri.

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.

Gaia DR2 has revealed new small-scale and large-scale patterns in the phase-space distribution of stars in the Milky Way. In cylindrical Galactic coordinates $(R,\phi,z)$, ridge-like structures can be seen in the \vphiR{} plane and asymmetric arch-like structures in the \vphivR{} plane. We show that the ridges are also clearly present when the third dimension of the \vphiR{} plane is represented by $\langle z \rangle$, $\langle V_z \rangle$, $\langle V_R \rangle$, $\langle$[Fe/H]$\rangle$ and $\langle[\alpha/{\rm Fe}]\rangle$. The maps suggest that stars along the ridges lie preferentially close to the Galactic midplane ($|z|<0.2$ kpc), and have metallicity and $\alpha$ elemental abundance similar to that of the Sun. We show that phase mixing of disrupting spiral arms can generate both the ridges and the arches. It also generates discrete groupings in orbital energy $-$ the ridges and arches are simply surfaces of constant energy. We identify 8 distinct ridges in the \gaia{} data: six of them have constant energy while two have constant angular momentum. Given that the signature is strongest for stars close to the plane, the presence of ridges in $\langle z \rangle$ and $\langle V_z \rangle$ suggests a coupling between planar and vertical directions. We demonstrate, using N-body simulations that such coupling can be generated both in isolated discs and in discs perturbed by an orbiting satellite like the Sagittarius dwarf galaxy.

We present the determination of stellar parameters and individual elemental abundances for 6 million stars from $\sim$8 million low-resolution ($R\sim1800$) spectra from LAMOST DR5. This is based on a modeling approach that we dub $The$ $Data$--$Driven$ $Payne$ ($DD$--$Payne$), which inherits essential ingredients from both {\it The Payne} \citep{Ting2019} and $The$ $Cannon$ \citep{Ness2015}. It is a data-driven model that incorporates constraints from theoretical spectral models to ensure the derived abundance estimates are physically sensible. Stars in LAMOST DR5 that are in common with either GALAH DR2 or APOGEE DR14 are used to train a model that delivers stellar parameters ($T_{\rm eff}$, $\log g$, $V_{\rm mic}$) and abundances for 16 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, Co, Ni, Cu, and Ba) when applied to LAMOST spectra. Cross-validation and repeat observations suggest that, for ${\rm S/N}_{\rm pix}\ge 50$, the typical internal abundance precision is 0.03--0.1\,dex for the majority of these elements, with 0.2--0.3\,dex for Cu and Ba, and the internal precision of $T_{\rm eff}$ and $\log g$ is better than 30\,K and 0.07\,dex, respectively. Abundance systematics at the $\sim$0.1\,dex level are present in these estimates, but are inherited from the high-resolution surveys' training labels. For some elements, GALAH provides more robust training labels, for others, APOGEE. We provide flags to guide the quality of the label determination and to identify binary/multiple stars in LAMOST DR5. The abundance catalogs are publicly accessible via \href{url}{http://dr5.lamost.org/doc/vac}.

A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR.

High-multiplex and deep spectroscopic follow-up of upcoming panoramic deep-imaging surveys like LSST, Euclid, and WFIRST is a widely recognized and increasingly urgent necessity. No current or planned facility at a U.S. observatory meets the sensitivity, multiplex, and rapid-response time needed to exploit these future datasets. FOBOS, the Fiber-Optic Broadband Optical Spectrograph, is a near-term fiber-based facility that addresses these spectroscopic needs by optimizing depth over area and exploiting the aperture advantage of the existing 10m Keck II Telescope. The result is an instrument with a uniquely blue-sensitive wavelength range (0.31-1.0 um) at R~3500, high-multiplex (1800 fibers), and a factor 1.7 greater survey speed and order-of-magnitude greater sampling density than Subaru's Prime Focus Spectrograph (PFS). In the era of panoramic deep imaging, FOBOS will excel at building the deep, spectroscopic reference data sets needed to interpret vast imaging data. At the same time, its flexible focal plane, including a mode with 25 deployable integral-field units (IFUs) across a 20 arcmin diameter field, enables an expansive range of scientific investigations. Its key programmatic areas include (1) nested stellar-parameter training sets that enable studies of the Milky Way and M31 halo sub-structure, as well as local group dwarf galaxies, (2) a comprehensive picture of galaxy formation thanks to detailed mapping of the baryonic environment at z~2 and statistical linking of evolving populations to the present day, and (3) dramatic enhancements in cosmological constraints via precise photometric redshifts and determined redshift distributions. In combination with Keck I instrumentation, FOBOS also provides instant access to medium-resolution spectroscopy for transient sources with full coverage from the UV to the K-band.

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.

The problem of the origin of the elements is a fundamental one in astronomy and one that has many open questions. Prominent examples include (1) the nature of Type Ia supernovae and the timescale of their contributions; (2) the observational identification of elements such as titanium and potassium with the $\alpha$-elements in conflict with core-collapse supernova predictions; (3) the number and relative importance of r-process sites; (4) the origin of carbon and nitrogen and the influence of mixing and mass loss in winds; and (5) the origin of the intermediate elements, such as Cu, Ge, As, and Se, that bridge the region between charged-particle and neutron-capture reactions. The next decade will bring to maturity many of the new tools that have recently made their mark, such as large-scale chemical cartography of the Milky Way and its satellites, the addition of astrometric and asteroseismic information, the detection and characterization of gravitational wave events, 3-D simulations of convection and model atmospheres, and improved laboratory measurements for transition probabilities and nuclear masses. All of these areas are key for continued improvement, and such improvement will benefit many areas of astrophysics.

Aims We study the effects of non-local thermodynamic equilibrium (NLTE) on the determination of stellar parameters and abundances of Fe, Mg, and Ti from the medium-resolution spectra of FGK stars.} Methods We extend the Payne fitting approach to draw on NLTE and LTE spectral models. These are used to analyse the spectra of the Gaia-ESO benchmark stars and the spectra of 742 stars in 13 open and globular clusters in the Milky Way: NGC 3532, NGC 5927, NGC 2243, NGC 104, NGC 1851, NGC 2808, NGC 362, M 2, MGC 6752, NGC 1904, NGC 4833, NGC 4372 and M 15. Results Our approach accurately recovers effective temperatures, surface gravities, and abundances of the benchmark stars and clusters members. The differences between NLTE and LTE are significant in the metal-poor regime, [Fe/H] $\lesssim -1$. The NLTE [Fe/H] values are systematically higher, whereas the average NLTE [Mg/Fe] abundance ratios are $\sim 0.15$ dex lower, compared to LTE. Our LTE measurements of metallicities and abundances of stars in Galactic clusters are in a good agreement with the literature. Yet, for most clusters, our study yields the first estimates of NLTE abundances of Fe, Mg and Ti. Conclusions All clusters investigated in this work are homogeneous in Fe and Ti, with the intra-cluster abundance variations of less than 0.04 dex. NGC 2808, NGC 4833, M 2 and M 15 show significant dispersions in [Mg/Fe]. Contrary to common assumptions, the NLTE analysis changes the mean abundance ratios in the clusters, but it does not influence the intra-cluster abundance dispersions.

The latest Gaia data release enables us to accurately identify stars that are more luminous than would be expected on the basis of their spectral type and distance. During an investigation of the 329 best Solar twin candidates uncovered among the spectra acquired by the GALAH survey, we identified 64 such over-luminous stars. In order to investigate their exact composition, we developed a data-driven methodology that can generate a synthetic photometric signature and spectrum of a single star. By combining multiple such synthetic stars into an unresolved binary or triple system and comparing the results to the actual photometric and spectroscopic observations, we uncovered 6 definitive triple stellar system candidates and an additional 14 potential candidates whose combined spectrum mimics the Solar spectrum. Considering the volume correction factor for a magnitude limited survey, the fraction of probable unresolved triple stars with long orbital periods is ~2 %. Possible orbital configurations of the candidates were investigated using the selection and observational limits. To validate the discovered multiplicity fraction, the same procedure was used to evaluate the multiplicity fraction of other stellar types.

We present The Payne, a general method for the precise and simultaneous determination of numerous stellar labels from observed spectra, based on fitting physical spectral models. The Payne combines a number of important methodological aspects: it exploits the information from much of the available spectral range; it fits all labels (stellar parameters and element abundances) simultaneously; it uses spectral models, where the atmosphere structure and the radiative transport are consistently calculated to reflect the stellar labels. At its core The Payne has an approach to accurate and precise interpolation and prediction of the spectrum in high-dimensional label-space, which is flexible and robust, yet based on only a moderate number of ab initio models (O(1000) for 25 labels). With a simple neural-net-like functional form and a suitable choice of training labels, this interpolation yields a spectral flux prediction good to $10^{-3}$ rms across a wide range of $T_{\rm eff}$ and log g (including dwarfs and giants). We illustrate the power of this approach by applying it to the APOGEE DR14 data set, drawing on Kurucz models with recently improved line lists: without recalibration, we obtain physically sensible stellar parameters as well as 15 element abundances that appear to be more precise than the published APOGEE DR14 values. In short, The Payne is an approach that for the first time combines all these key ingredients, necessary for progress towards optimal modelling of survey spectra; and it leads to both precise and accurate estimates of stellar labels, based on physical models and without re-calibration. Both the codes and catalog are made publicly available online.

It has long been known that the vertical motions of Galactic disk stars increase with stellar age, commonly interpreted as vertical heating through orbit scattering. Here we map the vertical actions of disk stars as a function of age ($\tau\le$ 8 Gyrs) and across a large range of Galactocentric radii, $\overline{R}_{\rm GC}$, drawing on APOGEE and Gaia data. We fit $\widehat{J_z}(\overline{R}_{\rm GC},\tau )$ as a combination of the vertical action at birth, $\widehat{J_{z,0}}$, and subsequent heating $\widehat{\Delta J_z}_{\, {\rm 1Gyr}}(\overline{R}_{\rm GC})$ that scales as $\tau^{\gamma(\overline{R}_{\rm GC})}$. The inferred birth temperature, $\widehat{J_{z,0}}(\overline{R}_{\rm GC})$ is 1 kpc km/s for 3 kpc < $\overline{R}_{\rm GC}$ < 10 kpc, consistent with the ISM velocity dispersion; but it rapidly rises outward, to 8 kpc km/s for $\overline{R}_{\rm GC}$ = 14 kpc, likely reflecting the stars' birth in a warped or flared gas disk. We find the heating rate $\widehat{\Delta J_z}_{\, {\rm 1Gyr}}$ to be modest and nearly constant across all radii, 1.6 kpc km/s Gyr$^{-1}$. The stellar age dependence $\gamma$ gently grows with Galactocentric radius, from $\gamma \simeq$ 1 for $\overline{R}_{\rm GC}\lesssim R_\odot$ to $\gamma \simeq$ 1.3 at $\overline{R}_{\rm GC}$ =14 kpc. The observed $J_z - \tau$ relation at all radii is considerably steeper ($\gamma\gtrsim$ 1) than the time dependence theoretically expected from orbit scattering, $J_z\propto t^{0.5}$. We illustrate how this conundrum can be resolved if we also account for the fact that at earlier epochs the scatterers were more common, and the restoring force from the stellar disk surface mass density was low. Our analysis may re-instate gradual orbital scattering as a plausible and viable mechanism to explain the age-dependent vertical motions of disk stars.

Asteroseismology is a promising tool to study Galactic structure and evolution because it can probe the ages of stars. Earlier attempts comparing seismic data from the {\it Kepler} satellite with predictions from Galaxy models found that the models predicted more low-mass stars compared to the observed distribution of masses. It was unclear if the mismatch was due to inaccuracies in the Galactic models, or the unknown aspects of the selection function of the stars. Using new data from the K2 mission, which has a well-defined selection function, we find that an old metal-poor thick disc, as used in previous Galactic models, is incompatible with the asteroseismic information. We show that spectroscopic measurements of [Fe/H] and [$\alpha$/Fe] elemental abundances from the GALAH survey indicate a mean metallicity of $\log (Z/Z_{\odot})=-0.16$ for the thick disc. Here $Z$ is the effective solar-scaled metallicity, which is a function of [Fe/H] and [$\alpha$/Fe]. With the revised disc metallicities, for the first time, the theoretically predicted distribution of seismic masses show excellent agreement with the observed distribution of masses. This provides an indirect verification of the asteroseismic mass scaling relation is good to within five percent. Using an importance-sampling framework that takes the selection function into account, we fit a population synthesis model of the Galaxy to the observed seismic and spectroscopic data. Assuming the asteroseismic scaling relations are correct, we estimate the mean age of the thick disc to be about 10 Gyr, in agreement with the traditional idea of an old $\alpha$-enhanced thick disc.

We investigate aluminum abundance variations in the stellar populations of globular clusters using both literature measurements of sodium and aluminum and APOGEE measurements of nitrogen and aluminum abundances. For the latter, we show that the Payne is the most suitable of the five available abundance pipelines for our purposes. Our combined sample of 42 globular clusters spans approximately 2 dex in [Fe/H] and 1.5 dex in $\log{M_{GC}/M_{\odot}}$. We find no fewer than five globular clusters with significant internal variations in nitrogen and/or sodium with little-to-no corresponding variation in aluminum, and that the minimum present-day cluster mass for aluminum enrichment in metal-rich systems is $\log{M_{GC}/M_{\odot}} \approx 4.50 + 2.17(\rm{[Fe/H]}+1.30)$. We demonstrate that the slopes of the [Al/Fe] vs [Na/Fe] and [Al/Fe] vs [N/Fe] relations for stars without field-like abundances are approximately log-linearly dependent on both the metallicity and the stellar mass of the globular clusters. In contrast, the relationship between [Na/Fe] and [N/Fe] shows no evidence of such dependencies. This suggests that there were (at least) two classes of non-supernovae chemical polluters that were common in the early universe, and that their relative contributions within globular clusters somehow scaled with the metallicity and mass of globular clusters. The first of these classes is predominantly responsible for the CNO and NeNa abundance variations, and likewise the second for the MgAl abundance variations. 47 Tuc and M4 are particularly striking examples of this dichotomy. As an auxiliary finding, we argue that abundance variations among Terzan 5 stars are consistent with it being a normal globular cluster.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

It is challenging to reliably identify stars that were born together outside of actively star-forming regions and bound stellar systems. However, co-natal stars should be present throughout the Galaxy, and their demographics can shed light on the clustered nature of star formation and the dynamical state of the disk. In previous work we presented a set of simulations of the Galactic disk that followed the clustered formation and dynamical evolution of 4 billion individual stars over the last 5 Gyr. The simulations predict that a high fraction of co-moving stars with physical and 3D velocity separation of $\Delta r < 20$ pc and $\Delta v < 1.5$ km s$^{-1}$ are co-natal. In this \textit{Letter}, we use \textit{Gaia} DR2 and LAMOST DR4 data to identify and study co-moving pairs. We find that the distribution of relative velocities and separations of pairs in the data is in good agreement with the predictions from the simulation. We identify 111 co-moving pairs in the Solar neighborhood with reliable astrometric and spectroscopic measurements. These pairs show a strong preference for having similar metallicities when compared to random field pairs. We therefore conclude that these pairs were very likely born together. The simulations predict that co-natal pairs originate preferentially from high-mass and relatively young ($< 1$ Gyr) star clusters. \textit{Gaia} will eventually deliver well-determined metallicities for the brightest stars, enabling the identification of thousands of co-natal pairs due to disrupting star clusters in the solar neighborhood.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

The distribution of gas-phase abundances in galaxy disks encodes the history of nucleosynthesis and transport through the interstellar medium (ISM) over cosmic time. Multi-object and high resolution integral-field spectroscopy have started to measure these distributions across hundreds of HII regions individually resolved at $\lesssim 100$ pc scales in a handful of objects, but in the coming decade these studies will expand to larger samples of galaxies. This will allow us to understand the role of feedback and turbulence in driving the mixing and diffusion of metals in the ISM, and statistically assess the role of galaxy environment and disk dynamics in modifying how mixing proceeds. Detailed searches for over- and under-enriched regions can address to what extent star formation is triggered by previous generations of star formation and by pristine and recycled gas flows. Local galaxies, for which these detailed measurements will be possible within the next decade, will inform the interpretation of integrated measurements at high-z, where very different dynamical gas-rich environments are found in early disk galaxies. Currently, progress in the field is severely hampered by the 0.2-0.3 dex level systematic uncertainties plaguing nebular abundance diagnostics. Improving our detailed understanding of ionized nebulae at $<$20 pc scales will help us find a solution to this problem, which will prove key to the study of metal enrichment and mixing across the galaxy population in the next decade.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

The clustered nature of star formation should produce a high degree of structure in the combined phase and chemical space in the Galactic disk. To date, observed structure of this kind has been mostly limited to bound clusters and moving groups. In this paper we present a new dynamical model of the Galactic disk that takes into account the clustered nature of star formation. This model predicts that the combined phase and chemical space is rich in substructure, and that this structure is sensitive to both the precise nature of clustered star formation and the large-scale properties of the Galaxy. The model self-consistently evolves 4 billion stars over the last 5 Gyr in a realistic potential that includes an axisymmetric component, a bar, spiral arms, and giant molecular clouds (GMCs). All stars are born in clusters with an observationally-motivated range of initial conditions. As direct \textit{N}-body calculations for billions of stars is computationally infeasible, we have developed a method of initializing star cluster particles to mimic the effects of direct \textit{N}-body effects, while the actual orbit integrations are treated as test particles within the analytic potential. We demonstrate that the combination of chemical and phase space information is much more effective at identifying truly co-natal populations than either chemical or phase space alone. Furthermore, we show that co-moving pairs of stars are very likely to be co-natal if their velocity separation is $< 2$ km s$^{-1}$ and their metallicity separation is $< 0.05$ dex. The results presented here bode well for harnessing the synergies between \textit{Gaia} and spectroscopic surveys to reveal the assembly history of the Galactic disk.

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 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.

Swan bands - characteristic molecular absorption features of the C$_2$ molecule - are a spectroscopic signature of carbon-enhanced stars. They can also be used to identify carbon-enhanced metal-poor (CEMP) stars. The GALAH (GALactic Archaeology with Hermes) is a magnitude-limited survey of stars producing high-resolution, high signal-to-noise spectra. We used 627,708 GALAH spectra to search for carbon-enhanced stars with a supervised and unsupervised classification algorithm, relying on the imprint of the Swan bands. We identified 918 carbon-enhanced stars, including 12 already described in the literature. An unbiased selection function of the GALAH survey allows us to perform a population study of carbon-enhanced stars. Most of them are giants, out of which we find 28 CEMP candidates. A large fraction of our carbon-enhanced stars with repeated observations show variation in radial velocity, hinting that there is a large fraction of variables among them. 32 of the detected stars also show strong Lithium enhancement in their spectra.

If the Galaxy is axisymmetric and in dynamical equilibrium, we expect negligible fluctuations in the residual line-of-sight velocity field. Recent results using the \apg{} survey find significant fluctuations in velocity for stars in the midplane ($|z|<$0.25 kpc) out to 5 kpc, suggesting that the dynamical influence of non-axisymmetric features i.e., the Milky Way's bar, spiral arms and merger events extends out to the Solar neighborhood. Their measured power spectrum has a characteristic amplitude of 11 \kms{} on a scale of 2.5 kpc. The existence of such large-scale streaming motions has important implications for determining the Sun's motion about the Galactic Centre. Using Red Clump stars from \glh{} and \apg{}, we map the line-of-sight velocities around the Sun (d$<$5 kpc), and $|z|<$1.25 kpc from the midplane. By subtracting a smooth axisymmetric model for the velocity field, we study the residual fluctuations and compare our findings with mock survey generated by \glx{}. We find negligible large-scale fluctuations away from the plane. In the mid-plane, we reproduce the earlier \apg{} power spectrum but with 20\% smaller amplitude (9.3 \kms{}) after taking into account a few systematics (e.g., volume completeness). Using a flexible axisymmetric model the power-amplitude is further reduced to 6.3 \kms{}. Additionally, our simulations show that, in the plane, distances are underestimated for high-mass Red Clump stars which can lead to spurious power-amplitude of about 5.2 \kms{}. Taking this into account, we estimate the amplitude of real fluctuations to be $<$4.6 \kms{}, about a factor of three less than the \apg{} result.

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.

Open cluster members are coeval and share the same initial bulk chemical compositions. Consequently, differences in surface abundances between members of a cluster that are at different evolutionary stages can be used to study the effects of mixing and internal chemical processing. We carry out an abundance analysis of seven elements (Li, O, Na, Mg, Al, Si, Fe) in 66 stars belonging to the open cluster M67, based on high resolution GALAH spectra, 1D MARCS model atmospheres, and, for the first time for a large sample of M67 stars, non-local thermodynamic equilibrium (non-LTE) radiative transfer. From the non-LTE analysis, we find a typical star-to-star scatter in the abundance ratios of around 0.05 dex; this scatter is slightly but systematically larger when LTE is assumed instead. We find trends in the abundance ratios with effective temperature, indicating systematic differences in the surface abundances between turn-off and giant stars; these trends are more pronounced when LTE is assumed. However, in the non-LTE analysis, most of the element trends have been flattened. Two species are exceptions to this behaviour, namely Al and Si, which both clearly display remaining trends in the non-LTE analysis. We comment on the possible origin of these trends, by comparing them with recent stellar models that include atomic diffusion.

We use the framework developed as part of the MESA Isochrones and Stellar Tracks (MIST) project to assess the utility of several types of observables in jointly measuring the age and 1D stellar model parameters in star clusters. We begin with a pedagogical overview summarizing the effects of stellar model parameters, such as the helium abundance, mass-loss efficiency, and the mixing length parameter, on observational diagnostics including the color-magnitude diagram, mass-radius relation, and surface abundances, amongst others. We find that these parameters and the stellar age influence observables in qualitatively distinctive, degeneracy-breaking ways. To assess the current state of affairs, we use the recent Gaia Data Release 2 (DR2) along with data from the literature to investigate three well-studied old open clusters---NGC6819, M67, NGC6791---as case studies. Although there is no obvious tension between the existing observations and the MIST models for NGC6819, there are interesting discrepancies in the cases of M67 and NGC6791. At this time, parallax zero point uncertainties in Gaia DR2 remain one of the limiting factors in the analysis of these clusters. With a combination of exquisite photometry, parallax distances, and cluster memberships from Gaia at the end of its mission, we anticipate precise and accurate ages for these and other star clusters in the Galaxy.

Sparse open clusters can be found at high galactic latitudes where loosely populated clusters are more easily detected against the lower stellar background. Because most star formation takes place in the thin disk, the observed population of clusters far from the Galactic plane is hard to explain. We combined spectral parameters from the GALAH survey with the Gaia DR2 catalogue to study the dynamics and chemistry of five old sparse high-latitude clusters in more detail. We find that four of them (NGC 1252, NGC 6994, NGC 7772, NGC 7826) - originally classified in 1888 - are not clusters but are instead chance projections on the sky. Member stars quoted in the literature for these four clusters are unrelated in our multi-dimensional physical parameter space; the published cluster properties are therefore irrelevant. We confirm the existence of morphologically similar NGC 1901 for which we provide a probabilistic membership analysis. An overdensity in three spatial dimensions proves to be enough to reliably detect sparse clusters, but the whole 6-dimensional space must be used to identify members with high confidence, as demonstrated in the case of NGC 1901.

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 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.

Oxygen is the most abundant "metal" element in stars and in the cosmos. But determining oxygen abundances in stars has proven challenging, because of the shortage of detectable atomic oxygen lines in their optical spectra as well as observational and theoretical complications with these lines (e.g., blends, 3D, non-LTE). Nonetheless, Ting et al. (2017) were recently able to demonstrate that oxygen abundances can be determined from low-resolution (R$\simeq$2000) optical spectra. Here we investigate the physical processes that enable such a measurement for cool stars, such as K-giants. We show that the strongest spectral diagnostics of oxygen come from the CNO atomic-molecular network, but are manifested in spectral features that do not involve oxygen. In the outer atmosphere layers most of the carbon is locked up in CO, and changes to the oxygen abundance directly affect the abundances of all other carbon-bearing molecules, thereby changing the strength of CH, CN, and C$_2$ features across the optical spectrum. In deeper atmosphere layers most of the carbon is in atomic form, and any change in the oxygen abundance has little effect on the other carbon-bearing molecules. The key physical effect enabling such oxygen abundance measurements is that spectral features in the optical arise from both the CO-dominant and the atomic carbon-dominant regions, providing non-degenerate constraints on both C and O. Beyond the case at hand, the results show that physically sound abundances measurements need not be limited to those elements that have observable lines themselves.

The effective temperature (Teff) distribution of stellar evolution models along the red giant branch (RGB) is sensitive to a number of parameters including the overall metallicity, elemental abundance patterns, the efficiency of convection, and the treatment of the surface boundary condition. Recently there has been interest in using observational estimates of the RGB Teff to place constraints on the mixing length parameter, a_MLT, and possible variation with metallicity. Here we use 1D MESA stellar evolution models to explore the sensitivity of the RGB Teff to the treatment of the surface boundary condition. We find that different surface boundary conditions can lead to +/- 100 K metallicity-dependent offsets on the RGB relative to one another in spite of the fact that all models can reproduce the properties of the Sun. Moreover, for a given atmosphere T-tau relation, we find that the RGB Teff is also sensitive to the optical depth at which the surface boundary condition is applied in the stellar model. Nearly all models adopt the photosphere as the location of the surface boundary condition but this choice is somewhat arbitrary. We compare our models to stellar parameters derived from the APOGEE-Kepler sample of first ascent red giants and find that systematic uncertainties in the models due to treatment of the surface boundary condition place a limit of ~100 K below which it is not possible to make firm conclusions regarding the fidelity of the current generation of stellar models.

We report the serendipitous observations of 571 luminous supergiants in the Magellanic Clouds by the spectroscopic GALAH and TESS-HERMES surveys: 434 stars in the Large Magellanic Cloud and 137 in the Small Magellanic Cloud. We also find one star that appears associated with structured star formation in the Magellanic Bridge. Both of these surveys are aimed at the local volume of the Galaxy but have simple, magnitude-limited selection functions that mean they include some observations of luminous extra-Galactic stars. The surveys determine stellar parameter and abundances using The Cannon, a data-driven generative modelling approach. In this work, we explore the results from The Cannon when it is fed the spectra of these intrinsically luminous supergiants in the Magellanic Clouds, which are well outside the normal bounds of The Cannon's training set. We find that, although the parameters are astrophysically incorrect, the $v\sin i$ and the abundances of lithium, barium, and magnesium are excellent discriminants of these stars. It shows that in the future, with an expanded training set, it should be possible to determine accurate values for these types of stars.

Core helium-burning red clump (RC) stars are excellent standard candles in the Milky Way. These stars may have more precise distance estimates from spectrophotometry than from Gaia parallaxes beyond 3 kpc. However, RC stars have $T_{\rm eff}$ and log g very similar to some red giant branch (RGB) stars. Especially for low-resolution spectroscopic studies where $T_{\rm eff}$, log g, and [Fe/H] can only be estimated with limited precision, separating RC stars from RGB through established method can incur ~20% contamination. Recently, Hawkins et al. (2018) demonstrated that the additional information in single-epoch spectra, such as the C/N ratio, can be exploited to cleanly differentiate RC and RGB stars. In this second paper of the series, we establish a data-driven mapping from spectral flux space to independently determined asteroseismic parameters, the frequency and the period spacing. From this, we identify 210,371 RC stars from the publicly available LAMOST DR3 and APOGEE DR14 data, with ~9% of contamination. We provide an RC sample of 92,249 stars with a contamination of only ~3%, by restricting the combined analysis to LAMOST stars with ${\rm S/N}_{\rm pix}$ > 75. This demonstrates that high-S/N, low-resolution spectra covering a broad wavelength range can identify RC samples at least as pristine as their high-resolution counterparts. As coming and ongoing surveys such as TESS, DESI, and LAMOST will continue to improve the overlapping training spectroscopic-asteroseismic sample, the method presented in this study provides an efficient and straightforward way to derive a vast yet pristine RC stars to reveal the 3D structure of the Milky Way.

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.

The primary goal of Galactic archaeology is to learn about the origin of the Milky Way from the detailed chemistry and kinematics of millions of stars. Wide-field multi-fibre spectrographs are increasingly used to obtain spectral information for huge samples of stars. Some surveys (e.g. GALAH) are attempting to measure up to 30 separate elements per star. Stellar abundance spectroscopy is a subtle art that requires a very high degree of spectral uniformity across each of the fibres. However wide-field spectrographs are notoriously non-uniform due to the fast output optics necessary to image many fibre outputs onto the detector. We show that precise spectroscopy is possible with such instruments across all fibres by employing a photonic comb -- a device that produces uniformly spaced spots of light on the CCD to precisely map complex aberrations. Aberrations are parametrized by a set of orthogonal moments with $\sim100$ independent parameters. We then reproduce the observed image by convolving high resolution spectral templates with measured aberrations as opposed to extracting the spectra from the observed image. Such a forward modeling approach also trivializes some spectroscopic reduction problems like fibre cross-talk, and reliably extracts spectra with a resolution $\sim2.3$ times above the nominal resolution of the instrument. Our rigorous treatment of optical aberrations also encourages a less conservative spectrograph design in the future.

Using GALAH survey data of nearby stars, we look at how structure in the planar (u,v) velocity distribution depends on metallicity and on viewing direction within the Galaxy. In nearby stars, with distance d < 1 kpc, the Hercules stream is most strongly seen in higher metallicity stars [Fe/H] > 0.2. The Hercules stream peak v value depends on viewed galactic longitude, which we interpret as due to the gap between the stellar stream and more circular orbits being associated with a specific angular momentum value of about 1640 km/s kpc. The association of the gap with a particular angular momentum value supports a bar resonant model for the Hercules stream. Moving groups previously identified in Hipparcos observations are easiest to see in stars nearer than 250 pc, and their visibility and peak velocities in the velocity distributions depends on both viewing direction (galactic longitude and hemisphere) and metallicity. We infer that there is fine structure in local velocity distributions that varies over distances of a few hundred pc in the Galaxy.

We present a Bayesian implementation of isochrone fitting in deriving stellar ages and masses, incorporating absolute K magnitude ($\rm M_K$) derived from 2MASS photometry and Gaia DR1 parallax and differentiation between initial bulk metallicity and present day surface metallicity, with allowance for incorporating further constraints (e.g., asteroseismology) when available. As a test, we re-computed stellar ages and masses of $\sim 4000$ stars in the solar neighbourhood from six well-studied literature samples using both Hipparcos and TGAS parallaxes. Our ages are found to be compatible with literature values but with reduced uncertainties in general. The inclusion of parallax-based $\rm M_K$ serves as an additional constraint on the derived quantities, especially when systematic errors in stellar parameters are underestimated. We reconstructed the age-metallicity relationship in the solar neighbourhood by re-analysing the Geneva-Copenhagen Survey with the inclusion of TGAS-parallaxes and initial bulk metallicity sampling. We found a flat trend for disk stars with ages $<$11\,Gyr but with smaller scatter at all ages compared to literature.

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.

We develop a data-driven spectral model for identifying and characterizing spatially unresolved multiple-star systems and apply it to APOGEE DR13 spectra of main-sequence stars. Binaries and triples are identified as targets whose spectra can be significantly better fit by a superposition of two or three model spectra, drawn from the same isochrone, than any single-star model. From an initial sample of $\sim$20,000 main-sequence targets, we identify $\sim$2,500 binaries in which both the primary and secondary star contribute detectably to the spectrum, simultaneously fitting for the velocities and stellar parameters of both components. We additionally identify and fit $\sim$200 triple systems, as well as $\sim$700 velocity-variable systems in which the secondary does not contribute detectably to the spectrum. Our model simplifies the process of simultaneously fitting single- or multi-epoch spectra with composite models and does not depend on a velocity offset between the two components of a binary, making it sensitive to traditionally undetectable systems with periods of hundreds or thousands of years. In agreement with conventional expectations, almost all the spectrally-identified binaries with measured parallaxes fall above the main sequence in the color-magnitude diagram. We find excellent agreement between spectrally and dynamically inferred mass ratios for the $\sim$600 binaries in which a dynamical mass ratio can be measured from multi-epoch radial velocities. We obtain full orbital solutions for 64 systems, including 14 close binaries within hierarchical triples. We make available catalogs of stellar parameters, abundances, mass ratios, and orbital parameters.

Gas-rich minor mergers contribute significantly to the gas reservoir of early-type galaxies (ETGs) at low redshift, yet the star formation efficiency (SFE; the star formation rate divided by the molecular gas mass) appears to be strongly suppressed following some of these events, in contrast to the more well-known merger-driven starbursts. We present observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of six ETGs, which have each recently undergone a gas-rich minor merger, as evidenced by their disturbed stellar morphologies. These galaxies were selected because they exhibit extremely low SFEs. We use the resolving power of ALMA to study the morphology and kinematics of the molecular gas. The majority of our galaxies exhibit spatial and kinematical irregularities, such as detached gas clouds, warps, and other asymmetries. These asymmetries support the interpretation that the suppression of the SFE is caused by dynamical effects stabilizing the gas against gravitational collapse. Through kinematic modelling we derive high velocity dispersions and Toomre Q stability parameters for the gas, but caution that such measurements in edge-on galaxies suffer from degeneracies. We estimate merger ages to be about 100~Myr based on the observed disturbances in the gas distribution. Furthermore, we determine that these galaxies lie, on average, two orders of magnitude below the Kennicutt-Schmidt relation for star-forming galaxies as well as below the relation for relaxed ETGs. We discuss potential dynamical processes responsible for this strong suppression of star formation surface density at fixed molecular gas surface density.

The distributions of a galaxy's gas and stars in chemical space encodes a tremendous amount of information about that galaxy's physical properties and assembly history. However, present methods for extracting information from chemical distributions are based either on coarse averages measured over galactic scales (e.g., metallicity gradients) or on searching for clusters in chemical space that can be identified with individual star clusters or gas clouds on $\sim 1$ pc scales. These approaches discard most of the information, because in galaxies gas and young stars are observed to be distributed fractally, with correlations on all scales, and the same is likely to be true of metals. In this paper we introduce a first theoretical model, based on stochastically-forced diffusion, capable of predicting the multi-scale statistics of metal fields. We derive the variance, correlation function, and power spectrum of the metal distribution from first principles, and determine how these quantities depend on elements' astrophysical origin sites and on the large-scale properties of galaxies. Among other results, we explain for the first time why the typical abundance scatter observed in the interstellar media of nearby galaxies is $\approx 0.1$ dex, and we predict that this scatter will be correlated on spatial scales of $\sim 0.5-1$ kpc, and over time scales of $\sim 100-300$ Myr. We discuss the implications of our results for future chemical tagging studies.

Core helium burning primary red clump (RC) stars are evolved red giant stars which are excellent standard candles. As such, these stars are routinely used to map the Milky Way or determine the distance to other galaxies among other things. However distinguishing RC stars from their less evolved precursors, namely red giant branch (RGB) stars, is still a difficult challenge and has been deemed the domain of asteroseismology. In this work, we use a sample of 1,676 RGB and RC stars which have both single epoch infrared spectra from the APOGEE survey and asteroseismic parameters and classification to show that the spectra alone can be used to (1) predict asteroseismic parameters with precision high enough to (2) distinguish core helium burning RC from other giant stars with less than 2% contamination. This will not only allow for a clean selection of a large number of standard candles across our own and other galaxies from spectroscopic surveys, but also will remove one of the primary roadblocks for stellar evolution studies of mixing and mass loss in red giant stars.

The observable spectrum of an unresolved binary star system is a superposition of two single-star spectra. Even without a detectable velocity offset between the two stellar components, the combined spectrum of a binary system is in general different from that of either component, and fitting it with single-star models may yield inaccurate stellar parameters and abundances. We perform simple experiments with synthetic spectra to investigate the effect of unresolved main-sequence binaries on spectral fitting, modeling spectra similar to those collected by the APOGEE, GALAH, and LAMOST surveys. We find that fitting unresolved binaries with single-star models introduces systematic biases in the derived stellar parameters and abundances that are modest but certainly not negligible, with typical systematic errors of $300\,\rm K$ in $T_{\rm eff}$, 0.1 dex in $\log g$, and 0.1 dex in $[\rm Fe/H]$ for APOGEE-like spectra of solar-type stars. These biases are smaller for spectra at optical wavelengths than in the near-infrared. We show that biases can be corrected by fitting spectra with a binary model, which adds only two labels to the fit and includes single-star models as a special case. Our model provides a promising new method to constrain the Galactic binary population, including systems with single-epoch spectra and no detectable velocity offset between the two stars.

The LAMOST survey has acquired low-resolution spectra (R=1,800) for 5 million stars across the Milky Way, far more than any current stellar survey at a corresponding or higher spectral resolution. It is often assumed that only very few elemental abundances can be measured from such low-resolution spectra, limiting their utility for Galactic archaeology studies. However, Ting et al. (2017) used ab initio models to argue that low-resolution spectra should enable precision measurements of many elemental abundances, at least in theory. Here we verify this claim in practice by measuring the relative abundances of 14 elements from LAMOST spectra with a precision of $\lesssim$ 0.1 dex for objects with ${\rm S/N}_{\rm LAMOST}$ > 30 (per pixel). We employ a spectral modeling method in which a data-driven model is combined with priors that the model gradient spectra should resemble ab initio spectral models. This approach assures that the data-driven abundance determinations draw on physically sensible features in the spectrum in their predictions and do not just exploit astrophysical correlations among abundances. Our analysis is constrained to the number of elemental abundances measured in the APOGEE survey, which is the source of the training labels. Obtaining high quality/resolution spectra for a subset of LAMOST stars to measure more elemental abundances as training labels and then applying this method to the full LAMOST catalog will provide a sample with more than 20 elemental abundances that is an order of magnitude larger than current high-resolution surveys, substantially increasing the sample size for Galactic archaeology.

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 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.

Understanding the evolution of the Milky Way calls for the precise abundance determination of many elements in many stars. A common perception is that deriving more than a few elemental abundances ([Fe/H], [$\alpha$/Fe], perhaps [C/H], [N/H]) requires medium-to-high spectral resolution, $R \gtrsim 10,000$, mostly to overcome the effects of line blending. In recent work (Rix et al. 2016; Ting et al. 2016) we presented an efficient and practical way to model the full stellar spectrum, even when fitting a large number of stellar labels simultaneously. In this paper we quantify to what precision the abundances of many different elements can be recovered, as a function of spectroscopic resolution and wavelength range. In the limit of perfect spectral models and spectral normalization, we show that the precision of elemental abundances is nearly independent of resolution, for a fixed exposure time and number of detector pixels; low-resolution spectra simply afford much higher S/N per pixel and generally larger wavelength range in a single setting. We also show that estimates of most stellar labels are not strongly correlated with one another once $R \gtrsim 1,000$. Modest errors in the line spread function, as well as small radial velocity errors, do not affect these conclusions, and data driven models indicate that spectral (continuum) normalization can be achieved well enough in practice. These results, to be confirmed with an analysis of observed low-resolution data, open up new possibilities for the design of large spectroscopic stellar surveys and for the re-analysis of archival low-resolution datasets.

Galah is an ongoing high-resolution spectroscopic survey with the goal of disentangling the formation history of the Milky Way, using the fossil remnants of disrupted star formation sites which are now dispersed around the Galaxy. It is targeting a randomly selected, magnitude limited ($V \leq 14$) sample of stars, with the goal of observing one million objects. To date, 300,000 spectra have been obtained. Not all of them are correctly processed by parameter estimation pipelines and we need to know about them. We present a semi-automated classification scheme which identifies different types of peculiar spectral morphologies, in an effort to discover and flag potentially problematic spectra and thus help to preserve the integrity of the survey's results. To this end we employ a recently developed dimensionality reduction technique t-SNE (t-distributed Stochastic Neighbour Embedding), which enables us to represent the complex spectral morphology in a two-dimensional projection map while still preserving the properties of the local neighbourhoods of spectra. We find that the majority (178,483) of the 209,533 Galah spectra considered in this study represents normal single stars, whereas 31,050 peculiar and problematic spectra with very diverse spectral features pertaining to 28,579 stars are distributed into 10 classification categories: Hot stars, Cool metal-poor giants, Molecular absorption bands, Binary stars, H$\alpha$/H$\beta$ emission, H$\alpha$/H$\beta$ emission superimposed on absorption, H$\alpha$/H$\beta$ P-Cygni, H$\alpha$/H$\beta$ inverted P-Cygni, Lithium absorption, and Problematic. Classified spectra with supplementary information are presented in the catalogue, indicating candidates for follow-up observations and population studies of the short-lived phases of stellar evolution.

We examine the processes triggering star formation and Active Galactic Nucleus (AGN) activity in a sample of 25 low redshift ($z<0.13$) gas-rich galaxy mergers observed at milli-arcsecond resolution with Very Long Baseline Interferometry as part of the mJy Imaging VLBA Exploration at 20cm (mJIVE-20) survey. The high ($>10^7$ K) brightness temperature required for an mJIVE-20 detection allows us to unambiguously identify the radio AGN in our sample. We find three such objects. Our VLBI AGN identifications are classified as Seyferts or LINERs in narrow line optical diagnostic plots; mid-infrared colours of our targets and the comparison of H$\alpha$ star formation rates with integrated radio luminosity are also consistent with the VLBI identifications. We reconstruct star formation histories in our galaxies using optical and UV photometry, and find that these radio AGN are not triggered promptly in the merger process, consistent with previous findings for non-VLBI samples of radio AGN. This delay can significantly limit the efficiency of feedback by radio AGN triggered in galaxy mergers. We find that radio AGN hosts have lower star formation rates than non-AGN radio-selected galaxies at the same starburst age. Conventional and VLBI radio imaging shows these AGN to be compact on arcsecond scales. Our modeling suggests that the actual sizes of AGN-inflated radio lobes may be much larger than this, but these are too faint to be detected in existing observations. Deep radio imaging is required to map out the true extent of the AGN, and to determine whether the low star formation rates in radio AGN hosts are a result of the special conditions required for radio jet triggering, or the effect of AGN feedback.

Stellar spectra depend on the stellar parameters and on dozens of photospheric elemental abundances. Simultaneous fitting of these $\mathcal{N}\sim 10-40$ model labels to observed spectra has been deemed unfeasible, because the number of ab initio spectral model grid calculations scales exponentially with $\mathcal{N}$. We suggest instead the construction of a polynomial spectral model (PSM) of order $\mathcal{O}$ for the model flux at each wavelength. Building this approximation requires a minimum of only ${\mathcal{N}+\mathcal{O}\choose\mathcal{O}}$ calculations: e.g. a quadratic spectral model ($\mathcal{O}=2$) to fit $\mathcal{N}=20$ labels simultaneously, can be constructed from as few as $231$ ab initio spectral model calculations; in practice, a somewhat larger number ($\sim 300-1000$) of randomly chosen models lead to a better performing PSM. Such a PSM can be a good approximation only over a portion of label space, which will vary case by case. Yet, taking the APOGEE survey as an example, a single quadratic PSM provides a remarkably good approximation to the exact ab initio spectral models across much of this survey: for random labels within that survey the PSM approximates the flux to within $10^{-3}$, and recovers the abundances to within $\sim 0.02$ dex rms of the exact models. This enormous speed-up enables the simultaneous many-label fitting of spectra with computationally expensive ab initio models for stellar spectra, such as non-LTE models. A PSM also enables the simultaneous fitting of observational parameters, such as the spectrum's continuum or line-spread function.

Stellar spectra are often modeled and fit by interpolating within a rectilinear grid of synthetic spectra to derive the stars' labels: stellar parameters and elemental abundances. However, the number of synthetic spectra needed for a rectilinear grid grows exponentially with the label space dimensions, precluding the simultaneous and self-consistent fitting of more than a few elemental abundances. Shortcuts such as fitting subsets of labels separately can introduce unknown systematics and do not produce correct error covariances in the derived labels. In this paper we present a new approach -- CHAT (Convex Hull Adaptive Tessellation) -- which includes several new ideas for inexpensively generating a sufficient stellar synthetic library, using linear algebra and the concept of an adaptive, data-driven grid. A convex hull approximates the region where the data lie in the label space. A variety of tests with mock datasets demonstrate that CHAT can reduce the number of required synthetic model calculations by three orders of magnitude in an 8D label space. The reduction will be even larger for higher-dimensional label spaces. In CHAT the computational effort increases only linearly with the number of labels that are fit simultaneously. Around each of these grid points in label space an approximate synthetic spectrum can be generated through linear expansion using a set of "gradient spectra" that represent flux derivatives at every wavelength point with respect to all labels. These techniques provide new opportunities to fit the full stellar spectra from large surveys with 15-30 labels simultaneously.

Stars born from the same molecular cloud should be nearly homogeneous in their element abundances. The concept of chemical tagging is to identify members of disrupted clusters by their clustering in element abundance space. Chemical tagging requires large samples of stars with precise abundances for many individual elements. With uncertainties of $\sigma_{[X/{\rm Fe}]}$ and $\sigma_{\rm [Fe/H]} \simeq 0.05$ for 10 elements measured for $>10^4$ stars, the APOGEE DR12 spectra may be the first well-suited data set to put this idea into practice. We find that even APOGEE data offer only $\sim 500$ independent volume elements in the 10-dimensional abundance space, when we focus on the $\alpha$-enhanced Galactic disk. We develop and apply a new algorithm to search for chemically homogeneous sets of stars against a dominant background. By injecting star clusters into the APOGEE data set we show that chemically homogeneous clusters with masses $\gtrsim 3 \times 10^7 \, {\rm M}_\odot$ would be easily detectable and yet no such signal is seen in the data. By generalizing this approach, we put a first abundance-based constraint on the cluster mass function for the old disk stars in the Milky Way.

In this work we present IRAM-30m telescope observations of a sample of bulge-dominated galaxies with large dust lanes, which have had a recent minor merger. We find these galaxies are very gas rich, with H2 masses between 4x10^8 and 2x10^10 Msun. We use these molecular gas masses, combined with atomic gas masses from an accompanying paper, to calculate gas-to-dust and gas-to-stellar mass ratios. The gas-to-dust ratios of our sample objects vary widely (between ~50 and 750), suggesting many objects have low gas-phase metallicities, and thus that the gas has been accreted through a recent merger with a lower mass companion. We calculate the implied minor companion masses and gas fractions, finding a median predicted stellar mass ratio of ~40:1. The minor companion likely had masses between ~10^7 - 10^10 Msun. The implied merger mass ratios are consistent with the expectation for low redshift gas-rich mergers from simulations. We then go on to present evidence that (no matter which star-formation rate indicator is used) our sample objects have very low star-formation efficiencies (star-formation rate per unit gas mass), lower even than the early-type galaxies from ATLAS3D which already show a suppression. This suggests that minor mergers can actually suppress star-formation activity. We discuss mechanisms that could cause such a suppression, include dynamical effects induced by the minor merger.

It is now well-established that the elemental abundance patterns of stars holds key clues not only to their formation but also to the assembly histories of galaxies. One of the most exciting possibilities is the use of stellar abundance patterns as "chemical tags" to identify stars that were born in the same molecular cloud. In this paper we assess the prospects of chemical tagging as a function of several key underlying parameters. We show that in the fiducial case of $10^4$ distinct cells in chemical space and $10^5-10^6$ stars in the survey, one can expect to detect $\sim 10^2-10^3$ groups that are $\ge 5\sigma$ overdensities in the chemical space. However, we find that even very large overdensities in chemical space do not guarantee that the overdensity is due to a single set of stars from a common birth cloud. In fact, for our fiducial model parameters, the typical $5\sigma$ overdensity is comprised of stars from a wide range of clusters with the most dominant cluster contributing only 25% of the stars. The most important factors limiting the identification of disrupted clusters via chemical tagging are the number of chemical cells in the chemical space and the survey sampling rate of the underlying stellar population. Both of these factors can be improved through strategic observational plans. While recovering individual clusters through chemical tagging may prove challenging, we show, in agreement with previous work, that different CMFs imprint different degrees of clumpiness in chemical space. These differences provide the opportunity to statistically reconstruct the slope and high mass cutoff of CMF and its evolution through cosmic time.

In a paper by Villata (2011), the possibility of a repulsive gravitational interaction between antimatter and ordinary matter was discussed. The author argued that this anti-gravity can be regarded as a prediction of general relativity, under the assumption of CPT symmetry. Stringent experimental constraints have been established against such a suggestion. The measurement of free-fall accelerations of various nuclei by the Eot-Wash group and searches for equivalence principle violation through the gravitational splitting in kaon physics consistently establish null results on any difference between the gravitational behaviour of antimatter and ordinary matter. The original arguments against antigravity were questioned by Nieto & Goldman (1991). In the light of new experiments as well as theoretical developments in the past 20 years, some of Nieto & Goldman's concerns have been addressed. While a precise measurement of the free-fall acceleration of antihydrogen will eventually lay this issue to rest, the purpose of this short letter is to argue that the substitution of dark energy with anti-gravity by antimatter, as suggested by Villata, is highly unlikely.

The GALAH survey is a large high-resolution spectroscopic survey using the newly commissioned HERMES spectrograph on the Anglo-Australian Telescope. The HERMES spectrograph provides high-resolution (R ~28,000) spectra in four passbands for 392 stars simultaneously over a 2 degree field of view. The goal of the survey is to unravel the formation and evolutionary history of the Milky Way, using fossil remnants of ancient star formation events which have been disrupted and are now dispersed throughout the Galaxy. Chemical tagging seeks to identify such dispersed remnants solely from their common and unique chemical signatures; these groups are unidentifiable from their spatial, photometric or kinematic properties. To carry out chemical tagging, the GALAH survey will acquire spectra for a million stars down to V~14. The HERMES spectra of FGK stars contain absorption lines from 29 elements including light proton-capture elements, alpha-elements, odd-Z elements, iron-peak elements and n-capture elements from the light and heavy s-process and the r-process. This paper describes the motivation and planned execution of the GALAH survey, and presents some results on the first-light performance of HERMES.

We use multi-wavelength (0.12 - 500 micron) photometry from Herschel-ATLAS, WISE, UKIDSS, SDSS and GALEX, to study 23 nearby spheroidal galaxies with prominent dust lanes (DLSGs). DLSGs are considered to be remnants of recent minor mergers, making them ideal laboratories for studying both the interstellar medium (ISM) of spheroids and minor-merger-driven star formation in the nearby Universe. The DLSGs exhibit star formation rates (SFRs) between 0.01 and 10 MSun yr^-1, with a median of 0.26 MSun yr^-1 (a factor of 3.5 greater than the average SG). The median dust mass, dust-to-stellar mass ratio and dust temperature in these galaxies are around 10^7.6 MSun yr^-1, ~0.05% and ~19.5 K respectively. The dust masses are at least a factor of 50 greater than that expected from stellar mass loss and, like the SFRs, show no correlation with galaxy luminosity, suggesting that both the ISM and the star formation have external drivers. Adopting literature gas-to-dust ratios and star formation histories derived from fits to the panchromatic photometry, we estimate that the median current and initial gas-to-stellar mass ratios in these systems are ~4% and ~7% respectively. If, as indicated by recent work, minor mergers that drive star formation in spheroids with (NUV-r)>3.8 (the colour range of our DLSGs) have stellar mass ratios between 1:6 and 1:10, then the satellite gas fractions are likely >50%.

We present a rigorous and practical way of constraining the Galactic potential based on the phase-space information for many individual stars. Such an approach is needed to dynamically model the data from ongoing spectroscopic surveys of the Galaxy and in the future Gaia. This approach describes the orbit distribution of stars by a family of parametrized distribution function (DF) proposed by McMillan and Binney, which are based on actions. We find that these parametrized DFs are flexible enough to capture well the observed phase-space distributions of individual abundance-selected Galactic subpopulations of stars (`mono-abundance populations') for a disc-like gravitational potential, which enables independent dynamical constraints from each of the Galactic mono-abundance populations. We lay out a statistically rigorous way to constrain the Galactic potential parameters by constructing the joint likelihood of potential and DF parameters, and subsequently marginalizing over the DF parameters. This approach explicitly incorporates the spatial selection function inherent to all Galactic surveys, and can account for the uncertainties of the individual position--velocity observations. On that basis, we study the precision of the parameters of the Galactic potential that can be reached with various sample sizes and realistic spatial selection functions. By creating mock samples from the DF, we show that, even under a restrictive and realistic spatial selection function, given a two-parameter gravitational potential, one can recover the true potential parameters to a few per cent with sample sizes of a few thousands. The assumptions of axisymmetry, of DFs that are smooth in the actions and of no time variation remain important limitations in our current study.

We present detailed elemental abundances of 12 subgiants in the open cluster IC 4756 including Na, Al, Mg, Si, Ca, Ti, Cr, Ni, Fe, Zn and Ba. We measure the cluster to have [Fe/H] = -0.01 +/- 0.10. Most of the measured star-to-star [X/H] abundance variation is below sigma < 0.03, as expected from a coeval stellar population preserving natal abundance patterns, supporting the use of elemental abundances as a probe to reconstruct dispersed clusters. We find discrepancies between Cr I and Cr II abundances as well as between Ti I and Ti II abundances, where the ionized abundances are larger by about 0.2 dex. This follows other such studies which demonstrate the effects of overionization in cool stars. IC 4756 are supersolar in Mg, Si, Na and Al, but are solar in the other elements. The fact that IC 4756 is supersolar in some alpha-elements (Mg, Si) but solar in the others (Ca, Ti) suggests that the production of alpha-elements is not simply one dimensional and could be exploited for chemical tagging.

We present the second of two papers concerning the origin and evolution of local early-type galaxies exhibiting dust features. We use optical and radio data to examine the nature of active galactic nucleus (AGN) activity in these objects, and compare these with a carefully constructed control sample. We find that dust lane early-type galaxies are much more likely to host emission-line AGN than the control sample galaxies. Moreover, there is a strong correlation between radio and emission-line AGN activity in dust lane early-types, but not the control sample. Dust lane early-type galaxies show the same distribution of AGN properties in rich and poor environments, suggesting a similar triggering mechanism. By contrast, this is not the case for early-types with no dust features. These findings strongly suggest that dust lane early-type galaxies are starburst systems formed in gas-rich mergers. Further evidence in support of this scenario is provided by enhanced star formation and black hole accretion rates in these objects. Dust lane early-types therefore represent an evolutionary stage between starbursting and quiescent galaxies. In these objects, the AGN has already been triggered but has not as yet completely destroyed the gas reservoir required for star formation.

We study dust and associated molecular gas in 352 nearby early-type galaxies (ETGs) with prominent dust lanes. 65% of these `dusty ETGs' (D-ETGs) are morphologically disturbed, suggesting a merger origin. This is consistent with the D-ETGs residing in lower density environments compared to the controls drawn from the general ETG population. 80% of D-ETGs inhabit the field (compared to 60% of the controls) and <2% inhabit clusters (compared to 10% of the controls). Compared to the controls, D-ETGs exhibit bluer UV-optical colours (indicating enhanced star formation) and an AGN fraction that is more than an order of magnitude greater (indicating higher incidence of nuclear activity). The clumpy dust mass residing in large-scale features is estimated, using the SDSS r-band images, to be 10^{4.5}-10^{6.5} MSun. A comparison to the total (clumpy + diffuse) dust masses- calculated using the far-IR fluxes of 15% of the D-ETGs that are detected by the IRAS- indicates that only ~20% of the dust resides in these large-scale features. The dust masses are several times larger than the maximum value expected from stellar mass loss, ruling out an internal origin. The dust content shows no correlation with the blue luminosity, indicating that it is not related to a galactic scale cooling flow. No correlation is found with the age of the recent starburst, suggesting that the dust is accreted directly in the merger rather than being produced in situ by the triggered star formation. Using molecular gas-to-dust ratios of ETGs in the literature we estimate that the median current and initial molecular gas fraction are ~1.3% and ~4%, respectively. Recent work suggests that the merger activity in nearby ETGs largely involves minor mergers (mass ratios between 1:10 and 1:4). If the IRAS-detected D-ETGs form via this channel, then the original gas fractions of the accreted satellites are 20%-44%. [Abridged]

In preparation for the HERMES chemical tagging survey of about a million Galactic FGK stars, we estimate the number of independent dimensions of the space defined by the stellar chemical element abundances [X/Fe]. [...] We explore abundances in several environments, including solar neighbourhood thin/thick disk stars, halo metal-poor stars, globular clusters, open clusters, the Large Magellanic Cloud and the Fornax dwarf spheroidal galaxy. [...] We find that, especially at low metallicity, the production of r-process elements is likely to be associated with the production of alpha-elements. This may support the core-collapse supernovae as the r-process site. We also verify the over-abundances of light s-process elements at low metallicity, and find that the relative contribution decreases at higher metallicity, which suggests that this lighter elements primary process may be associated with massive stars. [...] Our analysis reveals two types of core-collapse supernovae: one produces mainly alpha-elements, the other produces both alpha-elements and Fe-peak elements with a large enhancement of heavy Fe-peak elements which may be the contribution from hypernovae. [...] The extra contribution from low mass AGB stars at high metallicity compensates the dimension loss due to the homogenization of the core-collapse supernovae ejecta. [...] the number of independent dimensions of the [X/Fe]+[Fe/H] chemical space in the solar neighbourhood for HERMES is about 8 to 9. Comparing fainter galaxies and the solar neighbourhood, we find that the chemical space for fainter galaxies such as Fornax and the Large Magellanic Cloud has a higher dimensionality. This is consistent with the slower star formation history of fainter galaxies. [...]