Issue |
A&A
Volume 616, August 2018
|
|
---|---|---|
Article Number | A161 | |
Number of page(s) | 8 | |
Section | Stellar atmospheres | |
DOI | https://doi.org/10.1051/0004-6361/201732222 | |
Published online | 06 September 2018 |
Long-term variations in the X-ray activity of HR 1099
1
Hamburger Sternwarte, Universität Hamburg, Gojenbergsweg 112, 21029, Hamburg, Germany
e-mail: vperdelwitz@hs.uni-hamburg.de
2
Departamento de Astronomía, Facultad Ciencias Físicas y Matemáticas, Universidad de Concepción, Av. Esteban Iturra s/n Barrio Universitario, Casilla 160-C, Concepción, Chile
Received:
1
November
2017
Accepted:
5
June
2018
Context. Although timing variations in close binary systems have been studied for a long time, their underlying causes are still unclear. A possible explanation is the so-called Applegate mechanism, where a strong, variable magnetic field can periodically change the gravitational quadrupole moment of a stellar component, thus causing observable period changes. One of the systems exhibiting such strong orbital variations is the RS CVn binary HR 1099, whose activity cycle has been studied by various authors via photospheric and chromospheric activity indicators, resulting in contradicting periods.
Aims. We aim at independently determining the magnetic activity cycle of HR 1099 using archival X-ray data to allow for a comparison to orbital period variations.
Methods. Archival X-ray data from 80 different observations of HR 1099 acquired with 12 different X-ray facilities and covering almost four decades were used to determine X-ray fluxes in the energy range of 2−10 keV via spectral fitting and flux conversion. Via the Lomb-Scargle periodogram we analyze the resulting long-term X-ray light curve to search for periodicities.
Results. We do not detect any statistically significant periodicities within the X-ray data. An analysis of optical data of HR 1099 shows that the derivation of such periods is strongly dependent on the time coverage of available data, since the observed optical variations strongly deviate from a pure sine wave. We argue that this offers an explanation as to why other authors derive such a wide range of activity cycle periods based on optical data. We furthermore show that X-ray and optical variations are correlated in the sense that the star tends to be optically fainter when it is X-ray bright.
Conclusions. We conclude that our analysis constitutes, to our knowledge, the longest stellar X-ray activity light curve acquired to date, yet the still rather sparse sampling of the X-ray data, along with stochastic flaring activity, does not allow for the independent determination of an X-ray activity cycle.
Key words: stars: magnetic field / binaries: spectroscopic / stars: activity / stars: coronae / X-rays: binaries
© ESO 2018
1. Introduction
Stellar activity cycles, first discovered for our Sun by Schwabe (1844), have been observed in various stars via photospheric and chromospheric activity indicators, both in single and binary systems (e.g., Berdyugina 2005; Hall et al. 2007). With the advent of space-borne X-ray observatories it has become possible to extend this research to the X-ray range, which has been shown to be a robust tracer of activity cycles (Favata et al. 2004; Robrade et al. 2012; Lalitha & Schmitt 2013; Wargelin et al. 2017).
In addition to enabling the study of the underlying dynamo, magnetic activity cycles have been suggested to be the cause of timing variations in close binary systems. These have been observed over a broad range of such systems, including Algol and RS CVn systems, with relative period modulation amplitudes of ΔP/P ~ 10−5 and typical cycle periods of 30−50 yrs, as well as cataclysmic variables and W UMa systems with ΔP/P ~ 10−6 and cycles of 5−30 yrs (e.g., Warner 1988; Hall 1989). Post-common-envelope binaries (PCEBs) exhibit a very similar be-havior, with ΔP/P ~ 10−6 and modulation periods of 10−30 yrs (Zorotovic & Schreiber 2013; Bours et al. 2016).
A simple explanation of these period variations is to invoke the presence of a third body, giving rise to the variation via the light travel effect (LTE) effect (see e.g., Pribulla et al. 2012). However, while the existence of such circumbinary companions has been postulated to explain timing variations in various systems, to date none of them has been confirmed with an independent method; for example, the third body causing timing variations in V471 Tau first proposed by Beavers et al. (1986) could not be confirmed with a direct imaging attempt with VLT-SPHERE (Hardy et al. 2015).
It is therefore essential to assess whether other mechanisms can produce the observed period variations. Such a mechanism could be a powerful magnetic activity cycle as expected in the presence of rapidly rotating and convective stars (Baliunas et al. 1996). During the course of a cycle, the stellar angular momentum could be redistributed, thus leading to quasi-periodic changes of the stellar quadrupole moment and orbital period. This model, originally proposed by Applegate (1992), has been improved by Lanza & Rodonò (1999) by adopting a consistent formulation for virial equilibrium. Brinkworth et al. (2006) extended the model by introducing a finite-shell formalism, showing that the latter substantially increases the energy required to drive the Applegate mechanism. Völschow et al. (2016) applied this model to a total of 16 close binary systems and showed that the mechanism is clearly feasible in four of the systems, while it can be ruled out on energetic grounds in eight cases; an additional amount of four systems remained unclear. Further, improved models for the Applegate mechanism have been put forward by Lanza (2005, 2006). The authors found that if the observed period variations were driven by the Applegate mechanism, the energy dissipation in the convection zone would exceed the star’s luminosity. As shown recently, the resulting changes of the quadrupole moment may also affect the mass transfer rate between the stars, and potentially give rise to the long cycles in Double Periodic Variables (Schleicher & Mennickent 2017).
One of the systems exhibiting substantial period variations is HR 1099 (or V711 Tau), a RS CVn binary of spectral type K1 IV+G5 V with an orbital period of 2.84 days first identified by Bopp & Fekel (1976). Its proximity of 30.7 pc (van Leeuwen 2007) and its very high level of magnetic activity in all available indicators have made it one of the prime candidates for the study of stellar dynamos (see e.g., Donati 1999, García-Alvarez et al. 2002, Osten et al. 2004 and Cao & Gu 2015). Several authors have searched for magnetic activity periods in HR 1099 via long-term photometry, which is sensitive to detecting variations in the spot coverage. Such activity cycles have been reported by Lanza et al. (2006), who obtain a 19.5 ± 2 yr period, by Berdyugina & Henry (2007), who derive a period of 15−16 yr, and by Muneer et al. (2010) who report a 14.1 ± 0.3 yr period. An inspection of the light curves shows clear periodicities despite a very large scatter in the optical photometry. In addition to their activity period derivation, Muneer et al. (2010) show that the O−C diagram of HR 1099 can be improved considerably by fitting the orbital period of the system to the O–C data, resulting in a well-defined sinusoidal shape with a period of 36.3 ± 1.9 yrs.
As a consequence, the question arises whether magnetic activity is also periodic in other proxy indicators. Drake et al. (2014) have carried out study into the long-term X-ray variability of AR Lac found that it is constant and displays no signs of a coronal activity cycle. In this paper we therefore independently investigate the magnetic activity period of HR 1099 using an X-ray light curve spanning almost 40 yrs in order to test whether this method is a valid means to test the Applegate mechanism. A theoretical motivation is given in Sect. 2 and the observational data are described in Sect. 3. A summary and discussion of our results is given in Sect. 4.
2. Theoretical motivation
We show in the following that a long-term magnetic activity cycle can be qualitatively understood through simple dynamo models. For this purpose, we employ the relation between rotation period and activity cycle given as (Soon et al. 1993; Baliunas et al. 1996)
where D is the dynamo number and α is a parameter of order 0.5. Furthermore, D is related to the Rossby number, Ro, via
with Ro ∝ P/τ, where P is the period and τ is the convective turnover time, can also be obtained from (Soker 2000)
with vc denoting the convective velocity, Hp the pressure scale height, and ω, ωKep, and PKep the angular velocity, Keplerian angular velocity, and Keplerian orbital period, respectively. The convective velocity and the pressure scale height are properties inside the star, which we calculate using the stellar evolution code MESA (Paxton et al. 2011). We first evolve a Sun-like star with solar-like rotation up to 4.7 Gyr, which we use to calibrate α in Eq. (1). The model employs a solar metallicity, and we choose mixing length theory MLT = ‘ML1’ with αMLT = 1.5. Considering that the solar dynamo is expected to take place at the bottom of the convection zone with R ~ 0.7 R⊙, we use the properties at that point and find that α ~ 0.86. We subsequently run a simulation for the subgiant in HR 1099 assuming that it is tidally locked, so that its rotation corresponds to a fixed fraction of the critical angular velocity with ΩZAMS/Ωcrit = 0.28 following from the binary rotation period1. The mass of the subgiant is set to M = 1.3 M⊙ based on Bopp & Fekel (1976). The run is stopped when the radius is R = 4.0 R⊙ . We calculate the dynamo number D as a function of radius from the convection velocity vc and the pressure scale height Hp, allowing us to compute the expected dynamo activity cycle as a function of radius. In Fig. 1, the black curve shows the fraction Pobs/Pcalc as a function of the radial coordinate r/R and the horizontal red line is Pobs/Pcalc = 1, where Pobs is assumed here to be of the order of 20 yrs. They intersect at r/R ~ 0.8.
While based on simplifying assumptions, the result suggests that a dynamo cycle of ~20 yrs would be predominantly driven in the outer parts of the shell. We further make a preliminary assessment if such a dynamo cycle may be responsible for the O–C variations observed in the system through the Applegate mechanism. As a remark, we note that the ratio between activity cycle and O–C variation cycle is approximately 1:2, as is the case in previously observed systems (e.g., Lanza et al. 2001). The expected ratio is however not well explored in theoretical models. To assess whether the Applegate mechanism is energetically feasible, we employ the formalism developed by Völschow et al. (2016) using the two-zone model, which has been implemented in a publicly available Applegate calculator2. In this calculation, we neglect the mass contribution from the radiative core, where dynamo effects and related changes of the angular momentum cannot be expected to be relevant (see also discussion in Lanza 2006). We first recalculate the k1 and k2 coefficients defined by Völschow et al. (2016) based on the MESA output neglecting the central core. The lowest energy is found when we consider the exchange of angular momentum between the shells with 0.2 < r/Rin < 0.65 and 0.65 < r/Rout < 1. The result suggests that the Applegate mechanism may be energetically feasible, and would be driven somewhat more in the interior of the star.
![]() |
Fig. 1. Radial profile of Pobs/Pcalc obtained from Eq. (1) (black curve) and Pobs/Pcalc = 1 (horizontal red line). The lines meet at r/R ~ 0.8, suggesting that the observed activity cycle is driven in this region. |
3. Observations
HR 1099 is one of the brightest known coronal X-ray sources. Therefore HR 1099 has been observed – starting in 1978 – by virtually all X-ray facilities, even those with small effective areas and lower sensitivity, thus forty years of X-ray observations are available. Specifically, we have identified a total of 78 data sets, 60 in the form of spectra and 18 as count rates, from 12 different facilities. Many of these observations have been carried out for calibration purposes or with the goal of studying stellar flares (Kashyap & Drake 1999; Brinkman et al. 2001; Osten et al. 2004; Tsuboi et al. 2016), yet, to our knowledge, no long-term activity study of HR 1099 based on several data sets has been carried out so far. In Table 1 we provide a detailed list of these X-ray observations with the appropriate references.
For the derivation we have to distinguish between two types of data, namely (i) spectra and (ii) data without spectral information, which we converted into fluxes assuming a model based on XMM-Newton spectra.
Instruments in use.
3.1. Spectral data
The majority of the spectral X-ray data of HR 1099 are available as spectral pipeline products from the NASA HEASARC3 data base; we used the HEASARC data whenever available. Specifically, we custom processed the XMM-Newton data (acquired from the ESA archive) with the Science Analysis System (Gabriel et al. 2004) and downloaded Chandra spectra from the TGCat database (Huenemoerder et al. 2011). MAXI data were obtained via the on demand tool on the project website4, with a binning of 2 yrs due to the small effective area of the instrument and target and background apertures of 2 and 3 deg. Swift data were generated with the Build Swift-XRT products tool (Evans et al. 2009) on the Swift homepage5.
For the following spectral analysis we used XSPEC version 12.7.1 (Arnaud 1996). The 60 spectra were grouped with the grppha task to a minimum of 15 cts per spectral bin and clipped to the range of 2–10 keV, a spectral range covered by all but one instrument, the Einstein SSS. Due to an erroneous normalization of the pipeline products of the Einstein SSS and MPC, it was not possible to fit spectra from both instruments simultaneously and thereby cover the entire wavelength range. This and the presence of an instrumental artifact at ≈3.5 keV led us to utilize only the SSS data and extrapolate the spectra beyond 3 keV by using the XSPEC energies and adding 100 additional linear bins up to 10 keV.
The first spectral fit was performed with a cflux✻(mekal + mekal) model, where the model abundances were linked to each other, leaving six free parameters: flux, kT1, norm1, kT2, norm2 and abundance. All spectra with a reduced χ2 > 2 were then refitted after adding a third mekal component to the model. We were merely interested in obtaining a good flux estimate in the energy range in question, and not in any coronal parameters such as plasma temperatures or abundances. For the majority of the spectra our fits yielded a reduced χ2 in the range of 0.63–1.8 (see Table 1), the exceptions being one Einstein SSS and two RXTE spectra.
The output fluxes from the XSPEC cflux model were used for further analysis along with flux uncertainties obtained with the XSPEC error function for a 90% confidence level.
In the case of multi-instrument missions with more than one instrument covering 2–10 keV, all spectra were fitted simultaneously, with the exception of the photon-starved BeppoSAX data, where a simultaneous fit yielded a reduced χ2 > 3.
3.2. Flux conversion
We converted the 18 data sets from those instruments, for which only count rates were available (ROSAT, Einstein IPC, Ginga) into physical fluxes by computing conversion factors in the following manner. Two XMM-Newton spectra were fitted with a 3-T apec model in the spectral range 0.2–12.0 keV, one with a medium flux level as determined in the analysis described in the previous Section (Obs.-ID 0116890901), and one with the lowest flux of the entire XMM-Newton data set (Obs.-ID 0116200701). After the first fit, the plasma temperatures and abundances were frozen to the closest value available in the HEASARC WebPIMMS tool, after which the fit is repeated. The determined plasma temperatures and relative emission measures, along with the model component normalizations, were used as input for WebPIMMS, and a conversion factor was determined for all instruments in question. The conversion factors derived via the XMM-Newton Obs.-ID 0116890901 were then used to convert instrument count rates into physical fluxes, while the difference between the mid- and low-state conversion factors was used as a systematic error. The inclusion of this systematic error was necessary due to the fact that a flux conversion without any a priori knowledge of the coronal properties during a given observation must consider the entire range of system parameters representative of the system. All conversion factors are displayed in Table 1 along with the other data.
3.3. Light curve analysis
3.3.1. Period analysis
In this Section we analyze possible periodicities in the X-ray data alone. The resulting X-ray light curve comprising of 80 observations is displayed in the top panel of Fig. 2.
We excluded the Swift exposure ID 582894000 from further analysis on the grounds that the derived flux is two orders of magnitude larger compared to the rest of the data; this data set is the only Swift observation of the target connected to a TDRSS message, and therefore was likely triggered by an extreme flare of the system.
The remaining 79 data points were analyzed with the generalized Lomb-Scargle algorithm included in the PyAstronomy package with a normalization following Zechmeister & Kürster (2009) while weighting each observation with the determined flux error; the upper left panel of Fig. 3 shows the resulting periodogram. In order to assess the significance derived via the Lomb-Scargle periodogram, we ran a bootstrap approach with 104 repetitions where, in every repetition, the fluxes were randomly drawn with replacement and distributed to the observation epochs. The FAP is then the ratio between all repetitions with a higher LS power than that of the unperturbed analysis and the total number of repetitions. In the top left panel Fig. 3 we marked the FAP levels of 50 and 80% as dashed blue lines, and in the top right panel the maximum period distribution of all repetitions is displayed as a histogram. While the periodogram exhibits a peak at P = 21.2 yrs, the false alarm probability is higher than 50%, indicating that no significant period could be detected. A folded light curve assuming the putative period of 21.2 yrs (Fig. 4) confirms this finding, that is, at any given phase the entire range of fluxes is covered.
In order to check the influence of flaring events on the period analysis, we checked all available light curves for flares, and identified 33 observations without obvious large flaring events (see last column of Table 1). However, a period analysis of this reduced data set yielded no significant period, most likely due to insufficient temporal coverage.
![]() |
Fig. 2. Light curves of HR 1099 in the X-ray and optical regime. Top panel: X-ray flux of HR 1099 in the energy range 2–10 keV as a function of time. For clarity, we omitted the error bars in this plot, but the values can be found in Table 1. Bottom panel: change in optical flux (V-band) as a function of time. Blue data are taken from Berdyugina & Henry (2007), green data from Muneer et al. (2010) and red data from Jetsu et al. (2017). The black dots represent an annual binning of the data. |
![]() |
Fig. 3. Lomb-Scargle periodograms and distribution of periods in the bootstrap approach of the X-ray and optical data. Upper left panel: Lomb-Scargle periodogram of the X-ray data. The FAP levels of 50 and 80% resulting from the bootstrap analysis described in Sect. 3.3 are marked as blue lines. Upper right panel: distribution of maximum periods determined for all 104 repetitions. Lower panels: results of a similar analysis of the optical data. No FAP levels are displayed here, since none of the bootstrap repetitions yielded a higher Lomb-Scargle power than the unperturbed data. |
![]() |
Fig. 4. Folded X-ray light curve under the assumption of the 21.2 yr period determined in the LS analysis. |
3.3.2. Comparison to optical data
We performed a Lomb-Scargle analysis similar to the one described in Sect. 3.3 on an optical data set comprising of archival data provided by Berdyugina & Henry (2007), Muneer et al. (2010) and Jetsu et al. (2017). The resulting periodogram and distribution of periods in the bootstrap approach is displayed in the lower panels of Fig. 3. Since none of the 104 repetitions yielded a higher LS power than the maximum-likelihood period of 14.8 yrs, no FAP could be determined, but rather an upper limit of 0.01%.
We tested the correlation between optical and X-ray data by binning each data set into bins of one year and then performed a cross-correlation between the two. The data yielded only a marginal correlation (see Fig. 5) with a Pearson coefficient of 0.29. Computing a periodogram for different sets of archival optical data, starting with a window of ten years from the first data set and sequentially adding a year of data, we can show the effects of temporal coverage on derived periods. Fig.6 shows the complete photometric light curve (top panel), as well as the period of maximum Lomb-Scargle power as a function of the upper limit of the time window. While the analysis of a sinusoidal data set with this method should result in near-constant values, we observe a large spread in the derived periods, which clearly depend on the time coverage of available data (as also evidenced by the previously mentioned disagreement between other authors).
![]() |
Fig. 5. Correlation between the optical and X-ray light curve. The best fit linear curve is marked with a red line. |
![]() |
Fig. 6. Dependence of the maximum-likelihood period on temporal coverage of available data. Top: photometric light curve in V-band. Bottom: maximum periods derived via Lomb-Scargle analysis for data sets with increasing temporal coverage. The x-value of each point is the upper limit of the analyzed time window. The range of derived periods has a strong dependence on the time window. |
4. Discussion and conclusions
Our study of the longterm X-ray evolution of the active RS CVn binary HR 1099 has produced the first X-ray light curve of the system spanning several decades. We stress that most data used for this work are available as standard pipeline products.
We must nevertheless add some caveats: The X-ray instrumentation used substantially differs, some of the instruments are imaging, others are non-imaging, some instruments offer high spectral resolution, others none and the covered spectral bands differ. Some instruments cover only the energy band above 2 keV (such as RXTE or MAXI); for those instruments HR 1099 appears as a “soft” stellar source and the derived fluxes depend sensitively on the temperatures used. Other instruments, such as XMM-Newton cover the whole band pass, yet others, like the ROSAT PSPC have little to no sensitivity at energies above 2 keV and therefore the flux needs to be extrapolated, which again depends on the spectral model used. We argue that our choices for the flux conversion are reasonable, however, we also admit that they are likely affected by systematic uncertainties that are beyond our control such as the strength of the soft component during the RXTE observations and the strength of the hard component during the ROSAT observations.
Only a marginal correlation between the optical and X-ray data could be determined and we could not detect any significant periodicity in the X-ray data. Furthermore, the bootstrap method applied to estimate the FAP only accounts for white noise, and an approach including red noise would further decrease the significance of the detected peak (Vaughan 2005; VanderPlas 2018). One possible explanation for the large FAP is the substantial scatter within the light curve, which is most likely dominated by stochastic X-ray flaring activity.
Another likely contributor is the shape of the activity modulation itself. Even when taking into account that the peak in the Lomb-Scargle periodograms of both X-ray and optical data is broad, indicating an error ≥5 yrs, the maximum-likelihood periods are only in partial agreement with the photometric cycles derived by other authors, namely Lanza et al. (2006; 19.5 ± 2 yrs), Berdyugina & Henry (2007; 15 − 16) yrs and Muneer et al. (2010; 14.1 ± 0.3 yrs), as well as our own analysis of the complete data set (14.8 yrs). We thus conclude that the activity cycle of HR 1099 is more complicated than a simple sinusoidal, and we could only determine a characteristic timescale of 14.8 yrs based on the optical data.
Lanza et al. (2006) claim a correlation of the magnetic cycle of HR 1099 with the orbital period variations, arguing that the orbital period variations display a period twice as long as that of the activity cycle, hinting at the Applegate mechanism. Muneer et al. (2010) do, in fact, determine the period variations to be 36.3 ± 1.9 yr. However, in order to test this hypothesis, a continued long-term monitoring of HR 1099 and related objects using multiple activity indicators is required.
Here, we use the the notation of Paxton et al. (2013).
Acknowledgments
VP acknowledges funding through the DLR and DFG. VP, DRGS, FHN and JAZF thank for funding through Fondecyt regular (project code 1161247). DRGS is grateful for funding through the “Concurso Proyectos In-ternacionales de Investigación, Convocatoria 2015” (project code PII20150171) and ALMA-Conicyt (project code 31160001). REM acknowledges funding via VRID-ENLACE 214.016.002-1.0 and VRID-ENLACE 218.016.004-1.0. DRGS, FHN and REM further thank for funding through the BASAL Centro de Astrofísica y Tecnologías Afines (CATA) PFB-06/2007. This research has made extensive use of data and/or software provided by the High Energy Astrophysics Science Archive Research Center (HEASARC), which is a service of the Astrophysics Science Division at NASA/GSFC and the High Energy Astrophysics Division of the Smithsonian Astrophysical Observatory. This research has made use of MAXI data provided by RIKEN, JAXA and the MAXI team. It is partly based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. The scientific results reported in this article are based in part on data obtained from the Chandra Data Archive. This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester. We would like to thank Steve Drake, Michael Corcoran and Bryan Irby (HEASARC) for help with the Einstein data set.
References
- Applegate, J. H. 1992, ApJ, 385, 621 [NASA ADS] [CrossRef] [Google Scholar]
- Arnaud, K. A. 1996, in Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby, & J. Barnes, Astronomical Society of the Pacific Conference Series, 101, 17 [NASA ADS] [Google Scholar]
- Baliunas, S. L., Nesme-Ribes, E., Sokoloff, D., & Soon, W. H. 1996, ApJ, 460, 848 [NASA ADS] [CrossRef] [Google Scholar]
- Beavers, W. I., Lui, A., & Herczeg, T. J. 1986, ApJ, 300, 785 [NASA ADS] [CrossRef] [Google Scholar]
- Berdyugina, S. V. 2005, Liv. Rev. Sol. Phys., 2, 8 [Google Scholar]
- Berdyugina, S. V., & Henry, G. W. 2007, ApJ, 659, L157 [NASA ADS] [CrossRef] [Google Scholar]
- Boella, G., Butler, R. C., Perola, G. C., et al. 1997, A&AS, 122, 199 [Google Scholar]
- Bopp, B. W., & Fekel, Jr., F. 1976, AJ, 81, 771 [NASA ADS] [CrossRef] [Google Scholar]
- Bours, M. C. P., Marsh, T. R., Parsons, S. G., et al. 2016, MNRAS, 460, 3873 [NASA ADS] [CrossRef] [Google Scholar]
- Brinkman, A. C., Behar, E., Güdel, M., et al. 2001, A&A, 365, L324 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Brinkworth, C. S., Marsh, T. R., Dhillon, V. S., & Knigge, C. 2006, MNRAS, 365, 287 [NASA ADS] [CrossRef] [Google Scholar]
- Canizares, C. R., Davis, J. E., Dewey, D., et al. 2005, PASP, 117, 1144 [NASA ADS] [CrossRef] [Google Scholar]
- Cao, D., & Gu, S. 2015, MNRAS, 449, 1380 [NASA ADS] [CrossRef] [Google Scholar]
- Donati, J.-F. 1999, MNRAS, 302, 457 [NASA ADS] [CrossRef] [Google Scholar]
- Drake, J. J., Ratzlaff, P., Kashyap, V., et al. 2014, ApJ, 783, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2009, MNRAS, 397, 1177 [NASA ADS] [CrossRef] [Google Scholar]
- Favata, F., Micela, G., Baliunas, S. L., et al. 2004, A&A, 418, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gabriel, C., Denby, M., Fyfe, D. J., et al. 2004, in Astronomical Data Analysis Software and Systems (ADASS) XIII, ed. F. Ochsenbein, M. G. Allen, & D. Egret, Astronomical Society of the Pacific Conference Series, 314, 759 [Google Scholar]
- García-Alvarez, D., Foing, B. H., Montes, D., et al. 2002, A&A, 397, 285 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gehrels, N., Chincarini, G., Giommi, P., et al. 2004, ApJ, 611, 1005 [NASA ADS] [CrossRef] [Google Scholar]
- Giacconi, R., Branduardi, G., Briel, U., et al. 1979, ApJ, 230, 540 [NASA ADS] [CrossRef] [Google Scholar]
- Hall, D. S. 1989, Space Sci. Rev., 50, 219 [NASA ADS] [CrossRef] [Google Scholar]
- Hall, J. C., Lockwood, G. W., & Skiff, B. A. 2007, AJ, 133, 862 [NASA ADS] [CrossRef] [Google Scholar]
- Hardy, A., Schreiber, M. R., Parsons, S. G., et al. 2015, ApJ, 800, L24 [NASA ADS] [CrossRef] [Google Scholar]
- Huenemoerder, D. P., Mitschang, A., Dewey, D., et al. 2011, AJ, 141, 129 [NASA ADS] [CrossRef] [Google Scholar]
- Jahoda, K., Swank, J. H., Giles, A. B., et al. 1996, in EUV, X-Ray, and Gamma-Ray Instrumentation for Astronomy VII, ed. O. H. Siegmund, & M. A. Gummin, Proc. SPIE, 2808, 59 [NASA ADS] [CrossRef] [Google Scholar]
- Jansen, F., Lumb, D., Altieri, B., et al. 2001, A&A, 365, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jetsu, L., Henry, G. W., & Lehtinen, J. 2017, ApJ, 838, 122 [NASA ADS] [CrossRef] [Google Scholar]
- Kashyap, V., & Drake, J. J. 1999, ApJ, 524, 988 [NASA ADS] [CrossRef] [Google Scholar]
- Lalitha, S., & Schmitt, J. H. M. M. 2013, A&A, 559, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lanza, A. F. 2005, MNRAS, 364, 238 [NASA ADS] [CrossRef] [Google Scholar]
- Lanza, A. F. 2006, MNRAS, 369, 1773 [NASA ADS] [CrossRef] [Google Scholar]
- Lanza, A. F., & Rodonò, M. 1999, A&A, 349, 887 [Google Scholar]
- Lanza, A. F., Rodonò, M., Mazzola, L., & Messina, S. 2001, A&A, 376, 1011 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lanza, A. F., Piluso, N., Rodonò, M., Messina, S., & Cutispoto, G. 2006, A&A, 455, 595 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Matsuoka, M., Kawasaki, K., Ueno, S., et al. 2009, PASJ, 61, 999 [NASA ADS] [CrossRef] [Google Scholar]
- Muneer, S., Jayakumar, K., Rosario, M. J., Raveendran, A. V., & Mekkaden, M. V. 2010, A&A, 521, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Osten, R. A., Brown, A., Ayres, T. R., et al. 2004, ApJS, 153, 317 [NASA ADS] [CrossRef] [Google Scholar]
- Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [Google Scholar]
- Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4 [NASA ADS] [CrossRef] [Google Scholar]
- Pribulla, T., Vaňko, M., Ammler-von Eiff, M., et al. 2012, Astron. Nachr., 333, 754 [NASA ADS] [CrossRef] [Google Scholar]
- Robrade, J., Schmitt, J. H. M. M., & Favata, F. 2012, A&A, 543, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rothschild, R., Boldt, E., Holt, S., et al. 1979, Space Sci. Instrum., 4, 269 [NASA ADS] [Google Scholar]
- Schleicher, D. R. G., & Mennickent, R. E. 2017, A&A, 602, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schwabe, M. 1844, Astron. Nachr., 21, 233 [NASA ADS] [CrossRef] [Google Scholar]
- Soker, N. 2000, ApJ, 540, 436 [NASA ADS] [CrossRef] [Google Scholar]
- Soon, W. H., Baliunas, S. L., & Zhang, Q. 1993, ApJ, 414, L33 [NASA ADS] [CrossRef] [Google Scholar]
- Tanaka, Y., Inoue, H., & Holt, S. S. 1994, PASJ, 46, L37 [NASA ADS] [Google Scholar]
- Taylor, B. G., Andresen, R. D., Peacock, A., & Zobl, R. 1981, Space Sci. Rev., 30, 479 [NASA ADS] [CrossRef] [Google Scholar]
- Truemper, J. 1982, Adv. Space Res., 2, 241 [NASA ADS] [CrossRef] [Google Scholar]
- Tsuboi, Y., Yamazaki, K., Sugawara, Y., et al. 2016, PASJ, 68, 90 [NASA ADS] [CrossRef] [Google Scholar]
- Turner, M. J. L., Thomas, H. D., Patchett, B. E., et al. 1989, PASJ, 41, 345 [NASA ADS] [Google Scholar]
- van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- VanderPlas, J. T. 2018, ApJS, 236, 16 [NASA ADS] [CrossRef] [Google Scholar]
- Vaughan, S. 2005, A&A, 431, 391 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Völschow, M., Schleicher, D. R. G., Perdelwitz, V., & Banerjee, R. 2016, A&A, 587, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wargelin, B. J., Saar, S. H., Pojmański, G., Drake, J. J., & Kashyap, V. L. 2017, MNRAS, 464, 3281 [NASA ADS] [CrossRef] [Google Scholar]
- Warner, B. 1988, Nature, 336, 129 [NASA ADS] [CrossRef] [Google Scholar]
- Weisskopf, M. C., Brinkman, B., Canizares, C., et al. 2002, PASP, 114, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zorotovic, M., & Schreiber, M. R. 2013, A&A, 549, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A
Data resulting from the X-ray analysis.
Notes. Column one states the X-ray facility, and, in the case of multi-instrument missions, the instrument. The observing date is given as MJD and year in cols. three and four. The mission-specific Observation ID and duration of the exposure are given in cols. two and five. The resulting X-ray flux and error in the energy range 2–10 keV are given in cols. six and seven. The flux conversion factor (if data is derived from count rates) is given in col. eight, the reduced χ2 for spectral fits in col. nine and the flare flag in col. ten indicates whether visual inspection of the light curve indicated a major flare during the observation, where “y” denotes an obvious flaring event, “n” stands for no large flare and “–” marks all observations for which a visual inspection was not possible.
All Tables
All Figures
![]() |
Fig. 1. Radial profile of Pobs/Pcalc obtained from Eq. (1) (black curve) and Pobs/Pcalc = 1 (horizontal red line). The lines meet at r/R ~ 0.8, suggesting that the observed activity cycle is driven in this region. |
In the text |
![]() |
Fig. 2. Light curves of HR 1099 in the X-ray and optical regime. Top panel: X-ray flux of HR 1099 in the energy range 2–10 keV as a function of time. For clarity, we omitted the error bars in this plot, but the values can be found in Table 1. Bottom panel: change in optical flux (V-band) as a function of time. Blue data are taken from Berdyugina & Henry (2007), green data from Muneer et al. (2010) and red data from Jetsu et al. (2017). The black dots represent an annual binning of the data. |
In the text |
![]() |
Fig. 3. Lomb-Scargle periodograms and distribution of periods in the bootstrap approach of the X-ray and optical data. Upper left panel: Lomb-Scargle periodogram of the X-ray data. The FAP levels of 50 and 80% resulting from the bootstrap analysis described in Sect. 3.3 are marked as blue lines. Upper right panel: distribution of maximum periods determined for all 104 repetitions. Lower panels: results of a similar analysis of the optical data. No FAP levels are displayed here, since none of the bootstrap repetitions yielded a higher Lomb-Scargle power than the unperturbed data. |
In the text |
![]() |
Fig. 4. Folded X-ray light curve under the assumption of the 21.2 yr period determined in the LS analysis. |
In the text |
![]() |
Fig. 5. Correlation between the optical and X-ray light curve. The best fit linear curve is marked with a red line. |
In the text |
![]() |
Fig. 6. Dependence of the maximum-likelihood period on temporal coverage of available data. Top: photometric light curve in V-band. Bottom: maximum periods derived via Lomb-Scargle analysis for data sets with increasing temporal coverage. The x-value of each point is the upper limit of the analyzed time window. The range of derived periods has a strong dependence on the time window. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.