EDP Sciences
Home All issues Volume 593 (September 2016) A&A, 593 (2016) A123 Full HTML
Free Access
Issue
A&A
Volume 593, September 2016
Article Number A123
Number of page(s) 12
Section Stellar atmospheres
DOI https://doi.org/10.1051/0004-6361/201628545
Published online 04 October 2016

© ESO, 2016

1. Introduction

Well-sampled photometric monitoring of spotted cool stars allows tracing photospheric brightness variations that are caused by rotational modulation. Long-term photometric data moreover allow investigating spot activity cycles and relate inherent changes of the surface spot distribution to photometric period variations, both of which are an indicator of the underlying surface differential rotation. Phenomena such as active longitudes (Henry at al. 1995; Jetsu 1996) and flip-flops of active longitudes (Jetsu et al. 1991; Berdyugina & Tuominen 1998) were revealed by such long-term photometry, but they are far from being fully understood. Holzwarth & Schüssler (2003a,b) demonstrated on theoretical grounds that the tidal effects of a companion star in a binary system may alter the surface distribution of spots by breaking the axial symmetry, which could lead to preferred longitudes. Magnetic dynamo simulations (Fluri & Berdyugina 2004; Moss 2005; Elstner & Korhonen 2005) of the flip-flop phenomenon basically concluded that the axisymmetric and the non-axisymmetric (magnetic) field components must coexist to produce a flip-flop in the first place. These studies highlighted that meridional circulation and differential rotation are a prerequisite for this type of dynamo but do not yet allow predicting the relation between differential rotation and flip-flop period, for example. Recently, Roettenbacher et al. (2015a,b) showed that there might be ambiguity in interpreting light curves of RS CVn binaries in terms of persistent active longitudes and ellipsoidal distortion. In this case, some part of the light curve might be caused by binarity effects and not by cool spots.

For the characterization of differential surface rotation, time-series Doppler imaging became the most powerful method (see, e.g., Donati & Collier Cameron 1997; Strassmeier 2009; Künstler et al. 2015; Kővári et al. 2016). The differential rotation signature is either based on the longitudinal cross-correlation of images from successive rotation cycles or on the addition of the image shear as an additional free parameter in the line-profile inversion itself. Neither of these variants is free of problems, see Dunstone et al. (2008) for a comparison, but both were successfully applied to a few single and binary stars in different evolutionary stages (e.g., Donati et al. 2003; Dunstone et al. 2008; Kővári et al. 2012). For main-sequence stars, only solar-type differential rotation (equator rotates faster than poles) was found so far, while anti-solar differential rotation (poles rotate faster than equator) was found for a few fast-rotating giants (Weber et al. 2005). Kitchatinov & Rüdiger (2004) have shown that fast meridional circulation in a very deep convection zone with magnetically induced thermal inhomogeneities would lead to anti-solar differential rotation. However, at least in some cases, different authors found contradictory results for one and the same target, for example, for the K sub-giant of the RS CVn binary HR 1099, for which Vogt et al. (1999) analyzed high-resolution unpolarized spectroscopic data and found weak anti-solar differential rotation, while Petit et al. (2004) used high-resolution polarized spectra and found solar-type differential rotation. It is mostly assumed that the spots follow the rotation of the surface, while differential rotation is investigated by photometric spot modeling and/or time-series Doppler imaging. However, Korhonen & Elstner (2011) suggested that large spot structures could reflect geometric properties of the large-scale dynamo and not the surface rotation. Still, the spots are the tracers we rely on.

HD 208472 (V2075 Cyg, HIP 108198, V = 74, Pphot = 22.4 d) is a conspicuous target that seems to show active spot longitudes and an activity cycle with a period of about six years Özdarcan et al. (2010; hereafter Paper I). In 1991, W. Bidelman discovered its strong Ca ii H and K emission, and a number of follow-up studies revealed and confirmed that the system is a single-lined spectroscopic binary with a G8III primary component on a nearly circular orbit that exhibits rotationally modulated light variations with variable light-curve amplitude (Henry at al. 1995; Fekel et al. 1999; Strassmeier et al. 1999; Koen & Eyer 2002). It was also suggested that the system may be a good candidate for Doppler imaging (Henry at al. 1995). First such attempts for Doppler imaging were made by Weber et al. (2001), Weber (2004) and Weber et al. (2005), whose initial study was based on 70 consecutive nights of spectroscopic observations with the NSO McMath-Pierce telescope between 1996 and 1997. They found that spots at that time were concentrated at low and intermediate latitudes without any structures close to the one visible pole of the star. They applied the sheared-image method and found a weak anti-solar differential rotation of α = −0.04 ± 0.02, which, however, was not considered a conclusive result. Erdem et al. (2009) estimated an α parameter of the same amount, but relative to the orbital period, through photometric spot modeling based on BVRcIc photometry in 2006. More comprehensive photometry was presented in our Paper I, where we analyzed 17 yrs of data from Automatic Photoelectric Telescopes (APT). A spot-modeling and photometric-period analysis led us to a differential rotation coefficient of α = 0.004 ± 0.010, which practically indicated a non-detection of differential rotation. Furthermore, we found that the longitudinal position of spots varied coherently with the mean brightness and seemed to exhibit a 6.28 yr period. This cyclic variation was interpreted as a stellar analog of the solar 11-yr sunspot cycle.

In the current study, we present new high-resolution time-series spectroscopy of the system from 2009–2011 together with contemporaneous V-band photometry. We redetermine precise atmospheric parameters, obtain an improved spectroscopic orbit ,and investigate the evolutionary status of the primary star. The unique time-series data from STELLA enables us to obtain Doppler images of the star for many successive rotational cycles, and for three consecutive years. We first summarize the instrumental setup and the data collection and reduction procedures in Sect. 2. In Sect. 3 we derive stellar atmospheric parameters, spectroscopic orbital elements and physical properties of the system. Section 4 comprises the Doppler imaging of the primary star, including data preparation and a brief description of the iMap inversion code. We present a total of 13 surface maps and their analysis in Sect. 5. In the final section, we summarize and discuss our results.

2. Observations and data reductions

Spectroscopic observations were carried out with the robotic STELLA facility on Tenerife, Spain, from March 22, 2009 to July 20, 2011. All spectra were obtained with the fiber-fed STELLA Echelle Spectrograph (SES) with a resolving power of R = 55 000 corresponding to a two-pixel resolution of 0.12 Å at 6400 Å. The spectra cover the wavelength range 3900–8800 Å. In most cases, a single exposure was taken per available night. On some occasions, two successive spectra were obtained. Exposure time was set to 7200 s during the 2009 observing season, and 4500 s thereafter. The average signal-to-noise ratio (S/N) varied between 240 and 160 for the three observing seasons. By 2010, the efficiency for coupling the starlight into the fiber was significantly improved, and equally good S/N spectra were obtained with a shorter exposure time. In total, 229 usable spectra were acquired during three years.

Raw spectra were reduced with the SES data reduction pipeline, which is mainly based on IRAF1 command-language scripts (Weber et al. 2008). CCD images are corrected for bad pixel and cosmic-ray impacts. Average overscan is calculated and subtracted from each image. An already overscan-subtracted master bias frame is applied to the images to remove bias levels. Nightly obtained flat-field images (around dusk, dawn, and midnight) are used to produce a master flat-field image of the related night, and this master image is normalized to unity. Then target spectra are divided by this master flat. Scattered-light correction is done from the inter-order space, and one-dimensional spectra are extracted with the standard IRAF optimal-extraction routine. After this, a wavelength calibration is applied by using consecutively recorded Th-Ar spectra. Finally, continuum normalization of the extracted spectral orders is done by dividing with a flux-normalized synthetic spectrum of spectral classification just like HD 208472.

Photometric observations were conducted with the T7 0.75 m APT Amadeus at Fairborn Observatory in southern Arizona (Strassmeier et al. 1997). Its data were extensively used and described in Paper I, and we refer to this paper for instrumental and observational details.

3. Astrophysical parameters

3.1. Unspotted brightness

Figure 1a shows all available Johnson V-band data of HD 208472. With the addition of our new APT data for the time range from 2010 to 2015, the photometric coverage now extends continuously over 23 yrs. We determined the peak-to-peak light curve amplitudes and minimum brightness levels for each observing season and found that the system was never brighter than by the end of 2014, reaching a V magnitude of 713 ± 0.01. This is even brighter than the predicted unspotted (apparent) brightness of 7323 previously estimated in Paper I. Using the minimum brightness vs. amplitude method of Oláh et al. (1997), we find a new value for the unspotted V brightness of 704 (Fig. 1b). Using the interstellar reddening obtained in Paper I (E(B−V) = 0127), we revise the intrinsic unspotted V brightness of HD 208472 to 665 ± 0.01.

This updated value together with the Hipparcos distance of 153 pc (van Leeuwen 2007) modifies the absolute V magnitude of HD 208472 to +073 ± 0.18. Its luminosity is therefore 64 ± 5 L⊙ (adopting Mbol, ⊙ = +474).

thumbnail Fig. 1

a) Long-term V-band photometry of HD 208472. The dashed line denotes the updated unspotted brightness level. The blue filled dots show data from Paper I, while the red dots denote the recent data. Plus signs show differential magnitudes in the sense of check-minus-comparison star. These measurements are shifted by a constant amount and plotted on a common scale with measurements of the variable star. b) Minimum brightness vs. V amplitude diagram. The straight line is a linear representation of the upper envelope of the scatter.
Open with DEXTER

3.2. Orbital elements

Radial velocities (RV) from STELLA-SES spectra were derived from an order-by-order cross correlation with a synthetic template spectrum and then averaged. A total of 60 orders out of the 80 available were used, and a synthetic template spectrum of a G8 giant was adopted. The standard (external) error of a single STELLA-SES observation of HD 208472 is ≈30 m s-1. No systemic zero-point shift was added to the data in this paper, but see Strassmeier et al. (2012) for its determination with respect to the CORAVEL reference system. The complete data set of 229 RVs is available in electronic form at the CDS Strasbourg.

For an initial orbital solution we adopted the elements from Fekel et al. (1999) as starting values. The preliminary solution showed clear and systematic residuals with a full amplitude of up to 1 km s-1 (Fig. 2a). Two anti-phase waves can be identified (lower panel of the figure), indicating that two spots or spot groups traversed the stellar disk during the time of observation. We applied a discrete Fourier transform to the residuals using the Period04 (Lenz & Breger 2004) software package and searched for significant periods. This was done separately for the individual observing season. In 2010, we subdivided the time coverage into two data groups because the variation pattern had changed between the first and second half of the season. The resulting periods with significance of 5σ or more in amplitude are tabulated in Table 1 together with their amplitudes. The dominant period is close to the photometric (=rotational) period in 2010 and 2011, while the other periods are just the first, second, and third harmonics of it. In 2009, half of the photometric period is the most dominant one and the other periods are the photometric period itself and its other harmonics. This confirms that the RV scatter is due to starspots migrating in and out of sight and not due to an as yet undiscovered planet.

In the next step, we pre-whitened the RV measurements of every individual observing sequence with a least-squares fit from the periods in Table1. These data are called the spot-corrected velocities. We then repeated the orbital solution and obtained the final elements given in Table 2. Figure 2b shows the best orbital fit with respect to the spot-corrected RVs and its residuals. We note that the orbital elements remain almost identical, but the errors were reduced on average by a factor of 4.

thumbnail Fig. 2

a) Radial velocities (dots) before spot correction and the representation of the orbital fit (line). The bottom panel shows the residuals. b) Spot-corrected radial velocities and the final fit for the orbital elements. The rms of the fit decreased by a factor four with respect to the uncorrected velocities.
Open with DEXTER

Table 1

Periodogram results of radial velocity residuals.

Table 2

Spectroscopic orbital elements.

3.3. Astrophysical parameters

Effective temperature, surface gravity, metallicity and microturbulence were determined from each individual SES spectrum with the PARSES pipeline (Allende Prieto 2004; Allende Prieto et al. 2006; Jovanovic et al. 2013). This fits precalculated grids of synthetic spectra to selected wavelength regions of the observed spectra. Grids of synthetic spectra were calculated from MARCS model atmospheres in local thermodynamic equilibrium (LTE) (Gustafsson et al. 2008). The pipeline searches for the best-fit synthetic spectrum by minimizing the difference between the observed and computed spectrum and then adopts its parameters as its best solution. In addition to the above four atmospheric parameters, a line-broadening parameter consisting of the combined effect of vsini and macroturbulence was adopted as a free parameter during the fitting process. According to the atmospheric properties of the star, an assumed macroturbulence of 5 km s-1 was adopted from Gray (2005) and kept fixed.

The PARSES results for 218 SES spectra are plotted as a function of Julian date in Fig. 3. An effective temperature of 4,630 K, log g = 2.25, [Fe / H] = −0.49, and a microturbulence of 1.7 km s-1 compromise the best-fit averages. We excluded 11 of the 229 spectra that grossly deviated from the bulk, possibly because of problems arising from the continuum level. From the remaining individual values, we took unweighted averages and adopted these as our final atmospheric parameters in Table 3. Their rms deviations are estimated to be the respective internal errors.

thumbnail Fig. 3

PARSES spectrum-synthesis results of four atmospheric parameters for each individual SES observation.
Open with DEXTER

The PARSES pipeline returns average vsini values between 18 and 24 km s-1, depending on the assumed radial-tangential macro- and microturbulence and the wavelength range selected. A much more accurate vsini value of 19.4 ± 0.2 km s-1 is possible from our Doppler imaging analysis (see Sect. 4), in good agreement with 19.7 km s-1 given by Fekel et al. (1999). The radius of the star derived with the Stefan-Boltzmann law is then 12.5 ± 0.5 R⊙, a typical value for an early-K giant (solar Teff and Mbol were adopted as 5780 K and 474, respectively). The bolometric correction is adopted from Flower (1996). Combination of this radius and the vsini leads to an inclination angle of 44°. Simple propagation of the internal errors of the measured and calculated parameters results in an unrealistically small uncertainty for the inclination (just a few degrees). In particular for active giant stars, uncertainties on calculated radius and luminosity are expected to be much larger, hence we estimate ±10° for the range of the most likely inclination.

Table 3

Astrophysical parameters of HD 208472. Errors quoted are internal errors.

Figure 4 compares the position of the star in the log Teff vs. log L/L⊙ plane with evolutionary tracks and isochrones from Bertelli et al. (2008). We chose the Y = 0.30 and Z = 0.008 tracks, which are closest to but somewhat less metal deficient than our PARSES value (depending on the adopted metallicity of the Sun). The inset in Fig. 4 provides a detailed zoom into the red-giant branch (RGB). The position of HD 208472 suggests a 1.2 M⊙ mass and an age of ~3.5 Gyr. It also indicates that the star is past the RGB bump. Error ranges of ~0.1 M⊙ in mass and an age range of 2.8–4.2 Gyr are estimated from the error bars in Fig. 4. A summary of the astrophysical parameters is given in Table 3.

Because we do not see evidence of the secondary component at optical wavelengths, we expect at least a brightness difference between the components of ≈25. In this case, the secondary component must be fainter than MV≈ 32. The low-mass function already suggests a low mass for the secondary component. With our best-estimate for the inclination and a 1.2 M⊙ mass for the primary, we estimate ~0.4 M⊙ for the secondary mass. This indicates either a white dwarf or a normal M0 dwarf. Unfortunately, there are no ultraviolet observations for HD 208472, and therefore we cannot decide between these two possibilities.

thumbnail Fig. 4

Position of HD 208472 in the H-R diagram. The tracks (black dashed lines) and isochrones (red continuous lines) are adopted from Bertelli et al. (2008) for Y = 0.30 and Z = 0.008.
Open with DEXTER

4. Doppler imaging

4.1. Technical setup

We used our (Zeeman) Doppler-imaging code iMap (Carroll et al. 2007, 2008, 2009, 2012) to reconstruct the surface temperature distribution from time-series STELLA spectra. The code calculates local line profiles by solving the radiative transfer equation under LTE assumption and inverts very many photospheric absorption lines simultaneously. The iteratively regularized Landweber algorithm (Engl et al. 1996; Carroll et al. 2012) is implemented. During calculations, the surface of the star is divided into 5° × 5° segments; 2592 segments in total. Local line-profile calculations are made for each segment with ATLAS-9 (Castelli & Kurucz 2004) model atmospheres for the temperature range between 3500 K and 6000 K in steps of 250 K and interpolated to the gravity, metallicity, and microturbulence value of HD 208472. For further details about the code iMap, we refer to Carroll et al. (2012).

4.2. Line selection and data preparation

We selected 34 absorption lines that were relatively blend-free for the given temperature and gravity and had a depth larger than 20% of the continuum. These lines were taken into account simultaneously during the inversion. We applied the multiline technique as described in Carroll et al. (2012). Spectral line synthesis was calculated between ±50 km s-1in the velocity domain, however, the velocity range of the line inversion was restricted to ± 35 km s-1, where the spectral line wings reach the continuum. The lines and their standard parameters from VALD-3 (Kupka et al. 2011) are listed in Table 4.

Table 4

Absorption lines used for the mean profile. Atomic data are taken from VALD-3.

We divided the spectra into 13 separate data sets, where each set covers one and a half rotational cycles at most. No set has any spectrum from the previous or the following data set in common. This allowed us to trace the spot evolution and differential rotation without a biasing correlation between subsets that would be due to shared data. We summarize the observing log of all 13 data sets in Table 5.

Table 5

Observing log of Doppler-imaging sets.

The rotational phase for the spectra is based on our new orbital elements by using the orbital period and a modified epoch as zero point such that one-fourth of the orbital period was subtracted from T0 in Table 2 (“spot corrected” column); (1)At this modified epoch, the giant component is in inferior conjunction, that is, zero phase is 0° longitude. We assumed that the star rotates in the counter-clockwise direction and that the longitude increases in the same direction so that phase 0.25 (φ = 0.25) corresponds to 270° longitude on the stellar surface, for instance.

4.3. v sin i validation

Before calculating surface maps, we searched for the vsini value that provided the best representation for the rotationally broadened line profiles. For this purpose, we chose the third data set from the 2009 observing season, which is one of the best sets with very good sampling and high S/N. We reconstructed surface images for this set by iterating vsini values between 17.7 and 21.2 km s-1 while keeping the remaining parameters fixed. A parabolic fit to the distribution of χ2 in Fig. 5 indicates a best-match vsini of 19.4 ± 0.2 km s-1. The error comes from the width of the parabola, ±0.2 km s-1steps in the sampling, where the deviation is 1σ. We adopted this value as our best estimate and kept it fixed for the remaining imaging.

thumbnail Fig. 5

Quality-of-fit (χ2) for a series of trial vsini values. The red dashed curve is a quadratic fit to the data points. The arrow indicates the position of the best-fit vsini of 19.4 km s-1.
Open with DEXTER

5. Results

5.1. Surface maps

Figures 6 to 8 present our Doppler images for 2009 to 2011. We plot observed and inverted line profiles of each set in Fig. A.1. None of the images shows a dominant polar spot, in agreement with earlier maps by Weber (2004). Spots on HD 208472 are generally not very large and on average seem cooler than the undisturbed photosphere by just ~500 K. We did not recover any bright features warmer than the photosphere.

thumbnail Fig. 6

Doppler images of HD 208472 for the 2009 observing season.
Open with DEXTER

thumbnail Fig. 7

Doppler images of HD 208472 for the 2010 observing season.
Open with DEXTER

thumbnail Fig. 8

Doppler images of HD 208472 for the 2011 observing season.
Open with DEXTER

thumbnail Fig. 9

a) Longitudinal distribution of spots on HD 208472 from photometric two-spot models. The shaded regions are the time ranges where spots change their preferential positions (hemispheres) on the surface of the star. The continuous line represents the longitudes of the apsidal line of the binary system, while the two dashed lines restrict the stellar surface. The vertical dot-dashed lines correspond to times when changes occurred in light-curve amplitude and photometric period. We overplot the longitudes of the main spot of the Doppler images with cross symbols in red. b) Accumulated V brightness phased with the ephemeris in Eq. (1), and separated into the corresponding time intervals in panel a).
Open with DEXTER

Spots typically appear at latitudes higher than ≈30°, except for a few cases where an extension of a larger spot may bridge that latitude (e.g., in 2010, sets 5 and 9). The longitudinal concentration, on the other hand, had changed abruptly at one point. While in 2009, spots concentrated around zero phase (φ = 0.00; facing the observer at inferior conjunction), in the following season in 2010 the main concentration occurred near φ = 0.50 and remained so into 2011 until the end of our data set. We focus on this change in Sect. 5.2, but first we describe the results from the individual observing seasons.

2009 observing season. We have complete phase coverage for each image reconstructed from 2009. In the first image (set 1), we reconstruct three individual spots with latitudes between 45° and 75° and separated by ~0.40 in phase. The second image (set 2) shows the larger of the three spots at ~60° latitude and at φ = 0.05. Two of the spots seem to be marginally connected with each other, and this region appears to evolve systematically. Because there are almost four rotations between sets 1 and 2, some spot evolution is even expected. As time progresses from set 2 to set 4 (i.e., until the last image in 2009), the main spot at φ = 0.05 still exists but appears to have drifted toward decreasing phases with respect to the epoch and orbital period in Eq. (1). Spot migration like this is also typical for the shifts expected from differential rotation. One of the two originally connected spots had continuously drifted toward the larger spot. It had completely disappeared by the end of the season.

2010 observing season. In 2010 the phase coverage is not as good as in 2009, but still quite good for reconstructing reliable surface maps. The time gap between successive maps is typically a few days, except between sets 7 and 8, where it is ten days. Nevertheless, spot evolution can be traced safely through six successive stellar rotations. In the beginning of 2010 (set 5), we see the main spot at φ = 0.83 and two weaker spots at φ = 0.50 and φ = 0.18 at lower latitudes than the main spot. As time progresses, the spot at φ = 0.18 fades, while the spot at φ = 0.50 cools and becomes the most dominant feature. The spot at φ = 0.83 still survives but warms up and eventually approaches the main spot at φ = 0.50. We also observe that the two spots are connected to each other with a pale structure around ~75° latitude that is slightly cooler than the photosphere. All spots appear paler in set 7, while these spots appeared darker in the previous and in the following sets. This is most likely related to the data quality of set 7, where the average S/N is just around 120, while it is usually between 170–200. Moreover, the visible surface from sets 9 and 10 at φ = 0.00 appears to be almost unspotted except for the high-latitude features close to the rim of the visible surface. We emphasize that sets 9 and 10 have considerable phase gaps (listed in Table 5). We have no observation between φ = 0.75 and φ = 0.25 in set 10 and φ = 0.70 and φ = 0.00 in set 9. The consequence is that whatever weak structure may have existed at these longitudes, it cannot be reconstructed. The sudden disappearance of spots within these phases is related to these phase gaps and not to spot decay, for example.

2011 observing season. We have three maps for 2011 covering four consecutive stellar rotations. There is a time gap of one rotation between the end of set 11 and the beginning of set 12. We note that set 11 comprises only seven spectra with a maximum phase gap between two successive observations of 0.38 between φ = 0.28 and φ = 0.66. Care must be taken to avoid overinterpreting spots at these phases. Set 11 shows one cool, dominant spot at φ = 0.40 and two warmer spots at φ = 0.70 and φ = 0.00. After two rotations (i.e., in set 12), the main spot appears to have broken into two smaller and warmer spots, while the spot at φ = 0.70 became the most dominant spot. However, this evolution scenario is hampered by the phase gaps in set 11 and is not conclusive. The spot at φ = 0.00 seems to have moved toward decreased phases by the time of set 12 and then has vanished in set 13. At the same time, a small appendage of the originally cooler spot has migrated to φ = 0.10.

5.2. Active longitudes

In Paper I we showed that an activity cycle repeats almost regularly in time. A preliminary 6.28-yr activity cycle was reconstructed from long-term photometry. It predicted a new activity cycle to begin around 2009/2010. In the previous section, we noted that the main spot activity was located around φ = 0.00 in 2009, while it was around φ = 0.50 in 2010 and 2011. There were also almost no spots at φ = 0.50 in the 2009 maps, while the main spot concentration was observed at exactly that phase in 2010 and 2011. This means that we apparently witnessed a flip-flop event. The new longitudes of the main activity zone appear to correspond to either the sub-stellar point (i.e., φ = 0.00) or its antipode (φ = 0.50). We interpret this as evidence for two persistent active longitudes separated by 180° on the surface of HD 208472, which might be related with the binary nature of the star.

We also applied the same two-spot model technique as in Paper I to our recent V-band photometry from 2010–2015. In this way, we obtained the longitudinal distribution of spots, but over a longer period in time than the Doppler imaging. For details of the technique we refer to Paper I and references therein. The results are shown in Fig. 9a, which is the updated version of the Fig. 10 in Paper I. The spot positions show the continuing systematic drifts and verify the beginning of the new activity cycle around 2009/2010. Furthermore, the most recent photometric data show some evidence for yet another phase jump in 2015, in agreement with the prediction from the ~6 yr cycle length. The light curves in the orbital phase frame are shown in Fig. 9b. The active longitudes clearly show the alternating locations of minima, but the figure also emphasizes the existence of the ever increasing brightness of the system (see Fig. 1).

5.3. Differential rotation

We traced the position of spots by cross-correlating successive Doppler maps in longitudinal direction and examining the amount of longitudinal shift in a given latitude range where spots occur (e.g., Kővári et al. 2015; Künstler et al. 2015). This is basically equivalent to tracking the motion of sunspots as applied to solar-disk images (e.g., Howard 1984; Brajša et al. 2002).

In our present application, we are able to produce three cross-correlation-function (ccf) maps for the observing season 2009, five ccf maps for 2010, and two ccf maps for 2011. Altogether, we have ten ccf maps at our disposal. Because the Doppler images are based on a 5° × 5° equal-degree segmentation of the surface, we also calculated the cross correlation for latitudinal bins of 5° width. A Gaussian profile was used for fitting and identifying the CCF peaks. We present the grand average ccf map in Fig. 10a, which was obtained by averaging all ten ccf maps from three consecutive years. It shows that the available spot-displacement information extends from the visible pole down to ≈10° latitude. Thus, HD 208472 provides a wide latitude range with spots and indeed allows for a good estimation of its differential rotation. In Fig. 10b we plot the positions of the ccf peaks along with their estimated errors defined as the full-width at half maximum of the respective Gaussian fits. We then fit a solar-like differential rotation law to the global ccf pattern of the star. Considering the angular velocities at a latitude b and at the equator (Ω(b) and Ωeq, respectively), and the difference between the angular velocities at the pole and the equator (ΔΩ = Ωeq−Ωpole), we adopt a typical differential rotation law as (2)where the surface shear parameter, α, is defined by α = ΔΩ / Ωeq. Using Eq. (2) to fit the differential rotation pattern in Fig. 10b, we find α = + 0.015 ± 0.003 and Ωeq = 0.2801 ± 0.0004 rad/day, which corresponds to an equatorial period Peq = 22.431 ± 0.031 d.

A better representation of the differential rotation pattern can be achieved by adding another term to the differential rotation law, which is commonly adopted for the Sun, (3)where Ωpole is given by the sum of Ωeq, Ω1, and Ω2. In this case, we find a shear parameter α of + 0.012 ± 0.008, Ωeq = 0.2777 ± 0.0006 rad/day and Peq = 22.628 ± 0.048 d, where Ω1 = 0.0101 ± 0.0016 rad/day and Ω2 = −0.0134 ± 0.0015 rad/day.

Both obtained values of Peq differ by more than their individual errors. This discrepancy is caused by the lack of data points below a latitude of ±10° deg. The two α values agree to within their errors and suggest weak but solar-type differential rotation with a surface shear about one order of magnitude weaker than for the Sun. Such weak differential rotation is a common result for RS CVn binary stars with late-type cool giants, such as for XX Tri (Künstler et al. 2015) or ζ And (Kővári et al. 2012), for which both studies were based on time-series Doppler imaging.

thumbnail Fig. 10

a) Average ccf map from ten seasonal ccf maps. b) Average differential rotation pattern. The dashed line denotes the best fit using the single-term law in Eq. (2), while the solid line shows the best fit from the double-term law in Eq. (3).
Open with DEXTER

6. Summary and conclusions

Three years of high-resolution time-series STELLA spectra enabled us to not only refine precise spectroscopic orbital elements and atmospheric parameters of HD 208472, but reconstruct 13 separate surface images from three observing seasons through multiple-line inversions. Surface reconstructions were made by tracing the distortions in the absorption line profiles with time, which is the main principle of our time-series approach of Doppler imaging. The line-profile distortions are also the cause of systematic shifts in radial velocities, known as spot jitter in case of chromospherically active stars. After proper removal, it decreased the RV residuals for our new spectroscopic orbital solution by a factor of four from 380 m s-1 to 88 m s-1. Fourier analysis of the RV residuals indicated that the dominant periods are always very close to the orbital period and its harmonics, which is the typical signature of spot jitter.

HD 208472 exhibits active longitudes on adjacent hemispheres and occasional flip-flops between them that overlap the continuous drift in longitude with respect to the orbital reference frame. Similar active longitudes and flip-flops were observed in several RS CVn binaries (Jetsu 1996; Berdyugina & Tuominen 1998; Künstler et al. 2015), FK Com-type stars (Jetsu et al. 1991, 1999), and young solar-like stars (Berdyugina & Järvinen 2005). There is some evidence that flip-flops occur more or less regularly with timescales of a few years up to a decade, sometimes referred to as flip-flop cycle (see the references above), but the evidence remains inconclusive. To detect these cycle lengths, we still need long-term photometry with phase-resolved sampling or, even better, modern time-series Doppler imaging.

A 6.28 yr flip-flop cycle was discovered for HD 208472 in our Paper I, where we predicted the begin of a new cycle for 2009/2010. Our Doppler imagery in the present paper caught this event. We witnessed the change of the main activity zone from φ = 0.00 to φ = 0.50 between the 2009 and 2010 observing seasons over a time range of approximately one year. Thus, the location of the main activity zone corresponds to the sub-stellar point and its antipode, as has been suggested by Oláh (2006). Theoretically, coexistence of axisymmetric and non-axisymmetric dynamo modes may produce active longitudes and also give rise to the flip-flop phenomenon (Fluri & Berdyugina 2004; Moss 2005; Elstner & Korhonen 2005). However, why these are aligned and influenced by the secondary star and why the brightness continues to increase over several activity cycle is yet to be understood.

From cross-correlating successive Doppler maps, we found a surface shear of α = + 0.015 ± 0.003 for a sin2-latitude differential rotation law and α = + 0.012 ± 0.008 for a sin2 + sin4-latitude differential rotation law. This differential rotation is in solar-like direction, meaning that poles rotate slower than the equator, but ≈15 times weaker than on the Sun. The average value is close to the α value of XX Tri (Künstler et al. 2015), which is a similar RS CVn system in terms of physical properties and spectroscopic orbit, but maybe about twice as old as HD 208472. Comparing other stars whose α values were determined from Doppler imaging, our target star belongs to a group of weak solar-type differential rotators (AB Dor, Donati & Collier Cameron 1997; ζ And, Kővári et al. 2012; IL Hya, Kővári et al. 2014). On the other hand, anti-solar differential rotation was found (or at least claimed) in some other active stars, such as σ Gem (Kővári et al. 2015), UZ Lib (Vida et al. 2007), HU Vir (Strassmeier 1994; Harutyunyan et al. 2016), and HD 31993 (Strassmeier et al. 2003). Tidal force in a binary system is considered as a possible explanation for the existence of two different directions of differential rotation, depending on the orbital and physical properties of the related system. Holzwarth & Schüssler (2002) showed that the tidal force in an RS CVn binary may lead to preferred longitudes. It may also affect the angular momentum distribution in the convective envelope of the giant component and alter its differential rotation pattern (Kővári et al. 2012). Another effect of the tidal force on differential rotation might be the suppression of the strength of the differential rotation through tidal locking (Collier Cameron 2007). The RS CVn binaries mentioned above have α values of just a few of per cent (solar or anti-solar). In comparison, the single K2 giant HD 31993 (Strassmeier et al. 2003), where no tidal interaction is present, has a dramatically higher shear value of α = −0.125. Another comparison can be made with the single K1 giant DP CVn (Kővári et al. 2013), which has a surface shear value of α = −0.035, which is a stronger surface shear in magnitude than HD 208472. This shows that fast-rotating single giants with α values even comparable to the solar value have been reported, while no such RS CVn stars have been found. This supports the picture that differential rotation is suppressed by tidal forces as suggested by Scharlemann (1982), for example.

1

The Image Reduction and Analysis Facility is hosted by the National Optical Astronomy Observatories in Tucson, Arizona at URL iraf.noao.edu.

Acknowledgments

O.Ö. thanks the Scientific and Technological Research Council of Turkey (TÜBİTAK) for funding this study under “2219 post-doc research scholarship program”. O.Ö. appreciates friendly and warm hospitality of AIP Stellar Activity Group during his post-doc research at AIP. K.G.S. thanks the State of Brandenburg and the German BMBF for the continuous funding of the STELLA facility in Tenerife and the APTs in southern Arizona.

References

Appendix A: Line profiles of Doppler images

Figure A.1 shows observed and inverted line profiles of each Doppler image given in Figs. 6 to 8.

thumbnail Fig. A.1

Observed (dots) and inverted (continuous curve) line profiles of 13 sets. For a given set, the corresponding orbital phase of each line profile is shown at the right side of the related panel. Observing season and set numbers are shown at the top of the panels.
Open with DEXTER

All Tables

Table 1

Periodogram results of radial velocity residuals.

In the text
Table 2

Spectroscopic orbital elements.

In the text
Table 3

Astrophysical parameters of HD 208472. Errors quoted are internal errors.

In the text
Table 4

Absorption lines used for the mean profile. Atomic data are taken from VALD-3.

In the text
Table 5

Observing log of Doppler-imaging sets.

In the text

All Figures

thumbnail Fig. 1

a) Long-term V-band photometry of HD 208472. The dashed line denotes the updated unspotted brightness level. The blue filled dots show data from Paper I, while the red dots denote the recent data. Plus signs show differential magnitudes in the sense of check-minus-comparison star. These measurements are shifted by a constant amount and plotted on a common scale with measurements of the variable star. b) Minimum brightness vs. V amplitude diagram. The straight line is a linear representation of the upper envelope of the scatter.
Open with DEXTER
In the text
thumbnail Fig. 2

a) Radial velocities (dots) before spot correction and the representation of the orbital fit (line). The bottom panel shows the residuals. b) Spot-corrected radial velocities and the final fit for the orbital elements. The rms of the fit decreased by a factor four with respect to the uncorrected velocities.
Open with DEXTER
In the text
thumbnail Fig. 3

PARSES spectrum-synthesis results of four atmospheric parameters for each individual SES observation.
Open with DEXTER
In the text
thumbnail Fig. 4

Position of HD 208472 in the H-R diagram. The tracks (black dashed lines) and isochrones (red continuous lines) are adopted from Bertelli et al. (2008) for Y = 0.30 and Z = 0.008.
Open with DEXTER
In the text
thumbnail Fig. 5

Quality-of-fit (χ2) for a series of trial vsini values. The red dashed curve is a quadratic fit to the data points. The arrow indicates the position of the best-fit vsini of 19.4 km s-1.
Open with DEXTER
In the text
thumbnail Fig. 6

Doppler images of HD 208472 for the 2009 observing season.
Open with DEXTER
In the text
thumbnail Fig. 7

Doppler images of HD 208472 for the 2010 observing season.
Open with DEXTER
In the text
thumbnail Fig. 8

Doppler images of HD 208472 for the 2011 observing season.
Open with DEXTER
In the text
thumbnail Fig. 9

a) Longitudinal distribution of spots on HD 208472 from photometric two-spot models. The shaded regions are the time ranges where spots change their preferential positions (hemispheres) on the surface of the star. The continuous line represents the longitudes of the apsidal line of the binary system, while the two dashed lines restrict the stellar surface. The vertical dot-dashed lines correspond to times when changes occurred in light-curve amplitude and photometric period. We overplot the longitudes of the main spot of the Doppler images with cross symbols in red. b) Accumulated V brightness phased with the ephemeris in Eq. (1), and separated into the corresponding time intervals in panel a).
Open with DEXTER
In the text
thumbnail Fig. 10

a) Average ccf map from ten seasonal ccf maps. b) Average differential rotation pattern. The dashed line denotes the best fit using the single-term law in Eq. (2), while the solid line shows the best fit from the double-term law in Eq. (3).
Open with DEXTER
In the text
thumbnail Fig. A.1

Observed (dots) and inverted (continuous curve) line profiles of 13 sets. For a given set, the corresponding orbital phase of each line profile is shown at the right side of the related panel. Observing season and set numbers are shown at the top of the panels.
Open with DEXTER
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.