No Arabic abstract
Novae are the observational manifestations of thermonuclear runaways on the surface of accreting white dwarfs (WDs). Although novae are an ubiquitous phenomenon, their properties at low metallicity are not well understood. Using the publicly-available stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA), we model the evolution of accreting carbon-oxygen WDs and consider models which accrete matter with metallicity Z=0.02 or $10^{-4}$. We consider both models without mixing and with matter enriched by CO-elements assuming that mixing occurs in the process of accretion (with mixing fraction 0.25). We present and contrast ignition mass, ejected mass, recurrence period and maximum luminosity of novae for different WD masses and accretion rates for these metallicities and mixing cases. We find that models with Z = 0.02 have ignition masses and recurrence periods smaller than models with low Z, while the ejected mass and maximum luminosity are larger. Retention efficiency during novae outbursts decreases with increasing metallicity. In our implementation, inclusion of mixing at the H/He interface reduces accreted mass, ejected mass and recurrence period as compared to the no-mixing case, while the maximum luminosity becomes larger. Retention efficiency is significantly reduced, becoming negative in most of our models. For ease of use, we provide a tabular summary of our results.
The production of the neutron-capture isotopes beyond iron that we observe today in the solar system is the result of the combined contribution of the r-process, the s- process and possibly the i-process. Low-mass AGB (2 < M/Msun < 3) and massive (M >10 Msun ) stars have been identified as the sites of the s-process. In this work we consider the evolution and nucleosynthesis of low-mass AGB stars. We provide an update of the NuGrid Set models, adopting the same general physics assumptions but using an updated convective-boundary mixing model accounting for the contribution from internal gravity waves. The combined data set includes the initial masses Mzams/Msun = 2, 3 for Z = 0.03, 0.02, 0.01. These models are computed with the MESA stellar code and the evolution is followed up to the end of the AGB phase. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid mppnp code. The convective boundary mixing model leads to the formation of a 13C-pocket three times wider compared to the one obtained in the previous set of models, bringing the simulation results now in closer agreement with observations. We also discuss the potential impact of other processes inducing mixing, like rotation, adopting parametric models compatible with theory and observations. Complete yield data tables, derived data products and online analytic data access are provided.
The remarkable HST datasets from the CANDELS, HUDF09, HUDF12, ERS, and BoRG/HIPPIES programs have allowed us to map out the evolution of the UV LF from z~10 to z~4. We have identified 5859, 3001, 857, 481, 217, and 6 galaxy candidates at z~4, z~5, z~6, z~7, z~8, and z~10, respectively from the ~1000 arcmin**2 area probed. The selection of z~4-8 galaxies over the five CANDELS fields allows us to assess the cosmic variance; the largest variations are apparent at z>=7. Our new LF determinations at z~4 and z~5 span a 6-mag baseline (-22.5 to -16 AB mag). These determinations agree well with previous estimates, but the larger samples and volumes probed here result in a more reliable sampling of >L* galaxies and allow us to reassess the form of the UV LFs. Our new LF results strengthen our earlier findings to 3.4 sigma significance for a steeper faint-end slope to the UV LF at z>4, with alpha evolving from alpha=-1.64+/-0.04 at z~4 to alpha=-2.06+/-0.13 at z~7 (and alpha = -2.02+/-0.23 at z~8), consistent with that expected from the evolution of the halo mass function. With our improved constraints at the bright end, we find less evolution in the characteristic luminosity M* over the redshift range z~4 to z~7; the observed evolution in the LF is now largely represented by changes in phi*. No evidence for a non-Schechter-like form to the z~4-8 LFs is found. A simple conditional LF model based on halo growth and evolution in the M/L ratio of halos ((1+z)**-1.5) provides a good representation of the observed evolution.
The effects of rotation on stellar evolution are particularly important at low metallicity, when mass loss by stellar winds diminishes and the surface enrichment due to rotational mixing becomes relatively more pronounced than at high metallicities. Here we investigate the impact of rotation and metallicity on stellar evolution. Using a similar physics as in our previous large grids of models at Z=0.002 and Z=0.014, we compute stellar evolution models with the Geneva code for rotating and nonrotating stars with initial masses (Mini) between 1.7 and 120 Msun and Z=0.0004 (1/35 solar). This is comparable to the metallicities of the most metal poor galaxies observed so far, such as I Zw 18. Concerning massive stars, both rotating and nonrotating models spend most of their core-helium burning phase with an effective temperature higher than 8000 K. Stars become red supergiants only at the end of their lifetimes, and few RSGs are expected. Our models predict very few to no classical Wolf-Rayet stars as a results of weak stellar winds at low metallicity. The most massive stars end their lifetimes as luminous blue supergiants or luminous blue variables, a feature that is not predicted by models with higher metallicities. Interestingly, due to the behavior of the intermediate convective zone, the mass domain of stars producing pair-instability supernovae is smaller at Z=0.0004 than at Z=0.002. We find that during the main sequence phase, the ratio between nitrogen and carbon abundances (N/C) remains unchanged for nonrotating models. However, N/C increases by factors of 10-20 in rotating models at the end of the MS. Cepheids coming from stars with Mini > 4-6 Msun are beyond the core helium burning phase and spend little time in the instability strip. Since they would evolve towards cooler effective temperatures, these Cepheids should show an increase of the pulsation period as a function of age.
We compute the escape of ionizing radiation from galaxies in the redshift interval z=4-10, i.e., during and after the epoch of reionization, using a high-resolution set of galaxies, formed in fully cosmological simulations. The simulations invoke early, energetic feedback, and the galaxies evolve into a realistic population at z=0. Our galaxies cover nearly four orders of magnitude in masses (10^{7.8}-10^{11.5}msun) and more than five orders in star formation rates (10^{-3.5}-10^{1.7}msunyr^{-1}), and we include an approximate treatment of dust absorption. We show that the source-averaged Lyman-limit escape fraction at z=10.4 is close to 80% declining monotonically with time as more massive objects build up at lower redshifts. Although the amount of dust absorption is uncertain to 1-1.5 dex, it is tightly correlated with metallicity; we find that dust is unlikely to significantly impact the observed UV output. These results support reionization by stellar radiation from low-luminosity dwarf galaxies and are also compatible with Lyman continuum observations and theoretical predictions at zsim3-4.
High-resolution spectroscopy has revealed large concentrations of CNO and sometimes other intermediate-mass elements in the shells ejected during nova outbursts, suggesting that the solar composition material transferred from the secondary mixes with the outermost layers of the underlying white dwarf during the thermonuclear runaway. Multidimensional simulations have shown that Kelvin-Helmholtz instabilities provide self-enrichment of the accreted envelope with material from the outermost layers of the white dwarf, at levels that agree with observations. However, the Eulerian and time-explicit nature of most multidimensional codes used to date and the overwhelming computational load have limited their applicability, and no multidimensional simulation has been conducted for a full nova cycle. This paper explores a new methodology that combines 1-D and 3-D simulations. The early stages of the explosion (i.e., mass-accretion and initiation of the runaway) have been computed with the 1-D hydrodynamic code SHIVA. When convection extends throughout the entire envelope, the structures for each model were mapped into 3-D Cartesian grids and were subsequently followed with the multidimensional code FLASH. Two key physical quantities were extracted from the 3-D simulations and subsequently implemented into SHIVA, which was used to complete the simulation through the late expansion and ejection stages: the time-dependent amount of mass dredged-up from the outer white dwarf layers, and the time-dependent convective velocity profile throughout the envelope. More massive envelopes than those reported from previous models with pre-enrichment have been found. This results in more violent outbursts, characterized by higher peak temperatures and greater ejected masses, with metallicity enhancements in agreement with observations.