The R Coronae Borealis (RCB) stars are extremely hydrogen-deficient carbon stars which produce large amounts of dust, causing sudden deep declines in brightness. They are believed to be formed primarily through white dwarf mergers. In this paper, we use MESA to investigate how post-merger objects with a range of initial He-burning shell temperatures from 2.1 - 5.4 × 108 K with solar and subsolar metallicities evolve into RCB stars. The most successful model of these has subsolar metallicity and an initial temperature near 3 × 108 K. We find a strong dependence on initial He-burning shell temperature for surface abundances of elements involved in the CNO cycle, as well as differences in effective temperature and radius of RCBs. Elements involved in nucleosynthesis present around 1 dex diminished surface abundances in the 10% solar metallicity models, with the exception of carbon and lithium which are discussed in detail. Models with subsolar metallicities also exhibit longer lifetimes than their solar counterparts. Additionally, we find that convective mixing of the burned material occurs only in the first few years of post-merger evolution, after which the surface abundances are constant during and after the RCB phase, providing evidence for why these stars show a strong enhancement of partial He-burning products.

R Coronae Borealis (RCB) stars are rare hydrogen-deficient carbon-rich variable supergiants thought to be the result of dynamically unstable white dwarf mergers. We attempt to model RCB stars through all the relevant timescales by simulating a merger event in Octo-tiger, a 3D adaptive mesh refinement (AMR) hydrodynamics code, and mapping the post-merger object into MESA, a 1D stellar evolution code. We then post-process the nucleosynthesis on a much larger nuclear reaction network to study the enhancement of s-process elements. We present models that match observations or previous studies in most surface abundances, isotopic ratios, early evolution, and lifetimes. We also observe similar mixing behavior to previous modeling attempts that result in the partial He-burning products visible on the surface in observations. However, we do note that our subsolar models lack any enhancement in s-process elements, which we attribute to a lack of hydrogen in the envelope. We also find that the 16O/18O isotopic ratio is very sensitive to initial hydrogen abundance and increases outside of the acceptable range with a hydrogen mass fraction greater than 10−4.

We present an improved numerical method to model subsolar He+CO-WD merger progenitors of R Corona Borealis stars that builds on our previous work. These improvements include a smooth entropy transition from the core to the envelope of the post-merger, inclusion of single-zone nucleosynthesis to mimic the effects of burning during the merger event, and post-processing the models with a larger nuclear network for analysis of s-process nucleosynthesis. We perform a parameter study to understand the effects of the entropy transition, peak temperature, and overshooting on our models. The models that best agree with observations of R Corona Borealis stars are processed with a much larger nuclear network to investigate s-process nucleosynthesis and the dredge-up of s-process products into the outer envelope in detail. We present a model with a significant enhancement in s-process elements, which also agrees with observed surface abundances and isotopic ratios of 16O/16O and C/O between 1 and 10. Finally, we find that the neutron exposure and initial neutron densities this model requires to obtain such an enhancement are much more consistent with i-process nucleosynthesis.

Context. R Coronae Borealis (RCB) variables and their non-variable counterparts, the dustless Hydrogen-deficient Carbon (dLHdC) stars have been known to exhibit enhanced s-processed material on their surfaces, especially Sr, Y, and Ba. No comprehensive work has been done to explore the s-process in these types of stars, however one particular RCB star, U Aqr, has been under scrutiny for its extraordinary Sr enhancement. Aims. We aim to identify RCB and dLHdC stars that have significantly enhanced Sr abundances, such as U Aqr, and use stellar evolution models to begin to estimate the type of neutron exposure that occurs in a typical HdC star. Methods. We compared the strength of the Sr II 4077 Å spectral line to Ca II H to identify the new subclass of Sr-rich HdCs. We additionally used the structural and abundance information from existing RCB MESA models to calculate the neutron exposure parameter, τ. Results. We identified six stars in the Sr-rich class. Two are RCBs, and four are dLHdCs. We additionally found that the preferred RCB MESA model has a neutron exposure τ 0.1 mb −1 , which is lower than the estimated τ between 0.15 and 0.6 mb −1 for the Sr-rich star U Aqr found in the literature. We found trends in the neutron exposure corresponding to He-burning shell temperature, metallicity, and assumed s-processing site. Conclusions. We have found a sub-class of six HdCs known as the Sr-rich class, which tend to lie in the halo, outside the typical distribution of RCBs and dLHdCs. We found that dLHdC stars are more likely to be Sr-rich than RCBs, with an occurrence rate of ∼13% for dLHdCs and ∼2% for RCBs. This is one of the first potential spectroscopic differences between RCBs and dLHdCs, along with dLHdCs having stronger surface abundances of 18 O. We additionally found neutron exposure trends in our RCB models that will aide in understanding the interplay between model parameters and surface s-process elements.

No abstract