Thermohaline Mixing Research Articles

Hybrid C–O–Ne white dwarfs as progenitors of Type Ia supernovae: dependence on Urca process and mixing assumptions

When carbon is ignited off-centre in a CO core of a super-asymptotic giant branch star, its burning in a convective shell tends to propagate to the centre. Whether the C flame will actually be able to reach the centre depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the centre. As a result, a C–O–Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a ‘hybrid’ WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations, the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-centre C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of Type Ia supernova (SN Ia) progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the centre. Our hybrid SN Ia progenitor models have much lower C-to-O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.

On the numerical treatment and dependence of thermohaline mixing in red giants

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.

Spectroscopic study on the beryllium abundances of red giant stars

Abstract An extensive spectroscopic study was carried out for the beryllium abundances of 200 red giants (mostly of late G- and early K-type), which were determined from the near-UV Be ii 3131.066 line based on high-dispersion spectra obtained by Subaru/High Dispersion Spectrograph, with an aim of investigating the nature of surface Be contents in these evolved giants; e.g., dependence upon stellar parameters, degree of peculiarity along with its origin and build-up timing. We found that Be is considerably deficient (to widely different degrees from star to star) in the photosphere of these evolved giants by ∼ 1–3 dex (or more) compared to the initial abundance. While the resulting Be abundances [A(Be)] appear to depend weakly upon Teff, log g, [Fe/H], M, age, and ve sin i, this may be attributed to the metallicity dependence of A(Be) coupled with the mutual correlation between these stellar parameters, since such tendencies almost disappear in the metallicity-scaled Be abundance ([Be/Fe]). By comparing the Be abundances (as well as their correlations with Li and C) to the recent theoretical predictions based on sophisticated stellar evolution calculations, we concluded that such a considerable extent/diversity of Be deficit is difficult to explain only by the standard theory of first dredge-up in the envelope of red giants, and that some extra mixing process (such as rotational or thermohaline mixing) must be responsible, which presumably starts to operate already in the main-sequence phase. This view is supported by the fact that appreciable Be depletion is seen in less evolved intermediate-mass B–A type stars near to the main sequence.

Pre-main-sequence accretion and the formation of multiple populations in globular clusters

We investigate the viability of a model in which the chemical anomalies among Globular Cluster stars are due to accretion of gas onto the protostellar discs of low mass stars. This model has been suggested as a way to reduce the large initial cluster masses required by other models for the formation of multiple stellar generations. We numerically follow the evolution of the accreting stars, and we show that the structure of the seed star does not remain fully convective for the whole duration of the accretion phase. Stellar populations showing discrete abundances of helium in the core, that seem to be present in some clusters, might be formed with this mechanism only if accretion occurs before the core of the stars become radiative (within 2-3Myr) or if a thermohaline instability is triggered, to achieve full mixing after the accretion phase ends. We also show that the lithium abundances in accreted structures may vary by orders of magnitude in equal masses obtained by accreting different masses. In addition, the same thermohaline mixing which could provide a homogeneous helium distribution down to the stellar center, would destroy any lithium surviving in the envelope, so that both helium homogeneity and lithium survival require that the accretion phase be limited to the first couple of million years of the cluster evolution. Such a short accretion phase strongly reduces the amount of processed matter available, and reintroduces the requirement of an extremely large initial mass for the protocluster.

Which physics determines the location of the mean molecular weight minimum in red giants?

Stars ascending the red giant branch develop an inversion in mean molecular weight (mu) owing to the burning of He-3 in the region immediately above their hydrogen-burning shells. This inversion may drive thermohaline mixing and thereby be responsible for the extra mixing which is observationally indicated on the red giant branch. In this paper, we investigate the physical influences that determine the mass and temperature at which the inversion in mu develops. We find that it depends most strongly on the thermal structure of the envelope - the profiles of density and temperature in the region of the star immediately above the shell - and is otherwise relatively insensitive to abundances and nuclear reaction rates. The changes in the effects of thermohaline mixing as stars proceed up the giant branch can mostly be understood in terms of their changing thermal structure, driven by their increasing core mass.

MODULES FOR EXPERIMENTS IN STELLAR ASTROPHYSICS (MESA): PLANETS, OSCILLATIONS, ROTATION, AND MASSIVE STARS

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

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.

THE C-FLAME QUENCHING BY CONVECTIVE BOUNDARY MIXING IN SUPER-AGB STARS AND THE FORMATION OF HYBRID C/O/Ne WHITE DWARFS AND SN PROGENITORS

After off-center C ignition in the cores of super-AGB stars the C flame propagates all the way down to the center, trailing behind it the C-shell convective zone, and thus building a degenerate ONe core. This standard picture is obtained in stellar evolution simulations if the bottom C-shell convection boundary is assumed to be a discontinuity associated with a strict interpretation of the Schwarzschild condition for convective instability. However, this boundary is prone to additional mixing processes, such as thermohaline convection and convective boundary mixing. Using hydrodynamic simulations we show that, contrary to previous results, thermohaline mixing is too inefficient to interfere with the C-flame propagation. However, even a small amount of convective boundary mixing removes the physical conditions required for the C-flame propagation all the way to the center. This result holds even if we allow for some turbulent heat transport in the CBM region. As a result, super AGB stars build in their interiors hybrid C-O-Ne degenerate cores composed of a relatively large CO core (M_CO ~ 0.2 M_sun) surrounded by a thick ONe zone (M_ONe ~ 0.85 M_sun) with another thin CO layer above. If exposed by mass loss, these cores will become hybrid C-O-Ne white dwarfs. Otherwise, the ignition of C-rich material in the central core, surrounded by the thick ONe zone, may trigger a thermonuclear supernova explosion. The quenching of the C-flame may have implications for the ignition mechanism of SN Ia in the double-degenerate merger scenario.

The thermohaline, Richardson, Rayleigh-Taylor, Solberg–Høiland, and GSF criteria in rotating stars

Aims. We examine the interactions of various instabilities in rotating stars, which usually are considered as independent. Methods. An analytical study of the problem is performed, account is given to radiative losses, mu-gradients and horizontal turbulence. Results. The diffusion coefficient for an ensemble of instabilities is not given by the sum of the specific coefficients for each instability, but by the solution of a general equation. We find that thermohaline mixing is possible in low-mass red giants only if the horizontal turbulence is very weak. In rotating stars the Rayleigh-Taylor and the shear instabilities need simultaneous treating. This has for consequence that rotation laws of the form 1/r^(alpha) are predicted to be unstable for alpha > 1.6568, while the usual Rayleigh criterion predicts instability only for alpha > 2. Also, the shear instabilities are somehow reduced in Main Sequence stars by the effect of the Rayleigh-Taylor criterion. Various instability criteria should be expressed differently in rotating stars than in simplified geometries.

Red clump stars of the Milky Way – laboratories of extra-mixing

In this work we present the main atmospheric parameters, carbon, nitrogen and oxygen abundances, and 12C/13C ratios determined in a sample of 28 Galactic clump stars. Abundances of carbon were studied using the C2 band at 5086.2 {\AA}. The wavelength interval 7980-8130 {\AA} with strong CN features was analysed in order to determine nitrogen abundances and 12C/13C isotope ratios. The oxygen abundances were determined from the [O I] line at 6300 {\AA}. The mean abundances of C, N and O abundances in the investigated clump stars support our previous estimations that, compared to the Sun and dwarf stars of the Galactic disc, carbon is depleted by about 0.2 dex, nitrogen is enhanced by 0.2 dex and oxygen is close to abundances in dwarfs. The 12C/13C and C/N ratios for galactic red clump stars analysed were compared to the evolutionary models of extra-mixing. The steeper drop of 12C/13C ratio in the model of thermohaline mixing by Charbonnel & Lagarde better reflects the observational data at low stellar masses than the more shallow model of cool bottom processing by Boothroyd & Sackmann. For stars of about 2 Msun masses a modelling of rotationally induced mixing should be considered with rotation of about 250 km s-1 at the time when a star was at the hydrogen-core-burning stage.

Accretion of planetary matter from debris disks around white dwarfs: the fate of planetary systems

Heavy elements are observed in the atmospheres of many DA and DB white dwarfs, and their presence is attributed to the accretion of matter coming from debris disks. Several authors have deduced accretion rates from the observed abundances, taking into account the mixing induced by the convective zones and the gravitational settling. The obtained values are different for DA and DB white dwarfs. Here we show that an important process was forgotten in all these computations: thermohaline mixing, induced by the inverse μ-gradient built during the accretion process. Taking this mixing into account leads to an increase of the derived accretion rates, specially for DA white dwarfs, and modifies the conclusions.

Thermohaline instability and rotation-induced mixing in low and intermediate mass stars: Consequences on global asteroseismic quantities

Thermohaline mixing has been recently identified as the probable dominating process that governs the photospheric composition of low-mass bright red giant stars. Here, we present the predictions of stellar models computed with the code STAREVOL including this process together with rotational mixing. We compare our theoretical predictions with recent observations, and discuss the effects of both mechanisms on asteroseismic diagnostics.

Thermohaline instability and rotation-induced mixing

The availability of asteroseismic constraints for a large sample of stars from the missions CoRoT and Kepler paves the way for various statistical studies of the seismic properties of stellar populations. In this paper, we evaluate the impact of rotation-induced mixing and thermohaline instability on the global asteroseismic parameters at different stages of the stellar evolution from the Zero Age Main Sequence to the Thermally Pulsating Asymptotic Giant Branch to distinguish stellar populations. We present a grid of stellar evolutionary models for four metallicities (Z = 0.0001, 0.002, 0.004, and 0.014) in the mass range between 0.85 to 6.0 Msun. The models are computed either with standard prescriptions or including both thermohaline convection and rotation-induced mixing. For the whole grid we provide the usual stellar parameters (luminosity, effective temperature, lifetimes, ...), together with the global seismic parameters, i.e. the large frequency separation and asymptotic relations, the frequency corresponding to the maximum oscillation power {\nu}_{max}, the maximal amplitude A_{max}, the asymptotic period spacing of g-modes, and different acoustic radii. We discuss the signature of rotation-induced mixing on the global asteroseismic quantities, that can be detected observationally. Thermohaline mixing whose effects can be identified by spectroscopic studies cannot be caracterized with the global seismic parameters studied here. But it is not excluded that individual mode frequencies or other well chosen asteroseismic quantities might help constraining this mixing.

THERMOHALINE INSTABILITIES INSIDE STARS: A SYNTHETIC STUDY INCLUDING EXTERNAL TURBULENCE AND RADIATIVE LEVITATION

We have derived a new expression for the thermohaline mixing coefficient in stars, including the effects of radiative levitation and external turbulence, by solving Boussinesq equations in a quasi-incompressible fluid with a linear approximation. It is well known that radiative levitation of individual elements can lead to their accumulation in specific stellar layers. In some cases, it can induce important effects on the stellar structure. Here we confirm that this accumulation is moderated by thermohaline convection due to the resulting inverse $\mu$-gradient. The new coefficient that we have derived shows that the effect of radiative accelerations on the thermohaline instability itself is small. This effect must however be checked in all computations. We also confirm that the presence of large horizontal turbulence can reduce or even suppress the thermohaline convection. These results are important as they concern all the cases of heavy element accumulation in stars. The computations of radiative diffusion have to be revisited including thermohaline convection and its consequences. It may be one of the basic reasons for the fact that the observed abundances are always smaller than those predicted by pure atomic diffusion. In any case, these processes have to compete with rotation-induced mixing, but this competition is more complex than previously thought due to their mutual interaction.

Effects of thermohaline instability and rotation-induced mixing on the evolution of light elements in the Galaxy: D,3He and4He

Recent studies of low- and intermediate-mass stars show that the evolution of the chemical elements in these stars is very different from that proposed by standard stellar models. Rotation-induced mixing modifies the internal chemical structure of main sequence stars, although its signatures are revealed only later in the evolution when the first dredge-up occurs. Thermohaline mixing is likely the dominating process that governs the photospheric composition of low-mass red giant branch stars and has been shown to drastically reduce the net 3He production in these stars. The predictions of these new stellar models need to be tested against galaxy evolution. In particular, the resulting evolution of the light elements D, 3He and 4He should be compared with their primordial values inferred from the Wilkinson Microwave Anisotropy Probe data and with the abundances derived from observations of different Galactic regions. We study the effects of thermohaline mixing and rotation-induced mixing on the evolution of the light elements in the Milky Way. We compute Galactic evolutionary models including new yields from stellar models computed with thermohaline instability and rotation-induced mixing. We discuss the effects of these important physical processes acting in stars on the evolution of the light elements D, 3He, and 4He in the Galaxy. Galactic chemical evolution models computed with stellar yields including thermohaline mixing and rotation fit better observations of 3He and 4He in the Galaxy than models computed with standard stellar yields. The inclusion of thermohaline mixing in stellar models provides a solution to the long-standing "3He problem" on a Galactic scale. Stellar models including rotation-induced mixing and thermohaline instability reproduce also the observations of D and 4He.

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.

Weak G-band stars on the H-R diagram: clues to the origin of the Li anomaly

Weak G-band (WGB) stars are a rare class of cool luminous stars that present a strong depletion in carbon, but also lithium abundance anomalies that have been little explored in the literature since the first discovery of these peculiar objects in the early 50's. Here we focus on the Li-rich WGB stars and report on their evolutionary status. We explore different paths to propose a tentative explanation for the lithium anomaly. Using archive data, we derive the fundamental parameters of WGB (Teff, log g, log(L/Lsun)) using Hipparcos parallaxes and recent temperature scales. From the equivalent widths of Li resonance line at 6707 {\AA}, we uniformly derive the lithium abundances and apply when possible NLTE corrections following the procedure described by Lind et al. (2009). We also compute dedicated stellar evolution models in the mass range 3.0 to 4.5 Msun, exploring the effects of rotation-induced and thermohaline mixing. These models are used to locate the WGB stars in the H-R diagram and to explore the origin of the abundance anomalies. The location of WGB stars in the H-R diagram shows that these are intermediate mass stars of masses ranging from 3.0 to 4.5 Msun located at the clump, which implies a degeneracy of their evolutionary status between subgiant/red giant branch and core helium burning phases. The atmospheres of a large proportion of WGB stars (more than 50%) exhibit lithium abundances A(Li) \geq 1.4 dex similar to Li-rich K giants. The position of WGB stars along with the Li-rich K giants in the H-R diagram however indicates that both are well separated groups. The combined and tentatively consistent analysis of the abundance pattern for lithium, carbon and nitrogen of WGB stars seems to indicate that carbon underabundance could be decorrelated from the lithium and nitrogen overabundances.

METAL-RICH ACCRETION AND THERMOHALINE INSTABILITIES IN EXOPLANET-HOST STARS: CONSEQUENCES ON THE LIGHT ELEMENTS ABUNDANCES

The early evolution of planetary systems is expected to depend on various periods of disk matter accretion onto the central star, which may include the accretion of metal-rich matter after the star settles on the main sequence. When this happens, the accreted material is rapidly mixed within the surface convective zone and induces an inverse mean-molecular-weight gradient, unstable for thermohaline convection. The induced mixing, which dilutes the metal excess, may also have important consequences on the light elements abundances. We model and analyze this process, and present the results according to various possible accretion scenarios. We give a detailed discussion of the different ways of treating thermohaline mixing, as proposed by previous authors, and we converge on a consistent view, including the most recent numerical simulations. We show how the observations of light elements in stars can be used as tracers of such events.

Thermohaline instability and rotation-induced mixing

Context. The 3He content of Galactic HII regions is very close to that of the Sun and the solar system, and only slightly higher than the primordial 3He abundance as predicted by the standard Big Bang nucleosynthesis. However, the classical theory of stellar evolution predicts a high production of 3He by low-mass stars, implying a strong increase of 3He with time in the Galaxy. This is the well-known "3He problem". Aims. We study the effects of thermohaline and rotation-induced mixings on the production and destruction of 3He over the lifetime of low- and intermediate-mass stars at various metallicities. Methods. We compute stellar evolutionary models in the mass range 1 to 6M\odot for four metallicities, taking into account thermohaline instability and rotation-induced mixing. For the thermohaline diffusivity we use the prescription based on the linear stability analysis, which reproduces Red Giant Branch (RGB) abundance patterns at all metallicities. Rotation-induced mixing is treated taking into account meridional circulation and shear turbulence. We discuss the effects of these processes on internal and surface abundances of 3He and on the net yields. Results. Over the whole mass and metallicity range investigated, rotation-induced mixing lowers the 3He production, as well as the upper mass limit at which stars destroy 3He. For low-mass stars, thermohaline mixing occuring beyond the RGB bump is the dominant process in strongly reducing the net 3He yield compared to standard computations. Yet these stars remain net 3He producers. Conclusions. Overall, the net 3He yields are strongly reduced compared to the standard framework predictions.

CLOSE TO THE DREDGE: PRECISE X-RAY C AND N ABUNDANCES IN λ ANDROMEDA AND ITS PRECOCIOUS RED GIANT BRANCH 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.