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.

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.

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 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.

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.

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.

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 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.

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.

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.