The nucleus of the planetary nebula Hen2-428 is a short orbital-period (4.2h), double-lined spectroscopic binary, whose status as a potential supernova type Ia progenitor has raised some controversy in the literature. With the aim of resolving this debate, we carried out an in-depth reanalysis of the system. Our approach combines a refined wavelength calibration, thorough line-identifications, improved radial-velocity measurements, non-LTE spectral modeling, as well as multi-band light-curve fitting. Our results are then discussed in view of state-of-the-art stellar evolutionary models. Besides systematic zero-point shifts in the wavelength calibration of the OSIRIS spectra which were also used in the previous analysis of the system, we found that the spectra are contaminated with diffuse interstellar bands. Our Voigt-profile radial velocity fitting method, which considers the additional absorption of these diffuse interstellar bands, reveals significantly lower masses ($M_1=0.66\pm0.11M_\odot$ and $M_2=0.42\pm0.07M_\odot$) than previously reported and a mass ratio that is clearly below unity. Our spectral and light curve analyses lead to consistent results, however, we find higher effective temperatures and smaller radii than previously reported. Moreover, we find that the red-excess that was reported before to prove to be a mere artifact of an outdated reddening law that was applied. Our work shows that blends of HeII 5412A with diffuse interstellar bands have led to an overestimation of the previously reported dynamical masses of Hen2-428. The merging event of Hen2-428 will not be recognised as a supernova type Ia, but most likely leads to the formation of a H-deficient star. We suggest that the system was formed via a first stable mass transfer episode, followed by common envelope evolution, and it is now composed of a post-early AGB star and a reheated He-core white dwarf.

The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approach using detached eclipsing binaries. We then move to more indirect and model-dependent methods, such as the quite commonly used isochrone or stellar track fitting. The arrival of quantitative asteroseismology has opened a completely new approach to determine stellar masses and to complement and improve the accuracy of other methods. We include methods for different evolutionary stages, from the pre-main sequence to evolved (super)giants and final remnants. For all methods uncertainties and restrictions will be discussed. We provide lists of altogether more than 200 benchmark stars with relative mass accuracies between $[0.3,2]\%$ for the covered mass range of $M\in [0.1,16]\,M_\odot$, $75\%$ of which are stars burning hydrogen in their core and the other $25\%$ covering all other evolved stages. We close with a recommendation how to combine various methods to arrive at a "mass-ladder" for stars.

Recent observations of protoplanetary disks have revealed ring-like structures that can be associated to pressure maxima. Pressure maxima are known to be dust collectors and planet migration traps. Most of planet formation works are based either on the pebble accretion model or on the planetesimal accretion model. However, recent studies proposed the possible formation of Jupiter by the hybrid accretion of pebbles and planetesimals. We aim to study the full process of planet formation consisting of dust evolution, planetesimal formation and planet growth at a pressure maximum in a protoplanetary disk. We compute, through numerical simulations, the gas and dust evolution, including dust growth, fragmentation, radial drift and particle accumulation at a pressure bump. We also consider the formation of planetesimals by streaming instability and the formation of a moon-size embryo that grows into a giant planet by the hybrid accretion of pebbles and planetesimals. We find that pressure maxima in protoplanetary disks are efficient collectors of dust drifting inwards. The condition of planetesimal formation by streaming instability is fulfilled due to the large amount of dust accumulated at the pressure bump. Then, a massive core is quickly formed (in $\sim 10^4$ yr) by the accretion of pebbles. After the pebble isolation mass is reached, the growth of the core slowly continues by the accretion of planetesimals. The energy released by planetesimal accretion delays the onset of runaway gas accretion, allowing a gas giant to form after $\sim$1 Myr of disk evolution. The pressure maximum also acts as a migration trap. Pressure maxima in protoplanetary disks are preferential locations for dust traps, planetesimal formation by streaming instability and planet migration traps. All these conditions allow the fast formation of a giant planet by the hybrid accretion of pebbles and planetesimals.

Planetary nebulae represent a potential late stage of stellar evolution, however the central stars (CSPNe) are relatively faint and therefore pertinent information is merely available for <20% of the Galactic sample. Consequently, the literature was surveyed to construct a new catalogue of 620 CSPNe featuring important spectral classifications and information. The catalogue supersedes the existing iteration by 25%, and includes physical parameters such as luminosity, surface gravity, temperature, magnitude estimates, and references for published spectra. The marked statistical improvement enabled the following pertinent conclusions to be determined: the H-rich/H-poor ratio is 2:1, there is a deficiency of CSPNe with types [WC 5-6], and nearly 80% of binary central stars belong to the H-rich group. The last finding suggests that evolutionary scenarios leading to the formation of binary central stars interfere with the conditions required for the formation of H-poor CSPN. Approximately 50% of the sample with derived values of log L, log Teff, and log g, exhibit masses and ages consistent with single stellar evolutionary models. The implication is that single stars are indeed able to form planetary nebulae. Moreover, it is shown that H-poor CSPNe are formed by higher mass progenitors. The catalogue is available through the Vizier database.

All evolved stars with masses $M_\star\lesssim 2M_\odot$ undergo a helium(He)-core flash at the end of their first stage as a giant star. Although theoretically predicted more than 50 years ago, this core-flash phase has yet to be observationally probed. We show here that gravity modes (g modes) stochastically excited by He-flash driven convection are able to reach the stellar surface, and induce periodic photometric variabilities in hot-subdwarf stars with amplitudes of the order of a few mmag. As such they can now be detected by space-based photometry with the Transiting Exoplanet Survey Satellite (TESS) in relatively bright stars (e.g. magnitudes $I_C\lesssim 13$). The range of predicted periods spans from a few thousand seconds to tens of thousand seconds, depending on the details of the excitation region. In addition, we find that stochastically excited pulsations reproduce the pulsations observed in a couple of He-rich hot subdwarf stars. These stars, and in particular the future TESS target Feige 46, are the most promising candidates to probe the He-core flash for the first time.

Context. The large quantity of high-quality asteroseismic data that obtained from space-based photometric missions and the accuracy of the resulting frequencies motivate a careful consideration of the accuracy of computed oscillation frequencies of stellar models, when applied as diagnostics of the model properties. Aims. Based on models of red-giant stars that have been independently calculated using different stellar evolution codes, we investigate the extent to which the differences in the model calculation affect the model oscillation frequencies. Methods. For each of the models, which cover four different masses and different evolution stages on the red-giant branch, we computed full sets of low-degree oscillation frequencies using a single pulsation code and, from these frequencies, typical asteroseismic diagnostics. In addition, we carried out preliminary analyses to relate differences in the oscillation properties to the corresponding model differences. Results. In general, the differences in asteroseismic properties between the different models greatly exceed the observational precision of these properties, in particular for the nonradial modes whose mixed acoustic and gravity-wave character makes them sensitive to the structure of the deep stellar interior. In some cases, identifying these differences led to improvements in the final models presented here and in Paper I; here we illustrate particular examples of this. Conclusions. Further improvements in stellar modelling are required in order fully to utilise the observational accuracy to probe intrinsic limitations in the modelling. However, our analysis of the frequency differences and their relation to stellar internal properties provides a striking illustration of the potential of the mixed modes of red-giant stars for the diagnostics of stellar interiors.

The treatment of convective boundaries remains an important source of uncertainty within stellar evolution, with drastic implications for the thermally-pulsing stars on the Asymptotic Giant Branch (AGB). Various sources are taken as motivation for the incorporation of convective boundary mixing during this phase, from s-process nucleosynthesis to hydrodynamical models. In spite of the considerable evidence in favour of the existence of convective boundary mixing on the pre-AGB evolution, this mixing is not universally included in models of TP-AGB stars. The aim of this investigation is to ascertain the extent of convective boundary mixing, which is compatible with observations when considering full evolutionary models. Additionally, we investigate a theoretical argument that has been made that momentum-driven overshooting at the base of the pulse-driven convection zone should be negligible. We show that, while the argument holds,estimations based on the picture of turbulent entrainment suggest that mixing is possible at both convective boundaries. We demonstrate that additional mixing at convective boundaries during core-burning phases prior to the thermally-pulsing Asymptotic Giant Branch has an impact on the later evolution, changing the mass range at which the third dredge-up and hot-bottom burning occur, and thus also the final surface composition. In addition, an effort has been made to constrain the efficiency of convective boundary mixing at the different convective boundaries, using observational constraints. Our study suggests a strong tension between different constraints that makes it impossible to reproduce all observables simultaneously within the framework of an exponentially decaying overshooting. This result calls for a reassessment of both the models of convective boundary mixing and the observational constraints.

(Abridged). We introduce the Aarhus Red Giants Challenge, a series of detailed comparisons between widely used stellar evolution and oscillation codes aiming at establishing the minimum level of uncertainties in properties of red giants arising solely from numerical implementations. Using 9 state-of-the-art stellar evolution codes, we defined a set of input physics and physical constants for our calculations and calibrated the convective efficiency to a specific point on the main sequence. We produced evolutionary tracks and stellar structure models at fixed radius along the red-giant branch for masses of 1.0 M$_\odot$, 1.5 M$_\odot$, 2.0 M$_\odot$, and 2.5 M$_\odot$, and compared the predicted stellar properties. Once models have been calibrated on the main sequence we find a residual spread in the predicted effective temperatures across all codes of ~20 K at solar radius and ~30-40 K in the RGB regardless of the considered stellar mass. The predicted ages show variations of 2-5% (increasing with stellar mass) which we track down to differences in the numerical implementation of energy generation. The luminosity of the RGB-bump shows a spread of about 10% for the considered codes, which translates into magnitude differences of ~0.1 mag in the optical V-band. We also compare the predicted [C/N] abundance ratio and found a spread of 0.1 dex or more for all considered masses. Our comparisons show that differences at the level of a few percent still remain in evolutionary calculations of red giants branch stars despite the use of the same input physics. These are mostly due to differences in the energy generation routines and interpolation across opacities, and call for further investigations on these matters in the context of using properties of red giants as benchmarks for astrophysical studies.

We describe a revised procedure for the numerical simulation of planetary nebulae luminosity functions (PNLF), improving on previous work (M\'endez & Soffner 1997). The procedure now is based on new H-burning post-AGB evolutionary tracks (Miller Bertolami 2016). For a given stellar mass, the new central stars are more luminous and evolve faster. We have slightly changed the distribution of the [OIII] 5007 intensities relative to those of H$\beta$ and the generation of absorbing factors, while still basing their numerical modeling on empirical information extracted from studies of galactic planetary nebulae (PNs) and their central stars. We argue that the assumption of PNs being completely optically thick to H-ionizing photons leads to conflicts with observations and show that to account for optically thin PNs is necessary. We then use the new simulations to estimate a maximum final mass, clarifying its meaning, and discuss the effect of internal dust extinction as a possible way of explaining the persistent discrepancy between PNLF and surface brightness fluctuation (SBF) distances. By adjusting the range of minimum to maximum final mass, it is also possible to explain the observed variety of PNLF shapes at intermediate magnitudes. The new PN formation rates are calculated to be slightly lower than suggested by previous simulations based on older post-AGB evolutionary tracks.

The transition from the asymptotic giant branch (AGB) to the final white dwarf (WD) stage is arguably the least understood phase in the evolution of single low- and intermediate-mass stars ($0.8 \lesssim M_{\rm ZAMS}/M_\odot\lesssim 8...10$). Here we briefly review the progress in the last 50 years of the modeling of stars during the post-AGB phase. We show that although the main features, like the extreme mass dependency of post-AGB timescales were already present in the earliest post-AGB models, the quantitative values of the computed post-AGB timescales changed every time new physics was included in the modeling of post-AGB stars and their progenitors. Then we discuss the predictions and uncertainties of the latest available models regarding the evolutionary timescales of post-AGB stars.

We present Gemini-South observations of nine faint and extended planetary nebulae. Using direct images taken with the spectrograph GMOS, we built the $(u' - g')$ vs. $(g' - r')$ diagrams of the stars in the observed areas which allowed us, also considering their geometrical positions, to identify the probable central stars of the nebulae. Our stellar spectra of seven stars, also taken with GMOS, indicate that four (and probably two more) objects are white dwarfs of the DAO subtype. Moreover, the white dwarf status of the four stars is confirmed by the parameters $ T_{\mathrm{eff}}$ and $ \log g$ derived with the help of theoretical stellar spectra. Given this evidence, we propose that these hot stars are the central ionizing sources of the nebulae. With this work we hope to help improve the current scarce statistics on central white dwarfs in planetary nebulae.

Asteroseismology of white dwarf (WD) stars is a powerful tool that allows to reveal the hidden chemical structure of WD and infer details about their evolution by comparing the observed periods with those obtained from stellar models. A recent asteroseismological study has reproduced the period spectrum of the helium rich pulsating WD KIC 08626021 with an unprecedented precision. The chemical structure derived from that analysis is notably different from that expected for a WD according to currently accepted formation channels, thus posing a challenge to the theory of stellar evolution. We explore the relevant micro- and macro-physics processes acting during the formation and evolution of KIC 08626021 that could lead to a chemical structure similar to that found through asteroseismology. We quantify to which extent is necessary to modify the physical processes that shapes the chemical structure, in order to reproduce the most important features of the asteroseismic model. We model the previous evolution of KIC 08626021 by exploring specific changes in the 12C+alpha reaction rate, screening processes, microscopic diffusion, as well as convective boundary mixing during core-He burning. We find that, in order to reproduce the core chemical profile derived for KIC 0862602, the 12C+alpha nuclear reaction rate has to be increased by a factor of $\sim$ 10 during the helium-core burning, and reduced by a factor of $\sim$ 1000 during the following helium-shell burning, as compared with the standard predictions for this rate. In addition, the main chemical structures derived for KIC 0862602 cannot be reconciled with our present knowledge of white dwarf formation. We find that within our current understanding of white dwarf formation and evolution, it is difficult to reproduce the most important asteroseismologically-derived features of the chemical structure of KIC 08626021.

Stellar post asymptotic giant branch (post-AGB) evolution can be completely altered by a final thermal pulse (FTP) which may occur when the star is still leaving the AGB (AFTP), at the departure from the AGB at still constant luminosity (late TP, LTP) or after the entry to the white-dwarf cooling sequence (very late TP, VLTP). Then convection mixes the He-rich material with the H-rich envelope. According to stellar evolution models the result is a star with a surface composition of $\mathrm{H}\approx\,20\,$% by mass (AFTP), $\approx 1\,$% (LTP), or (almost) no H (VLTP). Since FTP stars exhibit intershell material at their surface, spectral analyses establish constraints for AGB nucleosynthesis and stellar evolution. We performed a spectral analysis of the so-called hybrid PG 1159-type central stars (CS) of the planetary nebulae Abell 43 and NGC7094 by means of non-local thermodynamical equilibrium models. We confirm the previously determined effective temperatures of $T_\mathrm{eff} = 115\,000\pm 5\,000\,$K and determine surface gravities of $\log (g\,/\,\mathrm{cm/s^2}) = 5.6\pm 0.1$ for both. From a comparison with AFTP evolutionary tracks, we derive stellar masses of $0.57^{+0.07}_{-0.04}\,M_\odot$ and determine the abundances of H, He, and metals up to Xe. Both CS are likely AFTP stars with a surface H mass fraction of $0.25 \pm 0.03$ and $0.15 \pm 0.03$, respectively, and a Fe deficiency indicating subsolar initial metallicities. The light metals show typical PG 1159-type abundances and the elemental composition is in good agreement with predictions from AFTP evolutionary models. However, the expansion ages do not agree with evolution timescales expected from the AFTP scenario and alternatives should be explored.

White dwarf stars constitute the final evolutionary stage for more than 95 per cent of all stars. The Galactic population of white dwarfs conveys a wealth of information about several fundamental issues and are of vital importance to study the structure, evolution and chemical enrichment of our Galaxy and its components ---including the star formation history of the Milky Way. In addition, white dwarfs are tracers of the evolution of planetary systems along several phases of stellar evolution. Also, white dwarfs are used as laboratories for astro-particle physics, being their interest focused on physics beyond the standard model. The last decade has witnessed a great progress in the study of white dwarfs. In particular, a wealth of information of these stars from different surveys has allowed us to make meaningful comparison of evolutionary models with observations. While some information like surface chemical composition, temperature and gravity of isolated white dwarfs can be inferred from spectroscopy, and the total mass and radius can be derived as well when they are in binaries, the internal structure of these compact stars can be unveiled only by means of asteroseismology, an approach based on the comparison between the observed pulsation periods of variable stars and the periods predicted by appropriate theoretical models. The asteroseismological techniques allow us to infer details of the internal chemical stratification, the total mass, and even the stellar rotation profile. In this review, we first revise the evolutionary channels currently accepted that lead to the formation of white-dwarf stars, and then, we give a detailed account of the different sub-types of pulsating white dwarfs known so far, emphasizing the recent observational and theoretical advancements in the study of these fascinating variable stars.

We have performed a census of the UV-bright population in 78 globular clusters using wide-field UV telescopes. This population includes a variety of phases of post-horizontal branch (HB) evolution, including hot post-asymptotic giant branch (AGB) stars, and post-early AGB stars. There are indications that old stellar systems like globular clusters produce fewer post-(early) AGB stars than currently predicted by evolutionary models, but observations are still scarce. We obtained FORS2 spectroscopy of eleven of these UV-selected objects (covering a range of -2.3<[Fe/H]<-1.0), which we (re-)analysed together with previously observed data. We used model atmospheres of different metallicities, including super-solar ones. Where possible, we verified our atmospheric parameters using UV spectrophotometry and searched for metal lines in the optical spectra. We calculated evolutionary sequences for four metallicity regimes and used them together with information about the HB morphology of the globular clusters to estimate the expected numbers of post-AGB stars. Seven of the eleven new luminous UV-bright stars are post-AGB or post-early AGB stars, two are evolving away from the HB, one is a foreground white dwarf, and one is a white dwarf merger. So spectroscopy is clearly required to identify the evolutionary status of hot UV-bright stars. For hotter stars, metal-rich model spectra are required to reproduce their optical and UV spectra, which may affect the flux contribution of hot post-AGB stars to the UV spectra of evolved populations. Adding published information on other hot UV-bright stars in globular clusters, we find that the number of observed hot post-AGB stars generally agrees with the predicted values, although the numbers are still low.

As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that were unable to reproduce the observed properties of extrasolar planets. In the last years, several improvements have been made in order to extend the study of type I migration rates to non-isothermal discs. Moreover, it was recently shown that if the planet's luminosity due to solid accretion is taken into account, inward migration could be slowed down and even reversed. In this work, we study the planet formation process incorporating, and comparing, updated type I migration rates for non-isothermal discs and the role of planet's luminosity over such rates. We find that the latter can have important effects on planetary evolution, producing a significant outward migration for the growing planets.

We have analyzed FUSE, COS, GHRS, and Keck HIRES spectra of the UV-bright star Barnard 29 in M13 (NGC 6205). By comparing the photospheric abundances derived from multiple ionization states of C, N, O, Si, and S, we infer an effective temperature T_eff = 21,400 +/- 400 K. Balmer-line fits yield a surface gravity log g = 3.10 +/- 0.03. We derive photospheric abundances of He, C, N, O, Mg, Al, Si, P, S, Cl, Ar, Ti, Cr, Fe, Ni, and Ge. Barnard 29 exhibits an abundance pattern typical of the first-generation stars in M13, enhanced in oxygen and depleted in aluminum. An underabundance of C and an overabundance of N suggest that the star experienced nonconvective mixing on the RGB. We see no evidence of significant chemical evolution since the star left the RGB; in particular, it did not undergo third dredge-up. Previous workers found that the star's FUV spectra yield an iron abundance about 0.5 dex lower than its optical spectrum, but the iron abundances derived from all of our spectra are consistent with the cluster value. We attribute this difference to our use of model atmospheres without microturbulence, which is ruled out by careful fits to optical absorption features. We derive a mass M_*/M_sun = 0.45 - 0.55 and luminosity log (L_*/L_sun) = 3.26 - 3.35. Comparison with stellar-evolution models suggests that Barnard 29 evolved from a ZAHB star of mass M_*/M_sun between 0.50 and 0.55, near the boundary between the extreme and blue horizontal branches.

Ultra-massive white dwarfs are powerful tools to study various physical processes in the Asymptotic Giant Branch (AGB), type Ia supernova explosions and the theory of crystallization through white dwarf asteroseismology. Despite the interest in these white dwarfs, there are few evolutionary studies in the literature devoted to them. Here, we present new ultra-massive white dwarf evolutionary sequences that constitute an improvement over previous ones. In these new sequences, we take into account for the first time the process of phase separation expected during the crystallization stage of these white dwarfs, by relying on the most up-to-date phase diagram of dense oxygen/neon mixtures. Realistic chemical profiles resulting from the full computation of progenitor evolution during the semidegenerate carbon burning along the super-AGB phase are also considered in our sequences. Outer boundary conditions for our evolving models are provided by detailed non-gray white dwarf model atmospheres for hydrogen and helium composition. We assessed the impact of all these improvements on the evolutionary properties of ultra-massive white dwarfs, providing up-dated evolutionary sequences for these stars. We conclude that crystallization is expected to affect the majority of the massive white dwarfs observed with effective temperatures below $40\,000\, \rm K$. Moreover, the calculation of the phase separation process induced by crystallization is necessary to accurately determine the cooling age and the mass-radius relation of massive white dwarfs. We also provide colors in the GAIA photometric bands for our H-rich white dwarf evolutionary sequences on the basis of new models atmospheres. Finally, these new white dwarf sequences provide a new theoretical frame to perform asteroseismological studies on the recently detected ultra-massive pulsating white dwarfs.

Identifying progenitor systems for the double-degenerate scenario is crucial to check the reliability of type Ia super-novae as cosmological standard candles. Santander-Garcia et al. (2015) claimed that Hen 2-428 has a double-degenerate core whose combined mass significantly exceeds the Chandrasekhar limit. Together with the short orbital period (4.2 hours), the authors concluded that the system should merge within a Hubble time triggering a type Ia supernova event. Garcia-Berro et al. (2016) explored alternative scenarios to explain the observational evidence, as the high mass conclusion is highly unlikely within predictions from stellar evolution theory. They conclude that the evidence supporting the supernova progenitor status of the system is premature. Here we present the first quantitative spectral analysis of Hen 2-428 which allows us to derive the effective temperatures, surface gravities and helium abundance of the two CSPNe based on state-of-the-art, non-LTE model atmospheres. These results provide constrains for further studies of this particularly interesting system.

Hot subdwarfs represent a poorly understood late phase of stellar evolution. While binary evolution plays an important role for B-subdwarfs (sdB), the origin of the He dominated O-subdwarfs (He-sdO) is unknown. We search for chemical signatures of their genesis by means of quantitative analyses of visual and FUV spectra. Using Tlusty/Synspec to compute line blanketed NLTE model atmospheres, Teff/log g and abundances have been derived for 4 prototypical He-sdO stars. Final models included H,He,C,N,O,Ne,Mg,Al,Si,P,S,Fe&Ni. Because of the enrichment of either N or C, models including these elements at the appropriate high abundance provide sufficiently accurate approximations to the temperature stratification of full models. No indications for binarity were found, neither radial velocity variations nor photometric evidence. All stars have atmospheres almost free of H and Teff =42000K to 47000K, and log g= 5.4 to 5.7. CD-31 4800 displays the signatures of CNO burning, while heavier elements are subsolar by ~0.4 dex, except Ne & Si being ~solar. The C-rich He-sdOs show a slightly subsolar metallicity accompanied by N-enrichment and O-deficiency. Ne is enriched up to a factor of ten with respect to the sun in LS IV+10 9. Stellar masses scatter around the canonical helium core flash mass. The abundance pattern of CD-31 4800 is consistent with predictions of models for cold mergers of He-WDs except for the low O abundance observed. Models for composite mergers predict abundance pattern dissimilar to those determined for the C-rich stars. [CW83]0904-02, though, may be a candidate for a composite He-WD merger, as it rotates and is more massive than the other stars. New evolutionary models for the hot flasher scenario predict abundance patterns similar to those determined for the C-rich stars. Hence, C-rich He-sdO may result from late He flashes with deep-mixing episodes. (abbreviated)

The nucleus of Hen 2-428 is a short orbital period (4.2 h) spectroscopic binary, whose status as potential supernovae type Ia progenitor has raised some controversy in the literature. We present preliminary results of a thorough analysis of this interesting system, which combines quantitative non-local thermodynamic equilibrium spectral modelling (non-LTE), radial velocity analysis,multi-band light curve fitting, and state-of-the art stellar evolutionary calculations. Importantly, we find that the dynamical system mass that is derived by using all available He II lines does not exceed the Chandrasekhar mass limit. Furthermore, the individual masses of the two central stars are too small to lead to an SN Ia in case of a dynamical explosion during the merger process.

Planetary nebulae are ionized clouds of gas formed by the hydrogen-rich envelopes of low- and intermediate-mass stars ejected at late evolutionary stages. The strong UV flux from their central stars causes a highly stratified ionization structure, with species of higher ionization potential closer to the star. Here we report on the exceptional case of HuBi 1, a double-shell planetary nebula whose inner shell presents emission from low-ionization species close to the star and emission from high-ionization species farther away. Spectral analysis demonstrates that the inner shell of HuBi 1 is excited by shocks, whereas its outer shell is recombining. The anomalous excitation of these shells can be traced to its low-temperature [WC10] central star whose optical brightness has declined continuously by 10 magnitudes in a period of 46 years. Evolutionary models reveal that this star is the descendent of a low-mass star ($\simeq$1.1 $M_\odot$) that has experienced a born-again event whose ejecta shock-excite the inner shell. HuBi 1 represents the missing link in the formation of metal-rich central stars of planetary nebulae from low-mass progenitors, offering unique insight regarding the future evolution of the born-again Sakurai's object. Coming from a solar-mass progenitor, HuBi 1 represents a potential end-state for our Sun.

Planetary Nebulae (PN) emit enormous amount of energy in several emission lines. Measuring the line-flux for PNe in a given stellar population, the Planetary Nebula Luminosity Function (PNLF) can be compiled. Surveys of PNe revealed that the faint-end of the PNLF can be approximated by a simple exponential dependency expected for an expanding spherical shell. However at the bright-end there exists a steep cut-off which was unexpected and remains unexplained. Interestingly, the cut-off value appears to be nearly the same for different stellar populations as young spiral galaxies and old elliptical galaxies and, despite the lack of understanding, became an extragalactic distance estimator. Here we show that the recently computed post-AGB evolutionary tracks are capable to explain the decades old mystery. All new models with ages between 1 and 7 Gyr (progenitor masses between 2.0 and 1.1 of solar mass) evolve fast enough to ionize the PN, and at similar post-AGB luminosity which allows the [O III] 500.7nm line to reach nearly the same magnitude. The new models predict that the Sun at the end of its life will form a rather faint PN.

The $\epsilon$ mechanism is a self-excitation mechanism of stellar pulsations which acts in regions inside the star where nuclear burning takes place. It has been shown that the $\epsilon$ mechanism can excite pulsations in models of hot pre-horizontal branch stars, and that the shortest periods of LS IV-14$^{\circ}$116 could be explained that way. We aim to study the $\epsilon$ mechanism in stellar models appropriate for hot pre-horizontal branch stars to predict their pulsational properties and the instability domain in the $\log g-\log T_{\rm eff}$ plane. We perform detailed computations of non-adiabatic non-radial pulsations on stellar models during the helium subflashes just before the helium-core burning phase. We find an instability domain of long-period gravity modes due to the $\epsilon$ mechanism in the $\log g-\log T_{\rm eff}$ plane at roughly $22000\,{\rm K} \lesssim T_{\rm eff}\lesssim 50000\,$K and $4.67 \lesssim \log g \lesssim 6.15$. Consequently, we find instabilities due to the $\epsilon$ mechanism on pre-extreme horizontal branch stellar models ($T_{\rm eff}\gtrsim 22000\,$K), but not on pre-blue horizontal branch stellar models ($T_{\rm eff}\lesssim 21000\,$K). The periods of excited modes range between $\sim 200$ and $\sim2000\,$s. Comparison with the three known pulsating He-rich subdwarfs shows that $\epsilon$ mechanism can excite gravity modes in stars with similar surface properties, but in our models it is only able to excite modes in the range of the shortest observed periods. We predict a new instability strip for hot-subdwarf stars of which LS IV-14$^{\circ}$116 could be the first inhabitant. Based on simple estimates we expect 1 to 10 stars in the current samples of hot-subdwarf stars to be pulsating by the $\epsilon$ mechanism. Our results could constitute a theoretical basis for future searches of pulsators in the Galactic field.

We report deep spectroscopy of ten planetary nebulae (PNe) in the Andromeda Galaxy (M31) using the 10.4m GTC. Our targets reside in different regions of M31, including halo streams and dwarf satellite M32, and kinematically deviate from the extended disk. The temperature-sensitive [O III] 4363 auroral line is observed in all targets. For four PNe, the GTC spectra extend beyond 1 micron, enabling explicit detection of the [S III] 6312 and 9069,9531 lines and thus determination of the [S III] temperature. Abundance ratios are derived and generally consistent with AGB model predictions. Our PNe probably all evolved from low-mass (<2 M_sun) stars, as analyzed with the most up-to-date post-AGB evolutionary models, and their main-sequence ages are mostly ~2-5 Gyr. Compared to the underlying, smooth, metal-poor halo of M31, our targets are uniformly metal-rich ([O/H]> -0.4), and seem to resemble the younger population in the stream. We thus speculate that our halo PNe formed in the Giant Stream's progenitor through extended star formation. Alternatively, they might have formed from the same metal-rich gas as did the outer-disk PNe, but was displaced into their present locations as a result of galactic interactions. These interpretations are, although speculative, qualitatively in line with the current picture, as inferred from previous wide-field photometric surveys, that M31's halo is the result of complex interactions and merger processes. The behavior of N/O of the combined sample of the outer-disk and our halo/substructure PNe signifies that hot bottom burning might actually occur at <3 M_sun, but careful assessment is needed.

The primary goal of this paper is to make a direct comparison between the measured and model-predicted abundances of He, C and N in a sample of 35 well-observed Galactic planetary nebulae (PN). All observations, data reductions, and abundance determinations were performed in house to ensure maximum homogeneity. Progenitor star masses (M < 4M_sun) were inferred using two published sets of post-AGB model tracks and L and T_eff values. We conclude the following: 1) the mean values of N/O across the progenitor mass range exceeds the solar value, indicating significant N enrichment in the majority of our objects; 2) the onset of hot bottom burning appears to begin around 2 solar masses, i.e., lower than ~5 M_sun implied by theory; 3) most of our objects show a clear He enrichment, as expected from dredge-up episodes; 4) the average sample C/O value is 1.23, consistent with the effects of third dredge-up; and 5) model grids used to compare to observations successfully span the distribution over metallicity space of all C/O and many He/H data points but mostly fail to do so in the case of N/O. The evident enrichment of N in PN and the general discrepancy between the observed and model-predicted N/O abundance ratios signal the need for extra-mixing as an effect of rotation and/or thermohaline mixing in the models. The unexpectedly high N enrichment that is implied here for low mass stars, if confirmed, will likely impact our conclusions about the source of N in the Universe.

We present a spectral analysis of the UV-bright star Y453 in M4. Model fits to the star's optical spectrum yield T_eff ~ 56,000 K. Fits to the star's FUV spectrum, obtained with the Cosmic Origins Spectrograph (COS) on board the Hubble Space Telescope, reveal it to be considerably hotter, with T_eff ~ 72,000 K. We adopt T_eff = 72,000 +/- 2000 K and log g = 5.7 +/- 0.2 as our best-fit parameters. Scaling the model spectrum to match the star's optical and near-infrared magnitudes, we derive a mass M_* = 0.53 +/- 0.24 M_sun and luminosity log L/L_sun = 2.84 +/- 0.05, consistent with the values expected of an evolved star in a globular cluster. Comparing the star with post-horizontal branch evolutionary tracks, we conclude that it most likely evolved from the blue horizontal branch, departing the AGB before third dredge-up. It should thus exhibit the abundance pattern (O-poor and Na-rich) characteristic of the second-generation (SG) stars in M4. We derive the star's photospheric abundances of He, C, N, O, Si, S, Ti, Cr, Fe, and Ni. CNO abundances are roughly 0.25 dex greater than those of the cluster's SG stars, while the Si and S abundances agree match the cluster values. Abundances of the iron-peak elements (except for iron itself) are enhanced by 1 to 3 dex. Rather than revealing the star's origin and evolution, this pattern reflects the combined effects of diffusive and mechanical processes in the stellar atmosphere.

The discovery of giant planets in wide orbits represents a major challenge for planet formation theory. In the standard core accretion paradigm planets are expected to form at radial distances $\lesssim 20$ au in order to form massive cores (with masses $\gtrsim 10~\textrm{M}_{\oplus}$) able to trigger the gaseous runaway growth before the dissipation of the disc. This has encouraged authors to find modifications of the standard scenario as well as alternative theories like the formation of planets by gravitational instabilities in the disc to explain the existence of giant planets in wide orbits. However, there is not yet consensus on how these systems are formed. In this letter, we present a new natural mechanism for the formation of giant planets in wide orbits within the core accretion paradigm. If photoevaporation is considered, after a few Myr of viscous evolution a gap in the gaseous disc is opened. We found that, under particular circumstances planet migration becomes synchronised with the evolution of the gap, which results in an efficient outward planet migration. This mechanism is found to allow the formation of giant planets with masses $M_p\lesssim 1 M_{\rm Jup}$ in wide stable orbits as large as $\sim$130 au from the central star.

The post-asymptotic giant branch (AGB) phase is arguably one of the least understood phases of the evolution of low- and intermediate- mass stars. The recent post-AGB evolutionary sequences computed by Miller Bertolami (2016) are at least three to ten times faster than those previously published by Vassiliadis & Wood (1994) and Bloecker (1995) which have been used in a large number of studies. This is true for the whole mass and metallicity range. The new models are also $\sim$0.1-0.3 dex brighter than the previous models with similar remnant masses. In this short article we comment on the main reasons behind these differences, and discuss possible implications for other studies of post-AGB stars or planetary nebulae.

For 32 central stars of PNe we present their parameters interpolated among the new evolutionary sequences. The derived stellar final masses are confined between 0.53 and 0.58 $M_\odot$ in good agreement with the peak in the white dwarf mass distribution. Consequently, the inferred star formation history of the Galactic bulge is well restricted between 3 and 11 Gyr and is compatible with other published studies. The new evolutionary tracks proved a very good as a tool for analysis of late stages of stars life. The result provide a compelling confirmation of the accelerated post-AGB evolution.

Systematic photometric and asteroseismological studies in the last decade support the belief that white dwarfs in the solar neighborhood harbor a broad range of hydrogen-layer contents. The reasons behind this spread of hydrogen-layer masses are not understood and usually misunderstood. In this work we present, and review, the different mechanisms that can (or cannot) lead to the formation of white dwarfs with a broad range hydrogen contents.

SAO244567 is a rare example of a star that allows us to witness stellar evolution in real time. Between 1971 and 1990 it changed from a B-type star into the hot central star of the Stingray Nebula. This observed rapid heating has been a mystery for decades, since it is in strong contradiction with the low mass of the star and canonical post-asymptotic giant branch (AGB) evolution. We speculated that SAO244567 might have suffered from a late thermal pulse (LTP) and obtained new observations with HST/COS to follow the evolution of the surface properties of SAO244567 and to verify the LTP hypothesis. Our non-LTE spectral analysis reveals that the star cooled significantly since 2002 and that its envelope is now expanding. Therefore, we conclude that SAO244567 is currently on its way back towards the AGB, which strongly supports the LTP hypothesis. A comparison with state-of-the-art LTP evolutionary calculations shows that these models cannot fully reproduce the evolution of all surface parameters simultaneously, pointing out possible shortcomings of stellar evolution models. Thereby, SAO244567 keeps on challenging stellar evolution theory and we highly encourage further investigations.

We employ an asteroseismic model of L19-2, a relatively massive ($M_{\star} \sim 0.75 M_{\sun}$) and hot ($T_{\rm eff} \sim 12\,100$ K) pulsating DA (H-rich atmosphere) white dwarf star (DAV or ZZ Ceti variable), and use the observed values of the temporal rates of period change of its dominant pulsation modes ($\Pi \sim 113$ s and $\Pi \sim 192$ s), to derive a new constraint on the mass of the axion, the hypothetical non-barionic particle considered as a possible component of the dark matter of the Universe. If the asteroseismic model employed is an accurate representation of L19-2, then our results indicate hints of extra cooling in this star, compatible with emission of axions of mass $m_{\rm a} \cos^2 \beta \lesssim 25$ meV or an axion-electron coupling constant of $g_{\rm ae} \lesssim 7 \times 10^{-13}$.

The Peccei-Quinn mechanism proposed to solve the CP problem of Quantum Chromodynamics has as consequence the existence of axions, hypothetical weakly interacting particles whose mass is constrained to be on the sub-eV range. If these particles exist and interact with electrons, they would be emitted from the dense interior of white dwarfs, becoming an important energy sink for the star. Due to their well known physics, white dwarfs are good laboratories to study the properties of fundamental particles such as the axions. We study the general effect of axion emission on the evolution of helium-rich white dwarfs and on their pulsational properties. To this aim, we calculate evolutionary helium-rich white dwarf models with axion emission, and asses the pulsational properties of this models. Our results indicate that the rates of change of pulsation periods are significantly affected by the existence of axions. We are able for the first time to independently constrain the mass of the axion from the study of pulsating helium-rich white dwarfs. To do this, we use an estimation of the rate of change of period of the pulsating white dwarf PG 1351+489 corresponding to the dominant pulsation period. From an asteroseismological model of PG 1351+489 we obtain $g_{ae}<3.3\times10^{-13}$ for the axion-electron coupling constant, or $m_a\cos^2{\beta}\lesssim$ 11.5 meV for the axion mass. This constraint is relaxed to $g_{ae}<5.5\times10^{-13}$ ($m_a\cos^2{\beta}\lesssim$ 19.5 meV), when no detailed asteroseismological model is adopted for the comparison with observations.

Because of the large neutron excess of $^{22}$Ne, this isotope rapidly sediments in the interior of the white dwarfs. This process releases an additional amount of energy, thus delaying the cooling times of the white dwarf. This influences the ages of different stellar populations derived using white dwarf cosmochronology. Furthermore, the overabundance of $^{22}$Ne in the inner regions of the star, modifies the Brunt-V\"ais\"al\"a frequency, thus altering the pulsational properties of these stars. In this work, we discuss the impact of $^{22}$Ne sedimentation in white dwarfs resulting from Solar metallicity progenitors ($Z=0.02$). We performed evolutionary calculations of white dwarfs of masses $0.528$, $0.576$, $0.657$ and $0.833$ M$_{\sun}$, derived from full evolutionary computations of their progenitor stars, starting at the Zero Age Main Sequence all the way through central hydrogen and helium burning, thermally-pulsing AGB and post-AGB phases. Our computations show that at low luminosities ($\log(L/L_{\sun})\la -4.25$), $^{22}$Ne sedimentation delays the cooling of white dwarfs with Solar metallicity progenitors by about 1~Gyr. Additionally, we studied the consequences of $^{22}$Ne sedimentation on the pulsational properties of ZZ~Ceti white dwarfs. We find that $^{22}$Ne sedimentation induces differences in the periods of these stars larger than the present observational uncertainties, particularly in more massive white dwarfs.

Recently Pasetto et al. have proposed a new method to derive a convection theory appropriate for the implementation in stellar evolution codes. Their approach is based on the simple physical picture of spherical bubbles moving within a potential flow in dynamically unstable regions, and a detailed computation of the bubble dynamics. Based on this approach the authors derive a new theory of convection which is claimed to be parameter free, non-local and time-dependent. This is a very strong claim, as such a theory is the holy grail of stellar physics. Unfortunately we have identified several distinct problems in the derivation which ultimately render their theory inapplicable to any physical regime. In addition we show that the framework of spherical bubbles in potential flows is unable to capture the essence of stellar convection, even when equations are derived correctly.

The Post Asymptotic Giant Branch (AGB) phase is arguably one of the least understood phases of the evolution of low- and intermediate- mass stars. The two grids of models presently available are based on outdated micro- and macro-physics and do not agree with each other. We study the timescales of post-AGB and CSPNe in the context of our present understanding of the micro- and macro-physics of stars. We want to assess whether new post-AGB models, based on the latter improvements in TP-AGB modeling, can help to understand the discrepancies between observation and theory and within theory itself. We compute a grid of post-AGB full evolutionary sequences that include all previous evolutionary stages from the Zero Age Main Sequence to the White Dwarf phase. Models are computed for initial masses between 0.8 and 4 $M_\odot$ and for a wide range of initial metallicities ($Z_0=$0.02, 0.01, 0.001, 0.0001), this allow us to provide post-AGB timescales and properties for H-burning post-AGB objects with masses in the relevant range for the formation of planetary nebulae ($\sim$ 0.5 - 0.8, $M_\odot$). We find post-AGB timescales that are at least $\sim 3$ to $\sim 10$ times shorter than those of old post-AGB stellar evolution models. This is true for the whole mass and metallicity range. The new models are also $\sim$ 0.1 - 0.3 dex brighter than the previous models with similar remnant masses. Post-AGB timescales show only a mild dependence on metallicity. The shorter post-AGB timescales derived in the present work are in agreement with recent semiempirical determinations of the post-AGB timescales from the CSPNe in the Galactic Bulge. Due to the very different post-AGB crossing times, initial-final mass relation and luminosities of the present models, they will have a significant impact in the predictions for the formation of planetary nebulae and the planetary nebulae luminosity function.

We present new white dwarf evolutionary sequences for low-metallicity progenitors. White dwarf sequences have been derived from full evolutionary calculations that take into account the entire history of progenitor stars, including the thermally-pulsing and the post-asymptotic giant branch phases. We show that for progenitor metallicities in the range 0.00003--0.001, and in the absence of carbon enrichment due to the occurrence of a third dredge-up episode, the resulting H envelope of the low-mass white dwarfs is thick enough to make stable H burning the most important energy source even at low luminosities. This has a significant impact on white dwarf cooling times. This result is independent of the adopted mass-loss rate during the thermally-pulsing and post-AGB phases, and the planetary nebulae stage. We conclude that in the absence of third dredge-up episodes, a significant part of the evolution of low-mass white dwarfs resulting from low-metallicity progenitors is dominated by stable H burning. Our study opens the possibility of using the observed white dwarf luminosity function of low-metallicity globular clusters to constrain the efficiency of third dredge up episodes during the thermally-pulsing AGB phase of low-metallicity progenitors.

Preliminary results from an ongoing project to compute a grid of post-AGB models is presented. Our preliminary results show that stellar evolution computations that include an updated treatment of the microphysics predict post-AGB timescales that are several times shorter that predicted by older models. Also the mass-luminosity relation of post-AGB models deviates from that of older grids. In addition, our results suggest only a slight metallicity dependence of the post-AGB timescales. We expect these results to have significant consequences for models of the formation of planetary nebulae and their luminosity function.

Recent determinations of the white dwarf luminosity function (WDLF) from very large surveys have extended our knowledge of the WDLF to very high luminosities. It has been shown that the shape of the luminosity function of white dwarfs (WDLF) is a powerful tool to test the possible properties and existence of fundamental weakly interacting subelectronvolt particles. This, together with the availability of new full evolutionary white dwarf models that are reliable at high luminosities, have opened the possibility of testing particle emission in the core of very hot white dwarfs. We use the available WDLFs from the Sloan Digital Sky Survey and the SuperCOSMOS Sky Survey to constrain the values of the neutrino magnetic dipole moment ($\mu_\nu$) and the axion-electron coupling constant ($g_{ae}$) of DFSZ-axions.

It has been shown that the shape of the luminosity function of white dwarfs (WDLF) is a powerful tool to check for the possible existence of DFSZ-axions, a proposed but not yet detected type of weakly interacting particles. With the aim of deriving new constraints on the axion mass, we compute in this paper new theoretical WDLFs on the basis of WD evolving models that incorporate for the feedback of axions on the thermal structure of the white dwarf. We find that the impact of the axion emission into the neutrino emission can not be neglected at high luminosities ($M_{\rm Bol}\lesssim 8$) and that the axion emission needs to be incorporated self-consistently into the evolution of the white dwarfs when dealing with axion masses larger than $m_a\cos^2\beta\gtrsim 5$ meV (i.e. axion-electron coupling constant $g_{ae}\gtrsim 1.4\times 10^{-13}$). We went beyond previous works by including 5 different derivations of the WDLF in our analysis. Then we have performed $\chi^2$-tests to have a quantitative measure of the assessment between the theoretical WDLFs ---computed under the assumptions of different axion masses and normalization methods--- and the observed WDLFs of the Galactic disk. While all the WDLF studied in this work disfavour axion masses in the range suggested by asteroseismology ($m_a\cos^2\beta\gtrsim 10$ meV; $g_{ae}\gtrsim 2.8\times 10^{-13}$) lower axion masses can not be discarded from our current knowledge of the WDLF of the Galactic Disk. A larger set of completely independent derivations of the WDLF of the galactic disk as well as a detailed study of the uncertainties of the theoretical WDLFs is needed before quantitative constraints on the axion-electron coupling constant can be made.

47 Tucanae is one of the most interesting and well observed and theoretically studied globular clusters. This allows us to study the reliability of our understanding of white dwarf cooling sequences, to confront different methods to determine its age, and to assess other important characteristics, like its star formation history. Here we present a population synthesis study of the cooling sequence of the globular cluster 47 Tucanae. In particular, we study the distribution of effective temperatures, the shape of the color-magnitude diagram, and the corresponding magnitude and color distributions. We do so using an up-to-date population synthesis code based on Monte Carlo techniques, that incorporates the most recent and reliable cooling sequences and an accurate modeling of the observational biases. We find a good agreement between our theoretical models and the observed data. Thus, our study, rules out previous claims that there are still missing physics in the white dwarf cooling models at moderately high effective temperatures. We also derive the age of the cluster using the termination of the cooling sequence, obtaining a good agreement with the age determinations using the main-sequence turn-off. Finally, we find that the star formation history of the cluster is compatible with that btained using main sequence stars, which predict the existence of two distinct populations. We conclude that a correct modeling of the white dwarf population of globular clusters, used in combination with the number counts of main sequence stars provides an unique tool to model the properties of globular clusters.

Pulsating white dwarf stars can be used as astrophysical laboratories to constrain the properties of weakly interacting particles. Comparing the cooling rates of these stars with the expected values from theoretical models allows us to search for additional sources of cooling due to the emission of axions, neutralinos, or neutrinos with magnetic dipole moment. In this work, we derive an upper bound to the neutrino magnetic dipole moment using an estimate of the rate of period change of the pulsating DB white dwarf star PG 1351+489. By comparing the theoretical rate of change of period expected for this star with the rate of change of period with time of PG 1351+489, we assess the possible existence of additional cooling by neutrinos with magnetic dipole moment. Our models suggest the existence of some additional cooling in this pulsating DB white dwarf, consistent with a non-zero magnetic dipole moment. Our upper limit for the neutrino magnetic dipole moment is somewhat less restrictive than, but still compatible with, other limits inferred from the white dwarf luminosity function or from the color-magnitude diagram of the Globular cluster M5. Further improvements of the measurement of the rate of period change of the dominant pulsation mode of PG 1351+489 will be necessary to confirm our bound.

Recent determinations of the white dwarf luminosity function (WDLF) from very large surveys have extended our knowledge of the WDLF to very high luminosities. This, together with the availability of new full evolutionary white dwarf models that are reliable at high luminosities, have opened the possibility of testing particle emission in the core of very hot white dwarfs, where neutrino processes are dominant. We use the available WDLFs from the Sloan Digital Sky Survey and the SuperCOSMOS Sky Survey to constrain the value of the neutrino magnetic dipole moment ($\mu_\nu$). We constructed theoretical WDLFs for different values of $\mu_\nu$ and performed a $\chi^2$-test to derive constraints on the value of $\mu_\nu$. We also constructed a unified WDLF by averaging the SDSS and SSS and estimated the uncertainties by taking into account the differences between the WDLF at each magnitude bin. Then we compared all WDLFs with theoretical WDLFs.Comparison between theoretical WDLFs and both the SDSS and the averaged WDLF indicates that $\mu_\nu$ should be $\mu_\nu<5\times 10^{-12}\, e\hbar/(2m_e c)$. In particular, a $\chi^2$-test on the averaged WDLF suggests that observations of the disk WDLF exclude values of $\mu_\nu>5\times 10^{-12}e\hbar/(2m_e c)$ at more than a 95\% confidence level, even when conservative estimates of the uncertainties are adopted. Our study shows that modern WDLFs, which extend to the high-luminosity regime, are an excellent tool for constraining the emission of particles in the core of hot white dwarfs. However, discrepancies between different WDLFs suggest there might be some relevant unaccounted systematic errors. A larger set of completely independent WDLFs, as well as more detailed studies of the theoretical WDLFs and their own uncertainties, is desirable to explore the systematic uncertainties behind this constraint.

We discuss the impact of residual nuclear burning in the cooling sequences of hydrogen-rich DA white dwarfs with very low metallicity progenitors ($Z=0.0001$). These cooling sequences are appropriate for the study of very old stellar populations. The results presented here are the product of self-consistent, fully evolutionary calculations. Specifically, we follow the evolution of white dwarf progenitors from the zero-age main sequence through all the evolutionary phases, namely the core hydrogen-burning phase, the helium-burning phase, and the thermally pulsing asymptotic giant branch phase to the white dwarf stage. This is done for the most relevant range of main sequence masses, covering the most usual interval of white dwarf masses --- from $0.53\, M_{\sun}$ to $0.83\, M_{\sun}$. Due to the low metallicity of the progenitor stars, white dwarfs are born with thicker hydrogen envelopes, leading to more intense hydrogen burning shells as compared with their solar metallicity counterparts. We study the phase in which nuclear reactions are still important and find that nuclear energy sources play a key role during long periods of time, considerably increasing the cooling times from those predicted by standard white dwarf models. In particular, we find that for this metallicity and for white dwarf masses smaller than about $0.6\, M_{\sun}$, nuclear reactions are the main contributor to the stellar luminosity for luminosities as low as $\log(L/L_{\sun})\simeq -3.2$. This, in turn, should have a noticeable impact in the white dwarf luminosity function of low-metallicity stellar populations.

We provide a fine and homogeneous grid of evolutionary sequences for He-core white dwarfs with masses 0.15-0.45 Msun, including the mass range for ELM white dwarfs (<0.20Msun). The grid is appropriate for mass and age determination, and to study their pulsational properties. White dwarf sequences have been computed by performing full evolutionary calculations that consider the main energy sources and processes of chemical abundance changes during white dwarf evolution. Initial models for the evolving white dwarfs have been obtained by computing the non-conservative evolution of a binary system consisting of a Msun ZAMS star and a 1.4 Msun neutron star for various initial orbital periods. To derive cooling ages and masses for He-core white dwarf we perform a least square fitting of the M(Teff, g) and Age(Teff, g) relations provided by our sequences by using a scheme that takes into account the time spent by models in different regions of the Teff-g plane. This is useful when multiple solutions for cooling age and mass determinations are possible in the case of CNO-flashing sequences. We also explore the adiabatic pulsational properties of models near the critical mass for the development of CNO flashes (~0.2 Msun). This is motivated by the discovery of pulsating white dwarfs with stellar masses near this threshold value. We obtain reliable and homogeneous mass and cooling age determinations for 58 very low-mass white dwarfs, including 3 pulsating stars. Also, we find substantial differences in the period spacing distributions of g-modes for models with stellar masses ~ 0.2 Msun, which could be used as a seismic tool to distinguish stars that have undergone CNO flashes in their early cooling phase from those that have not. Finally, for an easy application of our results, we provide a reduced grid of values useful to obtain masses and ages of He-core white dwarf.

It has been shown that the shape of the luminosity function of white dwarfs can be a powerful tool to check for the possible existence of DFSZ-axions. In particular, Isern et al. (2008) showed that, if the axion mass is of the order of a few meV, then the white dwarf luminosity function is sensitive enough to detect their existence. For axion masses of about $m_a > 5$ meV the axion emission can be a primary cooling mechanism for the white dwarf and the feedback of the axion emission into the thermal structure of the white dwarf needs to be considered. Here we present computations of white dwarf cooling sequences that take into account the effect of axion emission in a self consistent way by means of full stellar evolution computations. Then, we study and discuss the impact of the axion emission in the white dwarf luminosity function.

Helium-rich subdwarfs are a rare subclass of hot subdwarf stars which constitute a small and inhomogeneous group showing varying degrees of helium enrichment. Only one star, LS IV $^\circ$14 116 has been found to show multiperiodic luminosity variations. The variability of LS IV $^\circ$14 116 has been explained as the consequence of nonradial g-mode oscillations, whose excitation is difficult to understand within the frame of the standard $\kappa$-mechanism driving pulsations in sdBV stars. In a recent study, we have proposed that the pulsations of LS IV $^\circ$14 116 might be driven through the $\epsilon$-mechanism acting in unstable He-burning zones in the interior of the star, that appear before the quiescent He-burning phase. One of the few accepted scenarios for the formation of He-rich subdwarfs is the merger of two He-core white dwarfs. As part of this project, we present a study of the $\epsilon$-mechanism in post-merger remnants, and discuss the results in the light of the pulsations exhibited by LS IV $^\circ$14 116.

The development of surveys which will be able to cover a large region of the sky several times per year will allow the massive detection of transient events taking place in timescales of years. In addition, the projected full digitalization of the Harvard plate collection will open a new window to identify slow transients taking place in timescales of centuries. In particular, these projects will allow the detection of stars undergoing slow eruptions as those expected during late helium flashes in the post-AGB evolution. In order to identify those transients which correspond with late helium flashes the development of synthetic light curves of those events is mandatory. In this connection we present preliminary results of a project aimed at computing grids of theoretical light curves of born again stars.

White dwarf evolution is essentially a gravothermal cooling process, which,for cool white dwarfs, sensitively depends on the treatment of the outer boundary conditions. We provide detailed outer boundary conditions appropriate for computing the evolution of cool white dwarfs employing detailed non-gray model atmospheres for pure H composition. We also explore the impact on the white dwarf cooling times of different assumptions for energy transfer in the atmosphere of cool white dwarfs. Detailed non-gray model atmospheres are computed taken into account non-ideal effects in the gas equation of state and chemical equilibrium, collision-induced absorption from molecules, and the Lyman alpha quasi-molecular opacity. Our results show that the use of detailed outer boundary conditions becomes relevant for effective temperatures lower than 5800 and 6100K for sequences with 0.60 and 0.90 M_sun, respectively. Detailed model atmospheres predict ages that are up to approx 10% shorter at log L/L_sun=-4 when compared with the ages derived using Eddington-like approximations at tau_Ross=2/3. We also analyze the effects of various assumptions and physical processes of relevance in the calculation of outer boundary conditions. In particular, we find that the Ly_alpha red wing absorption does not affect substantially the evolution of white dwarfs. White dwarf cooling timescales are sensitive to the surface boundary conditions for T_eff < 6000K. Interestingly enough, non-gray effects have little consequences on these cooling times at observable luminosities. In fact, collision-induced absorption processes, which significantly affect the spectra and colors of old white dwarfs with hydrogen-rich atmospheres, have not noticeable effects in their cooling rates, except throughout the Rosseland mean opacity.

We employ a state-of-the-art asteroseismological model of G117-B15A, the archetype of the H-rich atmosphere (DA) white dwarf pulsators (also known as DAV or ZZ Ceti variables), and use the most recently measured value of the rate of period change for the dominant mode of this pulsating star to derive a new constraint on the mass of axion, the still conjectural non-barionic particle considered as candidate for dark matter of the Universe. Assuming that G117-B15A is truly represented by our asteroseismological model, and in particular, that the period of the dominant mode is associated to a pulsation g-mode trapped in the H envelope, we find strong indications of the existence of extra cooling in this star, compatible with emission of axions of mass m_a \cos^2 \beta = 17.4^{+2.3}_{-2.7} meV.

We present a large bank of chemical profiles and pulsation periods suited for asteroseismological studies of ZZ Ceti (or DAV) variable stars. Our background equilibrium DA white dwarf models are the result of fully evolutionary computations that take into account the complete history of the progenitor stars from the ZAMS. The models are characterized by self-consistent chemical structures from the centre to the surface, and cover a wide range of stellar masses, effective temperatures, and H envelope thicknesses. We present dipole and quadrupole pulsation $g$-mode periods comfortably covering the interval of periods observed in ZZ Ceti stars. Complete tabulations of chemical profiles and pulsation periods to be used in asteroseismological period fits, as well as other quantities of interest, can be freely downloaded from our website (\url{http://www.fcaglp.unlp.edu.ar/evolgroup}

DBV stars are pulsating white dwarfs with atmospheres rich in He. Asteroseismology of DBV stars can provide valuable clues about the origin, structure and evolution of hydrogen-deficient white dwarfs, and may allow to study neutrino and axion physics. Recently, a new DBV star, KIC 8626021, has been discovered in the field of the \emph{Kepler} spacecraft. It is expected that further monitoring of this star in the next years will enable astronomers to determine its detailed asteroseismic profile. We perform an asteroseismological analysis of KIC 8626021 on the basis of fully evolutionary DB white-dwarf models. We employ a complete set of evolutionary DB white-dwarf structures covering a wide range of effective temperatures and stellar masses. They have been obtained on the basis of a complete treatment of the evolutionary history of progenitors stars. We compute g-mode adiabatic pulsation periods for this set of models and compare them with the pulsation properties exhibited by KIC 8626021. On the basis of the mean period spacing of the star, we found that the stellar mass should be substantially larger than spectroscopy indicates. From period-to-period fits we found an asteroseismological model characterized by an effective temperature much higher than the spectroscopic estimate. In agreement with a recent asteroseismological analysis of this star by other authors, we conclude that KIC 8626021 is located near the blue edge of the DBV instability strip, contrarily to spectroscopic predictions. We also conclude that the mass of KIC 8626021 should be substantially larger than thought.

Asteroseismology of DAV stars (ZZ Ceti variables) can provide valuable clues about the origin, structure and evolution of DA (atmospheres rich in H) white dwarfs. Recently, a new DAV star, WD J191643.83+393849.7, has been discovered in the field of the Kepler spacecraft. It is expected that further monitoring of this star in the next years will enable astronomers to obtain the best lightcurve of a pulsating DA white dwarf ever recorded, and thus to know with unprecedented precision the hidden details of the internal structure of this star. In this paper, we perform a first asteroseismological analysis of WD J191643.83+393849.7 on the basis of fully evolutionary DA white-dwarf models. Specifically, we employ a complete set of evolutionary DA white-dwarf structures covering a wide range of effective temperatures, stellar masses, and H envelope thicknesses. These models have been obtained on the basis of a complete treatment of the evolutionary history of progenitors stars. We compute g-mode adiabatic pulsation periods for this set of models and compare them with the pulsation periods exhibited by WD J191643.83+393849.7. Based on a tentative estimation of the mean period spacing of the star, we find that the stellar mass should be substantially large ($\gtrsim 0.80 M_{\odot}$), in agreement with the spectroscopically derived stellar mass. Also, from period-to-period fits we find an asteroseismological model characterised by a low effective temperature, rather high stellar mass and a thin H envelope. The possibility that this rather massive pulsating white dwarf can be further monitored with Kepler with a high degree of detail turns the star WD J191643.83+393849.7 into a promising and unique object to study the physics of crystallization and carbon/oxygen phase diagrams at high densities.

We argue that NSV 11749, an eruption observed in the early twentieth century, was a rare event known as "very late thermal pulse" (VLTP). To support our argument we compare the lightcurve of NSV 11749 with those of the two bonafide VLTP objects known to date, V4334 Sgr and V605 Aql, and with those predicted by state of the art stellar evolution models. Next, we explore the IPHAS and 2MASS catalogues for possible counterparts of the eruption. Our analysis shows that the VLTP scenario outperforms all other proposed scenarios as an explanation of NSV 11749. We identify an IPHAS/2MASS source at the eruption location of NSV 11749. The derived colors suggest that the object is not enshrouded in a thick dust shell as V605 Aql and V4334 Sgr. Also the absence of an apparent planetary nebula (PN) at the eruption location suggests differences with known VLTP objects which might be linked to the intensity of the eruption and the mass of the object. Further exploration of this source and scenario seems desirable. If NSV 11749 was a born again star, it would be the third event of its kind to have been observed and will strongly help us to increase our understanding on the later stages of stellar evolution and violent reactive convective burning.

We present results of an ongoing study of the pulsational properties of sdB stellar models aimed at exploring the consequences of detailed chemical transitions for radial, $p-$ and $g-$ modes. In particular, we focus on the effects of diffusion at the H-He transition and of He-burning at the convective cores. We find that diffusion of He and H has a strong impact on the period spectrum of sdBVs stars, leading to less efficient mode trapping. Our results also suggests that asteroseismology of sdBVs stars might offer a very good opportunity to constrain extramixing processes in the He-burning cores of horizontal branch stars.

ZZ Ceti stars form the most numerous group of degenerate variable stars. They are otherwise normal DA (H-rich atmospheres) white dwarfs that exhibit pulsations. Here, we present an asteroseismological analysis for 44 bright ZZ Ceti stars based on a new set of fully evolutionary DA white dwarf models characterized by detailed chemical profiles from the centre to the surface. One of our targets is the archetypal ZZ Ceti star G117$-$B15A, for which we obtain an asteroseismological model with an effective temperature and a surface gravity in excellent agreement with the spectroscopy. The asteroseismological analysis of a set of 44 ZZ Ceti stars has the potential to characterize the global properties of the class, in particular the thicknesses of the hydrogen envelope and the stellar masses. Our results support the belief that white dwarfs in the solar neighbourhood harbor a broad range of hydrogen-layer thickness.

Cool white dwarfs are reliable and independent stellar chronometers. The most common white dwarfs have carbon-oxygen dense cores. Consequently, the cooling ages of very cool white dwarfs sensitively depend on the adopted phase diagram of the carbon-oxygen binary mixture. A new phase diagram of dense carbon-oxygen mixtures appropriate for white dwarf interiors has been recently obtained using direct molecular dynamics simulations. In this paper, we explore the consequences of this phase diagram in the evolution of cool white dwarfs. To do this we employ a detailed stellar evolutionary code and accurate initial white dwarf configurations, derived from the full evolution of progenitor stars. We use two different phase diagrams, that of Horowitz et al. (2010), which presents an azeotrope, and the phase diagram of Segretain & Chabrier (1993), which is of the spindle form. We computed the evolution of 0.593 and 0.878M_sun white dwarf models during the crystallization phase, and we found that the energy released by carbon-oxygen phase separation is smaller when the new phase diagram of Horowitz et al. (2010) is used. This translates into time delays that are on average a factor about 2 smaller than those obtained when the phase diagram of Segretain & Chabrier (1993) is employed. Our results have important implications for white dwarf cosmochronology, because the cooling ages of very old white dwarfs are different for the two phase diagrams. This may have a noticeable impact on the age determinations of very old globular clusters, for which the white dwarf color-magnitude diagram provides an independent way of estimating their age.

We investigate the pulsation driving mechanism responsible for the long-period photometric variations observed in LS IV-14{\deg}116, a subdwarf B star showing a He-enriched atmospheric composition. To this end, we perform detailed nonadiabatic pulsation computations over fully evolutionary post He-core-flash stellar structure models, appropriate for hot subdwarf stars at evolutionary phases previous to the He-core burning stage. We found that the variability of LS IV-14{\deg}116 can be attributed to nonradial g-mode pulsations excited by the {\epsilon}-mechanism acting in the He-burning shells that appear before the star settles on the He-core burning stage. Even more interestingly, our results show that LS IV-14{\deg}116 could be the first known pulsating star in which the {\epsilon}-mechanism of mode excitation is operating. Last but not least, we find that the period range of destabilized modes is sensitive to the exact location of the burning shell, something that might help to distinguish between the different evolutionary scenarios proposed for the formation of this star.

We perform an asteroseismological study on the DAV star G117-B15A on the basis of a modern set of fully evolutionary DA white dwarf models that have consistent chemical profiles at the core and the envelope. We found an asteroseismological model for G117-B15A that closely reproduces its observed pulsation periods. Then, we use the most recently measured value of the rate of period change for the dominant mode of this pulsating star to impose a preliminary upper limit to the mass of the axion.

We investigate the internal rotation profile of the GW Vir (PG1159-type) star PG 0122+200 by employing an asteroseismological model that closely reproduces the observed periods of this star. We adopt a forward approach and two inversion methods based on the rotational splitting of the pulsation frequencies to explore the properties of the rotation of PG 0122+200. We found evidence for differential rotation in this star.

We put asteroseismological constraints on the internal rotation profile of the GW Vir (PG1159-type) star PG 0122+200. To this end we employ a state-of-the-art asteroseismological model for this star and we assess the expected frequency splittings induced by rotation adopting a forward approach in which we compare the theoretical frequency separations with the observed ones assuming different types of plausible internal rotation profiles. We also employ two asteroseismological inversion methods for the inversion of the rotation profile of PG 0122+200. We find evidence for differential rotation in this star. We demonstrate that the frequency splittings of the rotational multiplets exhibited by PG 0122+200 are compatible with a rotation profile in which the central regions are spinning about 2.4 times faster than the stellar surface.

We propose a scenario for the formation of DA white dwarfs with very thin helium buffers. For these stars we explore the possible occurrence of diffusion-induced CNO- flashes, during their early cooling stage. In order to obtain very thin helium buffers, we simulate the formation of low mass remnants through an AGB final/late thermal pulse (AFTP/LTP scenario). Then we calculate the consequent white dwarf cooling evolution by means of a consistent treatment of element diffusion and nuclear burning. Based on physically sounding white dwarf models, we find that the range of helium buffer masses for these diffusion-induced novas to occur is significantly smaller than that predicted by the only previous study of this scenario. As a matter of fact, we find that these flashes do occur only in some low-mass (M < 0.6M) and low metallicity (Z_ZAMS <0.001) remnants about 10^6 - 10^7 yr after departing from the AGB. For these objects, we expect the luminosity to increase by about 4 orders of magnitude in less than a decade. We also show that diffusion-induced novas should display a very typical eruption lightcurve, with an increase of about a few magnitudes per year before reaching a maximum of M_V ~ -5 to -6. Our simulations show that surface abundances after the outburst are characterized by logNH/NHe ~ -0.15...0.6 and N>C>O by mass fractions. Contrary to previous speculations we show that these events are not recurrent and do not change substantially the final H-content of the cool (DA) white dwarf. (Abridged)

We compute full evolutionary sequences of red giant branch stars close to the luminosity bump by including state of the art composition transport prescriptions for the thermohaline mixing regimes. In particular we adopt a self-consistent double-diffusive convection theory, that allows to handle the instabilities that arise when thermal and composition gradients compete against each other, and a very recent empirically motivated and parameter free asymptotic scaling law for thermohaline composition transport. In agreement with previous works, we find that during the red giant stage, a thermohaline instability sets in shortly after the hydrogen burning shell (HBS) encounters the chemical discontinuity left behind by the first dredge-up. We also find that the thermohaline unstable region, initially appearing at the exterior wing of the HBS, is unable to reach the outer convective envelope, with the consequence that no mixing of elements that produces a non-canonical modification of the stellar surface abundances occurs. Also in agreement with previous works, we find that by artificially increasing the mixing efficiency of thermohaline regions it is possible to connect both unstable regions, thus affecting the photospheric composition. However, we find that in order to reproduce the observed abundances of red giant branch stars close to the luminosity bump, thermohaline mixing efficiency has to be artificially increased by about 4 orders of magnitude from that predicted by recent 3D numerical simulations of thermohaline convection close to astrophysical environments. From this we conclude the chemical abundance anomalies of red giant stars cannot be explained on the basis of thermohaline mixing alone.

We compute new chemical profiles for the core and envelope of white dwarfs appropriate for pulsational studies of ZZ Ceti stars. These profiles are extracted from the complete evolution of progenitor stars, evolved through the main sequence and the thermally-pulsing asymptotic giant branch (AGB) stages, and from time-dependent element diffusion during white dwarf evolution. We discuss the importance of the initial-final mass relationship for the white dwarf carbon-oxygen composition. In particular, we find that the central oxygen abundance may be underestimated by about 15% if the white dwarf mass is assumed to be the hydrogen-free core mass before the first thermal pulse. We also discuss the importance for the chemical profiles expected in the outermost layers of ZZ Ceti stars of the computation of the thermally-pulsing AGB phase and of the phase in which element diffusion is relevant. We find a strong dependence of the outer layer chemical stratification on the stellar mass. In particular, in the less massive models, the double-layered structure in the helium layer built up during the thermally-pulsing AGB phase is not removed by diffusion by the time the ZZ Ceti stage is reached. Finally, we perform adiabatic pulsation calculations and discuss the implications of our new chemical profiles for the pulsational properties of ZZ Ceti stars. We find that the whole $g-$mode period spectrum and the mode-trapping properties of these pulsating white dwarfs as derived from our new chemical profiles are substantially different from those based on chemical profiles widely used in existing asteroseismological studies. Thus, we expect the asteroseismological models derived from our chemical profiles to be significantly different from those found thus far.

We present full evolutionary calculations appropriate for the study of hydrogen-rich DA white dwarfs. This is done by evolving white dwarf progenitors from the zero age main sequence, through the core hydrogen burning phase, the helium burning phase and the thermally pulsing asymptotic giant branch phase to the white dwarf stage. Complete evolutionary sequences are computed for a wide range of stellar masses and for two different metallicities: Z=0.01, which is representative of the solar neighborhood, and Z=0.001, which is appropriate for the study of old stellar systems, like globular clusters. During the white dwarf cooling stage we compute self-consistently the phase in which nuclear reactions are still important, the diffusive evolution of the elements in the outer layers and, finally, we also take into account all the relevant energy sources in the deep interior of the white dwarf, like the release of latent heat and the release of gravitational energy due to carbon-oxygen phase separation upon crystallization. We also provide colors and magnitudes for these sequences, based on a new set of improved non-gray white dwarf model atmospheres, which include the most up-to-date physical inputs like the Lyman alpha quasi-molecular opacity. The calculations are extended down to an effective temperature of 2,500 K. Our calculations provide a homogeneous set of evolutionary cooling tracks appropriate for mass and age determinations of old DA white dwarfs and for white dwarf cosmochronology of the different Galactic populations.

We present full evolutionary calculations appropriate for the study of hot hydrogen-deficent DO white dwarfs, PG 1159 stars, and DB white dwarfs. White dwarf sequences are computed for a wide range of stellar masses and helium envelopes on the basis of a complete treatment of the evolutionary history of progenitors stars, including the core hydrogen and helium burning phases, the thermally-pulsing AGB phase, and the born-again episode that is responsible for the hydrogen deficiency. We also provide colors and magnitudes for the new sequences for $T_{\rm eff} < 40 000$ K, where the NLTE effects are not dominant. These new calculations provide an homogeneous set of evolutionary tracks appropriate for mass and age determinations for both PG 1159 stars and DO white dwarfs. The calculations are extended down to an effective temperature of 7 000 K. We applied these new tracks to redetermine stellar masses and ages of all known DO white dwarfs with spectroscopically-determined effective temperatures and gravities, and compare them with previous results. We also compare for the first time consistent mass determinations for both DO and PG 1159 stars, and find a considerably higher mean mass for the DO white dwarfs. We discuss as well the chemical profile expected in the envelope of variable DB white dwarfs from the consideration of the evolutionary history of progenitor stars. Finally, we present tentative evidence for a different evolutionary channel, other than that involving the PG 1159 stars, for the formation of hot, hydrogen-deficient white dwarfs.

We present a pulsational stability analysis of hot post-AGB H-deficient pre-white dwarf stars with active He-burning shells. The stellar models employed are state-of-the-art equilibrium structures representative of PG1159 stars derived from the complete evolution of the progenitor stars. On the basis of fully nonadiabatic pulsation computations, we confirmed theoretical evidence for the existence of a separate PG1159 instability strip in the $\log T_{\rm eff} - \log g$ diagram characterized by short-period $g$-modes excited by the $\epsilon$-mechanism. This instability strip partially overlaps the already known GW Vir instability strip of intermediate/long period $g$-modes destabilized by the classical $\kappa$-mechanism acting on the partial ionization of C and/or O in the envelope of PG1159 stars. We found that PG1159 stars characterized by thick He-rich envelopes and located inside this overlapping region could exhibit both short and intermediate/long periods simultaneously. we study the particular case of VV 47, a pulsating planetary nebula nucleus that has been reported to exhibit a series of unusually short pulsation periods. We found that the long periods exhibited by VV 47 can be readily explained by the classical $\kappa$-mechanism, while the observed short-period branch below $\approx 300$ s could correspond to modes triggered by the He-burning shell through the $\epsilon$-mechanism, although more observational work is needed to confirm the reality of these short-period modes. Were the existence of short-period $g$-modes in this star convincingly confirmed by future observations, VV 47 could be the first known pulsating star in which both the $\kappa$-mechanism and the $\epsilon$-mechanism of mode driving are simultaneously operating.

Motivated by the recent detection of single and binary He-core white dwarfs in metal-rich clusters, we present a full set of evolutionary calculations and colors appropriate for the study of such white dwarfs. The paper is also aimed at investigating whether stable hydrogen burning may constitute a main source of energy for massive He-core white dwarfs resulting from high-metallicity progenitors. White dwarf sequences are derived by taking into account the evolutionary history of progenitor stars with supersolar metallicities. We also incorporate a self-consistent, time-dependent treatment of gravitational settling and chemical diffusion, as well as of the residual nuclear burning. We find that the influence of residual nuclear burning during the late stages of white dwarf evolution is strongly dependent on the occurrence of chemical diffusion at the base of the hydrogen-rich envelope. When no diffusion is considered, residual hydrogen burning strongly influences the advanced stages of white dwarf cooling, introducing evolutionary delays of several Gyr. By contrast, when diffusion is taken into account the role of residual nuclear burning is strongly mitigated, and the evolution is dictated only by the thermal content stored in the ions. In addition, for all of our sequences, we provide accurate color and magnitudes on the basis of new and improved non gray model atmospheres which explicitly include Ly$\alpha$ quasi-molecular opacity.

A seismological study of the pulsating PG1159 stars PG2131+066 and PG 1707+427 is presented. We perform extensive adiabatic computations of g-mode pulsation periods of PG1159 evolutionary models with stellar masses ranging from 0.530 to 0.741 Msun. We constrain the stellar mass of PG2131+066 and PG 1707+427 by comparing the observed period spacing of each star with the theoretical asymptotic period spacings and with the average of the computed period spacings. We also employ the individual observed periods to find representative seismological models for both stars.

We present an asteroseismological study on the two high-gravity pulsating PG1159 (GW Vir or DOV) stars, PG 2131+066 and PG 1707+427, and on the pulsating [WCE] star NGC 1501. We compute adiabatic $g$-mode pulsation periods on PG1159 evolutionary models with stellar masses ranging from 0.530 to 0.741 Msun. These models take into account the complete evolution of progenitor stars, through the thermally pulsing AGB phase, and born-again episode. We constrain the stellar mass of PG 2131+066, PG 1707+427 and NGC 1501 by comparing the observed period spacing with the asymptotic period spacing and with the average of the computed period spacings. We also employ the individual observed periods in search of representative seismological models for each star. This work closes our short series of asteroseismological studies on pulsating pre-white dwarf stars. Our results demonstrate the usefulness of asteroseismology for probing the internal structure and evolutionary status of pre-white dwarf stars. In particular, asteroseismology is able to determine stellar masses of PG1159 stars with an accuracy comparable or even better than spectroscopy.

We present the first full evolutionary calculations aimed at exploring the origin of hot DQ white dwarfs. These calculations consistently cover the whole evolution from the born-again stage to the white dwarf cooling track. Our calculations provide strong support to the diffusive/convective-mixing picture for the formation of hot DQs. We find that the hot DQ stage is a short-lived stage and that the range of effective temperatures where hot DQ stars are found can be accounted for by different masses of residual helium and/or different initial stellar masses. In the frame of this scenario, a correlation between the effective temperature and the surface carbon abundance in DQs should be expected, with the largest carbon abundances expected in the hottest DQs. From our calculations, we suggest that most of the hot DQs could be the cooler descendants of some PG1159 stars characterized by He-rich envelopes markedly smaller than those predicted by the standard theory of stellar evolution. At least for one hot DQ, the high-gravity white dwarf SDSS J142625.70+575218.4, an evolutionary link between this star and the massive PG1159 star H1504+65 is plausible.

We present 1D numerical simulations aimed at studying the hot-flasher scenario for the formation of He-rich subdwarf stars. Sequences were calculated for a wide range of metallicities and physical assumptions, such as the stellar mass at the moment of the helium core flash. This allows us to study the two previously proposed flavors of the hot-flasher scenario ("deep" and "shallow" mixing cases) and to identify a third transition type. Our sequences are calculated by solving simultaneously the mixing and burning equations within a diffusive convection picture, and in the context of standard mixing length theory. We are able to follow chemical evolution during deep-mixing events in which hydrogen is burned violently, and therefore able to present a homogeneous set of abundances for different metallicities and varieties of hot-flashers. We extend the scope of our work by analyzing the effects of non-standard assumptions, such as the effect of chemical gradients, extra-mixing at convective boundaries, possible reduction in convective velocities, or the interplay between difussion and mass loss. Particular emphasis is placed on the predicted surface properties of the models. We find that the hot-flasher scenario is a viable explanation for the formation and surface properties of He-sdO stars. Our results also show that, during the early He-core burning stage, element diffusion may produce the transformation of (post hot-flasher) He-rich atmospheres into He-deficient ones. If this is so, then we find that He-sdO stars should be the progenitors of some of the hottest sdB stars.

We present 1D numerical simulations aimed at studying the hot-flasher scenario for the formation of He-rich subdwarf stars. Sequences were calculated for a wide range of metallicities and with the He core flash at different points of the post-RGB evolution (i.e. different remnant masses). We followed the complete evolution from the ZAMS, through the hot-flasher event, and to the subdwarf stage for all kinds of hot-flashers. This allows us to present a homogeneous set of abundances for different metallicities and all flavors of hot-flashers. We extend the scope of our work by analyzing the effects in the predicted surface abundances of some standard assumptions in convective mixing and the effects of element diffusion. We find that the hot-flasher scenario is a viable explanation for the formation of He-sdO stars. Our results also show that element diffusion may produce the transformation of (post hot-flasher) He-rich atmospheres into He-deficient ones. If this is so, then the hot-flasher scenario is able to reproduce both the observed properties and distribution of He-sdO stars.

Using the most recent results about white dwarfs in 10 open clusters, we revisit semi-empirical estimates of the initial-final mass relation in star clusters, with emphasis on the use of stellar evolution models. We discuss the influence of these models on each step of the derivation. One intention of our work is to use consistent sets of calculations both for the isochrones and the white dwarf cooling tracks. The second one is to derive the range of systematic errors arising from stellar evolution theory. This is achieved by using different sources for the stellar models and by varying physical assumptions and input data. We find that systematic errors, including the determination of the cluster age, are dominating the initial mass values, while observational uncertainties influence the final mass primarily. After having determined the systematic errors, the initial-final mass relation allows us finally to draw conclusions about the physics of the stellar models, in particular about convective overshooting.

We reexamine the theoretical instability domain of pulsating DB white dwarfs (DBV or V777 Her variables). We performed an extensive $g$-mode nonadiabatic pulsation analysis of DB evolutionary models considering a wide range of stellar masses, for which the complete evolutionary stages of their progenitors from the ZAMS, through the thermally pulsing AGB and born-again phases, the domain of the PG1159 stars, the hot phase of DO white dwarfs, and then the DB white dwarf stage have been considered. We explicitly account for the evolution of the chemical abundance distribution due to time-dependent chemical diffusion processes. We examine the impact of the different prescriptions of the MLT theory of convection and the effects of small amounts of H in the almost He-pure atmospheres of DB stars on the precise location of the theoretical blue edge of the DBV instability strip.

We present evidence that PG1159 stars could harbour He--rich envelopes substantially thinner than those predicted by current evolutionary models with current estimates of mass loss, which may be attributable to an extensive mass--loss episode during the born--again AGB phase. Specifically, we show that the models with thin He--rich envelopes predict remarkably large magnitudes of the rates of period change of the trapped and untrapped modes observed in the pulsating star PG 1159$-$035. This is a consequence of the much shorter evolutionary timescale of the models with thin He--rich envelopes during the low--gravity PG1159 regime. Our findings are particularly interesting in view of the suggestion of an evolutionary link between the helium--deficient PG1159 star H1504+65 and the recently discovered white dwarfs with almost pure carbon atmospheres.

An asteroseismological study of PG 1159-035, the prototype of the GW Vir variable stars, has been performed on the basis of detailed and full PG1159 evolutionary models presented by Miller Bertolami & Althaus (2006). We carried out extensive computations of adiabatic g-mode pulsation periods on PG1159 evolutionary models with stellar masses spanning the range 0.530 to 0.741 Mo. We derive a stellar mass in the range 0.56-0.59 Mo from the period-spacing data alone. We also find, on the basis of a period-fit procedure, a seismic model representative of PG 1159-035 that reproduces the observed period pattern with an average of the period differences of 0.64-1.03 s, consistent with the expected model uncertainties. The results of the period-fit analysis carried out in this work suggest that the surface gravity of PG 1159-035 would be 1 sigma larger than the spectroscopically inferred gravity. For our best-fit model of PG 1159-035, all of the pulsation modes are characterized by positive rates of period changes, at odds with the measurements by Costa & Kepler (2007).

We present numerical simulations of the very late thermal pulse (VLTP) scenario for a wide range of remnant masses. We show that, by taking into account the different possible remnant masses, the fast outburst evolution of V4334 Sgr (a.k.a. Sakurai's Object) and V605 Aql can be reproduced within standard 1D stellar evolution models. A dichotomy in the born again timescales is found, with lower mass remnants evolving in a few years and higher mass remnants (M $\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$}}}\hbox{$>$}}$ 0.6 M$_{\odot}$) failing to expand due to the H-flash and, as a consequence, evolving in timescales typical of He-shell flash driven born agains ($\sim 100$ yr).

We present first results from evolutionary simulations aimed at exploring the Hot-Flasher scenario for the formation of H-deficient subdwarf stars. The two types of late hot flashers that lead to He-enriched surfaces, "deep" and "shallow" mixing cases, are investigated for different metallicities.

We analyze systematics in the asteroseismological mass determination methods in pulsating PG 1159 stars. We compare the seismic masses resulting from the comparison of the observed mean period spacings with the usually adopted asymptotic period spacings, and the average of the computed period spacings. Computations are based on full PG1159 evolutionary models with stellar masses ranging from 0.530 to 0.741 Mo that take into account the complete evolution of progenitor stars. We conclude that asteroseismology is a precise and powerful technique that determines the masses to a high internal accuracy, but it depends on the adopted mass determination method. In particular, we find that in the case of pulsating PG 1159 stars characterized by short pulsation periods, like PG 2131+066 and PG 0122+200, the employment of the asymptotic period spacings overestimates the stellar mass by about 0.06 Mo as compared with inferences from the average of the period spacings. In this case, the discrepancy between asteroseismological and spectroscopical masses is markedly reduced when use is made of the mean period spacing instead of the asymptotic period spacing.

We present an asteroseismological study on PG 0122+200, the coolest known pulsating PG1159 (GW Vir) star. Our results are based on an augmented set of the full PG1159 evolutionary models recently presented by Miller Bertolami & Althaus (2006). We perform extensive computations of adiabatic g-mode pulsation periods on PG1159 evolutionary models with stellar masses ranging from 0.530 to 0.741 Msun. We derive a stellar mass of 0.626 Msun from a comparison between the observed period spacing and the computed asymptotic period spacing, and a stellar mass of 0.567 Msun by comparing the observed period spacing with the average of the computed period spacing. We also find, on the basis of a period-fit procedure, an asteroseismological model representative of PG 0122+200 which is able to reproduce the observed period pattern with an average of the period differences of 0.88 s. The model has an effective temperature of 81500 K, a stellar mass of 0.556 Msun, a surface gravity log g= 7.65, a stellar luminosity and radius of log(L/Lsun)= 1.14 and log(R/Rsun)= -1.73, respectively, and a He-rich envelope thickness of Menv= 0.019 Msun. We derive a seismic distance of about 614 pc and a parallax of about 1.6 mas. The results of the period-fit analysis carried out in this work suggest that the asteroseismological mass of PG 0122+200 could be 6-20 % lower than thought hitherto and in closer agreement (to within 5 %) with the spectroscopic mass. This result suggests that a reasonable consistency between the stellar mass values obtained from spectroscopy and asteroseismology can be expected when detailed PG1159 evolutionary models are considered.

We present 1-D numerical simulations of the very late thermal pulse (VLTP) scenario for a wide range of remnant masses. We show that by taking into account the different possible remnant masses, the observed evolution of V4334 Sgr (a.k.a. Sakurai's Object) can be reproduced within the standard 1D-MLT stellar evolutionary models without the inclusion of any $ad-hoc$ reduced mixing efficiency. Our simulations hint at a consistent picture with present observations of V4334 Sgr. From energetics, and within the standard MLT approach, we show that low mass remnants \hbox{($M\lesssim0.6$\msun)} are expected to behave markedly different than higher mass remnants \hbox{($M\gtrsim0.6$\msun)} in the sense that the latter are not expected to expand significantly as a result of the violent H-burning that takes place during the VLTP. We also assess the discrepancy in the born again times obtained by different authors by comparing the energy that can be liberated by H-burning during the VLTP event.

We analyze the robustness of H--deficient post--AGB tracks regarding previous evolution of their progenitor stars and the constitutive physics of the remnants. Our motivation is a recent suggestion of Werner & Herwig (2006) that previous evolution should be important in shaping the final post--AGB track and the persisting discrepancy between asteroseismological and spectroscopical mass determinations. This work is thus complementary to our previous work (Miller Bertolami & Althaus 2006) and intends to shed some light on the uncertainty behind the evolutionary tracks presented there. We compute full evolutionary models for PG1159 stars taking into account different extramixing (overshooting) efficiencies and lifetimes on the TP-AGB during the progenitor evolution. We also assess the effect of possible differences in the opacities and equation of state by artificially changing them before the PG1159 stage. Also comparisons are made with the few H-deficient post--AGB tracks available in the literature. Contrary to our expectations, we found that previous evolution is not a main factor in shaping H--deficient post--AGB tracks. Interestingly enough, we find that only an increase of $\sim50%$ in the intershell opacities at high effective temperatures may affect the tracks as to reconcile spectroscopic and asteroseismologic mass determinations. This forces us to conclude that our previous tracks (Miller Bertolami & Althaus 2006) are robust enough as to be used for spectroscopic mass determinations, unless opacities in the intershell region are substantially different. Our results, then, call for an analysis of possible systematics in the usually adopted asteroseismological mass determination methods.

We examine a recently-proposed evolutionary scenario that could explain the existence of the low-mass, helium-enriched PG1159 stars. We assess the overstability of pulsation g-modes of stellar models as evolution proceeds along the PG1159 domain. Stellar models are extracted from the full evolution of a 1 Mo model star that experiences its first thermal pulse as a late thermal pulse (LTP) after leaving the AGB. The evolutionary stages corresponding to the born-again episode and the subsequent helium sub-flashes are taken into account in detail. Under reasonable mass-loss rate assumptions, the evolutionary scenario reproduces the high helium abundances observed in some PG1159 stars. We find that, despite the high helium abundance in the driving layers, there exists a narrow region in the log g-log Teff diagram for which the helium-enriched PG1159 sequence exhibits unstable pulsation modes with periods in the range 500 to 1600 s. In particular, the nonpulsating helium-enriched PG1159 star, MCT 0130-1937, is located outside the theoretical instability domain. Our results suggest that MCT 0130-1937 is a real non-pulsating star and that the lack of pulsations should not be attributed to unfavorable geometry. Our study hints at a consistent picture between the evolutionary scenario that could explain the existence of helium-enriched PG1159 stars and the nonvariable nature of MCT 0130-1937. We also present theoretical support for the unusually high helium abundance observed in the nonpulsating PG1159 star HS 1517+7403. We suggest that HS 1517+7403 could be a transition object linking the low-mass helium-rich O(He) stars with the helium-enriched PG1159 stars via the evolutionary connection K1-27-->HS 1517+7403-->MCT 0130-1937.

We present asteroseismological inferences on RX J2117.1+3412, the hottest known pulsating PG1159 star. Our results are based on full PG1159 evolutionary models recently presented by Miller Bertolami & Althaus (2006). We performed extensive computations of adiabatic g-mode pulsation periods on PG1159 evolutionary models with stellar masses ranging from 0.530 to 0.741 Mo. PG1159 stellar models are extracted from the complete evolution of progenitor stars started from the ZAMS, through the thermally pulsing AGB and born-again phases to the domain of the PG 1159 stars. We constrained the stellar mass of RX J2117.1+3412 by comparing the observed period spacing with the asymptotic period spacing and with the average of the computed period spacings. We also employed the individual observed periods to find a representative seismological model. We derive a stellar mass of 0.56-0.57 Mo from the period spacing data alone. In addition, we found a best-fit model representative for RX J2117.1+3412 with an effective temperature of 163,400 K, a stellar mass of 0.565 Mo, and a surface gravity log g= 6.61. The derived stellar luminosity and radius are log(L/Lo)= 3.36 and log(R/Ro)= -1.23, respectively, and the He-rich envelope thickness is Menv= 0.02 Mo. We derive a seismic distance of 452 pc and a linear size of the planetary nebula of 1.72 pc. These inferences seem to solve the discrepancy between the RX J2117.1+3412 evolutionary timescale and the size of the nebula. All of the seismological tools we use concur to the conclusion that RX J2117.1+3412 must have a stellar mass of 0.565 Mo much in agreement with recent asteroseismology studies and in clear conflict with the predictions of spectroscopy plus evolutionary tracks.

We reexamine the theoretical instability domain of pulsating PG1159 stars (GW Vir variables). We performed an extensive g-mode stability analysis on PG1159 evolutionary models with stellar masses ranging from 0.530 to 0.741 Mo for which the complete evolutionary stages of their progenitors from the ZAMS, through the thermally pulsing AGB and born-again phases to the domain of the PG1159 stars have been considered. We found that pulsations in PG1159 stars are excited by the kappa-mechanism due to partial ionization of carbon and oxygen, and that no composition gradients are needed between the surface layers and the driving region, much in agreement with previous studies. We show, for the first time, the existence of a red edge of the instability strip at high luminosities. We found that all of the GW Vir stars lay within our theoretical instability strip. Our results suggest a qualitative good agreement between the observed and the predicted ranges of unstable periods of individual stars. Finally, we found that generally the seismic masses (derived from the period spacing) of GW Vir stars are somewhat different from the masses suggested by evolutionary tracks coupled with spectroscopy. Improvements in the evolution during the thermally pulsing AGB phase and/or during the core helium burning stage and early AGB could help to alleviate the persisting discrepancies.

We present full evolutionary calculations appropriate to post-AGB PG1159 stars for a wide range of stellar masses. We take into account the complete evolutionary stages of PG1159 progenitors starting from the Zero Age Main Sequence. We consider the two kinds of Born Again Scenarios, the very late thermal pulse (VLTP) and the late thermal pulse (LTP), that give rise to hydrogen-deficient compositions. The location of our PG1159 tracks in the effective temperature - gravity diagram and their comparison with previous calculations as well as the resulting surface compositions are discussed at some length. Our results reinforce the idea that the different abundances of $^{14}$N observed at the surface of those PG1159 stars with undetected hydrogen is an indication that the progenitors of these stars would have evolved through a VLTP episode, where most of the hydrogen content of the remnant is burnt, or LTP, where hydrogen is not burnt but instead diluted to very low surface abundances. We derive new values for spectroscopical masses based on these new models. We discuss the correlation between the presence of planetary nebulae and the $^{14}$N abundance as another indicator that $^{14}$N-rich objects should come from a VLTP episode while $^{14}$N-deficient ones should be the result of a LTP. Finally, we discuss an evolutionary scenario that could explain the existence of PG1159 stars with unusually high helium abundances and a possible evolutionary connection between these stars and the low mass O(He) stars.

We present evolutionary calculations aimed at describing the born-again scenario for post-AGB remnant stars of 0.5842 and 0.5885 \msun. Results are based on a detailed treatment of the physical processes responsible for the chemical abundance changes. We considered two theories of convection: the standard mixing length theory (MLT) and the double-diffusive GNA convection developed by Grossman et al. The latter accounts for the effect of the chemical gradient ($\nabla\mu$) in the mixing processes and in the transport of energy. We also explore the dependence of the born-again evolution on some physical hypothesis, such as the effect of the existence of non-zero chemical gradients, the prescription for the velocity of the convective elements and the size of the overshooting zones. Attention is given to the behavior of the born-again times and to the chemical evolution during the ingestion of protons. We find that in our calculations born again times are dependent on time resolution. In particular when the minimum allowed time step is below $5 \times 10^{-5}$ yr we obtain, with the standard mixing length theory, born again times of 5-10 yr. This is true without altering the prescription for the efficiency of convective mixing during the proton ingestion. On the other hand we find that the inclusion of the chemical gradients in the calculation of the mixing velocity tend to increase the born again times by about a factor of two. In addition we find that proton ingestion can be seriously altered if the occurrence of overshooting is modified by the $\nabla\mu$-barrier at the H-He interface, strongly altering born again times.

We study the formation and evolution of DA white dwarfs, the progenitors of which have experienced a late thermal pulse (LTP) shortly after the departure from the thermally pulsing AGB. To this end, we compute the complete evolution of an initially 2.7 Mo star all the way from the zero-age main sequence to the white dwarf stage. We find that most of the original H-rich material of the post-AGB remnant is burnt during the post-LTP evolution, with the result that, at entering its white dwarf cooling track, the remaining H envelope becomes 10^-6 Mo in agreement with asteroseismological inferences for some ZZ Ceti stars.