Source author record

Simon W. Campbell

Simon W. Campbell 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

17works
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

17 published item(s)

preprint2022arXiv

Integrated Mass Loss of Evolved Stars in M4 using Asteroseismology

Mass loss remains a major uncertainty in stellar modelling. In low-mass stars, mass loss is most significant on the red giant branch (RGB), and will impact the star's evolutionary path and final stellar remnant. Directly measuring the mass difference of stars in various phases of evolution represents one of the best ways to quantify integrated mass loss. Globular clusters (GCs) are ideal objects for this. M4 is currently the only GC for which asteroseismic data exists for stars in multiple phases of evolution. Using K2 photometry, we report asteroseismic masses for 75 red giants in M4, the largest seismic sample in a GC to date. We find an integrated RGB mass loss of $Δ\bar{M} = 0.17 \pm 0.01 ~\mathrm{M}_{\odot}$, equivalent to a Reimers' mass-loss coefficient of $η_R = 0.39$. Our results for initial mass, horizontal branch mass, $η_R$, and integrated RGB mass loss show remarkable agreement with previous studies, but with higher precision using asteroseismology. We also report the first detections of solar-like oscillations in early asymptotic giant branch (EAGB) stars in GCs. We find an average mass of $\bar{M}_{\text{EAGB}}=0.54 \pm 0.01 ~\mathrm{M}_{\odot}$, significantly lower than predicted by models. This suggests larger-than-expected mass loss on the horizontal branch. Alternatively, it could indicate unknown systematics in seismic scaling relations for the EAGB. We discover a tentative mass bi-modality in the RGB sample, possibly due to the multiple populations. In our red horizontal branch sample, we find a mass distribution consistent with a single value. We emphasise the importance of seismic studies of GCs since they could potentially resolve major uncertainties in stellar theory.

preprint2020arXiv

Discovery of ubiquitous lithium production in low-mass stars

The vast majority of stars with mass similar to the Sun are expected to only destroy lithium over the course of their lives, via low-temperature nuclear burning. This has now been supported by observations of hundreds of thousands of red giant stars (Brown et al. 1989, Kumar et al. 2011, Deepak et al. 2019, Singh et al. 2019, Casey et al. 2019). Here we perform the first large-scale systematic investigation into the Li content of stars in the red clump phase of evolution, which directly follows the red giant branch phase. Surprisingly we find that all red clump stars have high levels of lithium for their evolutionary stage. On average the lithium content increases by a factor of 40 after the end of the red giant branch stage. This suggests that all low-mass stars undergo a lithium production phase between the tip of the red giant branch and the red clump. We demonstrate that our finding is not predicted by stellar theory, revealing a stark tension between observations and models. We also show that the heavily studied (Brown et al. 1989, Reddy et al. 2005, Kumar et al. 2011, Singh et al. 2019, Casey et al. 2019) very Li-rich giants, with A(Li) $> +1.5$ dex, represent only the extreme tail of the lithium enhancement distribution, comprising 3% of red clump stars. Our findings suggest a new definition limit for Li-richness in red clump stars, A(Li) $> -0.9$ dex, which is much lower than the limit of A(Li) $> +1.5$ dex used over many decades (Brown et al. 1989, Castilho et al. 1995, Reddy et al. 2005, Carlberg et al. 2016, Casey et al. 2019, Holanda et al. 2020).

preprint2020arXiv

Neutron-capture elements in dwarf galaxies II: Challenges for the s- and i-processes at low metallicity

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.

preprint2020arXiv

The GALAH Survey: A new constraint on cosmological lithium and Galactic lithium evolution from warm dwarf stars

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.

preprint2018arXiv

Fates of the oldest intermediate-mass stars - Primordial to Extremely Metal-Poor AGB and Super-AGB Stars: White Dwarf or Supernova progenitors

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.

preprint2015arXiv

Beyond Mixing-length Theory: a step toward 321D

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éclet 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.

preprint2015arXiv

The treatment of mixing in core helium burning models - I. Implications for asteroseismology

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 $ΔΠ_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 $ΔΠ_1$ reported. However, we also find another possible solution that is related to the method used to determine $ΔΠ_1$: mode trapping can raise the observationally inferred $ΔΠ_1$ well above its true value. Even after accounting for these two proposed resolutions to the discrepancy in average $ΔΠ_1$, models still predict more CHeB stars with low $ΔΠ_1$ ($ < 270$ s) than are observed. We establish two possible remedies for this: i) there may be a difficulty in determining $ΔΠ_1$ for early CHeB stars (when $ΔΠ_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.

preprint2015arXiv

The treatment of mixing in core helium burning models -- II. Constraints from cluster star counts

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 $Δ\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.

preprint2014arXiv

Super and massive AGB stars - III. Nucleosynthesis in metal-poor and very metal-poor stars - Z=0.001 and 0.0001

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.

preprint2013arXiv

On the internal pollution mechanisms in the globular cluster NGC 6121 (M4): heavy-element abundances and AGB models

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.

preprint2013arXiv

Rubidium abundances in the globular clusters NGC 6752, NGC 1904 and NGC 104 (47 Tuc)

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.

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.

preprint2013arXiv

Transition of the Stellar Initial Mass Function Explored with Binary Population Synthesis

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.

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.

preprint2010arXiv

Evolution and nucleosynthesis of extremely metal-poor and metal-free low- and intermediate-mass stars II. s-process nucleosynthesis during the core He flash

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.

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

The core helium flash revisited III. From Pop I to Pop III stars

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.