Search

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 Gaia-ESO Survey is an European Southern Observatory (ESO) public spectroscopic survey that targeted $10^5$ stars in the Milky Way covering the major populations of the disk, bulge and halo. The observations were made using FLAMES on the VLT obtaining both UVES high ($R\sim47,000$) and GIRAFFE medium ($R\sim20,000$) resolution spectra. The analysis of the Gaia-ESO spectra was the work of multiple analysis teams (nodes) within five working groups (WG). The homogenisation of the stellar parameters within WG11 (high resolution observations of FGK stars) and the homogenisation of the stellar parameters within WG10 (medium resolution observations of FGK stars) is described here. In both cases, the homogenisation was carried out using a bayesian Inference method developed specifically for the Gaia-ESO Survey by WG11. The WG10 homogenisation primarily used the cross-match of stars with WG11 as the reference set in both the stellar parameter and chemical abundance homogenisation. In this way the WG10 homogenised results have been placed directly onto the WG11 stellar parameter and chemical abundance scales. The reference set for the metal-poor end was sparse which limited the effectiveness of the homogenisation in that regime. For WG11, the total number of stars for which stellar parameters were derived was 6,231 with typical uncertainties for Teff, log g and [Fe/H] of 32~K, 0.05 and 0.05 respectively. One or more chemical abundances out of a possible 39 elements were derived for 6,188 of the stars. For WG10, the total number of stars for which stellar parameters were derived was 76,675 with typical uncertainties for Teff, log g and [Fe/H] of 64~K, 0.15 and 0.07 respectively. One or more chemical abundances out of a possible 30 elements were derived for 64,177 of the stars., Comment: 31 pages, 19 figures

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

We present an empirical model of age-dependent photospheric lithium depletion, calibrated using a large, homogeneously-analysed sample of 6200 stars in 52 open clusters, with ages from 2--6000 Myr and $-0.3<{\rm [Fe/H}]<0.2$, observed in the Gaia-ESO spectroscopic survey. The model is used to obtain age estimates and posterior age probability distributions from measurements of the Li I 6708A equivalent width for individual (pre) main sequence stars with $3000 < T_{\rm eff}/{\rm K} <6500$, a domain where age determination from the HR diagram is either insensitive or highly model-dependent. In the best cases, precisions of 0.1 dex in log age are achievable; even higher precision can be obtained for coeval groups and associations where the individual age probabilities of their members can be combined. The method is validated on a sample of exoplanet-hosting young stars, finding agreement with claimed young ages for some, but not others. We obtain better than 10 per cent precision in age, and excellent agreement with published ages, for seven well-studied young moving groups. The derived ages for young clusters ($<1$ Gyr) in our sample are also in good agreement with their training ages, and consistent with several published, model-insensitive lithium depletion boundary ages. For older clusters there remain systematic age errors that could be as large as a factor of two. There is no evidence to link these errors to any strong systematic metallicity dependence of (pre) main sequence lithium depletion, at least in the range $-0.29 < {\rm [Fe/H]} < 0.18$. Our methods and model are provided as software -- "Empirical AGes from Lithium Equivalent widthS" (EAGLES)., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. Code available at https://github.com/robdjeff/eagles

The Gaia-ESO Survey is a public spectroscopic survey that has targeted $\gtrsim10^5$ stars covering all major components of the Milky Way from the end of 2011 to 2018, delivering its public final release in May 2022. Unlike other spectroscopic surveys, Gaia-ESO is the only survey that observed stars across all spectral types with dedicated, specialised analyses: from O ($T_\mathrm{eff} \sim 30,000-52,000$~K) all the way to K-M ($\gtrsim$3,500~K). The physics throughout these stellar regimes varies significantly, which has previously prohibited any detailed comparisons between stars of significantly different type. In the final data release (internal data release 6) of the Gaia-ESO Survey, we provide the final database containing a large number of products such as radial velocities, stellar parameters and elemental abundances, rotational velocity, and also, e.g., activity and accretion indicators in young stars and membership probability in star clusters for more than 114,000 stars. The spectral analysis is coordinated by a number of Working Groups (WGs) within the Survey, which specialise in the various stellar samples. Common targets are analysed across WGs to allow for comparisons (and calibrations) amongst instrumental setups and spectral types. Here we describe the procedures employed to ensure all Survey results are placed on a common scale to arrive at a single set of recommended results for all Survey collaborators to use. We also present some general quality and consistency checks performed over all Survey results., Comment: A&A accepted, minor revision, 36 pages, 38 figures

Understanding which sources are present in an astronomical catalogue and which are not is crucial for the accurate interpretation of astronomical data. In particular, for the multidimensional Gaia data, filters and cuts on different parameters or measurements introduces a selection function that may unintentionally alter scientific conclusions in subtle ways. We aim to develop a methodology to estimate the selection function for different sub-samples of stars in the Gaia catalogue. Comparing the number of stars in a given sub-sample to those in the overall Gaia catalogue, provides an estimate of the sub-sample membership probability, as a function of sky position, magnitude and colour. This estimate must differentiate the stochastic absence of sub-sample stars from selection effects. When multiplied with the overall Gaia catalogue selection function this provides the total selection function of the sub-sample. We present the method by estimating the selection function of the sources in Gaia DR3 with heliocentric radial velocity measurements. We also compute the selection function for the stars in the Gaia-Sausage/Enceladus sample, confirming that the apparent asymmetry of its debris across the sky is merely caused by selection effects. The developed method estimates the selection function of the stars present in a sub-sample of Gaia data, given that the sub-sample is completely contained in the Gaia parent catalogue (for which the selection function is known). This tool is made available in a GaiaUnlimited Python package., Comment: 13 pages, 12 figures. Submitted to A&A

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

The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022., Comment: 38 pages. A&A in press

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

In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article (Gilmore et al.) reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. The GES has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 31 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. The final catalogue has been released through the ESO archive at the end of May 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come., Comment: Accepted for publication in Astronomy and Astrophysics. 30 pages, 30 figures, 4 tables

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

Rotationally-induced mode splitting frequencies of low-luminosity subgiants suggest that angular momentum transport mechanisms are 1-2 orders of magnitude more efficient in these stars than predicted by theory. Constraints on the rotation profile of low-luminosity subgiants could be used to identify the dominant mechanism for angular momentum transport. We develop a forward model for the rotation profile given observed rotational splittings, assuming a step-like rotation profile. We identify a consistent degeneracy between the position of the profile discontinuity and the surface rotation rate. We perform mock experiments that show the discontinuity position can be better constrained with a prior on the surface rotation rate, which is informed by star spot modulations. We finally apply this approach to KIC 12508433, a well-studied low-luminosity subgiant, as an example case. With the observed surface rotation prior, we obtain a factor of two increase in precision of the position of strong rotation gradient. We recover the literature values of the core and surface rotation rates and find the highest support for a discontinuity in the radiative zone. Auxiliary measurements of surface rotation could substantially improve inferences on the rotation profile of low-luminosity subgiants with already available data., Comment: 9 pages, 11 figures, 2 table

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

We present chemical abundances for 21 elements (from Li to Eu) in 150 metal-poor Galactic stars spanning $-$4.1 $<$ [Fe/H] $<$ $-$2.1. The targets were selected from the SkyMapper survey and include 90 objects with [Fe/H] $\le$ $-$3 of which some 15 have [Fe/H] $\le$ $-$3.5. When combining the sample with our previous studies, we find that the metallicity distribution function has a power-law slope of $\Delta$(log N)/$\Delta$[Fe/H] = 1.51 $\pm$ 0.01 dex per dex over the range $-$4 $\le$ [Fe/H] $\le$ $-$3. With only seven carbon-enhanced metal-poor stars in the sample, we again find that the selection of metal-poor stars based on SkyMapper filters is biased against highly carbon rich stars for [Fe/H] $>$ $-$3.5. Of the 20 objects for which we could measure nitrogen, 11 are nitrogen-enhanced metal-poor stars. Within our sample, the high NEMP fraction (55\% $\pm$ 21\%) is compatible with the upper range of predicted values (between 12\% and 35\%). The chemical abundance ratios [X/Fe] versus [Fe/H] exhibit similar trends to previous studies of metal-poor stars and Galactic chemical evolution models. We report the discovery of nine new r-I stars, four new r-II stars, one of which is the most metal-poor known, nine low-$\alpha$ stars with [$\alpha$/Fe] $\le$ 0.15 as well as one unusual star with [Zn/Fe] = +1.4 and [Sr/Fe] = +1.2 but with normal [Ba/Fe]. Finally, we combine our sample with literature data to provide the most extensive view of the early chemical enrichment of the Milky Way Galaxy., Comment: MNRAS in press (see source file for full versions of long tables)

Neutron-star mergers were recently confirmed as sites of rapid-neutron-capture (r-process) nucleosynthesis. However, in Galactic chemical evolution models, neutron-star mergers alone cannot reproduce the observed element abundance patterns of extremely metal-poor stars, which indicates the existence of other sites of r-process nucleosynthesis. These sites may be investigated by studying the element abundance patterns of chemically primitive stars in the halo of the Milky Way, because these objects retain the nucleosynthetic signatures of the earliest generation of stars. Here we report the element abundance pattern of the extremely metal-poor star SMSS J200322.54-114203.3. We observe a large enhancement in r-process elements, with very low overall metallicity. The element abundance pattern is well matched by the yields of a single 25-solar-mass magnetorotational hypernova. Such a hypernova could produce not only the r-process elements, but also light elements during stellar evolution, and iron-peak elements during explosive nuclear burning. Hypernovae are often associated with long-duration gamma-ray bursts in the nearby Universe. This connection indicates that similar explosions of fast-spinning strongly magnetized stars occurred during the earliest epochs of star formation in our Galaxy., Comment: Author's version of a Letter published in Nature on July 8th, 2021

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

The recently discovered Indus stellar stream exhibits a diverse chemical signature compared to what is found for most other streams due to the abundances of two outlier stars, Indus$\_$0 and Indus$\_$13. Indus$\_$13, exhibits an extreme enhancement in rapid neutron-capture ($r$-)process elements with $\mathrm{[Eu/Fe]} = +1.81$. It thus provides direct evidence of the accreted nature of $r$-process enhanced stars. In this paper we present a detailed chemical analysis of the neutron-capture elements in Indus$\_$13, revealing the star to be slightly actinide poor. The other outlier, Indus$\_0$, displays a globular cluster-like signature with high N, Na, and Al abundances, while the rest of the Indus stars show abundances compatible with a dwarf galaxy origin. Hence, Indus$\_0$ provides the first chemical evidence of a fully disrupted dwarf containing a globular cluster. We use the chemical signature of the Indus stars to discuss the nature of the stream progenitor which was likely a chemically evolved system, with a mass somewhere in the range from Ursa Minor to Fornax., Comment: 14 pages, 4 figures, and 5 tables, accepted for publication in ApJ

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