Free Access
Issue
A&A
Volume 598, February 2017
Article Number A74
Number of page(s) 21
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201629814
Published online 02 February 2017

© ESO, 2017

1. Introduction

Slowly pulsating B (SPB) stars are main-sequence stars of between B3 and B9 in spectral type, 11 000 K and 22 000 K in effective temperature, and 2.5 ℳ and 8 ℳ in mass, which show non-radial, heat-driven, gravity-mode oscillations (see e.g. Aerts et al. 2010, Chap. 2). Although it has been a quarter of a century since their discovery by Waelkens (1991), in-depth seismic modelling based on high-order g modes was only achieved recently for KIC 10526294 by Pápics et al. (2014). Even though SPB stars are not massive according to the classical definition of supernova progenitors, they share the same internal structure on the main sequence by having a convective core and a radiative envelope. This means that SPB stars are ideal asteroseismic probes of ill-understood internal mixing processes that have a significant influence on the lifetime of the metal factories of the Universe, such as core overshooting, diffusive mixing, or internal differential rotation. Therefore these stars hold the key to a precise calibration of stellar structure and evolution models of massive stars. Such calibrations are much needed because the observed distribution of massive stars does not agree with theoretical predictions, as shown by Castro et al. (2014), and internal mixing is likely to be one of the possible solutions.

Before the space revolution of asteroseismology began with the launch of the MOST (Microvariablity and Oscillations of STars; Walker et al. 2003) and CoRoT (Convection Rotation and planetary Transits; Auvergne et al. 2009) missions, it was very difficult to collect sufficient photometric data of SPB stars, mainly because their dominant oscillation modes have typical periods of 0.5–3 days (see e.g. De Cat & Aerts 2002), and also their mode amplitudes are relatively low (e.g. compared to the slightly hotter β Cephei stars that pulsate in p modes). While some extremely time- and resource-intensive multisite campaigns and the CoRoT observations have resulted in a few in-depth observational studies showing an unexpected richness in the pulsation behaviour of OB stars in general (for a historical overview see e.g. Pápics 2013), the number of stars with actual constraints on the extent of the core overshoot and/or the internal rotation profile remained very limited (Aerts 2013). The real breakthrough arrived with the Kepler Mission, thanks to its monitoring over four years.

The primary goal of Kepler was the detection of Earth-like exoplanets (Borucki et al. 2010), but in the process the initial mission also provided us with four years of virtually uninterrupted photometric time series of more than 150 000 stars with micromagnitude precision in a 105-square-degree field of view (FOV) fixed between the constellations of Cygnus and Lyra, which are a goldmine for asteroseismology (Gilliland et al. 2010). While the positioning of the FOV above the Galactic plane was not optimal in terms of the number of observable OB stars (Debosscher et al. 2011; Balona et al. 2011; McNamara et al. 2012), the long time base and high precision meant that previously inaccessible seismic diagnostic methods could be applied with success for the first time.

Forward seismic modelling is not possible without mode identification. Earlier this could be only achieved for a few selected modes from ground-based, follow-up, multicolour photometry and/or line-profile variation studies (e.g. De Cat et al. 2005); but now the full process can be reduced to a relatively simple pattern search in the periodogram combined with the approximate knowledge of fundamental parameters from spectroscopy. It has been known for a long time that there is structure in the dense frequency spectrum of SPB stars (Dziembowski et al. 1993; Gautschy & Saio 1993). A series of g modes of the same degree with consecutive radial orders n is not only predicted to be nearly equally spaced in period, but the periodic deviations from such an equal spacing carry information about the physical properties in the near-core region (Miglio et al. 2008). The first detections of g-mode period series came from CoRoT data of the hybrid SPB/β Cep pulsators HD 50230 (Degroote et al. 2010, 2012) and HD 43317 (Pápics et al. 2012), but these were not precise enough to provide input for in-depth modelling.

This changed with Kepler and the case of the extremely slow rotator KIC 10526294. Aided especially by the tenfold increase in frequency resolution compared to CoRoT, Pápics et al. (2014) managed to identify a clear series formed by the rotationally split multiplets of the 19 consecutive (in radial order) dipole modes of the star that are quasi-equally spaced. Based on the unprecedented observational data, Moravveji et al. (2015) found that a diffusive core-overshooting prescription that is exponentially decaying fits better than a step function formulation, and these authors also put stringent constraints on the value of the core overshoot; they derived fov ∈ [0.017,0.018] reaching a ~5% precision in fov compared to the typical ~50% achieved before. They showed that the inclusion of extra diffusive mixing (log D ∈ [1.75,2.00] cm2 s-1) in the radiative envelope is needed to explain the structure in the period spacing pattern. The relative difference between the observed and their model frequencies is better than 0.5%. As a result of having such a good seismic model, Triana et al. (2015) have shown that, by using different inversion techniques, it is possible to derive the internal rotation profile of an SPB star using rotationally-split multiplets of different modes that penetrate to different depths of the star. In this particular case, while the seismic data are consistent with rigid rotation, a profile with counterrotation within the envelope has a statistical advantage over constant rotation. This finding was immediately validated by the 2D numerical simulations of Rogers (2015). Such an inversion has never been performed before based only on gravity modes, and even overall there were only a few earlier studies of hydrogen-burning stars (besides the Sun) that resulted in constraints on the internal rotation profile (see Aerts et al. 2003; Pamyatnykh et al. 2004; Briquet et al. 2007; Kurtz et al. 2014). In the meantime a faster rotator analogue of this star was also found. KIC 7760680 exhibits a rotationally affected series, consisting of 36 consecutive high-order, sectoral dipole gravity modes (Pápics et al. 2015), which carries clear signatures of chemical mixing and rotation. This star was modelled thoroughly by Moravveji et al. (2016), who reconfirmed that the exponential overshoot prescription on top of the receding convective core is superior to the step function formulation. They constrained the overshooting parameter to fov = 0.024 ± 0.001, which is slightly larger than the value in its virtually non-rotating counterpart KIC 10526294. Furthermore, an extra diffusive mixing of log D = 0.75 ± 0.25 cm2 s-1 was necessary to explain the deviations of the period spacing pattern from the asymptotic value.

In this paper we present our discovery of five additional SPB stars from the initial Kepler Mission that show rotationally affected period series. Besides analysing the remaining four targets of our original Guest Observer (GO) sample of eight stars, we also revisit one of the primary components of the binary stars already presented in Pápics et al. (2013). We refer to the latter for details on the target selection. Future modelling based on this observational data will significantly extend the list of SPB stars that have well-constrained internal mixing parameters (Aerts 2015), and bring us closer to the seismic calibration of stellar structure and evolution models of massive stars.

2. Spectroscopy

2.1. Spectroscopic data

We gathered high-resolution spectra of our targets between 2010 and 2015 to confirm the initial classification of the SPB stars, obtain precise fundamental parameters, test for their binary nature, and cover the binary orbit if necessary. Depending on the apparent brightness, we used the Hermes spectrograph (Raskin et al. 2011) installed on the 1.2 m Mercator telescope or the Isis spectrograph mounted on the 4.2 m William Herschel Telescope, which are both located on La Palma (Spain). We refer to Table 1 for a summary of the spectroscopic observations of the new SPB stars.

Table 1

Logbook of the spectroscopic observations – grouped by stars and observing runs.

The raw Hermes exposures were reduced using the dedicated pipeline of the instrument, resulting in the extracted, cosmic-removed, log-resampled (which includes barycentric correction), order-merged spectra that we used afterwards. The Isis frames were reduced using standard STARLINK routines following the optimal extraction described by Marsh (1989). All spectra were rectified using the method described in Pápics et al. (2013).

We used the method of spectral disentangling (Simon & Sturm 1994) formulated in Fourier space (Hadrava 1995) and implemented in the FDBinary code (Ilijic et al. 2004) to separate the individual spectral contributions for the KIC 4930889 binary star (see below for details). The method bypasses the step of radial velocity (RV) determination but simultaneously optimises orbital elements of the system and individual spectra of the binary components instead. We start with the disentangling of several (typically 4 or 5) spectral regions where the contributions from both binary components can be clearly visible and determine orbital elements from each of the considered regions separately. The final set of the orbital parameters is then computed by taking average values from the above-mentioned individual solution, and errors bars are computed as the standard deviations of the mean. Once computed, the orbital elements are fixed and the disentangling is performed in the entire wavelength range, including the Balmer lines.

The spectra of all stars in this paper were analysed using the Grid Search in Stellar Parameters (GSSP; Tkachenko 2015) software package1. The method is based on a comparison of the observed spectra in an arbitrary wavelength range with a grid of synthetic spectra in six fundamental parameters (Teff, log g, micro- and macroturbulent velocities, vsini, and metallicity) and optionally light dilution factor (specifically for the disentangled spectra of binary stars). The latter is assumed to be constant over considered wavelength interval, which is justified in our case, given similar spectral types of the individual binary components. Synthetic spectra are computed using the SynthV spectrum synthesis code (Tsymbal 1996) based on the grid of atmosphere models pre-computed with the most recent version of the LLmodels code (Shulyak et al. 2004). Both codes operate in LTE mode and make use of the atomic data extracted from the Vienna Atomic Line Database (VALD; Kupka et al. 2000). Chi-square merit function is used to evaluate the quality of the fit between the observed spectrum and each of the synthetic spectra in a grid; the error bars on the fundamental parameters are computed from chi-square statistics taking all possible correlations between the parameters in question into account. Tkachenko (2015) provides a detailed description of the method and full capabilities of the GSSP code.

In the following subsections we give the basic observable parameters and discuss the results from the spectral analyses for all targets individually. Pápics et al. (Table 5, 2013) provide the results for KIC 6352430 A. A general overview of the position of the stars in the SPB instability strip among other CoRoT and Kepler targets is given in Fig. 1.

thumbnail Fig. 1

Kiel diagram of a sample of B-type stars near the main sequence, for which an in-depth seismic data analysis is available, observed by CoRoT (circles) and Kepler (squares). The new SPBs from this paper are plotted with diamonds. Stars analysed or revisited in this paper are labelled using their full names, while for the others we use the following annotations: C1 – HD 48977, C2 – HD 43317, C3 – HD 50230 A, C4 – HD 50846 A, C5 – HD 182198, C6 – HD 181440, C7 – HD 46179, C8 – HD 174648, C9 – HD 174884 AB, K1 – KIC 4931738 AB, K2 – KIC 10526294, and K3 – KIC 7760680. The dot-dashed line represents the zero-age main sequence (ZAMS), while the dashed line indicates the terminal-age main sequence (TAMS). The cool edge of the β Cep instability strip is plotted with a dotted line (for the = 0 modes). The thin grey solid lines denote evolutionary tracks for selected masses. Isochrones from 20 to 28 Myr and isoradii corresponding to integer multiples of the Solar radius are plotted using thin dashed and dotted grey lines, respectively. The instability strip of SPB stars is defined as the region where the sum of excited = 1 plus = 2 gravity modes is at least three, and it is shown in colours ranging from purple (3 excited modes) to bright yellow (63 excited modes). All model data were taken from Moravveji (2016) for Z = 0.014, the Asplund et al. (2009) mixture, an exponential core overshoot of fov = 0.02, and opacity enhancement factors for iron and nickel of βFe = 1.75 and βNi = 1.75. Different error bars reflect differences in data quality and methodological approach.

2.1.1. KIC 3459297

The fundamental parameters from the spectral synthesis of KIC 3459297 are listed in Table 2, while we show the quality of the fit between the synthetic and observed spectra in Fig. 2. The derived parameters suggest that KIC 3459297 is the most evolved of the four new SPBs of this paper.

Table 2

Fundamental parameters and basic observable properties of KIC 3459297.

thumbnail Fig. 2

Comparison of the rectified observed spectrum and the synthetic spectrum for the two observed wavelength regions of KIC 3459297. In each panel, the Isis spectrum is plotted with a grey solid line and the synthetic spectrum is plotted with a red solid line.

2.1.2. KIC 4930889

KIC 4930889 turned out to be a double-lined spectroscopic binary (SB2), thus we took a total of 26 exposures with Hermes to cover the binary orbit. The SB2 nature of the original spectra is visible in Fig. 3. The derived orbital parameters from the disentangling process are listed in Table 3, while the reconstructed orbit and radial velocities are shown in Fig. 4. The fundamental parameters, obtained via spectral synthesis from the disentangled spectra place both components in the SPB instability strip, are listed in Table 3, while the fits of the disentangled spectra are shown in Fig. 5.

Table 3

Fundamental parameters, orbital parameters, and basic observable properties of KIC 4930889.

2.1.3. KIC 9020774

Its low brightness (with Kepler mag. = 15.091, the dimmest source in our GO sample) limited the precision of our analysis, but the fundamental parameters derived from the spectral synthesis (see Table 4, and Fig. 6) strongly suggest that KIC 9020774 is a young B star.

2.1.4. KIC 11971405

Spectroscopy of KIC 11971405 revealed that the star is a fast rotator near the middle of the SPB instability strip (see Table 5, and Fig. 7) with a projected rotational velocity of ~240 km s-1, which is the highest from the stars in our GO sample. Furthermore, the stellar spectrum shows very weak emission features near the core of the Hα line (6563 Å), and changes over time (see Fig. 8). The emission component in Hα is extremely weak; it appears only as two small bumps inside the absorption line around the core and it is not visible in any other Balmer line. Nevertheless, this grants the spectral classification of Be (for a review, see Rivinius et al. 2013). We present, in Sect. 3.4.5, that the star showed small photometric outbursts during the nominal Kepler Mission, which strengthens this classification. The exposures taken in October 2010 (before the start of the photometric outbursts) show excess absorption around the Mg ii lines at 4481 Å and 4534 Å compared to the observations from 2015 (after the start of the outbursts). These features consist of a broader and narrower absorption component. We checked whether these features could be from external contamination (e.g. of telluric origin), but the examination of other exposures from the same and neighbouring nights did not show anything unexpected at these wavelengths. Similar sharp red-shifted transient absorption features were observed in the Be star ω CMa (Štefl et al. 2003), but in the outburst phase. In addition, the emission features in Hα also show some minor changes between the two epochs.

thumbnail Fig. 3

Phase-folded residual spectra of KIC 4930889 around the H i line at 4861 Å obtained after subtraction of the average. The colour scale represents the residual intensity at each wavelength point. White horizontal regions denote phases that are more than 0.05π away from an observed spectrum.

thumbnail Fig. 4

Spectroscopic orbital solution for KIC 4930889 (solid line: primary; dashed line: secondary) and the radial velocities from the disentangling process. Filled circles indicate RVs of the primary; filled triangles denote RVs of the secondary component.

Table 4

Fundamental parameters and basic observable properties of KIC 9020774.

3. The Kepler photometry

3.1. Description and reduction of the data

For the photometric analyses presented later in this paper we used all available Kepler data from the nominal mission. These stars were observed in long cadence mode (LC) with a cadence of 29.43 min, for ~4 yr between May 2009 and May 2013 (from the commissioning quarter Q0 to the end of the last science quarter Q17). There was no data taken of KIC 6352430 A and KIC 9020774 during Q0, and because of the loss of module 3 of the CCD array, there was also no photometry gathered during Q6, Q10, and Q14 of KIC 3459297 and KIC 4930889. We used custom masks (see example in Fig. 9) to construct quarterly light curves from the target pixel files (Data Release 24) to minimise instrumental trends by adding pixels around the default masks with significant flux in them, and to avoid possible contamination from nearby objects by removing pixels containing light from sources other than the target. We cleaned all quarters from clear outliers, then detrended them using a division with a low-order (typically 2nd order) polynomial. In some specific cases quarters were broken up into pieces for this if the signal levels before and after an anomaly (e.g. the large coronal mass ejection during Q12) were visibly different. Afterwards the detrended quarters were merged to a continuous light curve and the counts were converted to ppm (see Fig. 10 for an overview of the final light curves).

thumbnail Fig. 5

Comparison of the Hermes spectrum at the phase of the largest observed radial velocity difference, disentangled component spectra, and synthetic spectra for each component for selected wavelength regions of KIC 4930889. In each panel, the Hermes spectrum is plotted with a black solid line, the disentangled spectra are plotted with grey solid lines, and synthetic spectra (scaled to the disentangled spectra) are plotted with red solid lines. For clarity, the primary and secondary spectra are plotted with a vertical offset of + 0.1 and + 0.2 flux units, respectively.

thumbnail Fig. 6

Comparison of the synthetic spectrum and the average of the rectified observed spectra for the two observed wavelength regions of KIC 9020774. In each panel, the Isis spectrum is plotted with a grey solid line, and the synthetic spectrum is plotted with a red solid line.

Table 5

Fundamental parameters and basic observable properties of KIC 11971405.

thumbnail Fig. 7

Comparison of the median of the rectified observed spectra, and the synthetic spectrum for selected wavelength regions of KIC 11971405. In each panel, the Hermes spectrum is plotted with a grey solid line, and the synthetic spectrum is plotted with a red solid line.

3.2. Frequency analysis

To extract the pulsation modes (and describe any other periodic light variation), we followed a standard prewhitening procedure similar to that described in Pápics et al. (2014), with the small modifications that a) the prewhitening was carried out in the order of peak significance (instead of amplitude); and b) the prewhitening procedure was stopped when the signal-to-noise ratio (S/N; calculated in a 1 d-1 window) of the last prewhitened peak dropped below four (instead of using a p value in hypothesis testing and filtering the frequency set later based on the measured S/N values). This resulted in the set of Fourier parameters (Aj amplitudes, fj frequencies, and θj phases) that we used later for the analysis. An overview of the power spectra is shown in Fig. 10.

3.3. Searching for patterns

There are two different kind of frequency patterns that we want to distinguish when dealing with SPB stars: combination peaks that are often found in groups and modes of the same degree forming a period series (of consecutive radial orders n). The former are a nuisance for the particular purpose of this paper because they complicate the identification of independent modes. However, combination frequencies could be used in the future to study nonlinear coupling of the modes, despite current uncertainties for their astrophysical implications in upper main-sequence stars; see e.g. Kurtz et al. (2015). Frequencies of modes with the same degree and consecutive radial order, on the other hand, depend explicitly on the internal structure and rotation frequency of the star, hence carry information about the physical properties inside the star (as already discussed in Sect. 1); therefore it is very important to separate peaks that belong to these two very different groups (see e.g. Pápics 2012). Combination peaks are harmonics or linear combinations of the individual frequencies, but being able to explain a peak as a linear combination of other peaks does not imply that this peak cannot be an actual pulsation mode.

thumbnail Fig. 8

Temporal changes in the spectrum of KIC 11971405. The average of the two consecutive exposures from October 2010 is plotted on the top using a black solid line, while the average spectrum from June 2015 is shown below with a red solid line with an offset of −0.15 in flux units. For clarity, the average from 2010 is repeated in the background using a grey solid line. The sharp absorption features in the bottom panel are telluric lines.

thumbnail Fig. 9

Custom pixel mask for the Q8 data of KIC 4930889. Individual pixel-level light curves are plotted for each pixel that was downloaded from the spacecraft. The numerical values and the background colour indicate the flux levels in each pixel. The pixels with green and red borders were used to extract the standard Kepler light curves. We added yellow pixels (significant signal is present) and removed red pixels (non-significant signal or contamination from nearby stars) in our custom mask for the light curve extraction. This results in a light curve with significantly less instrumental effects and/or contamination than the standard extraction.

Table 6

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the main period series of KIC 3459297.

This remark gets especially important when searches for combinations are extended to higher orders (\hbox{$\mathcal{O}>2$} in Pápics 2012) allowing positive and negative combinations simultaneously (e.g. cifi + cjfjckfk, where c ∈ N+, referred to as a mixed-type combination hereafter). Indeed, even with the frequency precision of the four years of Kepler data, true individual modes might easily get matched as combinations in a blind search. The main reason behind this is that the frequency spectrum of SPB stars is very dense. For example, the slow rotator KIC 10526294 shows a series of 19 = 1 modes (56 counting the rotationally split components) within 0.6 d-1 (Pápics et al. 2014), and the moderate rotator KIC 7760680 shows a series of 36 = 1 modes within 0.5 d-1 (Pápics et al. 2015); but, most importantly, these observed values are all in agreement with numbers from excitation calculations (Moravveji 2016). To complicate matters further, rotation introduces rotationally split components (and distorts period spacing patterns; see e.g. Bouabid et al. 2013) contributing to an even denser pulsation spectrum. Combination frequencies can most often explain the frequency groups found at higher frequencies (Kurtz et al. 2015), but one should pay special attention during a combination search because, with an ill-considered parameter selection, it is possible to match real frequencies in the main group of gravity modes as a mixed-type combination of a few selected peaks within the same group just by chance. That is why we advise against the use of mixed-type combinations during the search. In our experience, real combination peaks are in most cases simple harmonics (cifi), or very simple positive or negative combinations (cifi ± cjfj) of only two strong parent frequencies. It is extremely rare that there is no better physical explanation for a peak than a mixed-type combination.

After we filter out the true combination peaks, we attempt to identify series of modes of the same degree among the remaining frequencies. Since internal mixing and rotation influence the shape of the period spacing function (Miglio et al. 2008; Bouabid et al. 2013), we cannot use the autocorrelation function or histograms with success because the consecutive (in radial order n) modes are not expected to be equally spaced in period anymore. Instead, we have to look for characteristic patterns (rotation-induced slopes, and/or periodic waves) that resemble the expected period spacing function manually.

3.4. Results for individual targets

3.4.1. KIC 3459297

The final processed light curve of KIC 3459297 spans 1470.5 d, contains 52 411 data points, and has a duty cycle of 72.8%. The frequency resolution of the periodogram (also known as the Rayleigh limit) is 1 /T = 0.00068 d-1. The prewhitening procedure resulted in 477 statistically significant model frequencies contributing to a variance reduction of 99.75%, and bringing down the average signal levels from 45.9–27.1–5.0–3.9–3.8 ppm to 4.3–2.3–1.3–1.1–1.0 ppm, measured in 2 d-1 windows centred on 1, 2, 5, 10, and 20 d-1, respectively.

The dominant modes, as expected for an SPB star, are in the g-mode frequency regime between 0.7 and 1.2 d-1, with the dominant pulsation mode at 0.990040(7) d-1, having an amplitude of 4667(88) ppm. The amplitude weighted average of the model frequencies is 1.34 d-1, where the shift to higher frequencies is because of the presence of a large number of combination frequencies.

To confirm this, we follow the same combination search method as in Pápics et al. (2015), restricting the search to combinations of maximum two parent frequencies. We find that all low frequency peaks are simple negative combinations of parent peaks in the main frequency regime, with c ∈ [−1,1]. Almost all peaks towards higher frequencies, however, are simple positive combinations, with c ∈ [0,8], except for a few peaks between 2.2 and 2.8 d-1, some of which stay independent even when three parent frequencies and mixed-type combinations are allowed. Moreover, we find that looking at the histogram of frequency differences in Fig. 11 there are three bins that stand out: 0.990, 1.086, and 0.096 d-1. These match the frequency values of the two highest amplitude pulsation modes and their difference (f1, f2, and f2f1), which implies that the high amplitude parent peaks contribute to a large number of combination frequencies. A similar behaviour has already been detected in KIC 6352430 A (Pápics et al. 2013). Observing these differences in the histogram, for example f2f1, does not require the presence of negative combination frequencies in the periodogram, as two combination frequencies that are formed as f1 + fi and f2 + fi also result in an observed separation of f2f1.

While searching for characteristic patterns of a possible period series, we restrict ourselves to peaks that meet the Loumos & Deeming (1978) criterion, which states that the minimal frequency separation that two close peaks must have to avoid influence on their apparent frequencies in the periodogram is 2.5 × the Rayleigh limit. In practice, given fi and fi + 1 frequencies of consecutive significant peaks, where fi + 1fi < 2.5 /T, we keep the frequency with the larger amplitude and discard the other frequency. We repeat this pairwise test until there are no more consecutive peaks in the periodogram that have a frequency difference below 2.5 × the Rayleigh limit. This lowers the total number of frequencies to 337. This number is significantly reduced further when we remove the simple combination frequencies. Then we take the ~5 highest amplitude peaks and check whether they outline a recognisable period spacing pattern. Afterwards we try to add additional strong frequencies one-by-one, every time checking whether they fit in the period spacing pattern.

We find a long period series (consisting of 43 observed and 2 missing radial orders, which are listed in Table 6) tilted by rotation, which exhibits short period variations on top of the general downwards trend (see Fig. 12). This points towards a star that is already relatively evolved on the main sequence (in good agreement with the observed log g), and requires the absence of a large amount of diffusive mixing that would otherwise smear out the observed sawtooth-like features. Assuming that we observe a prograde l = 1 series (based on the observed spacing, the tilt, and the spectroscopic fundamental parameters), it is possible that the peaks between 2.2 and 2.8 d-1 (listed in Table 7) that are not identified as simple combination frequencies might be members of a prograde (m = 1 or m = 2) l = 2 series (see Fig. 13), but this could only be confirmed after carrying out modelling of the dipole series.

3.4.2. KIC 4930889

The reduced light curve of KIC 4930889 spans 1470.5 days and consists of 52 878 data points. The data has a duty cycle of 73.5%. The Rayleigh limit of the periodogram is 1 /T = 0.00068 d-1. The prewhitening procedure delivered 380 statistically significant model frequencies, resulting in a variance reduction of 99.5% and bringing down the average signal levels from 81.2–66.4–7.5–5.9–5.6 ppm to 12.2–7.1–1.3–0.8–0.5 ppm, measured in 2 d-1 windows centred on 1, 2, 5, 10, and 20 d-1, respectively.

The dominant modes are g modes between 0.75 and 1.4 d-1. The strongest pulsation mode is at 1.240146(6) d-1 with an amplitude of 8093(136) ppm. The amplitude weighted average of the model frequencies is 1.28 d-1. Almost all power above this frequency range comes from combination frequencies (with a possible exception of the relatively strong peak at 1.71770(4) d-1), while a number of low frequency peaks cannot be accounted for as simple negative combinations. Again, as illustrated by the five marked bins (at 0.023, 0.045, 1.194, 1.216, and 1.240 d-1) in Fig. 14, the three dominant modes contribute to a large number of combination frequencies (with the bins corresponding to f1f2,f1f3,f3,f2, and f1, respectively).

Table 7

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the possible l = 2 series in KIC 3459297.

Our manual search for possible period series (after filtering for close frequencies and leaving only 297 model peaks) gives a positive result again. We identify a tilted period series consisting of 20 consecutive radial orders (see Table 8 and Fig. 15). This period series is much smoother than in the case of KIC 3459297, showing only three shallow dips (around 0.82, 0.89, and 0.96 d). These properties point towards an earlier evolutionary stage and possibly a higher diffusive mixing.

Table 8

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the main period series of KIC 4930889.

Since we find an unusually high amount of modes at lower frequencies that cannot be explained as simple combinations, we also searched for period series in this region. This results in two very similar period series, both showing an upwards tilt typical for retrograde modes (see Table 9 and Fig. 16). During the construction of these series we disregarded the results from the combination search, as some peaks might very easily be identified as combination frequencies simply by chance since the region of possible parent frequencies is very dense, and because the fact that none of the higher amplitude peaks in the low frequency region were identified as combinations to start with.

thumbnail Fig. 10

Light curves (plotted using grey dots, showing a zoom-in to slightly less than 1/10th of the full data set) and the power spectra (plotted using red solid lines) of the five SPB stars studied in this paper along with the ones of KIC 10526294 and KIC 7760680 from Pápics et al. (2014, 2015) ordered according to the frequency of maximum power. For further information see Sect. 3.5.

thumbnail Fig. 11

Histogram of the distribution of frequency pairs in the frequency analysis of KIC 3459297, calculated using a bin width of 2.5 /T taking into account all statistically significant frequencies (grey solid line), or only a subset of frequencies where weaker signal from within 2.5 /T of the stronger peaks was excluded (black solid line). The three values that are mentioned in the text are marked with black triangles.

thumbnail Fig. 12

Zoom-in of the period spectrum of KIC 3459297 (top panel) with the members of the period spacing series indicated with vertical dashed lines. Peaks that are discarded during our pattern search are plotted using thin grey lines (for details, we refer to the text). The observed period spacing is shown in the bottom panel.

Table 9

Fourier parameters (periods, frequencies, amplitudes, and phases) of the two extra possible period series of KIC 4930889.

thumbnail Fig. 13

Same as Fig. 12, but for the possible l = 2 series in KIC 3459297.

3.4.3. KIC 6352430 A

The full reduced Kepler photometry of KIC 6352430 A covers 1459.5 d (compared to the 1141.5 days that were available in Pápics et al. 2013) and contains 65 478 data points (with a duty cycle of 91.7%). The Rayleigh limit is 1 /T = 0.000685 d-1. The iterative prewhitening was stopped after the removal of 858 significant frequencies, which resulted in a variance reduction of 99.9%. The average signal levels in the original and residual light curve are 55.8–52.6–6.2–5.1–4.8 ppm and 5.4–3.6–0.9–0.6–0.3 ppm, measured in 2 d-1 windows centred on 1, 2, 5, 10, and 20 d-1, respectively.

thumbnail Fig. 14

Same as Fig. 11, but for KIC 4930889.

The dominant modes are in the g mode regime between 0.9 and 1.8 d-1 with the strongest mode at 1.361716(6) d-1, having an amplitude of 7406(111) ppm. The amplitude weighted average of the model frequencies is 1.84 d-1, where the shift towards higher frequencies is again a clear sign of a large number of combination frequencies. Unlike the previous stars, the transition between the dominant region and the region of combination frequencies is continuous; there are no apparent gaps in the power spectrum where no signal could be found. All peaks below 0.7 d-1 can be explained as simple negative combinations, except for the frequency at 0.03768(3) d-1, which can be identified as the spectroscopically derived binary orbital frequency of fSB2 = 0.03766(3) d-1 from Pápics et al. (2013). We note that the fact that one can identify a frequency as a simple combination does not directly mean that there is no better explanation for the peak in question. In fact, we also observe four harmonics of the orbital frequency up to 5 × fSB2, but some of these get matched as simple negative combinations thanks to the high power density in the g-mode region; however a much more logical explanation, especially given their sequentially smaller amplitudes towards higher harmonics, is that these are orbital harmonics. Above 1.9 d-1 a quickly growing percentage of peaks are identified as combination frequencies, and beyond 2.5 d-1 basically all power can be traced back to parent frequencies in the dominant g-mode region (with a few relatively large amplitude exceptions below 4 d-1). Again, large amplitude modes have a major contribution to combination frequencies, as the strongest bins in Fig. 17 fall at frequencies of the strongest modes and their differences (with the marked bins at 0.102, 0.157, 1.361, and 1.519 d-1 corresponding to f3f1, f2f1, f1, and f2, respectively).

thumbnail Fig. 15

Same as Fig. 12, but for KIC 4930889.

thumbnail Fig. 16

Same as Fig. 12, but for the proposed additional series of the low frequency peaks in KIC 4930889.

thumbnail Fig. 17

Same as Fig. 11, but for KIC 6352430 A.

The purpose of revisiting KIC 6352430 A is to look for tilted period series, and (after filtering for close peaks that left us with 584 model peaks) screening through the independent frequencies in the region of interest, we manage to detect one (see Table 10 and Fig. 18) that is very similar in topography to that of KIC 3459297.

thumbnail Fig. 18

Same as Fig. 12, but for KIC 6352430 A.

Table 10

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the period series detected in KIC 6352430 A.

3.4.4. KIC 9020774

The reduced light curve of KIC 9020774 contains 64 863 data points over 1450.1 d (with a duty cycle of 91.4%). The Rayleigh limit of the data set is 1 /T = 0.00069 d-1. The iterative prewhitening procedure resulted in a model with 281 significant frequencies and a variance reduction of 99.1%. The average signal levels before and after the removal of the model are 41.7–45.3–7.6–6.3–6.0 ppm and 7.0–4.5–3.9–3.8–3.7 ppm, measured in 2 d-1 windows centred on 1, 2, 5, 10, and 20 d-1, respectively.

The dominant modes in the periodogram are g modes between 1.2 and 2.3 d-1. The strongest modes are localised in a tight cluster, with the dominant frequency at 1.900723(3) d-1 with an amplitude of 7672(65) ppm. The amplitude weighted average of all model frequencies is 2.08 d-1. This star also shows the combination frequencies that we have seen in the previous cases, meaning that nearly all peaks below 0.75 d-1 are simple negative combinations (with a possible exception of a stronger peak at 0.381316(7) d-1), and that peaks beyond the g-mode region are simple harmonics and positive combinations of the independent g modes. Not surprisingly anymore, the strongest modes contribute to a large number of these frequencies, as illustrated by the strong bins at 0.059 and 1.901 d-1 in Fig. 19, corresponding to f1f2 and f1, respectively.

thumbnail Fig. 19

Same as Fig. 11, but for KIC 9020774.

Screening through the 239 peaks that remain after filtering the close peaks out, we manage to find another tilted period series (see Table 11 and Fig. 20). The series is steep and smooth, which is in good agreement with the faster rotation rate of the star and its position in the Kiel diagram in Fig. 1.

thumbnail Fig. 20

Same as Fig. 12, but for KIC 9020774.

Table 11

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the period series of KIC 9020774.

3.4.5. KIC 11971405

From the viewpoint of its high projected rotational velocity of 242 ± 14 km s-1 and peculiar, exceptionally weak Be spectra, KIC 11971405 is the most extreme pulsator in our sample. As Kurtz et al. (2015) have already pointed out while studying the occurrence of linear combination frequencies in a sample of gravity mode pulsators, KIC 11971405 shows several outburst recurring on an irregular timescale. The authors also hypothesised the Be nature of the star, which we have now confirmed in Sect. 2.1.4. This is also in good agreement with the high vsini value.

The light curve looks like a typical SPB light curve (but with complex beating patterns and a shorter than typical period) until the middle of Q11, when the first large outburst occurs at TBJD (BJD−2 454 833) 1053. The outburst can be described as a general growth of both the average brightness and amplitude of the light variability. The outburst phase has a faster rising and a longer decaying part; the total length of the first outburst is ~15 d and the maximum brightness is reached ~5 d after the start. The second outburst, which is extremely similar to the first one both in morphology and temporal evolution, starts at TBJD 1129. These are followed by much shorter and shallower mini-outbursts centred at TBJD 1304, 1424, 1517, and 1547 (with a few possible smaller outbursts in between). Then just at the beginning of Q17, starting at TBJD 1559, there is another larger outburst, but since it starts exactly at the beginning of a quarter, it might be slightly influencing our detrending and vice versa.

Similar behaviour has been found in the hybrid B0.5 IVe CoRoT target HD 49330 by Huat et al. (2009). Likely because the CoRoT observation was 10 times shorter in time frame, only one outburst was recorded; but the magnitude of that event was very similar to the brightenings observed here, while its duration was 3 to 4 times longer. Furthermore, although not to the same extent, HD 49330 also exhibited enhanced pulsation amplitudes during the outburst. However, in contrast to the weak Be signatures observed in KIC 11971405, the CoRoT target showed a very strong (and varying) emission component (Floquet et al. 2009).

thumbnail Fig. 21

Short-time Fourier transform of the full light curve of KIC 11971405, showing the regions of dominant pulsation power.

We employ a short-time Fourier transformation (STFT) of the full data set to illustrate the long-term changes in the amplitude spectrum. Afterwards, we apply our normal iterative prewhitening procedure on the Q0-Q11.1 data (until the first monthly downlink within Q11) to derive model frequencies that we can use in our combination and period series search. We do not use the full data set for this procedure, as our initial attempts have shown that it is not possible to adequately model the full light curve with a limited set of frequencies in a meaningful way (which is not surprising when trying to reproduce aperiodic changes with a sum of sines).

The full light curve itself spans 1470.5 days and contains 65 968 data points (91.7% duty cycle and 0.00068 d-1 Rayleigh limit). Because of the complex beating patterns originating in many closely spaced frequencies and the limited duration of the outbursts, it is impossible to use a window width for the STFT that is long enough to avoid visual artefacts, but short enough so that it recovers some of the temporal changes that occur during an outburst. When the temporal resolution is chosen to be too small, such as in the wavelet analysis performed in Rivinius et al. (2016), dense frequency groups remain unresolved. This implies that the observed changes in amplitude are simply a result of the beating of the unresolved modes. Furthermore, because of the merging or smearing effect introduced by the degraded resolution in frequency, the beating among frequencies can lead to artificial changes in amplitudes and/or frequencies, turning a dense group of modes that are otherwise stable in frequency into an unresolved, periodically varying pseudo-mode. In Rivinius et al. (2016, Figs. 6. and 7.), g2, g3, g4, and g5 are suffering from this artefact, which explains the observed short-term correlations between the amplitudes of neighbouring frequency groups and hides the actual long-term changes. Here, we take a careful approach to avoid such artefacts.

From our tests it seems that all dominant frequencies get a boost in amplitude during the outburst events (although the one at 3.73380(1) d-1 to a lesser extent). This explains the visual broadening of the light curve around these events without an apparent change in the pattern of the light variation, which one would expect if the outburst were connected only to a given frequency. What is more interesting is that there is a clear change in the amplitude of two of the most dominant modes on a much longer timescale. The STFT of the three dominant frequency regions of KIC 11971405 are shown in Fig. 21 (calculated using 180 day wide windows that were slid with 30 days). Tracing the changes in amplitude and frequency (Fig. 22), we find that while the peak at 3.73380(1) d-1 continuously grows in time (from around 1000 to 4000 ppm), the peak at 0.27644(2) d-1 loses power (from around 1800 to 900 ppm). The smaller modes are much more difficult to track in the STFT, but is seems that there are some other peaks with smaller, long-term amplitude changes. Moreover, looking at the frequency values, the dominant peaks between 3.6 and 4.2 d-1 (in the 3rd panel of Fig. 21) show a small long-term decreasing trend in frequency, while the dominant modes in the other two frequency regions are stable over time; the low frequency peak at 0.27644(2) d-1 is strongly influenced by the presence of outbursts, which is why there are apparent variations after the first outburst. These changes in frequency are marginally larger than the Rayleigh limit of the full data set, but the trends are very clear and free of scatter, especially for the peak at 4.01030(1) d-1, which is the dominant frequency of the Q0-Q11.1 data (with an apparent amplitude of 3507(56) ppm).

thumbnail Fig. 22

Frequency (left panels) and amplitude (right panels) variations of the four strongest pulsation modes traced from the STFT of KIC 11971405, which was calculated using 360 d wide windows that were stepped with 30 d.

The Q0-Q11.1 data set that we use later on consists of 41 175 data points and covers 911.7 days (92.3% duty cycle and 0.0011 d-1 Rayleigh limit). The iterative prewhitening procedure resulted in 1541 model frequencies; the amplitude weighted average of all model frequencies is 3.14 d-1. While the data set is significantly shorter than those of the four other stars that we have analysed above, the number of frequencies is much higher, which raises some concerns. Although the frequency spectrum seems to be very rich in combination frequencies even beyond the three dominant frequency regions of 0 to 0.7 d-1, 1.3 to 2.3 d-1, and 3.1 to 4.5 d-1, there is no easy way to explain the presence of this many stable modes from a theoretical point of view. Moreover, there is also an excessive amount of peaks in this set of model frequencies that are so closely spaced that they must influence each other during the prewhitening, which is known to create spurious model frequencies. This is why we always filter these close frequencies before searching for period series, but in this star, the filtering removes half the frequencies (leaving 765 of them). Even though a rich and dense frequency spectrum in a fast rotator SPB is not at all impossible given the predicted presence of the numerous rotationally split components of various series of modes of different degrees and radial orders, it is more likely that a large amount of the very closely spaced frequency groups (where a larger central peak is surrounded by a few peaks that have lower amplitudes as we go further from the larger peak) are simply spurious peaks caused by the long-term frequency shifts observed in the STFT of the data. These structures are not to be confused with the narrowly spaced rotationally split triplets or quintuplets, which can have very similar appearance, since those only occur in extremely slow rotator stars.

While this model provides a variance reduction of 99.8% in the Q0-Q11.1 data, it fails miserably in reproducing the outbursts or the light curve after the outbursts. At this point we have no ways to model these features in the light curve, so we proceed with the analysis of the frequency spectrum of the Q0-Q11.1 data.

The frequency spectrum of KIC 11971405 is extremely structured. The histogram of the frequency pairs in Fig. 23 also reflects this (as frequency groups in the periodogram result in groups in the histogram too), but unlike the previous cases, there is not one bin or a few bins that would strongly rise above the underlying general structure. The bin corresponding to the dominant mode is stronger, but not as significant as in the other stars. There is a clear presence of possible combination frequencies, but owing to the complexity of the spectrum, at first sight it is much less clear which peaks could be parent modes and which are the combinations. While in our experience parent peaks are always stronger than the combinations, which is well illustrated by the previous four stars in this paper, it has been shown recently that it is also possible to have combination frequencies that are stronger than their parent peaks (Kurtz et al. 2015). Our standard combination search shows that the independent modes are situated in the 1.3 to 2.3 d-1 region, while as usual, the low frequency peaks are simple negative combinations and the peaks towards higher frequencies are low order positive combinations, even without omitting our standard amplitude criterion (Pápics 2012). There are some large amplitude modes (basically all strong frequencies – including the dominant one – outside the main 1.3 to 2.3 d-1 region), which cannot be explained as combinations without looking over this criterion. But when we only look at frequency values, they are also explained simply, as already shown in Kurtz et al. (2015). Except for these few large amplitude peaks, the amplitudes in the 1.3 to 2.3 d-1 region are in general larger than in the other groups.

thumbnail Fig. 23

Same as Fig. 11, but for KIC 11971405.

There is an additional observation that strengthens our opinion of why the 1.3 to 2.3 d-1 region must contain the independent pulsation modes, and that is the presence of possible period series. We find two possible series at period spacing values which are compatible with the observed high projected rotational velocity and they are both in this region (see Table 12 and Fig. 24). The second series gets slightly more noisy towards higher periods, likely because of the density of the frequency spectrum and the low value of the period spacing, causing the peaks to influence each other in that region. The detection of g-mode period series along with the observed weak Be spectrum confirms that Be stars are simply fast rotator SPB stars.

thumbnail Fig. 24

Same as Fig. 12, but for KIC 11971405.

Table 12

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the two possible period series of KIC 11971405.

There is a possible scenario that could simultaneously explain the formation of unexpectedly large amplitude combination peaks an the observed changes in the amplitude and frequency of these modes (Fig. 22). Take the peak at 4.01030(1) d-1, which is explained by Kurtz et al. as a combination of the peaks 1.82267(2) d-1 and 2.18761(2) d-1. This peak can also be explained as a combination of, at least, two other strong peaks (with a match better than the Rayleigh limit), namely 1.83831(2) d-1 and 2.17192(2) d-1. The difference between the two possible combination peaks is 0.000044 d-1 or 1/25th of the Rayleigh limit. Assuming that both these combination frequencies are actually produced by the four strong parent peaks, we get two peaks that cannot be resolved by 3–4 years of Kepler observations. Thus we detect these as one peak that shows amplitude and frequency variations on the timescale of the beat period of the two unresolved peaks, which is of the order of 60 years. This perfectly explains the excessive amplitude at a given phase of the beat period and the change in amplitude, although the observed change in frequency is larger than what this model predicts. In general, faster rotator SPB stars have denser g-mode period series, thus the formation of similar unresolved combination peaks becomes easier.

For sufficiently fast rotation, the frequency ranges of prograde quadrupole modes and of combination frequencies of dipole modes, may overlap. Given the large number of frequencies, combination frequencies of dipole modes matching observed quadrupole modes can be expected. On the other hand, a series of pure combination frequencies is unlikely to be identifiable as quadrupole modes. In order to avoid an incorrect interpretation, any seismic modelling of such a series has to be carried out with caution. In this case, we actually manage to detect a period spacing pattern (see Table 13 and Fig. 25) among the possible combination frequencies around 3.75 d-1. In Sect. 3.5 we show that this series is consistent with the rotation rate that we derived from the prograde dipole series, assuming (ℓ,m) = (2,2). This indicates that this series could be genuine, and we might be observing a complex mix of beating combination frequencies and actual independent pulsation modes.

Table 13

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the possible series among the high frequency peaks in KIC 11971405.

thumbnail Fig. 25

Same as Fig. 12, but for the possible series among the high frequency peaks in KIC 11971405.

Photometric outbursts:

inspecting the first two large outbursts in the light curve (see Fig. 26, first panel), we find that there is a very noticeable similarity between the two, as already mentioned above. To illustrate this, first we used a cross-correlation function (CCF; see Fig. 26, second panel) to calculate the time shift between the two events, then plotted them on top of each other (to illustrate the similarity; see Fig. 26, third panel) and measured whether there is a good linear correlation between these outbursts (see Fig. 26, last panel). We find that the two outbursts are indeed very similar (with a correlation coefficient R2 = 0.79), and the second follows the first outburst with a dT = 76.81 d delay. The fact that there is such a striking similarity between the topography of the pulsation signal during these outbursts must mean that there is a connection between the outbursts and pulsations. Since the topography is related to the phases of the pulsation modes and the different apparent beatings, the observed strong correlation means that the outbursts happen at the same phase (or constellation of phases) regarding the stellar oscillations. It has been predicted that g modes could indeed be responsible for the outburst events in Be stars, and there is growing observational evidence supporting this idea (see e.g. Rivinius et al. 1998; Shibahashi 2014).

thumbnail Fig. 26

Upper panel: full detrended light curve of KIC 11971405 (in grey) with the regions of the two most prominent outburst indicated in black and red. Second panel: the cross-correlation function of the two outbursts. Third panel: the second outburst (red solid line) shifted back in time with dT = 76.81 d (the value of the highest peak in the CCF function on the second panel) and plotted over the first outburst (black solid line). The match between the two light curves, especially around the maximum amplitude phase, is impressive. Lower panel: a correlation diagram between the observed photometric signal of the two outbursts in the range shown on the third panel. A 1:1 line is plotted in the background as a grey solid line, and the best-fit linear through the scatter plot is plotted using a red solid line (y = 0.81 ± 0.01 × x, R2 = 0.79).

3.5. Period spacing patterns across the SPB instability strip: Initial pattern fitting

To give an overview of the observed period spacing patterns across the SPB instability to date, we show each of the longest main series detected in SPB stars observed by Kepler in Fig. 27. Just by looking at this figure, along with Fig. 10, it is clear that stars that rotate faster (or – to be more precise – have a higher projected rotational velocity), show their strongest excited modes and their main period series towards shorter periods (higher frequencies), and have smaller average spacing values due to a larger tilt.

thumbnail Fig. 27

Main observed period spacing patterns of the five SPB stars studied in this paper along with those of KIC 10526294 and KIC 7760680 from Pápics et al. (2014, 2015). The colours represent the projected rotation velocity of the star; slow rotators are plotted using dark symbols and fast rotators using bright symbols.

The detected period spacing patterns can then be analysed to carry out a mode identification and derive the approximate near-core rotation rates by fitting model patterns to these series. Because the interior chemical gradients do not affect the global slope of the observed period spacing patterns, we can ignore their effect for now and start from an asymptotic period spacing pattern, as described by Tassoul (1980). Realistic values for the asymptotic period spacing of SPB stars range from around 4000 s to 11 000 s. The influence of rotation is then introduced using the traditional approximation, where we assumed that, first, the star is spherically symmetric and, second, is rigidly rotating. To this end, we applied the traditional approximation module from the stellar pulsation code GYRE v4.3 (Townsend & Teitler 2013). We then fit the resulting model patterns to the observed series. The used analysis methods are described in more detail in Van Reeth et al. (2016).

thumbnail Fig. 28

Dominant pulsation frequency in the inertial frame as a function of the derived rotation frequency frot for the detected period series (along with the previously known series of KIC 10526294 and KIC 7760680 from Pápics et al. 2014, 2015). The thick vertical lines indicate the full extent of the observed period series. The dashed black line indicates where finert equals to frot.

thumbnail Fig. 29

Same as Fig. 28, but in the corotating frame. The dashed red line shows where the pulsations pass from the superinertial regime (above the line) into the subinertial regime (below the line).

The dominant period spacing pattern (in Fig. 27) of each of the stars was found to be a series of (ℓ,m) = (1, + 1) prograde dipole modes. This is consistent with the results obtained by Van Reeth et al. (2016) for γ Dor stars. Other derived mode identifications include (ℓ,m) = (2,−2) for one of the rising series in KIC 4930889, and (2, + 2) for the extra series in KIC 3459297 and for the possible series in the range of combination frequencies in KIC 11971405. However, contrary to the study of the lower mass γ Dor stars, none of the observed series were identified as Rossby waves. A possible explanation for this non-detection may be related to the absence of a thin convective envelope in SPBs or different properties in the interior rotation between these two groups of gravity-mode pulsators; a proper interpretation of this different behaviour requires future theoretical research. The results of this exercise are illustrated in Figs. 28 and 29, while the derived rotation frequencies and asymptotic period spacings are listed in Table, 14. Finally, the derived asymptotic spacing for KIC 4930889 is more compatible with the fundamental parameters of the primary component of this system than of the secondary.

4. Conclusions

By presenting four new SPBs and revisiting one of the previously published binary components, we completed the observational analyses of the stars in our Kepler Guest Observer Programme. We raised the number of in-depth analysed SPB stars in the initial Kepler field from four to eight, and the number of SPB stars with an observed period spacing from four to nine; seven of these, including the five new SPB stars in this paper, are from Kepler. Simultaneously we have shown that period series are not only common, but they – and their combination peaks – dominate the frequency spectrum of SPB stars. We can find a wealth of information about the internal physics of these objects imprinted in the topology of these series, thus they serve as foundations on which we can build an improved theory.

While this investigation along with previous studies (Pápics et al. 2013, 2014, 2015) have already provided previously unseen constraints for theoretical modelling efforts (Moravveji et al. 2015; Triana et al. 2015; Moravveji et al. 2016), the exploitation of Kepler data needs to continue. There must be still a few dozen SPB stars in the Kepler data archives. These light curves are ten times better than anything we had pre-Kepler, and there is no upcoming mission the coming decade that would deliver data with the same or superior properties given that PLATO (Rauer et al. 2014) will only get launched in 2025. Although TESS (Transiting Exoplanet Survey Satellite; Ricker et al. 2014, launch in 2017) will be an all-sky survey, it will at best provide only one-fourth of the frequency resolution of Kepler, and only in the very limited continuous viewing zones around the ecliptic poles. Therefore in the meantime we must continue with screening the Kepler data for similar period series in SPB stars covering various regions of the instability strip to contribute to a statistically meaningful sample of these intriguing objects via in-depth modelling efforts based on their seismic barcodes.

Table 14

Results from the initial pattern fitting of the observed period series.


Acknowledgments

The research leading to these results has received funding from The Research Foundation – Flanders (FWO), Belgium, from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 670519: MAMSIE), and from the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme FP7/2007–2013 under REA grant agreement No. 623303 (ASAMBA). T.V.R. acknowledges financial support from the Fund for Scientific Research of Flanders (FWO), Belgium, under grant agreement G.0B69.13. We are grateful to Bill Paxton and Richard Townsend for their valuable work on the stellar evolution code MESA and stellar pulsation code GYRE. E.N. acknowledges support from the NCN grant No. 2014/13/B/ST9/00902. Funding for the Kepler Mission is provided by the NASA Science mission directorate. Some of the data presented in this paper were obtained from the Multimission Archive at the Space Telescope Science Institute (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 NNX09AF08G and by other grants and contracts. This work is partly based on observations made with the William Herschel Telescope, which is operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos (ORM) of the Instituto de Astrofísica de Canarias (IAC). We made use of the matplotlib colormaps by Nathaniel J. Smith, Stefan van der Walt, and Eric Firing.

References

  1. Aerts, C. 2013, in EAS Pub. Ser., 64, 323 [Google Scholar]
  2. Aerts, C. 2015, in IAU Symp., 307, 154 [Google Scholar]
  3. Aerts, C., Thoul, A., Daszyńska, J., et al. 2003, Science, 300, 1926 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  4. Aerts, C., Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology (Heidelberg: Springer Berlin) [Google Scholar]
  5. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  6. Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Balona, L. A., Pigulski, A., De Cat, P., et al. 2011, MNRAS, 413, 2403 [NASA ADS] [CrossRef] [Google Scholar]
  8. Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  9. Bouabid, M.-P., Dupret, M.-A., Salmon, S., et al. 2013, MNRAS, 429, 2500 [NASA ADS] [CrossRef] [Google Scholar]
  10. Briquet, M., Morel, T., Thoul, A., et al. 2007, MNRAS, 381, 1482 [NASA ADS] [CrossRef] [Google Scholar]
  11. Castro, N., Fossati, L., Langer, N., et al. 2014, A&A, 570, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. De Cat, P., & Aerts, C. 2002, A&A, 393, 965 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. De Cat, P., Briquet, M., Daszyńska-Daszkiewicz, J., et al. 2005, A&A, 432, 1013 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Debosscher, J., Blomme, J., Aerts, C., & De Ridder, J. 2011, A&A, 529, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Degroote, P., Aerts, C., Baglin, A., et al. 2010, Nature, 464, 259 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  16. Degroote, P., Aerts, C., Michel, E., et al. 2012, A&A, 542, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Dziembowski, W. A., Moskalik, P., & Pamyatnykh, A. A. 1993, MNRAS, 265, 588 [NASA ADS] [Google Scholar]
  18. Floquet, M., Hubert, A., Huat, A., et al. 2009, A&A, 506, 103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Gautschy, A., & Saio, H. 1993, MNRAS, 262, 213 [NASA ADS] [CrossRef] [Google Scholar]
  20. Gilliland, R. L., Brown, T. M., Christensen-Dalsgaard, J., et al. 2010, PASP, 122, 131 [NASA ADS] [CrossRef] [Google Scholar]
  21. Hadrava, P. 1995, A&AS, 114, 393 [NASA ADS] [Google Scholar]
  22. Huat, A., Hubert, A., Baudin, F., et al. 2009, A&A, 506, 95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Ilijic, S., Hensberge, H., Pavlovski, K., & Freyhammer, L. M. 2004, in Spectroscopically and Spatially Resolving the Components of the Close Binary Stars, eds. R. W. Hilditch, H. Hensberge, & K. Pavlovski, ASP Conf. Ser., 318, 111 [Google Scholar]
  24. Kepler Mission Team 2009, VizieR Online Data Catalog: V/133 [Google Scholar]
  25. Kupka, F. G., Ryabchikova, T. A., Piskunov, N. E., Stempels, H. C., & Weiss, W. W. 2000, Balt. Astron., 9, 590 [NASA ADS] [Google Scholar]
  26. Kurtz, D. W., Saio, H., Takata, M., et al. 2014, MNRAS, 444, 102 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kurtz, D. W., Shibahashi, H., Murphy, S. J., Bedding, T. R., & Bowman, D. M. 2015, MNRAS, 450, 3015 [NASA ADS] [CrossRef] [Google Scholar]
  28. Loumos, G. L., & Deeming, T. J. 1978, Ap&SS, 56, 285 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Marsh, T. R. 1989, PASP, 101, 1032 [NASA ADS] [CrossRef] [Google Scholar]
  30. McNamara, B. J., Jackiewicz, J., & McKeever, J. 2012, AJ, 143, 101 [NASA ADS] [CrossRef] [Google Scholar]
  31. Miglio, A., Montalbán, J., Noels, A., & Eggenberger, P. 2008, MNRAS, 386, 1487 [NASA ADS] [CrossRef] [Google Scholar]
  32. Moravveji, E. 2016, MNRAS, 455, L67 [NASA ADS] [Google Scholar]
  33. Moravveji, E., Aerts, C., Pápics, P. I., Triana, S. A., & Vandoren, B. 2015, A&A, 580, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Moravveji, E., Townsend, R. H. D., Aerts, C., & Mathis, S. 2016, ApJ, 823, 130 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  35. Pamyatnykh, A. A., Handler, G., & Dziembowski, W. A. 2004, MNRAS, 350, 1022 [NASA ADS] [CrossRef] [Google Scholar]
  36. Pápics, P. I. 2012, Astron. Nach., 333, 1053 [Google Scholar]
  37. Pápics, P. I. 2013, Ph.D. Thesis, KU Leuven, Belgium [Google Scholar]
  38. Pápics, P. I., Briquet, M., Baglin, A., et al. 2012, A&A, 542, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Pápics, P. I., Tkachenko, A., Aerts, C., et al. 2013, A&A, 553, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Pápics, P. I., Moravveji, E., Aerts, C., et al. 2014, A&A, 570, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Pápics, P. I., Tkachenko, A., Aerts, C., et al. 2015, ApJ, 803, L25 [Google Scholar]
  42. Raskin, G., van Winckel, H., Hensberge, H., et al. 2011, A&A, 526, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Rauer, H., Catala, C., Aerts, C., et al. 2014, Exp. Astron., 38, 249 [NASA ADS] [CrossRef] [Google Scholar]
  44. Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, in Space Telescopes and Instrumentation 2014: Optical, Infrared, and Millimeter Wave, Proc. SPIE, 9143, 914320 [Google Scholar]
  45. Rivinius, T., Baade, D., Stefl, S., et al. 1998, in A Half Century of Stellar Pulsation Interpretation, eds. P. A. Bradley, & J. A. Guzik, ASP Conf. Ser., 135, 343 [Google Scholar]
  46. Rivinius, T., Carciofi, A. C., & Martayan, C. 2013, A&ARv, 21, 69 [Google Scholar]
  47. Rivinius, T., Baade, D., & Carciofi, A. C. 2016, A&A, 593, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Rogers, T. M. 2015, ApJ, 815, L30 [NASA ADS] [CrossRef] [Google Scholar]
  49. Schmidt-Kaler, T. 1982, Landolt-Börnstein: Numerical Data and Functional Relationships in Science and Technology, Vol. 2b, eds. K. Schaifers, & H. H. Vogt (Springer-Verlag) [Google Scholar]
  50. Shibahashi, H. 2014, in Precision Asteroseismology, eds. J. A. Guzik, W. J. Chaplin, G. Handler, & A. Pigulski, IAU Symp., 301, 173 [Google Scholar]
  51. Shulyak, D., Tsymbal, V., Ryabchikova, T., Stütz, C., & Weiss, W. W. 2004, A&A, 428, 993 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Simon, K. P., & Sturm, E. 1994, A&A, 281, 286 [NASA ADS] [Google Scholar]
  53. Tassoul, M. 1980, ApJS, 43, 469 [NASA ADS] [CrossRef] [Google Scholar]
  54. Tkachenko, A. 2015, A&A, 581, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Townsend, R. H. D., & Teitler, S. A. 2013, MNRAS, 435, 3406 [NASA ADS] [CrossRef] [Google Scholar]
  56. Triana, S. A., Moravveji, E., Pápics, P. I., et al. 2015, ApJ, 810, 16 [NASA ADS] [CrossRef] [Google Scholar]
  57. Tsymbal, V. 1996, in M.A.S.S., Model Atmospheres and Spectrum Synthesis, eds. S. J. Adelman, F. Kupka, & W. W. Weiss, ASP Conf. Ser. 108, 198 [Google Scholar]
  58. Štefl, S., Baade, D., Rivinius, T., et al. 2003, A&A, 402, 253 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Van Reeth, T., Tkachenko, A., & Aerts, C. 2016, A&A, 593, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Waelkens, C. 1991, A&A, 246, 453 [NASA ADS] [Google Scholar]
  61. Walker, G., Matthews, J., Kuschnig, R., et al. 2003, PASP, 115, 1023 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Logbook of the spectroscopic observations – grouped by stars and observing runs.

Table 2

Fundamental parameters and basic observable properties of KIC 3459297.

Table 3

Fundamental parameters, orbital parameters, and basic observable properties of KIC 4930889.

Table 4

Fundamental parameters and basic observable properties of KIC 9020774.

Table 5

Fundamental parameters and basic observable properties of KIC 11971405.

Table 6

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the main period series of KIC 3459297.

Table 7

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the possible l = 2 series in KIC 3459297.

Table 8

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the main period series of KIC 4930889.

Table 9

Fourier parameters (periods, frequencies, amplitudes, and phases) of the two extra possible period series of KIC 4930889.

Table 10

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the period series detected in KIC 6352430 A.

Table 11

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the period series of KIC 9020774.

Table 12

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the two possible period series of KIC 11971405.

Table 13

Fourier parameters (periods, frequencies, amplitudes, and phases) of the modes of the possible series among the high frequency peaks in KIC 11971405.

Table 14

Results from the initial pattern fitting of the observed period series.

All Figures

thumbnail Fig. 1

Kiel diagram of a sample of B-type stars near the main sequence, for which an in-depth seismic data analysis is available, observed by CoRoT (circles) and Kepler (squares). The new SPBs from this paper are plotted with diamonds. Stars analysed or revisited in this paper are labelled using their full names, while for the others we use the following annotations: C1 – HD 48977, C2 – HD 43317, C3 – HD 50230 A, C4 – HD 50846 A, C5 – HD 182198, C6 – HD 181440, C7 – HD 46179, C8 – HD 174648, C9 – HD 174884 AB, K1 – KIC 4931738 AB, K2 – KIC 10526294, and K3 – KIC 7760680. The dot-dashed line represents the zero-age main sequence (ZAMS), while the dashed line indicates the terminal-age main sequence (TAMS). The cool edge of the β Cep instability strip is plotted with a dotted line (for the = 0 modes). The thin grey solid lines denote evolutionary tracks for selected masses. Isochrones from 20 to 28 Myr and isoradii corresponding to integer multiples of the Solar radius are plotted using thin dashed and dotted grey lines, respectively. The instability strip of SPB stars is defined as the region where the sum of excited = 1 plus = 2 gravity modes is at least three, and it is shown in colours ranging from purple (3 excited modes) to bright yellow (63 excited modes). All model data were taken from Moravveji (2016) for Z = 0.014, the Asplund et al. (2009) mixture, an exponential core overshoot of fov = 0.02, and opacity enhancement factors for iron and nickel of βFe = 1.75 and βNi = 1.75. Different error bars reflect differences in data quality and methodological approach.

In the text
thumbnail Fig. 2

Comparison of the rectified observed spectrum and the synthetic spectrum for the two observed wavelength regions of KIC 3459297. In each panel, the Isis spectrum is plotted with a grey solid line and the synthetic spectrum is plotted with a red solid line.

In the text
thumbnail Fig. 3

Phase-folded residual spectra of KIC 4930889 around the H i line at 4861 Å obtained after subtraction of the average. The colour scale represents the residual intensity at each wavelength point. White horizontal regions denote phases that are more than 0.05π away from an observed spectrum.

In the text
thumbnail Fig. 4

Spectroscopic orbital solution for KIC 4930889 (solid line: primary; dashed line: secondary) and the radial velocities from the disentangling process. Filled circles indicate RVs of the primary; filled triangles denote RVs of the secondary component.

In the text
thumbnail Fig. 5

Comparison of the Hermes spectrum at the phase of the largest observed radial velocity difference, disentangled component spectra, and synthetic spectra for each component for selected wavelength regions of KIC 4930889. In each panel, the Hermes spectrum is plotted with a black solid line, the disentangled spectra are plotted with grey solid lines, and synthetic spectra (scaled to the disentangled spectra) are plotted with red solid lines. For clarity, the primary and secondary spectra are plotted with a vertical offset of + 0.1 and + 0.2 flux units, respectively.

In the text
thumbnail Fig. 6

Comparison of the synthetic spectrum and the average of the rectified observed spectra for the two observed wavelength regions of KIC 9020774. In each panel, the Isis spectrum is plotted with a grey solid line, and the synthetic spectrum is plotted with a red solid line.

In the text
thumbnail Fig. 7

Comparison of the median of the rectified observed spectra, and the synthetic spectrum for selected wavelength regions of KIC 11971405. In each panel, the Hermes spectrum is plotted with a grey solid line, and the synthetic spectrum is plotted with a red solid line.

In the text
thumbnail Fig. 8

Temporal changes in the spectrum of KIC 11971405. The average of the two consecutive exposures from October 2010 is plotted on the top using a black solid line, while the average spectrum from June 2015 is shown below with a red solid line with an offset of −0.15 in flux units. For clarity, the average from 2010 is repeated in the background using a grey solid line. The sharp absorption features in the bottom panel are telluric lines.

In the text
thumbnail Fig. 9

Custom pixel mask for the Q8 data of KIC 4930889. Individual pixel-level light curves are plotted for each pixel that was downloaded from the spacecraft. The numerical values and the background colour indicate the flux levels in each pixel. The pixels with green and red borders were used to extract the standard Kepler light curves. We added yellow pixels (significant signal is present) and removed red pixels (non-significant signal or contamination from nearby stars) in our custom mask for the light curve extraction. This results in a light curve with significantly less instrumental effects and/or contamination than the standard extraction.

In the text
thumbnail Fig. 10

Light curves (plotted using grey dots, showing a zoom-in to slightly less than 1/10th of the full data set) and the power spectra (plotted using red solid lines) of the five SPB stars studied in this paper along with the ones of KIC 10526294 and KIC 7760680 from Pápics et al. (2014, 2015) ordered according to the frequency of maximum power. For further information see Sect. 3.5.

In the text
thumbnail Fig. 11

Histogram of the distribution of frequency pairs in the frequency analysis of KIC 3459297, calculated using a bin width of 2.5 /T taking into account all statistically significant frequencies (grey solid line), or only a subset of frequencies where weaker signal from within 2.5 /T of the stronger peaks was excluded (black solid line). The three values that are mentioned in the text are marked with black triangles.

In the text
thumbnail Fig. 12

Zoom-in of the period spectrum of KIC 3459297 (top panel) with the members of the period spacing series indicated with vertical dashed lines. Peaks that are discarded during our pattern search are plotted using thin grey lines (for details, we refer to the text). The observed period spacing is shown in the bottom panel.

In the text
thumbnail Fig. 13

Same as Fig. 12, but for the possible l = 2 series in KIC 3459297.

In the text
thumbnail Fig. 14

Same as Fig. 11, but for KIC 4930889.

In the text
thumbnail Fig. 15

Same as Fig. 12, but for KIC 4930889.

In the text
thumbnail Fig. 16

Same as Fig. 12, but for the proposed additional series of the low frequency peaks in KIC 4930889.

In the text
thumbnail Fig. 17

Same as Fig. 11, but for KIC 6352430 A.

In the text
thumbnail Fig. 18

Same as Fig. 12, but for KIC 6352430 A.

In the text
thumbnail Fig. 19

Same as Fig. 11, but for KIC 9020774.

In the text
thumbnail Fig. 20

Same as Fig. 12, but for KIC 9020774.

In the text
thumbnail Fig. 21

Short-time Fourier transform of the full light curve of KIC 11971405, showing the regions of dominant pulsation power.

In the text
thumbnail Fig. 22

Frequency (left panels) and amplitude (right panels) variations of the four strongest pulsation modes traced from the STFT of KIC 11971405, which was calculated using 360 d wide windows that were stepped with 30 d.

In the text
thumbnail Fig. 23

Same as Fig. 11, but for KIC 11971405.

In the text
thumbnail Fig. 24

Same as Fig. 12, but for KIC 11971405.

In the text
thumbnail Fig. 25

Same as Fig. 12, but for the possible series among the high frequency peaks in KIC 11971405.

In the text
thumbnail Fig. 26

Upper panel: full detrended light curve of KIC 11971405 (in grey) with the regions of the two most prominent outburst indicated in black and red. Second panel: the cross-correlation function of the two outbursts. Third panel: the second outburst (red solid line) shifted back in time with dT = 76.81 d (the value of the highest peak in the CCF function on the second panel) and plotted over the first outburst (black solid line). The match between the two light curves, especially around the maximum amplitude phase, is impressive. Lower panel: a correlation diagram between the observed photometric signal of the two outbursts in the range shown on the third panel. A 1:1 line is plotted in the background as a grey solid line, and the best-fit linear through the scatter plot is plotted using a red solid line (y = 0.81 ± 0.01 × x, R2 = 0.79).

In the text
thumbnail Fig. 27

Main observed period spacing patterns of the five SPB stars studied in this paper along with those of KIC 10526294 and KIC 7760680 from Pápics et al. (2014, 2015). The colours represent the projected rotation velocity of the star; slow rotators are plotted using dark symbols and fast rotators using bright symbols.

In the text
thumbnail Fig. 28

Dominant pulsation frequency in the inertial frame as a function of the derived rotation frequency frot for the detected period series (along with the previously known series of KIC 10526294 and KIC 7760680 from Pápics et al. 2014, 2015). The thick vertical lines indicate the full extent of the observed period series. The dashed black line indicates where finert equals to frot.

In the text
thumbnail Fig. 29

Same as Fig. 28, but in the corotating frame. The dashed red line shows where the pulsations pass from the superinertial regime (above the line) into the subinertial regime (below the line).

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.