Search

The innermost portions of the Milky Way's stellar halo have avoided scrutiny until recently. The lack of wide-area survey data, made it difficult to reconstruct an uninterrupted view of the density distribution of the metal-poor stars inside the Solar radius. In this study, we utilize red giant branch (RGB) stars from Gaia, with metallicities estimated using spectro-photometry from Gaia Data Release 3. Accounting for Gaia's selection function, we examine the spatial distribution of metal-poor ([M/H]<-1.3) RGB stars, from the Galactic centre (r~1 kpc) out to beyond the Solar radius (r~18 kpc). Our best-fitting single-component cored power-law model shows a vertical flattening of ~0.5 and a slope -3.4, consistent with previous studies. Motivated by the mounting evidence for two distinct stellar populations in the inner halo, we additionally test a range of two-component models. One of the components models the tidal debris from the Gaia Sausage/Enceladus merger, while the other captures the Aurora population -- stars that predate the Galactic disk formation. Our best-fit two-component model suggests that both populations contribute equally around the Solar radius, but Aurora dominates the inner halo with a steeper power-law index of -4.5, in agreement with the nitrogen-rich star distribution measured by Horta et al. (2021)., Comment: Submitted to MNRAS

Analysing stellar parameters and abundances from nearly one million Gaia DR3 Radial Velocity Spectrometer (RVS) spectra poses challenges due to the limited spectral coverage (restricted to the infrared Ca II triplet) and variable signal-to-noise ratios of the data. To address this, we use The Cannon, a data-driven method, to transfer stellar parameters and abundances from the GALAH Data Release 4 (DR4; R ~ 28,000) catalogue to the lower resolution Gaia DR3 RVS spectra (R ~ 11,500). Our model, trained on 14,484 common targets, predicts parameters such as Teff, log g, and [Fe/H], along with several other elements across approximately 800,000 Gaia RVS spectra. We utilise stars from open and globular clusters present in the Gaia RVS catalogue to validate our predicted mean [Fe/H] with high precision (~0.02-0.10 dex). Additionally, we recover the bimodal distribution of [Ti/Fe] versus [Fe/H], reflecting the high and low alpha-components of Milky Way disk stars, demonstrating The Cannon's capability for accurate stellar abundance determination from medium-resolution Gaia RVS spectra. The methodologies and resultant catalogue presented in this work highlight the remarkable potential of the RVS dataset, which by the end of the Gaia mission will comprise spectra of over 200 million stars., Comment: Submitted to MNRAS, 16 pages, 15 figures

The radial velocity catalog from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) is unique in its simultaneously large volume and high precision as a result of its decade-long survey duration, multiplexing (600 fibers), and spectral resolution of $R \sim 22,500$. However, previous data reductions of APOGEE have not fully realized the potential radial velocity (RV) precision of the instrument. Here we present an RV catalog based on a new reduction of all 2.6 million visits of APOGEE DR17 and validate it against improved estimates for the theoretical RV performance. The core ideas of the new reduction are the simultaneous modeling of all components in the spectra, rather than a separate subtraction of point estimates for the sky, and a marginalization over stellar types, rather than a grid search for an optimum. We show that this catalog, when restricted to RVs measured with the same fiber, achieves noise-limited precision down to 30 m/s and delivers well-calibrated uncertainties. We also introduce a general method for calibrating fiber-to-fiber constant RV offsets and demonstrate its importance for high RV precision work in multi-fiber spectrographs. After calibration, we achieve 47 m/s RV precision on the combined catalog with RVs measured with different fibers. This degradation in precision relative to measurements with only a single fiber suggests that refining line spread function models should be a focus in SDSS-V to improve the fiber-unified RV catalog., Comment: 22 pages, 13 figures, SDSS-V Technical Paper, Submitted to AJ

Context. The explosive burning that drives nova eruptions results in unique nucleosynthesis that heavily over-produces certain isotopes relative to the solar abundance. However, novae are often ignored when considering the chemical evolution of our Galaxy due to their low ejecta masses. Aims. In this work, we use previously computed synthetic nova populations and the galactic chemical evolution code OMEGA+ to assess the impact that novae have on the evolution of stable elemental and isotopic abundances. Methods. We combine populations of novae computed using the binary population synthesis code binary_c with the galactic chemical evolution code OMEGA+ and detailed, white dwarf mass-dependent nova yields to model the nucleosynthetic contributions of novae to the evolution of the Milky Way. We consider three different nova yield profiles, each corresponding to a different set of nova yield calculations. Results. Despite novae from low-mass white dwarfs (WDs) dominating nova ejecta contributions, we find that novae occurring on massive WDs are still able to contribute significantly to many isotopes, particularly those with high mass numbers. We find that novae can produce up to 35% of the Galactic 13C and 15N mass by the time the model Galaxy reaches [Fe/H] = 0, and earlier in the evolution of the Galaxy (between [Fe/H] = -2 and -1) novae may have been the dominant source of 15N. Predictions for [13C/Fe], [15N/Fe], 12C/13C, and 14N/15N abundances ratios vary by up to 0.2 dex at [Fe/H] = 0 and by up to 0.7 dex in [15N/Fe] and 14N/15N between [Fe/H] = -2 and -1 (corresponding approximately to Galactic ages of 170 Myr and 1 Gyr in our model). The Galactic evolution of other stable isotopes (excluding Li) is not noticeably affected by including novae.

Spectral synthesis codes are essential for inferring stellar parameters and detailed chemical abundances. These codes require many physical inputs to predict an emergent spectrum. Developers adopt the best measurements of those inputs at the time they release their code, but those measurements usually improve over time faster than the software is updated. In general, the impact of using incorrect or uncertain dissociation energies are largely unknown. Here we evaluate how incorrect dissociation energies impact abundances measured from C2, CN, CH, TiO, and MgO features. For each molecule we synthesised optical spectra of FGKM-type main-sequence and giant stars using the literature dissociation energy, and an incorrect (perturbed) dissociation energy. We find that the uncertainties in the dissociation energies adopted by spectral synthesis codes for CN, CH, TiO, and MgO lead to negligible differences in flux or abundances. C2 is the only diatomic molecule where the uncertainty of the inputted dissociation energy translates to a significant difference in flux, and carbon abundance differences of up to 0.2 dex. For Solar-like stars, the impact on carbon abundance is up to 0.09 dex. These large abundance differences demonstrate the importance of updating the inputs adopted by spectral synthesis codes, as well as a consensus on appropriate values between different codes., Comment: 7 pages, 4 figures

Stellar mergers and accretion events have been crucial in shaping the evolution of the Milky Way (MW). These events have been dynamically identified and chemically characterised using red giants and main-sequence stars. RR Lyrae (RRL) variables can play a crucial role in tracing the early formation of the MW since they are ubiquitous, old (t$\ge$10 Gyr) low-mass stars and accurate distance indicators. We exploited Data Release 3 of the GALAH survey to identify 78 field RRLs suitable for chemical analysis. Using synthetic spectra calculations, we determined atmospheric parameters and abundances of Fe, Mg, Ca, Y, and Ba. Most of our stars exhibit halo-like chemical compositions, with an iron peak around [Fe/H]$\approx -$1.40, and enhanced Ca and Mg content. Notably, we discovered a metal-rich tail, with [Fe/H] values ranging from $-$1 to approximately solar metallicity. This sub-group includes almost ~1/4 of the sample, it is characterised by thin disc kinematics and displays sub-solar $\alpha$-element abundances, marginally consistent with the majority of the MW stars. Surprisingly, they differ distinctly from typical MW disc stars in terms of the s-process elements Y and Ba. We took advantage of similar data available in the literature and built a total sample of 535 field RRLs for which we estimated kinematical and dynamical properties. We found that metal-rich RRLs (1/3 of the sample) likely represent an old component of the MW thin disc. We also detected RRLs with retrograde orbits and provided preliminary associations with the Gaia-Sausage-Enceladus, Helmi, Sequoia, Sagittarius, and Thamnos stellar streams., Comment: Accepted for publication in MNRAS. 29 pages, 20 figures

Gaia can not individually resolve very close binary systems, however, the collected data can still be used to identify them. A powerful indicator of stellar multiplicity is the sources reported Renormalized Unit Weight Error (ruwe), which effectively captures the astrometric deviations from single-source solutions. We aim to characterise the imprints left on ruwe caused by binarity. By flagging potential binary systems based on ruwe, we aim to characterise which of their properties will contribute the most to their detectability. We develop a model to estimate ruwe values for observations of Gaia sources, based on the biases to the single-source astrometric track arising from the presence of an unseen companion. Then, using the recipes from previous GaiaUnlimited selection functions, we estimate the selection probability of sources with high ruwe, and discuss what binary properties contribute to increasing the sources ruwe. We compute the maximum ruwe value which is compatible with single-source solutions as a function of their location on-sky. We see that binary systems selected as sources with a ruwe higher than this sky-varying threshold have a strong detectability window in their orbital period distribution, which peaks at periods equal to the Gaia observation time baseline. We demonstrate how our sky-varying ruwe threshold provides a more complete sample of binary systems when compared to single sky-averaged values by studying the unresolved binary population in the Gaia Catalogue of Nearby Stars. We provide the code and tools used in this study, as well as the sky-varying ruwe threshold through the GaiaUnlimited Python package, Comment: 11 pages, 10 figures. Accepted for publication in A&A

When there are many observations of an astronomical source - many images with different dithers, or many spectra taken at different barycentric velocities - it is standard practice to shift and stack the data, to (for example) make a high signal-to-noise average image or mean spectrum. Bound-saturating measurements are made by manipulating a likelihood function, where the data are treated as fixed, and model parameters are modified to fit the data. Traditional shifting and stacking of data can be converted into a model-fitting procedure, such that the data are not modified, and yet the output is the shift-adjusted mean. The key component of this conversion is a spectral model that is completely flexible but also a continuous function of wavelength (or position in the case of imaging) that can represent any signal being measured by the device after any reasonable translation (or rotation or field distortion). The benefits of a modeling approach are myriad: The sacred data never are modified. Noise maps, data gaps, and bad-data masks don't require interpolation. The output can take the form of an image or spectrum evaluated on a pixel grid, as is traditional. In addition to shifts, the model can account for line-spread or point-spread function variations, world-coordinate-system variations, and calibration or normalization variations. The noise in the output becomes uncorrelated across neighboring pixels as the shifts deliver good coverage in some sense. The only cost is a small increase in computational complexity over that of traditional methods. We demonstrate the method with a small data example and we provide open source sample code for re-use., Comment: submitted to AAS Journals

Detailed chemical studies of F/G/K -- or Solar-type -- stars have long been routine in stellar astrophysics, enabling studies in both Galactic chemodynamics, and exoplanet demographics. However, similar understanding of the chemistry of M and late-K dwarfs -- the most common stars in the Galaxy -- has been greatly hampered both observationally and theoretically by the complex molecular chemistry of their atmospheres. Here we present a new implementation of the data-driven \textit{Cannon} model, modelling $T_{\rm eff}$, $\log g$, [Fe/H], and [Ti/Fe] trained on low-medium resolution optical spectra ($4\,000-7\,000\,$\SI{}{\angstrom}) from 103 cool dwarf benchmarks. Alongside this, we also investigate the sensitivity of optical wavelengths to various atomic and molecular species using both data-driven and theoretical means via a custom grid of MARCS synthetic spectra, and make recommendations for where MARCS struggles to reproduce cool dwarf fluxes. Under leave-one-out cross-validation, our \textit{Cannon} model is capable of recovering $T_{\rm eff}$, $\log g$, [Fe/H], and [Ti/Fe] with precisions of 1.4\%, $\pm0.04\,$dex, $\pm0.10\,$dex, and $\pm0.06\,$dex respectively, with the recovery of [Ti/Fe] pointing to the as-yet mostly untapped potential of exploiting the abundant -- but complex -- chemical information within optical spectra of cool stars., Comment: 24 pages, 14 figures, 4 tables. Accepted for publication in MNRAS

The spatial distribution of Galactic stars with different chemical abundances encodes information on the processes that drove the formation and evolution of the Milky Way. Survey selection functions are indispensable for analysing astronomical catalogues produced by large-scale surveys. The use of these selection functions in data modelling is more complex when data from different surveys are to be modelled simultaneously. We introduce a procedure for constructing the selection function of a sample of red clump stars that have parallaxes and elemental abundances from the Gaia mission. We separately constructed the selection function of the APOGEE DR17 red clump stars, which depends on very different observables and has a very different spatial coverage. We combined the two surveys and accounted for their joint selection function to provide strong constraints on the radial and vertical density distribution of mono-abundance populations, with Gaia offering a dense coverage of the solar neighbourhood, while APOGEE reaches larger distances near the Galactic plane. We confirm that the radial density profile steepens with increasing metallicity. The combined sample also indicates a metallicity-dependent flaring of the alpha-poor disc. We provide the code for constructing the Gaia selection function we used in this study through the GaiaUnlimited Python package., Comment: 21 pages, 18 figures, accepted for publication in Astronomy & Astrophysics

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys., Comment: 11 pages + 22 page appendix, accepted to ApJL

We describe several updates to Korg, a package for 1D LTE spectral synthesis of FGKM stars. Built-in functions to fit observed spectra via synthesis or equivalent widths make it easy to take advantage of Korg's automatic differentiation. Comparison to a past analysis of 18 Sco shows that we obtain significantly reduced line-to-line abundance scatter with Korg. Fitting and synthesis are facilitated by a rigorously-tested model atmosphere interpolation method, which introduces negligible error to synthesized spectra for stars with $T_\mathrm{eff} \gtrsim 4000\,\mathrm{K}$. For cooler stars, atmosphere interpolation is complicated by the presence of molecules, though we demonstrate an adequate method for cool dwarfs. The chemical equilibrium solver has been extended to include polyatomic and charged molecules, extending Korg's regime of applicability to M stars. We also discuss a common oversight regarding the synthesis of hydrogen lines in the infrared, and show that Korg's Brackett line profiles are a much closer match to observations than others available. Documentation, installation instructions, and tutorials are available at https://github.com/ajwheeler/Korg.jl., Comment: updated to reflect published version

To accurately measure a star's atmospheric parameters and chemical abundances, it is crucial to have high-quality spectra. Analysing the detailed chemical abundances of groups of stars can help us better understand nucleosynthesis, galactic chemical enrichment, and stellar evolution. In this study, we explored whether stellar spots can affect a star's inferred metallicity and, if so, where the impact is the strongest. To investigate this, we created synthetic infrared spectra that included stellar spots for a sample of main-sequence stars younger than the sun. We then applied two models to the data: one that accounted for spots and one that did not. From this, we can determine the bias introduced when fitting spotted spectra with a non-spotted model and how this bias varies with different parameters. Our findings revealed that fitting spotted spectra with a non-spotted model can introduce a scatter of up to 0.05 dex in the inferred metallicity, especially for stars with high levels of spot coverage. This bias is similar in magnitude to other relevant effects, such as atomic diffusion, radiative levitation, or non-local thermodynamic equilibrium. We also found that the effect is most pronounced in young stars and decreases with age. These results suggest that stellar spots can introduce a systematic uncertainty in metallicity that is not currently accounted for in spectroscopic analysis. This could potentially limit scientific inferences for population-level studies of young stars and differential abundance analyses., Comment: 9 pages, 7 Figures

Stellar models predict that lithium (Li) inside a star is destroyed during the first dredge-up phase, yet 1.2% of red giant stars are Li-rich. We aim to uncover possible origins of this population, by analysing 1155 Li-rich giants (A(Li) $\geq$ 1.5) in GALAH DR3. To expose peculiar traits of Li-rich stars, we construct a reference sample of Li-normal (doppelg\"anger) stars with matched evolutionary state and fiducial supernova abundances. Comparing Li-rich and doppelg\"anger spectra reveals systematic differences in the H-$\alpha$ and Ca-triplet line profiles associated with the velocity broadening measurement. We also find twice as many Li-rich stars appear to be fast rotators (2% with $v_\textrm{broad} \gtrsim 20$ km s$^{-1}$) compared to doppelg\"angers. On average, Li-rich stars have higher abundances than their doppelg\"angers, for a subset of elements, and Li-rich stars at the base of RGB have higher mean $s-$process abundances ($\geq 0.05$ dex for Ba, Y, Zr), relative to their doppelg\"angers. External mass-transfer from intermediate-mass AGB companions could explain this signature. Additional companion analysis excludes binaries with mass ratios $\gtrsim$ 0.5 at $\gtrsim$ 7 AU. We also discover that highly Ba-enriched stars are missing from the Li-rich population, possibly due to low-mass AGB companions which preclude Li-enrichment. Finally, we confirm a prevalence of Li-rich stars on the red clump that increases with lithium, which supports an evolutionary state mechanism for Li-enhancement. Multiple culprits, including binary spin-up and mass-transfer, are therefore likely mechanisms of Li-enrichment., Comment: 29 pages, 19 figures, 6 tables. Submitted to ApJ

The progenitor systems and explosion mechanisms responsible for the thermonuclear events observationally classified as Type Ia supernovae are uncertain and difficult to uniquely constrain using traditional observations of Type Ia supernova host galaxies, progenitors, light curves, and remnants. For the subset of thermonuclear events that are prolific producers of iron, we use published theoretical nucleosynthetic yields to identify a set of elemental abundance ratios infrequently observed in metal-poor stars but shared across a range of progenitor systems and explosion mechanisms: [Na,Mg,Co/Fe]<0. We label stars with this abundance signature ``iron-rich metal-poor'' or IRMP stars. We suggest that IRMP stars formed in environments dominated by thermonuclear nucleosynthesis and consequently that their elemental abundances can be used to constrain both the progenitor systems and explosion mechanisms responsible for thermonuclear explosions. We identify three IRMP in the literature and homogeneously infer their elemental abundances. We find that the elemental abundances of BD+80 245, HE 0533--5340, and SMSS J034249.53--284216.0 are best explained by the (double) detonations of sub-Chandrasekhar mass CO white dwarfs. If our interpretation of IRMP stars is accurate, then they should be very rare in globular clusters and more common in the Magellanic Clouds and dwarf spheroidal galaxies than in the Milky Way's halo. We propose that future studies of IRMP stars will quantify the relative occurrences of different thermonuclear event progenitor systems and explosion mechanisms., Comment: Accepted for publication at AJ

The eighteenth data release of the Sloan Digital Sky Surveys (SDSS) is the first one for SDSS-V, the fifth generation of the survey. SDSS-V comprises three primary scientific programs, or "Mappers": Milky Way Mapper (MWM), Black Hole Mapper (BHM), and Local Volume Mapper (LVM). This data release contains extensive targeting information for the two multi-object spectroscopy programs (MWM and BHM), including input catalogs and selection functions for their numerous scientific objectives. We describe the production of the targeting databases and their calibration- and scientifically-focused components. DR18 also includes ~25,000 new SDSS spectra and supplemental information for X-ray sources identified by eROSITA in its eFEDS field. We present updates to some of the SDSS software pipelines and preview changes anticipated for DR19. We also describe three value-added catalogs (VACs) based on SDSS-IV data that have been published since DR17, and one VAC based on the SDSS-V data in the eFEDS field., Comment: Accepted to ApJS

We present Korg, a new package for 1D LTE (local thermodynamic equilibrium) spectral synthesis of FGK stars, which computes theoretical spectra from the near-ultraviolet to the near-infrared, and implements both plane-parallel and spherical radiative transfer. We outline the inputs and internals of Korg, and compare synthetic spectra from Korg, MOOG, Turbospectrum, and SME. The disagreements between Korg and the other codes are no larger than those between the other codes, although disagreement between codes is substantial. We examine the case of a C$_2$ band in detail, finding that uncertainties on physical inputs to spectral synthesis account for a significant fraction of the disagreement. Korg is 1-100 times faster than other codes in typical use, compatible with automatic differentiation libraries, and easily extensible, making it ideal for statistical inference and parameter estimation applied to large data sets. Documentation and installation instructions are available at https://ajwheeler.github.io/Korg.jl/stable/., Comment: to appear in ApJ; comments welcome

Interpreting and modelling astronomical catalogues requires an understanding of the catalogues' completeness or selection function: objects of what properties had a chance to end up in the catalogue. Here we set out to empirically quantify the completeness of the overall Gaia DR3 catalogue. This task is not straightforward because Gaia is the all-sky optical survey with the highest angular resolution to date and no consistent ``ground truth'' exists to allow direct comparisons. However, well-characterised deeper imaging enables an empirical assessment of Gaia's $G$-band completeness across parts of the sky. On this basis, we devised a simple analytical completeness model of Gaia as a function of the observed $G$ magnitude and position over the sky, which accounts for both the effects of crowding and the complex Gaia scanning law. Our model only depends on a single quantity: the median magnitude $M_{10}$ in a patch of the sky of catalogued sources with $\texttt{astrometric_matched_transits}$ $\leq 10$. $M_{10}$ reflects elementary completeness decisions in the Gaia pipeline and is computable from the Gaia DR3 catalogue itself and therefore applicable across the whole sky. We calibrate our model using the Dark Energy Camera Plane Survey (DECaPS) and test its predictions against Hubble Space Telescope observations of globular clusters. We find that our model predicts Gaia's completeness values to a few per cent across the sky. We make the model available as a part of the $\texttt{gaiasf}$ Python package built and maintained by the GaiaUnlimited project: $\texttt{https://github.com/gaia-unlimited/gaiaunlimited}$, Comment: submitted to A&A

Of all the light elements, the evolution of lithium (Li) in the Milky Way is perhaps the most difficult to explain. Li is difficult to synthesize and is easily destroyed, making most stellar sites unsuitable for producing Li in sufficient quantities to account for the proto-solar abundance. For decades, novae have been proposed as a potential explanation to this 'Galactic Li problem', and the recent detection of 7Be in the ejecta of multiple nova eruptions has breathed new life into this theory. In this work, we assess the viability of novae as dominant producers of Li in the Milky Way. We present the most comprehensive treatment of novae in a galactic chemical evolution code to date, testing theoretical- and observationally-derived nova Li yields by integrating metallicity-dependent nova ejecta profiles computed using the binary population synthesis code binary c with the galactic chemical evolution code OMEGA+. We find that our galactic chemical evolution models which use observationally-derived Li yields account for the proto-solar Li abundance very well, while models relying on theoretical nova yields cannot reproduce the proto-solar observation. A brief exploration of physical uncertainties including single-stellar yields, the metallicity resolution of our nova treatment, common-envelope physics, and nova accretion efficiencies indicates that this result is robust to physical assumptions. Scatter within the observationally-derived Li yields in novae is identified as the primary source of uncertainty, motivating further observations of 7Be in nova ejecta., Comment: Accepted for publication in ApJL 28/7/2022

Extremely Metal-Poor (EMP) stars provide a valuable probe of early chemical enrichment in the Milky Way. Here we leverage a large sample of $\sim600,000$ high-resolution stellar spectra from the GALAH survey plus a machine learning algorithm to find 54 candidates with estimated [Fe/H]~$\leq$~-3.0, 6 of which have [Fe/H]~$\leq$~-3.5. Our sample includes $\sim 20 \%$ main sequence EMP candidates, unusually high for \emp surveys. We find the magnitude-limited metallicity distribution function of our sample is consistent with previous work that used more complex selection criteria. The method we present has significant potential for application to the next generation of massive stellar spectroscopic surveys, which will expand the available spectroscopic data well into the millions of stars., Comment: 27 pages, 20 figures, accepted for publication in ApJ, candidate table available at this https://github.com/arvhug/GALAH---TSNE_EMP

Pioneering photometric, astrometric, and spectroscopic surveys are helping exoplanetary scientists better constrain the fundamental properties of stars within our galaxy, and the planets these stars host. In this study, we use the third data release from the stellar spectroscopic GALAH Survey, coupled with astrometric data of eDR3 from the \textit{Gaia} satellite, and other data from NASA's Exoplanet Archive, to refine our understanding of 279 confirmed and candidate exoplanet host stars and their exoplanets. This homogenously analysed data set comprises 105 confirmed exoplanets, along with 146 K2 candidates, 95 TESS Objects of Interest (TOIs) and 52 Community TOIs (CTOIs). Our analysis significantly shifts several previously (unknown) planet parameters while decreasing the uncertainties for others; Our radius estimates suggest that 35 planet candidates are more likely brown dwarfs or stellar companions due to their new radius values. We are able to refine the radii and masses of WASP-47 e, K2-106 b, and CoRoT-7 b to their most precise values yet, to less than 2.3\% and 8.5\% respectively. We also use stellar rotational values from GALAH to show that most planet candidates will have mass measurements that will be tough to obtain with current ground-based spectrographs. With GALAH's chemical abundances, we show through chemo-kinematics that there are five planet-hosts that are associated with the galaxy's thick disc, including NGTS-4, K2-183 and K2-337. Finally, we show there is no statistical difference between the chemical properties of hot Neptune and hot rocky exoplanet hosts, with the possibility that short-period rocky worlds might be the remnant cores of hotter, gaseous worlds., Comment: 21 pages with 11 figures and 5 tables. Accepted into MNRAS with appendix tables available within this arXiv submission

The metallicity of a star affects its evolution in a variety of ways, changing stellar radii, luminosities, lifetimes, and remnant properties. In this work, we use the population synthesis code binary_c to study how metallicity affects novae in the context of binary stellar evolution. We compute a 16-point grid of metallicities ranging from $Z=10^{-4}$ to 0.03, presenting distributions of nova white dwarf masses, accretion rates, delay-times, and initial system properties at the two extremes of our 16-point metallicity grid. We find a clear anti-correlation between metallicity and the number of novae produced, with the number of novae at $Z=0.03$ roughly half that at $Z=10^{-4}$. The white dwarf mass distribution has a strong systematic variation with metallicity, while the shape of the accretion rate distribution is relatively insensitive. We compute a current nova rate of approximately 33 novae per year for the Milky Way, a result consistent with observational estimates relying on extra-Galactic novae but an under-prediction relative to observational estimates relying on Galactic novae. However, the shape of our predicted Galactic white dwarf mass distribution differs significantly to existing observationally derived distributions, likely due to our underlying physical assumptions. In M31, we compute a current nova rate of approximately 36 novae per year, under-predicting the most recent observational estimate of $65^{+15}_{-16}$. Finally, we conclude that when making predictions about currently observable nova rates in spiral galaxies, or stellar environments where star formation has ceased in the distant past, metallicity can likely be considered of secondary importance compared to uncertainties in binary stellar evolution.

Orbital characteristics based on Gaia Early Data Release 3 astrometric parameters are analyzed for ${\sim} 4000$ metal-poor stars ([Fe/H] $\leq -0.8$) compiled from the Best $\&$ Brightest survey. Selected as metal-poor candidates based on broadband near- and far-IR photometry, $43\%$ of these stars had medium-resolution ($1200 \lesssim R \lesssim 2000$) validation spectra obtained over a seven-year campaign from $2014$ to $2020$ with a variety of telescopes. The remaining stars were chosen based on photometric metallicity determinations from the Huang et al. recalibration of the Sky Mapper Southern Survey. Dynamical clusters of these stars are obtained from the orbital energy and cylindrical actions using the \HDBSCAN ~unsupervised learning algorithm. We identify $52$ Dynamically Tagged Groups (DTGs) with between $5$ and $21$ members; $18$ DTGs have at least $10$ member stars. Milky Way (MW) substructures such as Gaia-Sausage-Enceladus, the Metal-Weak Thick-Disk, Thamnos, the Splashed Disk, and the Helmi Stream are identified. Associations with MW globular clusters are determined for $8$ DTGs; no recognized MW dwarf galaxies were associated with any of our DTGs. Previously identified dynamical groups are also associated with our DTGs, with emphasis placed on their structural determination and possible new identifications. Chemically peculiar stars are identified as members of several DTGs, with $6$ DTGs that are associated with \textit{r}-process-enhanced stars. We demonstrate that the mean carbon and $\alpha$-element abundances of our DTGs are correlated with their mean metallicity in an understandable manner. Similarly, we find that the mean metallicity, carbon, and $\alpha$-element abundances are separable into different regions of the mean rotational-velocity space.

'Publish or perish' is an expression describing the pressure on academics to consistently publish research to ensure a successful career in academia. With a global pandemic that has changed the world, how has it changed academic productivity? Here we show that academics are posting just as many publications on the arXiv pre-print server as if there were no pandemic: 168,630 were posted in 2020, a +12.6% change from 2019 and $+1.4\sigma$ deviation above the predicted 162,577 $\pm$ 4,393. However, some immediate impacts are visible in individual research fields. Conference cancellations have led to sharp drops in pre-prints, but laboratory closures have had mixed effects. Only some experimental fields show mild declines in outputs, with most being consistent on previous years or even increasing above model expectations. The most significant change is a 50% increase ($+8\sigma$) in quantitative biology research, all related to the COVID-19 pandemic. Some of these publications are by biologists using arXiv for the first time, and some are written by researchers from other fields (e.g., physicists, mathematicians). While quantitative biology pre-prints have returned to pre-pandemic levels, 20% of the research in this field is now focussed on the COVID-19 pandemic, demonstrating a strong shift in research focus., Comment: Submitted to RSOS

Stellar abundances and ages afford the means to link chemical enrichment to galactic formation. In the Milky Way, individual element abundances show tight correlations with age, which vary in slope across ([Fe/H]-[$\alpha$/Fe]). Here, we step from characterising abundances as measures of age, to understanding how abundances trace properties of stellar birth-environment in the disk over time. Using measurements from $\sim$27,000 APOGEE stars (R=22,500, SNR$>$200), we build simple local linear models to predict a sample of elements (X = Si, O, Ca, Ti, Ni, Al, Mn, Cr) using (Fe, Mg) abundances alone, as fiducial tracers of supernovae production channels. Given [Fe/H] and [Mg/H], we predict these elements, [X/H], to about double the uncertainty of their measurements. The intrinsic dispersion, after subtracting measurement errors in quadrature is $\approx 0.015-0.04$~dex. The residuals of the prediction (measurement $-$ model) for each element demonstrate that each element has an individual link to birth properties at fixed (Fe, Mg). Residuals from primarily massive-star supernovae (i.e. Si, O, Al) partially correlate with guiding radius. Residuals from primarily supernovae Ia (i.e. Mn, Ni) partially correlate with age. A fraction of the intrinsic scatter that persists at fixed (Fe, Mg), however, after accounting for correlations, does not appear to further discriminate between birth properties that can be traced with present-day measurements. Presumably, this is because the residuals are also, in part, a measure of the typical (in)-homogeneity of the disk's stellar birth environments, previously inferred only using open-cluster systems. Our study implies at fixed birth radius and time, there is a median scatter of $\approx 0.01-0.015$ dex in elements generated in supernovae sources., Comment: submitted to ApJ, 16 pages, 12 figures

Since the advent of $Gaia$ astrometry, it is possible to identify massive accreted systems within the Galaxy through their unique dynamical signatures. One such system, $Gaia$-Sausage-Enceladus (GSE), appears to be an early "building block" given its virial mass $> 10^{10}\,\mathrm{M_\odot}$ at infall ($z\sim1-3$). In order to separate the progenitor population from the background stars, we investigate its chemical properties with up to 30 element abundances from the GALAH+ Survey Data Release 3 (DR3). To inform our choice of elements for purely chemically selecting accreted stars, we analyse 4164 stars with low-$\alpha$ abundances and halo kinematics. These are most different to the Milky Way stars for abundances of Mg, Si, Na, Al, Mn, Fe, Ni, and Cu. Based on the significance of abundance differences and detection rates, we apply Gaussian mixture models to various element abundance combinations. We find the most populated and least contaminated component, which we confirm to represent GSE, contains 1049 stars selected via [Na/Fe] vs. [Mg/Mn] in GALAH+ DR3. We provide tables of our selections and report the chrono-chemodynamical properties (age, chemistry, and dynamics). Through a previously reported clean dynamical selection of GSE stars, including $30 < \sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} < 55$, we can characterise an unprecedented 24 abundances of this structure with GALAH+ DR3. Our chemical selection allows us to prevent circular reasoning and characterise the dynamical properties of the GSE, for example mean $\sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} = 26_{-14}^{+9}$. We find only $(29\pm1)\%$ of the GSE stars within the clean dynamical selection region. Our methodology will improve future studies of accreted structures and their importance for the formation of the Milky Way., Comment: 29 pages, 18 figures, 6 tables; publicshed in MNRAS; all code to recreate results (based on public data described in arXiv:2011.02505) and figures available at https://github.com/svenbuder/Accreted-stars-in-GALAH-DR3

The Phoenix stellar stream has a low intrinsic dispersion in velocity and metallicity that implies the progenitor was probably a low mass globular cluster. In this work we use Magellan/MIKE high-dispersion spectroscopy of eight Phoenix stream red giants to confirm this scenario. In particular, we find negligible intrinsic scatter in metallicity ($\sigma(\mathrm{[Fe~II/H]}) = 0.04^{+0.11}_{-0.03}$) and a large peak-to-peak range in [Na/Fe] and [Al/Fe] abundance ratios, consistent with the light element abundance patterns seen in the most metal-poor globular clusters. However, unlike any other globular cluster, we also find an intrinsic spread in [Sr II/Fe] spanning $\sim$1 dex, while [Ba II/Fe] shows nearly no intrinsic spread ($\sigma(\mathrm{[Ba~II/H]}) = {0.03}^{+0.10}_{-0.02}$). This abundance signature is best interpreted as slow neutron capture element production from a massive fast-rotating metal-poor star ($15-20 \mathrm{M}_\odot$, $v_\mathrm{ini}/v_\mathrm{crit} = 0.4$, $[\mathrm{Fe/H}] = -3.8$). The low inferred cluster mass suggests the system would have been unable to retain supernovae ejecta, implying that any massive fast-rotating metal-poor star that enriched the interstellar medium must have formed and evolved before the globular cluster formed. Neutron capture element production from asymptotic giant branch stars or magneto-rotational instabilities in core-collapse supernovae provide poor fits to the observations. We also report one Phoenix stream star to be a lithium-rich giant ($A(\mathrm{Li}) = 3.1 \pm 0.1$). At $[\mathrm{Fe/H}] = -2.93$ it is among the most metal-poor lithium-rich giants known., Comment: Accepted to ApJ 2021-07-09

Stellar members of binary systems are formed from the same material, therefore they should be chemically identical. However, recent high-precision studies have unveiled chemical differences between the two members of binary pairs composed by Sun-like stars. The very existence of these chemically inhomogeneous binaries represents one of the most contradictory examples in stellar astrophysics and source of tension between theory and observations. It is still unclear whether the abundance variations are the result of chemical inhomogeneities in the protostellar gas clouds or instead if they are due to planet engulfment events occurred after the stellar formation. While the former scenario would undermine the belief that the chemical makeup of a star provides the fossil information of the environment where it formed, a key assumption made by several studies of our Galaxy, the second scenario would shed light on the possible evolutionary paths of planetary systems. Here, we perform a statistical study on 107 binary systems composed by Sun-like stars to provide - for the first time - unambiguous evidence in favour of the planet engulfment scenario. We also establish that planet engulfment events occur in stars similar to our own Sun with a probability ranging between 20 and 35$\%$. This implies that a significant fraction of planetary systems undergo very dynamical evolutionary paths that can critically modify their architectures, unlike our Solar System which has preserved its planets on nearly circular orbits. This study also opens to the possibility of using chemical abundances of stars to identify which ones are the most likely to host analogues of the calm Solar System., Comment: This is an authors' version of the manuscript prior to formal peer review at Nature Astronomy

The chemical abundances of a galaxy's metal-poor stellar population can be used to investigate the earliest stages of its formation and chemical evolution. The Magellanic Clouds are the most massive of the Milky Way's satellite galaxies and are thought to have evolved in isolation until their recent accretion by the Milky Way. Unlike the Milky Way's less massive satellites, little is know about the Magellanic Clouds' metal-poor stars. We have used the mid-infrared metal-poor star selection of Schlaufman & Casey (2014) and archival data to target nine LMC and four SMC giants for high-resolution Magellan/MIKE spectroscopy. These nine LMC giants with $-2.4\lesssim[\text{Fe/H}]\lesssim-1.5$ and four SMC giants with $-2.6\lesssim[\text{Fe/H}]\lesssim-2.0$ are the most metal-poor stars in the Magellanic Clouds yet subject to a comprehensive abundance analysis. While we find that at constant metallicity these stars are similar to Milky Way stars in their $\alpha$, light, and iron-peak elemental abundances, both the LMC and SMC are enhanced relative to the Milky Way in the $r$-process element europium. These abundance offsets are highly significant, equivalent to $3.9\sigma$ for the LMC, $2.7\sigma$ for the SMC, and $5.0\sigma$ for the complete Magellanic Cloud sample. We propose that the $r$-process enhancement of the Magellanic Clouds' metal-poor stellar population is a result of the Magellanic Clouds' isolated chemical evolution and long history of accretion from the cosmic web combined with $r$-process nucleosynthesis on a timescale longer than the core-collapse supernova timescale but shorter than or comparable to the thermonuclear (i.e., Type Ia) supernova timescale., Comment: Accepted for publication in The Astronomical Journal. 30 pages, 7 figures, and 5 tables

Measuring the collapse time of a binary neutron star merger remnant can inform the physics of extreme matter and improve modelling of short gamma-ray bursts and associated kilonova. The lifetime of the post-merger remnant directly impacts the mechanisms available for the jet launch of short gamma-ray bursts. We develop and test a method to measure the collapse time of post-merger remnants. We show that for a GW170817-like event at $\sim\!40\,$Mpc, a network of Einstein Telescope with Cosmic Explorer is required to detect collapse times of $\sim\!10\,$ms. For a two-detector network at A+ design sensitivity, post-merger remnants with collapse times of $\sim\!10\,\mathrm{ms}$ must be $\lesssim 10\,$Mpc to be measureable. This increases to $\sim\!18-26\,$Mpc if we include the proposed Neutron star Extreme Matter Observatory (NEMO), increasing the effective volume by a factor of $\sim\!30$., Comment: 9 pages, 3 figures

Novae are some of the most commonly detected optical transients and have the potential to provide valuable information about binary evolution. Binary population synthesis codes have emerged as the most effective tool for modelling populations of binary systems, but such codes have traditionally employed greatly simplified nova physics, precluding detailed study. In this work, we implement a model treating H and He novae as individual events into the binary population synthesis code \binaryc. This treatment of novae represents a significant improvement on the `averaging' treatment currently employed in modern population synthesis codes. We discuss the evolutionary pathways leading to these phenomena and present nova event rates and distributions of several important physical parameters. Most novae are produced on massive white dwarfs, with approximately 70 and 55 per cent of nova events occurring on O/Ne white dwarfs for H and He novae respectively. Only 15 per cent of H-nova systems undergo a common-envelope phase, but these systems are responsible for the majority of H nova events. All He-accreting He-nova systems are considered post-common-envelope systems, and almost all will merge with their donor star in a gravitational-wave driven inspiral. We estimate the current annual rate of novae in M31 (Andromeda) to be approximately $41 \pm 4$ for H novae, underpredicting the current observational estimate of $65^{+15}_{-16}$, and $0.14\pm0.015$ for He novae. When varying common-envelope parameters, the H nova rate varies between 20 and 80 events per year., Comment: Accepted, MNRAS. 7 Jun 2020: Minor correction regarding AM CVn masses at period bounce, courtesy of P. Neuteufel

The results from the ESA Gaia astrometric mission and deep photometric surveys have revolutionized our knowledge of the Milky Way. There are many ongoing efforts to search these data for stellar substructure to find evidence of individual accretion events that built up the Milky Way and its halo. One of these newly identified features, called Nyx, was announced as an accreted stellar stream traveling in the plane of the disk. Using a combination of elemental abundances and stellar parameters from the GALAH and APOGEE surveys, we find that the abundances of the highest likelihood Nyx members are entirely consistent with membership of the thick disk, and inconsistent with a dwarf galaxy origin. We conclude that the postulated Nyx stream is most probably a high-velocity component of the Milky Way's thick disk. With the growing availability of large data sets including kinematics, stellar parameters, and detailed abundances, the probability of detecting chance associations increases, and hence new searches for substructure require confirmation across as many data dimensions as possible., Comment: Accepted for publication by ApJ Letters

We present the results of spectroscopic follow-up for 1897 low-metallicity star candidates, selected from the Best & Brightest (B&B) Survey, carried out with the GMOS-N/S (Gemini North/South telescopes) and Goodman (SOAR Telescope) spectrographs. From these low-resolution ($R \sim 2000$) spectra, we estimate stellar atmospheric parameters, as well as carbon and magnesium (representative of $\alpha$ elements) abundance ratios. We confirm that $56\%$ of our program stars are metal-poor ([Fe/H] $< -1.0$), $30\%$ are very metal-poor (VMP; [Fe/H] $< -2.0$) and $2\%$ are extremely metal-poor (EMP; [Fe/H] $< -3.0$). There are 191 carbon-enhanced metal-poor (CEMP) stars, resulting in CEMP fractions of $19\%$ and $43\%$ for the VMP and EMP regimes, respectively. A total of 94 confirmed CEMP stars belong to Group I ($A({\rm C}) \gtrsim 7.25$) and 97 to Group II ($A({\rm C}) \lesssim 7.25$) in the Yoon-Beers $A$(C)$-$[Fe/H] diagram. Moreover, we combine these data with Gaia EDR3 astrometric information to delineate new target-selection criteria, which have been applied to the Goodman/SOAR candidates, to more than double the efficiency for identification of bona-fide VMP and EMP stars in comparison to random draws from the B&B catalog. We demonstrate that this target-selection approach can achieve success rates of $96\%$, $76\%$, $28\%$ and $4\%$ for [Fe/H] $\leq -1.5$, $\leq -2.0$, $\leq -2.5$ and $\leq -3.0$, respectively. Finally, we investigate the presence of dynamically interesting stars in our sample. We find that several VMP/EMP ([Fe/H] $\leq -2.5$) stars can be associated with either the disk system or halo substructures like Gaia-Sausage/Enceladus and Sequoia., Comment: Accepted for ApJ. Long tables at the end of the manuscript. Minor typos fixed and references updated to match the published version

GALAH+ is a magnitude-limited survey of high resolution stellar spectra obtained by the HERMES spectrograph at the Australian Astronomical Observatory. Its third data release provides reduced spectra with new derivations of stellar parameters and abundances of 30 chemical elements for 584,015 dwarfs and giants, 88% of them in the Gaia magnitude range 11 < G < 14. Here we use these improved values of stellar parameters to build a library of observed spectra which is useful to study variations of individual spectral lines with stellar parameters. This and other improvements are used to derive radial velocities with uncertainties which are generally within 0.1 km/s or ~25% smaller than in the previous release. Median differences in radial velocities measured here and by the Gaia DR2 or APOGEE DR16 surveys are smaller than 30 m/s, a larger offset is present only for Gaia measurements of giant stars. We identify 4483 stars with intrinsically variable velocities and 225 stars for which the velocity stays constant over >=3 visits spanning more than a year. The combination of radial velocities from GALAH+ with distances and sky plane motions from Gaia enables studies of dynamics within streams and clusters. For example, we estimate that the open cluster M67 has a total mass of ~3300 Msun and its outer parts seem to be expanding, though astrometry with a larger time-span than currently available from Gaia eDR3 is needed to judge if the latter result is real., Comment: 15 pages, 4 tables, 17 figures + online data files, accepted for publication in MNRAS

Previous studies have found that the elemental abundances of a star correlate directly with its age and metallicity. Using this knowledge, we derive ages for a sample of 250,000 stars taken from GALAH DR3 using only their overall metallicity and chemical abundances. Stellar ages are estimated via the machine learning algorithm $XGBoost$, using main sequence turnoff stars with precise ages as our input training set. We find that the stellar ages for the bulk of the GALAH DR3 sample are accurate to 1-2 Gyr using this method. With these ages, we replicate many recent results on the age-kinematic trends of the nearby disk, including the age-velocity dispersion relationship of the solar neighborhood and the larger global velocity dispersion relations of the disk found using $Gaia$ and GALAH. The fact that chemical abundances alone can be used to determine a reliable age for a star have profound implications for the future study of the Galaxy as well as upcoming spectroscopic surveys. These results show that the chemical abundance variation at a given birth radius is quite small, and imply that strong chemical tagging of stars directly to birth clusters may prove difficult with our current elemental abundance precision. Our results highlight the need of spectroscopic surveys to deliver precision abundances for as many nucleosynthetic production sites as possible in order to estimate reliable ages for stars directly from their chemical abundances. Applying the methods outlined in this paper opens a new door into studies of the kinematic structure and evolution of the disk, as ages may potentially be estimated for a large fraction of stars in existing spectroscopic surveys. This would yield a sample of millions of stars with reliable age determinations, and allow precise constraints to be put on various kinematic processes in the disk, such as the efficiency and timescales of radial migration., Comment: 13 pages, 15 figures, submitted to MNRAS

The ESA Gaia astrometric mission has enabled the remarkable discovery that a large fraction of the stars near the Solar neighbourhood appear to be debris from a single in-falling system, the so-called Gaia-Enceladus-Sausage (GSE). One exciting feature of this result is that it gives astronomers for the first time a large sample of easily observable unevolved stars that formed in an extra-Galactic environment, which can be compared to stars that formed within our Milky Way. Here we use these stars to investigate the "Spite Plateau" -- the near-constant lithium abundance observed in metal-poor dwarf stars across a wide range of metallicities (-3<[Fe/H]<-1). In particular our aim is to test whether the stars that formed in the GSE show a different Spite Plateau to other Milky Way stars that inhabit the disk and halo. Individual galaxies could have different Spite Plateaus -- e.g., the ISM could be more depleted in lithium in a lower galactic mass system due to it having a smaller reservoir of gas. We identified 76 GSE dwarf stars observed and analyzed by the GALactic Archeology with HERMES (GALAH) survey as part of its Third Data Release. Orbital actions were used to select samples of Gaia-Enceladus stars, and comparison samples of halo and disk stars. We find that the Gaia-Enceladus stars show the same lithium abundance as other likely accreted stars and in situ Milky Way stars, strongly suggesting that the "lithium problem" is not a consequence of the formation environment. This result fits within the growing consensus that the Spite Plateau, and more generally the "cosmological lithium problem" -- the observed discrepancy between the amount of lithium in warm, metal-poor dwarf stars in our Galaxy, and the amount of lithium predicted to have been produced by Big Bang Nucleosynthesis -- is the result of lithium depletion processes within stars., Comment: 12 pages, 7 figures, updated version resubmitted to MNRAS after responding to positive referee report. Also now has attached galah_spite_stars.csv file which is full version of Table 1

Open clusters are unique tracers of the history of our own Galaxy's disk. According to our membership analysis based on \textit{Gaia} astrometry, out of the 226 potential clusters falling in the footprint of GALAH or APOGEE, we find that 205 have secure members that were observed by at least one of the survey. Furthermore, members of 134 clusters have high-quality spectroscopic data that we use to determine their chemical composition. We leverage this information to study the chemical distribution throughout the Galactic disk of 21 elements, from C to Eu. The radial metallicity gradient obtained from our analysis is $-$0.076$\pm$0.009 dex kpc$^{-1}$, which is in agreement with previous works based on smaller samples. Furthermore, the gradient in the [Fe/H] - guiding radius (r$_{\rm guid}$) plane is $-$0.073$\pm$0.008 dex kpc$^{-1}$. We show consistently that open clusters trace the distribution of chemical elements throughout the Galactic disk differently than field stars. In particular, at given radius, open clusters show an age-metallicity relation that has less scatter than field stars. As such scatter is often interpreted as an effect of radial migration, we suggest that these differences are due to the physical selection effect imposed by our Galaxy: clusters that would have migrated significantly also had higher chances to get destroyed. Finally, our results reveal trends in the [X/Fe]$-$r$_{\rm guid}$$-$age space, which are important to understand production rates of different elements as a function of space and time., Comment: 18 pages, 10 figures, accepted for publication MNRAS

The ensemble of chemical element abundance measurements for stars, along with precision distances and orbit properties, provides high-dimensional data to study the evolution of the Milky Way. With this third data release of the Galactic Archaeology with HERMES (GALAH) survey, we publish 678 423 spectra for 588 571 mostly nearby stars (81.2% of stars are within <2 kpc), observed with the HERMES spectrograph at the Anglo-Australian Telescope. This release (hereafter GALAH+ DR3) includes all observations from GALAH Phase 1 (bright, main, and faint survey, 70%), K2-HERMES (17%), TESS-HERMES (5%), and a subset of ancillary observations (8%) including the bulge and >75 stellar clusters. We derive stellar parameters $T_\text{eff}$, $\log g$, [Fe/H], $v_\text{mic}$, $v_\text{broad}$ & $v_\text{rad}$ using our modified version of the spectrum synthesis code Spectroscopy Made Easy (SME) and 1D MARCS model atmospheres. We break spectroscopic degeneracies in our spectrum analysis with astrometry from $Gaia$ DR2 and photometry from 2MASS. We report abundance ratios [X/Fe] for 30 different elements (11 of which are based on non-LTE computations) covering five nucleosynthetic pathways. We describe validations for accuracy and precision, flagging of peculiar stars/measurements and recommendations for using our results. Our catalogue comprises 65% dwarfs, 34% giants, and 1% other/unclassified stars. Based on unflagged chemical composition and age, we find 62% young low-$\alpha$, 9% young high-$\alpha$, 27% old high-$\alpha$, and 2% stars with $\mathrm{[Fe/H]} \leq -1$. Based on kinematics, 4% are halo stars. Several Value-Added-Catalogues, including stellar ages and dynamics, updated after $Gaia$ eDR3, accompany this release and allow chrono-chemodynamic analyses, as we showcase., Comment: 48 pages, 27 figures, 5 pages appendix, accepted for publication in MNRAS, comments welcome, GALAH+ DR3 can be found at https://docs.datacentral.org.au/galah/dr3/overview/, Catalogs can be downloaded from https://cloud.datacentral.org.au/teamdata/GALAH/public/GALAH_DR3/

We describe the atomic and molecular data that were used for the abundance analyses of FGK-type stars carried out within the Gaia-ESO Survey. We present an unprecedented effort to create a homogeneous line list, which was used by several abundance analysis groups to calculate synthetic spectra and equivalent widths. The atomic data are accompanied by quality indicators and detailed references to the sources. The atomic and molecular data are made publicly available in electronic form. In general experimental transition probabilities were preferred but theoretical values were also used. Astrophysical gf-values were avoided due to the model-dependence of such a procedure. For elements whose lines are significantly affected by hyperfine structure or isotopic splitting a concerted effort has been made to collate the necessary data for the individual line components. We also performed a detailed investigation of available data for line broadening due to collisions with neutral hydrogen atoms. Synthetic spectra calculated for the Sun and Arcturus were used to assess the blending properties of the lines. Among a subset of over 1300 lines of 35 elements in the wavelength ranges from 475 nm to 685 nm and from 850 nm to 895 nm we identified about 200 lines of 24 species which have accurate gf-values and are free of blends in the spectra of the Sun and Arcturus. For the broadening due to collisions with neutral hydrogen we recommend data based on Anstee-Barklem-O'Mara theory, where available, and to avoid lines of neutral species otherwise. Theoretical broadening data by R.L. Kurucz should be used for Sc II, Ti II, and Y II lines. For ionised rare-earth species the Uns\"old approximation with an enhancement factor of 1.5 for the line width can be used. Desirable improvements in atomic data were identified for a number of species, including Al I, S I, Cr II, Na I, Si I, Ca II, and Ni I., Comment: Accepted for publication in section 12. Atomic, molecular, and nuclear data of Astronomy and Astrophysics; main part 25 pages, 8 tables, 5 figures; appendices 53 pages, 23 tables, 29 figures

We present high-resolution Magellan/MIKE spectroscopy of 42 red giant stars in seven stellar streams confirmed by the Southern Stellar Stream Spectroscopic Survey (S5): ATLAS, Aliqa Uma, Chenab, Elqui, Indus, Jhelum, and Phoenix. Abundances of 30 elements have been derived from over 10,000 individual line measurements or upper limits using photometric stellar parameters and a standard LTE analysis. This is currently the most extensive set of element abundances for stars in stellar streams. Three streams (ATLAS, Aliqa Uma, and Phoenix) are disrupted metal-poor globular clusters, although only weak evidence is seen for the light element anticorrelations commonly observed in globular clusters. Four streams (Chenab, Elqui, Indus, and Jhelum) are disrupted dwarf galaxies, and their stars display abundance signatures that suggest progenitors with stellar masses ranging from $10^6-10^7 M_\odot$. Extensive description is provided for the analysis methods, including the derivation of a new method for including the effect of stellar parameter correlations on each star's abundance and uncertainty. This paper includes data gathered with the 6.5 meter Magellan Telescopes located at Las Campanas Observatory, Chile., Comment: Accepted to AJ. 39 pages, 13 figures, 7 tables

Globular clusters are some of the oldest bound stellar structures observed in the Universe. They are ubiquitous in large galaxies and are believed to trace intense star formation events and the hierarchical build-up of structure. Observations of globular clusters in the Milky Way, and a wide variety of other galaxies, have found evidence for a `metallicity floor', whereby no globular clusters are found with chemical (`metal') abundances below approximately 0.3 to 0.4 per cent of that of the Sun. The existence of this metallicity floor may reflect a minimum mass and a maximum redshift for surviving globular clusters to form, both critical components for understanding the build-up of mass in the universe. Here we report measurements from the Southern Stellar Streams Spectroscopic Survey of the spatially thin, dynamically cold Phoenix stellar stream in the halo of the Milky Way. The properties of the Phoenix stream are consistent with it being the tidally disrupted remains of a globular cluster. However, its metal abundance ([Fe/H] = -2.7) is substantially below that of the empirical metallicity floor. The Phoenix stream thus represents the debris of the most metal-poor globular cluster discovered so far, and its progenitor is distinct from the present-day globular cluster population in the local Universe. Its existence implies that globular clusters below the metallicity floor have probably existed, but were destroyed during Galactic evolution., Comment: Authors' version of an Article published in Nature on July 29th, 2020

The bulge is the oldest component of the Milky Way. Since numerous simulations of Milky Way formation have predicted that the oldest stars at a given metallicity are found on tightly bound orbits, the Galaxy's oldest stars are likely metal-poor stars in the inner bulge with small apocenters (i.e., $R_{\mathrm{apo}}\lesssim4$ kpc). In the past, stars with these properties have been impossible to find due to extreme reddening and extinction along the line of sight to the inner bulge. We have used the mid-infrared metal-poor star selection of Schlaufman & Casey (2014) on Spitzer/GLIMPSE data to overcome these problems and target candidate inner bulge metal-poor giants for moderate-resolution spectroscopy with AAT/AAOmega. We used those data to select three confirmed metal-poor giants ($[\mathrm{Fe/H}]=-3.15,-2.56,-2.03$) for follow-up high-resolution Magellan/MIKE spectroscopy. A comprehensive orbit analysis using Gaia DR2 astrometry and our measured radial velocities confirms that these stars are tightly bound inner bulge stars. We determine the elemental abundances of each star and find high titanium and iron-peak abundances relative to iron in our most metal-poor star. We propose that the distinct abundance signature we detect is a product of nucleosynthesis in the Chandrasekhar-mass thermonuclear supernova of a CO white dwarf accreting from a helium star with a delay time of about 10 Myr. Even though chemical evolution is expected to occur quickly in the bulge, the intense star formation in the core of the nascent Milky Way was apparently able to produce at least one Chandrasekhar-mass thermonuclear supernova progenitor before chemical evolution advanced beyond $[\mathrm{Fe/H}]\sim-3$., Comment: AAS Accepted in July 24, 2020. 25 pages, 9 figures

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., Comment: 5 pages, 2 figures, Accepted for publication by MNRAS

Detection and parameter estimation of binary neutron star merger remnants can shed light on the physics of hot matter at supranuclear densities. Here we develop a fast, simple model that can generate gravitational waveforms, and show it can be used for both detection and parameter estimation of post-merger remnants. The model consists of three exponentially-damped sinusoids with a linear frequency-drift term. The median fitting factors between the model waveforms and numerical-relativity simulations exceed 0.90. We detect remnants at a post-merger signal-to-noise ratio of $\ge 7$ using a Bayes-factor detection statistic with a threshold of 3000. We can constrain the primary post-merger frequency to $\pm_{1.2}^{1.4}\%$ at post-merger signal-to-noise ratios of 15 with an increase in precision to $\pm_{0.2}^{0.3}\%$ for post-merger signal-to-noise ratios of 50. The tidal coupling constant can be constrained to $\pm^{9}_{12}\%$ at post-merger signal-to-noise ratios of 15, and $\pm 5\%$ at post-merger signal-to-noise ratios of 50 using a hierarchical inference model.

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 complete results for planet-candidate hosts from the K2-HERMES survey, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra for more than 30,000 K2 stars. We present complete host-star parameters and planet-candidate radii for 224 K2 candidate planets from C1-C13. Our results cast severe doubt on 30 K2 candidates, as we derive unphysically large radii, larger than 2 $R_{Jup}$. This work highlights the importance of obtaining accurate, precise, and self-consistent stellar parameters for ongoing large planet search programs - something that will only become more important in the coming years, as TESS begins to deliver its own harvest of exoplanets., Comment: Accepted for publication in MNRAS. Full versions of online-only Tables 2 and 3 available on request. Same for Table 4, in exchange for the trick to make MNRAS style break long tables across pages

Magnetic fields and stellar spots can alter the equivalent widths of absorption lines in stellar spectra, varying during the activity cycle. This also influences the information that we derive through spectroscopic analysis. In this study we analyse high-resolution spectra of 211 Sun-like stars observed at different phases of their activity cycles, in order to investigate how stellar activity affects the spectroscopic determination of stellar parameters and chemical abundances. We observe that equivalent widths of lines can increase as a function of the activity index log R$^\prime_{\rm HK}$ during the stellar cycle, which also produces an artificial growth of the stellar microturbulence and a decrease in effective temperature and metallicity. This effect is visible for stars with activity indexes log R$^\prime_{\rm HK}$$\geq$$-$5.0 (i.e., younger than 4-5 Gyr) and it is more significant at higher activity levels. These results have fundamental implications on several topics in astrophysics that are discussed in the paper, including stellar nucleosynthesis, chemical tagging, the study of Galactic chemical evolution, chemically anomalous stars, the structure of the Milky Way disk, stellar formation rates, photoevaporation of circumstellar disks, and planet hunting., Comment: Accepted for publication by ApJ. 19 pages, 6 figures, 6 tables

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

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

Chemical tagging seeks to identify unique star formation sites from present-day stellar abundances. Previous techniques have treated each abundance dimension as being statistically independent, despite theoretical expectations that many elements can be produced by more than one nucleosynthetic process. In this work we introduce a data-driven model of nucleosynthesis where a set of latent factors (e.g., nucleosynthetic yields) contribute to all stars with different scores, and clustering (e.g., chemical tagging) is modelled by a mixture of multivariate Gaussians in a lower-dimensional latent space. We use an exact method to simultaneously estimate the factor scores for each star, the partial assignment of each star to each cluster, and the latent factors common to all stars, even in the presence of missing data entries. We use an information-theoretic Bayesian principle to estimate the number of latent factors and clusters. Using the second Galah data release we find that six latent factors are preferred to explain N = 2,566 stars with 17 chemical abundances. We identify the rapid- and slow-neutron capture processes, as well as latent factors consistent with Fe-peak and \alpha-element production, and another where K and Zn dominate. When we consider N ~ 160,000 stars with missing abundances we find another 7 factors, as well as 16 components in latent space. Despite these components showing separation in chemistry that is explained through different yield contributions, none show significant structure in their positions or motions. We argue that more data, and joint priors on cluster membership that are constrained by dynamical models, are necessary to realise chemical tagging at a galactic-scale. We release software that allows for model parameters to be optimised in seconds given a fixed number of latent factors, components, and $10^7$ abundance measurements., Comment: Accepted to ApJ

Observations of stellar remnants linked to Type Ia and Type Iax supernovae are necessary to fully understand their progenitors. Multiple progenitor scenarios predict a population of kicked donor remnants and partially-burnt primary remnants, both moving with relatively high velocity. But only a handful of examples consistent with these two predicted populations have been observed. Here we report the likely first known example of an unbound white dwarf that is consistent with being the fully-cooled primary remnant to a Type Iax supernova. The candidate, LP 93-21, is travelling with a galactocentric velocity of v$_{gal}$ $\simeq$ 605 km s$^{-1}$, and is gravitationally unbound to the Milky Way. We rule out an extragalactic origin. The Type Iax supernova ejection scenario is consistent with its peculiar unbound trajectory, given anomalous elemental abundances are detected in its photosphere via spectroscopic follow-up. This discovery reflects recent models that suggest stellar ejections likely occur often. Unfortunately the intrinsic faintness of white dwarfs, and the uncertainty associated with their direct progenitor systems, makes it difficult to detect and confirm such donors., Comment: 8 pages, 4 figures, 1 table