Low-Mass Close Binary: KIC 9761199
In this study, we present the results obtained from KIC 9761199’s the photometrical data acquired by the Kepler Mission. The light curve of the system, the sinusoidal variation out-of-eclipses and instant-short term flare events in the entire light curves were analyzed. The temperature of the secondary component was found to be 38911 K, while the mass ratio of the components () was found to be 0.690.01, and the orbital inclination () was computed as . The sinusoidal variation is caused by the stellar spots of two active regions separated by about longitudinally located around the latitudes of and . In addition, 94 flares were detected and their parameters were computed. The OPEA model was derived for these flares and its parameters were computed. The value as saturation level for the active component was found to be 1.9510.069 s, while the value was found to be 1014 s. The flare frequency was found to be 0.01351 , while the flare frequency was found to be 0.00006. Maximum flare rise time () was found to be 1118.098 s, while maximum flare total time () was found to be 6767.72 s. Comparing its analogue it is seen that the chromospheric activity level of KIC 9761199, which is a low-mass close binary system according to the light curve analyses, is an expected level according to the () color index of for the active component.
Low-Mass Close Binary: KIC 9761199
E. YOLDAŞ1 and H. A. DAL1,2
Draft version:11footnotetext: Department of Astronomy and Space Sciences, University of Ege, Bornova, 35100 İzmir, Turkey22footnotetext: Corresponding Author, Email: email@example.com
Favor de proporcionar un resumen en español. If you are unable to translate your abstract into Spanish, the editors will do it for you.
Key Words: techniques: photometric — methods: data analysis — stars: binaries: eclipsing — stars: low-mass — stars: flare — stars: individual: KIC 9761199
It is well known for several decades that the population rate of red dwarfs in our Galaxy is about 65, while seventy-five percent of them show flare activity. The stars exhibiting flare activity are known as UV Ceti type stars (Rod86). As it is seen that the population rate of UV Ceti type stars in our Galaxy is about 48.75, which means that one of each two stars in our galaxy shows the flare phenomenon. As it is expected, the general population rate of UV Ceti type stars is incredibly high in both the open clusters and the associations (Mir90; Pig90). However, the population rate of these stars has a reducing trend, while the age of the cluster gets older due to the Skumanich’s law (Sku72; Pet91; Sta91; Mar92). This is because the mass loss rate is so high due to the flare activity during the main sequence life of the stars that the evolution track of the stars are changed by these mass loss in their main sequence stages, and also in later stages (Mar92).
According to the recent studies, the mass loss rate of UV Ceti type stars is about per year due to flare like events, while the solar mass loss rate is about per year (Ger05). The difference between the mass loss ratios is also seen between the flare energy levels of the solar and stellar cases. As it is generally observed in the case of RS CVn type active binaries (Hai91), the highest energy detected from the most powerful flares, known as two-ribbon flares, occurring on the sun is found to be - erg (Ger05; Ben08). On the other hand, the flare energy level varies from erg to erg in the case of dMe stars (Hai91; Ger05), while some stars of young clusters such as the Pleiades cluster and Orion association exhibit some powerful flare events, which energies reach erg (Ger83).
Although both the mass loss ratios and the flare energy levels are remarkably different between the solar and stellar cases, the flare events occurring on a dMe star are generally tried to explain by the classical theory of solar flare. The main energy source in the flare events is magnetic reconnection processes in this theory (Ger05; Hud97). However, this theory does not figure out each flare phenomenon observed from the UV Ceti type stars. To reach the real solution, examining the flare events occurring on the different type stars, all the differences and similarities should be demonstrated. Then, it should be identified which parameters such as singularity, binarity, mass, age, ect, cause these differences and similarities. For example, some differences are seen between the flare frequencies from one star to the next. The similar case occurs for the flare energy spectra.
In this study, we figured out the nature of a low mass eclipsing binary, KIC 9761199, which is a bit different from the classical UV Ceti type stars from spectral type dMe due to being a binary system. We did complete light curve analyses of the system for the first time in the literature in order to find out the physical properties of the components, using the PHOEBE V.0.32 software (Prs05), whose method depends on the 2003 version of Wilson-Devinney Code (Wil71; Wil90). Then, the variations out-of-eclipses were analyzed, and the flares occurring on the chromospherically active component were detected to model the magnetic activity nature of the system, comparing the active component with its analogue.
KIC 9761199 (KOI 1459) was listed with the B band brightness of in the USNO ACT Catalog (Urb97) in the first time literature. Secondly, KIC 9761199 was listed as J19083435+4630290 in 2MASS All-Sky Survey Catalogue, in which the JHK brightness of the system were given as , , (Kha01; Cut03). Although the orbital period of the system () was determined as day for the first time by Wat06, analyzing the data obtained in the Kepler Mission, Cou11 indicated that the orbital period of the system () is day. There is no any light curve analysis of the system in detail, though Bor11 computed the semi-major axis () of the system as 0.013 , while Cou11 calculated the inclination () of the system as and noted that there is no third light excess. However, they also noted that the system should be examined by a complete light curve analysis in detail.
In the Kepler Database, the () color index of the system is given as , while the temperature ratio of the components is given 0.682 by Sla11. In addition, the inclination () of the system is listed as . However, this inclination value is not in agreement with given by Cou11. According to Wal13, the age of the system is 0.77 Gyr and the () color index of the system is .
The one of the spectral studies of the system is belong to Man13 and they gave the luminosity () of the system as 0.082 . Using the calibration of some spectral observations, Mui12 indicated that the system is a binary with components from the spectral type of M1, whose distance is about 198 pc. The total masses given in the literature for the components vary from 0.51 (Mui14) to 0.65 (Cou12), while the radius vary from 0.48 (Mui14) to 0.84 (Cou12). The temperature given in the literature for the system varies from 3742 K (Mui14) to 4060 K (Cou12). Using the calibrations, Cou12 calculated mass and radius of each component as , , , .
2 Data and Analyses
The Kepler Mission is one of the space missions to be aimed to find out exo-planet, in which More than 150.000 targets have been already observed (Bor10; Koc10; Cal10). The quality and sensitivity of these observations have the highest one ever reached in the photometry (Jen10a; Jen10b). In these observations, so many variable targets such as new eclipsing binaries, etc. have been also discovered apart from the exo-planets (Sla11; Mat12). There are lots of single or double stars, which some of them are the eclipsing binaries, exhibiting chromospheric activity among these newly discoveries (Bal15).
The data analyzed in this study were taken from the Kepler Mission Database (Sla11; Mat12). After removing the all the observations with large error due to the technical problems from the data, using the ephemeris taken from the Kepler Mission database, the phases were computed for all data used in the analysis, and obtained light curves were shown in Figure 1. Because of this study format, the detrended short cadence data were used in the analysis instead of the long cadence. The data were arranged as the suitable formats to analysis the flare events and also the light curve.
2.1 Light Curve Analysis
To compute the flare parameters, it needs the synthetic light curve as assumed the quiescent levels. In this purpose, the light curves shown in the upper panel of Figure 1 were analyzed, after removing the instant-short term light variations due to the flare activity. The initial analyses indicated that there are several solutions according to these data. Because of this, the averages of all the detrended short cadence data were computed phase by phase with interval of 0.001. Then, this averaged light curve shown in Figure 2 was analyzed.
Using the PHOEBE V.0.32 software (Prs05), which is employed in the 2003 version of the Wilson-Devinney Code (Wil71; Wil90), we analyzed the light curves obtained from the averages of all the detrended short cadence data. We attempted to analyze the light curves in various modes, including the detached system mode (Mod2), semi-detached system with the primary component filling its Roche-Lobe mode (Mod4), and semi-detached system with the secondary component filling its Roche-Lobe mode (Mod5). Our initial analyses demonstrated that an astrophysically reasonable solution was obtainable only in the detached system mode; no results that were statistically consistent with reasonable solutions could be obtained in any of the other modes.
There are several studies about KIC 9761199 in the literature, and lots of temperature values were given for the system varying from 3742 K (Mui14) to 4060 K (Cou12). To understand which one is right, in this study, we calculated both () and () color indexes from the JHK brightness (, , ) given by Kha01 and Cut03. Then, using the calibrations determined by Tok00, we derived de-reddened colors as and . In addition, using the same calibrations, we derived the temperature of the primary component as 404020 K depending on these de-reddened colors. In fact, the initial analyses with different temperature values between 3742 K - 4060 K indicated that an astrophysically acceptable solution could be obtained, if the temperature would be taken 4060 K for the primary component. Because of this, in the analysis, the temperature of the primary component was fixed to 4060 K, while the temperature of the secondary component was taken as adjustable parameter. Considering the spectral type corresponding to this temperature, the albedos ( and ) and the gravity-darkening coefficients ( and ) of the components were adopted for the stars with the convective envelopes (Luc67; Ruc69). The nonlinear limb-darkening coefficients ( and ) of the components were taken from Van93. In the analyses, their dimensionless potentials ( and ), the fractional luminosity () of the primary component, the inclination () of the system, the mass ratio of the system (), and the semi-major axis () were taken as the adjustable free parameters.
The sinusoidal variation out-of-eclipses were modelled with two cool spots on the primary component. In the analysis, although the third light was taken as the adjustable free parameters, it is seen that there is no third light excess in the total light.
All the parameters obtained from the light curve analysis are listed in Table 1. As it can be seen from the table, the error values of each parameter are saliently small. This is caused due to the averaged data used in the analyse. If the original data plotted in the upper panel of Figure 1 were used in the analyse, it would be waiting that the error values are more larger. In addition, apart from some theoretical parameters such as the nonlinear limb-darkening coefficients ( and ), the albedos ( and ) and the gravity-darkening coefficients ( and ), the temperature of the primary component is just one parameter to fixed in the analysis. However, the temperature of the primary component was not also arbitrary determined. This temperature was determined from the JHK data of system given in the 2MASS All-Sky Survey Catalogue. As a result, the parameters found from the light curve analysis in Table 1 are absolutely trustworthy to figure out the obvious nature of KIC 9761199.
The synthetic light curve obtained with these parameters is shown in Figure 2. In this figure, both primary and secondary minima of the light curve are also plotted in the wide plane to better view for the readers. In addition, the 3D model of Roche geometry depending on these parameters is shown in Figure 3. The Roche geometry and the cool spot locates depending on the parameters obtained from the light curve analysis is shown for different phases, such as the phases of 0.00, 0.25, 0.50, 0.75.
Although there is not any available radial velocity curve, we tried to estimate the absolute parameters of the components. Considering the calibrations given by Tok00, () color index was found to be for the primary component depending on its temperature value derived from the analysis, while it was found to be for the secondary components. These color indexes are in agreement with the value found by Wal13. Using the calibrations given by Tok00, the mass of the primary component must be 0.570.01 corresponding to its () color index. Considering the possible mass ratio of the system, the mass of the secondary component was found to be 0.390.02 . Using Kepler’s third law, we calculated the possible semi-major axis as a 5.160.02 . Considering this estimated semi-major axis, the radius of the primary component was computed as 0.620.05 , while that of the secondary component was computed as 0.560.05 .
2.2 Flare Activity and the OPEA Model
To demonstrate the chromospherically active nature of the system, it needs to reveal the light variations just due to flare events. For this aim, first of all, all the light variations due to both the geometrical effects and the rotational modulation caused by the spots were removed from the general light variation. For this purpose, the data of all the pre-whitened light curves were obtained. To acquire the pre-whitened data, we extracted the synthetic light curves from all the detrended short cadence data.
In the second step, to determine the basic flare parameters such as the first point of the flare beginning, the last point of the flare end and the flare energy, the quiescent levels for each flare should be derived. In this point, the synthetic model lead us to definite the quiescent levels for each flare at the same time of that flare. Using the synthetic model as the quiescent levels, the parameters of the flares were computed. Two samples of the flare light curves taken from observation data and the quiescent levels derived for these flares are shown in Figures 4.
Using the synthetic models assumed as the quiescent levels, both the beginning and the end of a flare for each one were defined, and then, some flare parameters, such as flare rise times (), decay times (), amplitudes of flare maxima, flare equivalent durations (), were computed. In total, 94 flare were detected from the available short cadence data in the Kepler Mission Database. All the computed parameters are listed in Table 2 for these 94 flares. In the table, flare maximum times, equivalent durations, rise times, decay times and amplitudes of flare maxima are listed from the first column to the last, respectively.
In the analysis, the equivalent durations of the flares were computed using Equation (1) taken from Ger72:
where is the flux of the star in the observing band while in the quiet state. In this study, we computed using by the synthetic models derived with the light curve analysis. is the intensity observed at the moment of the flare. is the flare-equivalent duration in the observing band. In this study, the flare energies () were not computed to be used in the following analyses due to the reasons described in detail by Dal10; Dal11.
Examining the relationships of the flare parameters among each other, it is seen that the distributions of flare equivalent durations on the logarithmic scale versus flare total durations are varying according to a rule. The distributions of flare equivalent durations on the logarithmic scale cannot be higher than a specific value for the star, and it is no matter how long the flare total duration is. Using the SPSS V17.0 (Gre99) and GrahpPad Prism V5.02 (Daw04) programs, Dal10; Dal11 indicated that the best function is the One Phase Exponential Association (hereafter OPEA) for the distributions of flare equivalent durations on the logarithmic scale versus flare total durations. The OPEA function has a term, and this makes it a special function in the analyses. The OPEA function is defined by Equation (2):
where the parameter is the flare equivalent duration on a logarithmic scale, the parameter is the flare total duration, according to the definition of Dal10, and the parameter is the flare-equivalent duration in on a logarithmic scale for the least total duration. In other words, the parameter is the least equivalent duration occurring in a flare for a star. Here is an important point that the parameter does not depends on only flare mechanism occurring on the star, but also depends on the sensitivity of the optical system used for the observations. The parameter value is upper limit for the flare equivalent duration on a logarithmic scale. Dal11 defined value as a saturation level for a star in the observing band.
Using the least-squares method, the OPEA model was derived for the distributions of flare equivalent durations on the logarithmic scale versus flare total durations. The derived model is shown in Figure 5 together with the observed flare equivalent durations, while the parameters computed from the model are listed in Table 3. The value listed in the table is difference between and values. According to the definition of the OPEA function, the parameter in Equation (2) is a constant for just this model, depending on the x values. The value is half of the first value, at which the model reaches the value. In other words, it is half of the minimum flare total time, which is enough to the maximum flare energy occurring in the flare mechanism.
It was tested by using three different methods, such as the D’Agostino-Pearson normality test, the Shapiro-Wilk normality test and also the Kolmogorov-Smirnov test, given by Dag86 to understand whether there are any other functions to model the distributions of flare equivalent durations on the logarithmic scale versus flare total durations. In these tests, the probability value called as was found to be and this means that there is no other function to model the distributions of flare equivalent durations (Mot07; Spa87).
KIC 9761199 was observed as long as 289.82931 day from JD 2455641.50631 to JD 2455931.33562 without any remarkable interruptions. In total, significant 94 flares were detected in these observations. The total flare equivalent duration computed from all the flares was found to be 628.11671 s (0.17448 hours). Ish91 described two frequencies for the stellar flare activity. These frequencies are defined as given by Equations (3) and (4):
where is the total flare number detected in the observations, and is the total observing duration, while is the total equivalent duration obtained from all the flares. In this study, frequency was found to be 0.01351 , while frequency was found to be .
3 Results and Discussion
The results obtained from the analyses of short cadence data in the Kepler Mission Database indicated that KIC 9761199 exhibits high level chromospheric activity. However, it needs to identify which component is active, then, it should be compared with its analogue to be certain about activity level of the system. All these need to determine the physical parameters of the components, firstly.
Although all the variations apart from the geometrical effect, such as flare activity, were removed from the short cadence data taken from the database, neither the eclipses nor any other variation could be seen clearly in these eliminated data, as it is seen from the upper panel of Figure 1. To use in the light curve analysis, we computed the averages phase by phase with interval of 0.001 for the eliminated short cadence data. As it can be seen from Figure 2, both the primary and secondary minima and also a sinusoidal variation due to the rotational modulation are clearly visible in these averaged data. Here is an important point that the amplitude difference between the primary and secondary minima can be visible. This is important, because the orbital period of the system is a controversial point in the literature (Wat06; Cou11). However, this figure demonstrates that the orbital period of the system is 1.3839980 day, not 0.692031 day. According to this result, first of all, KIC 9761199 is a double object. There are two different minima, because of this, the secondary object can not be a planet, it must be a star. If it was a planet, there would not be any secondary minimum in the light curve due to the contrast effect between a star and a planet.
In the literature, although there are a few approaches using some calibrations for the physical parameters of the components, there is no complete light curve analysis for KIC 9761199. In this purpose, the light curve obtained from the detrended short cadence data was analyzed, using the PHOEBE V.0.32 software (Prs05), which is employed in the 2003 version of the Wilson-Devinney Code (Wil71; Wil90). Although the spectral type was given as M1V for KIC 9761199 (Mui12), there are several temperature values mentioned for this system. According to the de-reddened colors and of the system, the temperature of the primary component was taken 4060 K as given by Cou12. In the light curve analysis with the PHOEBE V.0.32 software, the temperature of the secondary component was found to be 38911 K. The mass ratio of the system () was found to be 0.690.01, while the inclination () of the system was found to be .
In addition, as it is seen from Figure 1 and the upper panel of Figure 2, the light curve exhibits a remarkable sinusoidal variation at out-of-eclipses. This variation could be caused by the effects of tidal distortion or the mutual heating of the components themselves. All these cases were included to the analysis. For this aim, the albedos ( and ), the gravity-darkening coefficients ( and ) of the components were computed depending on their temperatures and their dimensionless potentials ( and ) were also derived in the analysis. However, it was seen from the initial results that the synthetic light curve obtained with these parameters did not fit the observations. In this point, considering both the temperatures and the fractional radii of the components, and also considering the flare activity, we assumed that this sinusoidal variation is caused by the rotational modulation due to the stellar cool spots. Because of this, the sinusoidal variation was modelled with two cool spots on the primary component in the light curve analysis. The spots are separated by about longitudinally, while one of them is located about the latitudes of and the other one is located about the latitudes of . Here, it must be noted that it was assumed the cool spots locate on the primary component. However, it is likely they can be located on the secondary component, too. Someone can easily reach an acceptable solution, astrophysically. A strong clue for the answer of this quandary is come from the discussion of flare activity in the later part of the text. The 3D model of Roche geometry for the components and also the cool spots on the primary component is shown in Figure 3 for different phases.
Considering both the component temperatures and fractional radii found from the light curve analysis, using the calibrations given by Tok00, and also Kepler’s third law, we tried to estimate the absolute parameters of the components. The mass of the primary component found be 0.570.01 while it was found to be 0.390.02 for the secondary component. In addition, the radius of the primary component was computed as 0.620.05 , while that of the secondary component was computed as 0.560.05 .
When it is considered the results of the analysis, it is seen that both the locations of the stellar spots are not changed on the primary component along 289.82931 day (9.66 months), during all the observing season. However, this is a serious moot point for the astrophysical perspective, because a spotted area evolves in 2 or 3 months at most in the solar case (Ger05). On the other hand, according to Hal89 and Ger05, it is well known that the spotted areas on the active components of some RS CVn binaries, such as V478 Lyr, can keep their shapes and locations about two years. Therefore, the behaviour of the cool spot activity observed in the case of KIC 9761199 is not inconsistent event in respect to the stellar spot activity phenomenon. It must be noted that although we stated the word spot, actually it refers to an active area, in which several spots can appear and disappear in time. Here, the question to be answered is why the locations of these areas are stable for 289.82931 days. This is a discussion about the differential rotation of the components.
We detected 94 flare events from KIC 9761199, and calculated some parameters for each flare, such as flare frequencies. Yol16 recently resolved the chromospheric activity nature of a different system, FL Lyr. They found both flare frequencies as 0.41632 and 0.00027 for that system. Comparing the flare frequencies of KIC 9761199 with those of FL Lyr, it is clearly seen that FL Lyr exhibits flare more frequently than KIC 9761199, and also its flares are more powerful than those of KIC 9761199. Comparing KIC 9761199’s flare frequency values with those found for UV Ceti type flare stars from spectral type dK5e to dM6e, the flare energies obtained from KIC 9761199 are remarkably lower than those obtained from UV Ceti stars. For instance, the observed flare number per hour for AD Leo was found to be 1.331 , while it was found to be 1.056 for EV Lac. Moreover, frequency was found to be for EQ Peg, while it was found to be 0.086 for AD Leo (Dal11). As it is clearly seen from these results, the flare frequencies of KIC 9761199 are remarkably small. However, it is well known from Dal11 that the flare frequency can dramatically changes from one season to the next for some stars, such as V1005 Ori, EV Lac, etc. In this case, there could be some changes in the flare frequency and the flare behaviour of KIC 9761199 in the next observing seasons. In addition, these results found from the flare frequency analyses reveal why any flare had not been detected from this system by any ground based telescope before the Kepler Mission. If frequency is especially considering, it will be understandable, because, frequency indicates that the flare events occurring on the active component of the system are so weak to observe by any ground based telescope.
The value was found to be 1.9510.069 s from the OPEA model derived from the variation of the flare equivalent duration distributions on the logarithmic scale versus flare total durations for 94 flares. However, Yol16 found the value as 1.2320.069 s for FL Lyr. Moreover, Dal11 computed the values for some UV Ceti type stars. They gave the values as 3.014 s for EV Lac () and 2.935 s for EQ Peg (), and also it is 2.637 s for V1005 Ori (). As it is seen that the maximum flare energy detected from KIC 9761199 is remarkably smaller than the maximum energy level obtained from UV Ceti type single flare stars. The value of this system a bit approaches just to the value obtained from V1005 Ori. It must be noted that Dal11 found that the value is always constant for a star, while it is changing from one star to the other depending on their color indexes. The authors defined the value as the energy saturation level for the flare mechanism occurring on the star. As a result, the flare activity should be occurred on the primary component due to its color index. If it would be the secondary component, the value of KIC 9761199 was seen a bit incompatible according to the analogue of the secondary component.
Using the regression calculations, the value was found to be 1014.0 s from the OPEA model for KIC 9761199. In the case of FL Lyr, it is 2291.7 s (Yol16). It means that a flare occurring on the FL Lyr can reach the maximum energy level when the flare total duration reaches minutes, while it takes minutes for KIC 9761199. Here is a constant depending on the OPEA function of a star (Spa87; Daw04; Mot07). In other words, the parameter refers a minimum duration limit for a flare reached maximum energy. In this perspective, the flares, whose total times are shorter than minutes, can never reach the value obtained from the OPEA model derived for the flares of KIC 9761199. This value is a few times higher than those obtained from single dMe stars. Because, in the case of single dMe stars, for example, it was found to be 433.10 s for DO Cep (), and 334.30 s for EQ Peg, while it is 226.30 s for V1005 Ori (Dal11). It means that in the case of the stars such as EQ Peg, V1005 Ori and DO Cep, the flares can reach the maximum energy level at their value, when their total durations reach about minutes, while it needs minutes for KIC 9761199.
On the other hand, maximum flare rise time () obtained from the flares of eclipsing binary KIC 9761199 was found to be 1118.098 s, while maximum flare total time () was found to be 6767.72 s. However, these values are s and s for FL Lyr. As a result, FL Lyr flare time scales are larger than those obtained from KIC 9761199, which is in agreement with the results found by Dal11 for the single flare stars from dMe. In the light curve analysis, the primary component is assumed as the chromospherically active component and its color index was found to be . However, the () color index of the active component of FL Lyr is . Consequently, the chromospherically active component of KIC 9761199 is cooler than that of FL Lyr. In this case, the flare time scales of FL Lyr must be larger than the other according to Dal11.
Finally considering the value, flare frequencies and also flare time scales, it is clearly seen that the flare activity level of KIC 9761199 is clearly lower than almost all the UV Ceti type stars, as it is in the case of FL Lyr. This result about the chromospheric activity level of the system is also in agreement with the result found by Dal11. The author indicated that all these parameters derived from the OPEA model of a star are depend on just the () color index of that star. In other word, the OPEA model parameters increase or decrease according to variation of the () color index from a star to the other. As a matter of course, the OPEA model parameters derived for KIC 9761199 were found as expected values for the () color index of the primary component, not the secondary component.
In this point, it could be better to recapitulate that comparing the flare activity nature seen from KIC 9761199 indicated that the assumption of ”the chromospherically active star is the primary component of the system” is correct. In addition, the statistical analyses of the OPEA model demonstrated that there is only one star exhibiting the flare activity in the system. Because, according to the statistical analyses, the probability value found to be . It means that there is no other function to model the distributions of flare equivalent durations (Mot07; Spa87). Therefore, there is only one star exhibiting flare activity, in front of us. Considering the classical theory of solar flare (Ger05), the star exhibiting flare activity and cool spot activity must be the same star.
According to Dal11, it is a contentious issue which parameter, the value or the flare frequencies (), is the best indicator for the activity level. In general, the chromospheric activity level of a system depends on some parameters such as especially stellar rotation velocity. The rotation velocity depends on generally the stellar age or being a component of a close binary. Because of this, we discussed KIC 9761199 for both cases. Firstly, the age of KIC 9761199 is given as 0.77 Gyr (Wal13), while the age of FL Lyr is given between 3.05 Gyr and 15.25 Gyr in the literature. Considering Sku72’s law, it would be expected that the activity level of KIC 9761199 is clearly higher than that of FL Lyr. Secondly, in this study we computed the radii of the components as and , and the semi-major axis as 5.16 . However, the radii of the FL Lyr components are given as , , and its semi-major axis is given as by Eke14. According to these absolute parameters of the systems, the activity level of KIC 9761199 should be higher than FL Lyr activity level. If the flare frequencies () are considered, the flare frequencies of FL Lyr are higher than KIC 9761199. However, if the values are considered, the value of KIC 9761199 is clearly higher than that of FL Lyr. According to the semi-major axis of the KIC 9761199, the components are closer to the each other than the status of FL Lyr’s components. Maybe this is why the value is higher in the case of KIC 9761199. However, there is still an unsolved point yet. Why the flare frequencies are higher in the case of FL Lyr? The similar phenomenon is also common among the UV Ceti type single stars (Dal10; Dal11).
As a result, according to possible masses and radii of the components, KIC 9761199 must be a low-mass and a close binary system. In addition, according to the OPEA model parameters, just one component of the system is chromospherically active star. Both the age given in the literature and the proximity of the components can help to keep the activity level as possible as high. However, the active component exhibits intense flare activity, while it exhibits quiescent spot activity. It does not mean that the cool spot activity level is low in this system. It is possible that the active component could be covered by some large spots spread to the all surface. In this case, there would be no rapid light variation in the light curve out-of-eclipses contrary to the case of FL Lyr. To reveal and understand the nature of the cool spot activity in this system, it needs to track the sinusoidal light variation due to the rotational modulation in the light curve. However, these observations will be very difficult for the ground based telescope due to the sensitivity problems. In addition, there is no detailed spectroscopic observation of the KIC 9761199 in the literature. Therefore, it is also needed for the future studies.
The authors thank the referee for useful comments that have contributed to the improvement of the paper.
- Balona (2015) Balona, L.A., 2015, MNRAS, 447, 2714
- Benz (2008) Benz, A.O., 2008, Living Rev. Solar Phys., 5, 1
- Borucki et al. (2010) Borucki, Wi.J., Koch, D.G., Basri, G., et al., 2010, Sci, 327, 977
- Borucki et al. (2011) Borucki, Wi.J., Koch, D.G., Basri, G., et al., 2011, ApJ, 736, 19
- Caldwell et al. (2010) Caldwell, D.A., Kolodziejczak, J.J., & Van Cleve, J. E., 2010, ApJL, 713, L92
- Coughlin et al. (2011) Coughlin, J.L., López-Morales, M., Harrison, T.E., et al., 2011, AJ, 141, 78
- Coughlin & López-Morales (2012) Coughlin, J.L., López-Morales, M., 2012, AJ, 143, 39
- Cutri et al. (2003) Cutri, R.M., Skrutskie, M.F., van Dyk, S., et al., 2003. The IRSA 2MASS all sky point source catalog. NASA IPAC Infrared Science Archive http://irsa.ipac.caltech.edu/applications/Gator
- D’Agostino & Stephens (1986) D’Agostino, R B. & Stephens, M.A., 1986, ”Tests for Normal Distribution” in Goodness-Of-Fit Techniques, Statistics: Textbooks and Monographs, New York: Dekker, edited by D’Agostino, R.B., Stephens, M.A.
- Dal & Evren (2010) Dal, H.A. & Evren, S., 2010, AJ, 140, 483
- Dal & Evren (2011) Dal, H.A. & Evren, S., 2011, AJ, 141, 33
- Dal et al. (2012) Dal, H.A., Sipahi, E., Özdarcan, O., 2012, PASA, 29, 150
- Dawson & Trapp (2004) Dawson, B. & Trapp, R.G., 2004, ”Basic and Clinical Biostatistics” (New York: McGraw-Hill), 61
- Eker et al. (2014) Eker, Z., Bilir, S., Soydugan, F., et al., 2014, PASA, 31, 24
- Gershberg & Shakhovskaya (1983) Gershberg, R. E. & Shakhovskaya, N. I., 1983, Astrophys. Space Sci., 95, 235
- Gershberg et al. (1972) Gershberg, R.E., 1972, Astrophys. Space Sci., 19, 75
- Gershberg (2005) Gershberg, R.E., 2005, ”Solar-Type Activity in Main-Secauence Stars”, Springer Berlin Heidelberg, New York, p.53, 191, 360
- Green et al. (1999) Green, S B., Salkind, N.J. & Akey, T.M., 1999, ”Using SPSS for Windows: Analyzing and Understanding Data” (Upper Saddle River, NJ: Prentice Hall), 50
- Haisch et al. (1991) Haisch, B., Strong, K.T., Rodonó, M., 1991, ARA&A, 29, 275
- Hall et al. (1989) Hall D.S., Henry G.W., Sowell J.R., 1989. AJ, 99, 396.
- Hudson & Khan (1997) Hudson, H.S. & Khan, J.I., 1997, in ASP Conf. Ser. 111, ”Magnetic Reconnection in the Solar Atmosphere”, ed. R. D. Bentley & J. T. Mariska (San Francisco, CA: ASP), 135
- Ishida et al. (1991) Ishida, K., Ichimura, K., Shimizu, Y., et al., 1991, Ap&SS, 182, 227
- Jenkins et al. (2010a) Jenkins, J.M., Caldwell, D.A., Chandrasekaran, H., et al., 2010a, ApJL, 713, L87
- Jenkins et al. (2010b) Jenkins, J.M., Chandrasekaran, H., McCauliff, S.D., et al., 2010b, Proc. SPIE, 7740, 77400
- Kharchenko (2001) Kharchenko, N.V., 2001, KFNT, 17, 409
- Koch et al. (2010) Koch, D.G., Borucki, W.J., Basri, G., et al. 2010, ApJL, 713, L79
- Lucy (1967) Lucy, L.B., 1967, Z. Astrophys, 65, 89
- Mann et al. (2013) Mann, A.W., Gaidos, E., Ansdell, M., 2013, ApJ, 779, 188
- Marcy & Chen (1992) Marcy, G.W. & Chen, G.H., 1992, ApJ, 390, 550
- Matijevič et al. (2012) Matijevič, G., Prša, A., Orosz, J.A., et al. 2012, AJ, 143, 123
- Mirzoyan (1990) Mirzoyan, L.V., 1990, in IAU Symp. 137, ”Flare stars in star clusters, associations and the solar vicinity”, proceedings of the 137th IAU Symposium, Byurakan, Armenian SSR, Oct. 23-27, 1989 (A91-55760 24-90). Dordrecht, Netherlands, Kluwer Academic Publishers, 1990, p.1
- Motulsky (2007) Motulsky, H., 2007, ”GraphPad Prism 5: Statistics Guide”, San Diego, CA: GraphPad Software Inc. Press, 94
- Muirhead et al. (2014) Muirhead, P.S., Becker, J., Feiden, G.A., et al., 2014, ApJS, 213, 5
- Muirhead et al. (2012) Muirhead, P.S., Hamren, K., Schlawin, E., et al., 2012, ApJL, 750, 37
- Pettersen (1991) Pettersen, B.R., 1991, Mem. Soc. Astron. Ital., 62, 217
- Pigatto (1990) Pigatto, L., 1990, in IAU Symp. 137, ”Flare stars in star clusters, associations and the solar vicinity”, proceedings of the 137th IAU Symposium, Byurakan, Armenian SSR, Oct. 23-27, 1989 (A91-55760 24-90). Dordrecht, Netherlands, Kluwer Academic Publishers, 1990, p.117
- Prša & Zwitter (2005) Prša, A. & Zwitter, T., 2005, ApJ, 628, 426
- Rodonó (1986) Rodonó, M., 1986, NASSP, 492, 409
- Rucinski (1969) Rucinski, S.M., 1969, AcA, 19, 245
- Skumanich (1972) Skumanich, A., 1972, ApJ 171, 565.
- Slawson et al. (2011) Slawson, R., Prša, A., Welsh, W.F., et al. 2011, AJ, 142, 160
- Spanier & Oldham (1987) Spanier, J. & Oldham, K.B., 1987, ”An Atlas of Function”, Washington, DC: Hemisphere Publishing Corporation Press, 233
- Stauffer (1991) Stauffer, J.R., 1991, in Proc. ”NATO Advanced Research Workshop on Angular Momentum Evolution of Young Stars”, ed. S. Catalano & J.R. Stauffer (Dordrecht: Kluwer), 117
- Tokunaga (2000) Tokunaga, A.T. 2000, in Allen’s Astrophysical Quantities, ed. A.N. Cox (Springer), p.143
- Urban et al. (1997) Urban, S., Corbin, T., Wycoff, G., 1997, AAS, 191, 5707
- Van Hamme (1993) Van Hamme, W., 1993, AJ, 106, 2096
- Walkowicz & Basri (2013) Walkowicz, L.M. & Basri, G.S., 2013, MNRAS, 436, 1883
- Watson (2006) Watson, C.L., 2006, SASS, 25, 47
- Wilson (1990) Wilson, R.E., 1990, ApJ, 356, 613
- Wilson & Devinney (1971) Wilson, R.E., Devinney, E.J., 1971, ApJ, 166, 605
- Yoldaş & Dal (2016) Yoldaş, E. and Dal, H.A., 2016, PASA, 33, 16