We study the evolution of close binary systems in order to account for the existence of the recently observed binary system containing the most massive millisecond pulsar ever detected, PSR J0740+6620, and its ultra-cool helium white dwarf companion. In order to find a progenitor for this object we compute the evolution of several binary systems composed by a neutron star and a normal donor star employing our stellar code. We assume conservative mass transfer. We also explore the effects of irradiation feedback on the system. We find that irradiated models also provide adequate models for the millisecond pulsar and its companion, so both irradiated and non irradiated systems are good progenitors for PSR J0740+6620. Finally, we obtain a binary system that evolves and accounts for the observational data of the system composed by PSR J0740+6620 (i.e. orbital period, mass, effective temperature and inferred metallicity of the companion, and mass of the neutron star) in a time scale smaller than the age of the Universe. In order to reach an effective temperature as low as observed, the donor star should have an helium envelope as demanded by observations.

We analyse the evolution of close binary systems containing a neutron star that lead to the formation of redback pulsars. Recently there has been some debate on the origin of such systems and the formation mechanism of redbacks may still be considered as an open problem. We show that the operation of a strong evaporation mechanism, starting from the moment when the donor star becomes fully convective (or alternatively since the formation of the neutron star by accretion induced collapse), produces systems with donor masses and orbital periods in the range corresponding to redbacks with donors appreciably smaller than their Roche lobes, i.e., they have low filling factors (lower than $0.75$). Models of redback pulsars can be constructed assuming the occurrence of irradiation feedback. They have been shown to undergo cyclic mass transfer during the epoch at which they attain donor masses and orbital periods corresponding to redbacks, and stay in quasi-Roche lobe overflow conditions with {\it high} filling factors. We show that, if irradiation feedback occurs and radio ejection inhibits further accretion onto the neutron star after the first mass transfer cycle, the redback systems feature {\it high} filling factors. We suggest that the filling factor should be considered as a useful tool for discriminating among those redback formation mechanisms. We compare theoretical results with available observations, and conclude that observations tend to favour models with high filling factors.

We study the present evolutionary status of the binary system containing the 2.66 ms pulsar PSR J1417-4402 in a 5.4 day orbit. This is the pulsar in the original source 3FGL J1417.5-4402, that has undergone a transition from X-ray state to a pulsar state, just like some redbacks did. The system has many characteristics similar to redback pulsars family, but is on a much wider orbit. We show that close binary evolution including irradiation feedback driven by the luminosity due to accretion onto the neutron star component of the pair, and evaporation due to pulsar emission, is able to account for the masses of the components and the photometric data of the donor star. The tracks leading to the present PSR J1417-4402 are degenerate within a range of parameters, suggesting that the {\it same} physics invoked to explain the redback/black widows groups leads to the formation of much wider orbit systems, outside the redback region limits.

It is difficult to establish the properties of massive stars that explode as supernovae. The electromagnetic emission during the first minutes to hours after the emergence of the shock from the stellar surface conveys important information about the final evolution and structure of the exploding star. However, the unpredictable nature of supernova events hinders the detection of this brief initial phase. Here we report the serendipitous discovery of a newly born, normal type IIb supernova (SN 2016gkg), which reveals a rapid brightening at optical wavelengths of about 40 magnitudes per day. The very frequent sampling of the observations allowed us to study in detail the outermost structure of the progenitor of the supernova and the physics of the emergence of the shock. We develop hydrodynamical models of the explosion that naturally account for the complete evolution of the supernova over distinct phases regulated by different physical processes. This result suggests that it is appropriate to decouple the treatment of the shock propagation from the unknown mechanism that triggers the explosion.

Binary evolution is investigated as the source for the extreme horizontal branch (EHB) stars in the old and metal rich open cluster NGC~6791. Employing an updated version of our binary stellar evolution code we demonstrate that EHB stars naturally emerge from the common envelope phase. In sum, the binary model reproduces the observed ($T_{\rm{eff}}$, $\log{g}$) and temporal properties of the EHB over-density tied to NGC 6971, without needing an ad-hoc and anomalous mass-loss prescription.

We study the evolution of close binary systems composed of a normal, intermediate mass star and a neutron star considering a chemical composition typical of that present in globular clusters (Z = 0.001). We look for similarities and differences with respect to solar composition donor stars, which we have extensively studied in the past. As a definite example, we perform an application on one of the redbacks located in a globular cluster. We performed a detailed grid of models in order to find systems that represent the so-called redback binary radio pulsar systems with donor star masses between 0.6 and 2.0 solar masses and orbital periods in the range 0.2 - 0.9 days. We find that the evolution of these binary systems is rather similar to those corresponding to solar composition objects, allowing us to account for the occurrence of redbacks in globular clusters, as the main physical ingredient is the irradiation feedback. Redback systems are in the quasi-RLOF state, that is, almost filling their corresponding Roche lobe. During the irradiation cycle the system alternates between semi-detached and detached states. While detached the system appears as a binary millisecond pulsar, called a redback. Circumstellar material, as seen in redbacks, is left behind after the previous semi-detached phase. The evolution of binary radio pulsar systems considering irradiation successfully accounts for, and provides a way for, the occurrence of redback pulsars in low-metallicity environments such as globular clusters. This is the case despite possible effects of the low metal content of the donor star that could drive systems away from redback configuration.

Two recently discovered very luminous supernovae (SNe) present stimulating cases to explore the extents of the available theoretical models. SN 2011kl represents the first detection of a supernova explosion associated with an ultra-long duration gamma ray burst. ASASSN-15lh was even claimed as the most luminous SN ever discovered, challenging the scenarios so far proposed for stellar explosions. Here we use our radiation hydrodynamics code in order to simulate magnetar powered SNe. To avoid explicitly assuming neutron star properties we adopt the magnetar luminosity and spin-down timescale as free parameters of the model. We find that the light curve (LC) of SN 2011kl is consistent with a magnetar power source, as previously proposed, but we note that some amount of 56^Ni (> 0.08 M_sun) is necessary to explain the low contrast between the LC peak and tail. For the case of ASASSN-15lh we find physically plausible magnetar parameters that reproduce the overall shape of the LC provided the progenitor mass is relatively large (a mass of the ejecta approx 6 M_sun). The ejecta hydrodynamics of this event is dominated by the magnetar input, while the effect is more moderate for SN 2011kl. We conclude that a magnetar model may be used for the interpretation of these events and that the hydrodynamic modeling is necessary to derive the properties of powerful magnetars and their progenitors.

Hubble Space Telescope observations of the site of the supernova (SN) 2008ax obtained in 2011 and 2013 reveal that the possible progenitor object detected in pre-explosion images was in fact multiple. Four point sources are resolved in the new, higher-resolution images. We identify one of the sources with the fading SN. The other three objects are consistent with single supergiant stars. We conclude that their light contaminated the previously identified progenitor candidate. After subtraction of these stars, the progenitor appears to be significantly fainter and bluer than previously measured. Post-explosion photometry at the SN location indicates that the progenitor object has disappeared. If single, the progenitor is compatible with a supergiant star of B to mid-A spectral type, while a Wolf-Rayet (WR) star would be too luminous in the ultraviolet to account for the observations. Moreover, our hydrodynamical modelling shows the pre-explosion mass was $4-5$ $M_\odot$ and the radius was $30-50$ $R_\odot$, which is incompatible with a WR progenitor. We present a possible interacting binary progenitor computed with our evolutionary models that reproduces all the observational evidence. A companion star as luminous as an O9-B0 main-sequence star may have remained after the explosion.

We study the occurrence of delayed SNe~Ia in the single degenerate (SD) scenario. We assume that a massive carbon-oxygen (CO) white dwarf (WD) accretes matter coming from a companion star, making it to spin at the critical rate. We assume uniform rotation due to magnetic field coupling. The carbon ignition mass for non-rotating WDs is M_{ig}^{NR} \approx 1.38 M_{\odot}; while for the case of uniformly rotating WDs it is a few percent larger (M_{ig}^{R} \approx 1.43 M_{\odot}). When accretion rate decreases, the WD begins to lose angular momentum, shrinks, and spins up; however, it does not overflow its critical rotation rate, avoiding mass shedding. Thus, angular momentum losses can lead the CO WD interior to compression and carbon ignition, which would induce an SN~Ia. The delay, largely due to the angular momentum losses timescale, may be large enough to allow the companion star to evolve to a He WD, becoming undetectable at the moment of explosion. This scenario supports the occurrence of delayed SNe~Ia if the final CO WD mass is 1.38 M_{\odot} < M < 1.43 M_{\odot}. We also find that if the delay is longer than ~3 Gyr, the WD would become too cold to explode, rather undergoing collapse.

We report on both high-precision photometry from the MOST space telescope and ground-based spectroscopy of the triple system delta Ori A consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7d, and a more distant tertiary (O9 IV P > 400 yrs). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for 3 weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy we have a well constrained radial velocity curve of Aa1. While we are unable to recover radial velocity variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries which show evidence for tidally induced pulsations.

We study the evolution of close binary systems formed by a normal (solar composition), intermediate mass donor star together with a neutron star. We consider models including irradiation feedback and evaporation. These non-standard ingredients deeply modify the mass transfer stages of these binaries. While models that neglect irradiation feedback undergo continuous, long standing mass transfer episodes, models including these effect suffer a number cycles of mass transfer and detachment. During mass transfer the systems should reveal themselves as low-mass X-ray binaries (LMXBs), whereas when detached they behave as a binary radio pulsars. We show that at these stages irradiated models are in a Roche lobe overflow (RLOF) state or in a quasi-RLOF state. Quasi-RLOF stars have a radius slightly smaller than its Roche lobe. Remarkably, these conditions are attained for orbital period and donor mass values in the range corresponding to a family of binary radio pulsars known as "redbacks". Thus, redback companions should be quasi-RLOF stars. We show that the characteristics of the redback system PSR J1723-2837 are accounted for by these models. In each mass transfer cycle these systems should switch from LMXB to binary radio pulsar states with a timescale of \sim million years. However, there is recent and fast growing evidence of systems switching on far shorter, human timescales. This should be related to instabilities in the accretion disc surrounding the neutron star and/or radio ejection, still to be included in the model having the quasi-RLOF state as a general condition.

We present Hubble Space Telescope (HST) observations of the field of the Type IIb supernova (SN) 2011dh in M51 performed at ~1161 rest-frame days after explosion using the Wide Field Camera 3 and near-UV filters F225W and F336W. A star-like object is detected in both bands and the photometry indicates it has negative (F225W - F336W) color. The observed object is compatible with the companion of the now-vanished yellow supergiant progenitor predicted in interacting binary models. We consider it unlikely that the SN is undergoing strong interaction and thus estimate that it makes a small contribution to the observed flux. The possibilities of having detected an unresolved light echo or an unrelated object are briefly discussed and judged unlikely. Adopting a possible range of extinction by dust, we constrain parameters of the proposed binary system. In particular, the efficiency of mass accretion onto the binary companion must be below 50%, if no significant extinction is produced by newly formed dust. Further multiband observations are required in order to confirm the identification of the object as the companion star. If confirmed, the companion star would already be dominant in the UV/optical regime, so it would readily provide a unique opportunity to perform a detailed study of its properties.

The recent detection in archival HST images of an object at the the location of supernova (SN) iPTF13bvn may represent the first direct evidence of the progenitor of a Type Ib SN. The object's photometry was found to be compatible with a Wolf-Rayet pre-SN star mass of ~11 Msun. However, based on hydrodynamical models we show that the progenitor had a pre-SN mass of ~3.5 Msun and that it could not be larger than ~8 Msun. We propose an interacting binary system as the SN progenitor and perform evolutionary calculations that are able to self-consistently explain the light-curve shape, the absence of hydrogen, and the pre-SN photometry. We further discuss the range of allowed binary systems and predict that the remaining companion is a luminous O-type star of significantly lower flux in the optical than the pre-SN object. A future detection of such star may be possible and would provide the first robust identification of a progenitor system for a Type Ib SN.

We calculate the evolution of close binary systems (CBSs) formed by a neutron star (behaving as a radio pulsar) and a normal donor star, evolving either to helium white dwarf (HeWD) or ultra short orbital period systems. We consider X-ray irradiation feedback and evaporation due to radio pulsar irradiation. We show that irradiation feedback leads to cyclic mass transfer episodes, allowing CBSs to be observed in-between as binary radio pulsars under conditions in which standard, non-irradiated models predict the occurrence of a low mass X-ray binary. This behavior accounts for the existence of a family of eclipsing binary systems known as redbacks. We predict that redback companions should almost fill their Roche lobe, as observed in PSR J1723-2837. This state is also possible for systems evolving with larger orbital periods. Therefore, binary radio pulsars with companion star masses usually interpreted as larger than expected to produce HeWDs may also result in such {\it quasi - Roche Lobe Overflow} states, rather than hosting a carbon-oxygen WD. We found that CBSs with initial orbital periods $\mathrm{P_{i}<1}$ day evolve into redbacks. Some of them produce low mass HeWDs, and a subgroup with shorter $\mathrm{P_{i}}$ become black widows (BWs). Thus, BWs descent from redbacks, although not all redbacks evolve into BWs. There is mounting observational evidence favoring that BW pulsars are very massive ($\mathrm{\gtrsim 2\; M_{\odot}}$). As they should be redback descendants, redback pulsars should also be very massive, since most of the mass is transferred before this stage.

We perform 1D calculations of neutrino opacities inside a young "strange star" assumed to be the result of the conversion process of a normal neutron object. We evaluate the deleptonization and cooling timescales, which happen to be longer than the proto-NS analogues, and preliminary address the features of the emerging neutrino signal.

We shall show in this report that the theoretical evolution of a particular class of systems containing a neutron star, the so-called "black widow" binaries, suggest masses above the $\sim 2 M_{\odot}$, a fact that would bring serious concerns about the right description of the equation of state above the saturation density. Moreover, so far the actual determinations of masses for these systems consistently give values above $2 M_{\odot}$, reinforcing the quandary. In this sense, and given that the confirmation of these ideas would create problems for the microphysical description, we argue that a crisis could be "in the works" in neutron star physics.

We present a new set of radial-velocity measurements of the spectroscopic binary HD 165052 obtained by disentangling of high-resolution optical spectra. The longitude of the periastron (60 +- 2 degrees) shows a variation with respect to previous studies. We have determined the apsidal motion rate of the system (12.1 +- 0.3 degree/yr), which was used to calculate the absolute masses of the binary components: M_1 = 22.5 +- 1.0 and M_2 = 20.5 +- 0.9 solar masses. Analysing the separated spectra we have re-classified the components as O7Vz and O7.5Vz stars.

The most recent member of the millisecond pulsar with very low-mass companions and short orbital periods class, PSR J1311-3430 (Pletsch et al. 2012) is a remarkable object in various senses. Besides being the first discovered in gamma-rays, its measured features include the very low or absent hydrogen content. We show in this Letter that this important piece of information leads to a very restricted range of initial periods for a given donor mass. For that purpose, we calculate in detail the evolution of the binary system self-consistently, including mass transfer and evaporation, finding the features of the new evolutionary path leading to the observed configuration. It is also important to remark that the detailed evolutionary history of the system naturally leads to a high final pulsar mass, as it seems to be demanded by observations.

We propose an alternative explanation for the post-breakout emission of SN 2008D associated with the X-ray transient 080109. Observations of this object show a very small contrast of 0.35 dex between the light-curve minimum occurring soon after the breakout, and the main luminosity peak that is due to radioactive heating of the ejecta. Hydrodynamical models show that the cooling of a shocked Wolf-Rayet star leads to a much greater difference (> 0.9 dex). Our proposed scenario is that of a jet produced during the explosion which deposits 56Ni-rich material in the outer layers of the ejecta. The presence of high-velocity radioactive material allows us to reproduce the complete luminosity evolution of the object. Without outer 56Ni it could be possible to reproduce the early emission purely from cooling of the shocked envelope by assuming a larger progenitor than a Wolf-Rayet star, but that would require an initial density structure significantly different from what is predicted by stellar evolution models. Analytic models of the cooling phase have been proposed reproduce the early emission of SN 2008D with an extended progenitor. However, we found that the models are valid only until 1.5 days after the explosion where only two data of SN 2008D are available. We also discuss the possibility of the interaction of the ejecta with a binary companion, based on published analytic expressions. However, the binary separation required to fit the early emission should be < 3 Rsun which is too small for a system containing two massive stars.

We perform binary stellar evolutionary calculations following the simultaneous evolution of both stars in the system to study a potential progenitor system for the Type IIb supernova 2011dh. Pre-explosion photometry as well as light-curve modeling have provided constraints on the physical properties of the progenitor system. Here we present a close binary system that is compatible with such constraints. The system is formed by stars of solar composition with 16 Msun + 10 Msun on a circular orbit with an initial period of 125 days. The primary star ends its evolution as a yellow supergiant with a mass of ~4 Msun, a final hydrogen content of ~3-5E-03 Msun and with an effective temperature and luminosity in agreement with the HST pre-explosion observations of SN 2011dh. These results are nearly insensitive to the adopted accretion efficiency factor beta. At the time of explosion, the companion star has an effective temperature of 22 to 40 thousand Kelvin, depending on the value of beta, and lies near the zero age main sequence. Considering the uncertainties in the HST pre-SN photometry the secondary star is only marginally detectable in the bluest observed band. Close binary systems, as opposed to single stars, provide a natural frame to explain the properties of SN 2011dh.

A set of hydrodynamical models based on stellar evolutionary progenitors is used to study the nature of SN 2011dh. Our modeling suggests that a large progenitor star ---with R ~200 Rsun---, is needed to reproduce the early light curve of SN 2011dh. This is consistent with the suggestion that the yellow super-giant star detected at the location of the SN in deep pre-explosion images is the progenitor star. From the main peak of the bolometric light curve and expansion velocities we constrain the mass of the ejecta to be ~2 Msun, the explosion energy to be E= 6-10 x 10^50 erg, and the 56Ni mass to be approximately 0.06 Msun. The progenitor star was composed of a helium core of 3 to 4 Msun and a thin hydrogen-rich envelope of ~0.1 M_sun with a main sequence mass estimated to be in the range of 12--15 Msun. Our models rule out progenitors with helium-core masses larger than 8 Msun, which correspond to M_ZAMS > 25 Msun. This suggests that a single star evolutionary scenario for SN 2011dh is unlikely.

The existence of millisecond pulsars with planet-mass companions (Bailes et al., 2011) in close orbits is challenging from the stellar evolution point of view. We calculate in detail the evolution of binary systems self-consistently, including mass transfer, evaporation and irradiation of the donor by X-rays feedback, demonstrating the existence of a new evolutionary path leading to short periods and compact donors as required by the observations of PSR J1719-1438. We also point out the alternative of an exotic nature of the companion planet-mass star.

In close binary systems composed of a normal, donor star and an accreting neutron star, the amount of material received by the accreting component is, so far, a real intrigue. In the literature there are available models that link the accretion disk surrounding the neutron star with the amount of material it receives, but there is no model linking the amount of matter lost by the donor star to that falling onto the neutron star. In this paper we explore the evolutionary response of these close binary systems when we vary the amount of material accreted by the neutron star. We consider a parameter \beta, which represents the fraction of material lost by the normal star that can be accreted by the neutron star. \beta is considered as constant throughout evolution. We have computed the evolution of a set of models considering initial donor star masses (in solar units) between 0.5 and 3.50, initial orbital periods (in days) between 0.175 and 12, initial masses of neutron stars (in solar units) of 0.80, 1.00, 1.20 and 1.40 and several values of beta. We assumed solar abundances. These systems evolve to ultracompact or to open binary systems, many of which form low mass helium white dwarfs. We present a grid of calculations and analyze how these results are affected upon changes in the value of \beta. We find a weak dependence of the final donor star mass with respect to \beta. In most cases this is also true for the final orbital period. The most sensitive quantity is the final mass of the accreting neutron star. As we do not know the initial mass and rotation rate of the neutron star of any system, we find that performing evolutionary studies is not helpful for determining \beta.

We construct a set of binary evolutionary sequences for systems composed by a normal, solar composition, donor star together with a neutron star. We consider a variety of masses for each star as well as for the initial orbital period corresponding to systems that evolve to ultra-compact or millisecond pulsar-helium white dwarf pairs. Specifically, we select a set of donor star masses of 0.50, 0.65, 0.80, 1.00, 1.25, 1.50, 1.75, 2.00, 2.25, 2.50, 3.00, and 3.50 solar masses, whereas for the accreting neutron star we consider initial masses values of 0.8, 1.0, 1.2, and 1.4 solar masses. The considered initial orbital period interval ranges from 0.5 to 12 days. It is found that the evolution of systems, with fixed initial values for the orbital period and the mass of the normal donor star, heavily depends upon the mass of the neutron star. In some cases, varying the initial value of the neutron star mass, we obtain evolved configurations ranging from ultra-compact to widely separated objects. We also analyse the dependence of the final orbital period with the mass of the white dwarf. In agreement with previous expectations, our calculations show that the final orbital period-white dwarf mass relation is fairly insensitive to the initial neutron star mass value. A new period-mass relation based on our own calculations is proposed, which is in good agreement with period-mass relations available in the literature. As consequence of considering a set of values for the initial neutron star mass, these models allow finding different plausible initial configurations (donor and neutron star masses and orbital period interval) for some of the best observed binary systems of the kind we are interested in here. We apply our calculations to analyse the case of PSR J0437-4715.

In the last few years, the so-called "Nice model" has got a significant importance in the study of the formation and evolution of the solar system. According to this model, the initial orbital configuration of the giant planets was much more compact than the one we observe today. We study the formation of the giant planets in connection with some parameters that describe the protoplanetary disk. The aim of this study is to establish the conditions that favor their simultaneous formation in line with the initial configuration proposed by the Nice model. We focus in the conditions that lead to the simultaneous formation of two massive cores, corresponding to Jupiter and Saturn, able to achieve the cross-over mass (where the mass of the envelope of the giant planet equals the mass of the core, and gaseous runway starts) while Uranus and Neptune have to be able to grow to their current masses. We compute the in situ planetary formation, employing the numerical code introduced in our previous work, for different density profiles of the protoplanetary disk. Planetesimal migration is taken into account and planetesimals are considered to follow a size distribution between $r_p^{min}$ (free parameter) and $r_p^{max}= 100$ km. The core's growth is computed according to the oligarchic growth regime. The simultaneous formation of the giant planets was successfully completed for several initial conditions of the disk. We find that for protoplanetary disks characterized by a power law ($\Sigma \propto r^{-p}$), smooth surface density profiles ($p \leq 1.5$) favor the simultaneous formation. However, for steep slopes ($p\sim 2$, as previously proposed by other authors) the simultaneous formation of the solar system giant planets is unlikely ...

It has been recently proposed that helium white dwarfs may provide promising conditions for the occurrence of the Bose-Einstein condensation. The argument supporting this expectation is that in some conditions attained in the core of these objects, the typical De Broglie wavelength associated with helium nuclei is of the order of the mean distance between neighboring nuclei. In these conditions the system should depart from classical behavior showing quantum effects. As helium nuclei are bosons, they are expected to condense. In order to explore the possibility of detecting the Bose-Einstein condensation in the evolution of helium white dwarfs we have computed a set of models for a variety of stellar masses and values of the condensation temperature. We do not perform a detailed treatment of the condensation process but mimic it by suppressing the nuclei contribution to the equation of state by applying an adequate function. As the cooling of white dwarfs depends on average properties of the whole stellar interior, this procedure should be suitable for exploring the departure of the cooling process from that predicted by the standard treatment. We find that the Bose-Einstein condensation has noticeable, but not dramatic effects on the cooling process only for the most massive white dwarfs compatible with a helium dominated interior (\approx 0.50 M_\odot) and very low luminosities (say, Log(L/L_\odot) < -4.0). These facts lead us to conclude that it seems extremely difficult to find observable signals of the Bose-Einstein condensation. Recently, it has been suggested that the population of helium white dwarfs detected in the globular cluster NGC 6397 is a good candidate for detecting signals of the Bose-Einstein condensation. We find that these stars have masses too low and are too bright to have an already condensed interior.

We present bolometric light curves of Type II-plateau supernovae (SNe II-P) obtained using a newly developed, one-dimensional Lagrangian hydrodynamic code with flux-limited radiation diffusion. Using our code we calculate the bolometric light curve and photospheric velocities of SN1999em obtaining a remarkably good agreement with observations despite the simplifications used in our calculation. The physical parameters used in our calculation are E=1.25 foe, M= 19 M_\odot, R= 800 R_\odot and M_{Ni}=0.056 M_\odot. We find that an extensive mixing of 56Ni is needed in order to reproduce a plateau as flat as that shown by the observations. We also study the possibility to fit the observations with lower values of the initial mass consistently with upper limits that have been inferred from pre-supernova imaging of SN1999em in connection with stellar evolution models. We cannot find a set of physical parameters that reproduce well the observations for models with pre-supernova mass of \leq 12 M_\odot, although models with 14 M_\odot cannot be fully discarded.

The core accretion mechanism is presently the most widely accepted cause of the formation of giant planets. For simplicity, most models presently assume that the growth of planetary embryos occurs in isolation. We explore how the simultaneous growth of two embryos at the present locations of Jupiter and Saturn affects the outcome of planetary formation. We model planet formation on the basis of the core accretion scenario and include several key physical ingredients. We consider a protoplanetary gas disk that exponentially decays with time. For planetesimals, we allow for a distribution of sizes from 100~m to 100~km with most of the mass in the smaller objects. We include planetesimal migration as well as different profiles for the surface density $\Sigma$ of the disk. The core growth is computed in the framework of the oligarchic growth regime and includes the viscous enhancement of the planetesimal capture cross-section. Planet migration is ignored. By comparing calculations assuming formation of embryos in isolation to calculations with simultaneous embryo growth, we find that the growth of one embryo generally significantly affects the other. This occurs in spite of the feeding zones of each planet never overlapping. The results may be classified as a function of the gas surface density profile $\Sigma$: if $\Sigma \propto r^{-3/2}$ and the protoplanetary disk is rather massive, Jupiter's formation inhibits the growth of Saturn. If $\Sigma \propto r^{-1}$ isolated and simultaneous formation lead to very similar outcomes; in the the case of $\Sigma \propto r^{-1/2}$ Saturn grows faster and induces a density wave that later acclerates the formation of Jupiter. Our results indicate that the simultaneous growth of several embryos impacts the final outcome and should be taken into account by planet formation models.

Giant planet formation process is still not completely understood. The current most accepted paradigm, the core instability model, explains several observed properties of the solar system's giant planets but, to date, has faced difficulties to account for a formation time shorter than the observational estimates of protoplanetary disks' lifetimes, especially for the cases of Uranus and Neptune. In the context of this model, and considering a recently proposed primordial solar system orbital structure, we performed numerical calculations of giant planet formation. Our results show that if accreted planetesimals follow a size distribution in which most of the mass lies in 30-100 meter sized bodies, Jupiter, Saturn, Uranus and Neptune may have formed according to the nucleated instability scenario. The formation of each planet occurs within the time constraints and they end up with core masses in good agreement with present estimations.

The equation of state calculated by Saumon and collaborators has been adopted in most core-accretion simulations of giant-planet formation performed to date. Since some minor errors have been found in their original paper, we present revised simulations of giant-planet formation that considers a corrected equation of state. We employ the same code as Fortier and collaborators in repeating our previous simulations of the formation of Jupiter. Although the general conclusions of Fortier and collaborators remain valid, we obtain significantly lower core masses and shorter formation times in all cases considered. The minor errors in the previously published equation of state have been shown to affect directly the adiabatic gradient and the specific heat, causing an overestimation of both the core masses and formation times.

Aims. In the context of the core instability model, we present calculations of in situ giant planet formation. The oligarchic growth regime of solid protoplanets is the model adopted for the growth of the core. Methods. The full differential equations of giant planet formation were numerically solved with an adaptation of a Henyey-type code. The planetesimals accretion rate was coupled in a self-consistent way to the envelope's evolution. Results. We performed several simulations for the formation of a Jupiter-like object by assuming various surface densities for the protoplanetary disc and two different sizes for the accreted planetesimals. We find that the atmospheric gas drag gives rise to a major enhancement on the effective capture radius of the protoplanet, thus leading to an average timescale reduction of 30% -- 55% and ultimately to an increase by a factor of 2 of the final mass of solids accreted as compared to the situation in which drag effects are neglected. With regard to the size of accreted planetesimals, we find that for a swarm of planetesimals having a radius of 10 km, the formation time is a factor 2 to 3 shorter than that of planetesimals of 100 km, the factor depending on the surface density of the nebula. Moreover, planetesimal size does not seem to have a significant impact on the final mass of the core.

Splaver and coworkers have measured the masses of the white dwarf and the neutron star components of the PSR J1713+0747 binary system pair by Shapiro Delay. We attempt to find the original configuration of this system performing a set of binary evolution calculations to simultaneously account for the masses of both stars and the orbital period. We considered initial masses of 1.5 and 1.4 \msun for the normal (donor) and the neutron star, respectively. We assumed two metallicity values (Z = 0.010 and 0.020), and an initial orbital period near 3 days. We assume that the neutron star is only able to retain \lesssim 0.10 of the matter transferred by the donor star. Calculations were performed employing our binary hydro code that handles the mass transfer rate in a fully implicit way together with state-of-the-art physical ingredients, diffusion and a non-grey atmospheres. We compare the structure of the resulting white dwarfs with the characteristic age of PSR J1713+0747 finding a nice agreement with observations by Lundgren et al. especially for the case of a donor star with Z= 0.010. This result indicates that the evolution of this kind of binary system is well understood. The models predict that, due to diffusion, the atmosphere of the white dwarf is an almost hydrogen-pure one. We find that such structures are unable to account for the colours measured by Lundgren et al. within their error bars. Thus, some discrepances in the white dwarf emergent radiation remain to be explained.

We study the evolution of a system composed by a 1.4 \msun neutron star and a normal, solar composition star of 2 \msun in orbit with a period of 1 day. Calculations were performed employing the binary hydro code presented in Benvenuto & De Vito (2003) that handle the mass transfer rate in a fully implicit way. Now we included the main standard physical ingredients together with diffusion processes and a proper outer boundary condition. We have assumed fully non - conservative mass transfer episodes. In order to study the interplay of mass loss episodes and diffusion we considered evolutionary sequences with and without diffusion in which all Roche lobe overflows (RLOFs) produce mass transfer. Another two sequences in which thermonuclearly-driven RLOFs are not allowed to drive mass transfer have been computed with and without diffusion. To our notice, this study represents the first binary evolution calculations in which diffusion is considered. The system produces a helium white dwarf of \sim 0.21 \msun in an orbit with a period of \sim 4.3 days for the four cases. We find that mass transfer episodes induced by hydrogen thermonuclear flashes drive a tiny amount of mass transfer. As diffusion produces stronger flashes, the amount of hydrogen - rich matter transferred is slightly higher than in models without diffusion. We find that diffusion is the main agent in determining the evolutionary timescale of low mass white dwarfs even in presence of mass transfer episodes.

We present a numerical code intended for calculating stellar evolution in close binary systems. In doing so, we consider that mass transfer episodes occur when the stellar size overflows the corresponding Roche lobe. In such situation we equate the radius of the star with the equivalent radius of the Roche lobe. This equation is handled implicitly together with those corresponding to the whole structure of the star. We describe in detail the necessary modifications to the standard Henyey technique for treating the mass loss rate implicitly together with thin outer layers integrations. We have applied this code to the calculation of the formation of low mass, helium white dwarfs in low mass close binary systems. We found that the global numerical convergence properties are fairly good. In particular, the onset and end of mass transfer episodes is computed automatically.

The present work is designed to explore the evolution of helium-core white dwarf (HeWD) stars for the case of metallicities much lower than the solar one (Z=0.001 and Z=0.0002). Evolution is followed in a self-consistent way with the predictions of detalied and new non-grey atmospheres, time-dependent element diffusion and the history of the white dwarf progenitor. Reliable initial models for low mass HeWDs are obtained by applying mass loss rates to a 1msun stellar model. The loss of angular momentum caused by gravitational wave emission and magnetic stellar wind braking are considered. Model atmospheres, based on a detailed treatment of the microphysics entering the WD atmosphere enable us to provide accurate colours and magnitudes at both early and advanced evolutionary stages. We find that most of our evolutionary sequences experience several episodes of hydrogen thermonuclear flashes. In particular, the lower the metallicity, the larger the minimum stellar mass for the occurrence fo flashes induced by CNO cycle reactions. The existence of a mass-threshold for the occurrence of diffusion-induced CNO flashes leadss to a marked dichotomy in the age of our models. Another finding of this study is that our HeWD models experience unstable hydrogen burning via PP nuclear reactions at late cooling stages as a result of hydrogen chemically diffusing inwards. Such PP flashes take place in models with very low metal content. We also find that models experiencing CNO flashes exhibit a pronouncede turn-off in most of their colours at M_V=16 approximately. Finally, colour-magnitude diagrams for our models are presented and compared with recent observational data of HeWD candidates in the globular clusters NGC 6397 and 47 Tucanae.

We examine the emergent fluxes from helium-core white dwarfs following their evolution from the end of pre-white dwarf stages down to advanced cooling stages. For this purpose, we include a detailed treatment of the physics of the atmosphere, particularly an improved representation of the state of the gas by taking into account non-ideal effects according to the so-called occupation probability formalism. The present calculations also incorporate hydrogen line opacity from Lyman, Balmer and Paschen series, pseudo-continuum absorptions and new updated induced-dipole absorption from H$_2$-H$_2$, H$_2$-He and H-He pairs. We find that the non-ideal effects and line absorption alter the appearance of the stellar spectrum and have a significant influence upon the photometric colours in the UBVRI-JHKL system. This occurs specially for hot models $T_{\rm eff}\ga 8000$ due to line and pseudo-continuum opacities, and for cool models $T_{\rm eff}\la 4000$ where the perturbation of atoms and molecules by neighbour particles affects the chemical equilibrium of the gas. In the present study, we also include new cooling sequences for helium-core white dwarfs of very low mass (0.160 and 0.148 M$_\odot$) with metallicity $Z=0.02$. These computations provide theoretical support to search for and identify white dwarfs of very low mass, specially useful for recent and future observational studies of globular cluster, where these objects have began to be detected.

The present work is designed to explore the effects of the time-dependent element diffusion on the mode trapping properties of DA white dwarf models with various thickness of the hydrogen envelope. Our predictions are compared with the standard assumption of diffusive equilibrium in the trace element approximation. We find that element diffusion markedly weakens the presence of mode trapping originated in the outer layers of the models, even for the case of thin hydrogen envelopes.

An adiabatic, non-radial pulsation study of a 0.563 M_sun DA white dwarf model is presented on the basis of new evolutionary calculations performed in a self-consistent way with the predictions of time dependent element diffusion, nuclear burning and the history of the white dwarf progenitor. Emphasis is placed on the role played by the internal chemical stratification of these new models in the behaviour of the eigenmodes, and the expectations for the full g-spectrum of periods. The implications for the mode trapping properties are discussed at length. In this regard, we find that, for high periods, the viability of mode trapping as a mode selection mechanism is markedly weaker for our models, as compared with the situation in which the hydrogen-helium transition region is treated assuming equilibrium diffusion in the trace element approximation.

We study the structural characteristic of the variable DA white dwarf G117B-15A by applying the methods of asteroseismology. For such a purpose, we construct white dwarf evolutionary models considering a detailed and up-to-date physical description as well as several processes responsible for the occurrence of element diffusion. We have considered several thickness for the outermost hydrogen layer, whereas for the inner helium-, carbon- and oxygen-rich layers we considered realistic profiles predicted by calculations of the white dwarf progenitor evolution. The evolution of each of the considered model sequences were followed down to very low effective temperature; in particular, from 12500K on we computed the dipolar, linear, adiabatic oscillations with low radial order. We find that asteroseismological results are not univocal regarding mode identification for the case of G117B-15A. However, our asteroseismological results are compatible with spectroscopical data only if the observed periods of 215.2, 271.0 and 304.4 s are due to dipolar modes with $k=2, 3, 4$ respectively. Our calculations indicate that the best fit to the observed period pattern of G117B-15A corresponds to a DA white dwarf structure with a stellar mass of 0.525 \msun, with a hydrogen mass fraction \lmh$\gtrsim$-3.83 at an effective temperature \teff$\approx$11800K.

We explore the effects of element diffusion due to gravitational settling and thermal and chemical diffusion on the pulsational properties of DA white dwarfs. To this end, we employ an updated evolutionary code coupled with a pulsational, finite difference code for computing the linear, non-radial g-modes in the adiabatic approximation. We follow the evolution of a 0.55 \msun white dwarf model in a self-consistent way with the evolution of chemical abundance distribution as given by time dependent diffusion processes. Results are compared with the standard treatment of diffusive equilibrium in the trace element approximation. Appreciable differences are found between the two employed treatments. We conclude that time dependent element diffusion plays an important role in determining the whole oscillation pattern and the temporal derivative of the periods in DAV white dwarfs. In addition, we discuss the plausibility of the standard description employed in accounting for diffusion in most of white dwarf asteroseismological studies.

We present a method to calculate masses for components of both eclipsing and non-eclipsing binary systems as long as their apsidal motion rates are available. The method is based on the fact that the equation that gives the rate of apsidal motion is a supplementary equation that allows the computation of the masses of the components, if the radii and the internal structure constants of them can be obtained from theoretical models. For this reason the use of this equation makes the method presented here model dependent. We apply this method to calculate the mass of the components of the non-eclipsing massive binary system HD 93205 (O3V+O8V), which is suspected to be a very young system. To this end, we computed a grid of evolutionary models covering the mass range of interest, and taking the mass of the primary (M_1) as the only independent variable, we solve the equation of apsidal motion for M_1 as a function of the age of the system. The mass of the primary we find ranges from M_1= 60+-19 msun for ZAMS models, which sets an upper limit for M_1, down to M_1= 40+-9 msun for an age of 2 Myr. Accordingly, the upper limit derived for the mass of the secondary (M_2= Q M_1) is M_2= 25 msun is in very good agreement with the masses derived for other O8V stars occurring in eclipsing binaries.

The purpose of this paper is to present new full evolutionary calculations for DA white dwarf stars with the major aim of providing a physically sound reference frame for exploring the pulsation properties of the resulting models in future communications. Here, white dwarf evolution is followed in a self-consistent way with the predictions of time dependent element diffusion and nuclear burning. In addition, full account is taken of the evolutionary stages prior to the white dwarf formation. In particular, we follow the evolution of a 3 \msun model from the zero-age main sequence (the adopted metallicity is Z=0.02) all the way from the stages of hydrogen and helium burning in the core up to the thermally pulsing phase. After experiencing 11 thermal pulses, the model is forced to evolve towards its white dwarf configuration by invoking strong mass loss episodes. Further evolution is followed down to the domain of the ZZ Ceti stars on the white dwarf cooling branch. Emphasis is placed on the evolution of the chemical abundance distribution due to diffusion processes and the role played by hydrogen burning during the white dwarf evolution. Furthermore, the implications of our evolutionary models for the main quantities relevant for adiabatic pulsation analysis are discussed. Interestingly, the shape of the Ledoux term is markedly smoother as compared with previous detailed studies of white dwarfs. This is translated into a different behaviour of the Brunt-Vaisala frequency.

In this letter we investigate the pulsational properties of ZZ Ceti stars on the basis of new white dwarf evolutionary models calculated in a self-consistent way with the predictions of time dependent element diffusion and nuclear burning. In addition, full account is taken of the evolutionary stages prior to the white dwarf formation. Emphasis is placed on the trapping properties of such models. By means of adiabatic, non-radial pulsation calculations, we find, as a result of time dependent diffusion, a much weaker mode trapping effect, particularly for the high-period regime of the pulsation g-spectrum. This result is valid at least for models with massive hydrogen-rich envelopes. Thus, mode trapping would not be an effective mechanism to explain the fact that all the high periods expected from standard models of stratified white dwarfs are not observed in the ZZ Ceti stars.

We present a finite difference code intended for computing linear, adiabatic, nonradial pulsations of spherical stars. This code is based on a general Newton - Raphson technique in order to handle the relaxation of the eigenvalue (square of the eigenfrequency) of the modes and their corresponding eigenfunctions. This code has been tested computing the pulsation spectra of polytropic spheres finding a good agreement with previous work. Then, we have coupled this code to our evolutionary code and applied it to the computation of the pulsation spectrum of a low mass, pure - helium white dwarf of 0.3 M_{sun} for a wide range of effective temperatures. In making this calculation we have taken an evolutionary time step short enough such that eigenmodes corresponding to a given model are used as initial approximation to those of the next one. Specifically, we have computed periods, period spacing, eigenfunctions, weight functions, kinetic energies and variational periods for a wide range of modes. To our notice this is the first effort in studying the pulsation properties of helium white dwarfs. The solution we have found working with these realistic white dwarf models are in good accord with the predictions of the asymptotic theory of Tassoul (1980) for high order modes. This indicates that the code presented here is able to work adequately also with realistic stellar models.

White dwarfs are well studied objects. The relative simplicity of their physics allows to obtain very detailed models which can be ultimately compared with their observed properties. Among white dwarfs there is a specific class of stars, known as ZZ-Ceti objects, which have a hydrogen-rich envelope and show periodic variations in their light curves. G117-B15A belongs to this particular set of stars. The luminosity variations have been successfully explained as due to g-mode pulsations. G117-B15A has been recently claimed to be the most stable optical clock ever found, being the rate of change of its 215.2 s period very small: \dot{P}= (2.3 +- 1.4)x10^{-15} s s^-1, with a stability comparable to that of the most stable millisecond pulsars. The rate of change of the period is closely related to its cooling timescale, which can be accurately computed. In this paper we study the pulsational properties of G117-B15A and we use the observed rate of change of the period to impose constraints on the axion emissivity and, thus, to obtain a preliminary upper bound to the mass of the axion. This upper bound turns out to be 4cos^{2}{\beta} meV at the 95% confidence level. Although there are still several observational and theoretical uncertainties, we conclude that G117-B15A is a very promising stellar object to set up constraints on particle physics.

The purpose of this work is to explore the evolution of helium-core white dwarf stars in a self-consistent way with the predictions of detailed non-gray model atmospheres and element diffusion. To this end, we consider helium-core white dwarf models with stellar masses of 0.406, 0.360, 0.327, 0.292, 0.242, 0.196 and 0.169 solar masses and follow their evolution from the end of mass loss episodes during their pre-white dwarf evolution down to very low surface luminosities. We find that when the effective temperature decreases below 4000K, the emergent spectrum of these stars becomes bluer within time-scales of astrophysical interest. In particular, we analyse the evolution of our models in the colour-colour and colour-magnitude diagrams and we find that helium-core white dwarfs with masses ranging from approx. 0.18 to 0.3 solar masses can reach the turn-off in their colours and become blue again within cooling times much less than 15 Gyr and then remain brighter than M_V approx. 16.5. In view of these results, many low-mass helium white dwarfs could have had time enough to evolve to the domain of collision-induced absorption from molecular hydrogen, showing blue colours.

We investigate the effects of element diffusion on the structure and evolution of low-mass helium white dwarfs (WD). Attention is focused on the occurrence of hydrogen shell flashes induced by diffusion processes during cooling phases. Initial models from 0.406 to 0.161 solar masses are constructed by applying mass loss rates at different stages of the RGB evolution of a solar model. The multicomponent flow equations describing gravitational settling, and chemical and thermal diffusion are solved and the diffusion calculations are coupled to an evolutionary code. In addition, the same sequences are computed but neglecting diffusion. We find that element diffusion strongly affects the structure and cooling history of helium WD. In particular, diffusion induces the occurrence of hydrogen shell flashes in models with masses ranging from 0.18 to 0.41 solar masses, which is in sharp contrast from the situation when diffusion is neglected. In connection with the further evolution, these diffusion-induced flashes lead to much thinner hydrogen envelopes, preventing stable nuclear burning from being an appreciable energy source at advanced stages of evolution. This implies much shorter cooling ages than in the case when diffusion is neglected. These new WD models are discussed in light of recent observational data of some millisecond pulsar systems with WD companions. We find that age discrepancies between the predictions of standard evolutionary models and such observations appear to be the result of ignoring element diffusion in such models. Indeed, such discrepancies vanish when account is made of diffusion.

A 0.242 Msun object that finally becomes a helium white dwarf is evolved from Roche lobe detachment down to very low luminosities. In doing so, we employ our stellar code to which we have added a set of routines that compute the effects due to gravitational settling, and chemical and thermal diffusion. Initial model is constructed by abstracting mass to a 1 Msun red giant branch model up to the moment at which the model begins to evolve bluewards. We find that element diffusion introduces noticeable changes in the internal structure of the star. In particular, models undergo three thermonuclear flashes instead of one flash as we found with the standard treatment. This fact has a large impact on the total mass fraction of hydrogen left in the star at entering the final cooling track. As a result, at late stages of evolution models with diffusion are characterized by a much smaller nuclear energy release, and they evolve significantly faster compared to those found with the standard treatment. We find that models in which diffusion is considered predict evolutionary ages for the white dwarf companion to the millisecond pulsar PSR B1855+09 in good agreement with the spin-down age of the pulsar.

Recent measurements made by Hipparcos (Provencal et al. 1998) present observational evidence supporting the existence of some white dwarf (WD) stars with iron - rich, core composition. In this connection, the present paper is aimed at exploring the structure and evolution of iron - core WDs by means of a detailed and updated evolutionary code. In particular, we examine the evolution of the central conditions, neutrino luminosity, surface gravity, crystallization, internal luminosity profiles and ages. We find that the evolution of iron - rich WDs is markedly different from that of their carbon - oxygen counterparts. In particular, cooling is strongly accelerated as compared with the standard case. Thus, if iron WDs were very numerous, some of them would have had time enough to evolve at lower luminosities than that corresponding to the fall - off in the observed WD luminosity function.

The purpose of this work is to present accurate and detailed mass-radius relations for white dwarf (WD) models with helium, carbon, oxygen, silicon and iron cores, by using a fully updated stellar evolutionary code. We considered masses from 0.15 to 0.5 Mo for the case of helium core, from 0.45 to 1.2 Mo for carbon, oxygen and silicon cores and from 0.45 to 1.0 Mo for the case of an iron core. In view of recent measurements made by Hipparcos that strongly suggest the existence of WDs with an iron-dominated core, we focus our attention mainly on the finite-temperature, mass-radius relations for WD models with iron interiors. Furthermore, we explore the effects of gravitational, chemical and thermal diffusion on low-mass helium WD models with hydrogen and helium envelopes.

There are at least three ways in which a varying gravitational constant $G$ could affect the interpretation of the recent high-redhisft Type Ia supernovae results. If the local value of $G$ at the space-time location of distant supernovae is different, it would change both the thermonuclear energy release and the time scale of the supernova outburst. In both cases the effect is related to a change in the Chandrasekhar mass $M_{\rm Ch}\propto G^{-3/2}$. Moreover the integrated variation of $G$ with time would also affect cosmic evolution and therefore the luminosity distance relation. Here we investigate in a consistent way how these different effects of a varying $G$ could change the current interpretation of the Hubble diagram of Type Ia supernovae. We parametrize the variation of $G$ using scalar-tensor theories of gravity, such as the Jordan-Brans-Dicke theory or its extensions. It is remarkable that Dirac's hypothesis that $G$ should decrease with time can qualitatively explain the observed $\Delta m \simeq 0.2$ mag decrease at $z\simeq 0.5$ (with respect to a decelerating universe) and, at the same time, reduce the duration of the risetimes of distant Type Ia supernovae as recently reported.

We present detailed evolutionary calculations for carbon - oxygen - and helium - core white dwarf (WD) models with masses ranging from M= 0.1 to M= 1.2 solar masses and for metallicities Z= 0.001 and Z= 0. The sequences cover a wide range of hydrogen envelopes as well. We employed a detailed WD evolutionary code. In particular, the energy transport by convectcion is treated within the formalism of the full spectrum turbulence theory. The set of models presented here is very detailed and should be valuable for the interpretation of the observational data on low - mass WDs recently discovered in numerous binary configurations and also for the general problem of determining the theoretical luminosity function for WDs. In this context, we compare our cooling sequences with the observed WD luminosity function recently improved by Leggett, Ruiz and Bergeron (1998) and we obtain an age for the Galactic disc of approximately 8 Gyr. Finally, we applied the results of this paper to derive stellar masses of a sample of low - mass white dwarfs.

In this study we compute the structure and evolution of carbon-oxygen DA white dwarfs by means of a detailed and updated evolutionary code. We treat the energy transport by convection within the formalism of the full spectrum turbulence theory, as given by the Canuto, Goldman and Mazzitelli (CGM) model. We explore the effect of various hydrogen layer masses on both the surface gravity and the hydrogen burning. Convective mixing at low luminosities is also considered. One of our main interests in this work has been to study the evolution of ZZ Ceti models, with the aim of comparing the CGM and mixing length theory (MLT) predictions. In this connection, we find that the temperature profile given by the CGM model is markedly different from that of the ML1 and ML2 versions of MLT. We have also computed approximate effective temperatures for the theoretical blue edge of the DA instability strip by using thermal timescale arguments for our evolving models. In this context, we found that the CGM theory leads to blue edges that are cooler than the observed ones. However, because the determination of the atmospheric parameters of ZZ Ceti stars is dependent on the assumed convection description in model atmosphere calculation, observed blue edges computed considering the CGM theory are required in order to perform a sef consistent comparison of our results with observations. Finally, detailed pulsation calculations of ZZ Ceti models considering the CGM convection would be necessary to place the results found in this paper on a firmer basis.