Lithium depletion and enrichment in the cosmos is not yet well understood. To help tighten constraints on stellar and Galactic evolution models, we present the largest high-resolution analysis of Li abundances A(Li) to date, with results for over 100 000 GALAH field stars spanning effective temperatures $5900\,\mathrm{K} \lesssim \rm{T_{eff}} \lesssim7000\,\mathrm{K}$ and metallicities $-3 \lesssim \rm[Fe/H] \lesssim +0.5$. We separated these stars into two groups, on the warm and cool side of the so-called Li-dip, a localised region of the Kiel diagram wherein lithium is severely depleted. We discovered that stars in these two groups show similar trends in the A(Li)-[Fe/H] plane, but with a roughly constant offset in A(Li) of 0.4 dex, the warm group having higher Li abundances. At $\rm[Fe/H]\gtrsim-0.5$, a significant increasing in Li abundance with increasing metallicity is evident in both groups, signalling the onset of significant Galactic production. At lower metallicity, stars in the cool group sit on the Spite plateau, showing a reduced lithium of around 0.4 dex relative to the primordial value predicted from Big Bang nucleosynthesis (BBN). However, stars in the warm group between [Fe/H] = -1.0 and -0.5, form an elevated plateau that is largely consistent with the BBN prediction. This may indicate that these stars in fact preserve the primordial Li produced in the early Universe.

We investigate the properties of 1306 red giant stars with high photospheric abundances of lithium observed by the GALAH, K2-HERMES and TESS-HERMES surveys, and discuss them in the context of proposed mechanisms for lithium enrichment in giant stars. We confirm that Li-rich giants are rare, making up only 1.1 per cent of our giant star sample. We use GALAH+ DR3 stellar parameters and a Bayesian isochrone analysis to divide the sample into first-ascent red giant branch and red clump stars, and confirm these classifications using asteroseismic data from K2. We find that red clump stars are more than three times as likely to be lithium-rich as red giant branch stars, and the occurrence rate of lithium richness with metallicity is quite different between the two populations. The probability for a rapidly rotating giant to be lithium-rich is distinctly higher as compared to a non rapidly rotating giant, independent of its evolutionary state. There is also a clear correspondence on the red clump between rapid rotation and extreme levels of lithium enrichment. There are almost no stars on the secondary red clump with lithium abundances above the primordial level, indicating a mass dependence in the red clump lithium enrichment mechanism(s). The complex distribution of lithium-rich giants across evolutionary phase, metallicity, rotation rate and mass implies multiple independent mechanisms for producing lithium enrichment.

The slow ($s$) and intermediate ($i$) neutron ($n$) capture processes occur both in asymptotic giant branch (AGB) stars, and in massive stars. To study the build-up of the $s$- and $i$-products at low metallicity, we investigate the abundances of Y, Ba, La, Nd, and Eu in 98 stars, at $-2.4<\text{[Fe/H]}<-0.9$, in the Sculptor dwarf spheroidal galaxy. The chemical enrichment from AGB stars becomes apparent at $\text{[Fe/H]}\approx-2$ in Sculptor, and causes [Y/Ba], [La/Ba], [Nd/Ba] and [Eu/Ba] to decrease with metallicity, reaching subsolar values at the highest $\text{[Fe/H]}\approx-1$. To investigate individual nucleosynthetic sites, we compared three $n$-rich Sculptor stars with theoretical yields. One carbon-enhanced metal-poor (CEMP-no) star with high $\text{[Sr, Y, Zr]}>+0.7$ is best fit with a model of a rapidly-rotating massive star, the second (likely CH star) with the $i$-process, while the third has no satisfactory fit. For a more general understanding of the build-up of the heavy elements, we calculate for the first time the cumulative contribution of the $s$- and $i$-processes to the chemical enrichment in Sculptor, and compare with theoretical predictions. By correcting for the $r$-process, we derive $\text{[Y/Ba]}_{s/i}=-0.85\pm0.16$, $\text{[La/Ba]}_{s/i}=-0.49\pm0.17$, and $\text{[Nd/Ba]}_{s/i}=-0.48\pm0.12$, in the overall $s$- and/or $i$-process in Sculptor. These abundance ratios are within the range of those of CEMP stars in the Milky Way, which have either $s$- or $i$-process signatures. The low $\text{[Y/Ba]}_{s/i}$ and $\text{[La/Ba]}_{s/i}$ that we measure in Sculptor are inconsistent with them arising from the $s$-process only, but are more compatible with models of the $i$-process. Thus we conclude that both the $s$- and $i$-processes were important for the build-up of $n$-capture elements in the Sculptor dwarf spheroidal galaxy.

We find an asymptotic limit for the integral scale dissipation length in strongly stratified stellar convection; this is not adjustable, and we identify it with the "mixing length parameter", $\alpha \sim 5/3$, a theoretical prediction which agrees to within 10% with the fitted empirical values quoted by eight independent stellar evolutionary groups, as well as several 3D stellar atmosphere groups. For strong stratification, the dissipation length approaches the density scale height; in weak stratification it shrinks to the depth for such thin convective regions. Simulations with different zoning give a monotonic decrease in resolution errors with added zones, down to barely detectable for the most difficult case (lower boundary layer), and an estimate of the dissipation due to the turbulent cascade. No explicit viscosity is required; stellar dissipation occurs from nonlinear fluid effects (Taylor, Onsager). These implicit large eddy simulations of the Euler equations, with simultaneous analysis by Reynolds averaging (RA-ILES), resolve the energy-containing eddies, develop a turbulent cascade down to the grid scale, and produce surfaces of separation at boundary layers. Like experiments (Warhaft), and direct numerical simulations of the Navier-Stokes equations (DNS, Sreenivasan), our simulations develop "anomalous terms" for dissipation in the turbulent cascade (due to intermittency, anisotropy, and interactions between coherent structures). They have more realistic behavior at high order than the K41 theory of Kolmogorov. Mixing length theory (MLT, Boehm-Vitense) overestimates the convective enthalpy flux for strong stratification, which has consequences for the theory of luminous blue variables, for the solar metal abundance, and deep solar convection.

The Andromeda Galaxy is the closest spiral galaxy to us and has been the subject of numerous studies. It harbors a massive dark matter (DM) halo which may span up to ~600 kpc across and comprises ~90% of the galaxy's total mass. This halo size translates into a large diameter of 42 degrees on the sky for an M31-Milky Way (MW) distance of 785 kpc, but its presumably low surface brightness makes it challenging to detect with gamma-ray telescopes. Using 7.6 years of Fermi Large Area Telescope (Fermi-LAT) observations, we make a detailed study of the gamma-ray emission between 1-100 GeV towards M31's outer halo, with a total field radius of 60 degrees centered at M31, and perform an in-depth analysis of the systematic uncertainties related to the observations. We use the cosmic ray (CR) propagation code GALPROP to construct specialized interstellar emission models (IEMs) to characterize the foreground gamma-ray emission from the MW, including a self-consistent determination of the isotropic component. We find evidence for an extended excess that appears to be distinct from the conventional MW foreground, having a total radial extension upwards of ~120-200 kpc from the center of M31. We discuss plausible interpretations of the excess emission but emphasize that uncertainties in the MW foreground, and in particular, modeling of the H I-related components, have not been fully explored and may impact the results.

The Andromeda Galaxy is the closest spiral galaxy to us and has been the subject of numerous studies. It harbors a massive dark matter (DM) halo which may span up to ~600 kpc across and comprises ~90% of the galaxy's total mass. This halo size translates into a large diameter of 42 degrees on the sky for an M31-Milky Way (MW) distance of 785 kpc, but its presumably low surface brightness makes it challenging to detect with gamma-ray telescopes. Using 7.6 years of Fermi Large Area Telescope (Fermi-LAT) observations, we make a detailed study of the gamma-ray emission between 1-100 GeV towards M31's outer halo, with a total field radius of 60 degrees centered at M31, and perform an in-depth analysis of the systematic uncertainties related to the observations. We use the cosmic ray (CR) propagation code GALPROP to construct specialized interstellar emission models (IEMs) to characterize the foreground gamma-ray emission from the MW, including a self-consistent determination of the isotropic component. We find evidence for an extended excess that appears to be distinct from the conventional MW foreground, having a total radial extension upwards of ~120-200 kpc from the center of M31. We discuss plausible interpretations of the excess emission but emphasize that uncertainties in the MW foreground, and in particular, modeling of the H I-related components, have not been fully explored and may impact the results. This study was first presented in a poster at the 8th International Fermi Symposium. The article describing the full analysis is under preparation and will be published elsewhere.

Renzini (1987) wrote an influential critique of mixing-length theory (MLT) as used in stellar evolution codes, and concluded that three-dimensional (3D) fluid dynamical simulations were needed to clarify several important issues. We have critically explored the limitations of the numerical methods and conclude that they are approaching the required accuracy. Implicit large eddy simulations (ILES) automatically connect large scale turbulence to a Kolmogorov cascade below the grid scale, allowing turbulent boundary layers to remove singularities that appear in the theory. Interactions between coherent structures give multi-modal behavior, driving intermittency and fluctuations. Reynolds averaging (RA) allows us to abstract the essential features of this dynamical behavior of boundaries which are appropriate to stellar evolution, and consider how they relate static boundary conditions (Richardson, Schwarzschild or Ledoux). We clarify several questions concerning when and why MLT works, and does not work, using both analytical theory and 3D high resolution numerical simulations. The composition gradients and boundary layer structure which are produced by our simulations suggest a self-consistent approach to boundary layers, removing the need for ad hoc procedures for 'convective overshooting' and `semi-convection'. In a companion paper we quantify the adequacy of our numerical resolution, determine of the length scale of dissipation (the `mixing length') without astronomical calibration, quantify agreement with the four-fifths law of Kolmogorov for weak stratification, and extend MLT to deal with strong stratification.

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.

Several recent spectroscopic investigations have presented conflicting results on the existence of Na-rich asymptotic giant branch (AGB) stars in the Galactic globular cluster M4 (NGC6121). The studies disagree on whether or not Na-rich red giant branch (RGB) stars evolve to the AGB. For a sample of previously published HER- MES/AAT AGB and RGB stellar spectra we present a re-analysis of O, Na, and Fe abundances, and a new analysis of Mg and Al abundances; we also present CN band strengths for this sample, derived from low-resolution AAOmega spectra. Following a detailed literature comparison, we find that the AGB samples of all studies consistently show lower abundances of Na and Al, and are weaker in CN, than RGB stars in the cluster. This is similar to recent observations of AGB stars in NGC 6752 and M 62. In an attempt to explain this result, we present new theoretical stellar evolutionary models for M 4; however, these predict that all stars, including Na-rich RGB stars, evolve onto the AGB. We test the robustness of our abundance results using a variety of atmospheric models and spectroscopic methods; however, we do not find evidence that systematic modelling uncertainties can explain the apparent lack of Na- rich AGB stars in M4. We conclude that an unexplained, but robust, discordance between observations and theory remains for the AGB stars in M 4.

The turbulent burning of nuclei is a common phenomenon in the evolution of stars. Here we examine a challenging case: the merging of the neon and oxygen burning shells in a 23 M$_{\odot}$ star. A previously unknown quasi-steady state is established by the interplay between mixing, turbulent transport, and nuclear burning. The resulting stellar structure has two burning shells within a single convection zone. We find that the new neon burning layer covers an extended region of the convection zone, with the burning peak occurring substantially below where the Damk\"ohler number first becomes equal to unity. These characteristics differ from those predicted by 1D stellar evolution models of similar ingestion events. We develop the mean-field turbulence equations that govern compositional evolution, and use them to interpret our data set. An important byproduct is a means to quantify sub-grid-scale effects intrinsic to the numerical hydrodynamic scheme. For implicit large eddy simulations, the analysis method is particularly powerful because it can reveal where and how simulated flows are modified by resolution, and provide straightforward physical interpretations of the effects of dissipation or induced transport. Focusing on the mean-field composition variance equations for our analysis, we recover a Kolmogorov rate of turbulent dissipation without it being imposed, in agreement with previous results which used the turbulent kinetic energy equation.

We examine a new method for peak localization, the centroid of the data raised to some power, which we call the generalized centroid. We derive the peak localization uncertainty for the generalized centroid and compare it with the Cramer-Rao lower bound for both Gaussian and quadratic fits (with Gaussian signal and noise). We find that the centroid of squares and the Gaussian fit yield the best results in both one and two dimensions. We perform similar analysis with a Lorentz-like signal and find that the centroid of cubes and the nonlinear least squares fit provide the best results. We support our derivations with simulations, and also show simulation results when the maximum function is used to initially estimate the peak location.

Context. Classical novae are thermonuclear explosions hosted by accreting white dwarfs in stellar binary systems. Material piles up on top of the white dwarf star under mildly degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope (mostly proton-capture reactions and beta-decays) heats the material up to peak temperatures ranging from 100 to 400 MK. In these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al), are ejected into the interstellar medium. Aims. To date, most of the efforts undertaken in the modeling of classical nova outbursts have focused on the early stages of the explosion and ejection, ignoring the interaction of the ejecta, first with the accretion disk orbiting the white dwarf, and ultimately with the secondary star. Methods. A suite of three-dimensional, SPH simulations of the interaction between the nova ejecta, the accretion disk, and the stellar companion have been performed to fill this gap, aimed at testing the influence of the different parameters (i.e., mass and velocity of the ejecta, mass and geometry of the accretion disk) on the dynamical and chemical properties of the system. Results. We discuss the conditions that lead to the disruption of the accretion disk and to mass loss from the binary system. In addition, we discuss the likelihood of chemical contamination of the stellar secondary induced by the impact with the nova ejecta and its potential effect on the next nova cycle.

It has been well established that Galactic Globular clusters (GCs) harbour more than one stellar population, distinguishable by the anti-correlations of light element abundances (C-N, Na-O, and Mg-Al). These studies have been extended recently to the asymptotic giant branch (AGB). Here we investigate the AGB of NGC 6397 for the first time. We have performed an abundance analysis of high-resolution spectra of 47 RGB and 8 AGB stars, deriving Fe, Na, O, Mg and Al abundances. We find that NGC 6397 shows no evidence of a deficit in Na-rich AGB stars, as reported for some other GCs - the subpopulation ratios of the AGB and RGB in NGC 6397 are identical, within uncertainties. This agrees with expectations from stellar theory. This GC acts as a control for our earlier work on the AGB of M 4 (with contrasting results), since the same tools and methods were used.

We present an analysis of the radial angular profile of the galacto-isotropic (GI) $\gamma$-ray flux--the statistically uniform flux in circular annuli about the Galactic center. Two different approaches are used to measure the GI flux profile in 85 months of Fermi-LAT data: the BDS statistic method which identifies spatial correlations, and a new Poisson ordered-pixel method which identifies non-Poisson contributions. Both methods produce similar GI flux profiles. The GI flux profile is well-described by an existing model of bremsstrahlung, $\pi^0$ production, inverse Compton scattering, and the isotropic background. Discrepancies with data in our full-sky model are not present in the GI component, and are therefore due to mis-modeling of the non-GI emission. Dark matter annihilation constraints based solely on the observed GI profile are close to the thermal WIMP cross section below 100 GeV, for fixed models of the dark matter density profile and astrophysical $\gamma$-ray foregrounds. Refined measurements of the GI profile are expected to improve these constraints by a factor of a few.

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.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

A recent study reported a strong apparent depression of Fe I, relative to Fe II, in the AGB stars of NGC 6752. This depression is much greater than that expected from the neglect of non-local thermodynamic equilibrium effects, in particular the dominant effect of overionisation. Here we attempt to reproduce the apparent Fe discrepancy, and investigate differences in reported sodium abundances. We compare in detail the methods and results of the recent study with those of an earlier study of NGC 6752 AGB stars. Iron and sodium abundances are derived using Fe I, Fe II, and Na I lines. Various uncertainties are explored. We reproduce the large Fe I depression found by the recent study, using different observational data and computational tools. Further investigation shows that the degree of the apparent Fe I depression is strongly dependent on the adopted stellar effective temperature. To minimise uncertainties in Fe I we derive temperatures for each star individually using the infrared flux method (IRFM). We find that the $T_{\rm{eff}}$ scales used by both the previous studies are cooler, by up to 100 K; such underestimated temperatures amplify the apparent Fe I depression. Our IRFM temperatures result in negligible apparent depression, consistent with theory. We also re-derived sodium abundances and, remarkably, found them to be unaffected by the new temperature scale. [Na/H] in the AGB stars is consistent between all studies. Since Fe is constant, it follows that [Na/Fe] is also consistent between studies, apart from any systematic offsets in Fe. We recommend the use of $(V-K)$ relations for AGB stars. We plan to investigate the effect of the improved temperature scale on other elements, and re-evaluate the subpopulation distributions on the AGB, in the next paper of this series. [abridged]

Galactic Globular clusters (GCs) are now known to harbour multiple stellar populations, which are chemically distinct in many light element abundances. It is becoming increasingly clear that asymptotic giant branch (AGB) stars in GCs show different abundance distributions in light elements compared to those in the red giant branch (RGB) and other phases, skewing toward more primordial, field-star-like abundances, which we refer to as subpopulation one (SP1). As part of a larger program targeting giants in GCs, we obtained high-resolution spectra for a sample of 106 RGB and 15 AGB stars in Messier 4 (NGC 6121) using the 2dF+HERMES facility on the Anglo-Australian Telescope. In this Letter we report an extreme paucity of AGB stars with [Na/O] > -0.17 in M4, which contrasts with the RGB that has abundances up to [Na/O] =0.55. The AGB abundance distribution is consistent with all AGB stars being from SP1. This result appears to imply that all subpopulation two stars (SP2; Na-rich, O-poor) avoid the AGB phase. This is an unexpected result given M4's horizontal branch morphology -- it does not have an extended blue horizontal branch. This is the first abundance study to be performed utilising the HERMES spectrograph.

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.

Stellar models provide a vital basis for many aspects of astronomy and astrophysics. Recent advances in observational astronomy -- through asteroseismology, precision photometry, high-resolution spectroscopy, and large-scale surveys -- are placing stellar models under greater quantitative scrutiny than ever. The model limitations are being exposed and the next generation of stellar models is needed as soon as possible. The current uncertainties in the models propagate to the later phases of stellar evolution, hindering our understanding of stellar populations and chemical evolution. Here we give a brief overview of the evolution, importance, and substantial uncertainties of core helium burning stars in particular and then briefly discuss a range of methods, both theoretical and observational, that we are using to advance the modelling.

We explore modifications to the current scenario for the slow neutron capture process in asymptotic giant branch (AGB) stars to account for the Pb deficiency observed in post-AGB stars of low metallicity ([Fe/H] ~ -1.2) and low initial mass (~ 1 - 1.5 Msun) in the Large and Small Magellanic Clouds. We calculated the stellar evolution and nucleosynthesis for a 1.3 Msun star with [Fe/H]=-1.3 and tested different amounts and distributions of protons leading to the production of the main neutron source within the 13C-pocket and proton ingestion scenarios. No s-process models can fully reproduce the abundance patterns observed in the post-AGB stars. When the Pb production is lowered the abundances of the elements between Eu and Pb, such as Er, Yb, W, and Hf, are also lowered to below those observed. Neutron-capture processes with neutron densities intermediate between the s and the rapid neutron-capture processes may provide a solution to this problem and be a common occurrence in low-mass, low-metallicity AGB stars.

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.

We examine the physical basis for algorithms to replace mixing-length theory (MLT) in stellar evolutionary computations. Our 321D procedure is based on numerical solutions of the Navier-Stokes equations. These implicit large eddy simulations (ILES) are three-dimensional (3D), time-dependent, and turbulent, including the Kolmogorov cascade. We use the Reynolds-averaged Navier-Stokes (RANS) formulation to make concise the 3D simulation data, and use the 3D simulations to give closure for the RANS equations. We further analyze this data set with a simple analytical model, which is non-local and time-dependent, and which contains both MLT and the Lorenz convective roll as particular subsets of solutions. A characteristic length (the damping length) again emerges in the simulations; it is determined by an observed balance between (1) the large-scale driving, and (2) small-scale damping. The nature of mixing and convective boundaries is analyzed, including dynamic, thermal and compositional effects, and compared to a simple model. We find that (1) braking regions (boundary layers in which mixing occurs) automatically appear {\it beyond} the edges of convection as defined by the Schwarzschild criterion, (2) dynamic (non-local) terms imply a non-zero turbulent kinetic energy flux (unlike MLT), (3) the effects of composition gradients on flow can be comparable to thermal effects, and (4) convective boundaries in neutrino-cooled stages differ in nature from those in photon-cooled stages (different P\'eclet numbers). The algorithms are based upon ILES solutions to the Navier-Stokes equations, so that, unlike MLT, they do not require any calibration to astronomical systems in order to predict stellar properties. Implications for solar abundances, helioseismology, asteroseismology, nucleosynthesis yields, supernova progenitors and core collapse are indicated.

The GALAH survey is a large high-resolution spectroscopic survey using the newly commissioned HERMES spectrograph on the Anglo-Australian Telescope. The HERMES spectrograph provides high-resolution (R ~28,000) spectra in four passbands for 392 stars simultaneously over a 2 degree field of view. The goal of the survey is to unravel the formation and evolutionary history of the Milky Way, using fossil remnants of ancient star formation events which have been disrupted and are now dispersed throughout the Galaxy. Chemical tagging seeks to identify such dispersed remnants solely from their common and unique chemical signatures; these groups are unidentifiable from their spatial, photometric or kinematic properties. To carry out chemical tagging, the GALAH survey will acquire spectra for a million stars down to V~14. The HERMES spectra of FGK stars contain absorption lines from 29 elements including light proton-capture elements, alpha-elements, odd-Z elements, iron-peak elements and n-capture elements from the light and heavy s-process and the r-process. This paper describes the motivation and planned execution of the GALAH survey, and presents some results on the first-light performance of HERMES.

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.

Angular anisotropy techniques for cosmic diffuse radiation maps are powerful probes, even for quite small data sets. A popular observable is the angular power spectrum; we present a detailed study applicable to any unbinned source skymap S(n) from which N random, independent events are observed. Its exact variance, which is due to the finite statistics, depends only on S(n) and N; we also derive an unbiased estimator of the variance from the data. First-order effects agree with previous analytic estimates. Importantly, heretofore unidentified higher-order effects are found to contribute to the variance and may cause the uncertainty to be significantly larger than previous analytic estimates---potentially orders of magnitude larger. Neglect of these higher-order terms, when significant, may result in a spurious detection of the power spectrum. On the other hand, this would indicate the presence of higher-order spatial correlations, such as a large bispectrum, providing new clues about the sources. Numerical simulations are shown to support these conclusions. Applying the formalism to an ensemble of Gaussian-distributed skymaps, the noise-dominated part of the power spectrum uncertainty is significantly increased at high multipoles by the new, higher-order effects. This work is important for harmonic analyses of the distributions of diffuse high-energy gamma-rays, neutrinos, and charged cosmic rays, as well as for populations of sparse point sources such as active galactic nuclei.

In this paper, we study the prospects for determining the nature of neutrinos in the context of a supersymmetric $B-L$ extension of the standard model by using dark matter indirect detection signals and bounds on $N_{\text{eff}}$ from the cosmic microwave background data. The model contains two new dark matter candidates whose dominant annihilation channels produce more neutrinos than neutralino dark matter in the minimal supersymmetric standard model. The photon and neutrino counts may then be used to discriminate between the two models. If the dark matter comes from the B-L sector, its indirect signals and impact on the cosmic microwave background can shed light on the nature of the neutrinos. When the light neutrinos are of Majorana type, the indirect neutrino signal from the Sun and the galactic center may show a prompt neutrino box-feature, as well as an earlier cut-off in both neutrino and gamma ray energy spectra. When the light neutrinos are of Dirac type, their contribution to the effective number of neutrinos $N_{\text{eff}}$ is at a detectable level.

The extragalactic dark matter (DM) annihilation signal depends on the product of the clumping factor, <\delta^2>, and the velocity-weighted annihilation cross section, \sigma v. This "clumping factor-\sigma v" degeneracy can be broken by comparing DM annihilation signals from multiple sources. In particular, one can constrain the minimum DM halo mass, M_min, which depends on the mass of the DM particles and the kinetic decoupling temperature, by comparing observations of individual DM sources to the diffuse DM annihilation signal. We demonstrate this with careful semi-analytic treatments of the DM contribution to the diffuse Isotropic Gamma-Ray Background (IGRB), and compare it with two recent hints of DM from the Galactic Center, namely, ~130 GeV DM annihilating dominantly in the \chi\chi\ to \gamma\gamma\ channel, and (10-30) GeV DM annihilating in the \chi\chi\ to b\bar{b} or \chi\chi\ to \tau^{+}\tau^{-} channels. We show that, even in the most conservative analysis, the Fermi IGRB measurement already provides interesting sensitivity. A more detailed analysis of the IGRB, with new Fermi IGRB measurements and modeling of astrophysical backgrounds, may be able to probe values of M_min up to 1 M_sun for the 130 GeV candidate and 10^{-6} M_sun for the light DM candidates. Increasing the substructure content of halos by a reasonable amount would further improve these constraints.

We present a new grid of stellar models and nucleosynthetic yields for super-AGB stars with metallicities Z=0.001 and 0.0001, applicable for use within galactic chemical evolution models. Contrary to more metal rich stars where hot bottom burning is the main driver of the surface composition, in these lower metallicity models the effect of third dredge-up and corrosive second dredge-up also have a strong impact on the yields. These metal-poor and very metal-poor super-AGB stars create large amounts of He4, C13 and N14, as well as the heavy magnesium isotopes Mg25 and Mg26. There is a transition in yield trends at metallicity Z approximately 0.001, below which we find positive yields of C12, O16, N15, Al27 and Si28, which is not the case for higher metallicities. We explore the large uncertainties derived from wind prescriptions in super-AGB stars, finding approximately 2 orders of magnitude difference in yields of Ne22, Na23, Mg24,25,26, Al27 and our s-process proxy isotope g. We find inclusion of variable composition low temperature molecular opacities is only critical for super-AGB stars of metallicities below Z approximately 0.001. We analyze our results, and those in the literature, to address the question: Are super-AGB stars the polluters responsible for extreme population in the globular cluster NGC 2808? Our results, as well as those from previous studies, seem unable to satisfactorily match the extreme population in this globular cluster.

The vital importance of composition-dependent low-temperature opacity in low-mass (M < 3Msun) asymptotic giant branch (AGB) stellar models of metallicity Z > 0.001 has recently been demonstrated (e.g. Marigo 2002; Ventura & Marigo 2010). Its significance to more metal-poor, intermediate mass (M > 2.5Msun) models has yet to be investigated. We show that its inclusion in lower-metallicity models ([Fe/H] < -2) is essential, and that there exists no threshold metallicity below which composition-dependent molecular opacity may be neglected. We find it to be crucial in all intermediate-mass models investigated ([Fe/H] < -2 and 2.5 < M/Msun < 5), because of the evolution of the surface chemistry, including the orders of magnitude increase in the abundance of molecule-forming species. Its effect on these models mirrors that previously reported for higher-metallicity models - increase in radius, decrease in Teff, faster mass loss, shorter thermally pulsing AGB lifetime, reduced enrichment in third dredge-up products (by a factor of three to ten), and an increase in the mass limit for hot bottom burning. We show that the evolution of low-metallicity models with composition-dependent low-temperature opacity is relatively independent of initial metal abundance because its contribution to the opacity is far outweighed by changes due to dredge-up. Our results imply a significant reduction in the expected number of nitrogen-enhanced metal-poor stars, which may help explain their observed paucity. We note that these findings are partially a product of the macrophysics adopted in our models, in particular the Vassiliadis & Wood (1993) mass loss rate which is strongly dependent on radius.

Searches for dark matter annihilation signals in the diffuse gamma-ray background are promising. We present the first comprehensive study utilizing both flux and anisotropy, using the example of a spectral line signal. Besides improving sensitivity, a combined analysis will help separately determine the particle properties of dark matter and the cosmological aspects of its clumping into halo substructure. The significance of a signal in a shot-noise-dominated anisotropy analysis increases linearly with time $t$, as opposed to $\sqrt{t}$ for a flux analysis, so a flux hint might be confirmed with an anisotropy signal. A first combined line search with Fermi-LAT would provide powerful new sensitivity; one with the proposed GAMMA-400 would be dramatically better.

Our main goals are to get a deeper insight into the evolution and final fates of intermediate-mass, extremely metal-poor (EMP) stars. We also aim to investigate their C, N, and O yields. Using the Monash University Stellar Evolution code we computed and analysed the evolution of stars of metallicity Z = 10^-5 and masses between 4 and 9 M_sun, from their main sequence until the late thermally pulsing (super) asymptotic giant branch, TP-(S)AGB phase. Our model stars experience a strong C, N, and O envelope enrichment either due to the second dredge-up, the dredge-out phenomenon, or the third dredge-up early during the TP-(S)AGB phase. Their late evolution is therefore similar to that of higher metallicity objects. When using a standard prescription for the mass loss rates during the TP-(S)AGB phase, the computed stars lose most of their envelopes before their cores reach the Chandrasekhar mass, so our standard models do not predict the occurrence of SNI1/2 for Z = 10^-5 stars. However, we find that the reduction of only one order of magnitude in the mass-loss rates, which are particularly uncertain at this metallicity, would prevent the complete ejection of the envelope, allowing the stars to either explode as an SNI1/2 or become an electron-capture SN. Our calculations stop due to an instability near the base of the convective envelope that hampers further convergence and leaves remnant envelope masses between 0.25 M_sun for our 4 M_sun model and 1.5 M_sun for our 9 M_sun model. We present two sets of C, N, and O yields derived from our full calculations and computed under two different assumptions, namely, that the instability causes a practically instant loss of the remnant envelope or that the stars recover and proceed with further thermal pulses. Our results have implications for the early chemical evolution of the Universe.

Large star-to-star variations of the abundances of proton-capture elements, such as Na and O, in globular clusters (GCs) are interpreted as the effect of internal pollution resulting from the presence of multiple stellar populations. To better constrain this scenario we investigate the abundance distribution of the heavy element rubidium (Rb) in NGC 6752, NGC 1904, and NGC 104 (47 Tuc). Combining the results from our sample with those in the literature, we found that Rb exhibits no star-to-star variations, regardless the cluster metallicity, with the possible intriguing, though very uncertain, exception of the metal-rich bulge cluster NGC 6388. If no star-to-star variations will be confirmed for all GCs, it implies that the stellar source of the proton-capture element variations must not have produced significant amounts of Rb. This element is observed to be enhanced at extremely high levels in intermediate-mass AGB (IM-AGB) stars in the Magellanic Clouds (i.e., at a metallicity similar to 47 Tuc and NGC 6388). This may present a challenge to this popular candidate polluter, unless the mass range of the observed IM-AGB stars does not participate in the formation of the second-generation stars in GCs. A number of possible solutions are available to resolve this conundrum, also given that the Magellanic Clouds observations are very uncertain and may need to be revised. The fast rotating massive stars scenario would not face this potential problem as the slow mechanical winds of these stars during their main-sequence phase do not carry any Rb enhancements; however, these candidates face even bigger issues such as the production of Li and the close over-imposition with core-collapse supernova timescales. Observations of Sr, Rb, and Zr in metal-rich clusters such as NGC 6388 and NGC 6441 are sorely needed to clarify the situation.

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.

Globular clusters display significant variations in their light-element content, pointing to the existence of a second stellar generation formed from the ejecta of an earlier generation. The nature of these internal polluters is still a matter of debate: the two most popular scenarios indicate intermediate-mass AGB stars and fast rotating massive stars. Abundances determination for some key elements can help distinguish between these competitor candidates. We present in this paper Y abundances for a sample of 103 red giant branch stars in NGC 6121. Within measurement errors, we find that the [Y/Fe] is constant in this cluster contrary to a recent suggestion. For a subsample of six stars we also find [Rb/Fe] to be constant, consistent with previous studies showing no variation in other s-process elements. We also present a new set of stellar yields for intermediate-mass AGB stellar models of 5 and 6 solar masses, including heavy element s-process abundances. The uncertainties on the mass-loss rate, the mixing-length parameter, and the nuclear reaction rates have a major impact on the stellar abundances. Within the IM-AGB pollution scenario, the constant abundance of heavy elements inside the cluster requires a marginal s-process efficiency in IM-AGB stars. Such a constrain could still be satisfied by the present models assuming a stronger mass-loss rate. The uncertainties mentioned above are limiting the predictive power of intermediate-mass AGB models. For these reasons, at the moment we are not able to clearly rule out their role as main polluters of the second population stars in globular clusters.

The stellar initial mass function (IMF) plays a crucial role in determining the number of surviving stars in galaxies, the chemical composition of the interstellar medium, and the distribution of light in galaxies. A key unsolved question is whether the IMF is universal in time and space. Here we use state-of-the-art results of stellar evolution to show that the IMF of our Galaxy made a transition from an IMF dominated by massive stars to the present-day IMF at an early phase of the Galaxy formation. Updated results from stellar evolution in a wide range of metallicities have been implemented in a binary population synthesis code, and compared with the observations of carbon-enhanced metal-poor (CEMP) stars in our Galaxy. We find that applying the present-day IMF to Galactic halo stars causes serious contradictions with four observable quantities connected with the evolution of AGB stars. Furthermore, a comparison between our calculations and the observations of CEMP stars may help us to constrain the transition metallicity for the IMF which we tentatively set at [Fe/H] = -2. A novelty of the current study is the inclusion of mass loss suppression in intermediate-mass AGB stars at low-metallicity. This significantly reduces the overproduction of nitrogen-enhanced stars that was a major problem in using the high-mass star dominated IMF in previous studies. Our results also demonstrate that the use of the present day IMF for all time in chemical evolution models results in the overproduction of Type I.5 supernovae. More data on stellar abundances will help to understand how the IMF has changed and what caused such a transition.

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.

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.

We describe cosmic gamma-ray and neutrino signals of dark matter annihilation, explaining how the complementarity of these signals provides additional information that, if observable, can enlighten the particle nature of dark matter. This is discussed in the context of exploiting the separate galactic and extragalactic components of the signal, using the spherical halo model distribution of dark matter. We motivate the discussion with supersymmetric extensions of the standard model of particle physics. We consider the minimal supersymmetric standard model (MSSM) where both neutrinos and gamma-rays are produced from annihilations. We also consider a gauged B-L, baryon number minus lepton number, extension of the MSSM, where annihilation can be purely to heavy right-handed neutrinos. We compare the galactic and extragalactic components of these signals, and conclude that it is not yet clear which may dominate when looking out of the galactic plane. To answer this question, we must have an understanding of the contribution of halo substructure to the annihilation signals. We find that different theories with indistinguishable gamma-ray signals can be distinguished in the neutrino signal. Gamma-ray annihilation signals are difficult to observe from the galactic center, due to abundant astrophysical sources; but annihilation neutrinos from there would not be so hidden, if they can be observed over the atmospheric neutrinos produced by cosmic rays.

We use the 3D stellar structure code DJEHUTY to model the ingestion of protons into the intershell convection zone of a 1 solar mass asymptotic giant branch star of metallicity Z=10^-4. We have run two simulations: a low resolution one of around 300,000 zones, and a high resolution one consisting of 2,000,000 zones. Both simulations have been evolved for about 4 hours of stellar time. We observe the existence of fast, downward flowing plumes that are able to transport hydrogen into close proximity to the helium burning shell before burning takes place. The intershell in the 3D model is richer in protons than the 1D model by several orders of magnitude and so we obtain substantially higher hydrogen-burning luminosities -- over 10^8 solar luminosities in the high resolution simulation -- than are found in the 1D model. Convective velocities in these simulations are over 10 times greater than the predictions of mixing length theory, though the 3D simulations have greater energy generation due to the enhanced hydrogen burning. We find no evidence of the convective zone splitting into two, though this could be as a result of insufficient spatial resolution or because the models have not been evolved for long enough. We suggest that the 1D mixing length theory and particularly the use of a diffusion algorithm for mixing do not give an accurate picture of these events. An advective mixing scheme may give a better representation of the transport processes seen in the 3D models.

We present a formalism for estimating the angular power spectrum of extragalactic gamma-rays produced by dark matter annihilating with any general velocity-dependent cross section. The relevant density and velocity distribution of dark matter is modeled as an ensemble of smooth, universal, rigid, disjoint, spherical halos with distribution and universal properties constrained by simulation data. We apply this formalism to theories of dark matter with p-wave annihilation, for which the relative-velocity-weighted annihilation cross section is \sigma v=a+bv^2. We determine that this significantly increases the gamma-ray power if b/a >> 10^6. The effect of p-wave annihilation on the angular power spectrum is very similar for the sample of particle physics models we explored, suggesting that the important effect for a given b/a is largely determined by the cosmic dark matter distribution. If the dark matter relic from strong p-wave theories is thermally produced, the intensities of annihilation gamma-rays are strongly p-wave suppressed, making them difficult to observe. If an angular power spectrum consistent with a strong p-wave were to be observed, it would likely indicate non-thermal production of dark matter in the early Universe.

Models of primordial and hyper-metal-poor stars with masses similar to the Sun experience an ingestion of protons into the hot core during the core helium flash phase at the end of their red giant branch evolution. This produces a concurrent secondary flash powered by hydrogen burning that gives rise to further nucleosynthesis in the core. We perform post-process nucleosynthesis calculations on a one-dimensional stellar evolution calculation of a star of 1 solar mass and metallicity [Fe/H] = -6.5 that suffers a proton ingestion episode. Our network includes 320 nuclear species and 2,366 reactions and treats mixing and burning simultaneously. The mixing and burning of protons into the hot convective core leads to the production of 13C, which then burns via the 13C(alpha,n)16O reaction releasing a large number of free neutrons. During the first two years of neutron production the neutron poison 14N abundance is low, allowing the prodigious production of heavy elements such as strontium, barium, and lead via slow neutron captures (the s process). These nucleosynthetic products are later mixed to the stellar surface and ejected via stellar winds. We compare our results with observations of the hyper-metal-poor halo star HE 1327-2326, which shows a strong Sr overabundance. Our model provides the possibility of self-consistently explaining the Sr overabundance in HE 1327-2326 together with its C, N, and O overabundances (all within a factor of ~4) if the material were heavily diluted, for example, via mass transfer in a wide binary system. The model produces at least 18 times too much Ba than observed, but this may be within the large modelling uncertainties. In this scenario, binary systems of low mass must have formed in the early Universe. If true then this puts constraints on the primordial initial mass function.

We calculate the effects of velocity-dependent dark matter annihilation cross sections on the intensity of the extragalactic gamma-ray background. Our formalism does not assume a locally thermal distribution of dark matter particles in phase space, and is valid for arbitrary velocity-dependent annihilation. As concrete examples, we calculate the effects of p-wave annihilation (with the $v$-weighted cross section of $\sigma v=a+bv^2$) on the mean intensity of extragalactic gamma rays produced in cosmological dark matter halos. This velocity variation makes the shape of the energy spectrum harder, but this change in the shape is too small to see unless $b/a\agt 10^6$. While we find no such models in the parameter space of the Minimal Supersymmetric Standard Model (MSSM), we show that it is possible to find $b/a\agt 10^6$ in the extension MSSM$\otimes U(1)_{B-L}$. However, we find that the most dominant effect of the p-wave annihilation is the suppression of the amplitude of the gamma-ray background. A non-zero $b$ at the dark matter freeze-out epoch requires a smaller value of $a$ in order for the relic density constraint to be satisfied, suppressing the amplitude by a factor as low as $10^{-6}$ for a thermal relic. Non-thermal relics will have weaker amplitude suppression. As another velocity-dependent effect, we calculate the spectrum for s-wave annihilation into fermions enhanced by the attractive Sommerfeld effect. Resonances associated with this effect result in significantly enhanced intensities, with a slightly softer energy spectrum.

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.

Degenerate ignition of helium in low-mass stars at the end of the red giant branch phase leads to dynamic convection in their helium cores. One-dimensional (1D) stellar modeling of this intrinsically multi-dimensional dynamic event is likely to be inadequate. Previous hydrodynamic simulations imply that the single convection zone in the helium core of metal-rich Pop I stars grows during the flash on a dynamic timescale. This may lead to hydrogen injection into the core, and a double convection zone structure as known from one-dimensional core helium flash simulations of low-mass Pop III stars. We perform hydrodynamic simulations of the core helium flash in two and three dimensions to better constrain the nature of these events. To this end we study the hydrodynamics of convection within the helium cores of a 1.25 \Msun metal-rich Pop I star (Z=0.02), and a 0.85 \Msun metal-free Pop III star (Z=0) near the peak of the flash. These models possess single and double convection zones, respectively. We use 1D stellar models of the core helium flash computed with state-of-the-art stellar evolution codes as initial models for our multidimensional hydrodynamic study, and simulate the evolution of these models with the Riemann solver based hydrodynamics code Herakles which integrates the Euler equations coupled with source terms corresponding to gravity and nuclear burning. The hydrodynamic simulation of the Pop I model involving a single convection zone covers 27 hours of stellar evolution, while the first hydrodynamic simulations of a double convection zone, in the Pop III model, span 1.8 hours of stellar life. We find differences between the predictions of mixing length theory and our hydrodynamic simulations. The simulation of the single convection zone in the Pop I model shows a strong growth of the size of the convection zone due to turbulent entrainment. Hence we predict that for the Pop I model a hydrogen injection phase (i.e. hydrogen injection into the helium core) will commence after about 23 days, which should eventually lead to a double convection zone structure known from 1D stellar modeling of low-mass Pop III stars. Our two and three-dimensional hydrodynamic simulations of the double (Pop III) convection zone model show that the velocity field in the convection zones is different from that predicted by stellar evolutionary calculations. The simulations suggest that the double convection zone decays quickly, the flow eventually being dominated by internal gravity waves.

(Abridged) We demonstrate that the amount of extra mixing required to fit the observed low C/N and 12C/13C ratios in first giant branch (FGB) stars is also sufficient to explain the C and N abundances of Galactic AGB stars. We simulate the effect of extra mixing on the FGB by setting the composition of the envelope to that observed in low-mass FGB stars, and then evolve the models to the tip of the AGB. The inclusion of FGB extra mixing compositional changes has a strong effect on the C and N abundance in our AGB models, leading to compositions consistent with those measured in Galactic C-rich stars. The composition of the models is also consistent with C abundances measured in mainstream silicon carbide grains. While our models cover the range of C abundances measured in C stars in NGC 1846, we cannot simultaneously match the composition of the O and C-rich stars. Our models only match the O isotopic composition of K and some M, MS giants, and are not able to match the O composition of C-rich AGB stars. By increasing the 16O intershell abundance (based on observational evidence) it is possible to reproduce the observed trend of increasing 16O/18O and 16O/17O ratios with evolutionary phase. We conclude 1) if extra mixing occurs during the AGB it likely only occurs efficiently in low metallicity objects, or when the stars are heavily obscured making spectroscopic observations difficult, and 2) that the intershell compositions of AGB stars needs further investigation.

[Abridged] We calculate the structural evolution and nucleosynthesis of a grid of models covering the metallicity range: -6.5 < [Fe/H] < -3.0 (plus Z=0), and mass range: 0.85 < M < 3.0 Msun, amounting to 20 stars in total. In this paper, the first of a series describing and analysing this large data set, we present the resulting stellar yields. Many of the models experience violent nuclear burning episodes not seen at higher metallicities. We refer to these events as `Dual Flashes'. These events have also been reported by previous studies. Some of the material processed by the Dual Flashes is dredged up causing significant surface pollution with a distinct chemical composition. We also analyse the yields in terms of C and N, comparing them to the observed CEMP abundances. At the lowest metallicities ([Fe/H] < -4.0) we find the yields to contain ~1 to 2 dex too much carbon, in agreement with all previous studies. At higher metallicities ([Fe/H] = -3.0), where the observed data set is much larger, all our models produce yields with [C/Fe] values consistent with those observed in the most C-rich CEMPs. However it is only the low-mass models that undergo the Dual Shell Flash (which occurs at the start of the TPAGB) that can best reproduce the C and N observations. Normal Third Dredge-Up can not reproduce the observations because at these metallicities intermediate mass models (M > 2 Msun) suffer HBB which converts the C to N thus lowering [C/N] well below the observations, whilst if TDU were to occur in the low-mass (M < 1 Msun) models (we do not find it to occur in our models), the yields would be expected to be C-rich only, which is at odds with the `dual pollution' of C and N generally observed in the CEMPs.

[Context] Water fountain stars are very young post-AGB stars with high velocity water maser jets. They are the best objects to study the onset of bipolar jets from evolved stars due to their young dynamical ages. [Methods] We use the Arizona Radio Observatory 10m telescope to observe the CO J=2-1 line and compare the line parameters with that of masers. [Results] We report the detection of 12CO and 13COJ=2-1 lines from IRAS 16342-3814. The inferred 12CO mass loss rate is an order of magnitude lower than the infrared and OH mass loss rates, indicating a very cold and thick O-rich circumstellar envelope around the star. We also find a 12CO expansion velocity of Vexp = 46 +- 1 km/s that is too high for an AGB wind and confirm the systemic velocity of 44 +- 1 km/s. In addition we measure a very low 12CO/13CO line ratio of 1.7. [Conclusions] The first detection of CO lines has provided a new way to investigate the water fountain stars. Given the high expansion velocity of the CO gas and its relation to maser velocities, we infer that the CO emission region is co-located with the OH mainline masers in the warm base of the optical bipolar lobes, while the high velocity OH1612MHz and H2O masers are located in the side walls and at the farthest ends of the bipolar lobes, respectively. Further observations are highly desired to understand the very low 12CO/13CO line ratio.

A super-solar fluorine abundance was observed in the carbon-enhanced metal-poor (CEMP) star HE 1305+0132 ([F/Fe] = +2.90, [Fe/H] = -2.5). We propose that this observation can be explained using a binary model that involve mass transfer from an asymptotic giant branch (AGB) star companion and, based on this model, we predict F abundances in CEMP stars in general. We discuss wether F can be used to discriminate between the formation histories of most CEMP stars: via binary mass transfer or from the ejecta of fast-rotating massive stars. We compute AGB yields using different stellar evolution and nucleosynthesis codes to evaluate stellar model uncertainties. We use a simple dilution model to determine the factor by which the AGB yields should be diluted to match the abundances observed in HE 1305+0132. We further employ a binary population synthesis tool to estimate the probability of F-rich CEMP stars. The abundances observed in HE 1305+0132 can be explained if this star accreted 3-11% of the mass lost by its former AGB companion. The primary AGB star should have dredged-up at least 0.2 Msun of material from its He-rich region into the convective envelope via third dredge-up, which corresponds to AGB models of Z ~ 0.0001 and mass ~ 2 Msun. Many AGB model uncertainties, such as the treatment of convective borders and mass loss, require further investigation. We find that in the binary scenario most CEMP stars should also be FEMP stars, that is, have [F/Fe] > +1, while fast-rotating massive stars do not appear to produce fluorine. We conclude that fluorine is a signature of low-mass AGB pollution in CEMP stars, together with elements associated with the slow neutron-capture process.

We report on an investigation into stellar evolution and nucleosynthesis in the low and extremely low metallicity regime, including models of stars with a pure Big Bang composition (i.e. Z=0). The metallicity range of the extremely metal poor (EMP) models we have calculated is -6.5 < [Fe/H] < -3.0, whilst our models are in the mass range 0.85 < M < 3.0 Msun. Many of the EMP and Z=0 models experience violent evolutionary episodes not seen at higher metallicities. We refer to these events as `Dual Flashes' since they are characterised by peaks in the hydrogen and helium burning luminosities occurring at roughly the same time. Some of the material processed by these events is later dredged up by the convective envelope, causing signifcant surface pollution. These events have been reported by previous studies, so our results reaffirm their occurrence -- at least in 1D stellar models. The novelty of this study is that we have calculated the entire evolution of the Z=0 and EMP models, from the ZAMS to the end of the TPAGB, including detailed nucleosynthesis. We have also calculated the nucleosynthetic yields, which will soon be available in electronic format. Although subject to many uncertainties these are, as far as we are aware, the only yields available in this mass and metallicity range. In this paper we briefy describe some of the results in the context of abundance observations of EMP halo stars. This work formed part of SWC's PhD thesis (completed in March 2007) and a series of subsequent papers will describe the results of the study in more detail.

We numerically investigate abundance properties of the Galactic globular clusters (GCs) by adopting a new ``external pollution'' scenario. In this framework, GCs are assumed to originate in forming low-mass dwarfs embedded in dark matter subhalos at very high redshifts (z) and thus be chemically influenced by field AGB stars of the dwarfs during early GC formation processes. In this external pollution scenario, the ratio of the total mass of infalling gas to that of AGB ejecta during GC formation in a dwarf (s) and the time scale of gas infall (sigma_I) are the most important key parameters that can determine abundance properties of GCs. We mainly investigate the abundance inhomogeneity among light elements (e.g., C, N, O, Na, and Al) of stars in GCs by using the latest stellar yield models of metal-poor AGB stars with and without third dredge-up. Our principal results for the models with no third dredge-up, which are more consistent with observations, are as follows. Both [N/Fe] and [C/Fe] can be diverse among stars within a GC owing to chemical pollution from field AGB stars. [N/Fe] distributions in some GCs can clearly show bimodality whereas [C/Fe] is monomodal in most models. [N/Fe] distributions depend on s such that models with smaller $s$ (i.e., larger mass fraction of AGB ejecta used for GC formation) show the [N/Fe] bimodality more clearly. N-rich, C-poor stars in GCs also have higher He abundances owing to pollution from massive AGB stars with He-rich ejecta. The number fraction of He-rich stars (Y >0.30) is higher for the models with smaller s and shorter sigma_I for 3 < s <24 and 10^5 < sigma_I < 10^7 yr.Our models are in strong disagreement with the observed O-Na anticorrelation.

We present raytracing computations for light emitted from the surface of a rapidly-rotating neutron star in order to construct light curves for X-ray pulsars and bursters. These calculations are for realistic models of rapidly-rotating neutron stars which take into account both the correct exterior metric and the oblate shape of the star. We find that the most important effect arising from rotation comes from the oblate shape of the rotating star. We find that approximating a rotating neutron star as a sphere introduces serious errors in fitted values of the star's radius and mass if the rotation rate is very large. However, in most cases acceptable fits to the ratio M/R can be obtained with the spherical approximation.

We investigate the suggestion that there are stellar populations in some globular clusters with enhanced helium (Y from 0.28 to 0.40) compared to the primordial value. We assume that a previous generation of massive Asymptotic Giant Branch (AGB) stars have polluted the cluster. Two independent sets of AGB yields are used to follow the evolution of helium and CNO using a Salpeter initial mass function (IMF) and two top-heavy IMFs. In no case are we able to produce the postulated large Y ~ 0.35 without violating the observational constraint that the CNO content is nearly constant.

Using a chemical evolution model we investigate the intriguing suggestion that there are populations of stars in some globular clusters (e.g. NGC 2808, omega Centauri) with enhanced levels of helium (Y from about 0.28 to 0.40) compared to the majority of the population that presumably have a primordial helium abundance. We assume that a previous generation of massive low-metallicity Asymptotic Giant Branch (AGB) stars has polluted the cluster gas via a slow stellar wind. We use two independent sets of AGB yields computed from detailed models to follow the evolution of helium, carbon, nitrogen and oxygen in the cluster gas using a Salpeter initial mass function (IMF) and a number of top-heavy IMFs. In no case were we able to fit the observational constraints, Y > 0.30 and C+N+O approximately constant. Depending on the shape of the IMF and the yields, we either obtained Y approximately greater than 0.30 and large increases in C+N+O or Y < 0.30 and C+N+O approximately constant. These results suggest that either AGB stars alone are not responsible for the large helium enrichment or that any dredge-up from this generation of stars was less than predicted by standard models.

[Abridged] Norris et al. (1981) found that there was a distinct lack of cyanogen-strong (CN-strong) stars in their sample of AGB stars, as compared to their sample of RGB stars (which had roughly equal numbers of CN-normal and CNstrong stars). Similar features have been discovered in the AGB populations of other clusters. Unfortunately all of these studies suffer from low AGB star counts so the conclusions are not necessarily robust -- larger, statistically signifcant, sample sizes are needed. We outline the results of a literature search for relevant CN observations and describe our observing proposal to test the suggestion that there are substantial abundance differences between the AGB and RGB in galactic globular clusters. The literature search revealed that the AGB star counts for all studies (which are not, in general, studies about AGB stars in particular) are low, usually being < or = 10. The search also revealed that the picture may not be consistent between clusters. Although most clusters appear to have CN-weak AGBs, at least two seem to have CN-strong AGBs (M5 & 47 Tuc). To further complicate the picture, clusters often appear to have a combination of both CN-strong and CN-weak stars on their AGBs -- although one population tends to dominate. Again, all these assertions are however based on small sample sizes.We aim to increase the sample sizes by an order of magnitude using existing high quality photometry in which the AGB and RGB can be reliably separated. For the observations we will use a wide-field, low- to mid-resolution multi-object spectroscope. With the new information we hope to ascertain whether significant abundance differences really exist.

We are in the process of testing a popular theory that the observed abundance anomalies in the Globular Cluster NGC 6752 are due to `internal pollution' from intermediate mass asymptotic giant branch stars. To this end we are using a chemical evolution model incorporating custom-made stellar evolution yields calculated using a detailed stellar evolution code. By tracing the chemical evolution of the intracluster gas, which is polluted by two generations of stars, we are able to test the internal pollution scenario in which the Na- and Al-enhanced ejecta from intermediate mass stars is either accreted onto the surfaces of other stars, or goes toward forming new stars. In this paper we focus mainly on the nucleosynthetic yields of the AGB stars and discuss whether these stars are the source of the observed Na-O anticorrelation. Comparing our preliminary results with observational data suggests that the qualitative theory is not supported by this quantitative study.

The main astrophysical factories of fluorine (19F) are thought to be Type II supernovae, Wolf-Rayet stars, and the asymptotic giant branch (AGB) of intermediate mass stars. We present a model for the chemical evolution of fluorine in the Milky Way using a semi-analytic multi-zone chemical evolution model. For the first time, we demonstrate quantitatively the impact of fluorine nucleosynthesis in Wolf-Rayet and AGB stars. The inclusion of these latter two fluorine production sites provides a possible solution to the long-standing discrepancy between model predictions and the fluorine abundances observed in Milky Way giants. Finally, fluorine is discussed as a possible probe of the role of supernovae and intermediate mass stars in the chemical evolution history of the globular cluster omega Centauri.

A self-consistent model of the chemical evolution of the globular cluster NGC 6752 is presented to test a popular theory that observed abundance anomalies are due to ``internal pollution'' from intermediate mass asymptotic giant branch stars. We simulated the chemical evolution of the intracluster medium under the assumption that the products of Type II SNe are completely expelled from the globular cluster, whereas material ejected from stars with m < 7 M_sun is retained, due to their weak stellar winds. By tracing the chemical evolution of the intracluster gas, we tested an internal pollution scenario in which the Na- and Al-enhanced ejecta from intermediate mass stars is either accreted onto the surfaces of other stars, or goes toward forming new stars. The observed spread in Na and Al was reproduced, but not the O-Na and Mg-Al anticorrelations. In particular, neither O nor Mg are sufficiently depleted to account for the observations. We predict that the Mg content of Na-rich cluster stars should be overwhelmingly dominated by the 25,26Mg isotopes, whereas the latest data shows only a mild 26Mg enhancement and no correlation with 25Mg. Furthermore, stars bearing the imprint of intermediate mass stellar ejecta are predicted to be strongly enhanced in both C and N, in conflict with the empirical data. We find that while standard AGB stellar models do show the hot H burning that seems required to explain the observations, this is accompanied by He burning, producing primary C, N, Mg and Na (via HBB) which do not match the observations. (Abridged)

We present preliminary results of stellar structure and nucleosynthesis calculations for some early stars. The study (still in progress) seeks to explore the expected chemical signatures of second generation low- and intermediate-mass stars that may have formed out of a combination of Big Bang and Pop III (Z=0) supernovae material. Although the study is incomplete at this stage, we find some important features in our models. The initial chemical composition of these early stars is found to be significantly different to that given by just scaling the solar composition. The most notable difference is the lack of nitrogen. This should not affect the structural evolution significantly as nitrogen will be quickly produced through the CNO cycle due to the presence of carbon (and oxygen). It should however effect the nucleosynthetic yields. We also find that our very low metallicity five solar-mass model, with [Fe/H]=-4.01, does not reach the RGB - it goes directly to the helium burning phase. It does not experience the first dredge-up either. This is not a new finding but it will have an effect on the surface chemical evolution of the stars and should alter the nucleosynthetic yields that we are currently calculating. Our higher metallicity stars, with a globular cluster composition at [Fe/H]= -1.40, do experience all the standard phases of evolution but also have significantly higher surface temperatures and luminosities compared to solar metallicity stars. Their internal temperatures are also higher which should again effect the final nucleosynthetic yields.