there may be around @xmath1 stellar - mass black holes ( bhs ) in the galaxy ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , but only two dozen of them have been dynamically confirmed in x - ray binary ( xrb ) systems , in which a bh is accreting from its companion star ( * ? ? ?
* ; * ? ? ?
* for recent reviews ) .
a large fraction of the bh xrbs are close binaries with orbital periods @xmath2 day .
the spectral types of the companion stars range from a2v to m1v , indicating that most of them are less massive than @xmath3 , so these xrbs are correspondingly classified as low - mass x - ray binaries ( lmxbs ) .
all known bh lmxbs are transient systems , spending most of their time in quiescence with x - ray luminosities below @xmath4 , and occasionally exhibiting outbursts when the x - ray luminosities can rise up to @xmath5 .
galactic bh lmxbs are believed to evolve from primordial binaries consisting of a massive star and a low- or intermediate - mass companion ( * ? ? ?
* for a review ) .
the standard formation scenario involves a common envelope ( ce ) phase , during which the binary s orbital separation reduces drastically due to the loss of the orbital kinetic energy @xcite , since the orbital separations of the present day binaries are orders - of - magnitude smaller than the radii that the massive primaries can reach during their evolution .
when the primary star evolves and overflows its roche lobe , mass transfer proceeds on a dynamical timescale because the primary is much more massive than the secondary .
the secondary can not accrete all the mass transferred onto it , and an envelope forms enshrouding the two stars .
then the secondary star spirals into the envelope , and the orbital energy is used to expel the envelope ( for reviews of ce evolution , see * ? ? ?
* ; * ? ? ?
if surviving the ce event , the binary can continue to evolve ; otherwise the two stars would merge to be a single star .
though the ce stage is too short - lived to be detected , indirect observational evidence has emerged ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
bhs are thought to be the evolutionary outcomes of sufficiently massive stars ( @xmath6 ; e.g. , * ? ? ?
* ; * ? ? ?
. however , the envelope of such stars are generally too tightly bound for a low - mass secondary to expel @xcite , thus a merger is more likely to take place , unless the primary has lost a significant fraction of its envelope before the ce phase through stellar winds @xcite , or abnormally high values for the ce ejection efficiency are adopted @xcite .
an alternative scenario is that the secondaries in bh lmxbs are initially intermediate - mass stars , which can eject the primary s envelope more easily .
observational evidence for this scenario comes from the cno - processed material in xte j1118@xmath7480 from its ultraviolet spectra @xcite .
it is , however , difficult for an intermediate - mass x - ray binary ( imxb ) to evolve to short - period systems , unless some unusual angular momentum loss mechanisms are invoked ( e.g. , * ? ? ?
* ; * ? ? ?
meanwhile , these binaries reveal a discrepancy between the calculated effective temperatures and the observed spectral types of the donor stars .
the reason lies in that nuclear burning in the center of an intermediate - mass star causes a higher effective temperature than an ordinary low - mass star , though it has lost a significant fraction of its mass during the evolution @xcite .
clues to the initial mass range for the donors in bh lmxbs can be derived from the spin evolution of the bhs , since the bh spins increase with accretion @xcite .
two methods , i.e. , x - ray continuum - fitting @xcite and modeling relativistic reflection @xcite have been employed to estimate the spins of about 20 bhs ( * ? ? ?
* ; * ? ? ?
* and references therein ) .
assuming that bhs in lmxbs are born with negligible spins , @xcite showed that the initial masses of the secondaries in short - period binaries are generally @xmath8 .    from another point of view
, approaches have been taken in the study of ce evolution by refining the energy budget for the ejection as well as the core - envelope binding energy .
for example , @xcite argued that nuclear fusion should be added into the energy sources that may contribute to the envelope ejection .
@xcite proposed that a thermonuclear runaway which is caused by the mixing hydrogen from the secondary into the helium - burning shell of the primary would assist the ejection . when treating the binding energy of the envelope , @xcite pointed out that there might be mass outflows during the slow spiral - in stage , when the orbital energy is balanced with the enthalpy rather than merely the internal energy of the envelope .
that treatment can lower the binding energy by a factor of @xmath9 , and help allow for a low - mass companion to survive the ce phase .
it is also noted that the value of the binding energy is sensitive to the definition of the boundary between the remnant core and the ejected envelope ( e.g. , * ? ? ?
@xcite advocated that , the he core of a giant would experience a thermal readjustment stage after the ce event .
thus the boundary between the remnant core and the ejected envelope can be determined as the place where the hydrogen burning shell had maximal compression prior to ce evolution .
another issue in the formation of bh lmxbs is the mass range for the bh progenitors .
current supernova ( sn ) theories can not reliably predict which stars form bhs instead of neutron stars ( nss ) , although they are generally thought to be more massive than @xmath10 .
the bh masses in the 24 xrbs have been dynamically measured to be in the range of @xmath11 to @xmath12 ( * ? ? ? * and references therein ) .
however , there is statistical evidence for the presence of a dearth of nss or bhs with masses @xmath13 @xcite .
this is in contrast with the traditional thought that the distribution of bh masses should decay with mass ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , suggesting that the physics of sn explosions that lead to the formation of bhs is still unclear @xcite .
recently , @xcite argued that zero - age main sequence stars with mass @xmath14 might experience an unsuccessful explosion and evolve to bhs , while stars with mass around @xmath15 experience a successful sn explosion and become nss .
though this bifurcation is not well understood and may be related to the compactness parameter that describes the density profile outside the iron core @xcite , the failed sn mechanism can naturally account for the observed @xmath16 gap between nss and bhs .
@xcite proposed that this mechanism can interpret the red supergiant problem and the sn rate problem , which concerns the absence of red supergiants with masses @xmath17 as the progenitors of type iip sne , and the deficit of the observed cosmic sn rate compared to the observed cosmic star formation rate , respectively .    in this paper
we investigate the formation of galactic bh lmxbs incorporating the influence of the aforementioned factors with both stellar evolution and binary population synthesis ( bps ) methods . in section 2 ,
we calculate and compare the binding energy parameters in various models , which are then used in our bps study . in section 3
, we introduce the assumptions and prescriptions adopted in our bps calculations .
the results in various kinds of models are presented and compared with observations in section 4 .
we summarize the results and discuss their implications in section 5 .
in the standard ce model the change in the binary s orbital energy @xmath18 is responsible for the ejection of the ce , which is described by the following equations @xcite , @xmath19 and @xmath20 where @xmath21 is the binding energy of the envelope , @xmath22 the efficiency of converting the orbital energy to kinetic energy to eject the ce , @xmath23 the gravitational constant , @xmath24 the primary s mass and @xmath25 its core mass , @xmath26 the mass of the secondary , and @xmath27 and @xmath28 the orbital separations just before and after the ce phase , respectively .
when only gravitational binding energy is considered , @xmath29 is given by @xmath30 if the internal energy of the stellar matter ( including both thermal energy and recombination energy ) also contributes to the binding energy , then we have @xmath31 dm,\ ] ] where @xmath32 is internal energy @xcite .
recently , @xcite proposed a modified energy criterion taking into account the mass outflows during the spiral - in stage . in this so - called enthalpy model ,
the binding energy is expressed as @xmath33 dm,\ ] ] where @xmath34 and @xmath35 are the pressure and the density of the gas , respectively .
since the @xmath36 term is non - negative , the absolute value of @xmath21 becomes smaller in this case .    for convenience
, @xcite introduced a parameter @xmath37 characterizing the central concentration of the donor s envelope to describe the binding energy : @xmath38 where @xmath39 is the mass of the donor star s envelope , and @xmath40 is the ratio of the roche lobe radius and the orbital separation at the onset of ce .
generally , once a star fills its roche lobe , @xmath41 is taken to be the stellar radius .
combining eqs .
( [ equation3])([equation6 ] ) , one can obtain different @xmath37s from the gravitational binding energy , the total energy and the enthalpy prescriptions , namely @xmath42 , @xmath43 and @xmath44 , respectively .
we calculate their values using an updated version ( 7624 ) of the modules for experiments in stellar astrophysics ( mesa ) @xcite evolution code .
we consider stars with initial masses in the range of @xmath45 , which securely covers the range of the stellar - mass bh progenitors ( more massive stars produce negligible bhs according to the initial mass function and will evolve to wolf - rayet ( wr ) stars with almost all the hydrogen envelopes lost due to strong stellar winds ) .
the chemical composition of the stars is taken to be @xmath46 and @xmath47 .
the mixing length parameter @xmath48 is set to be 2.0 . for massive stars
, stellar winds may have non - ignorable influence on the @xmath37 parameter @xcite , but their loss rates have not been accurately determined . here , we adopt two prescriptions for the wind loss rate .
the first one is taken from @xcite and @xcite for o and b stars in different stages , denoted as wind1 hereafter , while the other one , denoted as wind2 , adopts the maximum value of the prescriptions above in all the evolutionary stages , to be consistent with @xcite , and to set an upper limit for the influence of the stellar wind .
the effect of stellar rotation is ignored , because in our calculation for lmxb formation , ce evolution usually occurs during case c mass transfer when the primary star has evolved to be a supergiant star with very slow rotation due to stellar evolution and/or tidal synchronization .
our calculations show that the evolution of @xmath44 , @xmath43 , and @xmath42 generally have the similar trend in most of the evolutionary stages , i.e. , decreasing with increasing stellar radius .
however , for stars less massive than @xmath49 , @xmath37s increase again when they have ascended the asymptotic giant branch and developed a deep envelope ( see also * ? ? ?
as an illustration , fig .
[ figure1 ] shows different @xmath37 values as a function of the stellar radius @xmath50 for a @xmath51 star , and fig .
[ figure2 ] shows the @xmath44 values as a function of the stellar radius for stars with different masses .
the solid and dashed lines correspond to the wind1 and wind2 prescriptions , respectively . in fig .
[ figure1 ] the red , blue , and green lines describe the evolution of @xmath44 , @xmath43 , and @xmath42 respectively in each case .
it is seen that the values of @xmath37 with the wind1 prescription are about three times as large as those with the wind2 prescription ( the binding energy with the wind1 prescription is half of the value with the wind2 prescription considering different envelope masses ) , indicating that the binding energy parameter for massive stars is sensitive to wind loss @xcite .
we plot the core mass against the stellar radius in fig .
[ figure3 ] , and list the stellar parameters along the evolutionary tracks in table  [ table1 ] in the two cases .
we also use an updated stellar evolution code called ev originally developed by @xcite to calculate the values of @xmath37s , and find that they are several times smaller than those with the mesa code in the wind2 case ( see also * ? ? ?
* ; * ? ? ?
the reason for this difference is that stars are generally more compact ( especially for stellar structure within the hydrogen burning shell ) during late evolution when modeled with the ev code .
we choose to use the results with the mesa code in the following because it adopts much denser grids for stellar structure than the ev code .
@xcite pointed out that a low - mass ( @xmath52 ) star is unable to eject the bh progenitor s envelope unless @xmath53 . in our calculation , this condition could be achieved for both @xmath44 and @xmath43 if the masses of the bh progenitors are @xmath54 .
binaries with more massive progenitor stars are unable to survive ce evolution even if @xmath44 is adopted ( see also fig .  [ figure2 ] ) .
we use an updated version of the bps code developed by @xcite to study the formation of bh lmxbs .
the code originates from the rapid single - star evolution code written by @xcite .
it involves various kinds of processes such as mass transfer , accretion , ce evolution , sn kicks , tidal interactions , and angular momentum loss .
we have made quite a few modifications of the bps code to follow the formation and evolution of xrbs @xcite .
we briefly describe some considerations adopted in this work .      given the binding energy parameter @xmath37 , we can estimate the post - ce separation with the following equation : @xmath55 if neither the he core nor the secondary star fills its roche lobe , the binary is regarded to survive the ce phase and can continue to evolve .
generally the value of @xmath22 is less than unity , and there is indication that the ce efficiency could be low @xcite .
moreover , it may vary with stellar structure and the binary parameters .
indeed , it has been suggested that @xmath22 may depend on the component mass and the orbital period in the studies of the formation of post - ce white dwarf binaries ( e.g. , * ? ? ?
* ; * ? ? ?
* ) . however , it is not known whether this can be applied to more massive binaries , so we still adopt constant values ( 1 and 0.2 ) for @xmath22 .
it is noted that @xmath22 and @xmath37 are usually combined to affect ce evolution .      besides ce evolution , the sn explosions play an important role in the formation of bh lmxbs .
the sn mechanisms are relevant to both the formation rate of xrbs and the mass distribution of the remnant bhs .
when determining the bh masses @xmath56 , one needs to account for the measured masses of known bhs and the @xmath57 gap between the ns and bh masses without any fine - tuning of stellar mass loss .
we consider both the prescriptions based on the neutrino - driven convection - enhanced sn mechanism , which is called the rapid sn mechanism , and the failed sn mechanism . in
the former case the gravitational remnant mass is described by @xcite , @xmath58 here @xmath59 is the mass of the proto - compact object and @xmath60 with @xmath61 where @xmath62 is the co core mass , @xmath63 , and @xmath64 .
the latter prescription is under the assumption that the bh formation is controlled by the compactness of the stellar core at the time of collapse : low compactness stars are more likely to explode as sne and produce nss , while high - compactness stars are more likely to evolve to failed sne that produce bhs @xcite .
in this case we assume @xmath65 or @xmath62 , where @xmath66 is the he core masses prior to core collapse ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the value of @xmath66 is determined as follows ( * ? ? ? * for a detailed description ) .
first , the core mass @xmath67 at the end of the hertzsprung - gap is set according to the initial mass of the star ( see eq .  28 of * ? ? ?
* ) . then the evolution of the core mass during the giant branch is determined by combining the power - law core mass - luminosity relation @xmath68 where @xmath69 , @xmath70 , and @xmath71 and @xmath72 depend on the stellar mass , and the energy conservation equation for hydrogen burning @xmath73 where @xmath74 is the envelope mass fraction of hydrogen and @xmath75 is the specific energy release .
thus @xmath76 where @xmath77 represents hydrogen rate constant .
a similar procedure is adopted to calculate the co core mass , with @xmath78 replaced by @xmath79 myr@xmath80 .    in all of the cases we require that @xmath81 .
we do not consider the formation of bhs by accretion - induced collapse of nss in xrbs .
we also assume that a natal  kick " is imparted on the newborn bhs , similar as nss .
the kick velocity is set to be inversely proportional to the remnant mass , i.e. , @xmath82 , where @xmath83 is the kick velocity for nss , which follows a maxwellian distribution with @xmath84 @xcite .
we consider the bh lmxbs evolved from primordial binaries with both incipient low - mass and intermediate - mass secondaries .
we assume a constant star formation rate of @xmath85yr@xmath80 @xcite , and evolve @xmath86 primordial binaries with the primary mass ( @xmath24 ) ranging from @xmath87 to @xmath88 following the @xcite initial mass function , and the secondary mass ( @xmath26 ) from @xmath89 to @xmath90 .
for the initial mass ratio @xmath91 , we consider both flat distribution @xmath92 with @xmath93 randomly distributed between 0 and 1 , and a  twin "- like distribution with @xmath94 .
the initial separation @xmath95 of the binary components varies between @xmath96 and @xmath97 , with a flat distribution for @xmath98 .
all binaries are assumed to be initially in circular orbits .
we stop the calculation when the donor mass in lmxbs drops down to @xmath89 , or the evolution time exceeds 15 gyr .
we summarize the models and the prescriptions employed in table .
[ table2 ] .
we first consider the formation of bhs via the rapid sn mechanism . in this case the progenitor stars are usually more massive than @xmath99 in binaries . for such massive stars
the binding energy of the envelope is so large that with @xmath100 a low - mass secondary inevitably merges with the he core of the primary , and no successful bh lmxbs can form ( see also * ? ? ?
so in models a1 and a2 we set @xmath101 with @xmath102 and @xmath44 , respectively .
the distributions of the donor mass of the resultant bh xrbs at current epoch are plotted in the left and right panels of fig .
[ figure4 ] , respectively . to reveal the donor mass distribution in different scale ,
the donor masses are confined within @xmath103 and @xmath104 in the upper and lower panels of the figure , respectively .
the solid and dashed red lines represent the results with the @xmath37s obtained with the wind1 and wind2 prescriptions , respectively . also displayed is the distribution of the observed bh lmxbs with the black line ( data taken from * ? ? ?
* ; * ? ? ?
the number of the bh xrbs is calculated by multiplying the formation rate with the duration in a specific stage .
we only display the sources which are subject to the thermal - viscous instability in an x - ray irradiated accretion disk and  assumed " to be transient @xcite .
they include almost all bh lmxbs produced .
the exceptions are the systems comprised of intermediate - mass donors or with extremely short periods , which are not relevant here . in model a1 and a2 , the donors of persistent sources locate in the mass range of @xmath105 , and their total number is comparable with that of transient sources . when @xmath106 , most bh xrbs harbor an intermediate - mass donor star ( with mass @xmath107 ) , implying that lower - mass secondaries are unlikely to survive ce evolution . the corresponding formation rate of the bh xrbs is too low ( @xmath108 yr@xmath80 ) to account for the observed systems in the galaxy . when @xmath109 , the formation rate enhances a bit ( @xmath110 @xmath111 ) , and the donor mass distribution peaks at around @xmath112 .
it is seen that in both cases the distributions significantly deviate from the observed one .
the calculated orbital period distributions are too long ( @xmath113 1 day ) to be consistent with observations ( see fig .
[ figure12 ] below ) .
these results clearly disfavors the rapid sn mechanism , so in the following we only consider the failed sn mechanism .    in models
b1 and b2 , we assume @xmath114 and a flat mass ratio distribution with @xmath115 and @xmath116 , respectively .
we present the calculated distributions of the donor mass for the two models and compare them with observations in fig .
[ figure5 ] ( left : b1 ; right : b2 ) . here
the red and blue lines correspond to @xmath101 and 0.2 , respectively , and the line styles are same as in fig .
[ figure4 ] . with the failed sn mechanism the masses of the bh progenitors can be lower than @xmath117 , so it is easier for low - mass secondaries to survive the ce phase . from the left panel of fig .
[ figure5 ] we see that in the cases of wind1/@xmath101 , wind2/@xmath101 , and wind1/@xmath118 the donor mass peaks at @xmath119 , matching the observational data , while in the case of wind2/@xmath118 , the donor masses cluster around @xmath120 , larger than the observed peak .
this is simply because more massive donor stars are required to survive the ce evolution when the product @xmath121 is lower .
however , larger @xmath122 does not always necessarily lead to larger number of bh lmxbs , since the formation of a bh xrb strongly depends on the orbital separation just before and after the ce phase . with the same initial orbital periods ,
larger @xmath122 results in wider post - ce orbits , which may cause the subsequent roche lobe overflow to be postponed , reducing the x - ray lifetime and thus the number of xrbs .
figure  [ figure6 ] shows the birthrates of the bh xrbs as a function of the initial donor mass in each case .
we mention that , although the bh xrb birthrates can be as high as @xmath123 yr@xmath80 , the birthrates of those with incipient low - mass donors are several times lower .
we use the kolmogorov - smirnov ( ks ) test to quantitatively compare the the univariate distributions of the calculated and measured donor mass .
the ks statistic @xmath72 quantifies a distance between the empirical distribution functions of the two samples .
the null hypothesis that the two distributions are drawn from the same one is rejected at a significance level @xmath124 if @xmath125 @xcite . in table
[ table2 ] we list the calculated values of @xmath72 and @xmath126 for all models ( except models a1 and a2 because of too few data points ) in the case of wind2/@xmath101 at @xmath127 and 0.05 .
figure  [ figure7 ] shows the donor mass distributions for models c1 ( left panel ) and c2 ( right panel ) under the assumptions that @xmath128 and @xmath109 . in models c1 and c2 we adopt a flat distribution of the initial mass ratio and a  twin " distribution ( * ? ? ?
* ; * ? ? ?
* see however , cantrell & dougan 2014 ) , respectively . here
we confine the donor masses within the range of @xmath129 because of the small number of xrbs with more massive donor stars .
figure  [ figure7 ] shows that the peak of the donor mass distribution in the cases of wind1/@xmath101 , wind2/@xmath101 , and wind1/@xmath130 is @xmath131 , but shifts towards @xmath132 in the case of wind2/@xmath118 , roughly in accordance with fig .
[ figure5 ] .
this is also reflected by their similar @xmath72 values . to see whether the resultant bh lmxbs evolve from incipient low- or intermediate - mass binaries , we repeat the calculation with the initial donor mass confined to be @xmath133 and @xmath134 respectively . the results in model c1 ( with the wind2 prescription )
are plotted in fig .
[ figure8 ] .
it is clearly seen that the total number of the bh lmxbs evolved from the primordial binaries with the initial donor mass @xmath135 is almost the same as in the left panel of fig .
[ figure7 ] , so we can safely conclude that almost all the bh lmxbs in our calculation have incipient low - mass secondaries .
although the bh xrbs with low - mass donors have a lower birthrate than those with intermediate - mass donors , it is compensated by their long x - ray lifetimes .
figure  [ figure7 ] also indicates that assuming a @xmath94 distribution ( model c2 ) leads to less efficient formation of bh lmxbs than with a flat @xmath93 distribution because of less progenitor binaries with low - mass secondary stars .
figure  [ figure9 ] shows the distribution of the orbital period vs. the donor mass ( left panel ) , and the distribution of the mass transfer rate vs. the donor mass ( right panel ) in model c1 .
the color scale represents the relative number of the xrbs in a parameter space .
it can be seen that the majority of bh xrbs are short - period ( @xmath136 day ) systems with mass transfer rates below @xmath137 .
the above results are based on the hypothesis that the bh mass is equal to its progenitor s he core mass at the time of core bounce . to check the effect of another choice of the remnant mass ,
we construct models d1 and d2 , which are same as models c1 and c2 except with @xmath138 , respectively .
the donor mass distributions are shown in fig .
[ figure10 ] .
we find that the formation rate in mode d1 decreases to about half that in model c1 , and the peak of the donor mass distribution extends to a broader range compared with fig .
[ figure6 ] . the remnant mass influences the formation of the bh lmxbs by determining the roche lobe size after the sn explosion .
figure  [ figure11 ] shows the distributions of the bh masses under the two choices .
the red and blue lines represent the results in models c1 and d1 , which peak at @xmath139 and @xmath140 , respectively . @xcite
have analyzed the data of 16 bh x - ray transients , and found that the observed bh masses can be best described by a narrow distribution at @xmath141 , which is obviously more compatible with the result in model c1 .      in figs .
[ figure12 ] we show the calculated distributions of the orbital period and the effective temperature of the donor stars in model a. the color scale represents the relative number of the xrbs .
we also plot the measured orbital periods and effective temperatures of the donors for galactic bh lmxbs . in fig .
[ figure12 ] for models a1 ( left panel ) and a2 ( right panel ) , the observed distribution obviously deviates from the calculation , showing a tendency towards lower temperature .
the reason is that , as stated before , most of the bh lmxbs in models a1 and a2 descended from imxbs , so the donor star is likely to have a higher temperature than an originally low - mass star .
we get similar results for model b2 with @xmath142 . in other models the calculated orbital periods are more compatible with observations .
we also use the ks test to quantitatively compare the calculated and observed distributions , and the results for models b1 , b2 , c1 , and d1 are shown in table  [ table3 ] .
we can see that , except model b2 with @xmath118 , other models are all acceptable .
we show the distributions in the best - fit model c1 in fig .
[ figure13 ] .
previous studies have already demonstrated that it is hard to produce compact bhxrbs consisting of an initially low - mass companion with the standard formation scenario , since the companion star is unable to expel the ce ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* see however , yungelson & lasota 2008 ; wiktorowicz et al . 2013 ) .
however , the orbital periods , the donors spectral types , and the measured spins of quite a few bhs strongly suggest that in the incipient binaries the donor stars are likely to be of low - mass .
we construct a series of models to explore the plausible solutions to the above - mentioned puzzle .
the factors we have taken into account include the binding energy parameter @xmath37 , the ce efficiency @xmath22 , the sn mechanisms for bh formation , and the distribution of the mass ratio . by adopting different choices of these prescriptions ,
we examine their roles in the formation of bh lmxbs by comparing the calculated parameter distributions with observations .
our calculations show that the @xmath37-values strongly depend on the stellar mass and the evolutionary stage .
generally @xmath44 is several times larger than @xmath42 and @xmath43 , and this can help the binary to avoid merging during the ce evolution without invoking unphysically high values of @xmath22 .
however , we fail to reproduce the observational characteristics of galactic bh lmxbs based on the rapid sn mechanism even with @xmath44 adopted , since the bh progenitors are still too massive ( @xmath143 ) in this scenario .
we then consider the alternative failed sn mechanism , and assume that the remnant bh mass equals the he or co core mass of the primary .
this time we obtain more satisfactory results , with the distributions of the donor mass , temperature and orbital period for the current binaries largely consistent with observations .
in particular , most bh lmxbs in this case originate from primordial binaries with incipient low - mass ( @xmath133 ) secondary stars , compatible with the requirement derived from the spin estimates of several bhs in lmxbs . finally , all models predict that most bh lmxbs are transient sources due to thermal instability in the accretion disks .
thus , in the framework of the failed sn mechanism , it is possible to reproduce the galactic bh lmxbs with normal value of @xmath144 .
realizing the uncertainties in the treatment of ce evolution , we also check the results by adopting a small @xmath118 in each case .
our results demonstrate that , when @xmath122 is large enough for a low - mass star to survive the ce phase , the donor mass distribution is not sensitive to the choice of @xmath37 or @xmath22 , which affects the number and birthrate of the produced bh lmxbs .
however , for very small @xmath122 , the donor mass distribution strongly depends on the value of @xmath37 or @xmath22 , because surviving the ce phase becomes paramount and difficult .
in addition , the prescription of the bh mass with the he core mass leads to better match of the observations than with the co core mass .
our results provide useful constraints on the formation process of bh lmxbs , but are subject to many uncertainties in both observational and theoretical aspects .
when comparing the calculated results with observations , one needs to be cautious that this comparison is influenced by small number statistics as well as selection effects .
currently there are around 20 dynamically confirmed stellar - mass bhs . it is not known whether they are correctly representative of the overall population in the galaxy .
most bh lmxbs are transient sources , and we are lack of thorough understanding of their outburst characteristics ( including the peak luminosity and its relation with other parameters , and the duty cycles ) and the physical mechanisms @xcite . to make things more complicated , a new class of very faint x - ray transients were recently discovered with peak luminosities of only @xmath145 ergs@xmath80 @xcite . further more ,
considering the selection effect that luminous sources are more likely to be observed , we may expect that a large fraction of the lmxbs possess even lower - mass donor stars than the current sample , because they are relatively dim in x - rays even during outbursts .
there are also big uncertainties in modeling the evolutionary sequences of compact star binaries .
we take into account @xmath146 in calculating the binding energy of the envelope ( see also * ? ? ?
* ) , but its credibility is controversial and there are arguments both for and against it in the literature , let alone the definition of the core - envelope boundary from which the binding energy is integrated ( see * ? ? ? * for a discussion ) .
it is also highly uncertain whether the ce efficiency @xmath22 is constant or depends on other physical parameters . @xcite and @xcite suggested that it may decrease with increasing mass ratio , but this effect has not been confirmed in modeling the sdss post - ce binary population ( e.g. , * ? ? ? * ) and needs further verification .
even if true , the physics behind it is not clear , and numerical simulations in 3d of ce evolution are eagerly required to resolve this issue . as to the sn and the bh formation mechanisms , our results base on the suggestion that stars of mass @xmath147 die in failed sne leaving a bh @xcite .
the model still needs to be tested by both 3d numerical simulations of the collapse of massive stars and targeted field surveys for sne in nearby galaxies .
obviously a thorough investigation on the evolutionary history of bh lmxbs needs to incorporate bps with the aforementioned input physics in a self - consistent way .
we are grateful to an anonymous referee for helpful comments .
this work was supported by the natural science foundation of china under grant numbers 11133001 and 11333004 , and the strategic priority research program of cas ( under grant number xdb09000000 ) .
99 natexlab#1#1 bailyn , c. d. , jain , r. j. , coppi , p. , & orosz , j. a. 1998 , , 499 , 367 belczynski , k. , wiktorowicz , g. , fryer , c. l. , holz , d. e. , & kalogera , v. 2012 , , 757 , 91 brown , g. e. , & bethe , h. a. 1994 , , 423 , 659 cantrell , a. g. , & dougan , t. j. 2014 , , 445 , 2028 casares , j. & jonker , p. g. 2014 , spa .
, 183 , 223 chen , w .- c .
, & li , x .- d .
2006 , , 373 , 305 clausen , d. , piro , a. l. , & ott , c. d. 2015 , , 799 , 190 coriat , m. , fender , r. p. , & dubus , g. 2012 , , 424 , 1991 davis , s. w. , blaes , o. m. , hubeny , i. , & turner , n. j. 2005 , , 621 , 372 davis , p. j. , kolb , u. , & knigge , c. 2012 , , 419 , 287 de kool , m. 1990 , , 358 , 189 de marco , o. , passy , j .- c . , moe , m. , herwig , f. , mac low m .-
m . , paxton , b. 2011 , , 411 , 2277 drake , j. j. , & sarna , m. j. 2003 , , 594 , 55 eggleton , p.
, , 151 , 351 eggleton , p.
, , 156 , 361 fabian , a. c. , parker , m. l. , wilkins , d. r. , et al .
2014 , , 439 , 2307 fabian , a. c. , rees , m. j. , stella , l. , & white , n. e. 1989 , , 238 , 729 farr , w. m. , sravan , n. , cantrell , a. et al .
2011 , , 741 , 103 feigelson , e.  d. , & jogesh , b. g. 2012 , modern statistical methods for astronomy , cambridge university press ( cambridge , uk ) fragos , t. , & mcclintock , j. e. 2015 , , 800 , 17 fryer , c. l. 1999 , , 522 , 413 fryer , c. l. , belczynski , k. , wiktorowicz , g. et al .
2012 , , 749 , 91 fryer , c. l. , & kalogera , v. 2001 , , 554 , 548 han , z. , podsiadlowski , p. , & eggleton , p.
, , 270 , 121 han , z. , podsiadlowski , p. , maxted , p. f. l. , marsh , t. r. , & ivanova , n. 2002 , , 336 , 449 han , z. , podsiadlowski , p. , maxted , p. f. l. , & marsh , t. r. 2003 , , 341 , 669 haswell , c. a. , hynes , r. i. , king , a. r. , & schenker , k. 2002 , , 332 , 928 hobbs , g. , lorimer , d. r. , lyne , a. g. , & kramer , m. 2005 , , 360 , 974 horiuchi , s. , nakamura , k. , takiwaki , t. , kotake , k. , tanaka , m. 2014 , , 445 , l99 hurley , j. r. , pols , o. r. , & tout , c. a. 2000 , , 315 , 543 hurley , j. r. , tout , c. a. , & pols , o. r. 2002 , , 329 , 897 ivanova , n. 2002 , phd thesis , oxford university ivanova , n. 2006 , , 653 , l37 ivanova , n. 2011 , , 730,76 ivanova , n. , & chaichenets , s. 2011 , , 731 , l36 ivanova , n. , justham , s. , chen , x. , et al .
2013 , a&ar , 21 , 59 justham , s. , rappaport , s. , & podsiadlowski , p. 2006 , , 366 , 1415 kalogera , v. 1999 , , 521 , 723 kiel , p. d. , & hurley , j. r. 2006 , , 369 , 1152 kiziltan , b. , kottas , a. , de yoreo , m. , & thorsett , s. e. 2013 , , 778 , 66 kobulnicky , h. a. , & fryer , c. l. 2006 , , 670 , 747 kochanek , c. s. 2014 , , 785 , 28 kochanek , c. s. 2015 , , 446 , 1213 kreidberg , l. , bailyn , c. d. , farr , w. m. , & kalogera , v. 2012 , , 757 , 36 kroupa , p. , tout , c. a. , & gilmore , g. 1993 , , 262 , 545 lasota , j .- p .
2001 , , 45 , 449 lasota , j .-
p . , dubus , g. , & kruk , k. 2008 , , 486 , 523 li , x .- d .
2015 , , 64 - 66 , 1 mcclintock , j. e. , narayan , r. , & steiner , j. f. 2014 , space sci .
, 183 , 295 oconnor , e. , & ott , c. d. 2011 , , 730 , 70 ohlmann , s. t. , rpke , f. k. , pakmor , r. , & springel , v. 2016 , , 816 , l9 zel , f. , psaltis , d. , narayan , r. , & mcclintock , j. e. 2010 , , 725 , 1918 zel , f. , psaltis , d. , narayan , r. , & santos villarreal , a. 2012 , , 757 , 55 paczynski , b. 1976 , in iau symposium , vol . 73 , structure and evolution of close binary systems , ed .
p. eggleton , s. mitton , & j. whelan ( dordrecht , d. reidel publishing co. ) , p. 75
paxton , b. , bildsten , l. , dotter , a. , et al .
2011 , , 192 , 3 paxton , b. , cantiello , m. , arras , p. , et al .
2013 , , 208 , 4 paxton , b. , marchant , p. , schwab , j. , et al .
2015 , , 220 , 15 pinsonneault , m. h. , & stanek , k. z. 2006 , , 639 , l67 podsiadlowski , p. , rappaport , s. , & pfahl , e.d .
2002 , , 565 , 1107 podsiadlowski , p. , rappaport , s. , & han , z. 2003 , , 341 , 385 podsiadlowski , p. , ivanova , n. , justham , s. , & rappaport , s. 2010 , , 406 , 840 portegies zwart , s. f. , verbunt , f. , & ergma , e. 1997 , , 321 , 207 remillard , r. a. , & mcclintock , j. e. 2006 , ar , 44 , 49 shao , y. , & li , x .- d .
2014 , , 796 , 37 shiode , j. h. , & quataert , e. 2014 , , 780 , 96 sion , e. m. , bond , h. e. , lindler , d. , et al .
2012 , , 751 , 66 smith , l.  f. , biermann , p. , & mezger , p.  g.  1978 , , 66 , 65 smith , n. , & arnett , w. d. 2014 , , 785 , 82 smith , n. , li , w. , filippenko , a. v. , & chornock , r. 2011 , , 412 , 1522 sukhbold , t. & woosley , s. e. 2014 , , 783 , 10 taam , r. e. & sandquist , e. l. 2000 , ar , 38 , 113 tauris , t. m. , & dewi , j. d. m. 2001 , , 369 , 170 thorne , k. s. 1974 , , 191 , 507 timmes , f. x. , woosley , s. e. , weaver , t. a. 1996 , , 457 , 834 van den heuvel , e. p. j. 1983 , in accretion - driven stellar x - ray sources , ed . w. h. g. lewin & e. p. j. van den heuvel ( cambridge and new york , cambridge university press ) , p. 303
vink , j. s. , de koter , a. , & lamers , h. j. g. l. m. 2001 , , 369 , 574 webbink , r. f. 1984 , , 277 , 355 wijnands , r. , int zand , j. j. m. , rupen , m. , et al .
2006 , , 449 , 1117 wiktorowicz , g. , belczynski , k. , & maccarone , t. j. 2013 , arxiv:1312.5924v2 wong , t .- w . , valsecchi , f. , ansari , a. , et al .
2014 . , 790 , 119 woosley , s. e. , & weaver , t. a. 1995 , , 101 , 181 xu , x .- j .
, & li , x .- d .
2010a , , 722 , 1985 xu , x .-
j . , & li , x .- d .
2010b , , 716 , 114 yungelson , l. r. , & lasota , j .-
2008 , , 488 , 257 zhang , s. n. , cui , w. , & chen , w. 1997 , , 482 , l155 zorotovic , m. , schreiber , m. r. , gsicke , b. t. et al .
2011 , , 536 , l3     with the stellar radius @xmath50 for a @xmath148 star with initial mass of @xmath149 .
the the red , blue and green lines represent @xmath44 , @xmath43 , and @xmath42 , and the solid and dashed lines represent the results with the wind1 and wind2 prescriptions , respectively.[figure1],scaledwidth=100.0% ]     values as a function of the stellar radius for stars with different mass .
the solid and dashed lines represent the results with the wind1 and wind2 prescriptions , respectively .
[ figure2],scaledwidth=100.0% ]     star . here
, the boundary of the core and the envelope is determined to be the maximal compression point in the hydrogen burning shell .
the solid and dashed lines represent the results with the wind1 and wind2 prescriptions , respectively .
[ figure3],scaledwidth=100.0% ]                                                                              .the stellar parameters for a @xmath51 star at different evolutionary stages .
columns 2 - 4 and 5 - 7 correspond to the wind1 ( left ) and wind2 ( right ) prescriptions , respectively [ table1 ] [ cols="^,^,^,^,^,^,^",options="header " , ]