A Code for Stellar Binary Evolution and its Application to the Formation of Helium White Dwarfs

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.