Source author record

Richard J. Stancliffe

Richard J. Stancliffe appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

Radioactive nuclei in the early Solar System: analysis of the 15 isotopes produced by core-collapse supernovae

Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 to 100 Myr can be used to probe the origin of the Solar System. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: $^{26}$Al, $^{36}$Cl, $^{41}$Ca, $^{53}$Mn, $^{60}$Fe, $^{92}$Nb, $^{97}$Tc, $^{98}$Tc, $^{107}$Pd, $^{126}$Sn, $^{129}$I, $^{135}$Cs, $^{146}$Sm, $^{182}$Hf, and $^{205}$Pb. We probe the impact of the uncertainties of the core-collapse explosion mechanism by examining a collection of 62 core-collapse models with initial masses of 15, 20, and 25M$_{\odot}$, explosion energies between 3.4$\times$10$^{50}$ and 1.8$\times$10$^{52}$ ergs and compact remnant masses between 1.5M$_{\odot}$and 4.89M$_{\odot}$. We identify the impact of both explosion energy and remnant mass on the final yields of the SLRs. Isotopes produced within the innermost regions of the star, such as $^{92}$Nb and $^{97}$Tc, are the most affected by the remnant mass, $^{92}$Nb varying by five orders of magnitude. Isotopes synthesised primarily in explosive C-burning and explosive He-burning, such as $^{60}$Fe, are most affected by explosion energies. $^{60}$Fe increases by two orders of magnitude from the lowest to the highest explosion energy in the 15M$_{\odot}$model. The final yield of each examined SLR is used to compare to literature models.

preprint2019arXiv

Learning about the intermediate neutron-capture process from lead abundances

Lead (Pb) is predominantly produced by the slow neutron-capture process (s process) in asymptotic giant branch (AGB) stars. In contrast to significantly enhanced Pb abundances predicted by low-mass, low-metallicity AGB-models, observations of Magellanic post-AGB stars show incompatibly low Pb abundances. Observations of carbon-enhanced metal-poor (CEMP) stars whose s-process enrichments are accompanied by heavy elements traditionally associated with the rapid neutron-capture process (r process) have raised the need for a neutron-capture process operating at neutron densities intermediate to the s and r process: the so-called i process. We study i-process nucleosynthesis with single-zone nuclear-network calculations. Our i-process models can explain the heavy-element abundance patterns measured in Magellanic post-AGB stars including their puzzlingly low Pb abundances. Furthermore, the heavy-element enhancements in the post-AGB and CEMP-i stars, particularly their Pb abundance, allow us to characterise the neutron densities and exposures of the i process that produced the observed abundance patterns. We find that the lower-metallicity CEMP-i stars ($\left[ \mathrm{Fe} / \mathrm{H} \right] \approx -2.5$) have heavy-element abundances best matched by models with higher neutron densities and exposures ($τ> 2.0 \, \mathrm{mbarn}^{-1}$) compared to the higher-metallicity post-AGB stars ($\left[ \mathrm{Fe} / \mathrm{H} \right] \approx -1.3$, $τ< 1.3 \, \mathrm{mbarn}^{-1}$). This offers new constraints and insights regarding the properties of i-process sites and demonstrates that the responsible process operates on time scales of the order of a few years or less.

preprint2016arXiv

Constraining the Progenitor of the Type Ia Supernova SN 2012cg

The nature of the progenitors of SNe Ia is not yet fully understood. In the single-degenerate (SD) scenario, the collision of the SN ejecta with its companion star is expected to produce detectable UV emission in the first few days after the SN explosion within certain viewing angles. It was recently found that the $B-V$ colour of the nearby SN Ia SN 2012cg at about sixteen days before the maximum $B$-band brightness was about 0.2 mag bluer than those of other normal SNe Ia, which was reported as the first evidence for excess blue light from the interaction of normal SN Ia ejecta with its companion star. In this work, we compare current observations for SN 2012cg from its pre-explosion phase to the late-time nebular phase with theoretical predictions from binary evolution and population synthesis calculations for a variety of popular progenitor scenarios. We find that a main-sequence donor or a carbon-oxygen white dwarf donor binary system is more likely to be the progenitor of SN 2012cg. However, both scenarios also predict properties which are in contradiction to the observed features of this system. Taking both theoretical and observational uncertainties into account, we suggest that it might be too early to conclude that SN 2012cg was produced from an explosion of a Chandrasekhar-mass white dwarf in the SD scenario. Future observations and improved detailed theoretical modelling are still required to place a more stringent constraint on the progenitor of SN 2012cg.

preprint2016arXiv

How plausible are the proposed formation scenarios of CEMP-r/s stars?

CEMP-$r/s$ stars are metal-poor stars with enhanced abundances of carbon and heavy elements associated with the slow ($s$-) and rapid ($r$-) neutron-capture process. It is believed that carbon and $s$-elements were accreted from the wind of an AGB primary star, a scenario that is generally accepted to explain the formation of CEMP stars that are only enhanced in $s$-elements (CEMP-$s$ stars). The origin of $r$-element-enrichment in CEMP-$r/s$ stars is debated and many formation scenarios have been put forward. We aim to determine the likelihood of the scenarios proposed to explain the formation of CEMP-$r/s$ stars. We calculate the frequency of CEMP-$r/s$ stars among CEMP-$s$ stars for a variety of scenarios, and we compare it with that determined from an observed sample of CEMP-$r/s$ stars collected from the literature. The theoretical frequency of CEMP-$r/s$ stars predicted in most scenarios underestimates the observed ratio by at least a factor of 5. If the enrichments in $s$- and $r$-elements are independent, the model ratio of CEMP-$r/s$ to CEMP-$s$ stars is about 22%, that is approximately consistent with the lowest estimate of the observed ratio. However, this model predicts that about one third of all carbon-normal stars have [Ba/Fe] and [Eu/Fe] higher than 1, and that 40% of all CEMP stars have [Ba/Eu]$\le0$. Stars with these properties are at least ten times rarer in our observed sample. The $intermediate$ or $i$-process, which is supposedly active in some circumstances during the AGB phase, could provide an explanation of the origin of CEMP-$r/s$ stars, similar to that of CEMP-$s$ stars, in the context of wind mass accretion in binary systems. Further calculations of the nucleosynthesis of the $i$-process and of the detailed evolution of late AGB stars are needed to investigate if this scenario predicts a CEMP-$r/s$ star frequency consistent with the observations.

preprint2016arXiv

Models of low-mass helium white dwarfs including gravitational settling, thermal and chemical diffusion, and rotational mixing

A large number of extremely low-mass helium white dwarfs (ELM WDs) have been discovered in recent years. The majority of them are found in close binary systems suggesting they are formed either through a common-envelope phase or via stable mass transfer in a low-mass X-ray binary (LMXB) or a cataclysmic variable (CV) system. Here, we investigate the formation of these objects through the LMXB channel with emphasis on the proto-WD evolution in environments with different metallicities. We study, for the first time, the combined effects of rotational mixing and element diffusion (e.g. gravitational settling, thermal and chemical diffusion) on the evolution of proto-WDs and on the cooling properties of the resulting WDs. We present state-of-the-art binary stellar evolution models computed with MESA for metallicities between Z=0.0002 and Z=0.02, producing WDs with masses between 0.16-0.45 M$_{\odot}$. Our results confirm that element diffusion plays a significant role in the evolution of proto-WDs that experience hydrogen shell flashes. The occurrence of these flashes produces a clear dichotomy in the cooling timescales of ELM WDs, which has important consequences e.g. for the age determination of binary millisecond pulsars. Rotational mixing is found to counteract the effect of gravitational settling in the surface layers of young, bloated ELM proto-WDs and therefore plays a key role in determining their surface chemical abundances. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime, a crucial ingredient for understanding the newly observed ELM proto-WD pulsators. The hydrogen envelope at detachment, although small compared to the total mass of the WD, contains enough angular momentum such that the spin frequency of the resulting WD on the cooling track is well above the orbital frequency.

preprint2016arXiv

The intermediate neutron-capture process and carbon-enhanced metal-poor stars

Carbon-enhanced metal-poor (CEMP) stars in the Galactic Halo display enrichments in heavy elements associated with either the s (slow) or the r (rapid) neutron-capture process (e.g., barium and europium respectively), and in some cases they display evidence of both. The abundance patterns of these CEMP-s/r stars, which show both Ba and Eu enrichment, are particularly puzzling since the s and the r processes require neutron densities that are more than ten orders of magnitude apart, and hence are thought to occur in very different stellar sites with very different physical conditions. We investigate whether the abundance patterns of CEMP-s/r stars can arise from the nucleosynthesis of the intermediate neutron-capture process (the i process), which is characterised by neutron densities between those of the s and the r processes. Using nuclear network calculations, we study neutron capture nucleosynthesis at different constant neutron densities n ranging from $10^7$ to $10^{15}$ cm$^{-3}$. With respect to the classical s process resulting from neutron densities on the lowest side of this range, neutron densities on the highest side result in abundance patterns that show an increased production of heavy s-process and r-process elements but similar abundances of the light s-process elements. Such high values of n may occur in the thermal pulses of asymptotic giant branch (AGB) stars due to proton ingestion episodes. Comparison to the surface abundances of 20 CEMP-s/r stars show that our modelled i-process abundances successfully reproduce observed abundance patterns that could not be previously explained by s-process nucleosynthesis. Because the i-process models fit the abundances of CEMP-s/r stars so well, we propose that this class should be renamed as CEMP-i.

preprint2015arXiv

Confronting uncertainties in stellar physics: calibrating convective overshooting with eclipsing binaries

As part of a larger program aimed at better quantifying the uncertainties in stellar computations, we attempt to calibrate the extent of convective overshooting in low to intermediate mass stars by means of eclipsing binary systems. We model 12 such systems, with component masses between 1.3 and 6.2 solar masses, using the detailed binary stellar evolution code STARS, producing grids of models in both metallicity and overshooting parameter. From these, we determine the best fit parameters for each of our systems. For three systems, none of our models produce a satisfactory fit. For the remaining systems, no single value for the convective overshooting parameter fits all the systems, but most of our systems can be well described with an overshooting parameter between 0.09 and 0.15, corresponding to an extension of the mixed region above the core of about 0.1-0.3 pressure scale heights. Of the nine systems where we are able to obtain a good fit, seven can be reasonably well fit with a single parameter of 0.15. We find no evidence for a trend of the extent of overshooting with either mass or metallicity, though the data set is of limited size. We repeat our calculations with a second evolution code, MESA, and we find general agreement between the two codes. For the extension of the mixed region above the convective core required by the MESA models is about 0.15-0.4 pressure scale heights. For the system EI Cep, we find that MESA gives an overshooting region that is larger than the STARS one by about 0.1 pressure scale heights for the primary, while for the secondary the difference is only 0.05 pressure scale heights.

preprint2015arXiv

Early ultraviolet signatures from the interaction of Type Ia supernova ejecta with a stellar companion

The progenitors of SNe Ia are not yet fully understood. The two leading progenitor scenarios are the single-degenerate (SD) scenario and the double-degenerate scenario. In the SD scenario, the collision of the SN Ia ejecta with its companion star is expected to produce detectable ultraviolet (UV) emission in the first few days after the SN explosion within certain viewing angles. A strong UV flash has recently been detected in an SN 2002es-like peculiar SN Ia iPTF14atg by Cao et al., which is interpreted as evidence of an early-time UV signature due to SN ejecta interacting with its companion star, supporting the SD scenario. In this paper, we present the expected luminosity distributions of early-time UV emission arising from SN Ia ejecta-companion interaction by performing binary population synthesis calculations for different progenitor systems in the SD scenario. Our theoretical predictions will be helpful for future early-time observations of SNe Ia to constrain their possible progenitors. Assuming the observed strong UV pulse of iPTF14atg was indeed produced by the SN ejecta-companion interaction, our population synthesis model suggests that the progenitor system of iPTF14atg is most likely a red-giant donor binary system, and it is unlikely to have been a main-sequence or helium-star donor system.

preprint2015arXiv

Evolution and nucleosynthesis of helium-rich asymptotic giant branch models

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

preprint2015arXiv

Pre-explosion companion stars in Type Iax supernovae

Type Iax supernovae (SNe Iax) are proposed as one new sub-class of SNe Ia since they present observational properties that are sufficiently distinct from the bulk of SNe Ia. SNe Iax are the most common of all types of peculiar SNe by both number and rate, with an estimated rate of occurrence of about 5-30% of the total SN Ia rate. However, the progenitor systems of SNe Iax are still uncertain. Analyzing pre-explosion images at SN Iax positions provides a direct way to place strong constraints on the nature of progenitor systems of SNe Iax. In this work, we predict pre-explosion properties of binary companion stars in a variety of potential progenitor systems by performing detailed binary evolution calculations with the one-dimensional stellar evolution code STARS. This will be helpful for constraining progenitor systems of SNe Iax from their pre-explosion observations. With our binary evolution calculations, it is found that the non-degenerate helium (He) companion star to both a massive C/O WD (> 1.1 solar mass) and a hybrid C/O/Ne WD can provide an explanation for the observations of SN~2012Z-S1, but the hybrid WD+He star scenario is more favorable.

preprint2015arXiv

The interaction of core-collapse supernova ejecta with a companion star

The progenitors of many CCSNe are expected to be in binary systems. After the SN explosion, the companion may suffer from mass stripping and be shock heated as a result of the impact of the SN ejecta. If the binary system is disrupted, the companion is ejected as a runaway and hypervelocity star. By performing a series of 3D hydrodynamical simulations of the collision of SN ejecta with the companion star, we investigate how CCSN explosions affect their companions. We use the BEC code to construct the detailed companion structure at the time of SN explosion. The impact of the SN blast wave on the companion is followed by means of 3D SPH simulations using the Stellar GADGET code. For main-sequence (MS) companions, we find that the amount of removed mass, impact velocity, and chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Their relationship can be approximately fitted by power laws, which is consistent with the results obtained from impact simulations of SNe~Ia. However, we find that the impact velocity is sensitive to the momentum profile of the outer SN ejecta and, in fact, may decrease with increasing ejecta mass, depending on the modeling of the ejecta. Because most companions to Ib/c CCSNe are in their MS phase at the moment of the explosion, combined with the strongly decaying impact effects with increasing binary separation, we argue that the majority of these SNe lead to inefficient mass stripping and shock heating of the companion star following the impact of the ejecta. Our simulations show that the impact effects of Ib/c SN ejecta on the structure of MS companions, and thus their long-term post-explosion evolution, is in general not dramatic. We find that at most 10% of their mass is lost, and their resulting impact velocities are less than 100 km/s.

preprint2014arXiv

Carbon-Enhanced Metal-Poor Star Frequencies in the Galaxy: Corrections for the Effect of Evolutionary Status on Carbon Abundances

We revisit the observed frequencies of Carbon-Enhanced Metal-Poor (CEMP) stars as a function of the metallicity in the Galaxy, using data from the literature with available high-resolution spectroscopy. Our analysis excludes stars exhibiting clear over-abundances of neutron-capture elements, and takes into account the expected depletion of surface carbon abundance that occurs due to CN processing on the upper red-giant branch. This allows for the recovery of the initial carbon abundance of these stars, and thus for an accurate assessment of the frequencies of carbon-enhanced stars. The correction procedure we develope is based on stellar-evolution models, and depends on the surface gravity, log g, of a given star. Our analysis indicates that, for stars with [Fe/H]<=-2.0, 20% exhibit [C/Fe]>=+0.7. This fraction increases to 43% for [Fe/H]<=-3.0 and 81% for [Fe/H]<=-4.0, which is higher than have been previously inferred without taking the carbon-abundance correction into account. These CEMP-star frequencies provide important inputs for Galactic and stellar chemical-evolution models, as they constrain the evolution of carbon at early times and the possible formation channels for the CEMP-no stars. We also have developed a public online tool with which carbon corrections using our procedure can be easily obtained.

preprint2014arXiv

Constraints on single-degenerate Chandrasekhar mass progenitors of Type Iax supernovae

SNe Iax are proposed as one new sub-class of SNe Ia. SNe Iax have been estimated to account for ~5-30% of the total SN Ia rate, and most SNe Iax have been discovered in late-type galaxies. In addition, observations constrain the progenitor systems of some SN Iax progenitors have ages of <80Myr. Recently, the weak deflagration explosions of Ch-mass CO WDs seem to provide a viable physical scenario for SNe Iax. Here, comparing theoretical predictions from BPS calculations with SN Iax observations, we put constraints on the SD Ch-mass model as a possible Iax progenitor. A long delay times of >3Gyr and low SN rates of ~3e-5/yr are found in our RG donor channel, indicating that this channel is unlikely to produce SNe Iax. We predict that the Galactic rate from the MS (He) donor channel is ~1.5e-3/yr (3e-4/yr), which is consistent with the observed Iax rate. The short delay times in the He channel support the young host environments of SNe Iax. However, the relative long delay times in the MS donor channel are less favourable for the observational constraints on the Iax progenitor ages. Finally, we set an upper limit on the pre-SN mass-loss rate at < 10e-4 M_sun/yr (for V_wind=100 km/s). The delay times in the SD Ch-mass model do not account for that most SNe Iax are located in late-type galaxies. However, at least one Iax event (SN 2008ge) is hosted by a S0 galaxy with no signs of star formation. Moreover, current X-ray observations cannot rule out the SD Ch-mass model. Taking all these into account and considering the uncertainty of the observed rate for SNe Iax, we suggest that some SNe Iax may be produced from weak deflagrations of Ch-mass CO WDs in SD progenitors, especially in the He-star channel. This is consistent with recent analysis of HST observations for the SN Iax SN2012Z. However, this SD deflagration model is still unlikely to be the most common progenitor scenario for SNe Iax.

preprint2014arXiv

Seven New Carbon-Enhanced Metal-Poor RR Lyrae Stars

We report estimated carbon-abundance ratios, [C/Fe], for seven newly-discovered carbon-enhanced metal-poor (CEMP) RR Lyrae stars. These are well-studied RRab stars that had previously been selected as CEMP candidates based on low-resolution spectra. For this pilot study, we observed eight of these CEMP RR Lyrae candidates with the Wide Field Spectrograph (WiFeS) on the ANU 2.3m telescope. Prior to this study, only two CEMP RR Lyrae stars had been discovered: TY Gru and SDSS J1707+58. We compare our abundances to new theoretical models of the evolution of low-mass stars in binary systems. These simulations evolve the secondary stars, post accretion from an AGB donor, all the way to the RR Lyrae stage. The abundances of CEMP RR Lyrae stars can be used as direct probes of the nature of the donor star, such as its mass, and the amount of material accreted onto the secondary. We find that the majority of the sample of CEMP RR Lyrae stars is consistent with AGB donor masses of around 1.5 - 2.0 solar masses and accretion masses of a few hundredths of a solar mass. Future high-resolution studies of these newly-discovered CEMP RR Lyrae stars will help disentangle the effects of the proposed mixing processes that occur in such objects.

preprint2013arXiv

Metal-Poor Stars Observed with the Magellan Telescope I. Constraints on Progenitor Mass and Metallicity of AGB Stars Undergoing s-Process Nucleosynthesis

We present a comprehensive abundance analysis of two newly-discovered carbon-enhanced metal-poor (CEMP) stars. HE2138-3336 is a s-process-rich star with [Fe/H] = -2.79, and has the highest [Pb/Fe] abundance ratio measured thus far, if NLTE corrections are included ([Pb/Fe] = +3.84). HE2258-6358, with [Fe/H] = -2.67, exhibits enrichments in both s- and r-process elements. These stars were selected from a sample of candidate metal-poor stars from the Hamburg/ESO objective-prism survey, and followed up with medium-resolution (R ~ 2,000) spectroscopy with GEMINI/GMOS. We report here on derived abundances (or limits) for a total of 34 elements in each star, based on high-resolution (R ~ 30,000) spectroscopy obtained with Magellan-Clay/MIKE. Our results are compared to predictions from new theoretical AGB nucleosynthesis models of 1.3 Mo with [Fe/H] = -2.5 and -2.8, as well as to a set of AGB models of 1.0 to 6.0 Mo at [Fe/H] = -2.3. The agreement with the model predictions suggests that the neutron-capture material in HE2138-3336 originated from mass transfer from a binary companion star that previously went through the AGB phase, whereas for HE2258-6358, an additional process has to be taken into account to explain its abundance pattern. We find that a narrow range of progenitor masses (1.0 < M(Mo) < 1.3) and metallicities (-2.8 < [Fe/H] < -2.5) yield the best agreement with our observed elemental abundance patterns.

preprint2013arXiv

Modelling the nucleosynthetic properties of carbon-enhanced metal-poor RR Lyrae stars

Certain carbon-enhanced metal-poor stars likely obtained their composition via pollution from some of the earliest generations of asymptotic giant branch stars and as such provide important clues to early Universe nucleosynthesis. Recently, Kinman et al. discovered that the highly carbon- and barium-enriched metal-poor star SDSS J1707+58 is in fact an RR Lyrae pulsator. This gives us an object in a definite evolutionary state where the effects of dilution of material during the Main Sequence are minimised owing to the object having passed through first dredge-up. We perform detailed stellar modelling of putative progenitor systems in which we accreted material from asymptotic giant branch stars in the mass range 1-2 solar masses. We investigate how the surface abundances are affected by the inclusion of mechanisms like thermohaline mixing and gravitational settling. While we are able to find a reasonable fit to the carbon and sodium abundances of SDSS J1707+58, suggesting accretion of around 0.1 solar masses from a 2 solar mass companion, the strontium and barium abundances remain problematic and this object may have experienced something other than a standard $s$ process. We have more success in fitting the abundances of the mildly carbon-enriched, metal-poor RR Lyrae pulsator TY Gru (CS 22881-071), which we suggest received 0.1 solar masses of material from a companion of around 1 solar mass.

preprint2013arXiv

Sodium content as a predictor of the advanced evolution of globular cluster stars

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.

preprint2012arXiv

The occurrence of nitrogen-enhanced metal-poor stars: implications for the initial mass function in the early Galactic halo

Most carbon-enhanced metal-poor (CEMP) stars are thought to result from past mass transfer of He-burning material from an asymptotic giant branch (AGB) star to a low-mass companion star, which we now observe as a CEMP star. Because AGB stars of intermediate mass efficiently cycle carbon into nitrogen in their envelopes, the same evolution scenario predicts the existence of a population of nitrogen-enhanced metal-poor (NEMP) stars, with [N/Fe] > 1 and [N/C] > 0.5. Such NEMP stars are rare, although their occurrence depends on metallicity: they appear to be more common at [Fe/H] < -2.8 by about a factor of 10 compared to less metal-poor stars. We analyse the observed sample of metal-poor stars with measurements of both carbon and nitrogen to derive firm constraints on the occurrence of NEMP stars as a function of metallicity. We compare these constraints to binary population synthesis calculations in which we vary the initial distributions of mass, mass ratio and binary orbital periods. We show that the observed paucity of NEMP stars at [Fe/H] > -2.8 does not allow for large modifications in the initial mass function, as have been suggested in the literature to account for the high frequency of CEMP stars. The situation at lower metallicity is less clear, and we do not currently have stellar models to perform this comparison for [Fe/H] < -2.8. However, unless intermediate-mass AGB stars behave very differently at such low metallicity, the observed NEMP frequency at [Fe/H] < -2.8 appears incompatible with the top-heavy forms of the initial mass function suggested in the literature.

preprint2012arXiv

The Role of Thermohaline Mixing in Intermediate- and Low-Metallicity Globular Clusters

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.

preprint2011arXiv

3D hydrodynamical simulations of a proton ingestion episode in a low-metallicity asymptotic giant branch star

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.

preprint2011arXiv

Close to the Dredge: Precise X-ray C and N Abundances in lambda Andromeda and its Precocious RGB Mixing Problem

Chandra LETG+HRC-S and XMM-Newton RGS spectra of H-like C and N formed in the corona of the lambda And primary star, a mildly metal-poor G8 III-IV first ascent giant that completed dredge-up ~50 Myr ago, have been used to make a precise measurement of its surface C/N ratio. We obtain the formal result [C/N]=0.03+/-0.07, which is typical of old disk giants and in agreement with standard dredge-up theory for stars of about 1 M_sun or lower. In contrast, these stars as a group, including lambda And, have 12C/13C ~< 20, which is much lower than standard model predictions. We show that the abundances of the old disk giants are consistent with models including thermohaline mixing that begins at the red giant branch luminosity function "bump". Instead, lambda And indicates that the 12C/13C anomaly can be present immediately following dredge-up, contrary to current models of extra mixing on the red giant branch.

preprint2011arXiv

The s process in asymptotic giant branch stars of low metallicity and the composition of carbon-enhanced metal-poor stars

We present models for the slow neutron-capture process (s process) in asymptotic giant branch (AGB) stars of metallicity [Fe/H]=-2.3 and masses 0.9 Msun to 6 Msun. We encountered different regimes of neutron-capture nucleosynthesis increasing in importance as the stellar mass decreases: the 22Ne(alpha,n)25Mg reaction activated during the thermal pulses, the 13C(alpha,n)16O reaction activated in radiative conditions during the interpulse periods, and the 13C(alpha,n)16O reaction activated during the thermal pulses, also as a result of mild proton ingestion episodes. The models where the 13C burns radiatively (masses ~ 2 Msun) produce an overall good match to carbon-enhanced metal-poor (CEMP) stars showing s-process enhancements (CEMP-s), except they produce too much Na and F. On the other hand, none of our models can provide a match to the composition of CEMP stars also showing rapid-process enhancements (CEMP-s/r). The models fail to reproduce the observed Eu abundances, and they also fail to reproduce the correlation between the Eu and Ba abundances. They also cannot match the ratio of heavy to light s-process elements observed in many CEMP-s/r stars, which can be more than ten times higher than in the solar system. To explain the composition of CEMP-s/r stars we need to invoke the existence of a different "s/r" neutron-capture process either with features in-between the s and the r processes, or generated by superpositions of different neutron-capture processes in the same astrophysical site or in sites linked to each other - for example, in multiple stellar systems.

preprint2010arXiv

Does GD 356 have a Terrestrial Planetary Companion?

GD 356 is unique among magnetic white dwarfs because it shows Zeeman-split Balmer lines in pure emission. The lines originate from a region of nearly uniform field strength (delta B/B is approximately 0.1) that covers 10 per cent of the stellar surface in which there is a temperature inversion. The energy source that heats the photosphere remains a mystery but it is likely to be associated with the presence of a companion. Based on current models we use archival Spitzer IRAC observations to place a new and stringent upper limit of 12 Jupiter masses for the mass of such a companion. In the light of this result and the recent discovery of a 115 min photometric period for GD 356, we exclude previous models that invoke accretion and revisit the unipolar inductor model that has been proposed for this system. In this model a highly conducting planet with a metallic core orbits the magnetic white dwarf and, as it cuts through field lines, a current is set flowing between the two bodies. This current dissipates in the photosphere of the white dwarf and causes a temperature inversion. Such a planet is unlikely to have survived the RGB/AGB phases of evolution so we argue that it may have formed from the circumstellar disc of a disrupted He or CO core during a rare merger of two white dwarfs. GD 356 would then be a white dwarf counterpart of the millisecond binary pulsar PSR 1257+12 which is known to host a planetary system.

preprint2010arXiv

Is Extra Mixing Really Needed in Asymptotic Giant Branch Stars?

(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.

preprint2010arXiv

Thermohaline Mixing and its Role in the Evolution of Carbon and Nitrogen Abundances in Globular Cluster Red Giants: The Test Case of Messier 3

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.

preprint2009arXiv

Binary Populations of Carbon-Enhanced Metal-Poor Stars

We construct binary-star population nucleosynthesis models of carbon-enhanced metal poor (CEMP) stars. We compare the CEMP to EMP (extremely metal poor) ratio of our models to the observed ratio and find it is an order of magnitude too small. Through an increase in the efficiency of third dredge-up in low-mass, low-metallicity thermally-pulsing asymptotic-giant branch (TPAGB) stars our models better match the observations.

preprint2009arXiv

Modelling the binary progenitor of supernova 1993J

We have developed a detailed stellar evolution code capable of following the simultaneous evolution of both stars in a binary system, together with their orbital properties. To demonstrate the capabilities of the code we investigate potential progenitors for the Type IIb supernova 1993J, which is believed to have been an interacting binary system prior to its primary exploding. We use our detailed binary stellar evolution code to model this system to determine the possible range of primary and secondary masses that could have produced the observed characteristics of this system, with particular reference to the secondary. Using the luminosities and temperatures for both stars (as determined by Maund et al. 2004) and the remaining mass of the hydrogen envelope of the primary at the time of explosion, we find that if mass transfer is 100 per cent efficient the observations can be reproduced by a system consisting of a 15 solar mass primary and a 14 solar mass secondary in an orbit with an initial period of 2100 days. With a mass transfer efficiency of 50 per cent, a more massive system consisting of a 17 solar mass primary and a 16 solar mass secondary in an initial orbit of 2360 days is needed. We also investigate some of the uncertainties in the evolution, including the effects of tidal interaction, convective overshooting and thermohaline mixing.

preprint2009arXiv

Population Synthesis of Binary Carbon-enhanced Metal-poor Stars

The carbon-enhanced metal-poor (CEMP) stars constitute approximately one fifth of the metal-poor ([Fe/H] ~< -2) population but their origin is not well understood. The most widely accepted formation scenario, invokes mass-transfer of carbon-rich material from a thermally-pulsing asymptotic giant branch (TPAGB) primary star to a less massive main-sequence companion which is seen today. Recent studies explore the possibility that an initial mass function biased toward intermediate-mass stars is required to reproduce the observed CEMP fraction in stars with metallicity [Fe/H] < -2.5. These models also implicitly predict a large number of nitrogen-enhanced metal-poor (NEMP) stars which is not seen. We investigate whether the observed CEMP and NEMP to extremely metal-poor (EMP) ratios can be explained without invoking a change in the initial mass function. We confirm earlier findings that with current detailed TPAGB models the large observed CEMP fraction cannot be accounted for. We find that efficient third dredge up in low-mass (less than 1.25Msun), low-metallicity stars may offer at least a partial explanation to the large observed CEMP fraction while remaining consistent with the small observed NEMP fraction.

preprint2009arXiv

The depletion of carbon by extra mixing in metal-poor giants

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.

preprint2009arXiv

The effects of thermohaline mixing on low-metallicity asymptotic giant branch stars

We examine the effects of thermohaline mixing on the composition of the envelopes of low-metallicity asymptotic giant branch (AGB) stars. We have evolved models of 1, 1.5 and 2 solar masses from the pre-main sequence to the end of the thermally pulsing asymptotic giant branch with thermohaline mixing applied throughout the simulations. In agreement with other authors, we find that thermohaline mixing substantially reduces the abundance of helium-3 on the upper part of the red giant branch in our lowest mass model. However, the small amount of helium-3 that remains is enough to drive thermohaline mixing on the AGB. We find that thermohaline mixing is most efficient in the early thermal pulses and its efficiency drops from pulse to pulse. Nitrogen is not substantially affected by the process, but we do see substantial changes in carbon-13. The carbon-12 to carbon-13 ratio is substantially lowered during the early thermal pulses but the efficacy of the process is seen to diminish rapidly. As the process stops after a few pulses, the carbon-12 to carbon-13 ratio is still able to reach values of 10^3-10^4, which is inconsistent with the values measured in carbon-enhanced metal-poor stars. We also note a surprising increase in the lithium-7 abundance, with log epsilon(Li-7) reaching values of over 2.5 in the 1.5 solar mass model. It is thus possible to get stars which are both C- and Li-rich at the same time. We compare our models to measurements of carbon and lithium in carbon-enhanced metal-poor stars which have not yet reached the giant branch. These models can simultaneously reproduced the observed C and Li abundances of carbon-enhanced metal-poor turn-off stars that are Li-rich, but the observed nitrogen abundances still cannot be matched.

preprint2009arXiv

The evolution of low-metallicity asymptotic giant branch stars and the formation of carbon-enhanced metal-poor stars

We investigate the behaviour of asymptotic giant branch (AGB) stars between metallicities Z = 10-4 and Z = 10-8 . We determine which stars undergo an episode of flash-driven mixing, where protons are ingested into the intershell convection zone, as they enter the thermally pulsing AGB phase and which undergo third dredge-up. We find that flash-driven mixing does not occur above a metallicity of Z = 10-5 for any mass of star and that stars above 2 M do not experience this phenomenon at any metallicity. We find carbon ingestion (CI), the mixing of carbon into the tail of hydrogen burning region, occurs in the mass range 2 M to around 4 M . We suggest that CI may be a weak version of the flash-driven mechanism. We also investigate the effects of convective overshooting on the behaviour of these objects. Our models struggle to explain the frequency of CEMP stars that have both significant carbon and nitrogen enhancement. Carbon can be enhanced through flash-driven mixing, CI or just third dredge up. Nitrogen can be enhanced through hot bottom burning and the occurrence of hot dredge-up also converts carbon into nitrogen. The C/N ratio may be a good indicator of the mass of the primary AGB stars.

preprint2009arXiv

Why do low-mass stars become red giants?

We revisit the problem of why stars become red giants. We modify the physics of a standard stellar evolution code in order to determine what does and what does not contribute to a star becoming a red giant. In particular, we have run tests to try to separate the effects of changes in the mean molecular weight and in the energy generation. The implications for why stars become red giants are discussed. We find that while a change in the mean molecular weight is necessary (but not sufficient) for a 1 solar mass star to become a red giant, this is not the case in a star of 5 solar masses. It therefore seems that there may be more than one way to make a giant.