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

    Prediction of Short Stellar Activity Cycles using Derived and Established Empirical Relations between Activity and Rotation Periods

    2024-01-06 06:40:38AlthukairandTsiklauri
    Research in Astronomy and Astrophysics 2023年11期

    A.K.Althukair and D.Tsiklauri

    1 Department of Physics and Astronomy, School of Physical and Chemical Sciences, Queen Mary University of London, Mile End Road, London, E1 4NS, UK;a.k.althukair@qmul.ac.uk, d.tsiklauri@qmul.ac.uk

    2 Physics Department, College of Sciences, Princess Nourah Bint Abdulrahman University, Riyadh, P.O.Box 84428, Saudi Arabia

    Abstract In our previous work, we investigated the occurrence rate of super-flares on various types of stars and their statistical properties, with a particular focus on G-type dwarfs, using entire Kepler data.The said study also considered how the statistics change with stellar rotation period,which in turn,had to be determined.Using such new data, as a by-product, we found 138 Kepler IDs of F- and G-type main sequence stars with rotation periods less than a day (Prot<1 day).On one hand, previous studies have revealed short activity cycles in F-type and G-type stars and the question investigated was whether or not short-term activity cycles are a common phenomenon in these stars.On the other hand, extensive studies exist which establish an empirical connection between a star’s activity cycle and rotation periods.In this study, we compile all available Kepler data with Prot<1 day,and rely on an established empirical relation between Pcyc and Prot with the aim to provide predictions for very short 5.09 ≤Pcyc ≤38.46 day cases in a tabular form.We propose an observation to measure Pcyc using a monitoring program of stellar activity (e.g., activity-related chromospheric emission S-index) or a similar means for the Kepler IDs found in this study in order put the derived empirical relations between Pcyc and Prot derived here to the test.We also propose an alternative method for measuring very short Pcyc, using flare-detection algorithms applied to future space mission data.

    Key words: stars: activity – stars: flare – stars: rotation – stars: solar-type – stars: statistics – Sun: flares

    1.Introduction

    The 11 yr cycle of solar activity, discovered by Schwabe in 1844(Schwabe 1844),is a significant phenomenon in solar and stellar physics.The cycle is manifested by a periodic change in solar activity, including the appearance of sunspots and changes in the Sun’s magnetic field on this timescale.Smoothed sunspot numbers have been widely regarded as a proxy for solar activity over the past four centuries (Shepherd et al.2014).The idea of the sunspot number was first introduced by Waldmeier(1961)in the mid-20th century,and it has since become a standard measure for quantifying solar activity.These numbers reveal that there are almost regular cycles of about 11 yr, reflecting the Sun’s magnetic activity.

    During the course of a solar cycle, the Sun experiences alternating periods of strong and weak activity known as solar maximum and minimum(Hathaway et al.2002;Shepherd et al.2014; Reinhold et al.2017).As the solar cycle progresses, the magnetic field becomes more complex and twisted.This results in the emergence of sunspots, which are dark areas on the surface of the Sun with intense magnetic fields,and which vary in size and can last from days to several months (Petrovay &van Driel-Gesztelyi 1997), decaying into bright areas called faculae formed by smaller magnetic concentrations (Reinhold et al.2017).During the active phase of the solar cycle (solar maximum), the sunspot number increases, and their size becomes larger on the surface of the Sun.At the same time,bright faculae also become more prominent.As the cycle progresses, the number of sunspots decreases, the overall brightness of the Sun reduces and the Sun enters its least active phase of the solar cycle (solar minimum).These dark and bright features on the Sun’s surface contribute to variability in the total solar irradiance (TSI) (Marchenko et al.2022).Therefore,the TSI data can capture the combined effects of the evolving dark and bright features during the solar cycle(Domingo et al.2009; Reinhold et al.2017).

    Cyclic activity has been observed in stars other than the Sun through long-term brightness changes associated with increased occurrence of active regions on their surfaces or in their lower stellar atmospheres (Reinhold et al.2017).The Mount Wilson HK Project, which started in 1966 and lasted until the end of the 20th century, was the first to conduct a systematic search for activity cycles in main sequence stars(Wilson 1978; Baliunas et al.1995; Mittag et al.2019a).By examining chromospheric emission within the Ca II H&K spectral lines,the magnetic field associated with active regions on stellar surfaces is vital in conveying energy to the chromosphere.This heightened influx of energy into the chromosphere results in amplified chromospheric emission,notably observable in the central regions of the Ca II H&K spectral lines, as indicated by Reinhold et al.(2017).The measure of the chromospheric emission strength is described by the Mount Wilson S-index (Vaughan et al.1978) or by the quantityRH′K(Brandenburg et al.2017).Vaughan & Preston(1980) investigated the chromospheric activity levels in main sequence F-G-K-M stars by measuring the chromospheric Ca II H&K emission fluxes.They noted that these stars display varying degrees of chromospheric activity and observed a noticeable lack in the number of F-G stars displaying intermediate activity compared to both highly active and less active stars.They suggested that the absence of such stars could be attributed to a decline in chromospheric activity as the stars age.Noyes et al.(1984a) examined the relationship between chromospheric activity,specifically theRH′Kactivity index,and the Rossby number Ro=Prot/τcfor a sample of main sequence stars of spectral type F or later, where Protis the rotational period of the star and τcis a theoretically derived convective turnover time.They found a strong correlation between theRH′Kactivity index and the Rossby number.However,in contrast to the findings of Vaughan&Preston(1980),Noyes et al.(1984a)did not find any signs of the “Vaughan–Preston gap.” Noyes et al.(1984b) investigated the empirical relation between rotation period Prot,spectral type and activity cycle period Pcycfor 13 slowly rotating main-sequence stars.They found that the cycle period is related to the rotation period by a power law:(Brandenburg et al.2017; Mittag et al.2023).For stars of spectral type G0-K5, Baliunas et al.(1995)observed a pattern of variation in the rotation period and the measure of chromospheric activity(S-index).Their research revealed that the chromospheric activity levels were high in young stars with fast rotation periods.Chromospheric activity and rotation rates of stars in the intermediate age range were average.Alternatively, the chromospheric activity levels were low in old stars with slow rotation periods.This observation supports the existence of the Vaughan–Preston gap, indicating that chromospheric activity and rotation change over time as the stars age.The relation between rotation periods and activity cycles of a sample of stars was investigated by Baliunas et al.(1996), who discovered a correlation between the two variables.In particular, they observed that stars with slower rotation periods exhibit longer activity cycles, while stars with faster rotation periods tend to have shorter activity cycles.According to Oláh & Strassmeier (2002), the relation between rotation periods and cycle lengths is more evident for stars with shorter activity cycles.However, the association becomes less clear for longer cycle lengths when considering more recent findings on the time variability of solar cycles.This relationship can alternatively be expressed as

    In order to provide a background of results in this field of research,we now discuss previous literature.Vida et al.(2013)investigated the behavior and activity cycles of four fastrotating late-type stars with (Prot≤0.5 day), highlighting the presence of 1 yr cycles and the correlation between rotation rate and cycle length.Vida et al.(2014)used the short-term Fourier transform, a time-frequency analysis method, to examine the light curves of 39 fast-rotating late-type active stars with rotation periods of less than one day.Nine of the selected stars showed indications of activity cycles with periods between 300 and 900 days.These cycles were inferred based on the observed variations in the typical latitude of the starspots.These variations, along with the differential rotation of the stellar surface,result in changes in the observed rotation period during the activity cycle.This variation in the rotation period was attributed to the movement and evolution of starspots at different latitudes of the star.

    Reinhold et al.(2017) used four years of Kepler data to determine the cyclic variations in the amplitude of the light curve and the rotation period of stars by analyzing a sample of active stars and calculating the rotation period and variability amplitude for each star in each Kepler quarter.Then they searched for periodic variations in these time series using Lomb–Scargle periodograms and employed a false alarm probability (FAP) criterion for selection.The study’s findings indicate that amplitude periodicities, associated with underlying activity cycles, are detected in 3203 stars with cycle periods ranging from 0.5 to 6 yr and rotation periods ranging from 1 to 40 days.According to Brandenburg et al.(2017)’s analysis of new observations and previous data,the longer and shorter cycle periods closely match expectations based on the average activity levels and rotation periods, which indicate a connection between stellar activity and stellar rotation.

    Baliunas et al.(1995)reported an activity cycle of 11.6 yr in the F-type star τ Boo (HD 120136).However, the authors assigned an FAP “poor” grade to this finding.Mittag et al.(2017b) detected an activity cycle with a duration of 122 days in their analysis of the S-index data of τ Boo.This short activity cycle period suggests that τ Boo may exhibit variations on a relatively short timescale.Mittag et al.(2019a)focused on exploring the presence of short-term activity cycles in F-type stars, specifically using S-index time series data obtained with the Telescopio Internacional de Guanajuato Robótico Espectroscópico (TIGRE) facility.They utilized the generalized Lomb–Scargle periodogram method to analyze the data and search for periodic variations with a maximum length of 2 yr.Their sample of F-type stars identified four stars that exhibited cyclic variations with periods of less than a year.However,compared to solar-type stars with well-developed cyclic activity, the amplitudes of these short-term cyclic variations in F-type stars were smaller.Based on their findings, Mittag et al.(2019a) concluded that the activity behavior among F-type stars differs from that of the Sun and cooler main sequence stars.

    By studying 44 main sequence stars with confirmed activity cycles and rotation periods, Mittag et al.(2023) examined the relation between the length of the activity cycle and the Rossby number (Ro).They used empirical turnover periods based on the B ?V color index to calculate Rossby numbers,from which they deduced an empirical relationship between the Rossby number and the cycle duration.The study reported linear behavior in the double-logarithmic relationship between the Rossby number and cycle period.In addition, the relative convection zone depth was found to be correlated with cycle length and convective turnover time.

    Besides the 11 yr solar cycle,shorter cycles were discovered called Rieger cycles.The original Rieger cycles were first identified in the Sun by Rieger et al.(1984) with a specific periodicity of approximately 154 days for flare occurrences.The Rieger-type cycles(RTCs)encompass cycles with periods(PRTC) ranging from 109 to 276 days.These RTCs were observed in various phenomena beyond solar flares, such as solar magnetic field and sunspot indexes, indicating their widespread nature.The underlying nature of RTCs remains unclear.The RTCs become more pronounced during the solar activity maximum.There is a potential connection between RTCs and the modulation of the solar magnetic dynamo process,as discussed in Arkhypov&Khodachenko(2021)and references therein.Possible reasons encompass the role of inertial g- and r-waves, also known as Rossby waves, as modulators of the emergence of magnetic flux in the Sun.Arkhypov & Khodachenko (2021) analyze photometric data from 1726 main sequence stars with varying effective temperatures and rotation periods to study RTCs in other stars.Two types of RTCs are identified among the surveyed stars.The activity cycles with RTC periods (PRTC) are independent of the stellar rotation period and are suggested to be driven by Kelvin waves.The second type are activity cycles with PRTCproportional to the stellar rotation period and are suggested to be driven by Rossby waves.

    The Parker (1955) model of the α ?Ω dynamo introduced the concept of migratory dynamo waves, which play a crucial role in generating the observed solar cycle(Mittag et al.2023).The α-effect, arising from the twisting of rising magnetic field tubes due to Coriolis forces,creates the poloidal magnetic field required for the next sunspot cycle.This effect is responsible for the reversal of magnetic polarities between successive cycles(Parker 1955;Mittag et al.2023).On the other hand,the Ω-effect, resulting from the differential rotation of the star,generates a toroidal magnetic field by stretching the magnetic field lines in a longitudinal direction.The combination of the α-effect and the Ω-effect leads to the formation of migratory dynamo waves, where the toroidal field is periodically regenerated and transformed into the poloidal field through the action of the α-effect.These migratory dynamo waves propagate and interact within the star’s convective zone,causing the cyclic variations in the magnetic field(Mittag et al.2023).

    Now we describe existing theoretical knowledge about the possible relation between the magnetic cycle period and the rotation period of a star.In this context, according to Noyes et al.(1984b), the magnetic cycle period for G and K dwarfs,with convective turnover times(τc)between 11 and 26 days,is found to be proportional to the rotation period as follows:

    where n is 1.25.

    Simple dynamo models were discussed for understanding stellar magnetic activity and their implications for magnetic cycle periods in stars.Stix (1981) derived an equation to determine the critical dynamo number Dcritgiven by

    Parker (1955) provides a relation for the magnetic cycle frequency, ωmag_cyc, that involves the shear, H, and the αeffect.Stix (1976) presented an equivalent expression for the magnetic cycle frequency derived by Parker(1955),in terms of angular velocity gradient,Ω′, given by

    indicating a proportional relationship between the cycle frequency and rotation frequency.Based on the model’s assumptions, Equation (3) can be written as

    where l here is the length scale of turbulence and R?is the stellar radius.This equation indicates the theoretical prediction of the relation between the star’s activity cycle and its rotation period, which is Equation (6) in Mittag et al.(2023).

    According to the simple theoretical arguments quoted by Mittag et al.(2023), the magnetic cycle period Pmag_cycis proportional to the rotation period Prot.However, there is a modifying factor, l/R?the relative depth of turbulence, which depends on the stellar structure,which itself may depend on the effective temperature or B ?V color index of the star.This factor is expected to vary among different stars, especially those with different sizes, masses and ages.The smallness of the inverse relative depth of the turbulence ensures that the period of the magnetic activity cycle Pmag_cycis small.However, precisely what factors guarantee smallness of R?/l in a particular star is poorly understood.That is why it is unclear why stars with very short activity cycles,studied in thispaper, exist.All we can surmise is that the above theoretical arguments suggest Pmag_cycshould scale as

    Table 1 List of Star IDs with their Parameters, used in Previous Studies

    Table 1(Continued)

    Activity cycles, characterized by variations in magnetic activity over time, are essential for understanding the fundamental mechanisms that drive the magnetic fields of stars.A range of methodologies exist for the identification of activity cycles in stars.One such approach involves integrated flux measurements by continuously monitoring the total amount of energy emitted by a star, enabling the detection of variations in its magnetic activity (Kopp et al.2016; Reinhold et al.2020).Another approach is the analysis of chromospheric emission lines from the outer atmosphere of a star.Additionally, tracking of starspots by observing the movement and changes in starspots on a star’s surface serves as an indicator for fluctuations in magnetic activity (Montet et al.2017).Nevertheless, it is important to acknowledge that these techniques include limitations regarding photometric precision and the small sample size in spectroscopic observations(Scoggins et al.2019).An alternative approach involves detecting flares.Flares are a frequently observed phenomenon resulting from magnetic activity and are easier to detect,even at significant distances from stars.Wide-field photometric surveys allow for simultaneous monitoring of stars, making it possible to survey them for flare activity(Scoggins et al.2019).The Sun experiences variations in its flare rate by a factor of 10 between the solar maximum and minimum activity periods.Scoggins et al.(2019) focused on fluctuations in the frequency of flares from stars detected by the Kepler mission.The study examined a sample of 347 flare stars, which were selected based on having measured Kepler rotation periods, with a minimum of 100 candidate flare events.Scoggins et al.(2019) aimed to identify coherent variations in flare activity among these stars by computing the fractional luminosity emitted in flares.One star,KIC 8507979,was identified as the best candidate for flare activity variation.This star has a rotation period of 1.2 days and emits an average of 0.82 flares per day with energies exceeding 1032erg over the 18 Kepler quarters.The study observed a decline in flare activity from KIC 8507979 over time.Although the flare census derived from the Kepler light curve of this star did not provide definitive evidence for a stellar activity cycle,the observed variation of approximately 0.1 dex per year was consistent with cyclic behavior over ten years or more.

    The motivation for the current work is as follows:In Paper I(Althukair & Tsiklauri 2023a), we looked for super-flares on different types of stars and focused on G-type dwarfs using the entire Kepler data to study various aspects of statistical properties of the occurrence rate of super-flares.In Paper II(Althukair & Tsiklauri 2023b), as a by-product, we found 13 peculiar Kepler IDs that are Sun-like, slowly rotating cases with rotation periods of 24.5–44 days, and yet can produce a super-flare and six G-type and four M-type Kepler IDs with exceptionally large amplitude super-flares.As noted previously, these detections defy our current understanding of stars and hence deserve a further investigation.In this Paper III,the last in this series,we use the same data set as in Althukair&Tsiklauri (2023a) in order to study the empirical connection between a star’s activity cycle and rotation period for a sample of F and G main sequence stars with rotation periods of less than one day.Here our aim is to provide predictions for very short activity cycle cases in a tabular form and to investigate in the future whether these short activity cycles are a common phenomenon in these stars or not.Section 2 presents the method used in this work which includes a reproduction of the Mittag et al.(2023) fit, the data representation and fit and the target selection method.The main findings of the study are presented in Section 3,and Section 4 concludes this work with our main conclusions.

    2.Methods

    In our study,we adopt the terminology used by Brandenburg et al.(2017), Mittag et al.(2023) to categorize branches into two types: the “inactive” branch, referred to as the short-cycle branchPScyc, and the “active” branch, referred to as the longcycle branchPLcyc.These terms were introduced for the first time in Brandenburg et al.(2017).According to Mittag et al.(2023), this notation is more accurate and aligned with the actual characteristics of the branches.Therefore,they suggested that these terms should be used in future studies to refer to the two branches.

    Figure 1.Log-scale of rotation period vs.log-scale of observed activity cycle period(short cycle branch)for a sample of stars taken from Mittag et al.(2023).The deduced fits of the Prot vs.Pcyc relation are shown as solid lines.The blue line displays the fit when slope n is treated as an independent parameter while the red line delineates the fit with a fixed slope of n=1.

    2.1.Reproduction of Mittag et al.(2023) PScyc versus Prot Fit

    In this subsection,we reproduce the fit betweenPScycand Protdata from Mittag et al.(2023)to derive the fit parameters.First,we collected the data in Table 1, the first 32 rows are the observed activity cycle on the short-cycle branchPScycfrom Mittag et al.(2023), Table 1, along with the 32 corresponding rotation periods Prot.These cycle lengths and rotation periods can be found in Table 1.Then we plotted,in logarithmic scale,the rotation periods on the x-axis versus the observed cycle period on the y-axis as shown in Figure 1, using the empirical relation in Mittag et al.(2023) between the cycle periods and rotation periods in logarithmic terms that is given by

    Since the theoretical relation, Equation (4), implies a linear connection between Pcycand Prot, we fitted the data using a Python least-squares fit, a common technique for determining the best-fitting parameters for a given model, for two different slope adjustments as in Mittag et al.(2023).Also,we computed the R2coefficient of determination to measure how well the model fits the data.An R2value of 1 means that the predictions from the regression fit the data perfectly.First,we set the slope n to be 1 and deduced the value of a parameter as a=1.918±0.027 and the value of R2=0.87.The red line in Figure 1 illustrates this trend.Then we repeated the fit by treating slope n as an independent variable to derive a and n values as Equation (5) now becomes

    and the value of R2=0.93.The fit given by Equation (6) is identical to that of Mittag et al.2023.The blue line in Figure 1 represents this fit.It is obvious that the n=1 relation does not fit the short period data, as Mittag et al.(2023) pointed out.

    2.2.Data Representation and Fit

    In this subsection, we repeat the fit between ProtandPScycusing a larger data sample,taken from Mittag et al.(2023)and other previous studies.

    This sample,shown in Table 1,contains 92 Protand their 92 correspondingPScyc.In addition to the 32 observed activity cycles in Mittag et al.(2023),we aggregated 32 activity cycles on the short-cycle branchPScyccomputed by Mittag et al.(2023)together with the corresponding 32 rotation periods Prot.Furthermore, we included 28 activity cycles and their corresponding rotation periods that were collected from various other studies.ThesePScycwere taken from Ferreira Lopes et al.(2015), Egeland et al.(2015), Boro Saikia et al.(2016), Flores et al.(2016), Brandenburg et al.(2017), Mittag et al.(2019a).The star ID, effective temperature (Teff), color index (B ?V),Protand Pcycare shown in Table 1.Unavailable data are left blank in the table.It should be noted that we used in the fit twoPScycvalues for each of the 32 star IDs in Mittag et al.(2023),one was the observedPScycby a previous study, and the other was the calculatedPScycby Mittag et al.(2023), except for HD 16673 for which we collected threePScycdue to the multiple sources,as shown in Table 1.References for each ProtandPScycare shown in Table 1.

    In the same way as in Section 2.1, we utilized the empirical relation between Protand Pcycin logarithmic scale given by Equation(5)using the new data set in Table 1 to produce the fit parameters a and n.We performed a least-squares fit in Python to fit the data using two different slope adjustments again, one with a fixed slope n of 1 and another with the n treated as a free variable.This fit is depicted in Figure 2.For the fit with a fixed slope of 1, we determined the value for the parameter a=1.890±0.024 and R2=0.83.This trend is visualized by the red line in Figure 2.While for the fit with the slope n treated as a free variable,we deduced values for the parameters a and n as a=1.585±0.064, n=1.256±0.051 and R2=0.87.This fit is represented by the blue line in Figure 2, so that Equation (5) now becomes

    We note that our value of n=1.256±0.051 with the extended data set is closer to Noyes et al.(1984b)?s n=1.25 than Mittag et al.(2023)?s n=1.324±0.067.

    Figure 2.Log-scale of rotation period vs.log-scale of cycle period(short cycle branch) for a sample of 92 stars taken from previous studies in Table 1.The deduced fit of the Prot vs.Pcyc relation is displayed as solid lines.The blue line indicates the fit where slope n is treated as an independent parameter while the red line shows the fit with a fixed slope of n=1.

    2.3.Data Samples

    One of the main challenges in studying the relation between cycle length and rotation period is the lack of well-known and accurately measured activity cycles.This limitation introduces uncertainties in the derived empirical relations (Mittag et al.2023).To overcome these challenges, it is crucial to obtain more reliable cycle periods,particularly for long-period cycles.Achieving this requires long-term time series observations of stars to gather comprehensive and accurate data on their activity cycles (Mittag et al.2023).Therefore, when looking for activity cycles, it is more efficient to monitor fast-rotating objects, as cycles can be discovered within a few years of observation, as opposed to stars with longer rotation periods(Vida et al.2013).For this reason,we chose our sample for this study to include fast-rotating main sequence stars of type F and G from Kepler data with well-known rotation periods of less than one day.First, we collected all Kepler IDs which have well-known rotation periods.We then selected targets with rotation periods of less than a day.Using Gaia Data Release 2(Gaia DR2), we identified F- and G-type main sequence stars by their effective temperatures and radii based on the Harvard Spectral classification.The ranges of the effective temperature are 6000–7500 K and 5200–6000 K for F- and G-types,respectively.We thus obtained a total of 811 Kepler IDs of F- and G-type stars with less than one day rotation period.By applying the radius restriction of 1.15–1.4 R⊙a(bǔ)nd 0.96–1.15 R⊙to main sequence stars for F- and G-types,respectively, the final data sample reduced to 138 Kepler targets with a number of 83 F-type and 55 G-type mainsequence stars; 71.74% of the rotation periods for these stars were taken from McQuillan et al.(2014), 15.94% from Santos et al.(2021), 5.07% from Reinhold & Gizon (2015), 4.35%from Chowdhury et al.(2018) and 2.90% from Yang & Liu(2019).These 138 Kepler targets are listed in Table 2 with their effective temperature,radius,rotation period and the references for these rotation periods.

    Table 2 Lists of the 138 Kepler IDs with their Parameters and Predicted Pcyc

    Table 2(Continued)

    Figure 3.Prot vs.Pcyc using a simple linear regression without an intercept for a sample of stars whose Pcyc is less than 1 year.The determined fit of the Prot vs.Pcyc relation is shown as a solid green line.

    3.Results

    Using a data set of 138 Kepler IDs with Protranging from 0.202 to 0.997 day, we provide a prediction for the corresponding value of theirPScyc, by applying the empirical relation between Pcycand Protwith the derived parameters in Equation (7).Hence we obtained the predicted values of Pcycfrom

    From Equation(8),we calculated 138 Pcycfor 83 F-type and 55 G-type main sequence stars whose rotation period is less than a day.The shortest Pcycis equal to 5.09 days while the longest Pcycis equal to 38.46 days.All the 138 predicted Pcycare listed in Table 2.

    After predicting the values of the activity cycles for our extended, compared to Mittag et al.(2023), data sample, we wish to examine the theoretical prediction given by Equation (4) on short Pcyc<1 yr.This is because the latter equation is a theoretical prediction, based on first physical principles, as opposed to an empirical fit, which lacks any theoretical or conceptual justification.Therefore, we focus on the activity cycles derived from previous studies, as presented in Table 1.We choose 20 stars whose Pcycis less than a year and plot the fit between Protand Pcycas depicted in Figure 3 using a simple linear regression without an intercept given by

    We obtained the slope n=0.081±0.009 and the R2value is 0.997.Although the R2value for the fit is near unity,there is a large scatter indicating a poor quality fit.Note that Pcychere is in years,as in Figure 14 from Mittag et al.(2019a).Therefore,for the lower and upper bounds of our 138 Kepler IDs with Protranging from 0.202 day to 0.997 day, this simple theoretically justified equation predicts Pcyc=0.081×0.202×365.25=5.98 days and 0.081×0.997×365.25=29.50 days, which are not very different from applying the more accurate power law fit using Equation (8) of 5.09 days and 38.46 days,respectively.

    4.Conclusions

    In this work, we studied the empirical relation between star activity cycle and rotation period.First, we reproduced the fit between Protand Pcycusing Mittag et al.(2023) data and obtained the following fit parameters logPcyc≈(1.488 ±0.092) +(1.324 ± 0.067) logProt, which are the same parameters as in Mittag et al.(2023).Then, using a larger data set made up of Protand their associated Pcyctaken from prior studies, we again re-examined the fit between Protand Pcycand obtained the following fit parameters logPcyc≈(1.585 ±0.064) +(1.256 ± 0.051) logProt.Using these new parameters, we applied this relation to a sample of 83 F-type and 55 G-type main sequence stars with rotation periods of less than one day.The objective was to predict short activity cycles for these stars, aiming to ascertain, in future studies, if short activity cycles are a common occurrence in these stars or not.As a result, we derived 138 predicted Pcycranging from 5.09 to 38.46 days, which are listed in Table 2.

    The usefulness of measuring short stellar activity cycles hinges on two main general difficulties:

    (i)If the monitoring program of stellar activity(e.g.,activityrelated chromospheric emission S-index or similar) is used as in references such as Mittag et al.(2019a); or Baum et al.(2022), then cadence time of observations is too long, e.g.,according to Table 2 from the latter reference cadence could be 87 observations per year, i.e., 365/87=4 days.Resolving activity cycles with 5.09 ≤Pcyc≤38.46 days with such cadence would be nearly impossible.

    (ii) If Kepler data light curves are used for, e.g., plotting number of flares per day versus time then a large number of flare detections would be necessary to have reliable statistics.However, the problem is long cadence, 30 minutes, for the mainstream Kepler data.The photometer used by Kepler is sensitive to wavelengths ranging from 400 to 865 nm,covering the entire visible spectrum and a fraction of the infrared.The accuracy of the photometer on Kepler is approximately 0.01%or 0.1 mmag, when 30 minute integration times are used while considering stars with a magnitude of 12.Kepler’s 30 minute integration detected flare amplitudes are less than 0.1% of the stellar value and have energies of 2×1033erg.The duration of the flares ranged from one to three hours,with a rapid increase followed by a slow,exponential decline(Maehara et al.2012).When Kepler data are taken at a higher cadence or sampling rate of one minute, the accuracy of the measurements decreases.However, this higher cadence enables Kepler to detect flares that are too brief to be detected reliably using the main 30 minute integrations.With the one-minute cadence,Kepler can detect flares with energies as low as 1032erg(Maehara et al.2015).

    It is worth noting that earlier studies exist using different observations where the energy involved in the observed transient brightening is estimated to range from 1025to 1029erg (Shimizu 1995).Also, as far as the Sun is concerned,studies exist (Mason et al.2023) which consider flare frequency as a function of flare energy in the range 1027–1031erg, but this is only applicable to the Sun.

    In order to have good statistics for the Kepler IDs considered, we need to detect flares with energies 1027–32erg in order to see a variation in the number of flares per day on a timescale of 5.09 ≤Pcyc≤38.46 days.To achieve this goal, a new space mission is necessary with short time cadence(<1 minute) and photometric accuracy <0.01%.

    A typical example of such a proposed hypothetical space mission would record data on the number of flares per day for each target.These data can be presented in bins of, e.g., oneday width where the bin heights would show the number of flares detected in that bin.These bins would then exhibit a periodic variation over time.Fitting a sinusoidal curve then would enable deducing the activity cycle period.Thus,through this periodic variation, we could potentially detect the target’s magnetic activity cycle period.In some sense, our approach is similar to that of Scoggins et al.(2019).However, their observation was so short in duration that only decrease in the flare activity was seen.A longer duration of observations from a proposed new space mission would enable seeing periodic variation and hence deducing the activity cycle period.

    An alternative option could be conducting a shorter cadence ground-based S-index monitoring program of stellar activity with cadence ≈1 day or less.However it is unclear whether this is technically feasible.In any case, the present study provides predictions for 5.09 ≤Pcyc≤38.46 days and we hope that either future space or ground-based observational missions will put our predictions to the test.Until such time, the jury is still out.

    Acknowledgments

    Some of the data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST).STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555.Support for MAST for non-HST data is provided by the NASA Office of Space Science via grant NNX13AC07G and by other grants and contracts.

    The authors would like to thank Deborah Kenny of STScI for kind assistance in obtaining the data,and Cozmin Timis and Alex Owen of Queen Mary University of London for their assistance in data handling at the Astronomy Unit.

    A.K.Althukair wishes to thank Princess Nourah Bint Abdulrahman University, Riyadh, Saudi Arabia and Royal Embassy of Saudi Arabia Cultural Bureau in London, UK for the financial support of her PhD scholarship, held at Queen Mary University of London.

    Authors would like to thank an anonymous referee whose comments greatly improved this manuscript.

    Data Availability

    Some of the data underlying this article were accessed from the Mikulski Archive for Space Telescopes (MAST)https://mast.stsci.edu/portal/Mashup/Clients/Mast/Portal.html.This paper also has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium).Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.The derived data generated in this research will be shared on reasonable request to the corresponding author.

    ORCID iDs

    A.K.Althukair https://orcid.org/0000-0003-4075-4440

    D.Tsiklauri https://orcid.org/0000-0001-9180-4773

    欧美色视频一区免费| 性插视频无遮挡在线免费观看| 青春草国产在线视频 | 亚洲国产精品sss在线观看| 成人毛片60女人毛片免费| 国产伦理片在线播放av一区 | 精品人妻熟女av久视频| 91久久精品国产一区二区成人| а√天堂www在线а√下载| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 黄色配什么色好看| 久久久精品大字幕| 尾随美女入室| 国产综合懂色| 一级黄片播放器| 久久久久性生活片| 国产v大片淫在线免费观看| 午夜精品国产一区二区电影 | 给我免费播放毛片高清在线观看| 国产淫片久久久久久久久| 三级毛片av免费| 国产一级毛片在线| 日本一本二区三区精品| 久久久久久久久久成人| 熟妇人妻久久中文字幕3abv| 国产伦一二天堂av在线观看| 日韩欧美 国产精品| 亚洲国产欧洲综合997久久,| 国产真实乱freesex| 乱码一卡2卡4卡精品| 免费无遮挡裸体视频| 国产午夜精品久久久久久一区二区三区| 色综合站精品国产| 久久久久网色| 精品久久久久久久久久免费视频| 99久久久亚洲精品蜜臀av| 日韩欧美三级三区| 麻豆国产av国片精品| 男女那种视频在线观看| 91久久精品电影网| 淫秽高清视频在线观看| 青春草国产在线视频 | 黑人高潮一二区| 亚洲精品久久国产高清桃花| 欧美一级a爱片免费观看看| 熟女电影av网| 久久久久久久久大av| 一个人观看的视频www高清免费观看| 国产精品伦人一区二区| 岛国毛片在线播放| 国产精品久久久久久精品电影| 天天躁日日操中文字幕| 日本在线视频免费播放| 亚洲欧美精品自产自拍| 日韩欧美三级三区| 亚洲人与动物交配视频| 国产成人午夜福利电影在线观看| 日韩一本色道免费dvd| 2021天堂中文幕一二区在线观| 国产熟女欧美一区二区| 亚洲成人精品中文字幕电影| 国产成人精品久久久久久| .国产精品久久| 你懂的网址亚洲精品在线观看 | 干丝袜人妻中文字幕| 蜜桃亚洲精品一区二区三区| 国产色爽女视频免费观看| 别揉我奶头 嗯啊视频| 蜜桃亚洲精品一区二区三区| 精品99又大又爽又粗少妇毛片| 婷婷色av中文字幕| 色视频www国产| 嘟嘟电影网在线观看| 麻豆av噜噜一区二区三区| 亚洲精品乱码久久久v下载方式| 一个人看视频在线观看www免费| 久久久午夜欧美精品| 国产成人精品婷婷| 亚洲四区av| 精品人妻一区二区三区麻豆| 欧美最新免费一区二区三区| www.色视频.com| 91在线精品国自产拍蜜月| 夜夜夜夜夜久久久久| 亚洲人成网站在线观看播放| 久久精品人妻少妇| 久久久精品欧美日韩精品| 99久国产av精品国产电影| 极品教师在线视频| 中文字幕av在线有码专区| 老师上课跳d突然被开到最大视频| 少妇熟女欧美另类| 在线免费十八禁| 一本久久中文字幕| 级片在线观看| www.av在线官网国产| 国产精品福利在线免费观看| 国产黄色小视频在线观看| 最近2019中文字幕mv第一页| 99热这里只有是精品50| 99久久久亚洲精品蜜臀av| 国产国拍精品亚洲av在线观看| 1024手机看黄色片| 国产成人午夜福利电影在线观看| 国产午夜精品久久久久久一区二区三区| 欧美日韩乱码在线| 国产亚洲av片在线观看秒播厂 | 在现免费观看毛片| 欧美一级a爱片免费观看看| 寂寞人妻少妇视频99o| 精品一区二区免费观看| 日韩高清综合在线| 成人综合一区亚洲| 国产精品久久久久久久电影| av视频在线观看入口| 国产精品伦人一区二区| 国产亚洲精品久久久久久毛片| 亚洲av男天堂| 在线观看午夜福利视频| 精品久久久久久久久久免费视频| 欧美日韩在线观看h| а√天堂www在线а√下载| 国产单亲对白刺激| 欧美日韩综合久久久久久| av又黄又爽大尺度在线免费看 | 特大巨黑吊av在线直播| 国产精品久久久久久久电影| 18禁在线播放成人免费| 观看免费一级毛片| 蜜臀久久99精品久久宅男| 日本免费a在线| 亚洲av电影不卡..在线观看| 国产精品一区二区三区四区免费观看| 日韩制服骚丝袜av| 日韩 亚洲 欧美在线| 在线国产一区二区在线| 亚洲欧洲国产日韩| 国产精品久久久久久精品电影小说 | 麻豆一二三区av精品| 国产三级在线视频| 最后的刺客免费高清国语| 欧美潮喷喷水| 人人妻人人澡欧美一区二区| 在线观看午夜福利视频| 亚洲无线观看免费| 99久久中文字幕三级久久日本| 人妻少妇偷人精品九色| 亚洲经典国产精华液单| 国产一区二区亚洲精品在线观看| 午夜精品国产一区二区电影 | 国模一区二区三区四区视频| 少妇的逼水好多| 91av网一区二区| 午夜老司机福利剧场| 日韩视频在线欧美| 深夜a级毛片| 国产成人影院久久av| 亚洲丝袜综合中文字幕| 男人舔女人下体高潮全视频| 乱系列少妇在线播放| 欧美日韩国产亚洲二区| 美女黄网站色视频| 欧美丝袜亚洲另类| 夜夜看夜夜爽夜夜摸| 国产一区二区在线av高清观看| 99热这里只有是精品在线观看| 网址你懂的国产日韩在线| 国产一区亚洲一区在线观看| 国产精品伦人一区二区| 欧美一区二区精品小视频在线| 听说在线观看完整版免费高清| 99热这里只有精品一区| 欧美+亚洲+日韩+国产| 亚洲aⅴ乱码一区二区在线播放| 老女人水多毛片| 国产三级中文精品| 国产精品久久视频播放| 国产亚洲精品av在线| 国内揄拍国产精品人妻在线| 欧美xxxx性猛交bbbb| 国产精品不卡视频一区二区| www.色视频.com| 最近2019中文字幕mv第一页| 夫妻性生交免费视频一级片| 青春草亚洲视频在线观看| 内地一区二区视频在线| 亚洲精品日韩在线中文字幕 | 亚洲人成网站在线播放欧美日韩| 国产视频首页在线观看| 女的被弄到高潮叫床怎么办| 麻豆久久精品国产亚洲av| 成人特级av手机在线观看| 女的被弄到高潮叫床怎么办| 神马国产精品三级电影在线观看| 一区二区三区免费毛片| 在线播放无遮挡| 久久亚洲精品不卡| 欧美丝袜亚洲另类| 国产一级毛片七仙女欲春2| 激情 狠狠 欧美| 久久草成人影院| 国产精品久久久久久久电影| а√天堂www在线а√下载| 国产日韩欧美在线精品| 亚洲精品久久久久久婷婷小说 | 日本欧美国产在线视频| 少妇的逼水好多| 久久婷婷人人爽人人干人人爱| 日韩亚洲欧美综合| 精品少妇黑人巨大在线播放 | 插阴视频在线观看视频| 国产伦一二天堂av在线观看| 国产av不卡久久| 青春草国产在线视频 | 18禁在线播放成人免费| 午夜久久久久精精品| 亚洲真实伦在线观看| av卡一久久| 日韩大尺度精品在线看网址| 男女下面进入的视频免费午夜| 国产老妇伦熟女老妇高清| 久久鲁丝午夜福利片| 美女黄网站色视频| 99久国产av精品| 欧美日韩综合久久久久久| 国产视频首页在线观看| 亚洲在线自拍视频| 18禁黄网站禁片免费观看直播| 97热精品久久久久久| 一级毛片aaaaaa免费看小| 日本与韩国留学比较| 男女视频在线观看网站免费| 欧美日本视频| 我的女老师完整版在线观看| 老司机福利观看| 亚洲电影在线观看av| 我要看日韩黄色一级片| 日韩视频在线欧美| 国产精品,欧美在线| 99热全是精品| 变态另类成人亚洲欧美熟女| 日本爱情动作片www.在线观看| 欧美最黄视频在线播放免费| 久久99热这里只有精品18| 少妇人妻精品综合一区二区 | 欧美性感艳星| 国产大屁股一区二区在线视频| 中文欧美无线码| 99久国产av精品国产电影| 不卡视频在线观看欧美| 久久久久久久午夜电影| 黄色欧美视频在线观看| 亚洲av男天堂| 如何舔出高潮| 12—13女人毛片做爰片一| 亚洲aⅴ乱码一区二区在线播放| 看黄色毛片网站| 亚洲欧美成人综合另类久久久 | 色尼玛亚洲综合影院| 波多野结衣高清无吗| 国产成人影院久久av| 欧美人与善性xxx| 少妇的逼好多水| 高清毛片免费看| 日韩av不卡免费在线播放| 高清毛片免费观看视频网站| 国产白丝娇喘喷水9色精品| av免费在线看不卡| 成熟少妇高潮喷水视频| 日产精品乱码卡一卡2卡三| 美女cb高潮喷水在线观看| 欧美人与善性xxx| 欧美精品国产亚洲| 国产亚洲av片在线观看秒播厂 | 国产白丝娇喘喷水9色精品| 亚洲av成人av| 一卡2卡三卡四卡精品乱码亚洲| 久久热精品热| 最近手机中文字幕大全| 国产精品一区二区在线观看99 | 日本色播在线视频| 99在线视频只有这里精品首页| 亚洲国产精品合色在线| 亚洲最大成人av| 久久久a久久爽久久v久久| 国产三级在线视频| 午夜激情欧美在线| 12—13女人毛片做爰片一| 亚洲五月天丁香| 欧美日韩国产亚洲二区| 99九九线精品视频在线观看视频| 国产在线男女| 长腿黑丝高跟| 亚洲不卡免费看| 久久国产乱子免费精品| 精品无人区乱码1区二区| av.在线天堂| 亚洲高清免费不卡视频| 99久久九九国产精品国产免费| 国产一区二区三区在线臀色熟女| 中文字幕av在线有码专区| 免费不卡的大黄色大毛片视频在线观看 | 中国国产av一级| 真实男女啪啪啪动态图| 久久久久久伊人网av| 男人狂女人下面高潮的视频| 久久久午夜欧美精品| 欧美色视频一区免费| 美女脱内裤让男人舔精品视频 | 欧美日本视频| 亚洲精品粉嫩美女一区| 亚洲三级黄色毛片| 黑人高潮一二区| 亚洲国产精品久久男人天堂| 26uuu在线亚洲综合色| 99热6这里只有精品| 日本黄大片高清| videossex国产| 欧美不卡视频在线免费观看| 天美传媒精品一区二区| 91狼人影院| 99在线人妻在线中文字幕| 国产91av在线免费观看| 热99在线观看视频| 女的被弄到高潮叫床怎么办| 五月伊人婷婷丁香| 蜜桃亚洲精品一区二区三区| 深夜a级毛片| 亚洲欧美日韩卡通动漫| 久久人妻av系列| 日本与韩国留学比较| 午夜激情福利司机影院| 在线观看av片永久免费下载| 午夜激情欧美在线| 国产爱豆传媒在线观看| 人妻少妇偷人精品九色| 又粗又爽又猛毛片免费看| 国产亚洲av嫩草精品影院| av福利片在线观看| av在线播放精品| 亚洲精华国产精华液的使用体验 | 少妇的逼水好多| 天堂√8在线中文| av黄色大香蕉| 看非洲黑人一级黄片| 国产精品久久久久久久久免| 欧美日韩在线观看h| 最近的中文字幕免费完整| av在线蜜桃| 国产精品嫩草影院av在线观看| .国产精品久久| 一级毛片久久久久久久久女| 91精品一卡2卡3卡4卡| 国产一区二区激情短视频| 啦啦啦韩国在线观看视频| 菩萨蛮人人尽说江南好唐韦庄 | 成人综合一区亚洲| 97超视频在线观看视频| 99久久成人亚洲精品观看| 一区福利在线观看| 亚洲国产精品久久男人天堂| 一区福利在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产精品无大码| 欧美色欧美亚洲另类二区| av在线蜜桃| 午夜免费激情av| 亚洲经典国产精华液单| 欧美xxxx黑人xx丫x性爽| 网址你懂的国产日韩在线| 搡老妇女老女人老熟妇| 内射极品少妇av片p| 舔av片在线| 中国美白少妇内射xxxbb| 男女边吃奶边做爰视频| 中文字幕av在线有码专区| 又粗又硬又长又爽又黄的视频 | 国产高清激情床上av| 干丝袜人妻中文字幕| 欧美又色又爽又黄视频| 免费在线观看成人毛片| 久久久久久九九精品二区国产| 国产黄色小视频在线观看| 欧美色视频一区免费| 日韩高清综合在线| 真实男女啪啪啪动态图| 国产爱豆传媒在线观看| 有码 亚洲区| 精品午夜福利在线看| 老司机影院成人| 国产伦在线观看视频一区| 变态另类丝袜制服| 国产精品永久免费网站| 寂寞人妻少妇视频99o| 最近视频中文字幕2019在线8| 欧美成人精品欧美一级黄| 亚洲电影在线观看av| 国产成人a∨麻豆精品| 白带黄色成豆腐渣| 丝袜美腿在线中文| 三级经典国产精品| 国产精品久久视频播放| 国产一级毛片七仙女欲春2| 午夜精品在线福利| 国内少妇人妻偷人精品xxx网站| 少妇人妻精品综合一区二区 | 成人特级av手机在线观看| 波多野结衣巨乳人妻| 最近视频中文字幕2019在线8| 国产亚洲精品av在线| 日韩成人av中文字幕在线观看| 亚洲性久久影院| 一级二级三级毛片免费看| 国产伦在线观看视频一区| 九九热线精品视视频播放| 少妇的逼水好多| 舔av片在线| 免费看光身美女| 日日摸夜夜添夜夜爱| 九草在线视频观看| 变态另类丝袜制服| 极品教师在线视频| 91久久精品电影网| 色综合亚洲欧美另类图片| 天天躁日日操中文字幕| 亚洲三级黄色毛片| 国产午夜精品久久久久久一区二区三区| 免费搜索国产男女视频| 可以在线观看毛片的网站| 亚洲熟妇中文字幕五十中出| 欧美变态另类bdsm刘玥| 在线观看美女被高潮喷水网站| 欧美+日韩+精品| 校园人妻丝袜中文字幕| 国产精品一区二区三区四区免费观看| 最近手机中文字幕大全| 日本-黄色视频高清免费观看| av女优亚洲男人天堂| 91午夜精品亚洲一区二区三区| 久久久久性生活片| 日韩精品有码人妻一区| 真实男女啪啪啪动态图| 亚洲精品国产成人久久av| 中文资源天堂在线| 国产黄色小视频在线观看| 一个人看的www免费观看视频| 嫩草影院入口| 高清午夜精品一区二区三区 | 日本黄大片高清| 精品人妻熟女av久视频| 村上凉子中文字幕在线| 成人二区视频| 久久婷婷人人爽人人干人人爱| 国产精品不卡视频一区二区| 激情 狠狠 欧美| 国产精品乱码一区二三区的特点| 18禁在线无遮挡免费观看视频| 国产老妇女一区| 国产精华一区二区三区| 精品人妻偷拍中文字幕| 男人狂女人下面高潮的视频| 久久精品夜夜夜夜夜久久蜜豆| 日本在线视频免费播放| 一边亲一边摸免费视频| 久久久久久久久中文| 成人漫画全彩无遮挡| 亚洲最大成人中文| 一边摸一边抽搐一进一小说| 国内精品美女久久久久久| 国产精品一区www在线观看| 久久99精品国语久久久| 午夜福利成人在线免费观看| 搡女人真爽免费视频火全软件| 丝袜喷水一区| 久久精品影院6| 精品久久久久久久人妻蜜臀av| 一级毛片久久久久久久久女| 亚洲成人av在线免费| 老司机影院成人| 自拍偷自拍亚洲精品老妇| 中出人妻视频一区二区| 日日撸夜夜添| 天堂网av新在线| 美女xxoo啪啪120秒动态图| 老熟妇乱子伦视频在线观看| 亚洲最大成人手机在线| 日本五十路高清| 欧洲精品卡2卡3卡4卡5卡区| 欧美精品国产亚洲| 搞女人的毛片| av卡一久久| 一进一出抽搐gif免费好疼| 桃色一区二区三区在线观看| 色噜噜av男人的天堂激情| 国产日韩欧美在线精品| 黄色视频,在线免费观看| 人人妻人人澡人人爽人人夜夜 | 一边亲一边摸免费视频| 亚洲,欧美,日韩| 国产成人a∨麻豆精品| 亚洲成人久久性| 日本黄色片子视频| 高清日韩中文字幕在线| 赤兔流量卡办理| 亚洲欧美日韩卡通动漫| 亚洲成人久久性| 国产69精品久久久久777片| 国产伦精品一区二区三区视频9| 国产精品麻豆人妻色哟哟久久 | 久久精品影院6| 色5月婷婷丁香| 一级毛片我不卡| 91久久精品电影网| 国产精品久久久久久久久免| 亚洲国产欧洲综合997久久,| 国产伦一二天堂av在线观看| 久久久久九九精品影院| 天天躁夜夜躁狠狠久久av| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久久黄片| 尾随美女入室| 美女cb高潮喷水在线观看| 2022亚洲国产成人精品| 中文字幕人妻熟人妻熟丝袜美| 变态另类丝袜制服| 色综合色国产| 久久午夜亚洲精品久久| 亚洲成人精品中文字幕电影| 国产又黄又爽又无遮挡在线| 一级毛片久久久久久久久女| 国产熟女欧美一区二区| 麻豆乱淫一区二区| 日韩欧美一区二区三区在线观看| 亚洲av中文字字幕乱码综合| 亚洲高清免费不卡视频| 91午夜精品亚洲一区二区三区| 高清毛片免费看| 国产一区二区激情短视频| 白带黄色成豆腐渣| 乱系列少妇在线播放| 一区福利在线观看| 日韩一区二区视频免费看| 女人十人毛片免费观看3o分钟| 国内精品美女久久久久久| 女人被狂操c到高潮| 美女黄网站色视频| 搞女人的毛片| 男人和女人高潮做爰伦理| av免费观看日本| 国产亚洲av片在线观看秒播厂 | 精品人妻视频免费看| 国产色爽女视频免费观看| 日本免费a在线| 久久精品91蜜桃| 啦啦啦啦在线视频资源| 亚洲国产欧美在线一区| 日韩国内少妇激情av| 亚洲av免费高清在线观看| 91aial.com中文字幕在线观看| 真实男女啪啪啪动态图| 青青草视频在线视频观看| 成人特级黄色片久久久久久久| 99热这里只有是精品在线观看| 久久精品国产清高在天天线| 国产伦在线观看视频一区| 春色校园在线视频观看| 只有这里有精品99| 观看免费一级毛片| 青青草视频在线视频观看| 精品一区二区三区人妻视频| 特大巨黑吊av在线直播| 九色成人免费人妻av| 国产精品国产高清国产av| 亚洲欧美中文字幕日韩二区| 26uuu在线亚洲综合色| 中文字幕久久专区| 卡戴珊不雅视频在线播放| 一区二区三区高清视频在线| 久久草成人影院| 亚洲欧洲日产国产| www.av在线官网国产| 久99久视频精品免费| 99久久精品热视频| 成人欧美大片| 亚洲五月天丁香| 久久99蜜桃精品久久| 国产大屁股一区二区在线视频| 免费看美女性在线毛片视频| 99视频精品全部免费 在线| 麻豆乱淫一区二区| 美女黄网站色视频| 亚洲国产精品合色在线| 麻豆国产av国片精品| 婷婷六月久久综合丁香| 国产在线男女| 亚洲五月天丁香| 国产午夜精品久久久久久一区二区三区| 99热这里只有是精品50| 国产高清有码在线观看视频| 22中文网久久字幕| 九九爱精品视频在线观看| 亚洲精品国产成人久久av| 中文亚洲av片在线观看爽| 99热全是精品| 久久久久久久久久黄片| 婷婷亚洲欧美| 国产美女午夜福利| www.色视频.com| 亚洲,欧美,日韩| av女优亚洲男人天堂| 久久这里只有精品中国| 久久久久久久久中文| 亚洲最大成人手机在线| 精华霜和精华液先用哪个| 国产激情偷乱视频一区二区|