Large multi-object spectroscopic surveys require automated algorithms to optimise their observing strategy. One of the most ambitious upcoming spectroscopic surveys is the 4MOST survey. The 4MOST survey facility is a fibre-fed spectroscopic instrument on the VISTA telescope with a large enough field of view to survey a large fraction of the southern sky within a few years. Several Galactic and extragalactic surveys will be carried out simultaneously, so the combined target density will strongly vary. In this paper, we describe a new tiling algorithm that can naturally deal with the large target density variations on the sky and which automatically handles the different exposure times of targets. The tiling pattern is modelled as a marked point process, which is characterised by a probability density that integrates the requirements imposed by the 4MOST survey. The optimal tilling pattern with respect to the defined model is estimated by the tiles configuration that maximises the proposed probability density. In order to achieve this maximisation a simulated annealing algorithm is implemented. The algorithm automatically finds an optimal tiling pattern and assigns a tentative sky brightness condition and exposure time for each tile, while minimising the total execution time that is needed to observe the list of targets in the combined input catalogue of all surveys. Hence, the algorithm maximises the long-term observing efficiency and provides an optimal tiling solution for the survey. While designed for the 4MOST survey, the algorithm is flexible and can with simple modifications be applied to any other multi-object spectroscopic survey.

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

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.

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.

Recent and on-going large ground-based multi-object spectroscopic surveys allow to significantly increase the sample of spectroscopic binaries to get insight into their statistical properties. We investigate the repeated spectral observations of the Gaia-ESO Survey (GES) internal data release 5 to identify and characterize spectroscopic binaries with one visible component (SB1) in fields covering the discs, the bulge, the CoRot fields, and stellar clusters and associations. A statistical chi2-test is performed on spectra of the iDR5 sub-sample of approximately 43500 stars characterized by at least 2 observations and a S/N > 3. Our sample of RV variables is cleaned from contamination by pulsation/convection-induced variables using Gaia DR2 parallaxes and photometry. Monte-Carlo simulations using the SB9 catalogue of spectroscopic orbits allow to estimate our detection efficiency and to correct the SB1 rate to evaluate the GES SB1 binary fraction and its dependence with effective temperature and metallicity. We find 641 (resp., 803) FGK SB1 candidates at the 5 sigma (resp., 3 sigma) level. The orbital-period distribution is estimated from the RV standard-deviation distribution and reveals that the detected SB1 probe binaries with log(P[d]) < 4. We estimate the global GES SB1 fraction to be in the range 7-14% with a typical uncertainty of 4%. The GES SB1 frequency decreases with metallicity at a rate of -9+/-3%/dex in the metallicity range -2.7<FeH<+0.6. This anticorrelation is obtained with a confidence level higher than 93% on a homogeneous sample covering spectral types FGK and a large range of metallicities. When the present-day mass function is accounted for, this rate turns to 4+/-2%/dex with a confidence level higher than 88%. In addition we provide the variation of the SB1 fraction with metallicity separately for F, G, and K spectral types, as well as for dwarf and giant primaries.

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

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.

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.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

The signatures of the formation and evolution of a galaxy are imprinted in its stars. Their velocities, ages, and chemical compositions present major constraints on models of galaxy formation, and on various processes such as the gas inflows and outflows, the accretion of cold gas, radial migration, and the variability of star formation activity. Understanding the evolution of the Milky Way requires large observational datasets of stars via which these quantities can be determined accurately. This is the science driver of the 4MOST MIlky way Disc And BuLgE High-Resolution (4MIDABLE-HR) survey: to obtain high-resolution spectra at $R \sim 20\,000$ and to provide detailed elemental abundances for large samples of stars in the Galactic disc and bulge. High data quality will allow us to provide accurate spectroscopic diagnostics of two million stellar spectra: precise radial velocities; rotation; abundances of many elements, including those that are currently only accessible in the optical, such as Li, s-, and r-process; and multi-epoch spectra for a sub-sample of stars. Synergies with complementary missions like Gaia and TESS will provide masses, stellar ages and multiplicity, forming a multi-dimensional dataset that will allow us to explore and constrain the origin and structure of the Milky Way.

The mechanisms of the formation and evolution of the Milky Way are encoded in the orbits, chemistry and ages of its stars. With the 4MOST MIlky way Disk And BuLgE Low-Resolution Survey (4MIDABLE-LR) we aim to study kinematic and chemical substructures in the Milky Way disc and bulge region with samples of unprecedented size out to larger distances and greater precision than conceivable with Gaia alone or any other ongoing or planned survey. Gaia gives us the unique opportunity for target selection based almost entirely on parallax and magnitude range, hence increasing the efficiency in sampling larger Milky Way volumes with well-defined and effective selection functions. Our main goal is to provide a detailed chrono-chemo-kinematical extended map of our Galaxy and the largest Gaia follow-up down to $G = 19$ magnitudes (Vega). The complex nature of the disc components (for example, large target densities and highly structured extinction distribution in the Milky Way bulge and disc area), prompted us to develop a survey strategy with five main sub-surveys that are tailored to answer the still open questions about the assembly and evolution of our Galaxy, while taking full advantage of the Gaia data.

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

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.

Identifying and characterizing young stars in the Solar neighbourhood is essential to find and describe planets in the early stages of their evolution. This work seeks to identify nearby young stars showing a Lithium 6707.78$\,\AA$ absorption line in the GALAH survey. A robust, data-driven approach is used to search for corresponding templates in the pool of 434,215 measured dwarf spectra in the survey. It enables a model-free search for best-matching spectral templates for all stars, including M dwarfs with strong molecular absorption bands. 3147 stars have been found to have measurable Lithium: 1408 G and 892 K0-K5 dwarfs (EW(Li)$>$0.1$\,$\AA), 335 K5-K9 ($>$0.07$\,$\AA) and 512 M0$\sim$M4 dwarfs ($>$0.05$\,$\AA). Stars with such Lithium features are used to investigate the possibility of searching for young stars above the main sequence based merely on their parallaxes and broad-band photometry. Selection of young stars above the main sequence is highly effective for M dwarfs, moderately effective for K dwarfs and ineffective for G dwarfs. Using a combination of the Lithium information and the complete 6D kinematics from Gaia and GALAH, 305 new candidate moving group members have been found, 123 of which belong to the Scorpius-Centaurus association, 36 to the Pleiades and 25 to the Hyades clusters.

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.

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.

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.

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

V2492 Cyg is a young pre-main sequence star presenting repetitive brightness variations of significant amplitude (Delta R > 5 mag) whose physical origin has been ascribed to both extinction (UXor-type) and accretion (EXor-type) variability, although their mutual proportion has not been clarified yet. Recently, V2492 Cyg has reached a level of brightness ever registered in the period of its documented activity. Optical and near-infrared photometry and spectroscopy have been obtained in October 2016 and between March and July 2017. The source has remained bright until the end of May 2017, then it started to rapidly fade since the beginning of June at a rate of about 0.08 mag/day. On mid-July 2017 the source has reached the same low-brightness level as two years before. Extinction and mass accretion rate were derived by means of the luminosity of the brightest lines, in particular Halpha and Hbeta. A couple of optical high-resolution spectra are also presented to derive information on the gas kinematics. Visual extinction variations do not exceed a few magnitudes, while the mass accretion rate is estimated to vary from less than 10^-8 up to a few 10^-7 M_sun/yr. This latter is comparable to that estimated on the previous high-state in 2010, likely occurred under more severe extinction conditions. The combined analysis of the optical and near-infrared (NIR) observations extends to the present event the original suggestion that the V2492 Cyg variability is a combination of changing extinction and accretion.

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

The 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 apply the twin method to determine parallaxes to 232,545 stars of the RAVE survey using the parallaxes of Gaia DR1 as a reference. To search for twins in this large dataset, we apply the t-stochastic neighbour embedding t-SNE projection which distributes the data according to their spectral morphology on a two dimensional map. From this map we choose the twin candidates for which we calculate a chi^2 to select the best sets of twins. Our results show a competitive performance when compared to other model-dependent methods relying on stellar parameters and isochrones. The power of the method is shown by finding that the accuracy of our results is not significantly affected if the stars are normal or peculiar since the method is model free. We find twins for 60% of the RAVE sample which is not contained in TGAS or that have TGAS uncertainties which are larger than 20%. We could determine parallaxes with typical errors of 28%. We provide a complementary dataset for the RAVE stars not covered by TGAS, or that have TGAS uncertainties which are larger than 20%, with model-free parallaxes scaled to the Gaia measurements.

The Gaia-ESO Survey (GES) is a large spectroscopic survey that provides a unique opportunity to study the distribution of spectroscopic multiple systems among different populations of the Galaxy. We aim at detecting binarity/multiplicity for stars targeted by the GES from the analysis of the cross-correlation functions (CCFs) of the GES spectra with spectral templates. We develop a method based on the computation of the CCF successive derivatives to detect multiple peaks and determine their radial velocities, even when the peaks are strongly blended. The parameters of the detection of extrema (DOE) code have been optimized for each GES GIRAFFE and UVES setup to maximize detection. This code therefore allows to automatically detect multiple line spectroscopic binaries (SBn, n>1). We apply this method on the fourth GES internal data release and detect 354 SBn candidates (342 SB2, 11 SB3 and even one SB4), including only 9 SB2 known in the literature. This implies that about 98% of these SBn candidates are new (because of their faint visual magnitude that can reach V=19). Visual inspection of the SBn candidate spectra reveals that the most probable candidates have indeed a composite spectrum. Among SB2 candidates, an orbital solution could be computed for two previously unknown binaries: 06404608+0949173 (known as V642 Mon) in NGC 2264 and 19013257-0027338 in Berkeley 81. A detailed analysis of the unique SB4 (four peaks in the CCF) reveals that HD 74438 in the open cluster IC 2391 is a physically bound stellar quadruple system. The SB candidates belonging to stellar clusters are reviewed in detail to discard false detections. We warn against the use of atmospheric parameters for these system components rather than by SB-specific pipelines. Our implementation of an automatic detection of spectroscopic binaries within the GES has allowed an efficient discovery of many new multiple systems.

The orbit and physical parameters of the previously unsolved double-lined eclipsing binary Tyc 5227-1023-1, discovered during the search for RR Lyr variables candidate members of the Aquarius stream, are derived using high resolution \'echelle spectroscopy and $V$, $i^\prime$ photometry. Synthetic spectral analysis of both components has been performed, yielding metallicity [M/H] = $-0.63 \pm 0.11$ for both stars and the temperature of the secondary which is in close agreement with the one from the orbital solution, while the temperature of the primary is determined from photometry (T$_1 = 6350$ K). The masses and radii ($M_1 = 0.96 \pm 0.02, M_2 = 0.84 \pm 0.01$ M$_{\odot}$, $R_1 = 1.39 \pm 0.01, R_2 = 0.98 \pm 0.01$ R$_{\odot}$) reveal that both stars have already slightly evolved away from the Main Sequence band having an age of about $7$ Gyr, and the results of synthetic spectral analysis support the claim of co-rotation with the orbital motion. The radial velocity of the system is $-60 \pm 2$ km sec$^{-1}$ while its distance, computed from orbital parameters and the derived reddening $E_{B-V}$ = 0.053, is $496 \pm 35$ pc. Even though Tyc 5227-1023-1 was initially treated as a possible member of the Aquarius stream, the results presented here disagree with reported values for this ancient structure, and suggest a likely membership of the thick disk.

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 present the data reduction procedures being used by the GALAH survey, carried out with the HERMES fibre-fed, multi-object spectrograph on the 3.9~m Anglo-Australian Telescope. GALAH is a unique survey, targeting 1 million stars brighter than magnitude V=14 at a resolution of 28,000 with a goal to measure the abundances of 29 elements. Such a large number of high resolution spectra necessitates the development of a reduction pipeline optimized for speed, accuracy, and consistency. We outline the design and structure of the Iraf-based reduction pipeline that we developed, specifically for GALAH, to produce fully calibrated spectra aimed for subsequent stellar atmospheric parameter estimation. The pipeline takes advantage of existing Iraf routines and other readily available software so as to be simple to maintain, testable and reliable. A radial velocity and stellar atmospheric parameter estimator code is also presented, which is used for further data analysis and yields a useful verification of the reduction quality. We have used this estimator to quantify the data quality of GALAH for fibre cross-talk level ($\lesssim0.5$%) and scattered light ($\sim5$ counts in a typical 20 minutes exposure), resolution across the field, sky spectrum properties, wavelength solution reliability (better than $1$ $\mathrm{km\ s^{-1}}$ accuracy) and radial velocity precision.

Current ongoing stellar spectroscopic surveys (RAVE, GALAH, Gaia-ESO, LAMOST, APOGEE, Gaia) are mostly devoted to studying Galactic archaeology and structure of the Galaxy. But they allow for important auxiliary science: (i) Galactic interstellar medium can be studied in four dimensions (position in space + radial velocity) through weak but numerous diffuse insterstellar bands and atomic absorptions seen in spectra of background stars, (ii) emission spectra which are quite frequent even in field stars can serve as a good indicator of their youth, pointing e.g. to stars recently ejected from young stellar environments, (iii) astrometric solution of the photocenter of a binary to be obtained by Gaia can yield accurate masses when joined by spectroscopic information obtained serendipitously during a survey. These points are illustrated by first results from the first three surveys mentioned above. These hint at the near future: spectroscopic studies of the dynamics of the interstellar medium can identify and quantify Galactic fountains which may sustain star formation in the disk by entraining fresh gas from the halo; RAVE already provided a list of ~14,000 field stars with chromosperic emission in Ca II lines, to be supplemented by many more observations by Gaia in the same band, and by GALAH and Gaia-ESO observations of Balmer lines; several millions of astrometric binaries with periods up to a few years which are being observed by Gaia can yield accurate masses when supplemented with measurements from only a few high-quality ground based spectra.

We discuss the properties of H$\alpha$ emission stars across the sample of 22035 spectra from the Gaia-ESO Survey internal data release, observed with the GIRAFFE instrument and largely belonging to stars in young open clusters. Automated fits using two independent Gaussian profiles and a third component that accounts for the nebular emission allow us to discern distinct morphological types of H$\alpha$ line profiles with the introduction of a simplified classification scheme. All in all we find 3765 stars with intrinsic emission and sort their spectra into eight distinct morphological categories: single--component emission, emission blend, sharp emission peaks, double emission, P-Cygni, inverted P-Cygni, self--absorption, and emission in absorption. We have more than one observation for 1430 stars in our sample, thus allowing a quantitative discussion of the degree of variability of H$\alpha$ emission profiles, which is expected for young, active objects. We present a catalogue of stars with properties of their H$\alpha$ emission line profiles, morphological classification, analysis of variability with time and the supplementary information from the SIMBAD, VizieR, and ADS databases. The records in SIMBAD indicate the presence of H$\alpha$ emission for roughly 25% of all stars in our catalogue, while at least 305 of them have already been more thoroughly investigated according to the references in ADS. The most frequently identified morphological categories in our sample of spectra are emission blend (23%), emission in absorption (22%), and self--absorption (16%). Objects with repeated observations demonstrate that our classification into discrete categories is generally stable through time, but categories P-Cygni and self--absorption seem less stable, which is the consequence of discrete classification rules, as well as of the fundamental change in profile shape.

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 determine approximate numerical integrals of motion of 2D symmetric Hamiltonian systems. We detail for a few gravitational potentials the conditions under which quasi-integrals are obtained as polynomial series. We show that each of these potentials has a wide range of regular orbits that are accurately modelled with a unique approximate integral of motion.