The GALAH survey: An abundance, age, and kinematic inventory of the solar neighbourhood made with TGAS

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


Introduction
The first Gaia data release (Gaia Collaboration et al. 2016), which includes the Tycho-Gaia Astrometric Solution (TGAS), is a milestone of modern astronomy and has delivered positions, proper motions, and parallaxes for more than two million stars (Michalik et al. 2015;Lindegren et al. 2016). This has marked the beginning of the Gaia era, which will shed new light on our understanding of the formation and evolution of our Galaxy. The astrometric information delivered by Gaia will be particularly important when used in combination with spectroscopic quantities obtained using large ground-based surveys (e.g. Allende McMillan et al. 2018;Helmi et al. 2017;Kushniruk et al. 2017;Price-Whelan et al. 2017). Taken together, these data will prove particularly powerful in testing our models of Galactic assembly. The ensemble of GALAH (De Silva et al. 2015) and Gaia data have now begun to provide a high dimensionality mapping of the chemodynamical space of the nearby disk. In this study, we use the combination of GALAH and TGAS to describe the chemical, temporal, and kinematical distributions of nearby stars in the disk of the Milky Way.
The paper is organised as follows: We first introduce GALAH and discuss the strengths of combining this survey with the astrometric information provided by the Gaia satellite. In Sect. 2, we outline the observational strategy for GALAH. In Sect. 3 we explain our estimation of stellar properties, including stellar parameters, chemical composition, and ages. The results of our analysis are presented in Sect. 4 where we show the abundance and age trends for our sample and discuss the observed distribution of disk stars as a function of chemical composition, age and kinematics. In Sect. 5, we discuss the implications of our findings and make suggestions for further studies with GALAH (DR2) and Gaia (DR2) in the concluding section of the paper.

The GALAH survey
The GALAH survey is a ground-based, high-resolution stellar spectroscopic survey. It is executed with the High Efficiency and Resolution Multi-Element Spectrograph (HERMES) fed by the Two Degree Field (2dF) f/3.3 top end at the Anglo-Australian Telescope (Barden et al. 2010;Brzeski et al. 2011;Heijmans et al. 2012;Farrell et al. 2014;Sheinis et al. 2015). The overall scientific motivation for GALAH is presented in De Silva et al. (2015). The survey's primary goal is the chemical tagging experiment, as proposed by Freeman & Bland-Hawthorn (2002) and described in detail by Bland-Hawthorn et al. (2010). Chemical tagging offers the promise of linking stars that were born together via their chemical composition. As proposed by recent simulations (Ting et al. 2015(Ting et al. , 2016, this promise can be best explored with a high dimensionality in chemical space. Consequently, the spectrograph has been optimised to measure up to 30 different elements (more in very bright stars), covering a multitude of different nucleosynthesis channels, depending on the stellar type and evolution.
The GALAH survey selection function is both simple (magnitude limited) and will ensure that almost all stars observed by the GALAH survey are also measured by the Gaia satellite, see Sect. 2 and De . Our target selection limit of V ≤ 14 corresponds to Gaia 's peak performance with distance uncertainties for all GALAH stars expected to be better than 1%. Once GALAH is completed, this will provide both chemical and kinematical information for up to a million stars. These data will directly inform Galactic archaeology pursuits, enabling the empirical construction of the distribution function of stellar prop-erties and populations (chemical composition, age, position, orbits). The sample analysed here, is a first step in this direction.
For the analysis of the whole GALAH survey data, a combination of classical spectrum synthesis with Spectroscopy Made Easy (SME) by Piskunov & Valenti (2017) and with a datadriven propagation via The Cannon (Ness et al. 2015) is used. Prior to this study, data releases of GALAH (Martell et al. 2017), TESS-HERMES ) and K2-HERMES (Wittenmyer et al. 2018) were based only on spectroscopic input and provided the stellar properties obtained with the The Cannon. We stress that this work focuses on the first part of the usual GALAH analysis routine and is based only on a spectroscopic analysis with SME, but also includes photometric and astrometric input for the analysis. The stars used in this analysis will be a subset of the training set used for The Cannon in later data releases (GALAH DR2; Buder et al. 2018).

The strengths of combining GALAH and Gaia DR1 TGAS
Combining Gaia DR1 TGAS with high-resolution spectroscopy provides a variety of opportunities, ranging from an improved analysis of stellar parameters up to the expansion of the chemical space by the kinematical one.
Purely spectroscopic analyses may suffer from inaccuracies due to degeneracies in effective temperature (T eff ), surface gravity (log g), and chemical abundances. This is a consequence of simplified assumptions about stellar spectra and the subsequent construction of incomplete stellar models, for example assuming a 1D hydrostatic atmosphere and chemical compositions scaled with solar values and the metallicities. Previous high-resolution spectroscopic studies (see e.g. Bensby et al. 2014;Martell et al. 2017) have found inconsistencies between purely spectroscopic parameter estimates and those also based on photometric and asteroseismic information. Furthermore, many studies find that unphysical low surface gravities are estimated for G and K-type main sequence stars from spectroscopy alone Adibekyan et al. 2012). Cool dwarfs (T eff < 4500 K) are particularly challenging to study in the optical regime, because of the weakening of the singly ionised lines that are used to constrain the ionisation equilibrium, and due to the increasing influence of molecular blends as well as the failure of 1D LTE modelling; see e.g. Yong et al. (2004). Adding further (non-spectroscopic) information may alleviate these problems. Asteroseismic as well as interferometric and bolometric flux measurements for dwarf and turn-off stars are, however, still expensive and published values rare (especially for stars in the Southern hemisphere). Gaia DR2 and later releases will provide astrometric information for all observed GALAH stars and numerous stars that have been observed by other spectroscopic surveys, e.g. APOGEE (Majewski et al. 2017), RAVE (Kunder et al. 2017), Gaia-ESO (Gilmore et al. 2012;Randich et al. 2013), and LAMOST (Cui et al. 2012). The estimation of bolometric luminosities using both astrometric and photometric information will therefore be feasible for a large sample of stars in the near future.
Astrometric data, in combination with photometric and spectroscopic information, can also improve extinction estimates and narrow down uncertainties in the estimation of stellar ages from theoretical isochrones. Parallaxes and apparent magnitudes can also be used to identify binary systems with main sequence stars, which are not resolved by spectroscopy. This is important, because the companion contributes light to the spectrum and can hence significantly contaminate the analysis results (El-Badry et al. 2018a,b).
Prior to the first Gaia data release (DR1), the most notable observational chemodynamical studies were performed using the combination of the astrometric data from Hipparcos (ESA 1997;van Leeuwen 2007) and Tycho-2 (Høg et al. 2000) with additional observations by the Geneva-Copenhagen-Survey (Nordström et al. 2004;Casagrande et al. 2011) and highresolution follow-up observations (e.g. Bensby et al. 2014). Another approach, including post-correction of spectroscopic gravities has been adopted by Delgado Mena et al. (2017) for the HARPS-GTO sample. Large scale analyses have however been limited by the precision of astrometric measurements by Hipparcos to within the volume of a few hundred parsecs at most. With the new Gaia data, this volume is expanded to more than 1 kpc with DR1 and will be expanded even further with DR2, allowing the study of gradients or overdensities/groups in the chemodynamical space of the Milky Way disk and beyond.

Milky Way disk definitions and membership assignment
The disk is the most massive stellar component of the Milky Way. Numerous studies have observed the disk in the Solar neighbourhood. The pioneering studies by Yoshii (1982) and Gilmore & Reid (1983) found evidence for two thin and thick sub-populations in the disk based on stellar density distributions. Recent studies (e.g. Bovy et al. 2012Bovy et al. , 2016) find a structural continuity in thickness and kinematics, and the latter property has been shown to be a rather unreliable tracer of the disk subpopulations (Bensby et al. 2014). However, several seminal papers (e.g. Reddy et al. 2003;Fuhrmann 2011;Adibekyan et al. 2012;Bensby et al. 2014;Hayden et al. 2015) have established that the stellar disk consists of (at least) two major components in chemical space and age, commonly adopted as old, α-enhanced, metal-poor thick and the young thin disk with solar-like αenhancement at metallicities of −0.7 < [Fe/H] < +0.5. However, the bimodality between these two populations in the high-α metal-rich regime has been shown to become less or not significant and is still contentious, based on the chosen approaches and population cuts used for the definition of disk populations. (Adibekyan et al. 2011) even claimed a third sub-population in this regime. Some of the recent studies using chemistry assume the existence of two distinct populations in α-enhancement up to the most metal-rich stars. In these studies, the metal-rich stars are cut into high and low sequence memberships rather arbitrarily; either by eye or with rather fiducial straight lines (e.g. Adibekyan et al. (2012) Fig. 7, (Recio-Blanco et al. 2014 Fig. 12 or Hayden et al. (2017) Fig. 1). Other approaches to separate the α-sequences or stellar populations, e.g. kinematically (Bensby et al. 2014) or via age (Haywood et al. 2013) result in different separations. A consistent measure or definition to separate both α-sequences especially in the metal-rich regime remains elusive. For a discussion on combining chemistry and kinematics to separate the two α-sequences see e.g. Haywood et al. (2013) and Duong et al. (2018). For a more detailed overview regarding the definition of the stellar disks, we refer the reader to Martig et al. (2016b) and references therein.
Although element abundances are easier to determine than stellar ages, their surface abundances and measurements are subject to changes due to processes within the atmosphere of a star. Stellar ages are therefore the most promising tracer of stellar evolution and populations (Haywood et al. 2013;Bensby et al. 2014;Ness et al. 2016;Ho et al. 2017b). We note that, most recently, Hayden et al. (2017) investigated abundance sequences as a function of age, but with a representation as function of age ranges, starting with only old stars and subsequently including more younger ones (see their Fig. 3).
Until now, however, most authors suggested that this issue should be revisited when a larger, homogeneous, and less biased sample is available. With the observations obtained by the GALAH survey, we are now able to investigate the abundance sequences with a significantly larger and homogeneous data set with a rather simple selection function.

Observations
The GALAH survey collects data with HERMES, which can observe up to 360 science targets at the same time plus 40 fibres allocated for sky and guide stars (Sheinis et al. 2015;Heijmans et al. 2012;Brzeski et al. 2011;Barden et al. 2010). The selection of targets and observational setup are explained in detail by De Silva et al. (2015) and Martell et al. (2017). The observations used in this study were carried out between November 2013 and September 2016 with the lower of the two resolution modes (λ/∆λ ∼ 28000) with higher throughput, covering the four arms of HERMES, i.e. blue (4716 − 4896 Å including H β ), green (5650 − 5868 Å), red (6480 − 6734 Å including H α ) as well as the near infrared (7694 − 7876 Å, including the oxygen triplet).
The initial simple selection function of the GALAH survey was achieved with a random selection of stars within the limiting magnitudes 12 < V < 14 derived from 2MASS photometry (De Silva et al. 2015). To ensure a large overlap with TGAS (Michalik et al. 2015), our team added special bright fields (9 < V < 12) including a large number of stars in the Tycho-2 catalog (Høg et al. 2000) which were brighter than the nominal GALAH range (Martell et al. 2017). The exposure times were chosen to achieve a signal-to-noise ratio (S /N) of 100 per resolution element in the green channel / arm; 1 hour for main survey targets in optimal observing conditions, often longer in suboptimal observing conditions. The spectra are reduced with the GALAH pipeline (Kos et al. 2017), including initial estimates of T eff , log g, [Fe/H], and radial velocities (v rad ).
The GALAH+TGAS sample, observed until September 2016, consists of 23096 stars, covering mainly the spectral types F-K from pre-main sequence up to evolved asymptotic giant branch stars. For an overview of the sample, the spectroscopic parameters are depicted in Fig. 1, coloured by the parallax precision from TGAS. We note that the shown parameters are estimated as part of this study (see Sect. 3.1). The most precise parallaxes (σ( )/ ≤ 0.05) are available for main sequence stars cooler than 6000 K, decent parallaxes (σ( )/ ≤ 0.3) for most dwarfs and the lower luminosity end of the red giant branch. As expected by the magnitude constraints of TGAS as well as GALAH (De Silva et al. 2015), most of the overlap consists of dwarfs and turn-off stars (62 %) which also have smaller relative parallax uncertainties than the more distant giants, see Fig. 2. Cool evolved giants as well as hot turn-off stars have the least precise parallaxes of the GALAH+TGAS overlap because of their larger distances and are hence not included in the online tables.
In this work, we limit the sample for the analyses to dwarfs and turn-off stars (T eff ≥ 5500 K or log g ≥ 3.8 dex, see dashed line in Fig. 1) with relative parallax uncertainties smaller than 30 %. This allows the best estimation of ages from isochrones as well as reliable distance and kinematical information and avoids possible systematic differences in the analysis due to the different evolutionary stages of the stars. Evolutionary effects, such Article number, page 3 of 30 A&A proofs: manuscript no. Buder_GnT_v2 Kiel diagram with effective temperature T eff and surface gravity log g for the complete GALAH+TGAS overlap. The spectroscopic parameters are results of the analysis in Sect. 3.1. Colour indicates the relative parallax error. Most precise parallaxes are measured for the cool main sequence stars and the parallax precision decreases both towards the turn-off sequence and even more drastically towards evolved giants, which are the most distant stars in the sample. Dotted and dashed black lines indicate the limits to neglect hot stars (T eff > 6900 K) and giants (T eff < 5500 K and log g < 3.8 dex) for the subsequent analysis, respectively. See text for details on the exclusion of stars. Histograms of relative parallax uncertainties for both giants (top panel with T eff < 5500 K and log g < 3.8 dex) and dwarfs (lower panel, including main-sequence and turn-off stars). The majority of the GALAH+TGAS overlap consists of dwarfs. Their mean parallax precision is in the order of 10%, while giants parallaxes are less precise with most uncertainties above 20%. The histograms are truncated at σ( ) = for readability. All stars with parallax uncertainties larger than the parallax itself are hence contained in the last bin. as atomic diffusion, have been studied both with observations of clusters (see e.g. studies of the open cluster M67 by Önehag et al. (2014);Bertelli Motta et al. (2017); Gao et al. (2018)) as well as a theoretical predictions (Dotter et al. 2017) and are beyond the scope of this paper.
In addition to removing 8740 giant stars and 7674 stars with parallax uncertainties above 30%, we exclude some stars after a visual inspection of the spectra and using our quality analysis (explained in Sect. 3.1). We construct a final sample with reliable stellar parameters and element abundances. We neglect 54 stars with emission lines, 926 stars with bad spectra or reductions, 448 double-lined spectroscopic binary stars, 338 photometric binaries (see Sect. 3.4), 3429 stars with broadening velocities above 30 km/s (mostly hot turn-off stars with unbroken degeneracies of broadening velocity and stellar parameters with the GALAH setup), 1390 stars with T eff > 6900 K (for which we have not been able to measure element abundances) and 1048 stars with S /N below 25 in the green channel. Note that the different groups of excluded stars defined above are overlapping with each other.
For the final selected sample of 7066 stars, the majority of the individual S /N vary between 25 to 200, see Fig. 3. Most of the stars have a higher S /N than the targeted nominal survey value for the green channel. We note that the S /N of the blue arm is lower than in the others. For abundances measured within this arm, like Zn, we also estimate typically lower precision, see Sect. A.6.

Analysis
Our analysis combines the use of information derived from our GALAH spectra with additional photometric and astrometric measurements to achieve the best possible parameter estimation. We validate our analysis in a manner similar to other large-scale stellar surveys, such as APOGEE (Majewski et al. 2017;Abolfathi et al. 2018;García Pérez et al. 2016) or Gaia-ESO (Gilmore et al. 2012;Randich et al. 2013;Smiljanic et al. 2014;Pancino et al. 2017), using a set of well-studied stars, including the so called Gaia FGK benchmark stars (hereafter GBS, see Heiter et al. 2015a;Jofré et al. 2014). The stellar parameters T eff and log g of the GBS have been derived from direct observables: angular diameters, bolometric fluxes, and parallaxes, and are thus less model-dependent. They therefore provide reference parameters that do not suffer from the same model dependence as isolated spectroscopy. Among others, Schönrich & Bergemann (2014) and Bensby et al. (2014) showed the strength of combining spectroscopy and external information. The latter applied this approach for a sample of 714 nearby dwarfs with high accuracy astrometric parallaxes (van Leeuwen 2007) from the Hip-Article number, page 4 of 30 parcos mission. We will use their sample as a reference for this study, because the spectral analysis was performed in a similar way, including the anchoring of surface gravity to astrometric information. We stress that their study was performed with higher quality spectra (both regarding the spectral resolution and S /N) which allowed a higher precision on measurements to be achieved.

Stellar parameter determination
By using the fundamental relation between surface gravity, stellar mass, effective temperature, and bolometric luminosity the degeneracies with log g and other spectroscopically determined stellar parameters are effectively broken. The thereby improved values of T eff and log g leads to improved estimates of metallicities. Using broadband photometry (apparent magnitudes K S and inferred bolometric corrections BC K S as well as extinction A K S ) in combination with parallaxes or distances D , it is possible to estimate the bolometric magnitudes (M bol ) and luminosities (L bol ) to high precision and accuracy (see e.g. Alonso et al. 1995;Nissen et al. 1997;Bensby et al. 2014): The nominal values for the Sun (used in Eqs. 1 and 2) of T eff, = 5772 K, log(g ) = 4.438 dex, and M Bol, = 4.74 mag are taken from Prša et al. (2016).
Any filter with available bolometric corrections can be used for the computation of L bol ; the V band is commonly used for nearby stars. However, our GALAH data set also contains stars with substantial reddening and published catalogues of V band magnitudes, such as APASS (Henden et al. 2016), have multiple input sources, affecting the homogeneity of the data. We therefore decide in favour of using the K S band as given by 2MASS (Cutri et al. 2003), available for all our targets. Bolometric corrections BC = BC(T eff , log g, [Fe/H], E(B − V)) were interpolated with the grids from Casagrande & VandenBerg (2014). Distances were taken from Astraatmadja & Bailer-Jones (2016), using a Milky Way model as Bayesian prior. For attenuation, we use the RJCE method A K = A K (K S , W2) by Majewski et al. (2011);Zasowski et al. (2013). If K S or W2 could not be used, we use the approximation A K ∼ 0.38E(B − V) estimated by Savage & Mathis (1979) Schlegel et al. (1998). The reddening of our sample is on average E(B − V) = 0.12 ± 0.14 mag. For the nearby dwarfs, however, A K is very small and thus hard to estimate given the photometric uncertainties; hence it was set to 0 if the RJCE method yielded negative values.
With the exception of log g, stellar parameters and abundances are estimated using the spectrum synthesis code SME (Valenti & Piskunov 1996;Piskunov & Valenti 2017), which uses a Marquardt-Levenberg χ 2 -optimisation between the observed spectrum and synthetic spectra that are calculated on-thefly. As part of the GALAH+TGAS pipeline, SME version 360 is used, with marcs 1D model atmospheres (Gustafsson et al. 2008) and non-LTE-synthesis of iron from Lind et al. (2012). The chemical composition is assumed equal to the standard marcs composition, including gradual α-enhancement toward lower metallicity 1 . The pipeline is operated in the following way: 1. Stellar parameters are initialised from the analysis run used by Martell et al. (2017) if available and unflagged, otherwise the output from the reduction pipeline (Kos et al. 2017) is used and if these are flagged, we adopt generic starting values T eff = 5000 K, log g = 3.0, and [Fe/H] = −0.5. 2. Predefined 3 − 9 Å wide segments are normalised and unblended and well modelled Fe, Sc, and Ti lines within each segment are identified. Broader segments are used for the Balmer lines. The continuum shape is estimated by SME assuming a linear behaviour for each segment, and based on selected continuum points outside of the line masks. 3. Stellar parameters are iterated in two SME optimisation loops.
(a) SME parameters T eff , [Fe/H], v sin i ≡ v broad , and v rad are optimised by χ 2 minimisation using partial derivatives. (b) Whenever T eff or [Fe/H] change, log g and v mic are updated before the calculation of new model spectra and their χ 2 . We adjust log g according to Eq. 1 with isochrone-based masses M = M(T eff , log g, [Fe/H], M K S ) estimated by the Elli code, see Sect. 3.3. We adjust v mic following empirical relations estimated for GALAH 2 4. Each segment is re-normalised with a linear function while minimising the χ 2 distance for the chosen continuum points between observation and the synthetic spectrum created from the updated set of parameters (Piskunov & Valenti 2017) 5. The stellar parameters are iteratively optimised until the relative χ 2 -convergence criterium is reached.
During each optimisation iteration, a suite of synthetic spectra based on perturbed parameters and corresponding partial derivatives in χ 2 -space are computed to facilitate convergence. The parameters of the synthesis with lowest χ 2 are then either used as final parameters or as starting point of a new optimisation loop. For each synthesis, SME updates the line and continuous opacities and solves the equations of state and radiative transfer based on interpolated stellar model atmospheres (Piskunov & Valenti 2017). The optimisation has converged, when the fractional change in χ 2 is below 0.001. Non-converged optimisations after maximum 20 iteration are discarded. Figure 4 shows the final spectroscopic parameters T eff vs. log g of the final sample, colour coded by the fitted metallicities, masses, and ages from the Elli code, see Sect. 3.3.
We report the metallicity as [Fe/H] and base it on the iron scaling parameter of the best-fit model atmosphere (SME's internal parameter feh). It is mainly estimated from iron lines and hence traces to the true iron abundance, as our validation with [Fe/H] from the GBS shows.

Validation stars
To estimate the precision, we can rely on stars with multiple observations as part of the GALAH+TGAS sample: 334 stars have been observed twice and 44 stars have been observed three times.   Fig. 1 and excludes giants with T eff ≥ 5500 K and log g ≥ 3.8 dex. Stellar masses increase from the cool main sequence (∼ 0.8 M ) to the hottest turn-off stars (∼ 2.0 M ). Stellar ages decrease towards higher surface gravities on the cool main sequence and towards higher effective temperatures in the turn-off region. In contrast to this rather smooth trend, few metalpoor stars stand out with smaller stellar masses and higher stellar ages also at effective temperatures around 6000 K.
The individual differences of selected parameters are plotted in Fig. 5. Because we also use these multiple observations to assess the precision of the abundance estimates (see Sect. 3.5), we show the two element abundances Ti and Y as examples. We assume the uncertainties to be Gaussian and estimate the standard deviation of the multiple visits as a measure of precision. The resulting precisions based on the repeated observations are shown in Table 1. To estimate the accuracy, we use the GBS. These are, however, typically much brighter and closer than the survey targets and brighter than the bright magnitude limit of Gaia DR1 TGAS.  Histograms of parameter and abundance differences obtained from multiple observations of the same star. Shown are the absolute differences from two observations as well as from all three absolute difference combinations for three observations. A Gaussian distribution was fitted to the distributions (red curves). The obtained standard deviation is indicated in each panel.
Hence Hipparcos parallaxes (van Leeuwen 2007) are used additionally. New K S magnitudes are computed for GBS with 2MASS K S quality flag not equal to 'A', following the approach used by Heiter et al. (2015a) 3 . With this approach 22 GBS observations are analysed and compared to the estimates from the GALAH+TGAS pipeline, as depicted in Fig. 6. To have a statistically sufficient sample, we also include GBS giants in the analysis. We find a small bias of 51 ± 89 K in comparison to the systematic uncertainties present in both GBS and our parameters. We note, however, temperature-dependent biases of 110 ± 110 K for some stars around 5000 K. Towards higher temperatures, we also note an increasing disagreement, indicating that the temperatures of hotter stars are underestimated by our spectroscopic pipeline (−150 ± 130 K at 6600 K), a result likely to be caused by the application of 1D LTE atmospheres for hot stars (see e.g. Amarsi et al. 2018), where Balmer lines are the strongest or only contributor for the parameter estimation. For surface gravity, log g, and rotational/macroturbulence broadening, v broad , we find excellent agreement of 0.00 ± 0.05 dex and 0.9 ± 2.0 km/s respectively. The latter is computed as quadratic sum of v sin i and v mac for the GBS. For the metallicity, [Fe/H], we found a significant bias with respect to the GBS. Similar to previous studies of HERMES spectra (Martell et al. 2017;Sharma et al. 2018) we therefore shift the metallicity by +0.1 dex for our sample. The shift is chosen so that the overlap with GBS has consistent [Fe/H] in the solar regime. Two outliers for v mic can be seen to drive the bias of −0.14 ± 0.20 km/s, which we do not correct for, because the majority of the GBS sample agree well with our estimates and the two outliers are the most luminous giants, which are not representative of the final sample.
With these precision and accuracy estimates (the latter coming from the error-weighted standard deviation between GALAH and GBS estimates), we estimate the overall uncertainties of our parameters X (not mass and age, see Sect. 3.3) by summing them in quadrature to the formal covariance errors of SME (e 2 X,SME ): For element abundances, we estimate the overall uncertainties without the GBS term. In the case of log g, we replace e 2 log g,SME by the standard deviation of 10,000 Monte Carlo samples of Eq. 1. For this sampling, we use the uncertainties of e T eff ,final , the maximum likelihood masses as M with an error of 6% (based on mean mass uncertainties of an initial Elli run), e K S from 2MASS with mean uncertainties of 0.02 mag, and propagate this information to adjust BC (with typical changes below 0.07). Because Astraatmadja & Bailer-Jones (2016) only state the three quantiles, we sample two Gaussians with standard deviations estimated from the 5th and 95th distance percentile respectively. Because there are no Bayesian distance estimates for Hipparcos, we choose to sample parallaxes rather than distances D . For e A K we use the quadratically propagated uncertainties from the RJCE method (with mean uncertainties of 0.03 mag) or assume 0.05 mag for estimates based on E(B − V). We do not use Eq. 3 for age and mass, because they are estimated with the adjusted stellar parameters.

Mass and age determination
For the mass and age determination, we use the Elli code (Lin et al. 2018), employing a Bayesian implementation of fitting Dartmouth isochrones based on T eff , log g, [Fe/H], and absolute magnitude M K . M K is based on 2MASS K S , the distance estimates from Astraatmadja & Bailer-Jones (2016) and accounts for extinction A K (estimated as described in Sect. 3.1). The Dartmouth isochrones span ages from 0.25 to 15 Gyr and metallicities from −2.48 to +0.56 with α-enhancement analogous to the marcs atmosphere models 1 . Starting with a maximum likelihood mass and age estimation, MCMC samplers as part of the emcee package (Foreman-Mackey et al. 2013) are used to estimate masses and ages. Stellar ages and their uncertainties are estimated by computing the mean value and standard deviation of the posterior distribution. The stellar ages estimated with the Elli code have typical uncertainties of 1.6 Gyr (median of posterior standard deviations), which typically correspond to less than 30%, see Fig. 7. As pointed out for example by Feuillet et al. (2016), the posterior distribution does not necessary follow a Gaussian. Although this is the case for the large majority 3500 4000 4500 5000 5500 6000 6500 7000 incl. prec. and acc. of our stars, we also provide the 5th, 16th, 50th, 84th, and 95th percentiles to the community for follow-up studies. Because the results of this study did not change significantly with quality cuts for stellar ages, we do not apply them.

Binarity
The observational setup of the GALAH survey allocates one visit per observation (with exception of pilot and validation stars). Therefore, binaries or triples can usually not be identified via radial velocity changes.
Here, we use both the tSNE classifications by Traven et al. (2017), to identify obvious binaries, as well as visual inspection to identify double-line binaries which are less distinct from the tSNE classification. Within the sample, a binary fraction of 4% has been identified with high confidence from spectral peculiarities. Additionally, 338 probable photometric binaries on the main sequence are identified which show a significant deviation between spectroscopically determined log g or L bol with respect to photometrically determined ones. For these, the suspected secondary contributes significantly to the luminosity of the system without obvious features within the GALAH spectra. These stars lie above the main sequence within a colour-(absolute) magnitude diagram. We have identified the stars with photometric quantities beyond what is expected for a single star on the main sequence (shown as black dots in Fig. 8) by using a Dartmouth isochrone with the highest age (15 Gyr) and metallicity (+0.56 dex). We note that some of these stars show colour excesses. While these might have been mis-identified as binaries, they are definitely peculiar objects (e.g. pre-main-sequence stars), for which the pipeline is not adjusted and have subsequently been neglected. We want to stress again, that identified binaries have been excluded from the cleaned sample.

Abundance determination
With the stellar parameters estimated in Section 3.1, elemental abundances were calculated in the following way: 1. Predefined segments of the spectrum are normalised and the element lines chosen with two criteria. First, the lines have to have a certain depth, i.e., their absorption has to be significant. We use the internal SME parameter depth to assess this, see Piskunov & Valenti (2017). 2. The lines have to be unblended. This is tested by computing a synthetic spectrum of the segment with all lines and one only with the lines of the specific element. The χ 2 difference between the synthetic spectra for each point in the line mask has to be lower than 0.0005 or 0.01 (the latter for blended but The colour index is J − K S from 2MASS photometry and absolute magnitude for K S , inferred from 2MASS as well as distances D and extinction A K . A Dartmouth isochrone with age (15 Gyr) and metallicity (+0.56) is shown as white curve. This is used to identify 338 dwarfs with photometry outside of the expected range (above the white curve) for cool single main sequence stars (M K S > 2 mag), here shown in black. The identified stars are all nearby and their reddening is negligible, especially in the infrared. We note that for some stars, possibly mis-identified as binaries, the photometry indicates colour excesses or a pre-main-sequence stage, which is still an important reason to eliminate them from the subsequent analysis, as the pipeline is not adjusted for these stars.
indispensable lines), otherwise the specific point is neglected for the final abundance estimation. 3. The abundance for the measured element is optimised using up to 20 loops with the unblended line masks.
The selection of lines used for parameter and abundance analysis and their atomic data is a continuation of the work presented by Heiter et al. (2015b) and Heiter et al. (submitted). The complete linelist will be presented in Buder et al. (2018).
Abundances are estimated assuming LTE, with the exception of Li, O, Al, and Fe, for which we use corrections by Lind et al.
Nr. Stars 500 1000 1500 2000 2500 Fig. 9. The left three panels show the distribution of angular momentum L z , sorted by relative uncertainty and depict three close in views groups of with 75 stars with mean L z uncertainties of first 0.5%, second 5%, and third 30% in order to demonstrate the precision reached for different parallax qualities. Red is mean L z for each star and blue their 1σ area. A white dashed line indicated the Solar angular momentum. Average parallax precisions are indicated in the top of each panel. The best precision on parallaxes also lead to the most precise L z . The values of the least precise momenta (third panel) are significantly lower than those of the Sun, even when taking the standard deviation of the angular momentum estimates into account. Due to the selection of stars and the density structure of the disk, these stars are statistically further away and are expected to be at larger Galactic heights and closer to the Galactic centre. We note that their angular momenta are also different from the majority of stars, which have a Sun-like angular momenta, as shown in the fourth panel. For a discussion of the angular momenta of the stars in the chemodynamical context see Sect. 5.
Solar abundances are estimated based on a twilight flat in order to estimate the difference to the solar composition by Grevesse et al. (2007, G07). This difference for each element X, i.e. A(X) − A(X) G07 , is then subtracted from the element abundance of the stars of the sample.

Kinematic parameters
For our target stars, the space velocities U, V, and W were calculated using the galpy code by Bovy (2015), assuming (U , V , W ) = (9.58, 10.52, 7.01) km/s (Tian et al. 2015) relative to the local standard of rest.
We estimate kinematic probabilities of our sample stars to belong to the thin disk (D), thick disk (T D), and halo (H) following the approach by Bensby et al. (2014) (see their Appendix A) with adjusted solar velocities.
To estimate the Galactocentric coordinates and velocities as well as the action-angle coordinates of the sample, we use galpy. We choose the axisymmetric MWPotential2014 potential with a focal length of δ = 0.45 for the confocal coordinate system and the galpy length and velocity units 8 kpc and 220 km/s respectively. We place the Sun at a Galactic radius of 8 kpc and 25 pc above the Galactic plane. To speed up computations, we use the actionAngleStaeckel method. We estimate mean values and standard deviations of the action-angles per star from 1000 Monte Carlo samples of the 6D kinematical space randomly drawn within the uncertainties. We neglected the uncertainties of the 2D positions and estimated the standard deviation of the distances from the 5th and 95th percentiles given by Astraatmadja & Bailer-Jones (2016).
As shown in Fig. 9, the distance uncertainties are the dominant source of the action uncertainties. While for excellent parallaxes (left panel), the scatter in the action estimates is negligible, it becomes noticeable for parallaxes with uncertainties around 18%. For parallax uncertainties above 24%, the action uncertainties increase to as high as 31%. From the samples depicted in Fig. 9, one can see that these large uncertainties are particularly common for stars with low angular momentum. Because of the GALAH selection (observing in the Southern hemisphere and leaving out the Galactic plane) as well as the density structure of the disk with more stars towards the Galactic centre, we expect stars with larger distances (and hence larger distance uncertainties) to be situated at larger Galactic heights and smaller Galactic radii than the Sun. The right panels in Fig. 9 confirm this expectation. Stars with angular momenta comparable with the solar value have usually precisely estimated actions. The latter stars are also the majority of stars in the sample, as the histogram in the right panel shows.

Results
In Sect. 4.1, we describe the stellar age distribution, before presenting abundance and age trends in Sect. 4.2 and the kinematics of the sample in Sect. 4.3.
We note that the vast majority of the dwarfs from the GALAH+TGAS sample are more metal-rich than −0.5 dex, as seen in the metallicity distribution function in Fig. 10. These stars have no intrinsic selection bias in metallicity or kinematics.
From 10,000 Monte Carlo samples, we find the parameters of the metallicity distribution to be [Fe/H] = −0.04, σ [Fe/H] = 0.26, skewness = −0.667 ± 0.029, kurtosis 4 = −0.21 ± 0.23. The mean of our metallicity distribution is slightly lower but consistent within the uncertainties to the one estimated by Hayden et al. (2015) using APOGEE data for the same (solar) Galactic zone 5 . The APOGEE distribution also shows a narrower standard deviation (0.2 dex) around a mean value of +0.01 dex and  is less skewed (−0.53 ± 0.04) but more extended towards the metal-rich and metal-poor tail of the distribution (with a kurtosis of 0.86 ± 0.26). The kurtosis, a measure for the sharpness of the peak, indicates that the APOGEE distribution has a sharper peak than the GALAH distribution. The skewness indicates that the GALAH sample contains in general also relatively more metal-poor stars compared to the APOGEE sample. This is possibly caused by the different selection functions of the two surveys. GALAH avoids the Galactic plane (|b| ≤ 10 deg), whereas APOGEE targets the plane where we expect relatively more stars of the low-α-sequence that are more metal-rich than [Fe/H] = −0.7.

The age distribution
The age distribution of the GALAH+TGAS sample is shown in Fig. 11. It peaks between 3 to 3.5 Gyr, which is at an older age than estimated by the studies of Casagrande et al. (2011) andSilva Aguirre et al. (2018) who both placed the peak at approximately 2 Gyr. While this might be partially explained by a combination of both selection function, and target selection effects, we note that the exclusion of hot stars with effective temperatures above 6900 K in our sample, see Sect. 2, affects primarily stars with ages below the peak of the histogram. However, these hot stars have an average maximum likelihood age of 1.5 ± 0.8 Gyr and the location of the age peak does not change when including them.

The age-[α/Fe]-[Fe/H] distributions
We have detected abundances for up to 20 elements, which are presented in Sect. A. For an extended overview of abundance trends for the elements detectable across the whole GALAH range, we refer the reader to Buder et al. (2018). For this study, we focus on the α-elements and iron as well as their correlations with stellar age. The combination of these three parameters is shown in Fig. 12.
The abundance patterns of α-elements in the Galactic disks are expected to follow roughly a similar pattern according to the  Fig. 11. Distribution of stellar ages. The distribution peaks around 3 Gyr and decreases towards higher ages. We stress that the exclusion of stars with effective temperatures above 6900 K leads to fewer stars in the clean sample, with ages below 2 Gyr. The peak of the distribution is however not affected by this selection. stellar enrichment history by supernovae type Ia and II (see e.g. Gilmore et al. 1989). While both types of supernovae produce a variety of elements, there is a significant difference in the yields of iron and α-elements and the time in the Galactic evolution, when they each contribute to the chemical enrichment. Early in the chemical evolution of the Galaxy, SN Type II dominate the production of metals and large quantities of α-elements are then produced (e.g. Nomoto et al. 2013 The combined α-element abundance is estimated for 99% of our stars. For each of these stars, at least one α-process element has been detected and all significant measurements have been combined with their respective uncertainties as weight. Mg, Si, and Ti are the most precisely measured elements and have the highest weight. Hence, we note that the [α/Fe]-ratio, as defined here, is in practice very similar to the previously used error-weighted combination of Mg, Si, and Ti for GALAH DR1 (Martell et al. 2017) and for the study by Duong et al. (2018). We see overall good agreement in the [α/Fe] pattern with the stars in the solar vicinity analysed by the APOGEE survey (Hayden et al. 2015), i.e. predominantly solar ratios for −0.7 and +0.5 dex and fewer stars with increasing α-enhancement towards lower metallicity. We discuss this bimodality further in Sect. 4.2 by inspecting the quantitative distribution of [α/Fe] in several metallicity bins, see Fig. 13. We stress that there is no unambiguous or universal definition of α-enhancement, but studies estimate and define this parameter differently, which complicates comparison. Hereinafter we use an average, weighted with the inverse of the errors, of the four α-process elements (Mg, Si, Ca, and Ti) when we refer to [α/Fe] and α-enhancement. We note that because our definition of α-enhancement is driven by Ti as the most precisely determined element, our values are comparable with the study by Bensby et al. (2014) based on Ti. Fuhrmann (2011) used only Mg as tracer of the α-process ratio. All different definitions hence induce possible systematic trends.   Fig. 12 (a) to (c).

The [α/Fe]-[Fe/H] distribution is shown in
The pattern of our study is agrees very strongly with the results found by Ness et al. (2016) for APOGEE (see their Fig. 8) and Ho et al. (2017a) for LAMOST (see their Fig. 5), all three showing high age for the high-α sequence and younger ages for stars on the low-α sequence. We note that stars with larger ages usually have larger absolute age uncertainties. In contrast to the expected rather monotonic trend between α-enhancement Article number, page 11 of 30 A&A proofs: manuscript no. Buder_GnT_v2 and stellar age (especially at constant metallicity), we note that around −0.4 < [Fe/H] < 0, young and fast rotating stars are dominating the interim-[α/Fe] regime. For hotter stars with v sin i > 15 km/s, the estimated iron abundances A(Fe) are typically lower than the one of slow rotators. While this could be a trend introduced by the analysis approach that depends on sufficiently deep metal lines, another possibility is an actual correlation between [Fe/H] and rotation. An analysis of this correlation is complex and beyond the scope of this paper. When we neglect such stars (10% of the sample), the trend of stellar age and [α/Fe] is monotonic.
For the high-α metal-rich regime, a mix of different ages is noticeable, with an age spread up to 4 Gyr. We take a closer look at this region in Sect. 5.
To assess how distinct the two α-enhancement sequences are at different metallicities, we plot the histograms for five 0.15 dexwide metallicity bins in Fig. 13. By eye, two clear peaks can only be identified for the three lower metallicities with decreasing separation. However, the fit of two Gaussians recovers the two peaks for all five distributions. For the most metal-poor bin, the α-enhanced stars are more numerous with an enhancement of 0.25±0.03 dex, compared to the low-α stars at 0.13±0.06 dex. We note that even the low-α stars are slightly enhanced at these metallicities. At higher metallicities, the mean enhancement of the low-α sequence decreases gradually to become solar at solar metallicity.
The enhancement of the high-α sequence decreases more steeply down to 0.04 ± 0.05 dex at solar metallicity. The peaks of the two (forced) sequences are thus consistent within one sigma (indistinguishable) at solar metallicity. We stress that in our fit we forced two Gaussian distributions and the actual distribution looks like a positively skewed Gaussian distribution. This means that an assignment to the high-or low-α sequence based on a given [α/Fe] threshold Adibekyan et al. (e.g. 2012) is significantly less accurate or meaningful than in the metal-poor regime Duong et al. (see also 2018).
The widths of the Gaussian fits to the high and low-α sequence are of order 0.02 − 0.08 dex for [α/Fe] and similar to our measurement uncertainties and we note that the separation between the two sequences in the metal-poor regime is larger than this.
The [α/Fe]-age distribution is shown in Fig. 12 (d) to (f). The main findings from these panels are: 1. The mean α-enhancement stays rather constant at [α/Fe] ∼ 0.05 up until 8 Gyr and then increases with stellar age. A comparison of our relation with the one found for the stars analysed by Bensby et al. (2014), see left panel of Fig. 14, shows that the observed relations agree within their measurement uncertainties. 2. We find 6% of the stars below 8 Gyr with [α/Fe] > 0.125.
The corresponding fraction for stars older than 11 Gyr is 60%. This indicates a small jump around 8 to 10 Gyr (from mean [α/Fe] ∼ 0.05 to 0.09) between high and low-α enhancement, as also found by Haywood et al. (2013). At around 8 to 12 Gyr, we note a large range of metallicity in these coeval stars. 3. We find 67 stars among the young ones (i.e. 4978 stars with < 6 Gyr), that are significantly α-enhanced ([α/Fe] > 0.13) with normal rotational velocities (and another 59 with v sin i > 15 km/s). With ∼ 0.9% of our sample (∼ 1.8% when including the fast rotators), their ratio is in agreement with the sample analysed by Martig et al. (2015), who found 14 out of 1639 stars to be α-rich and similar to most of the ratios of other samples listed by Chiappini et al. (2015). Looking only at the young stars (< 6 Gyr) our ratio of 1.4 − 2.5% is however smaller than the one found by Martig et al. (2015) of 5.8%, pointing towards a different age distribution of the two different samples (containing either only giants or main sequences/turn-off stars). The 59 (47%) of the young α-rich stars in our sample with increased broadening are all hotter than 6000 K. We want to stress that for such stars the broadening in addition to the decreasing line strengths due to the atmosphere structure make the parameter estimation more uncertain than for the rest of our sample. 4. Among the old stars (> 11 Gyr), we find that a significant fraction of the sample (30%), are low-α stars ([α/Fe] < 0.125), in contradiction to Haywood et al. (2013). These stars are primarily cool main sequence or subgiant stars with metallicities above −0.5 dex, which causes their ages to have larger error bars. Silva Aguirre et al. (2018) also found such stars among APOGEE giants from the Kepler sample of stars, using asteroseismology to determine precise ages (see their Fig. 10). 5. Haywood et al. (2013) claimed a rather tight correlation between age and α-enhancement for the old high-α stars. However, Silva Aguirre et al. (2018) did not see evidence for such a tight relation. Our sample implies a tight trend, but is limited by the small number of these stars and we can not draw strong conclusions regarding the dispersion of chemistry and age. Fig. 12 (g) to (i) and shows:

The [Fe/H]-age distribution is depicted in
1. With increasing stellar age, the mean metallicity, indicated with a red line in panel (g), decreases steadily and non-linear from 0.04 at 1 Gyr to −0.56 at 13 Gyr respectively, but also with increasing dispersion. The recent study by Feuillet et al. (2018) has found that in their sample both the lowest [Fe/H] and highest [Fe/H] stars are older than the solar abundance stars and the iron abundance is hence less useful than [α/Fe] to predict age. 2. The dispersion increases steadily with stellar age from 0.17 at 2 Gyr to 0.35 at 13 Gyr, respectively, as the steadily increasing red dashed line (calculated in steps of 1.5 Gyr similar to the median age uncertainties and each containing at least 50 stars), which traces the dispersion of [Fe/H], shows in Fig. 12 (g). 3. Between ages of 0 and 8 Gyr, we see a rather flat mean metallicity with a spread of 0.5 dex around solar metallicity. In this age range, stars with lower metallicities (< −0.25) show an increase of α-enhancement of up to around 0.1 dex. 4. Above ages of 8 Gyr, metal-poor stars with [Fe/H] < −0.25 exhibit a decreasing trend of metallicity with increasing stellar age. We can not confirm the tight age-metallicity distribution for the oldest stars due to the small sample of these stars, although we note indications of a tight overdensity. Few old stars around solar metallicity stand out, as also found in previous studies by Casagrande et al. (2011, (see their Fig. 16)). These stars cause an increased dispersion, which is in agreement with the continuously increasing dispersion also seen at lower ages (see Fig. 12 (g)). Similar results were found by Haywood (2008b, see their Fig. 1), who interpreted this trend as an observational signature of radial migration. We will follow this up, when including kinematics, in the following section. We find that the old stars of the sample are more α-enhanced. The oldest stars (above 11 Gyr) are most metal-poor and αenhanced (around 0.25 dex at metallicity [Fe/H]∼ −0.5), while slightly younger stars (still above 8 Gyr) are less α-enhanced and more metal-rich (around 0.15 dex at solar metallicity). Most of the stars of the sample exhibit slightly increased or solar [α/Fe] and are on average younger than 6 to 8 Gyr, but the sample also contains a minority of old stars (around 8 Gyr) with low metallicity (around −0.6 dex) and only slight α-enhancement. We assess and discuss these results in more detail in Sect. 5

Kinematics
To get an overview of the kinematical content of the GALAH+TGAS overlap, we first examine the 3D velocities. We then use these velocities to assign membership probabilities to different Galactic components.
From the Toomre diagram in Fig. 15, we can deduce that most stars belong to the disk, because their total velocities are lower than 180 km/s, which is typically adopted as a good limit of halo kinematics (Venn et al. 2004;Nissen & Schuster 2010). This was expected from the target selection (De . Most of the sample shows an azimuthal relative velocity close to the local standard of rest, i.e. |V| < 50 km/s and also U and W are close to the local standard of rest, indicating a solarlike motion in the thin disk. However, there are stars that also show large deviations from the thin disk kinematics, with total velocities larger than 100 km/s relative to the local standard of rest. In previous studies, such stars have been identified as thick disk stars, based on a hard limit of 70 km/s (Fuhrmann 2004  Nr.

500
Nr. Adopting hard limits to separate populations is however not appropriate when it comes to kinematics, as both thin and thick disc populations show significant dispersions in their characteristic space velocities (Nordström et al. 2004). Several more sophisticated approaches have been implemented (see e.g. Reddy et al. 2006;Ruchti et al. 2011;Bensby et al. 2014). To begin with, we follow the approach by Bensby et al. (2014) to estimate the probability of each star to belong to one of the Milky Way components thin disk (D), thick disk (TD) or halo (H) by their kinematical information including the population velocity dispersions. The probability is influenced by the velocity distribution (assumed to be Gaussian), rotation velocities, as well as the expected ratio of stars among the components (Bensby et al. 2003). Similar to Bensby et al. (2014) we subsequently use the ratios of the membership probabilities.
For each component, it is possible to separate the sample into most likely thick disk stars (T D/D > 10) and most likely thin disk stars (T D/D < 0.1), as shown in Fig. 16. The 12 stars, that fit the kinematics of the halo best, are marked as big black circles and contribute 0.16 % to the sample. Almost all of these stars show a total space velocity of more than 180 km/s and would also be identified as halo stars with the simplified velocity criteria.
Most of the GALAH+TGAS dwarfs are most likely to be affiliated with the thin disk, while only a small fraction belongs to either thick disk or halo. For a very large fraction of the stars however, the probability ratio is indecisive.
While the analysis of kinematics with the classical Toomre diagram is a powerful tool for a small volume, the approximations (assuming similar positions and velocities) are less appropriate for larger volumes. We therefore include another way to interpret the kinematical information by calculating action-angle  coordinates and characterise orbits with integrals of motion as proposed by McMillan & Binney (2008). Contrary to U, V, and W, the three orbit labels, the actions J R , J Φ = L z , and J z , allow us to quantify and compare orbits of stars independent of their position relative to the Sun or the local standard of rest. The distribution of the three actions is shown in Fig. 17 and shows that most of the stars have similar orbits to the Sun, meaning low eccentricities and radial actions J R , low vertical oscillations and vertical actions J z , and azimuthal actions similar to the one of the Sun (with L Z, = R · v Φ, = 8 kpc · 220 km/s = 1760 kpc km/s). However there are several stars with J R > 75 kpc km/s on more eccentric orbits, which manifests in mean stellar radii either significantly closer to (L z L z, ) or further away from (L z L z, ) the Galactic center. We will follow this up in Sect. 5, when analysing angular momenta of the individual stars with different ages and chemistry, see Fig. 18.

Discussion of chemodynamics of the Galactic disks
The previous discussion highlighted that we report a significant change of chemistry around the stellar age of 8 to 10 Gyr. Haywood et al. (2013) and Bensby et al. (2014) used this to establish population assignments based on ages alone or ages and chemistry combined. Haywood et al. (2013)   In Fig. 18, we plot age bins; we can therefore focus on specific ages rather than age sequences. For each of the 1 Gyr bins, we show the [α/Fe]-[Fe/H] plane and colour the stars by their angular momentum, which is a function of Galactic radii. . Our main findings are: 1. The oldest stars (> 13 Gyr) are mostly α-enhanced and metal-poor and have not yet been significantly enriched by SN Ia because they were born before the delay time of SN Ia was reached. Their spread in metallicities can be explained by different SN II masses and and frequencies, as well as gas mixing efficiencies, in the progenitor clouds. Their angular momentum shows that stars older than 13 Gyr are usually located closer to the Galactic centre with mean L z = 1060 ± 600 km/s, which corresponds to an average Galactocentric distance of 4.8 ± 2.7 kpc. 2. We note however, four stars with roughly solar metallicities older than 13 Gyr, which could be unidentified binaries because of their proximity to the main sequence binary sequence. Such stars have also been found by Casagrande et al. (2011, see their  Their angular momenta and actions point towards solar-like orbits for half of them or eccentric orbits with significantly lower angular momenta than the Sun for the other half. The presence of these stars could be explained via chemical evolution, radial migration, but also influence of the bar. The recent study by Spitoni et al. (2018) showed, that a revision of the "two-infall" model can explain the presence of old stars with −0.5 < [Fe/H] < 0.25 and [α/Fe] < 0.05. Such stars have however also been found in the inner disk (Hayden et al. 2015) and even the bulge (see e.g. Bensby et al. 2017). The latter found even a high fraction of one-third of low-α stars among the old stars. However the models by Spitoni et al. (2018) are not able to fully recover e.g. the age distribution and additional dynamic process are needed to explain the observed data. If we assume radial migration of such stars, we expect also them to arrive in our solar neighborhood. Models for radial migration (see e.g. Frankel et al. 2018) however usually aim to model the secular evolution of stars and hence currently focus on stars younger than 8 Gyr. We do however not see a reason why the oldest stars, which clearly exist in the inner Galaxy (Hayden et al. 2015;Bensby et al. 2017) should not have migrated in a similar manor as the younger Fiducial lines indicate the low-α and high-α sequences to guide the eye. This plot shows, that stars with high-α enhancement have generally lower angular momentum, while stars in the low-α sequence, have mostly solar momentum. We note however that the metal-rich low-α stars have in general lower angular momenta than metal-poor stars with comparable α-enhancement. For a detailed discussion of the panels see Sect. 5. Article number, page 16 of 30 stars within the lookback time of the migration models. Further analyses beyond the scope of this paper will hopefully help to pin down the more likely reason for the presence of the old stars with solar [Fe/H] and [α/Fe]. Another explanation that can be tested with more extended data sets is the radial migration due to the influence of the Milky Way bar (see e.g. Grenon 1999;Minchev & Famaey 2010). 3. Stars between 12 and 13 Gyr exhibit more iron than the oldest stars although at similar [α/Fe]. In the metal-poor regime ([Fe/H] < 0), the high-α stars have lower angular momentum than the Sun ( L z = 1267 ± 389 kpc km/s with mean L z uncertainties of 63 kpc km/s), while more metalrich stars ([Fe/H] > 0) have angular momenta ( L z = 1513± 262 kpc km/s with mean L z uncertainties of 46 kpc km/s closer to the solar one (L z, = 1760 kpc km/s). 4. Below 12 Gyr, all stars except three outliers have metallicities above −1.0 dex. Between 12 and 9 Gyr, the stars on the high-α sequence become gradually less α enhanced and show increasing metallicities, hence indicating a continuous evolution of high-alpha stars along this sequence. This is consistent with the increasing enrichment of the ISM by SN Ia with a delay time distribution (see e.g. Maoz et al. 2012), producing significantly more iron than SN II. Their angular momenta are on average still significantly lower than the solar one. This indicates a continuous evolution of high-α stars along this sequence. We want to stress that these stars are still slightly α-enhanced, even at solar metallicities. Below 9 Gyr, only a few stars are on the high-α sequence and almost none below 7 Gyr. 5. Similar to previous studies (Lee et al. 2011;Haywood et al. 2013) we find a gradual increase of angular momentum and rotational velocity with metallicity among the alphaenhanced metal-poor stars, which are is also correlated with age above 10 Gyr. For stars below 10 Gyr (but even more pronounced for stars below 8 Gyr), the angular momentum decreases with metallicity, as shown in Fig. 19. 6. Around 10 Gyr and at younger times, stars at the metal-poor end of the low-α sequence appear. This finding is consistent with previous results by Haywood et al. (2013, see their Fig. 8). The angular momenta of many of them indicate an origin at larger Galactic radii, i.e. they are only visitors to the solar neighbourhood, in agreement with findings by Haywood (2008b) and Bovy et al. (2012). 7. Between 3 and 9 Gyr, the full range of metallicities from −0.7 up to 0.5 at the low-α sequence is covered. The stars at the low-α metal-poor end show on average significantly larger angular momenta than the Sun, an opposite trend with respect to the high-α stars. The angular momenta of metalrich stars (1640 ± 180 kpc km/s) are on average lower than the solar one. 8. At younger times than 3 Gyr, the spread in metallicities decreases and stars younger than 2 Gyr only cover metallicities between −0.3 and 0.3. Note that our cut for hot stars, see Sect.
2) has cut out most of the stars with ages below 2 Gyr, see Sect. 4.1. Casagrande et al. (2011) have found, however, that nearby stars younger than 1 Gyr also only cover the range of −0.3 < [Fe/H] < 0.2 (see their Fig. 16).
Similar conclusions of a gradual chemical enrichment of the high-α sequence with time can be drawn from the data by Hayden et al. (2017, see their Fig. 3), as well as from the analysis of Auriga simulations (Grand et al. 2018).
Finally, we revisit the α-enhancement of the sample, but dissect the sample by the kinematic probability of belonging to the thick or thin disk. These probabilities were estimated in Sect. 3.6 and we use them to recreate Fig. 19 from Bensby et al. (2014) with the GALAH+TGAS sample. From this dissection in Fig. 20 we conclude: 1. Stars that are kinematically > 10 times more likely to be part of the thick disc population (top panel), are mostly on the high-α sequence. We note that in contrast to Bensby et al. (2014), these stars are not all metal-poor, but also cover the high-α metal-rich regime. These stars are however almost exclusively older than 8 Gyr. The sample from Bensby et al. (2014) is expected to cover relatively more metal-poor stars, because they selected their sample specifically with the aim to trace the metal-poor limit of the thin disk, the metal-rich limit of the thick disk, the metal-poor limit of the thick disk. 2. Stars that are 2 − 10 times as likely to belong to the thick disk (2 < T D/D < 10) and those with inconclusive kinematics (0.5 < T D/D < 2) cover both the high-and low-α metalrich regime as well as a larger range of ages. 3. At the lowest probabilities of belonging to the thick disk (T D/D < 0.5), the stars are on average young (5.0 ± 2.2 Gyr) and cover the low-α sequence. However, there are still noticeable amounts of stars in the high-α metal-rich regime ([Fe/H] > −0.1 and above the dashed line, including stars older than 8 Gyr. This implies that there is no distinct kinematical separation of populations in the metal-rich regime and a significant overlap in kinematic properties is present.

Conclusions
The combination of spectroscopic data from the high-resolution GALAH survey and the astrometric data from the Tycho-Gaia Astrometric Solution (TGAS) spans a high-dimensional space of chemodynamical information. In this study, we have analysed 7066 dwarf and turn-off stars of this overlap. Our main results are summarised as follows: 6.1. Abundance and age trends 1. We show that our parameter and abundance estimates measured with the GALAH pipeline are accurate and precise.  Fig. 19 of Bensby et al. (2014), the abundance plateau of the thick disk and the decrease for thin disk have been plotted in black to guide the eye. While stars with highly thick-disk like kinematics (two upper panels) follow the high-α sequence, numerous stars with thin-disk like kinematics are also seen in the high-α metal-rich regime. However, most stars with thin disk like kinematics (T D D < 0.5 two lower panels) follow the low-α sequence. See Sect. 5 for further discussion.
2. Our selected stars are all within the solar vicinity and consist mostly of young stars with metallicities from −0.7 dex up to +0.5 dex as shown by the age and metallicity distributions. 3. We report stellar parameters (including iron in non-LTE), 18 element abundances in LTE, an error-weighted α-process element abundance in LTE, and Li as well as O abundances in non-LTE. We show that non-LTE corrections for the O triplet are a vital requirement for accurate abundance estimation of stars in different evolutionary stages. 4. We show in Fig. A.1, that the abundance trends estimated for the GALAH+TGAS sample agree well with previous studies with higher spectral quality. Due to the larger number statistics and the selection of our stars, we are able to assess abundance trends of the low-α sequence in more detail than previous studies. 5. Among the studied elements, we find the expected similarities of abundance trends among the α-elements Mg, Si, Ca, and Ti. O, however, shows a significantly steeper decrease with metallicity and no flattening at super-solar metallicities and was hence not treated as an α-element. The odd-Z elements Na and Al show similar increasing abundances towards the super-solar metallicities. We note that these trends also agree with those observed for the iron-peak elements Ni, Cu, and Zn.
6. We find significant trends of abundances with stellar age for αand s-process elements, see Fig. A.2, similar to previous studies of solar twins (Nissen 2015;Spina et al. 2016) and dwarfs in general (Bensby et al. 2014).

The age-[α/Fe]-[Fe/H] relationship
1. We recover the same pattern of the [α/Fe]-[Fe/H] in combination with age as previous studies (e.g. Ness et al. 2016;Ho et al. 2017a), namely that the high-αsequence is mainly consisting of old stars, while the low-α sequence covers ages usually below 8 Gyr. We note however, that at the the low-α metal-rich and low-α metal-poor regime, we also find stars with ages between 8 to 10 Gyr. 2. When using age and chemistry together, the high-α metalrich stars can be well explained to be part of a population formed from the same material as the canonical (lower metallicity) thick disk, although significantly enriched by both SN II and Ia. This is in contradiction to the claim by Adibekyan et al. (2012), but backed up by stellar kinematics. The angular momenta of the high-α metal-poor and high-α metal-rich stars are both significantly and consistently lower than the solar value. We conclude that these stars belong to the high-α sequence rather than representing a distinct population. 3. At solar and super-solar metallicity, both old and young stars as well as slightly α-enhanced and solar [α/Fe] stars are observed. We have shown in Fig. 13, that in this regime, the two distinct α-sequences (as seen for [Fe/H] < 0) become indistinguishable in [α/Fe] vs. [Fe/H]. We conclude, similarly to Haywood et al. (2013), that stellar age is a better identifier and discriminator of formation origin than metallicity, especially at high metallicity. Stellar age when available is more fundamental than an arbitrary cut in [α/Fe] for the purpose of understanding Milky Way populations. 4. The decreasing spread of metallicities for stars younger than 4 Gyr is strong support for radial migration. Radial migration predicts that stars will change their orbits after they have been born. The panel of Fig. 18 with stars below 2 Gyr shows a narrower range in metallicity compared to the panel for stars with ages of 3 − 4 Gyr. With increasing age, stars are observed at lower and higher metallicity ranges. Yet, there are no stars with metallicities at −0.5 and 0.5 dex for the youngest ages, although the GALAH selection is not biased against these stars. At increasingly younger ages we preferentially see stars formed in the solar neighbourhood, whereas older stars did not have time to migrate to this location. Therefore, we only see stars formed in the solar neighborhood. This observational support for radial migration in the thin disk has also been found in other studies Haywood (2008a); Feuillet et al. (2016, see e.g.) and tested by comparison with models including radial migration by Feuillet et al. (2018). We stress however, that we excluded stars with T eff > 6900 K from this study, which are mostly younger than 3 Gyr. While we advice caution that our result might be biased, Casagrande et al. (2011)  2. Stellar age does not change after the birth of a star. However, kinematic properties of stars from different populations overlap and might also change. Thus, age should be a more reliable definition of a population. 3. Independent of the population assignment of the high-αmetal-rich stars (see Sect. 6.2), their kinematics (e.g. their lower angular momentum compared to the Sun) indicates, that they have migrated towards us via blurring, meaning they are on eccentric orbits with mean radii closer to the Galactic center. 4. Around ages of 8 to 10, Gyr, stars at the metal-poor end of the thin disk are identified. When including kinematical information, we find a large number of these stars to be from outside the solar neighborhood (on eccentric orbits with mean radii further out). At the metal-rich end of our sample, the kinematical information points to a significant number of stars from the inner radii on eccentric orbits. This result is in agreement with the seminal paper by Edvardsson et al. (1993, see their Fig. 22) and confirms blurring as reason for a broad age-metallicity relation at different radii.

Future investigations
The follow-up of our study with more detailed studies focusing of certain aspects of the high-dimensional chemodynamical space is promising. Anders et al. (2018) have shown that the chemical space can be dissected by using t-distributed stochastic neighbour embedding. For the high-α metal-rich regime they find that these stars are different from the classical thick disk. They find however, that the stars in this regime are spread within the t-SNE distribution, indicating an evolution within this regime (see their Fig. 1). The application of their approach onto our data is promising but beyond the scope of this paper. In this study, we have shown that non-LTE corrections play a crucial role for several elements. We suggest a detailed further investigation of 3D and non-LTE effects for other elements. We are working on extending the non-LTE implementation and will apply them to 7 elements in GALAH DR2 .
A quantitative study of the correlation of element abundances and stellar ages for the stars of our sample could explore the influence of other parameters on the [Y/Mg]-age correlation beyond the solar twins (Nissen 2015;Spina et al. 2016), to test the hypothesis by Feltzing et al. (2017), that the [Y/Mg]-age relation is unique to solar analogues. Such analyses can be performed on all different elements, especially when using the second public Data Release from the GALAH survey ) and the second Gaia Data Release. With this large set of abundances, ages, and kinematics, it will also be possible to identify clumps in chemodynamical space or stellar streams (see e.g. Quillen et al. 2018). A central point of Galactic Archaeology in the future will also be the improvement of stellar age estimation.
ship. Parts of this research were conducted by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. We thank A. M. Amarsi for providing the grid of abundance corrections for oxygen. We thank U. Heiter for providing adjusted K S band magnitudes for the Gaia FGK benchmark stars. We thank H.-W. Rix, M. Fouesneau, A. M. Amarsi, M. Bergemann, and A. Quillen for useful discussions and comments. We would like to thank the anonymous referee for carefully reading our manuscript and for giving such constructive comments which substantially helped improving the quality of the paper. The computations were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC) at UPPMAX under project 2015/1-309 and 2016/1-400. This project was developed in part at the 2016 NYC Gaia Sprint, hosted by the Center for Computational Astrophysics at the Simons Foundation in New York City, and the 2017 Heidelberg Gaia Sprint, hosted by the Max-Planck-Institut für Astronomie, Heidelberg. This work has made use of data from the European Space Agency (ESA) mission Gaia (http://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, http://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France. The original description of the VizieR service was published in A&AS 143, 23. This research mad use of the TOPCAT tool, described in Taylor (2005). This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.

Appendix A: Abundance trends
In total, the GALAH survey wavelength range includes detectable atomic absorption lines of up to 30 elements of FKG stars. A subset of 20 of these can be measured in the dwarf and turn-off star spectra, which we examine in this work: Li, C, O, α-process elements (Mg, Si, Ca, and Ti), light odd-Z elements (Na, Al, K), iron-peak elements (Sc, V, Cr, Mn, Co, Ni, Cu, and Zn), and s-process neutron capture elements (Y and Ba). We list the mean precision (estimated from repeated observations), accuracy (inferred from the uncertainties of the oscillator strengths), measured fraction of the clean sample in percent, and numbers of measured lines in Table A.1. In Figures A.1 and A.2, we show our measured element abundances as a function of metallicity. We present the density distribution for these in Fig. A.1, coloured by stellar counts per bin. In Fig. A.2, we plot the same distribution, but colour each bin by the median age of its stars. When available, we include results from previous studies (Bensby et al. 2014;Battistini & Bensby 2015;Nissen et al. 2014;Zhao et al. 2016;Delgado Mena et al. 2017). These studies use spectra of higher quality (higher resolution and S /N), but are much smaller in sample size. Similar to the GALAH survey, these literature studies include dwarfs of both the low-and high-α sequences of the disk and span a large metallicity range (−2.8 < [Fe/H] < 0.4). We note that the selection function for these studies is different to that of GALAH and they contain much higher relative numbers of metal-poor stars.
For a quantitative discussion of age trends we refer the reader to other studies, e.g. for [Y/Mg] (Nissen 2015;Spina et al. 2016), [C/N] (Masseron & Gilmore 2015;Martig et al. 2016a;Ness et al. 2016) or the study of 17 abundance-age trends for APOGEE by Feuillet et al. (2018).

Appendix A.1: Lithium
Li is measured using 1D non-LTE corrections by Lind et al. (2009). In our sample, we can only detect Li in stars with a relatively large absolute Li abundance, of A(Li) > 2.0. In the metalpoor regime, the warm stars are situated on the Spite plateau (Spite & Spite 1982), around A(Li) = 2.3, as expected. Towards higher metallicities, A(Li) is mostly measured between 2.0 and 3.3. The latter value is close to the meteoritic A(Li) of 3.26±0.05 ). Li is expected to be a good tracer of the evolution of the star, because the initial composition is depleted by the proton bombardment processes at temperatures higher than 2.5 · 10 6 K (Pinsonneault 1997). From the respective panel in Fig. A.5, we see a strong correlation of higher Li for higher effective temperatures. The lines of Li in colder and evolved stars are below the 3σ detection limit. The strong anti-correlation between Li and effective temperature (or mass) is well known and previously depicted for thin/thick disk stars e.g. by Ramírez et al. (2012). This is due to the larger surface convective envelopes of cooler stars, extending to hotter layers in the stellar interior. When we look at the trend with ages, we find that among the stars with significant Li detections (all above the solar value of A(Li) = 1.05 dex), we see a tendency of lower Li with increasing age for a given metallicity, as expected (see e.g. the work for solar twins by Carlos et al. (2016) and references therein). The upper envelope for stars at −0.5 < [Fe/H] < 0.5 goes significantly above the Spite plateau due to interstellar and stellar production (see e.g. Prantzos et al. 2017) and even reaches the meteoritic values estimated by Asplund et al. (2009). While we see indications that Li is a good tracer of age for the younger dwarfs in our sample, we do not draw any conclusions because of the uncertain temperature-age causality as well as the influence of the metallicity, detection limits, and other potentially important factors such as rotational velocity.

Appendix A.2: Carbon
In the spectra, only atomic C lines with high excitation energies could be identified, which are strongest in hot and metalrich stars. Due to the high excitation energy of the C lines, in our sample only hot stars with metallicities above −0.3 dex have detectable line strengths. Starting from enhanced abundances at metallicities of −0.75 dex, we see a decreasing trend towards solar metallicities, which flattens at super-solar metallicities. Our results are consistent with the study of Nissen et al. (2014) who demonstrated a linear C-enhancement trend from [Fe/H] = 0 to −0.75 dex. Although it shows a behaviour like an α element, it is expected to follow the iron abundance more closely than these elements and the origin of C is still debated (see Nissen et al. 2014, and references therein). When we look at the trend with ages, we find no significant correlation between C and stellar age for our sample of dwarfs and turn-off stars. Several recent studies, e.g. Masseron & Gilmore (2015), have discovered significant age-trends for evolved stars, which is an observational manifestation of mass-dependent mixing during the dredge-up phase of the stellar evolution.(e.g. Feuillet et al. 2018). Since our study is limited to dwarf stars, we do not detect such correlations. Further, because of the detection limit, we can not draw conclusions regarding the stars with [Fe/H] < 0.

Appendix A.3: Oxygen
For O, we apply 1D non-LTE corrections based on the model atom and non-LTE radiative transfer code described in Amarsi et al. (2015Amarsi et al. ( , 2016a, but using 1D marcs model atmospheres.  Nr. density   Fig. A.5, indicating a large influence of non-LTE for the measured atomic resonance line K I 7699 Å. These effects are estimated to be of the order of −0.4 to −0.7 dex (Ivanova & Shimanskiȋ 2000). This line suffers from an interstellar contribution, which increases the line depth as a function of a reddening-dependent component, which has not been corrected for. It is therefore expected that the LTE trend differs from the non-LTE study by Zhao et al. (2016). K behaves similar to O both when assuming LTE and non-LTE, which indicates that K behaves like an α element. Due to the observational difficulties and strong expected non-LTE effects, we do not draw strong conclusions for this element from our results.
Appendix A.6: Iron-peak elements Scandium shows a similar behaviour as the α elements Si and Ti, i.e. a flat trend for super-solar metallicities and an increase of Sc abundances towards metal-poor stars, in agreement with Battistini & Bensby (2015). Adibekyan et al. (2012) found Sc trends similar to Al, with a significant increase of Al towards supersolar metallicities. All three studies have at least three Sc lines in common, hence it is unlikely that the difference originate in the chosen lines themselves. It is worth mentioning that Adibekyan et al. (2012), using the linelist by Neves et al. (2009), did not include hyper-fine structure splitting, contrary to this study. The difference in abundance measurements with or without the use of hyper-fine structure splitting is of a very complex nature, but has been shown to play an equally important role for abundance estimations as blending, microturbulence velocities, line choices, line centres, and oscillator strengths (Jofré et al. 2017).
Vanadium stays flat over most of the metallicity range, but shows a slight increase towards lower metallicities. A large number of lines are however too blended or weak to be used for the V abundance estimation in the metal-rich and poor regime. The estimated trend is consistent with those estimated by Adibekyan et al. (2012) and Bensby et al. (2014). Where we can detect V, [V/Fe] seems to be rather independent of stellar age, i.e. young and old stars span a range in [V/Fe] around ±0.15.
Chromium traces Fe along most of the metallicity range and our trends agree in general with Bensby et al. (2014) given our lower precision. In the metal-rich regime, super-solar [Cr/Fe] seems to be favoured. We note that Adibekyan et al. (2012) estimated a very slight anti-correlation of Cr with metallicity (around 0.03 dex [Cr/Fe] within 1 dex metallicity), which our measurements do not suggest.
Manganese measurements agree with the decreasing trend of Mn with metallicity found by Battistini & Bensby (2015) and Adibekyan et al. (2012) under the assumption of LTE. Battistini & Bensby (2015) also discuss the significant influence on this trend when considering non-LTE. We only use LTE in our study, but will include non-LTE in future GALAH analysis Cobalt was only detected for a very small number of stars with strong lines (4%), typically metal-rich stars in the sample. We therefore make no conclusions for Co.
Copper was not detected in the metal-poor regime; we observe solar [Cu/Fe] with a small spread at [Fe/H] < 0, an increasing spread of [Cu/Fe] towards solar [Fe/H], and increasing Cu abundance with metallicity in the super-solar regime. This trend agrees with previous studies by Delgado Mena et al. (2017), which was performed on a significantly smaller sample.
Zinc measurements in our sample follow the same trend as Ni and Cu, i.e. showing a rather flat behaviour at low metallicities and stronger increase at super-solar metallicities. Zn is however more scattered than Ni and Cu, although the mean uncertainties are comparable. Zn is not created in SN Type Ia according both to theoretical yields (Iwamoto et al. 1999;Kobayashi et al. 2006) and observations (Nissen & Schuster 2011;Mikolaitis et al. 2017;Skúladóttir et al. 2017). The yields of Zn are very metallicity-dependent (Kobayashi et al. 2006). We therefore do not expect Zn to behave exactly like α-elements. In our sample, stars with high [Zn/Fe] also show high [α/Fe] at and below solar metallicities. At super-solar metallicities, where the highand low-α sequence can not be distinguished, our measurements suggest a linear increase of Zn with a large spread or scatter, as also found by Bensby et al. (2014). Strong blending (especially towards cooler temperatures) decreases the precision of our measurements, but we note that Bensby et al. (2014);Delgado Mena et al. (2017) also estimated a significant spread of Zn with their high-resolution and high S /N data.
Appendix A.7: The neutron-capture elements Several neutron-capture elements have lines that are detectable in GALAH spectra only for giant stars. For this sample of dwarfs and subgiants, the line strength is too small to be significant at the survey, S /N. However, Y and Ba have strong singly ionised lines readily detectable also in the majority of our unevolved stars. We therefore only report the abundance of the s-process neutron capture elements Y and Ba relative to iron.
Yttrium shows a lens-like shape, with a lower almost flat lower limit at about −0.4 dex and a convex upper limit between metallicities of −0.7 and 0.4 dex with most abundances below +0.4 dex. While Bensby et al. (2014) also found some Yenhanced stars, the vast majority of their stars showed roughly solar Y abundances relative to iron. High abundance measurements in stars could be a result of non-LTE effects. However, our sample shows such elevated Y values across a wide range of temperatures and surface gravities, see Fig. A.6. We can not confirm a behaviour similar to O or K, for which LTE-based abundances are strongly overestimated in some stars. Additionally, a strong correlation between [Y/Fe] and age was found for solar twins by Nissen (2015) and Spina et al. (2016), which are insensitive to differential non-LTE effects because of the similar stellar parameters among solar twins. We find old stars of our sample to show depleted Y, while young stars are Y-enhanced. From both theory (Travaglio et al. 2004) and cluster observations (D'Orazi et al. 2009;Maiorca et al. 2011), this can be explained with the increasing contribution of s-process material from lowmass asymptotic giant branch stars to the Y and Ba abundances over time.
Barium was measured less frequently than Y, especially for metal-poor stars, but the overall lens-shaped trend of Ba with metallicity is similar and consistent with Bensby et al. (2014), who showed that a large fraction of their sample with T eff > 6100 K are Ba-enhanced (see their Fig. 16). We note that similar Article number, page 25 of 30 A&A proofs: manuscript no. Buder_GnT_v2 to Delgado Mena et al. (2017), the most metal-rich stars show lower [Ba/Fe] than [Y/Fe] by around 0.1 dex. Investigating the flagged stars, this trend seems to be caused by detection limits for Ba. When we look at the trend with ages, we find these hotter stars are however also younger and their Ba abundance correlates, similar to Y, with the stellar age, as shown by Nissen (2015) and Spina et al. (2016). When comparing [Y/Fe] and [Ba/Fe] in Fig. A.4, we find a reassuringly strong correlation of both, with a slightly steeper slope of Ba abundance compared to Y. In the right panels of Fig. A.4