On the Numerical Treatment and Dependence of Thermohaline Mixing in Red Giants

Abstract

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

arXiv comment19 pages, 18 figures, accepted for publication by MNRAS; revised version has corrected Fig 17

References in corpus (27)

• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars
• Approximate input physics for stellar modelling
• A Bitter Pill: The Primordial Lithium Problem Worsens
• Thermohaline mixing: A physical mechanism governing the photospheric composition of low-mass giants
• A consistent explanation for ${}^{12}$C/${}^{13}$C, ${}^7$Li, and ${}^3$He anomalies in red giant stars
• Thermohaline instability and rotation-induced mixing. I - Low- and intermediate-mass solar metallicity stars up to the end of the AGB
• Deep Mixing of He-3: Reconciling Big Bang and Stellar Nucleosynthesis
• Numerically determined transport laws for fingering ("thermohaline") convection in astrophysics
• Modelling the binary progenitor of supernova 1993J
• Chemical Transport and Spontaneous Layer Formation in Fingering Convection in Astrophysics
• Numerical Simulations of Thermohaline Convection: Implications for Extra-Mixing in Low-Mass RGB Stars
• Thermohaline mixing in evolved low-mass stars
• Thermohaline Mixing: Does It Really Govern the Atmospheric Chemical Composition of Low-Mass Red Giants?
• The depletion of carbon by extra mixing in metal-poor giants
• The effects of thermohaline mixing on low-metallicity asymptotic giant branch stars
• Lithium in the Intermediate-Age Open Cluster, NGC 3680
• Deep Mixing in Evolved Stars. II. Interpreting Li Abundances in RGB and AGB Stars
• The thermohaline, Richardson, Rayleigh-Taylor, Solberg-Hoiland, and GSF criteria in rotating stars
• The Role of Thermohaline Mixing in Intermediate- and Low-Metallicity Globular Clusters
• The solar, exoplanet and cosmological lithium problems
• Thermohaline Mixing and its Role in the Evolution of Carbon and Nitrogen Abundances in Globular Cluster Red Giants: The Test Case of Messier 3
• Lithium abundances in red giants of M4: evidence for asymptotic giant branch star pollution in globular clusters?
• The effect of the 19F(alpha,p)22Ne reaction rate uncertainty on the yield of fluorine from Wolf-Rayet Stars
• Effects of thermohaline instability and rotation-induced mixing on the evolution of light elements in the Galaxy : D, 3He and 4He
• Does simultaneous solution matter for stellar evolution codes?
• Double-diffusive mixing in stellar interiors in the presence of horizontal gradients

Cited by in corpus (48)

• Stellar yields from metal-rich asymptotic giant branch models
• NuGrid stellar data set I. Stellar yields from H to Bi for stars with metallicities Z = 0.02 and Z = 0.01
• On the sensitivity of gravito-inertial modes to differential rotation in intermediate-mass main-sequence stars
• Beyond Mixing-length Theory: a step toward 321D
• Tidal interactions between binary stars drives lithium production in low-mass red giants
• Lithium-Rich Giants in Globular Clusters
• On Carbon Burning in Super Asymptotic Giant Branch Stars
• Discovery of ubiquitous lithium production in low-mass stars
• The Gaia-ESO survey: impact of extra-mixing on C- and N-abundances of giant stars
• Abundances in the Local Region I: G and K Giants
• Uncertainties in Galactic Chemical Evolution Models
• BACCHUS Analysis of Weak Lines in APOGEE Spectra (BAWLAS)
• Lithium Enrichment Signatures of Planetary Engulfment Events in Evolved Stars
• ${}^7$Li evolution in the thin and thick disks of the Milky Way
• Light Elements in the Universe
• The AMBRE project: a study of Li evolution in the Galactic thin and thick discs
• The puzzle of the CNO isotope ratios in AGB carbon stars
• Models of red giants in the CoRoT asteroseismology fields combining asteroseismic and spectroscopic constraints
• Diagnostics of Stellar Modelling from Spectroscopy and Photometry of Globular Clusters
• Enhancement of Lithium in Red Clump Stars by the Additional Energy Loss Induced by New Physics
• Observing the metal-poor solar neighbourhood: a comparison of galactic chemical evolution predictions
• A phenomenological modification of thermohaline mixing in globular cluster red giants
• Multiple populations in massive star clusters under the magnifying glass of photometry: Theory and tools
• The Gaia-ESO survey: Lithium abundances in open cluster Red Clump stars
• The extent of mixing in stellar interiors: the open clusters Collinder 261 and Melotte 66
• 12C/13C isotopic ratios in red-giant stars of the open cluster NGC 6791
• WIYN Open Cluster Study: The Old Open Cluster, NGC 188, and a Re-evaluation of Lithium-Richness Among Red Giants
• On the serendipitous discovery of a Li-rich giant in the globular cluster NGC 362
• Is Thermohaline Mixing the Full Story? Evidence for Separate Mixing Events near the Red Giant Branch Bump
• The most metal-rich stars in the universe: chemical contributions of low and intermediate mass asymptotic giant branch stars with metallicities between $0.04\le Z\le \mathrm{0.10.}$
• Li-rich K giants, dust excess, and binarity
• Characterizing Observed Extra Mixing Trends in Red Giants using the Reduced Density Ratio from Thermohaline Models
• Bridging the Gap between Intermediate and Massive Stars I: Validation of MESA against the State-of-the-Art Monash Stellar Evolution Program for a 2$M_{\odot }$ AGB Star
• Fast and Slow Crystallization-driven Convection in White Dwarfs
• Chemical composition of evolved stars in the young open clusters NGC 4609 and NGC 5316
• Complex Organic and Inorganic Compounds in Shells of Lithium-rich K Giant Stars
• The Gaia-ESO Survey: Inhibited extra mixing in two giants of the open cluster Trumpler 20?
• Lithium in red giants: the roles of the He-core flash and the luminosity bump
• Li-rich and super Li-rich giants produced by element diffusion
• Chemical composition of giant stars in the open cluster IC 4756
• An unsupervised method for identifying $X$-enriched stars directly from spectra: Li in LAMOST
• Thermohaline mixing in extremely metal-poor stars
• Nucleosynthetic yields of intermediate-mass primordial to extremely metal-poor stars
• Convective Mixing: The Formation Channel of Li-rich Giants
• First Star Survivors as Metal-Rich Halo Stars that Experienced Supernova Explosions in Binary Systems
• Carbon Abundance of Globular Cluster M22 (NGC 6656) and the Surface Carbon Depletion Rates of the Milky Way Globular Clusters
• The Impact of Neutrino Magnetic Moments on the Evolution of the Helium Flash and Lithium-Rich Red Clump Stars
• Chemical element transport in stellar evolution models