• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Physical Parameters of W UMa Type Contact Binaries and Their Stability of Mass Transfer

    2023-05-29 10:12:42BerikolTekesteGebreyesusandSebluHumneNegu

    Berikol Tekeste Gebreyesus and Seblu Humne Negu

    1 Space Science and Geospatial Institute(SSGI),Entoto Observatory and Research Center(EORC)Astronomy and Astrophysics Department,P.O.Box 33679,Addis Ababa, Ethiopia; berikoltekeste@gmail.com

    2 Addis Ababa University, P.O.Box 1176, Addis Ababa Ethiopia; seblu1557@gmail.com

    Abstract In this study, we determined the physical parameters of W UMa type contact binaries and their stability of mass transfer with different stellar mass ranges over a broad space by applying the basic dynamical evolution equations of the W UMa type contact binaries using accretor and donor masses between 0.079 and 2.79 M⊙.In these systems,we have studied the three subclasses of W UMa systems of A-,B-and W-type contact binaries using the initial and final mass ranges and we investigated different stellar and orbital parameters for the subclasses of W UMa systems.We examined the stability of the W UMa type contact binaries using the orbital parameters such as critical mass ratio, Roche lobe radius of the donor star and mass ratio of these systems.Thus, we computed the observed and calculated physical parameters of A-,B-and W-type W UMa systems.Moreover,we determined the combined and color temperatures to classify the three subclasses of the systems.Also, we presented the result of the internal stellar structure and evolution of W UMa type contact binaries by using the polytropic model.

    Key words: (stars:) binaries (including multiple): close – (stars:) binaries: general – stars: evolution – stars:fundamental parameters –stars: early-type –stars: late-type –stars: low-mass –stars: luminosity function – mass function – stars: mass-loss

    1.Introduction

    A contact binary (CB) is a binary star whose two components are so close to each other that they merge to share their gaseous envelopes and they are gravitationally bound to each other (Duerbeck 1984; Muller & Kempf 1903).There are special types of CB stars, the most ample of which are W Ursae Majoris(W UMa)type(Mochnacki 1981).The W UMa type CBs are called low mass binaries and they include two stars bound together and orbiting around their barycenter(Boyajian et al.2016).All stars in the W UMa type CBs are dependent on the mass transfer between the components of the three subclasses (A-, B-, and W-type) of W UMa system via the physical principle of Roche lobe overflow (RLOF).

    In the subclass of W UMa type CBs, the components of the two stars share the same Roche lobe surface as a consequence of these systems.The stars in these systems fill their Roche lobe then gas flows from the outer layers of that star through the inner Lagrangian point that connects the two Roche lobes.The Roche lobe allows the mass transfer between components of W UMa type CBs(Eggleton 1983).The structures of stars in these systems have different properties from a single star,because of the mass transfer between the components of these systems.

    The physical parameters of W UMa type CBs were studied with the evolution and stellar structure of the stars using the initial and final mass ranges of the W UMa systems.Depending on the improvements in stellar theories and observational data of W UMa type CBs,there are the mass radius relation(MRR)and mass luminosity relation (MLR) which have been studied by Demircan & Kahraman (1991), Cester et al.(1983),Gimenez & Zamorano (1985), Yildiz (2014), Sun et al.(2020).But in this study we focus on the investigation of those relations and other relations such as luminosity temperature relation (LTR), radius temperature relation (RTR) and mass temperature relation(MTR)using the theoretical zero-age main sequence (ZAMS) and terminal-age main sequence (TAMS)models with empirical relations.

    Here, in this study we utilized the catalog of Yildiz (2014)and the late-type W UMa CB catalog of Sun et al.(2020).The catalog of Yildiz (2014) has different subclasses of W UMa type CBs such as A-type and W-type systems (Binnendijk 1970; Van Hamme 1982) whereas the catalog of Sun et al.(2020) has different subclasses of W UMa type CBs such as A-, B-, and W-type systems (Lucy & Wilson 1979).In these two catalogs, the primary components in the A-type are hotter than the secondaries, so an A-type binary has a hotter star and their orbital period is between 0.4 and 0.8 days (Binnendijk 1970)whereas the primary component of W-type is cooler than the secondary component and it has a shorter orbital period in between 0.22 and 0.4 days (Binnendijk 1970).The spectral types of A-and W-type W UMa type CBs are A-F and G-K (Webbink 2003; Rucinski 2004), respectively.Moreover,in the catalog of Sun et al.(2020)the B-type is in poor thermal contact, although it has different effective temperatures of the components which are significant (Lucy & Wilson 1979, Sun et al.2020).The B-type of late-type W UMa CBs comes from the similarity of their light curves (Lucy & Wilson 1979).

    In this study, we computed the observed and calculated parameters of 100 and 2314 stars from the catalog of Yildiz(2014) and Sun et al.(2020), respectively.In these two catalogs, we have studied orbital and stellar parameters to examine the stability of mass transfer and the evolution of the stars.

    In the present work,we have examined the orbital and stellar parameters of the three subclasses of W UMa type CBs using the initial and final mass range of these systems.Hence, we computed angular momentum loss (AML), critical mass ratio,Roche lobe radius of the donor star,mass ratio,semimajor axis for the accretor and donor stars, radius, luminosity and effective temperature, age of the stars, and combined and color temperature using the range of initial and final masses between 0.079 and 2.79 M⊙(Yildiz 2014; Sun et al.2020).Also, we investigated the internal stellar structure of W UMa type CBs using the polytropic model of the system which studies the internal part of the donor stars by introducing stellar structure models in W UMa type CBs (Chandrasekhar &Chandrasekhar 1957).

    The aim of this study is to make a statistical analysis of orbital and stellar parameters of the W UMa type CBs by using the catalog of Yildiz (2014), Sun et al.(2020).Also, we examined the internal stellar structure of W UMa type CBs employing the polytropic model.In these systems, we determined the evolution and stability of mass transfer by applying the ZAMS and TAMS models.

    This paper is organized as follows:In Section 2 we describe the orbital and stellar parameters of W UMa systems.In Section 3 we make a statistical analysis of W UMa type CBs’parameter distributions, construct distributions of the stellar and orbital parameters of the systems between the theoretical ZAMS and TAMS models, and compare them with the empirical relations.In Section 4, we present the results and also we parameterize the parameters of those selected systems from the catalog of W UMa type CBs, and then compare the results with the theoretical relations.The conclusions are drawn in Section 5.

    2.Orbital and Stellar Parameters of W UMa Type CBs

    2.1.Basic Equations of Orbital Parameters

    Variation of the orbital parameters with loss of mass from stars is investigated.We consider the accretor mass (Ma) and donor mass(Md)orbiting each other under the force of gravity.The orbital angular momentum of the W UMa system is given by Soberman et al.(1997), Siess et al.(2011), Eggen (1961)

    A system should satisfy Kepler’s third law to express the orbital period and semimajor axis of the W UMa type CBs.Then, the orbital period of a W UMa system is expressed as

    2.1.1.Roche Lobe Model

    The Roche lobe is a region where the material is bound to the star by gravity.The equivalent radius of the RLOF of W UMa type CBs is defined based on the radius of a sphere with the same volume and range of mass ratios.The radius of the Roche lobe is given by Eggleton (1983)

    where RLdis the Roche lobe radius of the donor star.

    It is difficult to distinguish between the two stars for the purpose of calculating the masses and radii.A simpler approximation of Roche lobe radius of the donor star in Equation (3) which works well for 0.1 <q <0.8 is stated by Paczynski (1971).

    2.1.2.Age of W UMa Type CBs

    The basic thing needed to derive the age of the stars is the process of mass transfer between the two components of the W UMa type CBs.In a W UMa system, the mass is being transferred from the more massive accretor star to the less massive donor star, so the age of the W UMa type CB will be changed (Bilir et al.2005).The star’s age cannot be measured,and it is only estimated based on the approximation of the main sequence(MS)relation(Soderblom 2010;Van Eylen et al.2016;Negu & Tessema 2018).

    2.1.3.Critical Mass Ratio of W UMa Type CBs

    In a W UMa system, the critical mass ratio mainly depends on the mass transfer,AML and stellar wind of the W UMa type CBs (Chen & Han 2008; Arbutina 2009).The critical mass ratio for W UMa type CBs is predicted as

    where Mcis core mass of the CB star as noted by Chen&Han(2008),and we used the core mass Mc=0.35 M⊙for the CBs.

    2.2.Stellar Parameters of the W UMa Type CBs

    The value of each mass in these systems strongly depends on the evolutionary phase of the mass losing star at the beginning of the mass loss transfer and the final stage process of A-, Band W-type W UMa CBs.

    Case 1: MRR and MLR- To determine the radius of the W UMa type CBs,we used the initial and final mass from the two catalogs.The radius of the star was increased or decreased based on the mass transfer between the components of these systems (Gimenez & Zamorano 1985).

    Thus,we investigate the radius based on the mass of A-,Wand B-type W UMa CBs.We illustrate the result of the initial and final luminosity for A-, B- and W-type W UMa CBs(Cester et al.1983) by calculating the luminosity of the stars and we compare them with the luminosity of the Sun(L⊙=3.846×1026W) (Demircan & Kahraman 1991,Harmanec 1988).

    Case 2: LTR and RTR- If the stars have low mass,the luminosity will be decreased, consequently the temperature will also be affected (Demircan & Kahraman 1991;Cester et al.1983; Gimenez & Zamorano 1985).To get the effective temperature of W UMa type CBs,we used the Stefan-Boltzmann law for the luminosity of stars.

    2.3.Combined and Color Temperature

    In the previous studies, Sun et al.(2020) made statements about the combined and color temperatures, but they did not calculate these parameters.Rather, they only give direction on how to calculate the two parameters.The combined and color temperatures are used to examine the subclasses of W UMa type CBs and their spectral type with the period luminosity relation.

    2.4.Stability of Mass Transfer in W UMa Type CBs

    The stability of the mass transfer depends on the donor mass loss and donor radius of the donor star and how the conservative process operates, and also depends on the AML process as well as Roche lobe of the donor star of W UMa type CBs.For the mass transfer to continue, the donor must remain within the Roche lobe and then the W UMa type CBs are stable(Hjellming & Webbink 1987; Webbink 1985).Two things are happening at the same time that affect the stability of mass transfer in W UMa type CBs:

    1.The mass-losing star changes its radius and,

    2.The AML could be transferred between the components of the subclass of W UMa type CBs via the donor response

    2.5.Internal Stellar Structure Equation for W UMa Type CBs

    In this section,we investigate the stellar structure equation of the donor stars in W UMa type CBs which is used to study the internal part of the stars with the polytropic model.In the polytropic model, different numbers of indexes are used to explain the internal stellar structure of the donor stars and the relation between the physical parameters of W UMa type CBs with stable mass transfer(Lane 1870;Emden 1907;Sirotkin&Kim 2009).A polytropic model of W UMa type CBs is applied to investigate how the pressure of the stars changes with density inside as one moves through the W UMa system.However, we describe the polytropic model as an approximation that represents the equation of state of the gas inside stars (Chandrasekhar & Chandrasekhar 1957).The polytropic stars have a relation between the pressure and density of W UMa type CBs.Using the central density of W UMa type CBs,ρc, we have the polytropic model, which is expressed as

    The appropriate boundary conditions from applying a polytropic model relation using Equation (14) are at the center of W UMa type CBs.

    When pressure and density go to zero in Equation (14), the surface of the polytropic model is located at ξ=ξ1, where θ(ξ1)=0.There is an analytical solution for different polytropic indexes n=0, 1, 5.

    3.Physical Parameter Distribution of W UMa Type CBs

    In this section, we are applying the basic equation of the mass transfer and stellar structure models, then we design the theoretical framework for the evolutionary computations of mass transfer in W UMa type CBs.The analysis relies on our main sample of data on W UMa type CBs from the catalogs of Yildiz (2014) and Sun et al.(2020), and we summarize the details of the observations and calculated parameters in Tables 2 and 3.

    3.1.Distribution of Orbital Parameters

    In the evolution of W UMa systems, there are different orbital parameter distributions of W UMa type CBs.

    Some orbital parameters for these systems are semimajor axis, orbital period, Roche lobe radius of the donor star and orbital angular momentum of the system.We determine the orbital parameters by using different initial and final mass ranges and comparing their results with the numerical solution of W UMa type CBs in Section 4.1.

    The left side of Figure 1 shows distribution of the semimajor axis in late-type W UMa CBs with the different ranges of final donor masses.The result of this plot shows the distribution of semimajor axis with variation of color lines.The variation of colors is used to identify the subtype of late-type W UMa CBs.In this result, the semimajor axis of late-type W UMa systems shows the expansion and shrinking of A-, B- and W-type latetype W UMa type CBs.The solid black line represents W-type stars with its peak at 172.0, and the corresponding range is between 1.936 R⊙and 2.103 R⊙.

    The right panel of Figure 1 shows the distribution of critical mass ratio for the final stage of donor mass of late-type W UMa CBs.Based on the mass of the stars from the catalog of Sun et al.(2020) and the core mass of CB stars, we calculated the critical mass ratio of late-type W UMa CBs utilizing Equation (6).The critical mass ratio is used to determine the stability of the mass transfer in W UMa systems.The peak point for the solid black line is 148 which indicates the number of W-type CBs that exist at that point and qcris between 0.625 and 0.6417.

    3.2.Distribution of Stellar Parameters

    In the evolution of W UMa systems there are different stellar parameter distributions of W UMa type CBs.Some stellar parameters of these systems are luminosity, radius, effective temperature, and combined and color temperature.

    Figure 1.The distribution of calculated semimajor axis with the different ranges of donor mass by applying Equation (2).The right side shows the distribution of critical mass ratio of A-, B- and W-type late-type W UMa CBs with the donor mass Md <0.7 M⊙using the catalog of Sun et al.(2020).

    Figure 2.The distribution of the accretor effective temperature with different ranges of donor masses for the left side of Figure 2 and distribution of the donor effective temperature using different ranges of initial accretor mass of W UMa systems based on the catalogs of Sun et al.(2020)and Yildiz(2014)by applying Equation(9).

    The left panel of Figure 2 displays the distribution of accretor effective temperature using the donor mass range Md<0.7 M⊙.The A-, B- and W-type late-type W UMa systems have different effective temperatures in the MS regions.The result, the solid black line, features a peak point at 189.0 which is in between 3.645 K and 3.650 K for W-type systems.The solid red line corresponds to B-type systems with high effective temperature and it is a hot star compared to the other types of W UMa CBs.

    The solid blue line shows the A-type star has moderate effective temperature and the solid black line of W-type of latetype W UMa CBs has cool effective temperature.The W UMa type CBs for A- and W-type have identical effective temperature, but there are slight differences between them(Yildiz 2014).The right panel of Figure 2 shows the variation of effective temperature of the donor stars for A-and W-type of W UMa CBs.The black dashed and solid black line distribution depicts the effective temperature for an A-type W UMa system and it has cool effective temperature compared to the corresponding W-type W UMa system.The W UMa type CBs,represented by the solid blue line of W-type,have a peak at 6, which is in between 3.539 K and 3.567 K.

    4.Result and Discussion

    4.1.Numerical Solution of W UMa Type CB Parameters

    Figure 3.The numerical result of the orbital period and AML of the late-type W UMa CBs using the catalog of Sun et al.(2020)and mass age relation with different range accretor mass of W UMa type CBs using the catalog of Yildiz (2014).

    The orbital angular momentum is used to determine the evolution of W UMa type CBs.The rate of change of AML increases as the orbital separation and orbital period decrease and then for AML decreases the total mass of the system also decreases (Smith 2006; Patterson 1984).The left side of Figure 3 shows the numerical solution of orbital period and orbital AML relation by applying Equations (1) and (2) using different ranges of orbital periods Porb=0.80,0.60 and 0.40 d.AML has direct or indirect relation with different orbital parameters, although for short orbital period the mass is high and the AML is also high, so orbital period and AML have an indirect relation.As we observed from the result, the short orbital period of the systems and the AML mechanism were governed by gravitational radiation (Patterson 1984; Taam &Spruit 2001).

    The right panel of Figure 3 features the numerical solution of the final mass distribution and age relation of the W UMa type CBs by applying Equation (5).The dashed black line with Md=0.1 M⊙has low mass and it has slow fusion to use in the core hydrogen burning.It has less pressure, low effective temperature and also low luminosity because these stars have a long lifetime on the MS range.The solid red line represents a more massive star with Ma=0.56 M⊙, so the solid red line corresponds to a short lifetime,with higher central temperature and high pressure to support itself against gravitational collapse and also has fast fusion to use in its core for hydrogen burning.Hence,the more massive stars have a short lifetime on the MS(Soderblom 2010).The A-type shows a long lifetime and the W-type also shows a short lifetime because of the mass difference The W-type is hotter, more massive and more luminous than the A-type W UMa CBs.

    4.2.The Comparison of Initial and Final Stellar Parameters

    Case 1:MLR and MRR-The left side of Figure 4 depicts the distribution of initial mass–luminosity relation of W UMa type CBs depending on the theoretical model of ZAMS and TAMS with an empirical relation.In these relations we incorporated the numerical solution oflogL(a,i) andlogL(d,i) with Ma,iand Md,i, where i is for initial.In the ZAMS model we used x=3.92 for Ma,i≥0.7 M⊙, and for the TAMS model we used y=2.62 and Ma,i<0.7 M⊙.The solid black and red lines representlogL(a,i) andlogL(d,i) respectively with x=3.92,although the solid blue and black dashed line corresponds to log(La,i) andlog(Ld,i) with y=2.62 for different initial mass ranges.

    The right side of Figure 4 shows the theoretical result of MRR using different initial accretor mass ranges.In W UMa systems the more massive stars have large radius and less massive stars have low radius,therefore,the radius and mass of W UMa stars have a direct relation.The radius of W UMa type CBs is utilized to examine the phase of stars in the MS region.Using Equation (9) we calculate the initial and final radius of the W UMa type CBs with the theoretical models of ZAMS and TAMS.

    Case 2:LTR and RTR-The left panel of Figure 5 shows the numerical solution of distribution final donor LTR.Hence,using different donor mass of W UMa systems Md=0.58,0.78 and 0.98 M⊙, we obtained a different luminosity distribution for W UMa type CBs.The solid red line shows a more massive star with high luminosity as well as high temperature.The right panel of Figure 5 displays results of accretor radius RTR using different range accretor masses, Ma=2.79 M⊙, Ma=2.29 M⊙and Ma=2.49 M⊙.However, as we observed from this result when the mass increases, the radius also increases, but the temperature will decrease.As we observed from the result, the solid black line represents a more massive star than the other and it has lower temperature but larger radius.The solid blue line shows low radius as well as high temperature.

    Figure 5.The two numerical results show the luminosity effective temperature relation, and radius effective temperature relation with different ranges by applying Equation (9) using the catalog of Yildiz (2014).

    The left side of Figure 6 shows a solid black, blue and red lines and a dashed black line to represent how the radius responds to the adiabatic (ζad) calculations using theoretical ZAMS and TAMS models with empirical relations by using the donor mass.The right side of Figure 6 displays solid black,red and blue lines to represent how the donor mass responds to the adiabatic (ζad) calculations with constant mass loss rates at Md=0.98, 0.78 and 0.58 M⊙.

    Comparing both plots, the star with lower mass will have a longer stability period than the star with greater mass.In a W UMa system,ζdon≥ζLd, implies stability andζdon<ζLdimplies instability of mass transfer in these systems(ζad=ζdon) (Eggleton 1983).These two plots depend on the donor mass of the star because ζaddepends on the mass loss for the donor star althoughζLdalso depends on the AML processes for W UMa type CBs.

    The left side of Figure 7 shows that the solid black and red line represents the relation between the semimajor axis and orbital period of A-, W- and B-type CBs.Hence, we classified the numerical solutions of the semimajor axis and orbital period based on the donor mass Md<0.7 M⊙of the late-type W UMa systems (Sun et al.2020).In these results the stars are compacted in one region;this shows less semimajor axis and a short orbital period for the W UMa systems.Kepler’s first law states the mass-orbital period relation of the star.

    Figure 6.The two numerical result plots demonstrate radius and donor mass response to the adiabatic(ζad)calculation,with the timescale of the W UMa type CBs by applying Equation (12) using the catalog of Yildiz (2014).

    Figure 7.The comparison of orbital period—semimajor axis relation and mass-orbital period relation for A-, B- and W-type of W UMa CBs with donor mass Md <0.7 M⊙by applying Equation (2) (Sun et al.2020).

    As Kepler’s law states,if a star is more massive it has a short orbital period and when a star has a low mass, it has a long orbital period.Based on Kepler’s laws,the orbital period of the late-type W UMa CBs are dependent on the mass of the stars.The right side of Figure 7 shows the final phase distribution of orbital period-donor mass of the late-type W UMa CBs compared with A-, B- and W-type CBs.The A-type of latetype W UMa CBs has a long orbital period and W-type has short orbital period (Van Hamme 1982).The solid blue line shows the numerical solution of these systems.

    The left side of Figure 8 shows the comparison between the distribution of the combined temperature and the calculated effective temperature of the late-type W UMa CBs that can be obtained from the catalog of Sun et al.(2020).In the result the blue-colors for B-type are dominated by the red-color of A-type, which happens because they have identical combined and effective temperature.

    The right side of Figure 8 shows the comparison between the color and effective temperature relation of the late-type W UMa CBs.However, using color temperature, we have illustrated either the stars are living in the MS region or they are out of the MS region (Pecaut & Mamajek 2013).Stars with high mass will have fast hydrogen burning with short lifetime, and they will have hotter temperature as well as being bluer (Bell et al.2012; Gullbring et al.1998).Some of the stars have red color because they are the stars that decrease in temperature and also have slow hydrogen burning in the core.This aspect makes the stars have red color in the MS(Bell et al.2012;Gullbring et al.1998).

    Figure 9.The result of orbital period and semimajor axis relation in the left panel and the right panel shows the mass and semimajor axis relation that can be obtained from the catalogs of Sun et al.(2020) and Yildiz (2014).In these results, the distribution shows the fitting function by applying the power-law function=a(x)b.

    Table 1 The Properties of the Lane–Emden Solution Using Different Polytropic Indexes

    The left side of Figure 9 shows the distribution of the orbital period and semimajor axis relation that can be obtained from the catalog of Sun et al.(2020) using different parameters of these systems.In the fitting function there are different types of relations.Some of them are linear, exponential, power-law,quadratic and others.In this thesis, we utilized a power-law fitting function.This function was used to investigate the error of the catalog of parameters from the catalog of Sun et al.(2020).

    The right side of Figure 9 displays the fitting function of semimajor axis and mass of the W UMa type CBs using the catalog of Yildiz (2014).To illustrate the error of the given data, we used a power-law fitting function with different fitting parameters“a”and“b.”The solid black line shows the relation between the semimajor axis and the mass of the W UMa type CBs and the solid red line indicates the error function.

    Figure 10.Numerical solution of the internal stellar structure of the W UMa type CBs for a donor star by the using Lane–Emden equation, Equation (14)with different polytropic indexes n=0, 1, 1.5, 2, 3 and 4 (Lane 1870;Emden 1907).

    4.2.1.Internal Stellar Structure of W UMa Type CBs

    We explain the internal stellar structure of W UMa type CBs using the polytropic model.In a polytropic model, different number of indexes are used to explain the internal stellar structure and the relation between pressure, radius, density,temperature and mass transfer in W UMa type CBs with a donor star (Lane 1870; Emden 1907; Sirotkin & Kim 2009).For n=0 the density is constant and the pressure is also constant,when n=1 the density of the star keeps the radius the same no matter what the mass is in the system and also the radius is independent of mass and the temperature is also constant;for n=5 the radius,mass,central density and central pressure are infinite, therefore, the binding energy is infinite.

    Table 1 shows the result of the polytropic case using the Lane–Emden equation for each polytropic index, where

    Figure 10 shows the numerical solution of the Lane–Emden equation by Lane(1870),Emden(1907).The polytropic index is zero and the pressure and density correspond to a constant, and the other analytical and numerical solutions exist for n=1,1.5,2,3,4 and also for 5.Hence,the n ≥5 solution contains infinite mass with the range of 0 <n <5 for the polytropic model.Each colored line shows a different polytropic index in Figure 10.

    5.Conclusions

    The reason for the W UMa type CBs exhibiting long-term evolution is because the W UMa type CBs are known as low mass stars.We have derived the basic evolution equations of mass transfer between the components of W UMa type CBs and we calculated the orbital and stellar parameters for these selected systems from the W UMa type CBs.

    We have presented the evolution of AML and effect of mass transfer in the subclasses of W UMa type CBs based on the catalogs of Yildiz(2014)and Sun et al.(2020).We determined the stability of stars based on the effect of mass transfer in the parameters of Roche lobe radius of the donor stars and critical mass ratio of W UMa systems.

    In this study, we have calculated some stellar and orbital parameters,which are not yet calculated in the catalog of Yildiz(2014) and Sun et al.(2020); some of them were critical mass ratio,Roche lobe radius of the donor star,mass ratio,AML and orbital angular velocity of the stars of the W UMa type CBs.We calculated the combined and color temperature based on the late-type W UMa CB catalog of Sun et al.(2020).These two parameters used to explain the color of stars are based on the temperature whether they are blue, red, yellow, orange or other colors.

    We explain the internal stellar structure of W UMa type CBs using a polytropic model.In the polytropic model, different numbers of indexes are used to explain the internal stellar structure of the donor stars and the relation between the physical parameters of W UMa type CBs with stable mass transfer (Lane 1870;Emden 1907;Sirotkin & Kim 2009).The final results are presented in Tables 2 and 3.

    Table 2 The Observed and Calculated Parameters for A- and W-type W UMa CBs Using Md ≥0.7 M⊙(Yildiz 2014)

    Table 3 The Observed and Calculated Parameters for A-, B- and W-type Late-type W UMa CBs Using Md ≥0.7 M⊙(Sun et al.2020)

    Acknowledgments

    We thank the Space Science and Geospatial Institute(SSGI)-Entoto Observatory and Research Center (EORC) Astronomy and Astrophysics Department for supporting this research.This research has made use of NASA’s Astrophysical Data System.

    美女脱内裤让男人舔精品视频| 精品人妻一区二区三区麻豆| 婷婷六月久久综合丁香| 亚洲av熟女| 国产伦精品一区二区三区视频9| 亚洲无线观看免费| 最近2019中文字幕mv第一页| 国产精品三级大全| 亚洲电影在线观看av| 欧美一级a爱片免费观看看| 好男人视频免费观看在线| 少妇丰满av| 午夜福利高清视频| 国产精品麻豆人妻色哟哟久久 | 我的老师免费观看完整版| 久久精品国产鲁丝片午夜精品| 亚洲五月天丁香| av在线观看视频网站免费| 尾随美女入室| 精品无人区乱码1区二区| 国产黄a三级三级三级人| 日本免费一区二区三区高清不卡| 日本一本二区三区精品| 一本久久精品| 亚洲精品日韩av片在线观看| 色综合色国产| 蜜桃久久精品国产亚洲av| 国产高清有码在线观看视频| 国产黄a三级三级三级人| 国产精品一区二区性色av| 联通29元200g的流量卡| 国产中年淑女户外野战色| 亚洲国产欧洲综合997久久,| 免费无遮挡裸体视频| 日韩 亚洲 欧美在线| 久久久久久久久大av| 午夜免费男女啪啪视频观看| 亚洲精品国产av成人精品| 最近最新中文字幕免费大全7| 日韩欧美 国产精品| 国产日韩欧美在线精品| 美女高潮的动态| 成人欧美大片| av在线蜜桃| 亚洲美女搞黄在线观看| 国产一区亚洲一区在线观看| 国产色婷婷99| 青春草亚洲视频在线观看| 一个人看视频在线观看www免费| av视频在线观看入口| 午夜福利在线观看免费完整高清在| 青青草视频在线视频观看| 国产免费男女视频| 嫩草影院入口| 我的老师免费观看完整版| 免费黄色在线免费观看| 国产精品久久久久久精品电影| 久久精品熟女亚洲av麻豆精品 | 亚洲一区高清亚洲精品| av免费观看日本| 国产又黄又爽又无遮挡在线| 国产黄片视频在线免费观看| 全区人妻精品视频| a级毛色黄片| 91在线精品国自产拍蜜月| 2021天堂中文幕一二区在线观| 免费av不卡在线播放| 青春草视频在线免费观看| 日韩 亚洲 欧美在线| 亚洲精品成人久久久久久| 成人午夜精彩视频在线观看| 亚洲va在线va天堂va国产| 日本黄大片高清| 日本爱情动作片www.在线观看| 国产午夜福利久久久久久| 欧美xxxx黑人xx丫x性爽| 亚洲欧美成人精品一区二区| 成人性生交大片免费视频hd| 日日撸夜夜添| 国产黄片美女视频| 最近的中文字幕免费完整| 最近最新中文字幕免费大全7| 一个人看视频在线观看www免费| 在线免费十八禁| 午夜a级毛片| 麻豆成人av视频| 精品人妻熟女av久视频| 亚洲国产精品成人综合色| 舔av片在线| 久久久亚洲精品成人影院| 亚洲三级黄色毛片| 国产黄色视频一区二区在线观看 | 久久久国产成人精品二区| 嫩草影院新地址| 国产一区二区在线观看日韩| av福利片在线观看| 国产乱人偷精品视频| 亚洲成人久久爱视频| 激情 狠狠 欧美| 97在线视频观看| 免费一级毛片在线播放高清视频| 麻豆成人午夜福利视频| 国产黄a三级三级三级人| 国产亚洲5aaaaa淫片| 免费av不卡在线播放| 夜夜爽夜夜爽视频| 女人被狂操c到高潮| 日韩av在线大香蕉| 精品久久久久久电影网 | 亚洲最大成人av| 欧美精品一区二区大全| 国产老妇伦熟女老妇高清| 亚洲精品影视一区二区三区av| 欧美精品国产亚洲| 黄片wwwwww| 欧美区成人在线视频| av专区在线播放| 99久久人妻综合| 午夜激情欧美在线| 国产免费视频播放在线视频 | 91aial.com中文字幕在线观看| 国产黄片美女视频| 久久久欧美国产精品| 国产高清三级在线| 麻豆av噜噜一区二区三区| 能在线免费看毛片的网站| 综合色丁香网| 黄片wwwwww| 精品酒店卫生间| 国产伦精品一区二区三区四那| 特大巨黑吊av在线直播| 久久久久九九精品影院| 午夜老司机福利剧场| 搡老妇女老女人老熟妇| 国产亚洲91精品色在线| 国产伦精品一区二区三区视频9| 日韩av不卡免费在线播放| 久久久久久九九精品二区国产| 亚洲国产高清在线一区二区三| 国产精品久久久久久久电影| 久久久久久久久久成人| 天堂av国产一区二区熟女人妻| 看片在线看免费视频| 国产激情偷乱视频一区二区| 国产淫语在线视频| 又黄又爽又刺激的免费视频.| 午夜福利成人在线免费观看| 亚洲精品乱码久久久v下载方式| 国产久久久一区二区三区| 非洲黑人性xxxx精品又粗又长| 国产精华一区二区三区| 搡老妇女老女人老熟妇| 99国产精品一区二区蜜桃av| av视频在线观看入口| 一夜夜www| 日本免费a在线| 国产一区二区在线av高清观看| 亚洲综合色惰| 男女啪啪激烈高潮av片| 一个人看视频在线观看www免费| 精品久久久久久久人妻蜜臀av| 精品久久久噜噜| 99热精品在线国产| 免费人成在线观看视频色| 人人妻人人看人人澡| 中文字幕久久专区| 国产成人精品一,二区| 欧美高清性xxxxhd video| 建设人人有责人人尽责人人享有的 | 国内揄拍国产精品人妻在线| 国产三级中文精品| 亚洲精品成人久久久久久| 最新中文字幕久久久久| 蜜臀久久99精品久久宅男| 久久精品人妻少妇| 两个人的视频大全免费| 国产免费一级a男人的天堂| 22中文网久久字幕| 熟女电影av网| 亚洲色图av天堂| 三级经典国产精品| 美女脱内裤让男人舔精品视频| 久久精品91蜜桃| 成人毛片a级毛片在线播放| 精品一区二区三区人妻视频| 寂寞人妻少妇视频99o| 久久久亚洲精品成人影院| 国产麻豆成人av免费视频| 少妇人妻精品综合一区二区| 极品教师在线视频| 亚洲美女搞黄在线观看| 色尼玛亚洲综合影院| 日韩视频在线欧美| 精品人妻一区二区三区麻豆| 久久精品国产鲁丝片午夜精品| 如何舔出高潮| 成人亚洲精品av一区二区| 级片在线观看| 99热这里只有是精品50| 精品久久久噜噜| 又黄又爽又刺激的免费视频.| 日韩欧美国产在线观看| 99久久成人亚洲精品观看| 亚洲三级黄色毛片| 免费人成在线观看视频色| 国产欧美日韩精品一区二区| 亚洲在线自拍视频| 国产日韩欧美在线精品| 亚洲天堂国产精品一区在线| 久久久a久久爽久久v久久| 午夜精品一区二区三区免费看| 国产色爽女视频免费观看| 高清在线视频一区二区三区 | 国产探花在线观看一区二区| 亚洲三级黄色毛片| 国产午夜精品一二区理论片| 小蜜桃在线观看免费完整版高清| 国产精品电影一区二区三区| 国产免费又黄又爽又色| 国产精品麻豆人妻色哟哟久久 | 免费无遮挡裸体视频| 久久精品国产自在天天线| 亚洲欧美日韩高清专用| 性色avwww在线观看| 精品欧美国产一区二区三| 特大巨黑吊av在线直播| 国产精品电影一区二区三区| 少妇高潮的动态图| 国产高潮美女av| av国产免费在线观看| 麻豆国产97在线/欧美| 99久久无色码亚洲精品果冻| 亚洲成人久久爱视频| 嘟嘟电影网在线观看| 久久这里有精品视频免费| 成年版毛片免费区| 91久久精品电影网| www日本黄色视频网| 高清av免费在线| 亚洲婷婷狠狠爱综合网| 亚洲精品乱码久久久v下载方式| 一级毛片aaaaaa免费看小| 国产精品99久久久久久久久| 九九在线视频观看精品| 久久99蜜桃精品久久| 国产av不卡久久| 卡戴珊不雅视频在线播放| 最近的中文字幕免费完整| 精品99又大又爽又粗少妇毛片| 青春草国产在线视频| 欧美高清性xxxxhd video| 久久久国产成人精品二区| 午夜久久久久精精品| 国产伦一二天堂av在线观看| 寂寞人妻少妇视频99o| 日本一本二区三区精品| 成人三级黄色视频| 热99在线观看视频| 性插视频无遮挡在线免费观看| 一夜夜www| 日韩高清综合在线| 少妇熟女欧美另类| 精品久久久噜噜| 在线观看66精品国产| 久久久亚洲精品成人影院| 最近最新中文字幕免费大全7| 成人漫画全彩无遮挡| 日本欧美国产在线视频| 久久精品国产99精品国产亚洲性色| 国产高清三级在线| 亚洲国产精品合色在线| 亚洲激情五月婷婷啪啪| 日韩精品有码人妻一区| 一级毛片久久久久久久久女| 99久久精品国产国产毛片| 欧美成人免费av一区二区三区| 中文字幕av在线有码专区| 国产黄片美女视频| 黄片无遮挡物在线观看| 欧美一区二区国产精品久久精品| av卡一久久| 欧美3d第一页| 天堂网av新在线| 三级经典国产精品| 深夜a级毛片| 国产精品一区二区三区四区久久| 两个人视频免费观看高清| 精品无人区乱码1区二区| 黄色欧美视频在线观看| 99视频精品全部免费 在线| 91av网一区二区| 综合色av麻豆| 国产淫片久久久久久久久| 国产成人精品一,二区| 免费无遮挡裸体视频| 亚洲精品日韩av片在线观看| 狂野欧美激情性xxxx在线观看| 中文字幕精品亚洲无线码一区| 午夜精品一区二区三区免费看| 欧美变态另类bdsm刘玥| 国产一级毛片七仙女欲春2| 色吧在线观看| 男女视频在线观看网站免费| 极品教师在线视频| 久久久久免费精品人妻一区二区| 天天躁夜夜躁狠狠久久av| 国产私拍福利视频在线观看| 日日摸夜夜添夜夜爱| 日本黄色片子视频| 波多野结衣巨乳人妻| 久久精品91蜜桃| 日韩欧美三级三区| 成人无遮挡网站| 欧美成人精品欧美一级黄| 国产精品不卡视频一区二区| 永久网站在线| 亚洲av男天堂| 国产在视频线精品| 亚洲精品乱码久久久v下载方式| 日本熟妇午夜| 日本午夜av视频| 最近最新中文字幕免费大全7| 亚洲av熟女| 精品99又大又爽又粗少妇毛片| 亚洲最大成人中文| 亚洲成人中文字幕在线播放| 国产免费视频播放在线视频 | 国产精品福利在线免费观看| 精品午夜福利在线看| 色播亚洲综合网| 成人亚洲精品av一区二区| 欧美人与善性xxx| 亚洲精品亚洲一区二区| 欧美性猛交黑人性爽| 国产精品久久久久久精品电影小说 | 亚洲在线自拍视频| 免费在线观看成人毛片| 一级毛片aaaaaa免费看小| 啦啦啦韩国在线观看视频| 国内精品一区二区在线观看| 人人妻人人澡欧美一区二区| 国产成人福利小说| 成人性生交大片免费视频hd| 欧美一级a爱片免费观看看| 国产熟女欧美一区二区| 菩萨蛮人人尽说江南好唐韦庄 | 美女cb高潮喷水在线观看| 蜜桃久久精品国产亚洲av| 看片在线看免费视频| 身体一侧抽搐| 国产精品一区二区三区四区久久| 18禁裸乳无遮挡免费网站照片| 亚洲精品乱码久久久v下载方式| 网址你懂的国产日韩在线| 日本色播在线视频| 午夜亚洲福利在线播放| 99久国产av精品| 日本熟妇午夜| 26uuu在线亚洲综合色| 欧美激情在线99| 18禁动态无遮挡网站| 看十八女毛片水多多多| av国产免费在线观看| 一二三四中文在线观看免费高清| 亚洲国产成人一精品久久久| 麻豆成人av视频| 看非洲黑人一级黄片| 偷拍熟女少妇极品色| 青春草国产在线视频| 国产成人freesex在线| 婷婷色综合大香蕉| 色播亚洲综合网| 亚洲精品乱码久久久v下载方式| 一个人看的www免费观看视频| www.色视频.com| 禁无遮挡网站| 别揉我奶头 嗯啊视频| 日本黄色视频三级网站网址| 99久久中文字幕三级久久日本| 国模一区二区三区四区视频| 色尼玛亚洲综合影院| 精品不卡国产一区二区三区| 免费av毛片视频| 性插视频无遮挡在线免费观看| 天堂影院成人在线观看| 老司机影院毛片| 99久久中文字幕三级久久日本| 中文乱码字字幕精品一区二区三区 | 日本av手机在线免费观看| 久热久热在线精品观看| 亚洲av二区三区四区| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲欧美精品综合久久99| 夫妻性生交免费视频一级片| 亚洲av男天堂| 日韩亚洲欧美综合| 精品国内亚洲2022精品成人| 亚洲精品乱码久久久久久按摩| 国产亚洲最大av| 在线免费观看不下载黄p国产| 亚洲国产日韩欧美精品在线观看| 欧美bdsm另类| 久久国内精品自在自线图片| 99热精品在线国产| 亚洲真实伦在线观看| 深爱激情五月婷婷| 欧美成人精品欧美一级黄| 国语自产精品视频在线第100页| av又黄又爽大尺度在线免费看 | 亚洲精品一区蜜桃| 精品国产一区二区三区久久久樱花 | 成人性生交大片免费视频hd| 国产亚洲精品av在线| 美女大奶头视频| 免费看美女性在线毛片视频| 国产精品伦人一区二区| 一个人看的www免费观看视频| 美女黄网站色视频| 简卡轻食公司| 欧美日韩国产亚洲二区| 哪个播放器可以免费观看大片| 亚洲精品影视一区二区三区av| 男女那种视频在线观看| 天堂中文最新版在线下载 | 能在线免费看毛片的网站| 国产成人精品一,二区| 国产三级中文精品| 欧美激情在线99| 大又大粗又爽又黄少妇毛片口| 国产三级中文精品| 一二三四中文在线观看免费高清| 真实男女啪啪啪动态图| 日韩一本色道免费dvd| 69人妻影院| 久久人妻av系列| 午夜福利网站1000一区二区三区| 成人漫画全彩无遮挡| 别揉我奶头 嗯啊视频| 精品人妻偷拍中文字幕| 岛国在线免费视频观看| 狠狠狠狠99中文字幕| 午夜福利视频1000在线观看| 少妇的逼水好多| 久久99蜜桃精品久久| 日韩欧美精品免费久久| 午夜久久久久精精品| 色网站视频免费| 日韩强制内射视频| 免费看日本二区| 夜夜看夜夜爽夜夜摸| 激情 狠狠 欧美| 久久精品国产亚洲av天美| 在线天堂最新版资源| 床上黄色一级片| 特大巨黑吊av在线直播| 亚洲电影在线观看av| 一级毛片电影观看 | 亚洲在线自拍视频| 久久久久久九九精品二区国产| 一个人免费在线观看电影| 一级爰片在线观看| 久久久午夜欧美精品| 卡戴珊不雅视频在线播放| 午夜精品在线福利| 午夜免费激情av| 国产亚洲最大av| 国产精品综合久久久久久久免费| 国产探花在线观看一区二区| 久99久视频精品免费| 久久久久久久久久成人| 最新中文字幕久久久久| 水蜜桃什么品种好| 日本熟妇午夜| 久久人人爽人人爽人人片va| 水蜜桃什么品种好| 成人亚洲欧美一区二区av| 99九九线精品视频在线观看视频| 国产精品麻豆人妻色哟哟久久 | 成人av在线播放网站| 高清视频免费观看一区二区 | 黄片wwwwww| 国产精品嫩草影院av在线观看| 一个人看视频在线观看www免费| 国产高清视频在线观看网站| 日本黄色片子视频| 欧美区成人在线视频| 日韩成人av中文字幕在线观看| 青春草视频在线免费观看| 男人舔女人下体高潮全视频| 在线免费观看不下载黄p国产| 你懂的网址亚洲精品在线观看 | 国产激情偷乱视频一区二区| 哪个播放器可以免费观看大片| 51国产日韩欧美| 国产白丝娇喘喷水9色精品| 青春草亚洲视频在线观看| 国产日韩欧美在线精品| 深夜a级毛片| 国产精品一区www在线观看| 免费无遮挡裸体视频| 欧美性猛交╳xxx乱大交人| 久久久精品大字幕| 久久婷婷人人爽人人干人人爱| 99热全是精品| 九九在线视频观看精品| 超碰av人人做人人爽久久| 联通29元200g的流量卡| 热99在线观看视频| 久久精品久久久久久久性| 精品一区二区三区人妻视频| 18禁在线播放成人免费| 国产久久久一区二区三区| 亚洲欧美精品专区久久| 午夜精品在线福利| 两性午夜刺激爽爽歪歪视频在线观看| 最近最新中文字幕免费大全7| 桃色一区二区三区在线观看| 亚洲中文字幕一区二区三区有码在线看| 欧美性猛交╳xxx乱大交人| 观看免费一级毛片| 成人三级黄色视频| 国产在视频线精品| 久久亚洲国产成人精品v| 色综合亚洲欧美另类图片| 国产真实乱freesex| h日本视频在线播放| 欧美一区二区国产精品久久精品| 精品久久久久久久久久久久久| 2022亚洲国产成人精品| 亚洲五月天丁香| 国产精品野战在线观看| 午夜a级毛片| 午夜激情福利司机影院| 你懂的网址亚洲精品在线观看 | 国产在线男女| 97超视频在线观看视频| 麻豆av噜噜一区二区三区| 边亲边吃奶的免费视频| 国产免费视频播放在线视频 | 18禁动态无遮挡网站| 久久欧美精品欧美久久欧美| 国产真实乱freesex| 国产老妇女一区| 91在线精品国自产拍蜜月| 成人美女网站在线观看视频| 黄片wwwwww| 桃色一区二区三区在线观看| 欧美又色又爽又黄视频| 久久久国产成人免费| 中文亚洲av片在线观看爽| 久久综合国产亚洲精品| 少妇的逼水好多| 精华霜和精华液先用哪个| 91精品国产九色| av播播在线观看一区| 国产真实乱freesex| 亚洲中文字幕日韩| 亚洲高清免费不卡视频| 日韩,欧美,国产一区二区三区 | 久久久久久久久久久丰满| 最新中文字幕久久久久| 成人亚洲精品av一区二区| 亚洲欧美日韩无卡精品| 桃色一区二区三区在线观看| 亚洲精品久久久久久婷婷小说 | 亚洲在线观看片| 国产精品三级大全| 看黄色毛片网站| 美女黄网站色视频| 真实男女啪啪啪动态图| 国产精品人妻久久久影院| 又黄又爽又刺激的免费视频.| 国产成人91sexporn| 又爽又黄a免费视频| 精品久久久久久电影网 | 成人三级黄色视频| 国产单亲对白刺激| 免费看a级黄色片| 99热全是精品| 亚洲精品,欧美精品| 国产69精品久久久久777片| 亚洲精品乱码久久久久久按摩| 国产美女午夜福利| 亚洲精品乱码久久久久久按摩| 国产成人aa在线观看| 日本免费一区二区三区高清不卡| 日韩亚洲欧美综合| 国产女主播在线喷水免费视频网站 | av在线观看视频网站免费| 国产人妻一区二区三区在| 精品人妻偷拍中文字幕| 亚洲欧美中文字幕日韩二区| ponron亚洲| 水蜜桃什么品种好| 神马国产精品三级电影在线观看| 久久人人爽人人片av| 亚洲精品,欧美精品| 免费看av在线观看网站| 男人舔奶头视频| 韩国高清视频一区二区三区| 国产精品国产三级专区第一集| 色视频www国产| 亚洲欧美日韩无卡精品| 一级av片app| 久久亚洲精品不卡| 中文字幕人妻熟人妻熟丝袜美| 久久这里只有精品中国| 尤物成人国产欧美一区二区三区| 亚洲国产精品专区欧美| 国产色爽女视频免费观看| 国产精品,欧美在线| 十八禁国产超污无遮挡网站| 久久人人爽人人爽人人片va| 成人午夜精彩视频在线观看| 男女国产视频网站|