cataclysmic variables ( cvs ) are close binary systems consisting of a low - mass main - sequence ( ms ) star that fills its roche lobe and a white dwarf .
main - sequence star ( mass - donating star , usually referred to as star 2 or secondary ) loses matter through the vicinity of the inner lagrangian point @xmath0 .
white dwarf ( accretor star , referred to as star 1 or primary ) accretes at least a fraction of this matter via accretion disk or via accretion columns in the polar zones in the case when white dwarf has a strong magnetic field .
physics and evolution of cvs as well as evolution of similar to them low - mass x - ray binaries were investigated starting from the late sixties ( see , e.g. , [ [ rjw82][fed97 ] ] and references therein ) .
these studies were motivated by the recognition of the fact that evolution of cvs is determined by losses of orbital angular momentum due to radiation of gravitational waves [ [ kmg62],[p67 ] ] and magnetically coupled stellar wind [ [ sch62][vz81 ] ] .
a number of studies investigated the influence of the loss of angular momentum associated with the material outflow from the system on the cv evolution ( see , e.g. , [ [ yu73][kk95 ] ] ) .
however , in the absence of gas dynamical simulations of the mass transfer in the binaries , these processes were considered under the parametric approach .
recent three - dimensional ( 3d ) gas dynamical simulations of the structure of gaseous flows in semi - detached binaries proved an important role of a circumbinary envelope ( see , e.g. [ [ bbkc98b],[bbkc98c ] ] ) .
these calculations also show that during mass exchange a significant fraction of matter leaves the system . the present work is mainly devoted to the numerical investigation of the evolution of cvs using the data of 3d gas dynamical calculations on the losses of mass and angular momentum from close binaries .
the main attention is paid to the stability of mass exchange against runaway mass loss and its dependence on the donor to accretor mass ratio @xmath1 .
it is known that the loss of matter by the donor star results in the violation of its hydrostatic and thermal equilibrium .
hydrostatic equilibrium is restored adiabatically , i.e. in dynamical time scale , while thermal equilibrium is restored in thermal time scale .
this transition to a new equilibrium state leads to a change of donor radius .
the sign of radius variation depends on the convective and radiative stability of outer envelope of the star . for stars with masses
@xmath2 with deep convective envelope as well as for white dwarfs mass loss results in increase of the radius , while for stars with radiative envelopes it leads to the shrinkage of the star .
mass exchange in close binary is unstable when in the course of evolution mass - losing star tends to overfill the roche lobe .
it can occur when radius of the donor @xmath3 increases faster ( or decreases slower ) than the effective radius of roche lobe @xmath4 .
is determined as the radius of a sphere with a volume equal to the volume of the roche lobe . ]
the case when the radius of the donor increases while the radius of the roche lobe decreases is also possible .
thus , the question of stability of mass exchange against runaway is determined by the balance of derivatives @xmath5  and @xmath6 .
it means that unstable mass exchange is possible even for stars which contract when losing matter .
the study of conditions of stable mass exchange is motivated by two related problems : i ) it is necessary to explain why @xmath7 10% of cvs with well - determined masses of components have combination of donor mass and mass ratio of components which is `` forbidden '' in the evolutionary models based on conventional assumptions on the conservation of angular momentum of a binary ; ii ) for population synthesis studies of cvs it is necessary to distinguish progenitors of cvs among all binaries containing a white dwarf and a low - mass companion ( i.e. , specify in which binaries of this type stable mass exchange is possible ) .
a problem similar to the last one appears in the studies of low - mass x - ray binaries .
the paper is organized as follows . in section  2
we describe major factors affecting the evolution of cvs and conditions of stable mass exchange in binaries .
the reasons which motivated us to rule out the conservative approximation for mass exchange are considered in section  3 .
sections  4 and  5 describe some results of 3d gas dynamical simulations of the flow structure in semi - detached binaries and introduce the model for description of angular momentum losses due to non - conservative mass exchange .
results of evolutionary calculations for conservative and `` non - conservative '' models are compared in section  6 .
criteria of stable mass exchange for `` non - conservative '' model with mass and angular momentum loss are determined in section  7 .
section  8 summarizes the main results of the work .
in the context of the present work the term `` stable mass exchange '' means dynamical stability . in other words ,
in the course of mass exchange the donor star can be out of thermal equilibrium and mass transfer rate may exceed the rate corresponding to the kelvin time scale ( @xmath8 , where radius @xmath9 , luminosity @xmath10 , and mass of the star @xmath11 are in solar units and mass transfer rate @xmath12 is in @xmath13 ) .
let us define unstable mass exchange as the situation when the radius of the donor changes faster than the effective radius of roche lobe : @xmath14 . using the derivatives of the radii w.r.t .
mass of the donor we can rewrite the condition of the stability of mass transfer as :    @xmath15    the effective radius of roche lobe can be estimated , for instance , using interpolation formula given by eggleton [ [ eggl83 ] ] :    @xmath16 where @xmath17 is the semimajor axis of the orbit .
for @xmath18 the approximation suggested by paczyski [ [ pacyn71 ] ] is more convenient :    @xmath19    derivative of the roche lobe radius can be written as    @xmath20    it is seen from ( [ rrp ] ) and ( [ drra ] ) that the variation of effective radius of roche lobe is defined by the change of donor and accretor masses @xmath21 and @xmath22 as well as by the change in the separation @xmath17 , and , hence , it is influenced by the loss of mass and angular momentum from the system . as a rule , only systemic angular momentum is considered as the orbital momentum , i.e. the sum of angular momenta of two stars , which are considered as point masses .
momenta of the spin of the components , momentum of the accretion disc ( if it s present ) , as well as momentum of gaseous flows inside the system as a rule are not included into the total momentum .
deviation of the momentum of the donor from that of the point mass object ( which can be significant , see [ [ kuz95 ] ] ) usually is not taken into account as well .
such a simplified approach is compelled , since it is very difficult to include all the abovementioned factors into calculations .
having in mind all these simplifications , one gets for the orbital momentum of a binary system with the circular orbit    @xmath23    where @xmath24 is the gravitational constant .    according to the widely accepted models ,
the evolution of cataclysmic binaries is determined by the loss of angular momentum from the system via gravitational waves radiation ( gwr ) and/or magnetic stellar wind ( msw ) of the donor as well as the mass exchange between components . in the standard models of evolution it s accepted that the mass exchange among components does not change the systemic angular momentum , and influences it indirectly  through a possible loss of mass from the system and angular momentum ablation by leaving matter ( term with subscript loss in eq .
[ dj ] below ) .
therefore , within conventional models , the equation for the change of orbital momentum can be written as :    @xmath25    let us consider the terms of expression ( [ dj ] ) : +    _ 1 ) loss of the systemic angular momentum due to gwr _ + the change of the systemic orbital momentum as a result of gwr is given by a formula ( see , e.g. , [ [ ll62 ] ] ) :    @xmath26    where @xmath27 is the speed of light .
intensity of gwr is characterized by a very strong dependence on the separation of components and respectively on the orbital period .
this process is essential for short - period systems with @xmath28 , since only then the time scale of angular momentum loss becomes shorter than the hubble time .
+    _ 2 ) loss of the systemic angular momentum due to msw of the donor _ + a mechanism for the loss of the angular momentum of the system by means of the msw was suggested by schatzman [ [ sch62 ] ] and by mestel [ [ mes68 ] ] .
if the star has a convective envelope and , subsequently , a surface magnetic field , its own axial rotation is braked by the magnetic stellar wind , and the angular momentum loss rate can be essential even for small mass loss rate . the subsequent synchronization of the spin of the donor and orbital rotation caused by tidal interaction between components results in the loss of the orbital momentum of the system and reduction of the separation of components @xmath17 . to take this effect into account , verbunt and
zwaan [ [ vz81 ] ] extrapolated the data on rotation of f - type stars [ [ skum72 ] ] to the k and m type components of cvs and suggested a widely used semi - empirical formula which allows to determine the temporal behavior of the orbital momentum :    @xmath29    were @xmath30 is the radius of gyration of the donor star ; @xmath31 is a numerical factor determined from the comparison of theoretical calculations and observational data .
well - defined reduction of magnetic field for stars with masses less than @xmath32 permits to suggest that vanishing of a radiative core of a star when its mass decreases to this threshold leads to an abrupt `` switch - off '' of the so - called `` @xmath33''-dynamo mechanism that is responsible for the generation of stellar magnetic fields [ [ sr83 ] ] . at the instant of cessation ( or sharp reduction ) of msw
the time scale of angular momentum loss [ defined by eq .
( [ djmsw ] ) ] is shorter than the thermal time scale of a donor with @xmath34 .
thus , the latter is out of thermal equilibrium and its radius is larger than the radius of thermally equilibrium star of the same mass.,[pol96 ] ] . ]
when the action of msw and corresponding loss of angular momentum stop , the rate of reduction of the semimajor axis declines , the donor shrinks down to the equilibrium radius and loses the contact with roche lobe [ [ sr83 ] ] .
meanwhile , the system continues to lose angular momentum due to gwr , components continue to approach each other and after some time the optical star fills its roche lobe again .
further evolution of the system is determined by the loss of angular momentum via gwr .
cessation of msw after the star becomes completely convective is widely adopted hypothesis explaining satisfactorily the so - called orbital period gap in cvs .
the value of coefficient @xmath31 in eq .
( [ djmsw ] ) can be determined by comparison of theoretical width of period gap and observational one . in the present study we adopted @xmath35 in accordance to [ [ ft94 ] ] .
+    _ 3 ) loss of the systemic angular momentum due to mass loss _ + mass loss from the system is mainly considered as a parameter fine tuning of which permits to get an agreement with observational minimal period of cvs ( e.g. [ [ rjw82 ] ] ) .
more definite assumptions are made only for specific cases when examining evolution of the systems , in which : i ) the rate of accretion on white dwarf is limited by the hydrogen burning rate ( @xmath36@xmath13 ) , and the mass loss rate by the donor is much higher than this limit , but does not exceed eddington limit for the dwarf ( that is close to @xmath37@xmath13  , see , e.g. , [ [ hkn96],[yl98 ] ] ) , the mass excess being lost by means of stellar wind [ [ kh94 ] ] ;  ii ) mass is lost due to outbursts occurring on white dwarf ( see , e.g. , [ [ yu73 ] , [ yuetal96][skr98 ] ] ) . in both cases
it is usually assumed that the specific angular momentum of the matter leaving the system is equal to the specific momentum of the accretor .    to describe the loss of mass and momentum out of the system
, it is convenient to introduce parameters describing the degree of non - conservativity of the evolution respective to the mass [ [ ty71 ] ] :    @xmath38 and to the angular momentum : sometimes specific angular momentum lost from the system is measured in units of @xmath39  and used instead of @xmath40 : @xmath41 ]    @xmath42    where @xmath43 .    with these parameters ,
the expression for the loss of momentum by the matter leaving the system gets the following form :    @xmath44    in the calculations of the cataclysmic binary evolution an equation for the variation of the semiaxis of the orbit is used directly , instead of the equation for the variation of momentum over time .
let us differentiate the expression ( [ j ] ) , substitute the result into equation ( [ dj ] ) , and use the relation @xmath45 .
then the equation for the variation of orbital semiaxis can be written as :    @xmath46    here , index ` exch ' corresponds to the change of @xmath17 resulting from the mass transfer between components .
note that while the term @xmath47 is absent in the equation ( [ dj ] ) , the appropriate term @xmath48 is present in the equation ( [ da ] ) , being a natural consequence of the dependence of the orbital momentum both on semimajor axis of the orbit and masses of components .
when mass exchange is conservative , the evolution of the orbital separation is determined by the assumption that the mass transfer does not change the orbital momentum of the system . in the course of mass exchange
the matter that initially had specific angular momentum of the donor is transferred onto accretor and finally gets the specific momentum of the accretor .
the assumption of conservation of the orbital momentum implies that if the mass of the donor is lower than the mass of accretor , the excess momentum of accreted gas should transform into the orbital one , and the mass exchange acts in the direction of increasing the orbital semiaxis .
if the donor is more massive than the accretor , the lacking momentum of the accreted gas should be taken from the orbital momentum , therefore , the mass exchange acts in the direction of shrinking of the binary system .
the numerical investigation of the evolution of cataclysmic binaries consists of simultaneous calculations of the evolution of the donor and of the variation of the orbital separation in time .
let us consider the processes that determine the evolution of the orbital semiaxis .
using parameters @xmath49 and @xmath40 , we can write formulae for the terms of equation ( [ da ] ) ( here masses and distances are given in solar units and time is given in years ) :    @xmath50    @xmath51    @xmath52    @xmath53    it is seen from eqs  ( [ daex])([damsw ] ) that changing of @xmath4 is mainly determined by the value of @xmath54 . on the other hand ,
changing of @xmath3 depends on @xmath21 and @xmath55 .
thus the most convenient way for studying of limits of stable mass exchange is analysis of `` ratio of donor - star mass to accretor mass @xmath54  donor - star mass @xmath21 '' diagram .    in a number of cases
one may simplify expression ( [ drra ] ) and obtain analytical formulae for the derivative of the roche lobe radius .
in particular , when the angular momentum changes due to gwr and msw can be neglected , one may deduce from ( [ daex ] ) and ( [ daloss ] ) the expression for the change of separation @xmath17 [ the first term in the formula ( [ drra ] ) ] as a function of mass ratio @xmath54 and the parameters of the loss of matter and momentum @xmath49 and @xmath40 :    @xmath56    using the expression ( [ rrp ] ) for @xmath4 we can express the second term of ( [ drra ] ) as :    @xmath57    using ( [ dabpsi ] ) and ( [ dabpsi1 ] ) , for the case of totally conservative mass exchange ( @xmath58 ) , we obtain :    @xmath59    from the same equations , for the so - called ` jeans mode ' of mass loss , when the mass exchange in the system does not occur , but the donor loses its mass due to stellar wind and matter leaves the system carrying away specific angular momentum of the donor ( @xmath60 , @xmath61 ) , we obtain :    @xmath62    for the case of stellar wind from the accretor when all matter lost by the donor is transferred onto accretor but then leaves the system carrying away specific angular moment of the accretor ( formally , one should take here @xmath60 and @xmath63 ) we obtain :    @xmath64    the comparison of the rates of variation of the radius of the donor and effective radius of the roche lobe ( [ drra ] ) determines the state of stability / instability of mass exchange in any phase of evolution .
it is possible to determine the boundaries of stable mass exchange region by two different methods :    * it is possible to calculate derivatives of the radius of the donor w.r.t .
its mass directly for different @xmath55 and to compare it with derivatives of the effective radius of roche lobe depending on @xmath22 , @xmath21 , and @xmath17 .
this method was used , for instance , in [ [ tfy82 ] ] and it permits to restrict a region of @xmath65 diagram where mass exchange is dynamically stable and also to detect the regions where mass exchange should occur in different time scales : thermal , nuclear , or in the time scale of angular momentum loss . * alternatively , a requirement that the rate of accretion should not exceed the eddington limit for the dwarf can be accepted as a criterion of stable mass exchange .
this assumption seems to be more reasonable since evolutionary calculations often suggest very high rate of mass loss during very initial moments of time after overfilling the roche lobe , while after that @xmath12 falls off rapidly .
in fact , this approach involves a restriction on dynamical stability of mass transfer but permits loss of matter in thermal time scale .
figure  [ ris1 ] depicts the boundary of the stable mass exchange region in conventional conservative model of the evolution of cv .
boundary 1 was found in [ [ tfy82 ] ] , and boundary 2 was found in [ [ kool92 ] ] .
these boundaries coincide well for of donors of low mass but differ a little for large ones .
unfortunately , in [ [ kool92 ] ] the procedure of the calculation of the boundary is nt described in detail .
we may assume that this discrepancy is due to the differences in the codes , in a different treatment of possible weakening of donor msw at masses greater than @xmath66 , and different techniques for computation of stellar radii derivatives in this region .
figure  [ ris1 ] also shows cvs with known masses of components from the catalog [ [ cat97 ] ] .
the total number of binaries in the diagram is 80 , and 11 of them are in the region of unstable mass exchange . for 8 of these 11 systems uncertainties of parameters
are known , they are shown in fig .
[ ris1 ] as error bars .
it is seen that the location of stars in a region of unstable mass exchange can not be explained by uncertainties in masses of components . in fig .
[ ris1]a we plot cvs with the best determined masses [ [ smith_dhillon ] ] , parameters for a number of systems being different compared to the data from [ [ cat97 ] ] . nevertheless , even for well - determined systems the problem of `` improper '' location of cvs remains  3 out of 22 systems lie in the `` forbidden '' region .
to explain this disagreement we suggest a model of cvs evolution which involves loss of mass and angular momentum from the system .
this model is based on 3d gas dynamical simulations of the flow structure in close binary systems .
the above - stated analysis proves that the rate of change of the distance between components @xmath17 in evolutionary calculations depends heavily on two parameters : the degree of non - conservativity of mass exchange w.r.t .
the matter @xmath49 and w.r.t . to the angular momentum @xmath40 .
these parameters can be determined only by means of 3d gas dynamical simulations of mass transfer in binary systems .
these investigations were made in [ [ bbkc98b ] ] for a binary system with mass ratio of components @xmath67 , and in [ [ bbkc98c ] ] for a binary with @xmath68 . besides that
, we have specially made a numerical simulation of a binary system with @xmath69 to study the influence of parameter @xmath54 on the flow structure .
analysis of the obtained results proves that in the binary systems considered by us for different values of parameter @xmath54 mass transfer is non - conservative and the non conservativity degree w.r.t .
the matter is @xmath70 .
simulations of flows in the 3d numerical models prove also that the matter leaves the system with a significant angular momentum .
the non conservativity degree of mass exchange w.r.t .
the matter can be easily determined directly from gas dynamical model , while the determination of non conservativity degree w.r.t .
the angular momentum is a sophisticated problem .
let us consider the equations which determines the transport of angular momentum in the system .
we should consider navier - stokes equations for a correct analysis because viscosity plays a significant role in the redistribution of angular momentum . from the steady - state gas dynamical equations in a rotating frame
we can deduce the equation for the angular momentum transfer    @xmath71 where , as usually , @xmath72 is the velocity vector , @xmath73 is the density , @xmath74 is the roche potential , @xmath75 is the orbital angular velocity , @xmath76 is the stress tensor , @xmath77 is the radius - vector of the center of masses of the system , and angular momentum @xmath78 ( in a laboratory , i.e. inertial or observer s frame ) is defined by expression :    @xmath79
let us write down the equation for angular momentum transfer in a divergent form :    @xmath80 from this equation , we can obtain the integral law of change of angular momentum as follows :    @xmath81 the term @xmath82 determines the change of the angular momentum due to non central nature of the force field , defined by the roche potential . by definition ,
the value    @xmath83    is the specific flux of angular momentum @xmath78 . as a result ,
for the steady - state case we obtain :    @xmath84 where @xmath85 and @xmath86 are the boundaries of the donor and the accretor , respectively , and @xmath87 is the outer boundary of the calculation domain . by analogy to the flux of matter    @xmath88 one may estimate the flux of angular momentum from the system as    @xmath89    and to use it in the expression ( [ eqpsi ] ) , determining the non conservativity parameter w.r.t . the angular momentum .
note that numerical simulations were conducted for eulerian equations for in - viscid gas and not for navier - stokes equations . respectively , to calculate integral ( [ error ] ) we use isotropic term describing gas dynamical pressure in the expression for the specific flux of angular momentum : @xmath90 .
is kronecker s symbol : @xmath91 ] this substitution is correct because on the outer border where the integral in ( [ error ] ) is estimated , viscosity does not play a significant role .    using expression ( [ error ] ) in the gas dynamical calculations we obtained @xmath92 corresponding to @xmath93 ( see definition of @xmath94 in the footnote on page ) for the binary system with @xmath67 [ [ bbkc98b ] ] , and @xmath95 that corresponds to @xmath96 [ [ bbkc98c ] ] for a binary with @xmath97 the last value is in good agreement with the results of two - dimensional calculations of binary system with equal masses of components [ [ sawada84 ] ] where @xmath98 was obtained .
however , the use of these estimates for @xmath40 in evolutionary calculations proves that this rate of an angular moment loss is so great that in the majority of binary systems the mass exchange rapidly becomes a runaway process and the rate of mass loss by the donor increases without any limit .
apparently , the application of formula ( [ error ] ) to the estimation of the angular momentum loss in the evolutionary calculations is not quite correct . since conjunction of gas dynamical and evolutionary models is not completely consistent .
the considered gas dynamical model does not take into account the change in the position of stars in time and , hence , the change in the angular momentum of the gas due to the non - central nature of the force field is not compensated by respective change in the orbital angular momentum .    in a more general gas dynamical model , in which the change in the position of stars is taken into account , one may correctly estimate the change of the systemic angular momentum in the form of angular momentum flux through the external boundary for the `` donor + accretor + gas '' system .
however , to take into account such gas dynamical calculations appears rather complicated .
such estimates are possible only for a specific stage in the life of a binary system . moreover , thus obtained gas dynamical results can not be used directly in the standard evolutionary models that ignore the presence of circumbinary envelope in the system . therefore , a simplified model for estimation of the angular momentum loss developed on the basis of gas dynamical calculations will be a promising approach in the nearest future .
gas dynamical simulations of mass transfer in semi - detached binaries confirm that outflow of matter from the donor occurs through a small vicinity of the inner lagrangian point @xmath0 . in this case
specific angular momentum of lost matter ( in corotating frame ) can be estimated as @xmath99 , where @xmath100 . respectively , the flow of angular momentum from the donor is equal to    @xmath101    due to non - conservativity of mass exchange , only a fraction @xmath49 of the matter outflowing through @xmath0 is accreted . specific angular momentum of accreted matter is equal to the angular momentum of accretor ( here we neglect the finite radius of accretor and/or residual angular momentum of accreted matter , and this approximation appears to be reasonable for cvs ) .
it means that accreted matter has zero angular momentum in the rotating frame , and , hence , does not accelerate the rotation of accretor .
thus the question on the efficiency of the angular momentum transfer from the spin rotation of accretor to the orbital rotation of the binary can be ignored . to summarize
, we can describe the flow of accreted angular momentum as    @xmath102    these expressions for the angular momentum flows do not contradict the general formula for specific flow of angular momentum [ see ( [ potok ] ) ] . indeed , viscosity does not play a significant role neither in the vicinity of the surface of the donor nor in the vicinity of the accretor where the matter has already lost its angular momentum and the flow has radial direction in the rotating frame .
it means , that in these cases the stress tensor @xmath76 is reduced to the isotropic pressure as well , and does not deposit into the integral of specific flow of angular momentum .
our assumptions permit to write the expression for the flow of angular momentum lost from the system in the following form :    @xmath103 where @xmath104 is a parameter defining the fraction of the circumbinary envelope angular momentum that is lost with the matter leaving the system .
then @xmath105 is the fraction of the circumbinary envelope angular momentum that returns back into the system via tidal interaction .
henceforth we adopt for the parameter @xmath104 a value of 1 that corresponds to the case when the angular momentum of the circumbinary envelope is totally lost with the matter leaving the system . in this case the specific angular momentum of the matter leaving the system is equal to    @xmath106 or in the units of systemic specific angular momentum ( @xmath107 ) :    @xmath108    where @xmath109 is dimensionless distance from the center of the donor to @xmath0 .
we would like to stress that formula ( [ xxxx ] ) was derived under the assumption of aligned synchronous rotation of components of the binary system .
it is known that the loss of angular momentum due to outflow of the matter from the vicinity of @xmath0 can result in the development of a misalignment of the vector of spin rotation of the donor and the vector of the orbital rotation of the binary [ [ matese83 ] ] .
the flow structure in the binaries with misaligned asynchronous rotation was studied by the authors in [ [ bbkc98d ] ] .
having in mind the simplified character of evolutionary models we do nt consider the effect of the misalignment in the present work .
the graph of @xmath110 is shown in fig .
[ ris2 ] for various values of non - conservativity degree w.r.t .
the matter @xmath49 .
the values @xmath111 and @xmath112 are also shown in the figure by dashed lines . in the region @xmath113 each unit of the matter leaving the system carries away an angular moment that is higher than the mean specific angular momentum of the system .
this results in diminishing of systemic specific angular momentum :    @xmath114    in the region @xmath115 each unit of the matter leaving the system carries away an angular moment that is lower than the mean specific angular momentum of the system , and , hence ,    @xmath116    figure  [ ris2 ] shows that the value of @xmath49 becomes negative in a certain range of @xmath54 for the majority of cases of non - conservative mass exchange ( except @xmath117 ) .
this means that the escape of the matter to infinity causes redistribution of the angular momentum in the `` donor + accretor + gas '' system .
this results in the growth of the total angular momentum of the binary : @xmath118 .
it is evident that such a `` pumping '' of angular momentum into the system can stabilize the mass exchange in cv .
indeed , consideration of orbit semiaxis change due to conservative mass exchange shows that the binary system widens for @xmath119 only . for the `` non - conservative '' model
the range of @xmath54 corresponding to the widening of binary due to mass and angular momentum losses is larger  for @xmath120 it happens for @xmath121 .
note that @xmath122 dependence has a curios feature [ see fig .
[ ris2 ] and formula ( [ xxxx ] ) ]  all curves @xmath110 for different values of @xmath49 pass through the same point .
this point corresponds to the case when the center of mass of the binary system is located exactly equidistantly between @xmath0 and accretor , and in this case @xmath40 does not depend on @xmath49 .
fortunately , this specific feature does not influence the analysis of solution .
three - dimensional gas dynamical simulation of mass transfer in cvs confirms conclusion that the outflow of matter from the system can lead both to decrease and to increase of the total angular momentum of a binary .
the typical flow structures [ see fig .
[ ris3]a , b ] for binaries with @xmath123 and @xmath124 demonstrate that in an inertial frame gas can have different directions of rotation ( in both figures orbital rotation corresponds to the counter - clockwise direction ) . respectively , the change of angular momentum has opposite signs in the obtained solutions  for @xmath67 the systemic angular momentum decreased , and for @xmath69 it increased .
it seems reasonable to explain the different types of gas dynamical solutions using the following qualitative speculation :    * the flow of the gas leaving the system consists of the matter of the circumbinary envelope that has an angular momentum large enough to overcome the gravitational attraction of the accretor .
initial velocity of this matter ( in rotational frame ) has the same direction as the orbital movement of the accretor .
* velocity of the gas leaving the system changes under the action of the gravitational attraction of both components of the binary system , centrifugal force , coriolis force , and pressure gradient .
the change of azimuthal velocity is determined mainly by the coriolis force and gravitation forces , since centrifugal force is axially - symmetrical and deviation of pressure gradient from the axial symmetry is small as well .
action of coriolis force results in deflecting of the stream in the direction opposite to that of the orbital movement .
consequently , this permits to subdivide the whole flow leaving the system through the vicinity of @xmath125 into two parts : the first part retains the direction of original movement , and the direction of movement of the second one is changed .
* it is evident that in an axially - symmetrical gravitational field the action of coriolis force can not produce a solution when the direction of gas movement in the inertial frame is opposite to the orbital one .
the situation is changed drastically for the case of non - axially - symmetrical gravitational field of a binary system . in this case gravitational attraction of the donor
can accelerate the gas of the stream and result in the solution with opposite direction of movement in the inertial frame .
* it follows from eq .
( [ xxxx ] ) and fig .
[ ris3 ] that there is a range of values of @xmath54 where gas moves in the direction opposite to the orbital one in the inertial frame .
these values of @xmath54 correspond to the cases when attraction of the donor is large enough to accelerate the gas properly ( left end of the range ) , but not too large to produce accretion of gas on the surface of the donor ( right end of the range ) .
the values of @xmath54 out of this range correspond to the case when the flow directed along the orbital movement prevails .
appearance of an additional range of values of @xmath54 where the `` non - conservative '' mass exchange is accompanied by increase of systemic angular momentum results in expansion of the stable mass exchange region . to illustrate this fact
let us consider a graph that is often used in the qualitative studies of the stability of mass transfer in cv .
this is a @xmath126 graph with curves corresponding to the analytical expressions for derivatives of effective radius of the roche lobe w.r.t .
the mass of the donor @xmath127 , and straight lines corresponding to the values of derivatives of the radius of the star w.r.t . to its mass @xmath128 ,
the latter being known from the investigations of the internal structure of star ( see , e.g. , [ [ tfy82 ] ] ) .
mass exchange is stable for the range of @xmath54 where curve @xmath127 passes below the straight line @xmath128 for the donor of corresponding type .
this graph is presented in fig .
it shows the derivative of effective radius of roche lobe @xmath129 as a function of @xmath54 for a scenario in which the loss of angular momentum is defined by formula ( [ xxxx ] ) for various values of the parameter @xmath49 .
dependencies @xmath130 for the scenario of the completely conservative mass exchange and for the case when the outflowing matter carries away specific angular momentum of the donor or specific angular momentum of accretor are shown as well . the same fig .
[ ris4 ] also shows the values of derivative @xmath131 for completely convective and degenerate stars ( @xmath132 , for thermally equilibrium stars of ms with @xmath133 ( @xmath134 , and for subgiants with degenerate low - mass helium cores ( @xmath135 . the analysis of fig .
[ ris4 ] proves that the outflow of matter and loss of angular momentum from the system in accordance to law ( [ xxxx ] ) stabilize mass exchange in a wider range of @xmath54 compared to the conventional models of conservative mass exchange .
it is also worth noting that for conservative scenario of mass exchange its time scale for stars with masses @xmath136 is determined by the time scale of the loss of angular momentum or the nuclear time scale for @xmath137 [ [ tfy82 ] ] , while for above - considered `` non - conservative '' model this limit on @xmath54 is sufficiently higher .
we would like to stress also that in the case under consideration there exists a region of unstable mass exchange where value of @xmath54 is very small .
the reason of instability [ i.e. the violation of criteria ( [ dr2r ] ) ] is fast increase of the ratio of the lost angular momentum to the mean specific momentum of binary system [ see eq .
( [ dabpsi ] ) and ( [ xxxx ] ) ] .
this circumstance can lead to the destruction of a donor of very low mass
. is also known .
it may be caused by low efficiency of interaction between accretion disk and orbital movement ( see , e.g. , [ [ rsh83],[hp84 ] ] ) .
] thus , cvs that begin their evolution in a stable mode finish it by a catastrophic disruption of the donor when @xmath54 becomes lower than a certain threshold .
to take into account results of 3d gas dynamical simulations of the mass transfer in cvs and to incorporate proper estimates of the loss of angular momentum we investigated the evolution of cvs on the basis of various assumptions on its conservativity .
evolution of the donor was calculated using a modified version of the code for numerical modeling of low - mass stars which was used in the previous works of the authors ( see , e.g. , [ [ ft94],[fed97],[fy84][tf89 ] ] ) .
the code uses opacity tables of huebner _
et al . _ [ [ hmma77 ] ] with addition of data of alexander _ et al . _
[ [ ajr83 ] ] for low - temperature region .
equation of state was adopted from fontaine _ et al .
_ [ [ fgh77 ] ] with modifications suggested by denisenkov [ [ den89 ] ] . rates of nuclear reaction were adopted from [ [ hfcz83],[cfhz85 ] ] .
all runs were carried out under assumption that the donor fills its roche lobe immediately after arriving on zams .
masses of the donors were adopted in the range from @xmath138 to @xmath139 and for various values of @xmath54 .
chemical composition was adopted as follows : @xmath140 , @xmath141 , @xmath142 . for the mixing length parameter
a value of @xmath143 was adopted .
we used the method of calculation of the donor mass loss rate for known stellar radius and effective radius of roche lobe suggested by kolb & ritter [ [ kr90 ] ] .
based on the results of gas dynamical calculations , we assumed that the fraction of matter accreted by white dwarf is @xmath120 , and that angular momentum carried away by the matter leaving the system is defined by expression ( [ xxxx ] ) .
we did nt include into the model mass and angular momentum loss due to outbursts .
let us consider evolution of cv with parameters corresponding to the stable mass exchange in conservative model . in fig .
[ ris5 ] ( `` orbital period  logarithm of donor - star mass loss rate '' diagram ) two evolutionary tracks are shown .
one of the tracks is calculated in the conservative approximation ( track 1 ) , and another one  in the `` non - conservative '' approximation ( track 2 ) .
initial masses of both donor and accretor were equal to @xmath144 .
it is seen that these two tracks are almost identical . during the early phase evolution
is driven by angular momentum loss due to the msw of the donor .
this results in decrease of the orbital semiaxis and orbital period in the course of evolution .
decrease of the mass of the donor as well as its radius lead to decrease of @xmath145 [ see eq .
( [ damsw ] ) ] , and consequently to decrease of mass loss rate by the donor @xmath55 .
for the `` non - conservative '' track , the value of @xmath55 is remarkably lower than for the `` conservative '' one , because in this phase the value of @xmath40 is negative ( see fig .
[ ris2 ] ) and net loss of orbital angular momentum in the `` non - conservative '' model is lower than in the conservative one .
the situation is inverted in later phases of evolution as @xmath54 reaches the values corresponding to large angular momentum loss .    after the star becomes completely convective , msw from the donor stops .
mass of the star at this instant is determined by the deviation from the thermal equilibrium .
single ms stars become completely convective at @xmath146 , while mass - losing components of binaries do at lower masses  @xmath147 ( under adopted assumption on chemical composition and opacity of stellar matter ) .    in the `` conservative '' track mixing occurs at @xmath148 ( when the orbital period is equal to @xmath149 ) , and in the `` non - conservative '' track  at @xmath150 ( @xmath151 ) . while in reality
in this case the evolution of binary has to be terminated by a thermonuclear explosion of accretor .
the continuation of calculation for @xmath152 can be advocated by the spin - up of accretor due to accretion because in this case the critical mass can remarkable exceed @xmath153 . ]
thus , values of donor mass and orbital period differ only a little for the `` conservative '' track and for the `` non - conservative '' one .
the rate of angular momentum loss decreases abruptly and contraction of binary systems becomes slower after the donor msw is stopped . as a result donor contracts ceases to lose mass . since at the moment of cessation of mass loss
the radius of the donor is larger than the radius of thermally equilibrium ms star , the star contracts to the equilibrium radius . as a result
the ratio of stellar radius to the effective radius of roche lobe decreases even more . in this phase the evolution of binary system
is determined by gwr only . during this detached phase the binary
can not be identified as cv , and certain range of the orbital periods becomes devoid of stars .
this explains the origin of so called period gap of cvs .
observed period gap of cvs lies between @xmath154 and @xmath155 [ [ cat97 ] ] .
the boundaries of the `` theoretical '' period gap depend on adopted rate of angular momentum loss due to msw and on input physics parameters for evolutionary code that determine theoretical radii of stars .
we used a code that gives somewhat overestimated value for the lower boundary .
gap boundaries for the `` conservative '' track are equal to @xmath156 , and for the `` non - conservative '' one they are equal to @xmath157 .
this discrepancy is related to the different rate of angular momentum loss just prior to gap , and with different masses of accretors after termination of semi - detached phase of evolution as well .
note that the period gap is characterized not by total absence of cvs but rather by deficiency of cvs .
the presence of cvs in the period gap can be explained as follows : i ) for cvs with masses of donors equal to @xmath158 mass exchange begins immediately in gap ; ii ) in the case when the donor fills its roche lobe not on zams but later , when hydrogen is almost burn out , then as mass decreases down to @xmath32 complete mixing does not occur due to presence a helium core in the star ( see , e.g. , [ [ tfey87],[fe89 ] ] ) .
though , it must be admitted that the number of such systems is very small .
the phase of evolution after period gap is characterized by lower rate of mass loss by the donor , because evolution now is driven by gwr . in this phase
the rate of angular momentum loss from the system is sufficiently lower compared to the donor msw driven phase of evolution .
another feature of this phase of evolution is existence of a minimum period of binary .
short - period cutoff is related to the onset of significant degeneracy of the material of the donor and consequently to the change of mass - radius dependence law [ [ f71],[p81 ] ] : for small masses of donors @xmath159 their radii increase as mass decreases and the exponent of the mass - radius law tends to @xmath160 .
minimal period corresponds to the value of exponent equal to @xmath161 . ) for the effective radius of roche lobe @xmath162 , condition @xmath163 , third kepler s law @xmath164 , and definition of @xmath128 ( [ dr2r ] ) . as a result
we have : @xmath165 @xmath166 @xmath167 and , finally @xmath168 this conclusion is independent of the assumption of conservativity or non - conservativity of mass exchange . ] for `` conservative '' track the minimum period is equal to 75 minutes and for `` non - conservative '' one it is equal to 81 minutes . the difference in the values of @xmath169 is related to the difference in the total mass of the system and to the difference in the deviations of the radii of donor stars from thermal equilibrium values .
the latter deviations are caused by different rates of mass loss .
after passing the minimum the orbital period begins to increase , and in this phase the rate of mass loss by the donor quickly decreases .
it is seen from fig .
[ ris5 ] that difference between `` conservative '' and `` non - conservative '' tracks becomes larger in this stage .
this can be naturally explained by increase of the rate of angular momentum loss with decrease of @xmath54 in the `` non - conservative '' case .
`` conservative '' and `` non - conservative '' tracks have appreciably different values of mass transfer rates after passing the minimum of orbital period .
this fact is of importance since the probability of discovery of a cv in a given range of @xmath170 is proportional to    @xmath171    where @xmath172 is a positive constant @xmath173  for visual magnitude limited sample [ [ kkr98 ] ] .    after passing the minimum of the orbital period
( i.e. in period range @xmath174 ) the mean value of the probability of discovery of cv for the `` non - conservative '' track is 1.2 times larger than for the `` conservative '' one .
thus , despite a substantial difference in the rate of mass exchange , the probability of discovery of cv differs only a little in two cases under consideration .
this is due to faster crossing of the corresponding period range in the `` non - conservative '' case . nevertheless , predicted difference in @xmath55 for `` conservative '' and `` non - conservative '' cases could be reflected in the distribution of cvs over types of variability since the character of variability of a cv depends on the rate of mass exchange .
figure  [ ris5 ] also shows the evolutionary tracks for stars that would be dynamically unstable within the conservative model .
initial masses of the components were equal to @xmath144 and @xmath175 .
of course , this combination of masses is rather extreme but we adopted these parameters for the sake of illustration .
track 3 was calculated under the conventional assumptions on the loss of angular momentum due to gwr and msw only .
this is a typical case of unstable mass exchange with unlimited growth of mass exchange rate and fast convergence of the time scale of mass loss to the dynamical one , i.e. runaway . strictly speaking this run
should be considered as purely illustrative because in this case a common envelope has to form that increases the rate of angular momentum loss even more ( it was not taken into account in calculations ) .
track 4 was calculated within the model in which the binary loses matter and angular momentum in accordance with eq .
( [ xxxx ] ) .
it resembles the tracks for conventional case of stable evolution of cv and differs slightly from track 2 by a higher rate of mass loss in the initial phase of evolution .
the latter is caused by some increase of @xmath176 at large @xmath54 [ see eq .
( [ damsw ] ) ] .
it results in more strong deviations of the donor from thermal equilibrium and corresponding change of the upper border of period gap : it is equal to @xmath177 for track 4 , and @xmath178 for track 2 .
the minimum of period for track 4 is equal to 76 minutes , i.e. it is nearly the same as for track 1 but slightly less than for track 2 .
boundaries of stable mass exchange region were determined as follows : several sets of evolutionary tracks were calculated for systems with predefined masses of donors and values of @xmath54 that were varied with step 0.1 .
the value of @xmath54 was considered as critical if for a given evolutionary track the mass accretion rate did nt exceed eddington s limit and if for the next evolutionary track ( with @xmath54 higher by 0.1 ) it did .
note that consequent increasing of @xmath54 leads to unlimited growth of mass loss rate ( at values of @xmath54 greater by 0.2  0.3 than critical ) .
note also , that results of such procedure are in good agreement with technique based on comparison of derivatives of stellar radius and effective roche lobe radius .    the new boundary of the stable mass exchange region is shown in fig .
[ ris1new ] ( under assumption of `` non - conservative '' evolution ) .
it is well seen that the new region of stable mass exchange contains all observed cvs with estimated masses of components .
thus , within the `` non - conservative '' model the problem of location of observed systems outside the region of stable mass exchange does not arise .
on the other hand , in the @xmath65 diagram all observed cvs are concentrated near @xmath179 , while `` non - conservative '' model supposes sufficiently higher values of @xmath54 .
to investigate this phenomenon , we have calculated a set of `` non - conservative '' tracks with initial donor masses from @xmath180 up to @xmath181 with step @xmath138 and with such values of @xmath54 that tracks start in the vicinity of the new boundary of stable mass exchange region .
these tracks are shown in fig .
[ ris7 ] ( @xmath65 diagram ) . based on calculations of these tracks , it is possible to get isochrones that are also shown in fig .
[ ris7 ] for the moments of time equal to @xmath182 , @xmath183 , @xmath184 , @xmath185 , @xmath186 , @xmath187 years ( time is measured from the onset of mass exchange between components ) .
years and @xmath185 years for tracks with initial donor mass @xmath188 merge , since this interval of time binaries spend in the period gap . ]
it is this `` extreme '' set of tracks that was chosen since these systems spent the majority of time between an old boundary of the stable mass exchange region and the new one .
it is clear , that systems with smaller values of @xmath54 ( tracks of which begin to the left from new boundary ) spent there less time .
thus , cvs with parameters that are close to the `` extreme '' ones have most chances to be found out in this region .    with respect to the location of isochrones in fig .
[ ris7 ] all tracks can be subdivided into two distinct groups : i ) tracks for systems with massive donors ( from @xmath189 up to @xmath139 ) , for which the evolution prior to the period gap is defined mainly by the loss of systemic angular momentum due to the donor msw ; ii ) tracks for systems of low - mass donors ( from @xmath138 up to @xmath190 ) , for which this stage is absent , because their donors are initially completely convective . for the tracks of the first group a high ( @xmath191@xmath13 ) rate of mass loss by the donor prior
the period gap is typical , and , hence , donor mass and mass ratio decrease rapidly .
this results in a fast evolution in the @xmath65 diagram . as a result ,
the part of the diagram near the new boundary is crossed very quickly .
this explains the absence of observed cvs near new boundary .
for the tracks of the second group ( dropping a short initial phase with large rate of mass exchange ) the evolution of binary system is characterized by low and nearly constant mass loss rate ( @xmath192@xmath13 ) for a long time until period minimum .
consequently , the mass of the donor and mass ratio of components change slowly .
hence , if beginning the evolution near the new boundary of stable mass exchange region , the system should remain there long enough . the absence of cvs in this region can be explained as follows : masses of dwarfs in zero - age cvs can not be higher than @xmath193 that is minimum mass of helium core for a star leaving ms .
results of three - dimensional gas dynamical simulations of the flow structure in semi - detached binaries show that during mass exchange a significant fraction of the matter leaves the system .
on the other hand , a number of observed cvs have a combination of donor mass and mass ratio of components forbidden in the evolutionary models based on conventional assumptions on the loss of angular momentum of binaries , since in these cases mass exchange should be unstable .
this contradiction can be resolved in the framework of the model suggested here .
this model takes into account losses of mass and angular momentum from the system according to the results of gas dynamical simulations .
it is shown that approximation for outflowing angular momentum as the difference of specific angular momentum of the donor matter in @xmath0 and specific angular momentum of accreted matter ( being multiplied by appropriate mass flows ) permits to explain observations satisfactorily .
note , that the new `` non - conservative '' model practically does not change such parameters of evolutionary tracks as limits of period gap and minimal period .
thus , the main conclusion of our study is : if we use a `` non - conservative '' approximation for cv evolution , we obtain a new boundary of the region of stable mass exchange in @xmath194 plane which explains the distribution of observed cvs better than conventional models .
this work was supported by the russian foundation for basic research grants 99 - 02 - 17619 and 99 - 02 - 16037 and russian federation presidential grant 99 - 15 - 96022 .
1 .   [ rjw82 ] _
rappaport , s. , joss , p.c . , & webbink , r.f . _ // astrophys .
2 .   [ rvj83 ] _ rappaport , s. , verbunt , f. , & joss , p.c . _ // astrophys .
[ it84 ] _ iben , i.,jr . , &
tutukov , a.v . _ // astrophys .
[ p84 ] _ patterson , j. _ //
[ w85 ] _ webbink , r.f . _
// interacting binary stars / eds w.h.g .
lewin , j. van paradijs , e.p.j .
van den heuvel .
cambridge : cambridge univ . press , 1985 .
[ kolb93 ] _ kolb , u. _ // astron . astrophys .
p.149 7 .
[ ft94 ] _ fedorova , a.v . , &
tutukov , a.v . _
// astron . zhurn .
p.431 ( astron .
reports v.38 .
p.377 ) . 8 .
[ rit95 ] _ ritter , h. _ //
evolutionary processes in close binary stars / eds r.a.m.j .
wijers , m.b .
davies , c.a . tout .
dordrecht : kluwer , 1992 .
[ pol96 ] _ politano , m. _ //
p.338 . 10 .
[ fed97 ] _ fedorova , a.v . _
// binary stars / ed .
moscow : cosmoinform , 1997 .
p.179 ( in russian ) . 11 .
[ kmg62 ] _ kraft , r.p . ,
mathews , j. , greenstein , j.l . _ //
p.312 . 12 .
[ p67 ] _ paczyski , b. _ // acta astron .
. p.267 . 13 .
[ sch62 ] _ schatzman , e. _ // ann .
p.18 . 14 .
[ mes68 ] _ mestel , l. _ //
mon . not .
r. astron . soc .
p.359 . 15 .
[ egg76 ] _ eggleton , p.p . _ // structure and evolution of close binary stars / eds p.p .
eggleton , s. mitton , j. whelan .
dordrecht : reidel , 1976 .
p.209 . 16 .
[ vz81 ] _ verbunt , f. , & zwaan , c. _ // astron . astrophys .
[ yu73 ] _ yungelson , l.r . _
// nauchn . inform .
p.71 ( in russian ) . 18 .
[ kieboom81 ] _ kieboom , k.h . , &
verbunt , f. _ // astron . astrophys .
p.l11 . 19 .
[ kk95 ] _ king , a. , & kolb , u. _ // astrophys .
p.330 . 20 .
[ bbkc98b ] _ bisikalo , d.v . ,
boyarchuk , a.a . ,
chechetkin , v.m . ,
kuznetsov , o.a . , & molteni , d. _ // mon . not .
r. astron .
p.39 . 21 .
[ bbkc98c ] _ bisikalo , d.v . ,
boyarchuk , a.a . , kuznetsov , o.a .
, & chechetkin , v.m . _ //
astron . zhurn .
1998 . v.75 .
p.706 ( astron .
reports v.42 .
p.621 ; preprint astro - ph/9806013 ) . 22 .
[ eggl83 ] _ eggleton , p.p . _ // astrophys .
v.268 . p.368 .
[ pacyn71 ] _ paczyski , b. _ //
astron . astrophys .
p.183 . 24 .
[ kuz95 ] _ kuznetsov , o.a . _ // astron .
p.508 ( astron .
reports v.39 .
p.450 ) 25 .
[ ll62 ] _ landau , l.d .
, & lifshitz , e.m . _ the classical theory of fields .
1971 , oxford : pergamon press 26 .
[ skum72]_skumanich a. _ //
p.565 . 27 .
[ sr83 ] _ spruit , h.c . , & ritter , h. _ // astron .
1983 , v.124 .
p.267 . 28 .
[ hkn96 ] _ hachisu , i. , kato , m. , & nomoto , k. _ //
[ yl98 ] _ yungelson , l. , & livio , m. _ // astrophys .
j. 1998 , v.497 . p.168 . 30 .
[ kh94 ] _ kato , m. , & hachisu , i. _ // astrophys .
p.802 . 31 .
[ yuetal96 ] _ yungelson , l. , livio , m. , truran , j. , tutukov , a. , & fedorova , a. _ //
. p.890 . 32 .
[ lgr91 ] _ livio , m. , govarie , a. , & ritter , h. _
// astron . astrophys .
v.246 . p.84
[ skr98 ] _ schenker , k. , kolb , u. , & ritter , h. _ //
mon . not .
r. astron . soc .
p.633 . 34 .
[ ty71 ] _ tutukov , a.v . , & yungelson , l.r .
// nauchn . inform .
p.88 ( in russian ) . 35 .
[ tfy82 ] _ tutukov , a.v . ,
fedorova , a.v .
, & yungelson , l.r .
pisma astron . zhurn .
p.365 ( sov .
p.198 ) . 36 .
[ kool92 ] _ de kool , m. _ // astron .
v.261 . p.188
[ cat97 ] _ ritter , h. , & kolb , u. _ //
astrophys . suppl .
p.83 . 38 .
[ smith_dhillon ] _
smith , d.a . , & dhillon , v.s .
// mon . not .
r. astron . soc .
v.301 . p.767 .
[ sawada84 ] _ sawada , k. , hachisu , i. , & matsuda , t. _ // mon . not .
r. astron . soc .
[ matese83 ] _ matese , j.j . , &
whitmire , d.p .
// astrophys .
v.266 . p.776
[ bbkc98d ] _ bisikalo , d.v . ,
boyarchuk , a.a . ,
kuznetsov , o.a . , & chechetkin , v.m . _ //
astron . zhurn .
1999 . v.76 .
p.270 ( astron .
reports v.43 .
p.229 ; preprint astro - ph/9812484 ) . 42 .
[ rsh83 ] _ ruderman , m.a . , & shaham , j. _ // nature . 1983 . v.304 .
p.425 . 43 .
[ hp84 ] _ hut , p. , & paczyski , b. _ //
astrophys j. 1984 .
v.284 . p.675
[ fy84 ] _ fedorova , a.v .
, & yungelson , l.r . _
// astrophys .
space sci .
p.125 . 45 .
[ tfey85 ] _ tutukov , a.v . ,
fedorova , a.v .
, ergma , e.v . , & yungelson , l.r .
pisma astron . zhurn .
p.123 ( sov .
p.52 ) . 46 .
[ tfey87 ] _ tutukov , a.v .
, fedorova , a.v . ,
ergma , e.v . , & yungelson , l.r . _ //
pisma astron . zhurn .
p.780 ( sov
p.328 ) . 47 .
[ fe89 ] _ fedorova , a.v . , & ergma , e.v .
// astrophys .
space sci .
[ tf89 ] _ tutukov , a.v . , &
fedorova , a.v . _ //
p.1172 ( sov
p.606 ) . 49 .
[ hmma77 ] _ huebner , w.f . ,
merts , a.l . , magee , n.h . , & argo , m.f .
astrophysical opacity library : los alamos sci . lab . rep . la-6760-m .
los alamos , 1977 . 50 .
[ ajr83 ] _ alexander , d.r . ,
johnson , h.r . , & rypma , r.l . _
// astrophys .
v.272 . p.773
[ fgh77 ] _ fontaine , g. , graboske , h.c .
, & van horn , h.m .
// astrophys .
v.35 . p.293 . 52 .
[ den89 ] _ denisenkov , p.a . _ // nauchn .
p.145 ( in russian ) . 53 .
[ hfcz83 ] _ harris , m.j . , fowler , w.a . ,
caughlan , g.r .
, & zimmerman , b.a .
astron . astrophys .
p.165 . 54 .
[ cfhz85 ] _ caughlan , g.r . ,
fowler , w.a .
, harris , m.j . , & zimmerman , b.a . _ // atom .
data and nucl .
data tabl .
p.197 . 55 .
[ kr90 ] _ kolb , u. , & ritter , h. _ //
astron . astrophys .
v.236 . p.385
[ f71 ] _ faulkner , j. _ //
1971 v.170 . p.l99 .
[ p81 ] _ paczyski , b. _ // acta astron .
[ kkr98 ] _ kolb , u. , king , a.r . , & ritter , h. _ // mon . not .
r. astron .
v.298 . p.l29 .