Mass Transfer in Binary Stellar Evolution and Its Stability

The evolution of a binary star system by various analytical and numerical approximations of mass transfer rate normalized to the equilibrium rate and its stability conditions are investigated. We present results from investigations of mass transfer and stability in close binary star systems using the different orbital parameters. The stability and instability of mass transfer in binary star evolution depends on the exchange of material which the response of the binary to the initial Roche lobe overflow causes the donor to loose even more material. Our work is mainly focused on basic mathematical derivations, analytical and numerical solutions in order to explain the mass transfer system in different orbital parameters as well as the results are compared with previous studies in both cases. Mass transfer is usually stable, as long as the winds specific angular momentum does not exceed the angular momentum per reduced mass of the system. This holds for both dynamical and thermal time scales. Those systems which are not stable will usually transfer mass on the thermal time scale. The variation of Roche lobe radius with mass ratio in the binary, for various orbital parameters in the conservative and non-conservative mass transfer, as well as the evolution equations, orbital angular momentum of the binary system and the corresponding analytical and numerical solutions for different cases, under certain restrictive approximations is derived, simulated and discussed.


Introduction
The dominant feature in the evolution of stars is in closed binaries, and the one which distinguishes it from the evolution of single stars is the presence of various forms of mass transfer between the two stars.Mass transfer occurs in many different types of systems, to widely varying effects [1].Contamination of the envelope of a less evolved star with chemically processed elements, as in Barium stars [2], and catastrophic mass transfer, by common envelope phase, as in W UMa systems; or a slow, steady mass transfer by Roche lobe overflow [2].
Close binary stars consisting of two compact stellar remnants (white dwarfs (WDs), neutron stars (NSs), or black holes (BHs)) are considered as primary targets of the forthcoming field of gravitational wave (GW) astronomy [3] since their orbital evolution has entirely controlled by the emission of gravitational waves and leads to ultimate coalescence (merger) of the components.GW emission is the sole factor responsible for the change of orbital parameters of a pair of compact (degenerate) stars.However, at the early stages of binary evolution, it is the mass transfer between the components and the loss of matter and its orbital angular momentum that play a dominant dynamical role.
Mass transfer is particularly interesting if one considers the evolution of a system with at least one degenerate star.In these cases, mass transfer produces spectacular effects, resulting in part from the intense magnetic and gravitational fields of the compact objects, pulsed X-ray emission, nuclear burning, novae outbursts, and so on.Also, since the mass transfer rates can be high and orbital period measurements accurate, one may see the dynamical effects of mass transfer on the binary orbit, as in Cygnus X-3 [4].
In the case of cataclysmic variables (CVs) and low mass X-ray binary system (LMXBs), one has a highly evolved, compact star (CVs and LMXBs contain white dwarfs and neutron stars, respectively) and a less evolved main sequence or red giant star.Mass transfer usually proceeds by accretion onto the compact object, and is secularly stable.Most of these objects involve one component filling its Roche lobe (the donor or secondary), transferring mass to a compact star (the accretor or primary) and the mass leaking out of the inner Lagrangian point forms an accretion disk around the accretor.A binary system starts out as detached [4] [5], with both stars in hydrostatic equilibrium and filling an equipotential surface inside their Roche lobes.In a binary, both stars (1 and 2 can fill their Roche lobes at subsequent evolutionary stages).We use the convention which 1 denotes the primary star, and 2 denotes the secondary star of the binary.In this convention, the mass ratio is q m m m m = = > .During phases of mass transfer, we denote the mass-losing star, the "donor", with subscript "d" and the companion star, the "accretor", with subscript "a".The transfer will be accompanied by a stellar wind from the mass-losing star, or ejection of matter from the accretor, as in novae and galactic jet sources.
There are various unanswered questions in the evolution of LMXBs into low mass binary pulsars (LMBPs, in which a millisecond radio pulsar is in a binary with a low mass white dwarf companion) and of CVs.Among these are the problem of the disparate birthrates of the LMXBs and the LMBPs (Kulkarni and Narayan1988), estimation of the strength of X-ray heating induced winds from the donor star.Now, if a star is a member of a binary star system, then it is possible that as the components evolving, one of them (or both) can fill up its Roche lobe.Thomas (1977) classified mass transferring binaries on the basis of what state of evolution the donor is in as follows: Case A: If the orbital separation is small enough (usually a few days), the star can fill its Roche lobe during its slow expansion through the main-sequence phase while still burning hydrogen in its core.
Case B: If the orbital period is less than about 100 days, but longer than a few days, the star will fill its Roche lobe during the rapid expansion to a red giant with a helium core.If the helium core ignites during this phase and the transfer is interrupted, the mass transfer is Case B.
Case C: If the orbital period is above 100 days, the star can evolve to the red supergiant phase before it fills its Roche lobe.In this case, the star may have a CO or ONe core.
Case A mass transfer occurs during the slow growth, Case B during the first rapid expansion, and Case C during the final expansion phase.The nature of the remnant depends upon the state of the primary during the onset of RLOF and the orbital properties of the resultant binary depend upon the details of the mass transfer.
The problem of mass transfer in binaries by Roche lobe overflow has received a good deal of attention in the literature over the past few decades, typically in investigations of one aspect or another of orbital evolution or stability.[6] examined mass transfer by isotropic winds and accretion, in investigating the evolution of short period binary X-ray sources with extreme mass ratios, and developed models with accretion and (typically) isotropic re-emission of transferred matter, in the context of the period gap in cataclysmic variables.
Stability has become an important in many area of astrophysical research.Many authors has been studied the stability in different models (for example, [7] [8]).Both authers didn't address stability and mass transfer in binary stellar evolution.[1] [2] [4] were tried to study stability criteria for mass transfer in binary stellar evolution using limits of pure models of mass transfer and extreme mass ratio.The present study will address the mass transfer in binary stellar evolution using analytical and numerical methods which will give brief relations about stability and mass transfer.
The aim of this paper is to study mass transfer in binary stellar evolution behavior using mass transfer and different orbital parameters.This will be used to study stability properties of binary stars using analytical and numerical solutions.We calculate the rate of change of orbital angular momentum, orbital periods, angular speed and stability analysis of mass transfer rates, with relevant evolution in general form.
The paper is organized as follows: In Section (2), we derive basic equations of mass transfer, mass transfer and evolution of orbital parameters, conservative and non-conservative mass transfer and evolution equations; in Section (3), stability analysis of mass transfer and its rate will be determined; in Section (4), analytical and numerical solutions will be presented.Finally our main conclusions are summarized in Section (5).

Mass Transfer and the Evolution of Orbital Parameters
In this section, expressions for the variation of orbital parameters with loss of mass from one of the stars are derived.In what follows, the two stars will be referred to as 1 m and 2 m , with the latter the mass losing star.The angular momentum (AM) of binary component i in a circular orbit is expressed as: where is the velocity of accretor and donor star; "i" stands for the accretor "a" and donor "d" stars respectively, i a is the position of the accretor and donor star, and ω is the orbital angular spin frequency.
Roche geometry changes are governed by the total orbital angular momentum J of the system which is given by ( ) where are the distance of the two stars from the center of mass (semi major-axis),

(
) is the total mass, and the mass ratio is Throughout this paper, we consider an eccentricity e will be zero.The period "P", the semi major axis "a" and the total mass transfer related through Kepler's law as: ( ) which implies that the orbital angular momentum J can be written as: ( ) The angular momentum for binary component "i" is given by and the specific angular momentum for binary component "i" can be written as: where i m is mass of the accretor and donor and i l is specific angular momentum of the accretor and donor star, respectively.From Kepler's 3 rd we relate orbital evolution as: Hence from Equation (4) the orbital angular momentum J of a binary star is given by ( ) In many cases the angular momentum stored in the rotation of the two stars is negligible compared to the orbital angular momentum, so that Equation (8) also represents the orbital angular momentum of the binary, to good approximation.By differentiating this expression we obtain a general equation for the evolution of orbital parameters: and ( ) In Equation ( 10), e is a parameter that determines the amount by which its orbit around another body deviates from a perfect circle.From Equation (7) and Equation (10) the rate of change of orbital period, semi major axis, and total mass transfer can be written as: which depends on time "t".By integrating Equation (11) we obtain: where 0 0 0 , , and P a M are initial orbital period, semi major axis, and total mass transfer of the binary stars, respectively.
In the case of Roche lobe overflow in an already circularized binary, the last term is zero.The J  term represents angular momentum loss from the binary, which can be due to spontaneous processes (such as gravitational wave radiation) or it can be associated with mass loss from the binary as a whole or from the component stars.

Conservative Mass Transfer
When no ejected matter leaves a binary system, the mass transfer is said to be conservative.During conservative mass transfer, the orbital elements of the binary can change due to transfer of angular momentum from one star to the companion.Consider a system with a total mass , semi-major axis, a, eccentricity, e, and the total orbital angular momentum, J is given by ( ) will also be conserved, where G, is universal gravitational constant.We first consider conservative mass transfer, in which the total mass and orbital angular momentum of the binary are conserved.In that case we can set as: In circular, 0 e e = =  and conservative mass exchange, 0 M J = =   and 0 a d m m = − > , then from Equation (13) we have: where ( ) , is a constant.Hence from Equation (10) we obtain: ( ) Using the Kepler's third law in Equation ( 15), the period variation with time P  due to mass transfer can be written as: and integrating Equation ( 16) which leads us 3 0 0 a P P a An explicit relation between the separation and the masses can be found by integrating Equation (15) with 0 e = ; 2 2 d a m m a = constant.Applying Kepler's law we obtain similar relation between the rate of change of the orbital period P  and the rate of mass transfer d m  , and between the period and masses directly from Equation ( 15): ( ) ( ) The importance of Equation (15) and Equation ( 18) is that it allows to determine the mass transfer rate of observed semi-detached binaries, if the masses and the period derivative can be measured.This is complicated by the fact that many binaries show short-term period fluctuations, while for the long-term average of the period derivative, the trend has been determined with reasonable conditions as follows: , mass transfer from the less massive star to the more massive one will cause the orbit to enlarge.
If 1 q > , 0 a <  , mass transfer from the more massive star to the less massive one will cause the orbit to shrink.

Non-Conservative Mass Transfer
In case of non-conservative mass transfer both mass and angular momentum can be removed from the system.Following [9] Orbital Angular Momentum (OAM) of a two body system is given by ( ) where is moment of inertia, and 2π is the angular frequency of the binary.
Assuming isotropic mass loss from the surface of the components, then, we have the rate of angular momentum: ( ) Since the dynamics of a two body system obey Kepler's third law therefore one expects transfer of mass would change the, a, P and M accordingly Assume q as a constant, differentiating Equation ( 19) with respect to time gives: Substituting Equation (21) into Equation ( 22) we get the following two equations: and Therefore, the two equations tell us how the loss of OAM and mass are caused the orbital period, P, and orbital radius, a, to change.The isotropic mass loss implies: This equation implies that the only source of OAM loss is mass loss, which were assumed to be isotropic.Substituting J J  from equation Equation (25) into Equation (23) will give: In Equation (26), since the P P  is measurable observationally, therefore, M M  can be calculated from Equation (26).Now if we use the following equations taken from Stepien (1995) to calculate relative angular momentum lost from a system as: ( ) ( ) and hence the relative OAM lost only by magnetized star wind: lost 1 3 where This equation (Equation (30)) tells us the total rate of change of the orbital angular frequency of the binary related to the mass lost from the system.

Evolution Equations
Now we have to enumerate some of the physical mechanisms that can influence the evolution of a binary, we can proceed to derive mathematical expressions in order to quantify these effects.To reiterate, we work under the following assumptions: i) Roche potential describes the gravitational field; ii) Kepler's laws are valid, and that the stars orbit around their common center of mass in circular orbits; iii) the spin and orbit axes are all parallel to one another and; iv) tidal effects are included even though we ignore the effects of the distortion of the stars, i.e. we assume spherically symmetric stars.Now, the total angular momentum of the system is given by where i k is dimensionless constants depending on the internal structure of the the angular spin frequencies( i ω ), and ( ) are the spins of the accretor and donor respectively.The first term in Equation (31) represents the orbital angular momentum of the components, and the two other terms represent the spin angular momenta of the stars.

(
) The rate of change of the spin momenta can be given as the sum of a consequential angular momentum term and a tidal term as: where, a l and d l indicate the specific angular momenta of the matter arriving at the accretor and the matter leaving the donor respectively.
The second term in the Equation (33) represents the tidal term, which is a function of the degree of ( ) and the tidal time scale ( ) Substituting Equation (33) and Equation (34) into Equation (31).
where ( ) Accordingly, from Equation (35) the rate of change of the orbital separation is given by which we re-write for the convenience as: The sign of the terms on the right hand side of Equation (37) determines if the orbit is shrinking or expanding.a q represents the net effect of the consequential transfer of angular momentum.In the absence of mass transfer, the orbit shrinks as a result of the systemic angular momentum loss.
If a q q > then the binary continues to shrink.On the other hand, for systems with a q q < , mass transfer leads to the orbit expanding or least not shrinking as rapidly as it would in the absence of mass transfer.
The tidal terms have a more complicated behavior because the sign of the tidal torque , i tid J is a function of the difference in the spin frequencies and the orbital frequency, i ω Ω − .If the spin frequency of a given component is higher than the orbital frequency, tides will tend to pump angular momentum into the orbit and thus help increase the separation.Conversely, if the orbital frequency is higher than the spins, it can lead to a runaway instability since the tides will tend to suck more angular momentum out from the orbit which in turn will increase the orbital frequency even more.
We now determine how the depth of contact evolves to complete the set of equations we need to specify the evolution of the binary system.We can write equation of the Roche lobe size as: where L R , L R  are the Roche lobe radius and the rate of change of the Roche Lobe radius and Symbolically, generalizing the meaning of the symbols introduced by [5] ( ) where the symbol L υ stands for driving terms and L ζ denotes logarithmic derivatives with respect to donor mass.We write the logarithmic time derivative of the donor radius ( ) as: where d υ represents the rate of change of the donor radius due to intrinsic processes such as thermal relaxation and nuclear evolution, whereas d ζ usually describes changes resulting from adiabatic variations of d m as described in Equation (44).We derive here a simple analytic approximation to the effective mass-radius exponent when the response of the donor is a combination of the adiabatic and thermal adjustments to mass loss.from Equation (44).
where d ζ is the effective mass-radius exponent and d υ stands for the thermal radial reaction rate, and s ζ is here the purely adiabatic mass-radius exponent.As a consequence of mass loss, the donors radius d R will differ from the equilibrium radius corresponding to its instantaneous mass ( ) . With these definitions we write: ( ) The secular evolution of the binary takes place on the mass transfer time scale Equation ( 46) with respect to time, we get: in Equation (47), and solving for d ζ , we obtain: This expression shows that if the evolution is much slower than the thermal relaxation ( ), the donor radius follows the equilibrium radius closely, whereas if mass transfer occurs rapidly, the donor reacts adiabatically.Equations ((34), (36), (40), and (44)) referred to as "the evolution equations".

Mass Transfer Rate
Mass transfer will proceed on a time-scale which depends critically on the changes in the radius of the donor star and that of its Roche lobe in response to the mass loss; if the star expands faster than its Roche lobe or shrinks less rapidly than its Roche lobe for a prolonged time, mass transfer will be unstable and the donor star may disintegrate.If the donor star expands less rapidly or shrinks faster than its Roche lobe, mass transfer will generally be stable and may continue for a long time.
We now introduce the concept of the equilibrium mass transfer rate, which is the mass transfer rate that a stable, semi-detached binary undergoes for a given rate of driving.The equilibrium mass transfer rate is a function of the driving rate, the consequential angular momentum loss mechanisms and the value of the mass ratio q.The mass transfer rate can be written generally as: Under the assumption that 0 M  is a function of the binary parameters, whilst f is a strong function of the depth of contact d R ∆ , we can write: For equilibrium mass transfer, we let 0 d m =  , i.e. the mass transfer rate itself does not evolve very fast.Thus, using the evolution equations and Equation (46) we have: where, stable 2 2 is the critical mass ratio for stability of mass transfer.For example, for a polytropic donor with 3 2 n = (which is representative of a white dwarf donor) 1 3 and so 2 3 In general one can write from Equation (51).
when stable q q < , the pre-factor on the right hand side of Equation ( 54) is negative and thus the mass transfer rate tends toward the equilibrium value implying that the mass transfer is stable.On the other hand, when stable q q > , we can see that the system cannot reach the equilibrium value and this implies that the mass transfer is unstable.It is possible during the course of the evolution that a system that initially has stable q q > can evolve into a system with stable q q < .One of the components of the binary gets into contact with its Roche lobe either due to the expansion of the star (the donor) and/or due to shrinking of the Roche lobe as a result of angular momentum loss by GWR.In order to determine exactly when a binary with given total mass and mass ratio q becomes semi-detached, one needs to specify the radii of the component stars and their corresponding Roche lobes.The condition for stability can be expressed as: which simply states that the donors radius must shrink at least as fast as its Roche lobe or that the donor can expand no faster than its Roche lobe expands.
The first factor we will consider is the change in the orbital separation resulting from a mass transfer event.The orbital angular momentum, orb J for a binary composed of two point masses is given by Equation ( 2).We can use Kepler's third law: to eliminate the angular frequency from Equation ( 2) to obtain: If we logarithmically differentiate Equation (57) with respect to time and from Equation (10) with 0 e = , we obtain: From Equation (15), Equation (58) becomes ( ) We can then see that if the donor is the more massive star (q > 1), the orbital separation will shrink upon mass transfer.Conversely, if q < 1, a will increase.
From Roche lobe geometry and following [10], the Roche lobe radius d L R can be crudely approximated as: A more accurate approximation for the Roche lobe radius was determined by Eggleton (1983) and is given by ( ) for 0 q < < ∞ Eggleton (1983).Differentiating Equation (60) with respect to time, we obtain: Combining Equation (62) with Equation (58) we arrive at For conservative mass transfer Equation (64) becomes 5 23 which parametrized as: In the conservative case Equation (65) indicates that if 5 6 the donor Roche lobe will be contract upon mass loss as d m  is negative.Finally, we consider the response of the donors radius to change in its mass.For stars like the sun it is well known that the mass and radius are approximately proportional to one another.Using the expression of (66), the mass radius relation for solar type stars implies If we substitute Equation (65) and Equation (67) into Equation (55) we obtain a limiting stable mass ratio for the binary star.If the mass ratio exceed this values the Roche lobe will be shrink faster than the star can contract and mass transfer will proceed on a dynamical time scale.If 4 3 q ≤ the star will, on a time scale set by the mass transfer rate detach from is Roche lobe.From Equation (64) any mechanism Schutz (1990) that removes orbital angular momentum from the binary will cause the Roche lobe of the donor to contract.

Analytical Solutions
Before we proceed upon studying the complete numerical solutions for the evolution of orbital parameters and the evolution equations we derived in Section (2) for different astrophysical scenarios, we obtain analytic solutions for the time evolution of the geometry of mass transfer rate.We can do this only on the assumption that most of the parameters characterizing the binary remain constant or evolve slowly as compared to the evolution of the mass transfer rate.[5] was among the first to attempt such analysis and following them, we generalize their results to an arbitrary polytropic index and to isothermal atmospheres.Analytic solutions are useful in providing physical insights into the expected behaviour of the binary system in the limit where the assumptions imposed to obtain the analytic solutions are valid approximately.To make Equations ((60), ( 61) and ( 68)), the two figures are plotted for different values of the mass ratio The simple form of Equation ( 60) is a reasonable approximation to the more accurate Eggleton [11] formula for mass ratios 0.9 q  . Thus we see that the Roche lobe radius is directly proportional to the separation for a fixed mass ratio.Given Roche geometry and the structure of the star, we can now specify the mass transfer rate.In general, we expect the mass transfer to be a strong function of the "depth of contact", defined as the amount by which the donor overflows its Roche lobe and of the structure of the star, which can be described in different ways depending on the problem at hand.Thus if 0, 0 which implies that the orbit will be enlarge and if , 0 that the orbit will be shirk.This means the radius of a Roche lobe filling star depends only on the binary separation and the mass ratio.The blue curve shows that the stable conservative mass transfer while the green one is the unstable mass transfer.
In Figure 2, we have plotted the equilibrium angular momentum values for a range of binary system with a time and for a mass ratio q = 1, taking the mass of the black hole 1.75 We here produce analytical solutions for the non conservative mass transfer in close binary system taking the initial masses of the primary and secondary as (years) and period p = 100 (days), the orbital angular momentum J continues to show a steady decreasing with the increasing time in Figure 3 and Figure 4.
Considering an interacting double degenerate system with initial masses of donor and accretor 0.24  (dotted line) and a "semi-degenerate" low-mass helium star donor plus white dwarf accretor of the same initial masses (dash or blue line), and the periods (P = 2, 3, 4) from Roche lobe overflow; lower limits of the disk is equilibrium instability region according to [12] for 0.5 q = and conventional evolutionary computations may be not adequate for description of mass transfer process.
From Equation ( 18), the change of orbital period is proportional to the change in dot M .Prolonged, conservative mass transfer is obtained when the initial mass of donor star is more massive, and the period decreases; hence star move close together.

The Binary Parameters
It is worthwhile to describe the binary parameters used by [4] [5] in their pioneering analysis in some detail, since we are using it as a comparison to our numerical solutions in Section (4.4).We shall thus obtain analytic as well as numerical solutions for a general polytropic index n to compare our results with those obtained from a full numerical evolution.The mass transfer rate for a polytropic donor is given by Jedrzejec (1969) assuming laminar flow, quoted by Paczynski & Sienkiewicz (1972) and from Equations.((51)-( 54)) we obtain:  ) corresponds to the case for which the donor has developed a 0.005M  helium core.
Raising both sides of the Equation (69) to the power ( ) and differentiating, we obtain: where we have set the factor L d R R to unity, given that in most situations Defining a positive dimensionless mass transfer ( ) ( ) , and a characteristic time scale The solution can be easily inverted to yield where is the initial mass transfer rate, normalized as Equation (72).Returning now to Equation (70), with the same definitions as Equation (72) for X and τ, we obtain for the general case in which driving is present ( ) )( ) is the equilibrium value normalized to 0 M .Note that in the stable case, this value is positive; while it is negative in the unstable case.Before we attempt to solve the above differential equation, it is clear from its form and the signs just discussed, that it describes a stable solution in which eq X X → when stable q q < .if, however, stable q q > , the right had side of Equation ( 73) is positive even if the mass transfer vanishes initially, and it just gets bigger as the mass transfer grows.Since X diverges for the non driving case in a finite time, the driven case diverges even sooner.
Considering Equation (70) again we define ( ) ) The differential equation for the evolution of mass transfer now becomes Thus, for the stable case ( * 0 X > , while ( * 0 X < for the unstable case, and * τ is defined positive).The ge- neral analytic solution comprising both the stable and the unstable case can be given in terms of the hyper geometric function, as follows: ( ) Though it is not possible to invert this general solution to obtain the mass transfer rate X as a function of time, particular solutions for specific values of n can be inverted to obtain simple solutions.For example, n = 3/2 yields the [5] solution: ( )

Isothermal Atmospheres
In case of isothermal atmospheres, the mass transfer rate [13] [14] is given by where, α is the scale height.This form of the mass transfer equation is much simpler to integrate than the one for the polytropes considered above.With the same approximations and notation as in the steps leading to Equation (73), and defining , we obtain: ( ) This is easily integrable and invertible to get the following: where ( ) is the time scale required for the driving to change the depth of contact by ~α , and 0 X is the initial value, always positive for physically meaningful cases.In the stable case 0 eq X > , and eq X X → , while 0 eq X < for the unstable case and X diverges in a finite time ( ) . If no driving is present, we may set 0 eq X = and integrate Equation (78) for an isothermal donor.The result is again simple and instructive.
( ) In the stable case, for any initial mass transfer, the system will detach and mass transfer will tend to zero.In the unstable case, any non-zero initial mass transfer will grow and diverge in a finite time.
We specify the parameters of the binary at its initial state; 0.62 Note that initially, the stable and the unstable solutions in Figure 5 almost overlap.This is true as long as the system is either detached or depth of contact is relatively small.As soon as the binary evolves to deeper contact, the unstable solution diverges rapidly and blows up in a finite time.This outcome is predicted by the unstable analytic solution which has been obtained by assuming that the basic primary parameters: the mass ratio q, driving rate L υ and the masses of the components do not change during the evolution.This is obvious- ly not true, especially in the unstable case, since by definition in the unstable case these parameters are evolving rather rapidly (Figure 5).

Numerical Solutions
In Section (4.2 and 4.3), we have derived the basic evolution equations; Equations ((72), ( 73), (75), and (76)) from Equations ((51)-( 54)) in Section (3.1) that for the numerical results and the corresponding analytic solutions for different cases, under certain restrictive approximations.These analytic solutions are a useful reference to compare our numerical results to, since in the limit that the assumptions used to derive the analytic solution are met, the numerical solutions should approach the analytic ones.However, in general, the analytic solutions cannot explain the behavior of a given system accurately and especially when the system evolves rapidly, solving the evolution equations numerically leads to a different outcome than what one would expect analytically.

Limitations of the Analytic Solutions
The analytic solutions of equations in Section (4.2 and 4.3) from Section (2.4) and through Section (3.1) are obtained under the following assumptions: 1) The driving rate given by L υ is constant throughout the evolution.2) The separation "a" of the binary is effectively constant even though the mass transfer rate changes, which is not true for "a" "real system" (see Equation (36)).3) The tidal effects are effectively ignored.4) A system that is initially unstable ( stable q q > ), remains so throughout the evolution, since q is assumed to be the same throughout and it is possible for an initially unstable system to evolve to a stable configuration.
In order to overcome some of these shortcomings one has to numerically integrate the evolution equations in a self consistent way.In particular, in what follows, we allow for the changes in the masses of the components assuming conservative mass transfer (and hence the mass ratio q), allow the binary separation to change as a result of any driving present, and compute the values of d ζ and L ζ as they evolve.The values of d ζ depend on the adopted mass-radius relationship for the donor.To calculate L ζ , we need to specify how the mass and angular momentum are redistributed in the binary during mass transfer which depends on the particular case at hand.For example, it depends on whether the stream impacts the accretor or if an accretion disk is present; if the mass transfer is sub-critical and conservative, or if mass and angular momentum are being ejected from the system following supercritical mass transfer (Figures 6-8). .The mass transfer rate in solar mass normalized to the initial equilibrium rate as a function of time in units of the initial τ: analytic (red curve) and numericalgreen (q = 0.985), blue (q = 0.875), magenta (q = 0.816), cyan (q = 0.656 same as [5]), Darkred (q = 0.625) and orange (q = 0.416).

Figure 7.
Comparison of numerical integrations with the analytic solutions for isothermal atmospheres.The natural logarithm of the mass transfer rate normalized to the initial equilibrium rate for the case with q = 0.663, is shown as a function of time in units of iso τ : analytic (black curve) and numerical green (q = 0.985), orange (0.613), cyan (q = 0.656 same as [5]), blue (0.875), magenta (0.816) and red (0.416).We simulate exactly the same initial conditions as the ones assumed to derive the analytic solution, a constant driving rate and the same mass-radius relationship.We then relax the constraints and let the system evolve in a self-consistent manner for different values of the mass ratio q.We observe that initially all the curves overlap, but as the depth of contact increases, the mass transfer rate increases.As the system evolves, the orbital separation, which is decreasing initially, begins to increase as per Equation (36) which in turn decreases the driving rate.At some point during the evolution q becomes less than stable q and the system evolves to stable mass transfer rates.From our analysis in Section (4.2), we expect the systems with stable q q > to be unstable, and the higher the mass ratio, the more unstable the system is.This is indeed what we observe in the numerical solutions.The numerical solutions reach a peak in the mass transfer rate in a finite time, but turn around and a stable mass transfer is established.In any case, our results suggest and demonstrate that at least some of the binaries that start off unstable do survive the mass transfer instability and evolve into systems like AM CVn-evolving to higher separations and diminishing mass transfer rates.This result, which is not predicted by the analytic solution, has important consequences for population synthesis, understanding AM CVn type systems.
In Figure 7, we plot the analytic and numerical solutions for an isothermal donor for different values of the mass ratio.Unlike the polytropic donor case, which has a fixed reference point where the donor fills its Roche lobe ( d L R R = ), the isothermal case does not have a fixed reference point to which we can normalize our results.We, rather arbitrarily start the integrations when the depth of contact is 11.
, corresponding to an initial mass transfer rate of 10 −5 of the reference rate.The behavior of the systems is qualitatively the same as the polytropic donors.Note that the y-axis in Figure 7 is the natural logarithm of x which from Section (4.3), is just the depth of contact d L R R − in units of the pressure scale height α, normalized to the reference depth corresponding to the equilibrium rate for q = 0.663.
Thus we here produce a numerical analysis for the rate of change of orbital angular momentum(green line), the rate of change of mass of the secondary star (blue line), and Increase in the mass of the primary with time (in year) in a close binary system taking the initial masses of the primary as , p = 30 (per year) and q = 0.8.But afterwards as the mass accretion rate by the primary varies slowly and consequently a m shows a very weak variation with time, the orbital angular momentum orb J continues to slow as steady decreasing with the decreasing in d m .

Conclusions
We have presented a simple, analytical and numerical approach to describe the mass transfer in the evolution of binary star systems in which we track the mass of each component, the total mass of the system, the orbital angular momentum and the spin angular momenta of each component.We have included the effects of mass and angular momentum loss from the system and also the exchange of mass and angular momentum between the components of the binary.The main reason why many binary systems transfer matter at some stage of their evolutionary lifetimes is during the system evolution.One of the stars in a binary system may increase in radius filling its Roche lobe, or the binary separation may shrink because orbital angular momentum is being lost from the system as a consequence of stellar wind mass loss, or gravitational radiation.
There are two general mechanisms for the transferring angular momentum out of the system: First, binary systems emit gravitational radiation, which carries away angular momentum.The rate of angular momentum loss increases as the orbital separation decreases, but it decreases as the total mass of the system decreases.Because both the secondary mass and the orbital separation decrease as CV systems evolving, the net effect is a loss mechanism that does not depend strongly on the orbital period.
Our investigation indicates that the stability of the mass transfer processes depends on: how the radius of the donor star responds to the mass loss, how the orbit responds to mass transfer, and how the mass-gainer responds to the matter that is being dumped on it.
The Roche lobe of the mass donor shrinks as a consequence of its mass loss, increasing the rate at which it loses mass is referred to unstable mass transfer and the Roche lobe of the mass donor grows as a consequence of its mass loss, stopping the mass transfer is referred to stable mass transfer.
In this paper we have made an effort to find a numerical and analytical solutions for conservative and non-conservative mass transfer in close binary systems which can address the gradually reducing rate of the mass accretion by the accretor from the donor with respect to time as well as with respect to the increase in mass of the accretor [15].In view of this we have presented the time dependent profiles of the masses of both the component stars in a close binary as well as the orbital angular momentum of the system.Our analytical solutions shows for the non-conservative mass transfer in close binary system taking the initial masses of the primary and secondary as corresponds to the case for which the donor has developed a 0.005M  helium core (see Figure 4), the orbital angular momentum J continues to show a steady decreasing with the increasing time (see Figure 2).
In conservative mass transfer and from Equation (15); if a m is the mass losing star, a m  will be negative, therefore we conclude that if the mass donor star is more massive, then the orbit will shrink and hence period decreases.But if the mass donor is less massive than the accretor component then orbit will be widened and period increases.
We have also presented a new finding (analytic and numeric solutions) for the evolution of binary systems, generalizing the results of [4] to binaries with any polytropic index and to binaries with isothermal atmospheres.The analytic solutions always predict that if a binary attains contact when the mass ratio q is such that the ensuing mass transfer is unstable, the binary merges in a finite time.On the other hand, the numerical solutions predict that a binary that undergoes unstable mass transfer at initial contact can evolve to a state of stable mass transfer, survive and evolve into an AM CVn type system.Moreover, the analytic solutions also discussed here assume that the driving rates d υ and L υ in Equation (70) remain constant while the depth of contact changes.This is only approximately true; and a self-consistent solution will require numerical integrations.This is a situation we encounter in some large-scale hydrodynamic simulations of mass transfer in poly tropic binaries.We note that from Equations ((54), and (72)) in the stable case, ( ) , the mass transfer decays asymptotically to zero over a characteristic time ( ) ( ) ( ) , it will blow up in finite time equal to chr τ .Thus, the essence of the stability of mass transfer in a binary is already contained in the sim- ple case of no driving.The presence of driving exacerbates the natural instability or, in the stable case, it settles asymptotically to a non-zero stable mass transfer, which we have been observed in AM CVns, CVs and LMXBs.In our study numerical solutions indicates more physical descriptions than analytical (Figures 1-5) one as shown in Figures 6-8.
Finally, the most useful thing, is that if one assumes only accretion and wind-like mass transfer, then most every binary in which mass is transferred from the less massive star is stable on both dynamical and thermal time scales.If the mass donor has a radiative envelope (not treated here), it will shrink in response to mass loss, and lose mass stable.If the donor has a convective envelope, a modestly sized core will stabilize it sufficiently to prevent mass loss on the dynamical time scale.
= for comparison in Figure 1 and Figure 2.

Figure 1 .Figure 2 .
Figure 1.Comparison of [10] (Eq.(60), blue curve) and [11] (Eq. (61)) formulae for the Roche radius, normalized to the binary separation vs. mass ratio change in mass for both the accretion and donor as well as the changing of the orbital angular frequency with time are derived in Equation (21) within the pre-assigned time scale of mass exchange as 4 10 t =

Figure 4 .
Figure 4.The total rate of change of the orbital angular frequency of the binary (Eq.(30)) related to the mass lost from the system of two point masses which will coalesce due to gravitational wave emission in a time interval shorter than 10 10 yr, as a function of the donor mass in solar unit ( ) d m M  .As we go from top to bottom on the right side; the lines (curves) are calculated for magenta ( ( ) 10 10 M M BH BH + +   Equations ((72), (73), (75), and (76)).We notice that since > d L d ζ ζ ; the ra- dius of the star increases at a rate faster than the Roche lobe radius, and the resulting mass transfer is unstable.The analytic solutions, both the stable and the unstable, are shown in Figure 5.

Figure 5 .
Figure 5.The solutions for the stable (dotted line) and unstable cases.

Figure 6 .
Figure 6.Numerical and Analytic Solutions for Polytropes.Comparison of numerical integrations with the analytic solution in Section (4.2).The mass transfer rate in solar mass normalized to the initial equilibrium rate as a function of time in units of the initial τ: analytic (red curve) and numericalgreen (q = 0.985), blue (q = 0.875), magenta (q = 0.816), cyan (q = 0.656 same as[5]), Darkred (q = 0.625) and orange (q = 0.416).

Figure 8 .
Figure 8. Loss of orbital angular momentum and decrease in the mass of the secondary with time (in year).
the present calculation, we consider the values of the parameters as 2.78 a = an initial 1M  donor whose central hydrogen abundance is 0.005 at the onset of mass transfer, donor whose central hydrogen abundance has just reached zero at the onset of mass transfer and 1