We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, dust-gas mixtures, H2 chemistry, physical viscosity, external forces including numerous galactic potentials, Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

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 ℓ = 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.

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.

Context. In the interior of stars, a convectively unstable zone typically borders a zone that is stable to convection. Convective motions can penetrate the boundary between these zones, creating a layer characterized by intermittent convective mixing, and gradual erosion of the density and temperature stratification. Aims. We examine a penetration layer formed between a central radiative zone and a large convection zone in the deep interior of a young low-mass star. Using the Multidimensional Stellar Implicit Code (MUSIC) to simulate two-dimensional compressible stellar convection in a spherical geometry over long times, we produce statistics that characterize the extent and impact of convective penetration in this layer. Methods. We apply extreme value theory to the maximal extent of convective penetration at any time. We compare statistical results from simulations which treat non-local convection, throughout a large portion of the stellar radius, with simulations designed to treat local convection in a small region surrounding the penetration layer. For each of these situations, we compare simulations of different resolution, which have different velocity magnitudes. We also compare statistical results between simulations that radiate energy at a constant rate to those that allow energy to radiate from the stellar surface according to the local surface temperature. Results. Based on the frequency and depth of penetrating convective structures, we observe two distinct layers that form between the convection zone and the stable radiative zone. We show that the probability density function of the maximal depth of convective penetration at any time corresponds closely in space with the radial position where internal waves are excited. We find that the maximal penetration depth can be modeled by a Weibull distribution with a small shape parameter. Using these results, and building on established scalings for diffusion enhanced by large-scale convective motions, we propose a new form for the diffusion coefficient that may be used for one-dimensional stellar evolution calculations in the large Péclet number regime. These results should contribute to the 321D link.

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 ± 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 AGB HB = 0.455 ± 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.

We study lithium depletion in low-mass and solar-like stars as a function of time, using a new diffusion coefficient describing extra-mixing taking place at the bottom of a convective envelope. This new form is motivated by multidimensional fully compressible, time-implicit hydrodynamic simulations performed with the MUSIC code. Intermittent convective mixing at the convective boundary in a star can be modeled using extreme value theory, a statistical analysis frequently used for finance, meteorology, and environmental science. In this Letter, we implement this statistical diffusion coefficient in a one-dimensional stellar evolution code, using parameters calibrated from multi-dimensional hydrodynamic simulations of a young low-mass star. We propose a new scenario that can explain observations of the surface abundance of lithium in the Sun and in clusters covering a wide range of ages, from ∼50 Myr to ∼4 Gyr. Because it relies on our physical model of convective penetration, this scenario has a limited number of assumptions. It can explain the observed trend between rotation and depletion, based on a single additional assumption, namely, that rotation affects the mixing efficiency at the convective boundary. We suggest the existence of a threshold in stellar rotation rate above which rotation strongly prevents the vertical penetration of plumes and below which rotation has small effects. In addition to providing a possible explanation for the long-standing problem of lithium depletion in pre-main-sequence and main-sequence stars, the strength of our scenario is that its basic assumptions can be tested by future hydrodynamic simulations.

Context. Studies of stellar convection typically use a spherical shell geometry. The radial extent of the shell and the boundary conditions applied are based on the model of the star investigated. We study the impact of different two-dimensional spherical shells on compressible convection. Realistic profiles for density and temperature from an established one-dimensional stellar evolution code are used to produce a model of a large stellar convection zone representative of a young low-mass star, such as our sun at 10 6 years of age. Aims. We analyze how the radial extent of the spherical shell changes the convective dynamics that result in the deep interior of the young sun model, far from the surface. In the near-surface layers, simple smaller-scale convection develops from the profiles of temperature and density. A central radiative zone below the convection zone provides a lower boundary on the convection zone. The inclusion of either of these physically distinct layers in the spherical shell can potentially affect the characteristics of deep convection. Methods. We perform hydrodynamic implicit large-eddy simulations of compressible convection using the MUltidimensional Stellar Implicit Code (MUSIC). Because MUSIC has been designed to use realistic stellar models produced from one-dimensional stellar evolution calculations, MUSIC simulations are capable of seamlessly modeling a whole star. Simulations in two-dimensional spherical shells that have different radial extents are performed over tens or even hundreds of convective turnover times, permitting the collection of well-converged statistics. Results. To measure the impact of the spherical shell geometry and our treatment of boundaries, we evaluate basic statistics of the convective turnover time, the convective velocity, and the overshooting layer. These quantities are selected for their relevance to one-dimensional stellar evolution calculations, so that our results are focused toward studies exploiting the so-called 321D link. We find that the inclusion in the spherical shell of the boundary between the radiative and convection zones decreases the amplitude of convective velocities in the convection zone. The inclusion of near-surface layers in the spherical shell can increase the amplitude of convective velocities, although the radial structure of the velocity profile established by deep convection is unchanged. The impact from including the near-surface layers depends on the speed and structure of small-scale convection in the near-surface layers. Larger convective velocities in the convection zone result in a commensurate increase in the overshooting layer width and decrease in the convective turnover time. These results provide support for non-local aspects of convection.

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 ≤ M/M ⊙ ≤ 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 Kepler field.