Search

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

Theoretical models of stellar evolution predict that most of the lithium inside a star is destroyed as the star becomes a red giant. However, observations reveal that about 1% of red giants are peculiarly rich in lithium, often exceeding the amount in the interstellar medium or predicted from the Big Bang. With only about 150 lithium-rich giants discovered in the past four decades, and no distinguishing properties other than lithium enhancement, the origin of lithium-rich giant stars is one of the oldest problems in stellar astrophysics. Here we report the discovery of 2,330 low-mass (1 to 3$\,M_\odot$) lithium-rich giant stars, which we argue are consistent with internal lithium production that is driven by tidal spin-up by a binary companion. Our sample reveals that most lithium-rich giants have helium-burning cores ($80^{+7}_{-6}\%$), and that the frequency of lithium-rich giants rises with increasing stellar metallicity. We find that while planet accretion may explain some lithium-rich giants, it cannot account for the majority that have helium-burning cores. We rule out most other proposed explanations as the primary mechanism for lithium-rich giants, including all stages related to single star evolution. Our analysis shows that giants remain lithium-rich for only about two million years. A prediction from this lithium depletion timescale is that most lithium-rich giants with a helium-burning core have a binary companion., Comment: ApJ accepted. Catalog of 2,330 lithium-rich giants available at https://zenodo.org/record/3515777

Abridged: Getting a better understanding of the evolution and nucleosynthetic yields of the most metal-poor stars (appr. Z<=10^-5) is critical because they are part of the big picture of the history of the primitive Universe. Yet many of the remaining unknowns of stellar evolution lie in the birth, life, and death of these objects. We review stellar evolution of intermediate-mass (IMS) Z<=10-5 models existing in the literature, with a focus on the problem of their final fates. The depth and efficiency of mixing episodes are critical to determine the mass limits for the formation of electron-capture supernovae, but our knowledge of these phenomena is not complete because they are strongly affected by the choice of input physics. We also consider the alternative SNI1/2 channel to form SNe out of the most metal-poor IMS. In this case, it is critical to understand the thermally-pulsing AGB evolution until the late stages. Efficient second dredge-up and, later, third dredge-up episodes could be able to pollute stellar envelopes enough for the stars to undergo thermal pulses in a way very similar to that of higher initial Z objects. Inefficient 2nd and/or 3rd dredge-up may leave an almost pristine envelope, unable to sustain strong stellar winds. This may allow the H-exhausted core to grow to M_Ch before the envelope is lost, and thus let the star explode as a SNI1/2. After reviewing the information available on these two possible channels for the formation of SNe, we discuss existing nucleosynthetic yields of stars of metallicity Z<=10^-5, and present an example of nucleosynthetic calculations for a thermally-pulsing Super-AGB star of Z=10^-5. We compare theoretical predictions with observations of the lowest [Fe/H] objects detected. The review closes by discussing current open questions as well as possible fruitful avenues for future research., Comment: 36 pages, 13 figures, 3 tables, PASA review, accepted

Extremely metal-poor (EMP) stars are an integral piece in the puzzle that is the early Universe, and although anomolous subclasses of EMP stars such as carbon-enhanced metal-poor (CEMP) stars are well-studied, they make up less than half of all EMP stars with [Fe/H] $\sim -3.0$. The amount of carbon depletion occurring on the red giant branch (carbon offset) is used to determine the evolutionary status of EMP stars, and this offset will differ between CEMP and normal EMP stars. The depletion mechanism employed in stellar models (from which carbon offfsets are derived) is very important, however the only widely available carbon offsets in the literature are derived from stellar models using a thermohaline mixing mechanism that cannot simultaneously match carbon and lithium abundances to observations for a single diffusion coeffcient. Our stellar evolution models utilise a modified thermohaline mixing model that can match carbon and lithium in the metal-poor globular cluster NGC 6397. We compare our models to the bulk of the EMP star sample at [Fe/H] $= -3$ and show that our modified models follow the trend of the observations and deplete less carbon compared to the standard thermohaline mixing theory. We conclude that stellar models that employ the standard thermohaline mixing formalism overestimate carbon offsets and hence CEMP star frequencies, particularly at metallicities where carbon-normal stars dominate the EMP star population., Comment: 8 pages, 4 figures, 1 table

Stars with unusual elemental abundances offer clues about rare astrophysical events or nucleosynthetic pathways. Stars with significantly depleted magnesium and enhanced potassium ([Mg/Fe] < -0.5; [K/Fe] > 1) have to date only been found in the massive globular cluster NGC 2419 and, to a lesser extent, NGC 2808. The origin of this abundance signature remains unknown, as does the reason for its apparent exclusivity to these two globular clusters. Here we present 112 field stars, identified from 454,180 LAMOST giants, that show significantly enhanced [K/Fe] and possibly depleted [Mg/Fe] abundance ratios. Our sample spans a wide range of metallicities (-1.5 < [Fe/H] < 0.3), yet none show abundance ratios of [K/Fe] or [Mg/Fe] that are as extreme as those observed in NGC 2419. If confirmed, the identified sample of stars represents evidence that the nucleosynthetic process producing the anomalous abundances ratios of [K/Fe] and [Mg/Fe] probably occurs at a wide range of metallicities. This would suggest that pollution scenarios that are limited to early epochs (such as Population III supernovae) are an unlikely explanation, although they cannot be ruled out entirely. This sample is expected to help guide modelling attempts to explain the origin of the Mg-K abundance signature.

A series of astronomical observations obtained over the period 1986 to 2018 supports the idea that life is a cosmic rather than a purely terrestrial or planetary phenomenon. These include (1) the detection of biologically relevant molecules in interstellar clouds and in comets, (2) mid-infrared spectra of interstellar grains and the dust from comets, (3) a diverse set of data from comets including the Rosetta mission showing consistency with biology and (4) the frequency of Earth-like or habitable planets in the Galaxy. We argue that the conjunction of all the available data suggests the operation of cometary biology and interstellar panspermia rather than the much weaker hypothesis of comets being only the source of the chemical building blocks of life. We conclude with specific predictions on the properties expected of extra-terrestrial life if it is discovered on Enceladus, Europa or beyond. A radically different biochemistry elsewhere can be considered as a falsification of the theory of interstellar panspermia., Comment: 13 pages, 5 figures, submitted to PASA

Theoretical predictions for the core helium burning phase of stellar evolution are highly sensitive to the uncertain treatment of mixing at convective boundaries. In the last few years, interest in constraining the uncertain structure of their deep interiors has been renewed by insights from asteroseismology. Recently, Spruit (2015) proposed a limit for the rate of growth of helium-burning convective cores based on the higher buoyancy of material ingested from outside the convective core. In this paper we test the implications of such a limit for stellar models with a range of initial mass and metallicity. We find that the constraint on mixing beyond the Schwarzschild boundary has a significant effect on the evolution late in core helium burning, when core breathing pulses occur and the ingestion rate of helium is fastest. Ordinarily, core breathing pulses prolong the core helium burning lifetime to such an extent that models are at odds with observations of globular cluster populations. Across a wide range of initial stellar masses ($0.83 \leq M/\text{M}_\odot \leq 5$), applying the Spruit constraint reduces the core helium burning lifetime because core breathing pulses are either avoided or their number and severity reduced. The constraint suggested by Spruit therefore helps to resolve significant discrepancies between observations and theoretical predictions. Specifically, we find improved agreement for $R_2$, the observed ratio of asymptotic giant branch to horizontal branch stars in globular clusters; the luminosity difference between these two groups; and in asteroseismology, the mixed-mode period spacing detected in red clump stars in the \textit{Kepler} field., Comment: Accepted for publication in MNRAS; 11 pages, 6 figures

Thermohaline mixing is a favoured mechanism for the so-called "extra mixing" on the red giant branch of low-mass stars. The mixing is triggered by the molecular weight inversion created above the hydrogen shell during first dredge-up when helium-3 burns via 3He(3He,2p)4He. The standard 1D diffusive mixing scheme cannot simultaneously match carbon and lithium abundances to NGC 6397 red giants. We investigate two modifications to the standard scheme: (1) an advective two stream mixing algorithm, and (2) modifications to the standard 1D thermohaline mixing formalism. We cannot simultaneously match carbon and lithium abundances using our two stream mixing approach. However we develop a modified diffusive scheme with an explicit temperature dependence that can simultaneously fit carbon and lithium abundances to NGC 6397 stars. Our modified diffusive scheme induces mixing that is faster than the standard theory predicts in the hotter part of the thermohaline region and mixing that is slower in the cooler part. Our results infer that the extra mixing mechanism needs further investigation and more observations are required, particularly for stars in different clusters spanning a range in metallicity., Comment: 15 pages, 15 figures

We review the lives, deaths and nucleosynthetic signatures of intermediate mass stars in the range approximately 6.5-12 Msun, which form super-AGB stars near the end of their lives. We examine the critical mass boundaries both between different types of massive white dwarfs (CO, CO-Ne, ONe) and between white dwarfs and supernovae and discuss the relative fraction of super-AGB stars that end life as either an ONe white dwarf or as a neutron star (or an ONeFe white dwarf), after undergoing an electron capture supernova. We also discuss the contribution of the other potential single-star channels to electron-capture supernovae, that of the failed massive stars. We describe the factors that influence these different final fates and mass limits, such as composition, the efficiency of convection, rotation, nuclear reaction rates, mass loss rates, and third dredge-up efficiency. We stress the importance of the binary evolution channels for producing electron-capture supernovae. We discuss recent nucleosynthesis calculations and elemental yield results and present a new set of s-process heavy element yield predictions. We assess the contribution from super-AGB star nucleosynthesis in a Galactic perspective, and consider the (super-)AGB scenario in the context of the multiple stellar populations seen in globular clusters. A brief summary of recent works on dust production is included. Lastly we conclude with a discussion of the observational constraints and potential future advances for study into these stars on the low mass/high mass star boundary., Comment: 28 pages, 11 figures. Invited review for Publications of the Astronomical Society of Australia, to be published in special issue on "Electron Capture Supernovae". Submitted

The treatment of convective boundaries during core helium burning is a fundamental problem in stellar evolution calculations. In Paper~I we showed that new asteroseismic observations of these stars imply they have either very large convective cores or semiconvection/partially mixed zones that trap g-modes. We probe this mixing by inferring the relative lifetimes of asymptotic giant branch (AGB) and horizontal branch (HB) from $R_2$, the observed ratio of these stars in recent HST photometry of 48 Galactic globular clusters. Our new determinations of $R_2$ are more self-consistent than those of previous studies and our overall calculation of $R_2 = 0.117 \pm 0.005$ is the most statistically robust now available. We also establish that the luminosity difference between the HB and the AGB clump is $\Delta \log{L}_\text{HB}^\text{AGB} = 0.455 \pm 0.012$. Our results accord with earlier findings that standard models predict a lower $R_2$ than is observed. We demonstrate that the dominant sources of uncertainty in models are the prescription for mixing and the stochastic effects that can result from its numerical treatment. The luminosity probability density functions that we derive from observations feature a sharp peak near the AGB clump. This constitutes a strong new argument against core breathing pulses, which broaden the predicted width of the peak. We conclude that the two mixing schemes that can match the asteroseismology are capable of matching globular cluster observations, but only if (i) core breathing pulses are avoided in models with a semiconvection/partially mixed zone, or (ii) that models with large convective cores have a particular depth of mixing beneath the Schwarzschild boundary during subsequent early-AGB `gravonuclear' convection., Comment: 22 pages, 14 figures. Accepted for publication in MNRAS

There is now strong evidence that some stars have been born with He mass fractions as high as $Y \approx 0.40$ (e.g., in $\omega$ Centauri). However, the advanced evolution, chemical yields, and final fates of He-rich stars are largely unexplored. We investigate the consequences of He-enhancement on the evolution and nucleosynthesis of intermediate-mass asymptotic giant branch (AGB) models of 3, 4, 5, and 6 M$_\odot$ with a metallicity of $Z = 0.0006$ ([Fe/H] $\approx -1.4$). We compare models with He-enhanced compositions ($Y=0.30, 0.35, 0.40$) to those with primordial He ($Y=0.24$). We find that the minimum initial mass for C burning and super-AGB stars with CO(Ne) or ONe cores decreases from above our highest mass of 6 M$_\odot$ to $\sim$ 4-5 M$_\odot$ with $Y=0.40$. We also model the production of trans-Fe elements via the slow neutron-capture process (s-process). He-enhancement substantially reduces the third dredge-up efficiency and the stellar yields of s-process elements (e.g., 90% less Ba for 6 M$_\odot$, $Y=0.40$). An exception occurs for 3 M$_\odot$, where the near-doubling in the number of thermal pulses with $Y=0.40$ leads to $\sim$ 50% higher yields of Ba-peak elements and Pb if the $^{13}$C neutron source is included. However, the thinner intershell and increased temperatures at the base of the convective envelope with $Y=0.40$ probably inhibit the $^{13}$C neutron source at this mass. Future chemical evolution models with our yields might explain the evolution of s-process elements among He-rich stars in $\omega$ Centauri., Comment: 21 pages, 16 figures, accepted for publication by MNRAS. Stellar yields included as online data tables

The detection of mixed oscillation modes offers a unique insight into the internal structure of core helium burning (CHeB) stars. The stellar structure during CHeB is very uncertain because the growth of the convective core, and/or the development of a semiconvection zone, is critically dependent on the treatment of convective boundaries. In this study we calculate a suite of stellar structure models and their non-radial pulsations to investigate why the predicted asymptotic g-mode $\ell = 1$ period spacing $\Delta\Pi_1$ is systematically lower than is inferred from Kepler field stars. We find that only models with large convective cores, such as those calculated with our newly proposed "maximal-overshoot" scheme, can match the average $\Delta\Pi_1$ reported. However, we also find another possible solution that is related to the method used to determine $\Delta\Pi_1$: mode trapping can raise the observationally inferred $\Delta\Pi_1$ well above its true value. Even after accounting for these two proposed resolutions to the discrepancy in average $\Delta\Pi_1$, models still predict more CHeB stars with low $\Delta\Pi_1$ ($ < 270$ s) than are observed. We establish two possible remedies for this: i) there may be a difficulty in determining $\Delta\Pi_1$ for early CHeB stars (when $\Delta\Pi_1$ is lowest) because of the effect that the sharp composition profile at the hydrogen burning shell has on the pulsations, or ii) the mass of the helium core at the flash is higher than predicted. Our conclusions highlight the need for the reporting of selection effects in asteroseismic population studies in order to safely use this information to constrain stellar evolution theory., Comment: 24 pages. 24 figures. Published in MNRAS

We have serendipitously identified the first lithium-rich giant star located close to the red giant branch bump in a globular cluster. Through intermediate-resolution FLAMES spectra we derived a lithium abundance of A(Li)=2.55 (assuming local thermodynamical equilibrium), which is extremely high considering the star's evolutionary stage. Kinematic and photometric analysis confirm the object as a member of the globular cluster NGC 362. This is the fourth Li-rich giant discovered in a globular cluster but the only one known to exist at a luminosity close to the bump magnitude. The three previous detections are clearly more evolved, located close to, or beyond the tip of their red giant branch. Our observations are able to discard the accretion of planets/brown dwarfs, as well as an enhanced mass-loss mechanism as a formation channel for this rare object. Whilst the star sits just above the cluster bump luminosity, its temperature places it towards the blue side of the giant branch in the colour-magnitude diagram. We require further dedicated observations to unambiguously identify the star as a red giant: we are currently unable to confirm whether Li production has occurred at the bump of the luminosity function or if the star is on the pre zero-age horizontal branch. The latter scenario provides the opportunity for the star to have synthesised Li rapidly during the core helium flash or gradually during its red giant branch ascent via some extra mixing process., Comment: Accepted for publication in The Astrophysical Journal Letters

We explore the final fates of massive intermediate-mass stars by computing detailed stellar models from the zero age main sequence until near the end of the thermally pulsing phase. These super-AGB and massive AGB star models are in the mass range between 5.0 and 10.0 Msun for metallicities spanning the range Z=0.02-0.0001. We probe the mass limits M_up, M_n and M_mass, the minimum masses for the onset of carbon burning, the formation of a neutron star, and the iron core-collapse supernovae respectively, to constrain the white dwarf/electron-capture supernova boundary. We provide a theoretical initial to final mass relation for the massive and ultra-massive white dwarfs and specify the mass range for the occurrence of hybrid CO(Ne) white dwarfs. We predict electron-capture supernova (EC-SN) rates for lower metallicities which are significantly lower than existing values from parametric studies in the literature. We conclude the EC-SN channel (for single stars and with the critical assumption being the choice of mass-loss rate) is very narrow in initial mass, at most approximately 0.2 Msun. This implies that between ~ 2-5 per cent of all gravitational collapse supernova are EC-SNe in the metallicity range Z=0.02 to 0.0001. With our choice for mass-loss prescription and computed core growth rates we find, within our metallicity range, that CO cores cannot grow sufficiently massive to undergo a Type 1.5 SN explosion., Comment: 15 pages, 7 figures, accepted for publication in MNRAS