arxiv: v2 [astro-ph.sr] 24 Feb 2013

Transcription

1 Mon. Not. R. Astron. Soc.,?? () Printed 6 February 3 (MN LATEX style file v.) On the origin of W UMa type contact binaries - a new method for computation of initial masses arxiv:3.635v [astro-ph.sr] 4 Feb 3 M. Yıldız and T. Doğan Department of Astronomy and Space Sciences, Science Faculty, Ege University, 35, Bornova, İzmir, Turkey. INTRODUCTION The evolution and structure of a star are primarily determined by its initial mass. Therefore, the structure of stars in contact binaries (CBs) can be very different from single stars because of mass transfer. The component stars of CBs evolve as single stars with initial mass (M i) until the binary interaction occurs. When the mass transfer starts, the initially more massive component (M i, the current secondary) loses mass that is transferred to the initially less massive component (M i, the current primary) and the rest leaves the system. Because of the interaction, the mass ratio of the system reverses. After mass transfer stops or sufficiently slows down, the secondary stars are found to be overluminous for a main-sequence (MS) star of the same mass. The reason for this luminosity excess is that the luminosity of such stars depends on their initial mass as much as their present mass. In this study, we develop a method to find initial masses of the secondary stars of W UMa type binaries in terms of the luminosity excess. W UMa type CBs are among the best studied stars in the literature. For about CBs (5 A-subtypes and 49 W-subtypes), accurate dimensions of component stars from photometric and spectroscopic data are successfully derived by different investigators. The subtypes are defined by their light-curves (Binnendijk 97). mutlu.yildiz@ege.edu.tr ABSTRACT W UMa type binaries have two defining characteristics. These are (i) the effective temperatures of both components are very similar, and (ii) the secondary (currently less massive) component is overluminous for its current mass. We consider the latter to be an indication of its mass before the mass-transfer event. For these stars, we define a mass difference (δm) between the mass determined from its luminosity and the present mass determined from fitting the binary orbit. We compare the observed values of the mass difference to stellar models with mass-loss. The range of initial secondary masses that we find for observed W UMa type binaries is.3.6m. We discover that the A- and the W-subtype contact binaries have different ranges of initial secondary masses. Binary systems with an initial mass higher than.8 ±.M become A-subtype while systems with initial masses lower than this become W-subtype. Only 6 per cent of systems violate this behavior. We also obtain the initial masses of the primaries using the following constraint for the reciprocal of the initial mass ratio: < /q i <. The range of initial masses we find for the primaries is..5m, except for two systems. Finally in comparing our models to observed systems, we find evidence that the mass transfer process is not conservative. We find that only 34 per cent of the mass from the secondary is transferred to the primary. The remainder is lost from the system. Key words: binaries: close-binaries: eclipsing-stars: evolution-stars: interior-stars: late-type. We compute initial masses of secondary and primary components of these systems. Our sample includes well-determined parameters of CBs from the compilations of Yakut & Eggleton (5), Gazeas et al. (5, 6) and Zola et al. (5, ). Beside the luminosity excess of the secondary components in the W UMa type binaries, another very important feature of these systems is that the effective temperatures, T e, of components are almost identical. However, the slight difference between effective temperatures of components is systematic between the two subtypes. For the A-subtype systems, the effective temperature of the massive component (T e) is higher than the effective temperature of the lower mass component (T e); the opposite is true for the W-subtype systems. There is debate in the literature concerning whether the subtypes represent an evolutionary sequence (Awadalla & Hanna 5; Gazeas & Niarchos 6; Eker et al. 6). According to Maceroni & van t Veer (996), however, there is no evolutionary link between most of the A-subtype systems and the W-subtype CBs. In this paper, we aim to confirm that the subtypes are determined by their initial parameters. If initial parameters of A- and W-subtypes are completely different, then they are not an evolutionary sequence. There are two evolutionary pathways that lead to a star filling its Roche lobe. These are nuclear and angular momentum evolution (Hilditch, King & McFarlane 989). It is thought that angular momentum evolution via magnetic braking (Schatzman 96; Mestel 968) is the primary mechanism for the formation and evolution of

2 M. Yıldız and T. Doğan W UMa type CBs (Okamoto & Sato 97; van t Veer 979; Vilhu 98). However, both these two mechanisms affect the kind of W UMa binaries that are formed. Also because nuclear and angular momentum evolution are very strongly dependent on stellar mass, one may expect that the initial masses of components of W UMa type binaries are the key parameters for our understanding of the formation and evolution of W UMa type binaries. Each component of a CB fills its Roche lobe and is under the gravitational interaction of the other. Apart from this, the more massive primary component is an MS star and the luminosity of the secondary component is higher than the expected MS luminosity (Yang & Liu ; Webbink 3; Li et al. 8). This excess may be due to large amount of energy transfer that occurs (see Section ) from the primary to the secondary (Lucy 968a,b; Webbink 3; Li et al. 8). However, the evolutionary status of CBs is still unclear. According to Hilditch, King & McFarlane (988), the primary components of the W-subtype CBs are not evolved MS stars, and the A-subtype CBs are close to the terminal-age of the main sequence (TAMS). Although there are many studies in the literature about contact systems, the contact models are still not able to satisfy all the observational constraints. It is still an open question whether CBs are in thermal equilibrium or not. It is generally accepted that a contact configuration in thermal equilibrium is not possible and therefore CBs suffer from thermal relaxation oscillation (Flannery 976; Qian ; Kähler a,b, 4; Kreiner et al. 3; Webbink 3; Li, Han & Zhang 4; Yakut & Eggleton 5; Paczynski et al. 6). However, according to some authors, a contact configuration can be in thermal equilibrium. They argue that contact systems can achieve thermal equilibrium if the secondary, currently less massive component, is more evolved than the massive component (Stȩpień 6a,b, 9, ). In such a model, CB stars originate from detached close binaries with initial orbital periods close to a few days and the system was exposed to mass exchange with a mass ratio reversal in the past. According to these assumptions, energy exchange does not influence stellar radii and occurs in thermal equilibrium (Kähler 4; Stȩpień ). Large-scale energy transfer from the primary to the secondary components solves several outstanding problems in CBs (Lucy 968a). However, the problem of energy transfer requires an intrinsic driving mechanism and its effect on the evolution of the component stars and the system is unknown. These processes are very poorly understood because of the complexity of detailed hydrodynamic treatment of the flow and also due to the unknown driving mechanism. This paper is organized as follows. We outline the physical properties of CBs in Section and we discuss the basic properties of CBs in the mass-luminosity (M L) and the mass-radius (M R) diagrams. In Section 3, we develop and apply a comprehensive method for computation of initial masses of CBs based on stellar modelling with mass-loss. In Section 4, we apply this method to the CBs and discuss the results. Finally, in Section 5, we draw our conclusions. THE PHYSICAL PROPERTIES OF CONTACT BINARIES AND THE PROBLEM OF ENERGY TRANSFER The most reliable accurate parameters (mass, radius and luminosity) of binary stars are computed by using both photometric and spectroscopic data. We compiled such well known CBs from the literature. Their basic properties are listed in Table A of the Appendix. These properties give clues to evolutionary phases of the primary components, at least. The primary components of CBs in M L and M R diagrams are much closer to that of the normal MS stars than the secondary components. Their luminosities and radii are plotted with respect to their masses in Figs (a) and (b), respectively. For comparison, the components of the welldetermined detached eclipsing binaries (+; hereafter DEBs) are also included (Torres, Andersen & Giménez ). The solid line in these figures shows zero-age MS (ZAMS), derived from models constructed with the MESA code (Paxton et al. ). There is an agreement between the luminosities of the primaries of CBs and that of the components of DEBs. If we compare the radii of the primary stars of CBs in Fig. (b), we can confirm that these two groups of stars are in very good agreement with the DEB stars. Fig. (b) also shows the TAMS line (as a dashed line). Most of the primaries of CBs occur between the ZAMS and the TAMS lines. TheM L andm R diagrams of the secondary components of W UMa type CBs are plotted in Figs (c) and (d), respectively. Unlike the primaries shown in Figs (a) and Figs (b), there is no agreement between the secondary components and the components of DEBs inm L andm R diagrams. Therefore, the secondary components of W UMa type CBs must have a very different structure or evolutionary path. All the secondary components are above the ZAMS line in both Figs (c) and (d). The mean masses of the secondary components of the A- and W-subtype CBs are.38 and.47m, respectively. Despite these low masses, the secondary stars are very bright and very large in comparison to their MS counterparts of the same mass. This excess in luminosity and radius may be explicable by the early evolution of the systems. In the literature, the close T es of components of CBs have been explained as a result of an energy transfer from the primary to the secondary star (Smith 984). According to this approach, the luminosity excess of the secondary components (δl = L M 4 ) must compensate the lack of luminosity of the primary components (δl = L M 4 ). If this is the case, there must be a correlation betweenδl andδl. In Fig.,δL is plotted with respect toδl. It should be noted that δl is always greater than zero and δl is in general (four fifths of systems) less than zero. If the luminosity excess of the secondary stars compensates the lack of luminosity of the primary stars, the position of the CB components should be on the δl = δl line (solid line). There is a large scatter for the A-subtype CBs (circles) and therefore no evidence for the energy transfer process is shown for this subtype. For the W-subtype CBs (filled circles), however, some of the systems with small values of δl and δl form a sequence on the line δl = δl. The sequence is very clear from the origin to the point with δl = L and δl = L. The details of this range are also shown in Fig.. The energy transfer process does not work, at least, for the A-subtype CBs. For the W-subtype CBs, however, there are some systems (about ) in which the transfer process may play a part for very close effective temperatures. For some binaries, however, both components have excess luminosity which the energy transfer process cannot explain. 3 METHOD FOR COMPUTATION OF INITIAL MASS FROM LUMINOSITY EXCESS As one of the essential parameters of a star, luminosity depends on the nuclear burning in the core. For a single star, this is deter-

3 On the origin of W UMa type contact binaries - a new method for computation of initial masses 3.5 A-subtype W-subtype DEB Primaries.6 Primaries.4 log L.5 log R (a) -. (b) log M Secondaries log M Secondaries.6.4 log L.5 log R (c) log M -.4 (d) log M Figure. Mass-luminosity and mass-radius relation for the primaries of the A- (circle) and the W- subtype (filled circle) CBs are plotted in panels (a) and (b), respectively. For comparison, the components of the well-known DEBs (+) are also plotted. Solar units are used. A solid line indicates the ZAMS line and a dashed line the TAMS line, obtained by using the MESA stellar evolution code. Panels (c) and (d) are same as (a) and (b), respectively, but for secondary components of W UMa type CBs. mined by the initial mass of the star. For a component star which has experienced mass transfer (and loss), M i = M f + M, () where M is the mass lost or transferred from the secondary, i.e. the mass difference between the initial (M i) and the final (M f ) masses. The physical conditions in the central regions of the secondary star, which yield luminosity higher than M 4 f, are similar to that of a star with a mass lower than M i but higher than M f. More precisely, the central physical conditions of the secondary star are very close to a single star of massm L L.5 (in solar units). For M L, which is the mass of an isolated star with the same luminosity, one can use a more realistic mass-luminosity relation. The M L relation we employ is L =.49M 4.6, () derived from the grids obtained by Yıldız () for TAMS models. The value of mass according to luminosity strongly depends on the evolutionary phase of a mass losing star at the beginning of the mass loss/transfer process. For the computations relating to the secondary stars of CBs, for now, we assume that mass transfer starts near the TAMS phase of the old massive component. Therefore, equation () is used throughout this paper [see Section 5.3 for correction due to deviation from equation ()]. For a typical TAMS star, the quantity M L defined by M L = (L /.49) /4.6 for the secondary star computed from equation () is equal to the observed M = M f. If there is a luminosity excess, then M L > M. We define the mass difference between these two masses as δm: δm = M L M. (3) According to our assumption, the luminosity excess of the secondary component of a CB is related to its initial mass as much as to its present mass, and therefore to the amount of total mass lost by the secondary ( M). This means that there must be a correlation between M and δm: M = f(δm). Then the initial mass of the secondary component of a CB is given as M i = M +f(δm) (4) We note thatm i is a function only of the observable quantitiesl andm. The masses, luminosities and radii are in solar units.

4 4 M. Yıldız and T. Doğan δl Table. Model results with mass-loss to obtain an expression for the difference between initial and final masses. The final mass is chosen as.5m and the mass transfer begins when central hydrogen abundance X C = 6. Model M i M f L f M L δm M No (M ) (M ) (L ) (M ) (M ) (M ) δl Figure. Excess in luminosity of the secondary component (δl = L M 4 ) is plotted with respect to luminosity lack (or excess) of the primary components (δl = L M 4 ). If the energy transfer is the dominant operating mechanism, all the CBs should make a sequence on the δl = δl line (solid line). This is the case for only some of W-subtype CBs. If one findsm i, thenm i can be computed by using the constraint on the reciprocal of initial mass ratio /q i: < /q i <. M i is a function of M = M i M and mass lost by the system (M lost ). The increase in the primary mass can be written as M M i = M M lost = M( γ). (5) where γ is the ratio M lost / M. γ = corresponds to the conservative case and γ = is the case in which all the mass lost by secondary component is lost by the system. The reciprocal of the initial mass ratio as a function of current masses, M i, and γ is given as /q i = Mi M (Mi M)( γ) =. (6) M i M i The unknown value of γ can be found from the constraint < /q i < (see section 4.). Our next task is to find an expression for f(δm). 3. Method derived from the MESA MODELS 3.. Initial masses of the secondary components As discussed above, the secondary components of CBs have a completely different interior from the normal MS stars because their initial mass is significantly higher than the present mass. This situation forces the secondary components to have a different evolutionary path from single MS stars without mass-loss. Change in mass essentially causes a breakdown in hydrostatic equilibrium, which is then restored by expansion of the central regions. We construct models in order to assess the ultimate effect of change in mass on the structure of stars by using the MESA stellar evolution code. From these models, we aim to find an expression for M as a function ofδm. Models with constant and variable masses are constructed by using the MESA code. The details of this code can be found in Paxton et al. (). The chemical composition and convective parameter for all the models are taken as X=.7, Z=., α=.89. Stellar models are constructed with the solar mixture given by Grevesse & Sauval (998). The radiative opacity is derived from Table. Same as Table, but for the final mass M f =.8M. Model M i M f L f M L δm M No (M ) (M ) (L ) (M ) (M ) (M ) recent OPAL tables (Iglesias & Rogers 996), implemented by the low-temperature tables of Ferguson et al. (5). The thermonuclear reactions rates are computed by using Angulo et al. (999) and Caughlan & Fowler (988), and the standard mixing length theory is employed for convection (Cox & Giuli 968). The MESA ρ T tables are based on the 5 update of the OPAL EOS tables (Rogers & Nayfonov ). The progenitors of the secondary stars of CBs should be evolved stars (Hilditch et al. 988), otherwise an excess in luminosity cannot be observed. Therefore, when we construct our stellar models we assume that mass transfer begins near the TAMS which we take to be when the central hydrogen abundance is 6. The basic properties of the constructed models with different initial masses are listed in Table. The range of initial masses is. 3.5M. For all of these models, the mass loss rate is 8 M yr and the final mass is.5m. After several tries for different mass-loss rates, we deduce that the results are independent of the mass-loss rate. The corresponding mass to the luminosity (M L) after mass transfer is computed from equation (). The mass difference δm is as given in equation (3) and M = M i M f. We also construct models withm f =.8M. The basic properties of these models are given in Table. The initial mass is plotted with respect to M L in Fig. 3 for M f =.5M and.8m. The relation between the initial and luminosity masses is almost linear. In both sets of models with either M f =.5M or.8m, we note that the higher the M, the higher is the δm. In Fig. 4, M is plotted with respect to δm (in units ofm ). Both sets of models withm f =.5M (filled circle) andm f =.8M (circle) have the same relation between M and δm. The solid line shows the fitted line[.5(δm.7).64 ] for both data sets. It is remarkable to confirm that the relation is not a linear one. Also shown in Fig. 4 is the dotted line for the linear relationship M =.δm. The initial mass of the secondary stars of CBs is then com-

5 On the origin of W UMa type contact binaries - a new method for computation of initial masses M f =.5M Sun M f =.8M Sun A-subtype W-subtype M i M,i M L M Figure 3. The relation between M L and M i for the final masses.5 (filled circle) and.8 (circle) M. M i -M f M f =.5M Sun M f =.8M Sun M L -M f Figure 4. M as a function of δm for the MESA models with final masses M f =.5 (filled circle) and M f =.8M (circle). The solid line represents the fitting curve M =.5(δM.7).64 and dotted line is for M =.δm. puted by using M i = M + M = M +.5(δM.7).64 (7) where the masses are in units of solar mass. 3.. Uncertainties in the initial masses of the secondary components The initial masses of the secondary components are computed from L and M. Therefore, the uncertainty inm i ( M i) arises from uncertainties M and L. Define f = M =.5(δM.7).64. Then, f is a function of M L and M. In terms of these quantities and their uncertainties, the uncertainty in M i can be written as M i = M + f dm L L + M L dl f M M. (8) Figure 5. Initial mass of the secondary components of CBs. Circles are for A-subtype and the filled circles are for W-subtype. The solid black line represents the derived transition mass with the dashed lines showing theσ error bars. Apart from a few systems, the A-subtype CBs have M i higher than.8 M and the W-subtype CBs have M i lower than.8m. Using the definition of δm given in equation (3), we find that f = f M L M =.59(δM.7).36. (9) From the employed M Lrelation, we obtain dm L =.. () dl L.76 The explicit form of M i as a function of M and L is (.59. L M i = M + (δm.7).36 4 RESULTS AND DISCUSSION L The initial masses of the secondary components + M ). () The initial mass of the secondary component of a CB is computed using the method from the previous section. The results are listed in Table A. If we compare an A-subtype with a W-subtype CB of the samem (present secondary mass), then we find that the initial mass of the secondary of an A-subtype is.5m higher than that of W-subtype. This means that the (present) secondary components of the A-subtype have lost/transferred much more mass than those of the W-subtype. In Fig. 5, the initial mass of the secondaries, derived by using the method based on MESA models, is plotted with respect to the present mass of the secondary components. The initial masses are plotted with their uncertainties computed by a method presented in the previous section. It can clearly be seen that subtypes are determined by the initial mass of the secondaries. If M i >.8 ±.M, the system becomes A-subtype. If M i <.8 ±.M, then the system exhibits the W-subtype properties. There are just a few systems that do not obey this prescription (see below). The mean initial mass of the secondaries is M ia =.97M for the A- and M iw =.56M for the W-subtype CBs. The underlying mechanism for this separation must be based on competition between the orbital and the nuclear time-scales.

6 6 M. Yıldız and T. Doğan There are only two W-subtype systems (Wa) having an initial mass higher than M ia. They are V78 Her and V4 Aur. These binaries are the hottest of the W-subtype CBs. Their M L relation and angular momentum properties are similar to that of the A-subtype rather than the W-subtype. The number of the A-subtype binaries (Aw) having an initial mass significantly less than M iw, however, is four. These binaries are EQ Tau, OO Aql, V58 Oph and TZ Boo. They have T e about 56 K [log(t e) = 3.75] and are the coolest of the A-subtype CBs. The O Connell effect may be responsible for this atypical situation if it is not due to uncertainties in the observed T es. However, in addition to the initial mass of the secondary components, another secondary effect may also play a role in the formation of W UMa subtypes. For example, initial orbital angular momentum may be important. Putting the atypical CBs aside, the initial mass range for the secondary components is.7.6m for A-subtype and.7.9m for W-subtype. However, for most of the CBs the range is.3.m. CC Com is the system with the minimum secondary initial mass,.86 M. The nuclear time-scale of such a star is about 6 Gyr, greater than the age of the galaxy. Therefore, this system (or most of the W-subtype CBs) must have some specific feature, which makes its fundamental properties unusual. As a matter of fact, both components of CC Com are brighter than ZAMS stars and fainter than TAMS stars. Formation of binaries with M i <.5M seems to be different from the other CBs. They have either a different value of γ or different structure for the progenitors of the secondary components during the mass transfer process. For example, the progenitors may not be as evolved (not TAMS) as we assume (see Section 5.3). The most realistic uncertainties in M i are computed from equation () and given in Table A. For some of the binaries, for example, HV UMa, the uncertainty M i is very high due to high uncertainties in M and L. However, this result does not change the situation much because the uncertainties are within an acceptable range for most of the binaries in our sample. 4. Initial masses of the primary components The reciprocal mass ratio/q i is the physical constraint for the determination of initial masses of the primary components. The fitting parameter γ (equation 5) is the mass lost by the system in units of M. First, we consider A-subtype CBs. If γ has a small value, then some A-subtype systems have an unreasonably negative value for M i. In order to keep M i > (/q i > ) for all A-subtype CBs, the minimum value of γ is.5. The maximum value of γ for /q i <, however, is.87. These two values of γ give mean value for /q i as.368 and.666, respectively. For /q i =.5, γ =.664. This value ofγ yields M i ranging from.3 to.5m for A-subtypes. This means that the mass lost by these subtypes is higher than the transferred mass: 66 per cent of mass lost by the secondary components is lost by the system and 34 per cent is transferred to the primary components. /q i of A-subtype CBs with γ =.664 is plotted with respect to M in Fig. 6. The reciprocal of the initial mass ratio is well within the expected range. For/q i of W-subtype CBs, we use the same value of γ. On average, the A-subtype CBs have higher /q i than the W-subtype CBs. There are only 3 violating systems. These binaries are the W-subtype CBs (see below) and have very lowδm values, less than.m. They may not be TAMS stars but are somehow in between the ZAMS and the TAMS lines. The value of γ we found is a mean value. Indeed, there is no /q i A-subtype W-subtype M Figure 6./q i computed from the expressions derived from the models by using MESA code is plotted with respect to M. All the A-subtype systems are in the range < /q i < ; just three W-subtype systems are greater than one. restriction for a single value of γ for all CBs. However, what is important is that there is a solution with a single value of γ for almost all CBs, and this is certainly true for the A-subtype CBs. It may be concluded that the value ofγ does not dramatically change much from system to system. The initial masses of the primary components are listed in Table A. The range of primary initial mass of the CBs is..5m, apart from two binaries (HV UMa and V376 And). This is the mass range in which magnetized stellar wind is an efficient mechanism for angular momentum loss (Tutukov, Dremova & Svechnikov 4). The mean value of the initial primary masses is M i =.99M for the A- and M i =.84M for the W-subtype. The uncertainties inm i are also given in Table A. M i is taken as M i/q i. 4.3 Correction of the initial masses of components of CBs withδm <.35M As previously stated, the initial mass of some secondaries is very small for some CBs, for example CC Com. In Fig. 7, the reciprocal of initial mass ratio /q i is plotted with respect to δm both for the A- and the W-subtype CBs. For A-subtype, /q i does not show any dependence onδm. For W-subtype, however,/q i is dependent on δm if δm <.4M. The solid line shows the curve fit( 3.5δM.5 +.8) for this range. This relation reflects our assumption about the evolutionary phase of secondary components at the beginning of mass transfer. It seems that for some of the W-subtype binaries the secondary components are not close to the TAMS line. Therefore, equation (7) is not valid for these CBs. Using the curve ( 3.5δM.5 +.8) of Fig. 7, we obtain a new equation for M i, and M i is computed as in Section 5.. This equation removes the slope in Fig. 7. The corrected initial masses of the secondary and the primary components of these CBs are listed in Table 3. The upper limit for the initial mass of the secondaries is.9m, consistent with the results given above. The minimum initial mass of the secondaries is with smallδm pertaining to BH Cas,.36M.

7 On the origin of W UMa type contact binaries - a new method for computation of initial masses A-subtype W-subtype Table 4. Cluster member CBs. Star M i M i t MS (Gyr) t cl (Gyr) P(day) Subtype Cluster /q i..8 TX Cnc W Praesepe AH Vir W Wolf 63 QX And A NGC 75 AH Cnc A M67.6 Table 5. Mass intervals for final products Figure 7./q i computed from the expressions derived from the models by using the MESA code is plotted with respect toδm. /q i for the A-subtype CBs has a very uniform distribution. For the W-subtype CBs, however, this is not the case. For small values ofδm,/q i is dependent onδm. The solid line shows the fitted curve for the W-subtype CBs ( 3.5δM.5 +.8). Table 3. The corrected masses for the secondary and the primary masses of the W-subtype CBs with δm <.35M. The correction is computed according to the expression obtained from the curve in Fig. 7 ( 3.5δM.5 +.8). Star δm M M M i M i /q i t MS (Gyr) DN Cam EF Boo ET Leo AA UMa YY Eri ER Ori BB Peg AE Phe W UMa AM Leo V78 Tau RZ Com AO Cam GZ And TW Cet BX Peg AB And SW Lac VW Cep V757 Cen CC Com XY Leo V53 Cas BH Cas RW Dor Some notes on evolutionary time-scales and cluster member CBs As mentioned above, the ultimate fate of a binary depends on its total mass and how this mass is shared between its components. If mass of the massive component (M i) is high enough, then nuclear evolution is very fast compared to orbital evolution. Therefore, the minm i (M ) maxm i (M ) M i (M ) t min (Gyr) Product.3±..8± W-subtype.8±..6± A-subtype lifetime of this star is so short that the system cannot be a CB. In this respect, a very interesting result we obtain is that there is an upper limit for the mass of the massive components of CBs,.6M (V9 Her). This result leads us to consider binaries with components having a mass higher than.6m as candidate systems for binaries with compact objects, for example. Ages of the W UMa type CBs are investigated by many different groups using different methods (Guinan & Bradstreet 988; van t Veer 996; Bilir at al. 5; Li, Han & Zhang 5; Li et al. 7) and range from.. to 7.Gyr. The life-time of a W UMa type CB can roughly be divided into three different phases. These are detached, semi-detached and contact phases. From the initial masses, at least for the A-subtype CBs, we find how long the detached phase is: t det = t MS. The MS age for a star with.6m is.46gyr. This time is likely to be very short for angular momentum evolution without mass-transfer. The MS age of a star with transition mass(.8m ) from the A- to the W-subtype CBs is.37gyr (EX Leo). For such systems, the nuclear and angular momentum time-scales are comparable. For some W-subtype systems, the nuclear time scale is liable to be so long that the components evolve into contact phase before significant nuclear evolution of the massive component. The MS lifetime for BH Cas, for example is 3.3Gyr. Such a time is very long for angular momentum evolution of BH Cas and similar CBs with two late-type components. According to the results, we deduce that ages of the A-subtype binaries (t s) are greater than the MS lifetime of the progenitor of their secondaries (t MS). For the W-subtype binaries, however, t s may also be less than t MS because it seems that the mass transfer for some CBs starts before the secondary component(m i) reaches the TAMS line. In order to test this point, we consider CBs in clusters. There are four such binaries in our CB sample. Two of them are the W-subtype CBs. The W-subtype binaries TX Cnc and AH Vir are members of Praesepe and Wolf 63, respectively. The A- subtype CBs are QX And of NGC 75 and AH Cnc of M67. The MS lifetime according to M i (t MS) of the binaries and age of clusters are listed in Table 4, respectively. While t MS is less than t s = t clus (cluster age) for A-subtype, t MS for the W-subtype binaries is less than t s for AH Vir and greater than for TX Cnc. Although there are only four well-known cluster member systems, these results are consistent with our expectations. Table 5 summarizes the results obtained for the initial masses of the A- and the W-subtype CBs. The MS lifetime according to M i is given in the fourth column. The age of the A-subtypes must be greater than this age. For the W-subtype CBs, however, age can be either greater or smaller than the given value.4gyr, depending

8 8 M. Yıldız and T. Doğan on the efficiency of the angular momentum loss mechanism via magnetic stellar wind. ACKNOWLEDGMENTS The anonymous referee and Peter P. Eggleton are acknowledged for their suggestions which improved presentation of the manuscript. 5 CONCLUSION In the present study we have investigated the physical properties of CBs. The approximate initial masses of the secondary and primary components are computed, for the first time ever, by using the basic parameters of W UMa type CBs (5 A- and 49 W-subtype). The initial mass of the secondary is computed from its luminosity excess by a method based on the expression derived from stellar models with mass-loss. We have also deduced that only one-third of the mass lost by the secondary components is transferred to the primary components. The remaining part leaves the systems. The range of initial masses of the secondaries of CBs is.3.6m. The upper part of this range is occupied by the A-subtype CBs and the W-subtype CBs are in the lower part. The transition occurs at about.8±.m. The main assumption in our method is that the mass transfer starts near or after the TAMS phase of the massive component (progenitor of the secondary component). This assumption seems fairly good for the A-subtypes. For the W-subtype, however, a correction is needed. From δm and the reciprocal mass-ratio diagram (Fig. 7), we develop a method and obtain corrected initial masses for the case of small δm. The corrected initial masses of the secondaries range from.3.9m. One of the important outcomes of the present study is that there is also a lower limit for the initial mass of the primary components. The reason for this result should be studied in the context of angular momentum evolution versus nuclear evolution. For the binary systems with an initial mass higher than.8m, although the angular momentum is lost by the less massive component the nuclear evolution is principally the mechanism responsible for the Roche lobe filling process. Therefore, semidetached systems withm i =.8.6M must yield A-subtype W UMa CBs. The binary systems withm i less than.8m evolve into the contact phase due to the rapid angular momentum evolution. The Roche lobes of their components contract so that they form the W-subtype contact or near-contact binaries despite the components not being much evolved in the MS phase. The primary components of W UMa type CBs have an initial mass range from. to.5m. A star with mass in this range loses its spin angular momentum relatively rapidly. If it has a tidally interacting companion, orbital angular momentum is lost by the system. The precursors of W UMa type CBs are the systems in which nuclear evolution of the massive components is in competition with the angular momentum loss from the less massive components. In the precursors of the W-subtype W UMa CBs, both components are effective in the angular momentum loss process. The method we develop for computation of the initial masses of the secondary components is highly suitable for a complete error analysis. From the derivatives of the function derived from the fitting curve for M as a function of L and M (equation 8), we successfully compute M i in terms of the observed uncertainties L and M (equation ). We also consider cluster member W UMa type CBs. Although the number of CBs with accurate dimensions available in the literature is very small (4 systems), the results support our findings on initial masses and time-scales for the A- and W-subtype CBs. REFERENCES Ahn, Y.S., Hill, G., Khalesseh, B., 99, A&A, 65, 597 Angulo, C., Arnould, M., Rayet, M., Descouvemont, P., Baye, D., Leclercq-Willain, C., Coc, A., et al., 999, Nucl. Phys. A, 656, 3 Awadalla, N.S., Yamasaki, A., 984, Ap&SS, 7, 347 Awadalla, N.S., Hanna, M.A., 5, J.Korean Ast. Soc., 38, 43 Baran, A., Zola, S., Rucinski, S.M., Kreiner, J.M., Siwak, M., Drozdz, M., 4, AcA, 54, 95 Barone, F., di Fiore, L., Milano, L., Russo, G., 993, ApJ, 47, 37 Bilir, S., Karatas, Y., Demircan, O., Eker, Z., 5, MNRAS, 357, 497 Binnendijk, L., 97, Vistas Astron.,, 7 Caughlan, G.R., Fowler, W.A., 988, At. Data Nucl. Data Tables, 4, 83 Chochol, D., van Houten, C.J., Pribulla, T., Grygar, J.,, Cont. Astron. Obs. Skalante Pleso, 3, 5 Cox, J.P., Giuli, R.T., 968, Principles of Stellar Structure (New York: Gordon & Breach) Csák, B., Kiss, L.L., Vinkó, J., Alfaro, E.J.,, A&A, 356, 63 Eker, Z., Demircan, O., Bilir, S., Karataş, Y., 6, MNRAS, 373, 483 Ferguson, J.W., Alexander, D.R., Allard, F., Barman, T., Bodnarik, J.G., Hauschildt, P.H., Heffner-Wong, A., Tamanai, A., 5, ApJ, 63, 585 Flannery, B.P., 976, ApJ,5, 7 Gazeas, K.D., Baran, A., Niarchos, P., Zola, S., Kreiner, J.M., Ogloza, W., Rucinski, S.M., Zakrzewski, B., Siwak, M., Pigulski, A., Dróżdż, M., 5, AcA, 55, 3 Gazeas, K.D., Niarchos, P., Zola, S., Kreiner, J.M., Rucinski, S.M., 6, AcA, 56, 7 Gazeas, K.D., Niarchos, P., 6, MNRAS, 37, L9 Goecking, K.-D., Duerbeck, H.W., 993, A&A, 78, 463 Goecking, K.-D., Duerbeck, H.W., Plewa, T., Kaluzny, J., Schertl, D., Weigelt, G., Flin, P., 994, A&A, 89, 87 Grevesse, N., Sauval, A.J., 998, Space Sci. Rev., 85, 6 Guinan, E.F., Bradstreet, D.H., 988, Formation and Evolution of Low Mass Stars, Dordrecht: Kluwer, p Gürol, B., Terzioğlu, Z., Gürsoytrak, S.H., Gökay, G., Derman, E., a, AN, 33, 69 Gürol, B., Derman, E., Saguner, T., Gürsoytrak, S.H., Terzioğlu, Z., Gökay, G., Demircan, Y., Okan, A., Saral, G., b, NewA, 6, 4 Hilditch, R.W., King, D.J., 986, MNRAS, 3, 58 Hilditch, R.W., King, D.J., McFarlane, T.M., 988, MNRAS, 3, 34 Hilditch, R.W., King, D.J., McFarlane, T.M., 989, MNRAS, 37, 447 Hrivank, B.J., 988, ApJ, 335, 39 Hrivank, B.J., Guinan E.F., DeWarf, L.E., Ribas, I.,, AJ,, 84 Iglesias, C.A., Rogers, F.J., 996, ApJ, 464, 943 Kallrath, J., Milone, E.F., Breinhorst, R.A., Wilson, R.E., Schnell, A., Purgathofer, A., 6, A&A, 45, 959 Kaluzny, J., Rucinski, S.M., 986, AJ, 9, 666

9 On the origin of W UMa type contact binaries - a new method for computation of initial masses 9 Kähler H., a, A&A, 395, 899 Kähler H., b, A&A, 395, 97 Kähler H., 4, A&A, 44, 37 Köse,., Kalomeni, B., Keskin, V., Ulaş, B., Yakut, K.,, AN, 33, 66 Kreiner J., Rucinski, S.M., Zola, S., Niarchos, P., Ogloza, W., Stachowski, G., Baran, A., Gazeas, K., Drozdz, M., Zakrzewski, B., Pokrzywka, B., Kjurkchieva, D., Marchev, D., 3, A&A, 4, 465 Li, L., Han, Z., Zhang, F., 4, MNRAS, 35, 37 Li, L., Han, Z., Zhang, F., 5, MNRAS, 36, 7 Li, L., Zhang, F., Han, Z., Jiang, D., 7, ApJ, 66, 596 Li, L., Zhang, F., Han, Z., Jiang, D., Jiang, T., 8, MNRAS, 387, 97 Lucy, L.B., 968a, ApJ, 5, 3 Lucy, L.B., 968b, ApJ, 53, 877 Maceroni, C., van t Veer, F., 996, A&A, 3, 53 Mestel, L., 968, MNRAS, 38, 359 Metcalfe, T.S., 999, AJ, 7, 53 Milone, E.F., Wilson, R.E., Hrivnak, B.J., 987, ApJ, 39, 35 Niarchos, P.G., Manimanis, V.N., 3, A&A, 45, 63 Okamoto, I., Sato, K., 97, PASJ,, 37 Özkardeş, B., Erdem, A., Bakış, V., 9, NewA, 4, 46 Özkardeş, B., Erdem, A.,, NewA, 7, 43 Paczynski, B., Szczygiel, D.M., Pilecki, B, Pojmanski, G., 6, MNRAS, 368, 3 Paxton, B., Bildsten, L., Dotter, A., Herwig, F., Lesaffre, P., Timmes, F.,, ApJS, 9, 3 Pribulla, T., Chochol, D., Rovithis-Livaniou, H., Rovithis, P., 999, A&A, 345, 37 Pribulla, T., Vaňko, M.,, Cont. Astron. Obs. Skalante Pleso, 3, 79 Qian, S.,, MNRAS, 38, 94 Qian, S.-B., Yang, Y.-G., 4, AJ, 8, 43 Rainger, P.P., Hilditch, R.W., Bell, S.A., 99, MNRAS, 46, 4 Ren, A.B., Zhang, X.B., Luo, C.Q., Luo, Y.P., Deng, L.C., Luo, Z.Q.,, NewA, 6, 94 Rogers, F.J., Nayfonov, A.,, ApJ, 576, 64 Samec, R.G., Hube, D.P., 99, AJ,, 7 Schatzman, E., 96, Annales dastrophysique, 5, 8 Smith, R.C., 984, QJRAS, 5, 45 Stȩpień, 6a, AcA, 56, 99 Stȩpień, 6b, AcA, 56, 347 Stȩpień, 9, MNRAS, 397, 857 Stȩpień,, AcA, 6, 39 Szalai, T., Kiss, L.L., Mészáros, Sz., Vinkó, J., Csizmadia, Sz., 7, A&A, 465, 943 Torres, G., Andersen, J., Giménez, A.,, A&ARv, 8, 67 Tutukov, A.V., Dremova, G.N., Svechnikov, M.A., 4, Ast.Rep., 48, 9 Ulaş, B., Kalomeni, B., Keskin, V., Köse, O., Yakut, K.,, NewA, 7, 46 Vaňko, M., Pribulla, T., Chochol, D., Parimucha, Š., Kim, C.H., Lee, J.W., Han, J.Y.,, Cont. Astron. Obs. Skalante Pleso, 3, 9 van t Veer, L., 979, A&A, 8, 87 van t Veer, F., 996, in Milone E. F., Mermilliod J.-C., eds, ASP Conf. Ser., The Origins, Evolutions and Destinies of Binary Stars., Astron. Soc. Pac., San Francisco, 9, 8 Vilhu, O., 98, A&A, 9, 7 Webbink, R.F., 3,in Turcotte, S., Keller, S.C., Cavallo, R.M., eds, ASP Conf. Ser., 93, 76, 3D Stellar Evolution Yakut, K., Kalomeni, B., İbanoğlu, C., 4, A&A, 47, 75 Yakut, K., Eggleton, P., 5, ApJ, 69, 55 Yakut, K., Ulaş, B., Kalomeni, B., Gülmen, Ö, 5, MNRAS, 363, 7 Yakut, K., Zima, W., Kalomeni, B., Van Winckel, H., Waelkens, C., De Cat, P., Bauwens, E., Vučković, M., Saesen, S., Le Guillou, L., Parmaksızoğlu, M., Uluç, K., Khamitov, I., Raskin, G., Aerts, C., 9, A&A, 53, 65 Yang, Y., Liu, Q.,, AJ,, 45 Yang, Y., Liu, Q., 3a, AJ, 6, 96 Yang, Y., Liu, Q., 3b, NewA, 8, 465 Yıldız, M.,, submitted Zhang, X.B., Zhang, R.X., 4, MNRAS, 347, 37 Zola, S., Rucinski, S.M., Baran, A., Ogloza, W., Pych, W., Kreiner, J.M., Stachowski, G., Gazeas, K., Niarchos, P., Siwak, M., 4, AcA, 54, 99 Zola S., Kreiner J.M., Zakrzewski B., Kjurkchieva D.P., Marchev D.V., Baran A., Rucinski S.M., Ogloza W., Siwak M., Koziel D., Drozdz M., Pokrzywka B., 5, AcA, 55, 389 Zola S., Gazeas K.D., Kreiner J.M., Ogloza W., Siwak M., Koziel- Wierzbowska D., Winiarski, M.,, MNRAS, 48, 464

10 M. Yıldız and T. Doğan APPENDIX A: BASIC PROPERTIES OF W UMA TYPE BINARIES Table A: Fundamental properties of the A- and W-subtype CBs. Columns are organized as name, subtype, period, mass ratio, effective temperatures of primary and secondary, primary and secondary radii, luminosities of primary and secondary, primary and secondary masses, initial masses of primary and secondary before the mass transfer (see Section 5), reciprocal of initial mass ratio, MS life time of secondary component according to initial mass and numbers for references. Star t P T T R R L L M M M,i M,i /q i t MS Ref d (K) (K) (R ) (R ) (L ) (L ) (M ) (M ) (M ) (M ) Gy HV UMa A.7 7± 7±.6±.5.8± ±.73.8±.6.5±.7.54±.69.± ,37 V376 And A ± 7335±.66±..55±. 3.4± ±.4.49±.6.76±.3.5±.3.9± RR Cen A ± ± 3.4±..7±.5.±.65.9±.35.9±.43.45±..±.38.5± V9 Her A ± 7346±.75±..4±. 3.5±.9 5.9±.4.7±.5.5±.3.6±.9.37± UZ Leo A ± 5 683± 5.9±..39±..±.8 3.7±.7.99±.4.6±..35±.7.4± V535 Ara A.69 8± 5 89± 5.9±.3.3±. 8.± 3. 6.±..94±.4.59±..64±.5.5± AQ Tuc A ±..3± ± ±.38.93±..69±.8.±.3.4± ,3,37 EF Dra A V388 Oph A.8 69± ± 3.6±..3±. 3.5±..43±.7.8±..34±..35±.4.± AW UMa A ±.5.66± ±.4.4±..±.3.3± ,37 V776 Cas A.44 67± 9 675± 9.8±..75±. 5.9±..±.6.75±.4.4±..5±.8.4± V59 Per A.76 68± 4 6± 4.5±.3.47±. 9.58±..5±.9.74±.6.68±.4.3±.5.9± XZ Leo A ± ± 3.69±..±. 6.93±.4.7±..74±.5.59±.3.3±..6± DK Cyg A NN Vir A.48 69± 5 695± 5.7±..5±. 5.9±.6 3.6±.6.73±..85±..95±.8.36± η CrA A ±.3.85± ±.4.±..8±.6.6± ,37 RZ Tau A.46 73± 8 794± 8.56±.7.4±.5 6.9±.6.6±.7.7±.6.64±.6.±.5.± V4 Cyg A AQ Psc A.476 6± 5 64± 5.75±..89±. 3.76±..98±..68±.3.39±..89±.7.8± AH Aur A ± 9 64± 9.85±..89± ±.3.68±.5.8±..±.4.± ,37 V839 Oph A FP Boo A ± ± 4.3±.3.77±..±..9±..6±.5.5±..9±.6.96± V73 Cyg A ± 659±.4±..33±. 8.9±.6.88±..6±..5±..9±.3.± EX Leo A ± 4 6± 4.56±..73±. 3.47±.9.66±.3.57±.3.3±..8±.8.7± AH Cnc A.36 63± 9 675± 9.4±.9.68±.5.78±.5.64±..47±.5.5±.3.85±.6.93± AP Leo A.43 65± 5 65± 5.48±.5.8±.3.8.9±.8.46±.4.43±..8±..± ,37 FG Hya A.38 6± 659±.44±..5±..7±.4.4±..45±.3.6±..79±.4.9± UX Eri A.445 6± 8 634± 8.47±..9±..64±.4.7±.7.43±.3.53±..8±..± YY CrB A ± 64±.43±..8±..59.8±.3.43±.3.35±..85±.5.93± ,37 CK Boo A V566 Oph A.4 7± 4 688± 4.47±..79±. 4.65±.4.5±..4±.3.33±..6±.3.8± TYC74 a A ± ± 36.45±..7±. 3.3±.4.74±..38±..6±..9±.4.83±..44. CN Hyi A ± 636± 6.6±..77±.3 4.±.6.87±..37±.6.5±..98±.6.79± GR Vir A ± 663± 9.4±.7.6±.4.87±.8.48±.6.37±.6.7±.6.83±.3.8± DZ Psc A.366 6± 687±.47±..6±..84±.8.49±..35±.6.8±..83±.8.8± RT LMi A.375 6± ± 58.9±..8±..8±.8.93±.6.3±.5.48±.3.75±.3.88± HN Uma A.383 6± 77 68± 77.44±..58±..55±.3.4±.3.8±.6.8±..76±.6.75± V4 Aur A ± 5983±.44±.3.59±..3±.5.39±..7±.6.7±.3.75±..74± SX Crv A ± 3 66± 3.3±..44±..5.5±..5±.4.±..68±.5.7± ,43 EQ Tau A ± 585±.4±..79± ±..±.3.54±..47±.9.9± ,37 Y Sex A ±..75±. 3.±.44.69±..±.5.±.3.9±.5.64± V357 Oph A ± 5 578± 5.39±..69±..78±.3.47±..9±..9±..68±.4.7± V44 Peg A QX And A.4 65± 5 64± 55.4±..7± ±.9.8±.7.4±.4.9±.7.6± VZ Lib A ± 598±.7±.5.7±.3.36±.9.59±.5.6±.6.35±.3.7±.3.6± OO Aql A.57 57± 3 568± 3.38±..8±..8±.38.55±.3.5±..88±..3±.33.9± V58 Oph A DX Tuc A ± 37 68± 37.±.4.7±..97±.5.66±.4.±.3.3±..8±.6.49± TV Mus A ± 3 688± 3.47±.7.67±.8.48±.4.55±..94±.4.3±..9±..34± XY Boo A ±.5.55± ±.34.5±.5.98±.5.3±..6.,37 TZ Boo A ± ± 65.96±.3.4±.5.98±.6.6±.3.7±.5.5±.4.47±.9.8± DN Cam W ± 3 67± 3.78±..±. 5.6±.3.67±.6.85±..8±..88±.9.5± V78 Her W V4 Aur W.64 67± ± 3.98±..9± ±.7.64±.5.33±..±.8.± ,43 EF Boo W ET Leo W.347 5± 8 55± 8.36±..84±..±..56±..59±..54±..39±.5.3± AA UMa W

11 On the origin of W UMa type contact binaries - a new method for computation of initial masses Table A: continued from previous page Star t P T T R R L L M M M,i M,i /q i t MS Ref d (K) (K) (R ) (R ) (L ) (L ) (M ) (M ) (M ) (M ) YY Eri W.3 537± 5 563± 6.±.6.77±.4.7±.7.5±.8.55±.4.6±.7.6±.45.37± ER Ori W ,37 V5 Oph W ± 7 666± 84.5±.4.93±.6.45±.4.±.8.5±.8.56±.8.73±.36.± V87 Ara W.4 586± 35 6± 35.67±..6±..96±.3.5±..5±..±..89±.4.9± BB Peg W ± 7 6± 7.8±..8±..6±.3.8±..4±..55±..57±.4.8± VY Sex W ± 5 596± 5.5±..86±..7±..83±..4±..45±..73±.4.99± V47 Aql W AE Phe W ± 7 66± 8.9±.3.8±..74±.8.83±.3.38±.6.63±..44±.8.± EZ Hya W UV Lyn W ± 66±.39±..89±..3.9±.3.36±..5±..8±.5.9± ,37 AH Vir W V84 Her W W UMa W ± 78 63± 9.7±.3.85±..6±.6.±.7.35±.9.69±.5.45±.3.9± QW Gem W ± 6±.6±.3.75±..7.69±.4.3±.4.44±..66±.6.9± ,37 SS Ari W V75 Cen W ,37 AM Leo W.366 6± 5 6± 5.3±..85±..84±..95±..9±..59±..59±.4.95± V78 Tau W ± 5 65± 5.±..85±..58±.4.93±.4.9±.7.57±.3.6±.3.94± ,38 V9 Cyg W ±.8.3±..85±.3.7± RZ Com W ± ± 98.±.3.78±..6±.6.93±.5.3±.9.55±.4.65±.3.86± GM Dra W ± ± 58.5±..6±..9±.7.56±..±.4.±..83±.7.67± FU Dra W.37 58± 8 633± 8.±..6±..7.46±..7±.4.9±..68±.5.7± ,37 U Peg W ± 584±.±..74± ±..5±..38±..66±.5.7± ,37 GZ And W.35 58± 5 6± 5.±..74±..3.73±.3.±..59±..44±..83± ,37 AO Cam W ± 33 59± 33.9±..73±..3±.7.58±..±..49±..5±.5.78± ,7 TW Cet W BV Dra W ± ± 3.±..76±..7.84±.4.4±..43±..76±.5.59± ,37 BX Peg W.8 53± 558±.97±..6±..66±..3±..±..38±.5.4±..68± OU Ser W ± 4 638± 4.9±..5± ±..±..8±..73±.5.5± ,37 AB And W ± ± 7.5±..76±..87±..59±..±..49±..5±.4.67± SW Lac W TY Boo W.37 58± 68±.±..69±..±.4.6±.9.93±..4±..66±..5± VW Cep W.78 5± 3 548± 4.9±..6±..47±.8.6±.4.93±..4±..7±..64± BW Dra W.9 598± 664±.98±..55±...39±..9±..6±..65±.5.45± ,37 TX Cnc W ± 7 665± 84.±.3.83±.4.35±..9±.5.9±..5±.6.7±.9.5± V757 Cen W V89 Her W ± 538±.7±..74±..5.4±..86±..37±..53±.5.47± ,43 CC Com W. 4± 8 43± 8.7±..53±..4±..9±..7±..38±..86±..56± XY Leo W ± 8 48± 8.83±.3.66±..±.8.±.7.76±.5.46±.6.3±.46.57± V53 Cas W.34 44± ± 5.74±.4.55±..8±..3±..75±.3.38±..±.3.54± BH Cas W ± 498±.5±..75±..5±..3±..74±.6.35±.3.43±..38±.3.6.8,37 RW Dor W RW Com W.37 5± 54±.7±..46±..3.6±..56±.6.±.3.4±..5± ,37 a) TYC Ref.: ; Ahn, Hill & Khalesseh (99), ; Awadalla & Yamasaki (984), 3; Baran et al. (4), 4; Barone et al. (993), 5; Chochol et al. (), 6; Csák et al. (), 7; Gazeas et al. (5), 8; Gazeas et al. (6), 9; Goecking & Duerbeck (993), ; Goecking et al. (994), ; Gürol et al. (a), ; Gürol et al. (b), 3; Hilditch & King (986), 4; Hilditch, King & McFarlane (988), 5; Hilditch, King & McFarlane (989), 6; Hrivank (988), 7; Hrivank et al. (), 8; Kallrath et al. (6), 9; Kaluzny & Rucinski (986), ; Köse et al. (),; Kreiner et al. (3), ; Metcalfe (999), 3; Milone et al. (987), 4; Niarchos & Manimanis (3), 5; Özkardeş, Erdem & Bak (9), 6; Özkardeş & Erdem (), 7; Pribulla et al. (999), 8; Pribulla & Vaňko (), 9; Qian & Yang (4), 3; Rainger, Hilditch & Bell (99), 3; Ren et al. (), 3; Samec & Hube (99), 33; Szalai et al. (7), 34; Ulaş et al. (), 35; Vaňko et al. (), 36; Yakut et al. (4), 37; Yakut & Eggleton (5), 38; Yakut et al. (5), 39; Yakut et al. (9), 4; Yang & Liu (3a), 4; Yang & Liu (3b), 4; Zhang & Zhang (4), 43; Zola et al. (4), 44; Zola et al. (5), 45; Zola et al. ()