The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approach using detached eclipsing binaries. We then move to more indirect and model-dependent methods, such as the quite commonly used isochrone or stellar track fitting. The arrival of quantitative asteroseismology has opened a completely new approach to determine stellar masses and to complement and improve the accuracy of other methods. We include methods for different evolutionary stages, from the pre-main sequence to evolved (super)giants and final remnants. For all methods uncertainties and restrictions will be discussed. We provide lists of altogether more than 200 benchmark stars with relative mass accuracies between $[0.3,2]\%$ for the covered mass range of $M\in [0.1,16]\,M_\odot$, $75\%$ of which are stars burning hydrogen in their core and the other $25\%$ covering all other evolved stages. We close with a recommendation how to combine various methods to arrive at a "mass-ladder" for stars.

Transport of angular momentum in stellar interiors is currently not well understood. Asteroseismology can provide us with estimates of internal rotation of stars and thereby advances our understanding of angular momentum transport. We can measure core-rotation rates in red-giant stars and we can place upper bounds on surface-rotation rates using measurements of dipole ($l=1$) modes. Here, we aim to determine the theoretical sensitivity of modes of different spherical degree towards the surface rotation. Additionally, we aim to identify modes that can potentially add sensitivity at intermediate radii. We used asteroseismic rotational inversions to probe the internal stellar rotation profiles in red-giant models from the base of the red-giant branch up to the luminosity bump. We used the inversion method of multiplicative optimally localised averages (MOLA) to assess how well internal and surface rotation rates can be recovered from different mode sets and different synthetic rotation profiles. We confirm that dipole mixed modes are sufficient to set constraints on the average core-rotation rates in red giants. However, surface-rotation rates estimated with only dipole mixed modes are contaminated by the core rotation. We show that the sensitivity to the surface rotation decreases from the base of the red-giant branch until it reaches a minimum at 0.6-0.8$L_\text{bump}$ due to a glitch in the buoyancy frequency. Thereafter a narrow range of increased surface sensitivity just below the bump luminosity exists. Quadrupole and octopole modes have more sensitivity in the outer parts of the star. If observed, quadrupole and octopole modes enable us to distinguish between differential and solid body rotation in the convection zone. To obtain accurate estimates of rotation rates at intermediate radii, acoustic oscillation modes with a spherical degree of $l\approx10$ are needed.

Convective boundary mixing (CBM) is ubiquitous in stellar evolution. It is a necessary ingredient in the models in order to match observational constraints from clusters, binaries and single stars alike. We compute `effective overshoot' measures that reflect the extent of mixing and which can differ significantly from the input overshoot values set in the stellar evolution codes. We use constraints from pressure modes to infer the CBM properties of Kepler and CoRoT main-sequence and subgiant oscillators, as well as in two radial velocity targets (Procyon A and $\alpha$ Cen A). Collectively these targets allow us to identify how measurement precision, stellar spectral type, and overshoot implementation impact the asteroseismic solution. With these new measures we find that the `effective overshoot' for most stars is in line with physical expectations and calibrations from binaries and clusters. However, two F-stars in the CoRoT field (HD 49933 and HD 181906) still necessitate high overshoot in the models. Due to short mode lifetimes, mode identification can be difficult in these stars. We demonstrate that an incongruence between the radial and non-radial modes drives the asteroseismic solution to extreme structures with highly-efficient CBM as an inevitable outcome. Understanding the cause of seemingly anomalous physics for such stars is vital for inferring accurate stellar parameters from TESS data with comparable timeseries length.

Context. The large quantity of high-quality asteroseismic data that obtained from space-based photometric missions and the accuracy of the resulting frequencies motivate a careful consideration of the accuracy of computed oscillation frequencies of stellar models, when applied as diagnostics of the model properties. Aims. Based on models of red-giant stars that have been independently calculated using different stellar evolution codes, we investigate the extent to which the differences in the model calculation affect the model oscillation frequencies. Methods. For each of the models, which cover four different masses and different evolution stages on the red-giant branch, we computed full sets of low-degree oscillation frequencies using a single pulsation code and, from these frequencies, typical asteroseismic diagnostics. In addition, we carried out preliminary analyses to relate differences in the oscillation properties to the corresponding model differences. Results. In general, the differences in asteroseismic properties between the different models greatly exceed the observational precision of these properties, in particular for the nonradial modes whose mixed acoustic and gravity-wave character makes them sensitive to the structure of the deep stellar interior. In some cases, identifying these differences led to improvements in the final models presented here and in Paper I; here we illustrate particular examples of this. Conclusions. Further improvements in stellar modelling are required in order fully to utilise the observational accuracy to probe intrinsic limitations in the modelling. However, our analysis of the frequency differences and their relation to stellar internal properties provides a striking illustration of the potential of the mixed modes of red-giant stars for the diagnostics of stellar interiors.

Since the onset of the `space revolution' of high-precision high-cadence photometry, asteroseismology has been demonstrated as a powerful tool for informing Galactic archaeology investigations. The launch of the NASA TESS mission has enabled seismic-based inferences to go full sky -- providing a clear advantage for large ensemble studies of the different Milky Way components. Here we demonstrate its potential for investigating the Galaxy by carrying out the first asteroseismic ensemble study of red giant stars observed by TESS. We use a sample of 25 stars for which we measure their global asteroseimic observables and estimate their fundamental stellar properties, such as radius, mass, and age. Significant improvements are seen in the uncertainties of our estimates when combining seismic observables from TESS with astrometric measurements from the Gaia mission compared to when the seismology and astrometry are applied separately. Specifically, when combined we show that stellar radii can be determined to a precision of a few percent, masses to 5-10% and ages to the 20% level. This is comparable to the precision typically obtained using end-of-mission Kepler data

The search for twins of the Sun and Earth relies on accurate characterization of stellar and exoplanetary parameters: i.e., ages, masses, and radii. In the modern era of asteroseismology, parameters of solar-like stars are derived by fitting theoretical models to observational data, which include measurements of their oscillation frequencies, metallicity [Fe/H], and effective temperature Teff. Combining this information with transit data furthermore yields the corresponding parameters for their exoplanets. While [Fe/H] and Teff are commonly stated to a precision of ~0.1 dex and ~100 K, the impact of errors in their measurement has not been studied in practice within the context of the parameters derived from them. Here we use the Stellar Parameters in an Instant (SPI) pipeline to estimate the parameters of nearly 100 stars observed by Kepler and Gaia, many of which are confirmed planet hosts. We adjust the reported spectroscopic measurements of these stars by introducing faux systematic errors and artificially increasing the reported uncertainties, and quantify the differences in the resulting parameters. We find that a systematic error of 0.1 dex in [Fe/H] translates to differences of only 4%, 2%, and 1% on average in the resulting stellar ages, masses, and radii, which are well within their uncertainties (~11%, 3.5%, 1.4%) as derived by SPI. We also find that increasing the uncertainty of [Fe/H] measurements by 0.1 dex increases the uncertainties by only 0.01 Gyr, 0.02 M_sun, and 0.01 R_sun, which are again well below their reported uncertainties (0.5 Gyr, 0.04 M_sun, 0.02 R_sun). The results for Teff at 100 K are similar. Stellar parameters from SPI are unchanged within uncertainties by errors of up to 0.14 dex or 175 K, and are even more robust to errors in Teff than the seismic scaling relations. Consequently, the parameters for their exoplanets are robust as well.

The evolution of low-mass stars into red giants is still poorly understood. During this evolution the core of the star contracts and, simultaneously, the envelope expands -- a process known as the `mirror'. Additionally, there is a short phase where the trend for increasing luminosity is reversed. This is known as the red-giant-branch bump. We explore the underlying physical reasons for these two phenomena by considering the specific entropy distribution in the star and its temporal changes. We find that between the luminosity maximum and luminosity minimum of the bump there is no mirror present and the star is fully contracting. The contraction is halted and the star regains its mirror when the hydrogen-burning shell reaches the mean molecular weight discontinuity. This marks the luminosity minimum of the bump.

(Abridged). We introduce the Aarhus Red Giants Challenge, a series of detailed comparisons between widely used stellar evolution and oscillation codes aiming at establishing the minimum level of uncertainties in properties of red giants arising solely from numerical implementations. Using 9 state-of-the-art stellar evolution codes, we defined a set of input physics and physical constants for our calculations and calibrated the convective efficiency to a specific point on the main sequence. We produced evolutionary tracks and stellar structure models at fixed radius along the red-giant branch for masses of 1.0 M$_\odot$, 1.5 M$_\odot$, 2.0 M$_\odot$, and 2.5 M$_\odot$, and compared the predicted stellar properties. Once models have been calibrated on the main sequence we find a residual spread in the predicted effective temperatures across all codes of ~20 K at solar radius and ~30-40 K in the RGB regardless of the considered stellar mass. The predicted ages show variations of 2-5% (increasing with stellar mass) which we track down to differences in the numerical implementation of energy generation. The luminosity of the RGB-bump shows a spread of about 10% for the considered codes, which translates into magnitude differences of ~0.1 mag in the optical V-band. We also compare the predicted [C/N] abundance ratio and found a spread of 0.1 dex or more for all considered masses. Our comparisons show that differences at the level of a few percent still remain in evolutionary calculations of red giants branch stars despite the use of the same input physics. These are mostly due to differences in the energy generation routines and interpolation across opacities, and call for further investigations on these matters in the context of using properties of red giants as benchmarks for astrophysical studies.

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.

Stellar models utilising one-dimensional (1D), heuristic theories of convection fail to adequately describe the energy transport in superadiabatic layers. The improper modelling leads to well-known discrepancies between observed and predicted oscillation frequencies for stars with convective envelopes. Recently, three-dimensional (3D) hydrodynamic simulations of stellar envelopes have been shown to facilitate a realistic depiction of superadiabatic convection in 1D stellar models. The resulting structural changes of the boundary layers have been demonstrated to impact not only the predicted oscillation spectra but evolution tracks as well. In this paper, we quantify the consequences that the change in boundary conditions has for stellar parameter estimates of main-sequence stars. For this purpose, we investigate two benchmark stars, Alpha Centauri A and B, using Bayesian inference. We show that the improved treatment of turbulent convection makes the obtained 1D stellar structures nearly insensitive to the mixing length parameter. By using 3D simulations in 1D stellar models, we hence overcome the degeneracy between the mixing length parameter and other stellar parameters. By lifting this degeneracy, the inclusion of 3D simulations has the potential to yield more robust parameter estimates. In this way, a more realistic depiction of superadiabatic convection has important implications for any field that relies on stellar models, including the study of the chemical evolution of the Milky Way Galaxy and exoplanet research.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

1D stellar evolution codes employ rudimentary treatments of turbulent convection. For stars with convective envelopes, this leads to systematic errors in the predicted oscillation frequencies needed for asteroseismology. One way of mending these structural inadequacies is through patching, whereby the outermost layers of 1D models are replaced by the mean stratifications from 3D simulations. In order to viably implement this approach in asteroseismic analysis, interpolation throughout precomputed 3D envelopes is required. We present a method that interpolates throughout precomputed 3D envelopes as a function of effective temperature, surface gravity, and metallicity. We conduct a series of validation tests that demonstrate that the scheme reliably and accurately reproduces the structures of stellar envelopes and apply our method to the Sun as well as two stars observed by Kepler. We parameterize the frequency shift that results from patching and show that the functional forms are evolutionary dependent. In addition we find that neglecting modal effects, such as non-adiabatic energetics, introduces systematic errors in asteroseimically obtained stellar parameters. Both these results suggest that a cautious approach is necessary when utilizing empirical surface corrections in lieu of patching models. Our results have important implications, particularly for characterizing exoplanet systems, where accuracy is of utmost concern.

Oscillation modes in fast-rotating stars can be split into several subclasses, each with their own properties. To date, seismology of these stars cannot rely on regular pattern analysis and scaling relations. However, recently there has been the promising discovery of large separations observed in spectra of fast-rotating $\delta$ Scuti stars: they were attributed to the island-mode subclass, and linked to the stellar mean density through a scaling law. In this work, we investigate the relevance of this scaling relation by computing models of fast-rotating stars and their oscillation spectra. In order to sort the thousands of oscillation modes thus obtained, we train a convolutional neural network isolating the island modes with 96\% accuracy. Arguing that the observed large separation is systematically smaller than the asymptotic one, we retrieve the observational $\Delta\nu - \overline{\rho}$ scaling law. This relation will be used to drive forward modelling efforts, and is a first step towards mode identification and inversions for fast-rotating stars.

The internal structures and properties of oscillating red-giant stars can be accurately inferred through their global oscillation modes (asteroseismology). Based on 1460 days of Kepler observations we perform a thorough asteroseismic study to probe the stellar parameters and evolutionary stages of three red giants in eclipsing binary systems. We present the first detailed analysis of individual oscillation modes of the red-giant components of KIC 8410637, KIC5640750 and KIC9540226. We obtain estimates of their asteroseismic masses, radii, mean densities and logarithmic surface gravities by using the asteroseismic scaling relations as well as grid-based modelling. As these red giants are in double-lined eclipsing binaries, it is possible to derive their independent dynamical masses and radii from the orbital solution and compare it with the seismically inferred values. For KIC 5640750 we compute the first spectroscopic orbit based on both components of this system. We use high-resolution spectroscopic data and light curves of the three systems to determine up-to-date values of the dynamical stellar parameters. With our comprehensive set of stellar parameters we explore consistencies between binary analysis and asteroseismic methods, and test the reliability of the well-known scaling relations. For the three red giants under study, we find agreement between dynamical and asteroseismic stellar parameters in cases where the asteroseismic methods account for metallicity, temperature and mass dependence as well as surface effects. We are able to attain agreement from the scaling laws in all three systems if we use $\Delta\nu_{\rm ref,emp} = 130.8 \pm 0.9\,\mu$Hz instead of the usual solar reference value.

The wealth of asteroseismic data for red-giant stars and the precision with which these data have been observed over the last decade calls for investigations to further understand the internal structures of these stars. The aim of this work is to validate a method to measure the underlying period spacing, coupling term and mode offset of pure gravity modes that are present in the deep interiors of red-giant stars. We subsequently investigate the physical conditions of the evanescent zone between the gravity mode cavity and the pressure mode cavity. We implement an alternative mathematical description to analyse observational data and to extract the underlying physical parameters that determine the frequencies of mixed modes. This description takes the radial order of the modes explicitly into account, which reduces its sensitivity to aliases. Additionally, and for the first time, this method allows us to constrain the gravity mode offset for red-giant stars. We determine the period spacing and the coupling term for the dipole modes within a few percent of literature values. Additionally, we find that the gravity mode offset varies on a star by star basis and should not be kept fixed in the analysis. Furthermore, we find that the coupling factor is logarithmically related to the physical width of the evanescent region normalised by the radius at which the evanescent zone is located. Finally, the local density contrast at the edge of the core of red giant branch models shows a tentative correlation with the offset. (abstract abriged)

The masses and radii of solar-like oscillators can be estimated through the asteroseismic scaling relations. These relations provide a direct link between observables, i.e. effective temperature and characteristics of the oscillation spectra, and stellar properties, i.e. mean density and surface gravity (thus mass and radius). These scaling relations are commonly used to characterize large samples of stars. Usually, the Sun is used as a reference from which the structure is scaled. However, for stars that do not have a similar structure as the Sun, using the Sun as a reference introduces systematic errors as large as 10\% in mass and 5\% in radius. Several alternatives for the reference of the scaling relation involving the large frequency separation (typical frequency difference between modes of the same degree and consecutive radial order) have been suggested in the literature. In a previous paper, we presented a reference function with a dependence on both effective temperature and metallicity. The accuracy of predicted masses and radii improved considerably when using reference values calculated from our reference function. However, the residuals indicated that stars on the red-giant branch possess a mass dependence that was not accounted for. Here, we present a reference function for the scaling relation involving the large frequency separation that includes the mass dependence. This new reference function improves the derived masses and radii significantly by removing the systematic differences and mitigates the trend with $\nu_{\rm max}$ (frequency of maximum oscillation power) that exists when using the solar value as a reference.

With the advent of dedicated photometric space missions, the ability to rapidly process huge catalogues of stars has become paramount. Bellinger and Angelou et al. (2016) recently introduced a new method based on machine learning for inferring the stellar parameters of main-sequence stars exhibiting solar-like oscillations. The method makes precise predictions that are consistent with other methods, but with the advantages of being able to explore many more parameters while costing practically no time. Here we apply the method to 52 so-called "LEGACY" main-sequence stars observed by the Kepler space mission. For each star, we present estimates and uncertainties of mass, age, radius, luminosity, core hydrogen abundance, surface helium abundance, surface gravity, initial helium abundance, and initial metallicity as well as estimates of their evolutionary model parameters of mixing length, overshooting coefficient, and diffusion multiplication factor. We obtain median uncertainties in stellar age, mass, and radius of 14.8%, 3.6%, and 1.7%, respectively. The source code for all analyses and for all figures appearing in this manuscript can be found electronically at: https://github.com/earlbellinger/asteroseismology

Astronomy is in an era where all-sky surveys are mapping the Galaxy. The plethora of photometric, spectroscopic, asteroseismic and astrometric data allows us to characterise the comprising stars in detail. Here we quantify to what extent precise stellar observations reveal information about the properties of a star, including properties that are unobserved, or even unobservable. We analyse the diagnostic potential of classical and asteroseismic observations for inferring stellar parameters such as age, mass and radius from evolutionary tracks of solar-like oscillators on the lower main sequence. We perform rank correlation tests in order to determine the capacity of each observable quantity to probe structural components of stars and infer their evolutionary histories. We also analyse the principal components of classic and asteroseismic observables to highlight the degree of redundancy present in the measured quantities and demonstrate the extent to which information of the model parameters can be extracted. We perform multiple regression using combinations of observable quantities in a grid of evolutionary simulations and appraise the predictive utility of each combination in determining the properties of stars. We identify the combinations that are useful and provide limits to where each type of observable quantity can reveal information about a star. We investigate the accuracy with which targets in the upcoming TESS and PLATO missions can be characterized. We demonstrate that the combination of observations from GAIA and PLATO will allow us to tightly constrain stellar masses, ages and radii with machine learning for the purposes of Galactic and planetary studies.

Owing to the remarkable photometric precision of space observatories like Kepler, stellar and planetary systems beyond our own are now being characterized en masse for the first time. These characterizations are pivotal for endeavors such as searching for Earth-like planets and solar twins, understanding the mechanisms that govern stellar evolution, and tracing the dynamics of our Galaxy. The volume of data that is becoming available, however, brings with it the need to process this information accurately and rapidly. While existing methods can constrain fundamental stellar parameters such as ages, masses, and radii from these observations, they require substantial computational efforts to do so. We develop a method based on machine learning for rapidly estimating fundamental parameters of main-sequence solar-like stars from classical and asteroseismic observations. We first demonstrate this method on a hare-and-hound exercise and then apply it to the Sun, 16 Cyg A & B, and 34 planet-hosting candidates that have been observed by the Kepler spacecraft. We find that our estimates and their associated uncertainties are comparable to the results of other methods, but with the additional benefit of being able to explore many more stellar parameters while using much less computation time. We furthermore use this method to present evidence for an empirical diffusion-mass relation. Our method is open source and freely available for the community to use. The source code for all analyses and for all figures appearing in this manuscript can be found electronically at https://github.com/earlbellinger/asteroseismology

We conduct a series of comparisons between spectroscopic and photometric observations of globular clusters and stellar models to examine their predictive power. Data from medium-to-high resolution spectroscopic surveys of lithium allow us to investigate first dredge-up and extra mixing in two clusters well separated in metallicity. Abundances at first dredge-up are satisfactorily reproduced but there is preliminary evidence to suggest that the models overestimate the luminosity at which the surface composition first changes in the lowest-metallicity system. Our models also begin extra mixing at luminosities that are too high, demonstrating a significant discrepancy with observations at low metallicity. We model the abundance changes during extra mixing as a thermohaline process and determine that the usual diffusive form of this mechanism cannot simultaneously reproduce both the carbon and lithium observations. Hubble Space Telescope photometry provides turnoff and bump magnitudes in a large number of globular clusters and offers the opportunity to better test stellar modelling as function of metallicity. We directly compare the predicted main-sequence turn-off and bump magnitudes as well as the distance-independent parameter $\Delta M_V ~^{\rm{MSTO}}_{\rm{bump}}$. We require 15 Gyr isochrones to match the main-sequence turn-off magnitude in some clusters and cannot match the bump in low-metallicity systems. Changes to the distance modulus, metallicity scale and bolometric corrections may impact on the direct comparisons but $\Delta M_V ~^{\rm{MSTO}}_{\rm{bump}}$, which is also underestimated from the models, can only be improved through changes to the input physics. Overshooting at the base of the convective envelope with an efficiency that is metallicity dependent is required to reproduce the empirically determined value of $\Delta M_V ~^{\rm{MSTO}}_{\rm{bump}}$.

The presence of multiple populations in globular clusters has been well established thanks to high-resolution spectroscopy. It is widely accepted that distinct populations are a consequence of different stellar generations: intra-cluster pollution episodes are required to produce the peculiar chemistry observed in almost all clusters. Unfortunately, the progenitors responsible have left an ambiguous signature and their nature remains unresolved. To constrain the candidate polluters, we have measured lithium and aluminium abundances in more than 180 giants across three systems: NGC~1904, NGC~2808, and NGC~362. The present investigation along with our previous analysis of M12 and M5 affords us the largest database of simultaneous determinations of Li and Al abundances. Our results indicate that Li production has occurred in each of the three clusters. In NGC~362 we detected an M12-like behaviour, with first and second-generation stars sharing very similar Li abundances favouring a progenitor that is able to produce Li, such as AGB stars. Multiple progenitor types are possible in NGC~1904 and NGC~2808, as they possess both an intermediate population comparable in lithium to the first generation stars and also an extreme population, that is enriched in Al but depleted in Li. A simple dilution model fails in reproducing this complex pattern. Finally, the internal Li variation seems to suggest that the production efficiency of this element is a function of the cluster's mass and metallicity - low-mass or relatively metal-rich clusters are more adept at producing Li.

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.

In recent years much interest has been shown in the process of thermohaline mixing in red giants. In low and intermediate mass stars this mechanism first activates at the position of the bump in the luminosity function, and has been identified as a likely candidate for driving the slow mixing inferred to occur in these stars. One particularly important consequence of this process, which is driven by a molecular weight inversion, is the destruction of lithium. We show that the degree of lithium destruction, or in some cases production, is extremely sensitive to the numerical details of the stellar models. Within the standard 1D diffusion approximation to thermohaline mixing, we find that dfferent evolution codes, with their default numerical schemes, can produce lithium abundances that differ from one another by many orders of magnitude. This disagreement is worse for faster mixing. We perform experiments with four independent stellar evolution codes, and derive conditions for the spatial and temporal resolution required for a converged numerical solution. The results are extremely sensitive to the timesteps used. We find that predicted lithium abundances published in the literature until now should be treated with caution.

Convergent lines of evidence suggest that globular clusters host multiple stellar populations. It appears that they experience at least two episodes of star formation whereby a fraction of first-generation stars contribute astrated ejecta to form the second generation(s). To identify the polluting progenitors we require distinguishing chemical signatures such as that provided by lithium. Theoretical models predict that lithium can be synthesised in AGB stars, whereas no net Li production is expected from other candidates. It has been shown that in order to reproduce the abundance pattern found in M4, Li production must occur within the polluters, favouring the AGB scenario. Here we present Li and Al abundances for a large sample of RGB stars in M12 and M5. These clusters have a very similar metallicity, whilst demonstrating differences in several cluster properties. Our results indicate that the first-generation and second-generation stars share the same Li content in M12; we recover an abundance pattern similar to that observed in M4. In M5 we find a higher degree of complexity and a simple dilution model fails in reproducing the majority of the stellar population. In both clusters we require Li production across the different stellar generations, but production seems to have occurred to different extents. We suggest that such a difference might be related to the cluster mass with the Li production being more efficient in less-massive clusters. This is the first time a statistically significant correlation between the Li spread within a GC and its luminosity has been demonstrated. Finally, although Li-producing polluters are required to account for the observed pattern, other mechanisms, such as MS depletion, might have played a role in contributing to the Li internal variation, though at relatively low level.

The asymptotic giant branch (AGB) phase is the final stage of nuclear burning for low-mass stars. Although Milky Way globular clusters are now known to harbour (at least) two generations of stars they still provide relatively homogeneous samples of stars that are used to constrain stellar evolution theory. It is predicted by stellar models that the majority of cluster stars with masses around the current turn-off mass (that is, the mass of the stars that are currently leaving the main sequence phase) will evolve through the AGB phase. Here we report that all of the second-generation stars in the globular cluster NGC 6752 -- 70 per cent of the cluster population -- fail to reach the AGB phase. Through spectroscopic abundance measurements, we found that every AGB star in our sample has a low sodium abundance, indicating that they are exclusively first-generation stars. This implies that many clusters cannot reliably be used for star counts to test stellar evolution timescales if the AGB population is included. We have no clear explanation for this observation.

The handful of available observations of AGB stars in Galactic Globular Clusters suggest that the GC AGB populations are dominated by cyanogen-weak stars. This contrasts strongly with the distributions in the RGB (and other) populations, which generally show a 50:50 bimodality in CN band strength. If it is true that the AGB populations show very different distributions then it presents a serious problem for low mass stellar evolution theory, since such a surface abundance change going from the RGB to AGB is not predicted by stellar models. However this is only a tentative conclusion, since it is based on very small AGB sample sizes. To test whether this problem really exists we have carried out an observational campaign specifically targeting AGB stars in GCs. We have obtained medium resolution spectra for about 250 AGB stars across 9 Galactic GCs using the multi-object spectrograph on the AAT (2df/AAOmega). We present some of the preliminary findings of the study for the second parameter trio of GCs: NGC 288, NGC 362 and NGC 1851. The results indeed show that there is a deficiency of stars with strong CN bands on the AGB. To confirm that this phenomenon is robust and not just confined to CN band strengths and their vagaries, we have made observations using FLAMES/VLT to measure elemental abundances for NGC 6752.We present some initial results from this study also. Our sodium abundance results show conclusively that only a subset of stars in GCs experience the AGB phase of evolution. This is the first direct, concrete confirmation of the phenomenon.

Observed chemical (anti)correlations in proton-capture elements among globular cluster stars are presently recognized as the signature of self-enrichment from now extinct, previous generations of stars. This defines the multiple population scenario. Since fluorine is also affected by proton captures, determining its abundance in globular clusters provides new and complementary clues regarding the nature of these previous generations, and supplies strong observational constraints to the chemical enrichment timescales. In this paper we present our results on near-infrared CRIRES spectroscopic observations of six cool giant stars in NGC 6656 (M22): the main objective is to derive the F content and its internal variation in this peculiar cluster, which exhibits significant changes in both light and heavy element abundances. We detected F variations across our sample beyond the measurement uncertainties and found that the F abundances are positively correlated with O and anticorrelated with Na, as expected according to the multiple population framework. Furthermore, our observations reveal an increase in the F content between the two different sub-groups, s-process rich and s-process poor, hosted within M22. The comparison with theoretical models suggests that asymptotic giant stars with masses between 4 and 5 Msun are responsible for the observed chemical pattern, confirming evidence from previous works: the difference in age between the two sub-components in M22 must be not larger than a few hundreds Myr.

The Galactic globular cluster NGC 1851 has raised much interest since HST photometry revealed that it hosts a double subgiant branch. Here we report on our homogeneous study into the cyanogen (CN) bandstrengths in the RGB population (17 stars) and AGB population (21 stars) using AAOmega/2dF spectra with R $\sim 3000$. We discover that NGC 1851 hosts a quadrimodal distribution of CN bandstrengths in its RGB and AGB populations. This result supports the merger formation scenario proposed for this cluster, such that the CN quadrimodality could be explained by the superposition of two `normal' bimodal populations. A small sample overlap with an abundance catalogue allowed us to tentatively explore the relationship between our CN populations and a range of elemental abundances. We found a striking correlation between CN and [O/Na]. We also found that the four CN peaks may be paired -- the two CN-weaker populations being associated with low Ba and the two CN-stronger populations with high Ba. If true then s-process abundances would be a good diagnostic for disentangling the two original clusters in the merger scenario. More observations are needed to confirm the quadrimodality, and also the relationship between the subpopulations. We also report CN results for NGC 288 as a comparison. Our relatively large samples of AGB stars show that both clusters have a bias towards CN-weak AGB populations.

It is now widely accepted that globular cluster red giant branch stars owe their strange abundance patterns to a combination of pollution from progenitor stars and in situ extra mixing. In this hybrid theory a first generation of stars imprint abundance patterns into the gas from which a second generation forms. The hybrid theory suggests that extra mixing is operating in both populations and we use the variation of [C/Fe] with luminosity to examine how efficient this mixing is. We investigate the observed red giant branches of M3, M13, M92, M15 and NGC 5466 as a means to test a theory of thermohaline mixing. The second parameter pair M3 and M13 are of intermediate metallicity and our models are able to account for the evolution of carbon along the RGB in both clusters. Although, in order to fit the most carbon-depleted main-sequence stars in M13 we require a model whose initial [C/Fe] abundance leads to a carbon abundance lower than is observed. Furthermore our results suggest that stars in M13 formed with some primary nitrogen (higher C+N+O than stars in M3). In the metal-poor regime only NGC 5466 can be tentatively explained by thermohaline mixing operating in multiple populations. We find thermohaline mixing unable to model the depletion of [C/Fe] with magnitude in M92 and M15. It appears as if extra mixing is occurring before the luminosity function bump in these clusters. To reconcile the data with the models would require first dredge-up to be deeper than found in extant models.

We review the observational evidence for extra mixing in stars on the red giant branch (RGB) and discuss why thermohaline mixing is a strong candidate mechanism. We recall the simple phenomenological description of thermohaline mixing, and aspects of mixing in stars in general. We use observations of M3 to constrain the form of the thermohaline diffusion coefficient and any associated free parameters. This is done by matching [C/Fe] and [N/Fe] along the RGB of M3. After taking into account a presumed initial primordial bimodality of [C/Fe] in the CN-weak and CN-strong stars our thermohaline mixing models can explain the full spread of [C/Fe]. Thermohaline mixing can produce a significant change in [N/Fe] as a function of absolute magnitude on the RGB for initially CN-weak stars, but not for initially CN-strong stars, which have so much nitrogen to begin with that any extra mixing does not significantly affect the surface nitrogen composition.

We provide a brief review of thermohaline physics and why it is a candidate extra mixing mechanism during the red giant branch (RGB). We discuss how thermohaline mixing (also called $\delta$ $\mu$ mixing) during the RGB due to helium-3 burning, is more complicated than the operation of thermohaline mixing in other stellar contexts (such as following accretion from a binary companion). We try to use observations of carbon depletion in globular clusters to help constrain the formalism and the diffusion coefficient or mixing velocity that should be used in stellar models. We are able to match the spread of carbon depletion for metal poor field giants but are unable to do so for cluster giants, which may show evidence of mixing prior to even the first dredge-up event.

A previously reported literature search suggested that the AGB stars in Galactic globular clusters may be showing different distributions of CN-strong and CN-weak stars as compared to their RGB stars. In most cases the second giant branches of GCs appeared to be deficient in stars with strong CN bands. However the sample sizes of AGB stars at that time were too small to give a definitive picture. Thus an observing campaign targeting GC AGB stars was proposed. We now have medium resolution spectral observations of about 250 GC AGB stars across 9 globular clusters, obtained with the 2dF/AAOmega instrument on the Anglo-Australian Telescope. In this paper we report some preliminary results regarding the distributions of CN-strong and CN-weak stars on the two giant branches of a selection of globular clusters. We find that some GCs show a total lack of CN-strong stars on the AGB, whilst some show a reduction in CN-strong stars as compared to the RGB. Standard stellar evolution does not predict this change in surface abundance between the two giant branches. We discuss some possible causes of this unexpected phenomenon.

There is an apparent dichotomy between the metal-poor ([Fe/H] <= -2) yet carbon-normal giants and their carbon-rich counterparts. The former undergo significant depletion of carbon on the red giant branch after they have undergone first dredge-up, whereas the latter do not appear to experience significant depletion. We investigate this in the context that the extra mixing occurs via the thermohaline instability that arises due to the burning of helium-3. We present the evolution of [C/Fe], [N/Fe] and carbon-12/carbon-13 for three models: a carbon-normal metal-poor star, and two stars that have accreted material from a 1.5 solar mass AGB companion, one having received 0.01 solar masses of material and the other having received 0.1 solar masses. We find the behaviour of the carbon-normal metal-poor stars is well reproduced by this mechanism. In addition, our models also show that the efficiency of carbon-depletion is significantly reduced in carbon-rich stars. This extra-mixing mechanism is able to reproduce the observed properties of both carbon-normal and carbon-rich stars.

Recent work has proposed that a merger event between a red-giant and a He white dwarf may be responsible for the production of R-stars (Izzard et al, 2007). We investigate the proposed evolution and nucleosynthesis of such a model. We simulate the hypothesized late ignition of the core flash by increasing the neutrino losses until the ignition occurs sufficiently far from the centre that the subsequent evolution produces dredge-up of carbon to the extent that the post-flash object is a carbon star. Detailed nucleosynthesis is performed within this approximation, and we show that the overall properties are broadly consistent with the observations. Details will depend on the dynamics of the merger event.