Asteroseismic modelling of the two F-type hybrid pulsators KIC 10080943A and KIC 10080943B

Asteroseismic modelling of the two F-type hybrid pulsators KIC 10080943A and KIC 10080943B

V. S. Schmid Aspirant PhD Fellow of the Research Foundation Flanders (FWO), Belgium1 Institute for Astronomy, KU Leuven, Celestijnenlaan 200D, B – 3001 Leuven, Belgium
   C. Aerts 1 Institute for Astronomy, KU Leuven, Celestijnenlaan 200D, B – 3001 Leuven, Belgium 2Department of Astrophysics/IMAPP, Radboud University Nijmegen, P.O. Box 9010, 6500 GL Nijmegen, The Netherlands 2
Received; accepted
Key Words.:
stars: individual: KIC 10080943 - binaries: spectroscopic - stars: oscillations - stars: evolution

Context:Pulsating binary stars are ideal targets for testing the theory of stellar structure and evolution. Fundamental parameters can be derived to high precision from binary modelling and provide crucial constraints for seismic modelling. High-order gravity modes are sensitive to the conditions near the convective core and therefore allow for a determination of parameters describing interior physics, especially the convective-core overshooting parameter. KIC 10080943 is a binary system that contains two gravity- and pressure-mode hybrid pulsators. A detailed observational study has provided fundamental and seismic parameters for both components.

Aims:We aim to find a model that is able to predict the observed g-mode period spacings and stellar parameters of both components of KIC 10080943.

Methods:By calculating model grids with the stellar evolution code MESA and the seismic code GYRE, we can compare theoretical properties to the observed mean period spacing and position in the Hertzsprung-Russell diagram.

Results:The masses of our best models are somewhat below the values estimated from binarity, which is a consequence of the low observed mean g-mode period spacing. We find that the amount of core overshooting and diffusive mixing can be well constrained by the equal-age requirement for the two stars, however, we find no significant difference for different shapes of core overshooting. The measured rotation rates are within the limit of validity for the first-order perturbation approximation. We can find a good fit by using the traditional approximation for the pulsations, when taking slightly younger models with a higher asymptotic period spacing. This is because the zonal modes experience a slight shift due to the Coriolis force, which the first-order perturbation approximation ignores.


1 Introduction

Stellar models play an important role in aiding astrophysical research in a wide variety of contexts, despite several shortcomings and uncertainties in the input physics. Fortunately, we can test and improve these models by comparing observed stellar pulsations, that is the asteroseismic fingerprints of the stellar interior, to theoretical predictions. It is imperative to choose the test cases wisely by the amount of observational constraints they hold. As such, gravity (g) and pressure (p) mode hybrid pulsators have a high potential for testing theories of stellar structure and evolution.

The periods of high-order g modes of a non-rotating, unevolved star with a convective core and a radiative envelope are predicted to be equally spaced (Tassoul, 1980). Yet, deviations from this equidistant spacing are expected as the star evolves or rotates (Miglio et al., 2008; Bouabid et al., 2013). The extent of the convective core and chemical composition gradients at the boundary of that core influence the mean value and structure of the period spacing, while rotation introduces a tilt in the period spacing pattern. This is of particular interest for stars of intermediate mass () that inhabit the transition region, where a radiative central region gives way to a convective core and the convective envelope becomes increasingly thinner. For these types of stars the g modes are thought to be excited by the convective flux blocking mechanism (Guzik et al., 2000; Dupret et al., 2005). The heat-driven p modes of hybrid pulsators can then provide an additional selection criterion based on the predictions of mode excitation. Very recently, there have been several studies exploiting the seismic potential of g-mode pulsators with carefully identified modes from four years of Kepler (Borucki et al., 2010) data (e.g. Kurtz et al., 2014; Moravveji et al., 2015; Pápics et al., 2015; Saio et al., 2015; Van Reeth et al., 2015a, b; Moravveji et al., 2016).

The main caveat of these studies is that they only concern single stars. Their derived fundamental parameters and ages are therefore still dependent on the accuracy of stellar models on which the seismic properties are based. An additional and independent test can be provided by pulsating binary stars, whose measurements of fundamental parameters rely only on Kepler’s laws and geometry. Double-lined spectroscopic binaries, where radial velocities of both components are measured, yield a mass ratio, while eclipsing binaries allow for the determination of absolute masses and radii. Such systems are rare and require a dedicated, detailed study to recover the seismic properties of the pulsating star and precise stellar parameters. There have only been a handful of published analyses so far (e.g. Maceroni et al., 2009, 2014; Welsh et al., 2011; Hambleton et al., 2013).

In this paper we present the seismic modelling of two hybrid pulsators; KIC 10080943A and KIC 10080943B. They reside in a double-lined spectroscopic binary system and have very similar seismic properties and stellar parameters, which are described in Sect. 2 and whose analysis has been published by Schmid et al. (2015). We use the one-dimensional (1D) stellar evolution code MESA (Paxton et al., 2011, 2013, 2015) and the stellar pulsation code GYRE (Townsend & Teitler, 2013) to find the best representation of both stars, assuming two coeval single stars and neglecting their tidal influence. The modelling process is explained in detail in Sect. 3. Finally, we present a discussion of our results in Sect. 4.

2 Observational constraints on KIC 10080943

Figure 1: Observed g modes in frequency (top) and period (bottom) as a function of degree are shown for the primary in red and the secondary in blue. The values are distinguished by the different symbols. The yellow stars denote the modes, which could not be assigned to either component.
Figure 2: Observed p-mode frequencies as a function of degree , which was estimated, assuming that all observed multiplets are complete. Modes of the primary are shown in red and those of the secondary are shown in blue. The values are distinguished by the different symbols.

KIC 10080943 is an eccentric, double-lined binary system, containing two  Dor/ Sct hybrid pulsators (for a detailed description of  Dor and  Sct type pulsators, see Aerts et al., 2010, chapt. 2). Keen et al. (2015) published the analysis of the period spacings of the g modes from Kepler photometry. Schmid et al. (2015) performed a detailed analysis of the system; these authors derived fundamental parameters of both components from modelling the binary signals of ellipsoidal variation and reflection in the Kepler light curve in addition to the radial velocity curves and disentangled component spectra from high-resolution follow-up spectroscopy covering the orbit. They were able to assign most period spacing series and high-amplitude p modes to either the primary or the secondary. This was achieved by analysing the phase modulations of the p modes that arise from the binary motion. Modes of the primary have a time delay of  s, while those of the secondary are delayed by  s. Combination frequencies between the p modes and two period spacing series, which form a series of rotationally split doublets, helped to determine that these belong to the secondary. Three other period spacing series form series of triplets with a different rotational splitting value and thus originate in the primary. Figure 1 shows the six period spacing series and which star they originate in. The series of modes could not be assigned to either star, as it does not show combination frequencies or rotational splitting. In Fig. 2 the p modes that can be connected to one star or the other are shown.

All derived stellar parameters including the mean period spacing values are listed in Table 1 and are used to place both stars in a Hertzsprung-Russell diagram (HRD) inside a 1- error box in Fig. 3. It is obvious that the parameters lack the precision typical for double-lined eclipsing binaries (Southworth, 2012), as no eclipses can be observed for KIC 10080943. The radii are then derived from the equipotential surfaces in Roche geometry, as opposed to the width of the eclipses. The ratio of the effective temperatures is somewhat constrained by both effects, while they are largely independent of the orbital inclination (for more information on the binary and a description of the modelling process, see Schmid et al., 2015). This high uncertainty of the orbital inclination propagates into the precision of the absolute masses. The mass ratio of the two components, however, was derived with great confidence from the ratio of the radial velocity amplitudes. Hence, the binarity of the system still yields crucial additional constraints for the theoretical stellar models compared to single star pulsators, as is discussed below.

Primary Secondary
-0.050.17 -0.090.30
Table 1: Fundamental and seismic parameters of the primary and secondary component of KIC 10080943.

The main diagnostic for finding a model representative of KIC 10080943A and KIC 10080943B is the period spacing of the high-order g modes. For both components, the mean period spacing value is rather low at  seconds for stars near . Indeed, comparing with the sample of  Dor stars in Van Reeth et al. (2015b), the measured values of both components hint at masses . Furthermore, the observed period spacing series show a significant wavy structure without clear dips and a slight tilt, which is caused by rotation. This points towards evolved stars with a low central hydrogen mass fraction, where the receding convective hydrogen-burning core has left behind a chemical composition gradient leading to a peak in the Brunt-Väisälä frequency and therefore a decrease in the asymptotic period spacing (e.g. Miglio et al., 2008).

In addition to the g modes, two p modes were identified as candidate radial modes and can help distinguish between possible models. To estimate their radial overtone, we calculate the pulsation constant (Breger, 1990) as


or, expressed in observable parameters


For the highest amplitude p mode () and for the parameters of the secondary, we find , which is indicative of the radial fundamental mode. The p mode with frequency belongs to the primary and has , which is consistent with (Stellingwerf, 1979).

In addition to these seismic diagnostics, the main constraint from the binarity is the requirement of equal age, equal composition, and the mass ratio. This means that the seismic properties of both models have to match the observations at the same age, while the ratio of their masses also has to be in agreement with the observational measurement.

3 Searching for the best seismic models

We used MESA (Paxton et al., 2011, 2013, 2015, version 7385) to calculate evolutionary tracks for different parameters, which are specified in the following subsection. For all tracks we used OPAL opacity tables (Iglesias & Rogers, 1993, 1996) constructed for the solar metal mixture provided by Asplund et al. (2009). Stellar layers are dynamically unstable against convection according to the Schwartzschild criterion and are modelled with the mixing length theory (MLT) by Cox & Giuli (1968, chap. 14), using . All equilibrium models along the evolutionary tracks are non-rotating, as both components are slow rotators (rotational period near the core d and d, Keen et al., 2015). Our MESA inlist is provided in Appendix A and is available for further choices of input physics.

As already mentioned in Sect. 2, the Brunt-Väisälä frequency increases with the steepness of the chemical composition gradient . The shape and strength of the mixing of the stellar material therefore influence the shape of and hence the periods and period spacings of high-order g modes. This can be illustrated by the asymptotic period spacing, which is defined by Tassoul (1980) as


where , and is the convective-core boundary. To probe the strength of mixing beyond the convective regions, we included convective-core overshooting and diffusive mixing as free parameters into our models. Convective-core overshooting mainly increases the size of the convective core, thus decreasing the integral . It is implemented in MESA using two different definitions that both depend on the MLT diffusion coefficient () and an adjustable parameter. In the exponential description, the region immediately outside the core has a diffusion coefficient as determined by


where is the local pressure scale height, is the distance from the convective-core boundary, and is a free parameter (Herwig, 2000). This definition results in an exponentially decaying . Alternatively overshooting can also be defined as a step function, where the extent of the overshooting layer is with the adjustable parameter . Typically and increasing either parameter leads to a wider mixed region beyond the core. In the radiative parts of the model, where neither convection nor overshooting is active, we set a minimum diffusion coefficient of . It is implemented such that the diffusion coefficient throughout the star does not fall below this value. The higher it is, the more the chemical composition gradient is washed out in the radial direction, leading to fewer dips in the period spacing pattern. Apart from affecting , both overshooting and diffusion also have the consequence of mixing more hydrogen into the central region and therefore extending the time the stars spend on the main sequence. In other words, otherwise equal stars have a higher at a given age when only the mixing strength is increased.

begin end step
coarse grid 704
() 1.7 2.2 0.05 11
0.012 0.018 0.002 4
0.0 0.03 0.01 4
0 3 1 4
fine grid (exponential overshooting) 3267
() 1.76 1.86 0.01 11
0.0105 0.0125 0.001 3
0.005 0.015 0.001 11
0 2 0.25 9
fine grid (step overshooting) 3267
() 1.76 1.86 0.01 11
0.0105 0.0125 0.001 3
0.05 0.15 0.01 11
0 2 0.25 9
111The overshooting parameters and are expressed in terms of local pressure scale height () and the diffusion coefficient has the unit cm s.
Table 2: Extents of the grid parameters.

Based on these stellar models, we calculate theoretical pulsation modes, using the stellar pulsation code GYRE (Townsend & Teitler, 2013, version 4.3). To take the effect rotation has on the observed frequencies into account, we perturbed the frequencies of the calculated zonal modes. In a first step we apply the first-order perturbative approach (Ledoux, 1951), where and is the Ledoux constant. We use as calculated by GYRE and the observed frequency splitting as a first guess for the rotational frequency . The validity of the perturbative approach for  Dor stars was tested by Ballot et al. (2010) for different rotation rates. For both stars and the observational estimates of the rotation rates (Schmid et al., 2015) , where is the rotation rate near the core and are the g-mode frequencies in the co-rotating frame of reference. This means that they are in the validity range of the perturbative approach to first order.

An alternative method to treat the coupling between rotation and pulsation is the traditional approximation of rotation (TAR; e.g. Townsend, 2003), which is also implemented in GYRE. This method is based on neglecting the horizontal component of the rotation vector and is valid when the horizontal displacement is larger than the radial displacement, which is the case for g modes in  Dor stars. Bouabid et al. (2013) studied the influence of the Coriolis force on high-order g modes, using the TAR, and compared their computations to predictions by the first-order perturbative approach. These authors showed that, even for slow rotation rates, the Coriolis force has an influence on the mode periods, which is stronger at higher radial orders. This force causes a departure from linearity in rotational splitting and also effects the zonal modes, which are unperturbed in the perturbative approach.

3.1 MESA model grids

Figure 3: Top panel: HRD with the observed positions and areas of both components of KIC 10080943, marked by the boxes (the primary is the cooler, more luminous star). The positions of the components in Model 2 are given by the diamonds and for Model 3 by the squares, while the respective evolutionary tracks are given by the red and black solid lines. Dark blue symbols mark the primary, while the light blue symbols indicate the secondary models. Parameters of these models are given in Table 3. Lower panel: Asymptotic period spacing as a function of central hydrogen mass fraction . The observed mean period spacing of the primary (secondary) is given as the horizontal solid (dashed) line and the areas are shown as the grey patch.
Model 1 Model 2 Model 3 Model 4 Model 5 Model 6
grid coarse exp. ov. (Ledoux) step ov. (Ledoux) exp. ov. (TAR) step ov. (TAR) morphology
() 1.85 1.82 1.81 1.82 1.81 1.67
() 1.75 1.76 1.76 1.76 1.76 1.60
0.012 0.0125 0.0125 0.0125 0.0125 0.010
or 0.01 0.008 0.11 0.008 0.11 0.007
or 0.0 0.005 0.05 0.005 0.05 0.006
0 0.25 0.25 0.25 0.25 0.5
2 1.5 1.75 1.5 1.75 0.75
(Gyr) 1.087 1.123 1.110 1.102 1.089 1.082
(Gyr) 1.082 1.127 1.110 1.100 1.081 1.085
0.09 0.1 0.11 0.12 0.13 0.22
0.24 0.19 0.21 0.22 0.23 0.29
(s) 2832 2806 2816 2873 2877 2829
(s) 2909 2905 2904 2939 2942 2851
(K) 7131 7074 7110 7157 7189 7583
(K) 7526 7308 7372 7387 7452 7597
() 2.92 2.81 2.72 2.74 2.66 2.12
() 2.29 2.43 2.39 2.38 2.33 1.90
1.30 1.25 1.23 1.25 1.23 1.13
1.18 1.18 1.18 1.18 1.18 1.04
222The overshooting parameters and are expressed in terms of local pressure scale height () and the diffusion coefficient has the unit cm s.
Table 3: Parameters of the best combination of primary and secondary models for each grid of Table 2.

From the results of the binary modelling by Schmid et al. (2015), we can draw an error box on the HRD, within which our best models should reside. We calculated a grid of models over the following parameters that influence the position of the model on the HRD: mass , metallicity , convective-core overshooting in the exponential description, and diffusive mixing (the extent of the grid is given in Table 2), while the initial hydrogen mass fraction was fixed at the solar value .

At every time step of each track, we calculated the asymptotic period spacing (Eq. 3). This parameter, as well as the fundamental model parameters , , , and were then interpolated on an equidistant age scale with steps of years. At every age step we evaluated models for the primary and the secondary simultaneously by comparing the observed to the asymptotic period spacing, allowing only combinations with equal composition at the zero-age main sequence (ZAMS; ) and which fulfil the observed mass ratio within . We calculated the as


where and refer to the mean of the observed period spacing of the zonal dipole modes for the primary and the secondary, respectively. For about of all 19 456 valid combinations and at the same age, meaning that the asymptotic period spacing of the model is within of the observations. When we compare the model tracks to the observations it becomes obvious that those with high diffusion () are too luminous and too old. Thus, we can limit the number of good models even further, when only considering those whose and lie within of the observed values, leaving 60 combinations of primary and secondary models fulfilling all observational constraints derived from the g-mode period spacings and the binarity of our target.

We find that the remaining 60 models have common characteristics. They have low masses (the highest primary mass ; the best models have and ) and low metallicity, with the majority having . While the core overshooting is limited to small values , diffusion with values are possible. As already predicted by the mean value and structure of the observed period spacing, both stars are old with  Gyr and low central hydrogen mass fractions and . The parameters of the best model are given in Table 3.

Based on these common characteristics, we computed two denser grids with refined steps in mass and the mixing parameters (see Table 2). The difference between those grids is the definition of core overshooting, where one grid uses the exponential description and the other uses step overshooting. We analysed them using the same method as described above and found 97 880 and 170 449 models fulfilling the requirements for the exponential-overshoot and step-overshoot grids, respectively. Their parameters agree well with the results of the coarse grid analysis and can be compared in Table 3 for the models of lowest . Figure 3 shows their position and evolutionary tracks in the HRD. In accordance with results from the previous grid, we find that higher values of overshooting combined with low values in envelope diffusion are favoured for the primary, while we find the opposite for the secondary (see Sect. 4 for a discussion on the relevance of these values). Furthermore, we again find that, based on this approach, the best models are old at around 1.1 Gyr and have a low (see Figs 4 and 5).

(a) Fine grid – exponential overshoot
(b) Fine grid – step overshoot
Figure 4: Chi-squared distribution of age for all coeval models in the two denser grids, which fall within of the HRD position of both stars and have equal composition.
(a) Fine grid – exponential overshoot
(b) Fine grid – step overshoot
Figure 5: Chi-squared distribution of central hydrogen mass fraction for the same models as in Fig. 4. The models for the primary are indicated by the dark blue points, while the models for the secondary are indicated by the light blue points.

So far, we have only taken into account the mean of the observed g-mode period spacing, neglecting the effect of rotation and the information that is held in the spacing morphology. In the next step, we therefore calculate the adiabatic frequencies of our models with GYRE and compare their seismic properties to the observations.

3.2 Period spacing with GYRE

We calculated the zonal dipole and quadrupole modes with GYRE for the coeval models with the lowest . As we clearly see the effect of rotation in the observed frequency spectrum (see Fig. 1), we apply Ledoux splitting to the zonal modes to obtain the theoretical prograde and retrograde modes and the rotational frequencies near the core for both components. We derive and and a negligible difference between the models with exponential and those with step overshooting. For these rotation rates and the co-rotating dipole-mode frequencies we find and . The comparison between the calculated modes of Models 2 and 3 and the observed period spacings is illustrated in Fig. 6, while Fig. 7 shows the comparison of the p modes. Our main focus lies in comparing the mean period spacing values and the structure of the period spacing pattern, rather than individual frequencies. We find a satisfactory fit for both models and cannot distinguish between exponential or step overshooting. Furthermore, we are able to assign the series (yellow stars in Fig. 1) to the primary, as there is a clear overlap with the retrograde sectoral modes of both models, while there is a gap for the models of the secondary. From Eq. (3), we expect that the ratio of the asymptotic period spacing of the quadrupole modes to the asymptotic period spacing of the dipole modes is . We use the derived to shift the observed quadrupole sectoral modes to their zonal counterparts, and find  s and . For the secondary and we find  s and . These results support that the series originates in the primary.

(a) Gravity modes of Model 2.
(b) Gravity modes of Model 3.
Figure 6: Comparison of theoretical (empty symbols) to observed (filled symbols) period spacing for the primary (top panels) and the secondary (bottom panels). Ledoux splitting has been applied to the theoretical zonal modes to fit the retrograde and prograde modes. The uncertainties of the observed period spacings are smaller than the symbols. The radial orders are indicated near selected modes. Dipole modes of radial order are shown for the primary and for the secondary, while the radial orders of the quadrupole modes are and , respectively.
(a) Pressure modes of Model 2.
(b) Pressure modes of Model 3.
Figure 7: Comparison of the theoretical zonal p modes (empty symbols) to the observed modes (filled symbols) as a function of degree . Modes of the primary are shown in the top panels and for the secondary in the bottom panels. The radial orders are given above the respective modes.

To test the accuracy of the rotation rates derived with the Ledoux splitting, i.e.  and , we recompute all modes using the TAR. At the rotation rates that were found from the first-order perturbations, we find a disagreement between the GYRE modes and the observations. However, we find that by altering the rotation rates, we cannot obtain a better fit. Rather, we find that the observed shift is due to a difference in the mean period spacing values. We selected Models 2 and 3 as having the lowest for the asymptotic period spacing, which is only valid for zonal modes in a non-rotating case. As the Coriolis force shifts the zonal modes to higher frequencies, we find a decreased period spacing. Therefore, we find a better fit, when selecting younger models (Models 4 and 5; see Table 3) where is about 64 s higher and is about 36 s higher, for both the exponential and step overshooting grids. Figure 8 shows how the period spacings of the best models from the grid search and these younger models compare, using the two different approximations (Ledoux versus TAR) for the puslations. We compare the period spacings of these younger models to the observations in Fig. 9. As can be seen, there is only a slight difference between Fig. 9 and Fig. 6, which mainly results from adopting younger models.

(a) Period spacing for Model 2 and Model 4.
(b) Period spacing for Model 3 and Model 5.
Figure 8: Difference between first-order perturbation approximation and TAR for the rotation rates and . GYRE zonal modes with a Ledoux shift applied are shown as the empty symbols and GYRE modes calculated using the TAR are shown as the filled symbols, where the primary is in dark blue and the secondary in light blue.
(a) Gravity modes of Model 4.
(b) Gravity modes of Model 5.
Figure 9: Same as Fig. 6, but using the TAR instead of the Ledoux splitting to estimate the influence of rotation.

4 Discussion

4.1 On the mixing parameters of the coeval models

Figure 10: Asymptotic period spacing as a function of central hydrogen mass fraction for different masses and with different values of (top panel) and (bottom panel). The red points denote ages of 0.6 Gyr.

During the analysis of all three grids of coeval models based on , , and , we found a clear trend for high overshooting and low diffusive mixing for the primary, and low overshooting and high diffusive mixing for the secondary. The main reasons why these parameters influence our model selection is that stronger mixing increases the main sequence lifetime of a star by providing fresh hydrogen to the nuclear-burning core. By directly affecting the stellar structure they also alter the Brunt-Väisälä frequency and thus . Both effects are illustrated in Fig. 10. In this figure it can also be seen that higher mass models have a higher period spacing at a given . Yet, the evolution of is faster for the primary and drops to the observed value in shorter time than for the secondary, when the strength of mixing is equal in both components. Thus, a higher value of for the primary is needed to slow the evolution of down. Additionally, it extends the evolutionary track to reach the observed , when it would otherwise turn off before the observed error box. For the secondary, additional mixing is mainly required to adjust the evolution of with that of the primary. As models with and already cross the observational error box on the HRD, increasing would result in a too cool and too big star. Diffusive mixing, on the other hand, has a weaker influence on the main sequence lifetime than overshooting and is thus the preferred mixing mechanism for the secondary. In summary we state that the equal age requirement and the stellar parameters from binary modelling, serve as a stringent constraint to derive the width of the core-overshooting region and the strength of radial envelope mixing in both stars at the current evolutionary stage. Variations throughout the stellar evolution of these parameters are likely to occur. Both stars could thus have had similar values near the ZAMS or started off with different values at birth.

4.2 The period spacing morphology

When focusing only on the morphology of the period spacing we find Models 2 and 3 to be less satisfying fits to the observations than suggested by the mean period spacing values. While the computed period spacing for the primary has variations that are too strong and too short, the structure of the secondary model is too smooth (see Figs. 6 and 9). Miglio et al. (2008) showed that the interval of the dips in period spacing depends on the steepness of the chemical composition gradient near the core, which grows with increasing age or decreasing as more and more hydrogen is burned into helium. In Fig. 11 we show how the period spacing evolves for different values of for a model of . For the secondary, the lack of structure in the computed period spacing can be explained by high diffusion, which diminishes near the core. In Fig. 12 it is clearly visible how the dips in period spacing disappear with increasing . Hence, the period spacing morphology points to stars with a higher than we found with the grid search, explained in Sect. 3. Indeed, we find a better fit with models where and , while we find that overshooting only has weak influence. The similarity of these values for both stars is a direct consequence of the similar morphology in the observed period spacings of both components, however, they are in contradiction to the other observational constraints. By increasing the stars move closer to the ZAMS and outside the error box defined by and from binary modelling. In even stronger disagreement, however, are the asymptotic period spacings, as they are several above the observed values.

(a) , , ,
(b) , , ,
Figure 11: Period spacing for different values of along two evolutionary tracks incorporating exponential core-overshooting (left panel) or step overshooting (right panel). The radial orders of the first and last mode in the period spacing is given. The filling of the symbols is alternated to enhance visibility.
(a) , , ,
(b) , , ,
Figure 12: Period spacing for different values of along two evolutionary tracks, one with exponential core-overshooting (left panel) and one with step overshooting (right panel). The radial orders of the first and last mode in the period spacing is given.

To obtain better agreement between the period spacing morphology and the mean period spacing values, as well as the stellar parameters, we exploit correlations between mass, metallicity, and . For a given , the asymptotic period spacing decreases with decreasing mass and decreasing metallicity. Moreover, models with low metal content are more luminous than those with high metal content and, therefore, mimic higher mass models with high metallicity in the HRD. To decrease as much as possible, we use a secondary mass of and determine the primary mass via the mass ratio as . We choose a metallicity of , which is still in agreement with the slightly sub-solar metallicity reported by Schmid et al. (2015). The mixing parameters were chosen to conform with the observed structure in the measured period spacing pattern. From the computed evolutionary tracks we determine the coeval Model 6, where is at a minimum (all grid and global parameters of Model 6 can be found in Table 3). For this model, we compute the zonal dipole and quadrupole modes with GYRE and apply the Ledoux splitting for the prograde and retrograde modes. The derived rotation rates for this model, are and . We show the result for this test in Fig. 13, where the computed g and p modes are compared to the observed modes. We find a much better fit of the period spacing morphology and even the asymptotic period spacing values lie within of the observed mean values (Fig. (a)a). Some room for improvement exists for the global parameters of the primary model, as can be seen in Fig. 14. The and of the model still lie within of the observed values, but its radius lies outside this range.

(a) Gravity modes of Model 6.
(b) Pressure modes of Model 6.
Figure 13: Comparison of theoretical modes (empty symbols) to observed modes (filled symbols) for the primary (top panels) and the secondary (bottom panels). The radial orders are given near selected modes.
Figure 14: Same as Fig. 3 but for Model 6, indicated by stars and its evolutionary track as the black solid line. Here the areas for and are also shown as dashed boxes.

4.3 Varying the He mass fraction

In the model grids discussed above, we have assumed an initial hydrogen mass fraction . Varying and thereby adapting the initial helium mass fraction influences the positions of the tracks on the HRD and the asymptotic period spacing. Increasing , while decreasing causes a star of the same mass and to become more luminous. Furthermore, for a given , increases. We therefore tested whether changing the helium mass fraction would lead to higher masses that are in better agreement with the values from binary modelling, following the same modelling procedure as carried out for fixed . We computed two additional small grids of models for and , using exponential overshooting. The best fitting model with has a primary mass of , which is slightly higher than the mass of Model 2. However, this increase in mass is accompanied with a lower metallicity of . The secondary mass is the same as for Model 2, while the overshooting and diffusion parameters only change slightly (, , , and ).

4.4 The p modes

The observed p modes provide weaker constraints for both components than the g modes because of their uncertain mode identification, which is based on the assumption that all observed multiplets are complete. Models 2 and 3, which best fit the observed mean period spacing and binary parameters, confirm the mode identification for the secondary (Fig. 7). The singlet and highest amplitude p mode fits the first overtone mode better than the radial fundamental. For the primary, the situation is somewhat more ambiguous. For the singlet and the triplet centre frequency , neither nor can be assigned without doubt. As for the g modes, it is not possible to distinguish between exponential and step overshooting using the p modes.

Schmid et al. (2015) used the rotational splittings observed in the p modes and the g modes to provide an estimate of the core-to-surface rotation rates. They showed that both stars might have a slight differential rotation in the radial direction. We now use our best models to compute the Ledoux constant and derive the envelope rotation rates based on the p-mode splittings. Since the splittings detected by Schmid et al. (2015) show a significant asymmetry, which cannot be reproduced by the first-order Ledoux approximation, we use the average per multiplet.

For the triplet associated with the secondary (see Fig. 2) we find that is close to 1, which is the expected value for pure p modes. This yields an envelope rotation rate of for the secondary. Thus, we observe differential rotation in the radial direction with a slower core than envelope, given the ratio of surface-to-core rotation rate . This value is equal for Models 2 and 3.

For the primary, the situation is slightly more complicated. Using for the modes of radial orders and for the dipole modes and for the quadrupole mode, we derive rotation rates in the range . The envelope of the primary can thus be either faster or slower than its core, depending on the p-mode multiplet. When we examine the values of the individual modes more closely, we find that some have values below 0.9. For the mode near we find and a rotation rate of , while and hence . This is because the mode has a mixed mode character and not only probes the envelope (like pure p modes) but also has a high amplitude near the convective core (like pure g modes). This is also the case for the quadrupole mode, where we find envelope rotation rates of . To obtain a better estimate of the surface-to-core rotation rate, we only use the rotation rates derived from the pure p-mode triplet near . For this triplet ranges from to , using either modes of or of both Models 2 and 3. Thus, we find almost uniform rotation with a slightly faster core than envelope, with up to .

Moreover, we use the p modes to compare the scenario that predicts the mean period spacing and the binary parameters (Models 2 and 3) to the scenario that predicts the period spacing morphology (Model 6). We find that the p modes of Model 6 have frequencies that are too high and therefore do not agree with the observations. Aside from the binary parameters, the p-mode frequencies are thus another indication that Models 2 and 3 are indeed better representations for the two stars of KIC 10080943 than Model 6.

5 Summary

In this paper we have presented detailed seismic modelling for two F-type g- and p-mode hybrid pulsators that reside in the binary system KIC 10089043 (Keen et al., 2015; Schmid et al., 2015). We calculated stellar model grids with MESA, covering the observed range of stellar parameters, and compared the observed pulsation modes to predictions by the stellar pulsation code GYRE. For the analysis, we were able to exploit observational constraints from the binarity and the pulsations. This was implemented by requiring an equal age and composition, the mass ratio and the position in the HRD within of the binary solution, together with asymptotic period spacing values within of the observed mean period spacing. With this approach, we provide the first such consistent seismic modelling of a binary F-type g-mode pulsator. We find that the stars have an age of  Gyr and have a low central hydrogen mass fraction of and . This is a consequence of the observed mean period spacing values, which are below 3000 s and, thus, rather low for stars. In fact, our best models have masses below the binary values. However, the age is in good agreement with the observed HRD positions of both stars with cool temperatures and big radii. Furthermore, we were able to determine the rotation rates of the regions near the core to be and . When the traditional approximation of rotation is used at these rotation rates instead of the first-order perturbative Ledoux approximation, a lower asymptotic period spacing is required for both stars. This is because the zonal modes are influenced by the Coriolis force and shifted to higher frequencies.

To test the shape of the convective-core overshooting we computed two grids that differ only in the overshooting prescription. One grid uses exponential core overshooting, while the other uses step overshooting. We could not find a significant difference between these two grids, in the sense that the difference in morphology of the period spacing pattern between these two descriptions is by far inferior to the discrepancy between the observed and theoretical morphology. We could show, however, that the stellar parameters derived from binary modelling and the requirement of equal age and composition for both stars, can constrain the amount of overshooting and diffusive mixing. For the primary, we find high values of overshooting and low values of diffusion and for the secondary we find the opposite. The best model for the exponential overshooting grid is , , , , , , and ; is expressed in local pressure scale height and has a unit of cm s. For the step overshooting grid, the best model is similar at , , , , , , and .

These best models fail to predict the detailed period spacing structure, which points to the limitations of 1D stellar models in terms of chemical mixing in the radial direction only. For these 1D stellar models, we find a better agreement with the observed structures for stars that have a higher and more similar mixing parameters, which contradicts the observed mean period spacing and stellar parameters. We are able to find a coeval model with low masses and metallicity that fits both the morphology and mean value of the period spacing. Yet, the p-mode frequencies of this model are too high to reproduce the observed frequencies. This points to mean densities that are too high or radii that are too small. This model, indeed, is discrepant with the binary modelling at level. Especially, the primary radius of the model is below the observed value. An independent measurement of the radius, for example through a precise distance provided by the ESA Gaia mission (Perryman et al., 2001), could help to firmly exclude such younger models. The fact that the variations of the period spacing of the primary are too strong clearly points towards missing extra mixing in the transversal direction, which is perhaps of rotational origin.

During the analysis for this paper, we ignored the influence both stars have on each other and modelled them as coeval single stars. This approach is justified, as the observed interaction is weak and we relied on the g modes, whose probing power is situated near the core. We can estimate the tidal influence by the tide-generating potential, which depends on the small parameter (e.g. Willems & Aerts, 2002) and find for KIC 10080943. The influence of the centrifugal force due to the rotation on the pulsations can be estimated by , which lies between 0.01 and 0.017 for the primary and between 0.005 and 0.014 for the secondary (see e.g. Aerts et al., 2002). Thus, the centrifugal force has a much stronger influence on the pulsations than the tides, but is so small that our approach based on the Coriolis force above is fully justified. Also, the brightening signal, which is detected in the high-precision Kepler light curve, has a small amplitude and the reflection effect outweighs the contribution of the tidal distortion. Some evidence of the tidal forces, however, could be detected in one g mode, which has a frequency of exactly , and in the rotational splitting of the p modes (Schmid et al., 2015). The next step towards improving the modelling of this binary would be to create a binary model that takes these interactions into account. Such models could also help gauge the effect the tidal forces have on the level of mixing or the convective-core size. The best models we obtained so far constitute a fruitful starting point for more detailed studies such as these. Other future work could include mode excitation. As of now, GYRE does not include a pulsation-convection interaction, which would be necessary to study the mode excitation by convective flux blocking (Guzik et al., 2000; Dupret et al., 2005). Therefore, KIC 10080943 is a rich target, which still holds potential for improving stellar and seismic models.

The research leading to these results received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement N 670519: MAMSIE). The computational resources and services used in this work were provided by the VSC (Flemish Supercomputer Center), funded by the Hercules Foundation and the Flemish Government – department EWI. We would like to thank Ehsan Moravveji, Péter I. Pápics, and Timothy Van Reeth for fruitful discussions and valuable input, and the anonymous referee for comments that helped improve this manuscript.


  • Aerts et al. (2010) Aerts, C., Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology, 1st edn. (Springer), 866
  • Aerts et al. (2002) Aerts, C., Handler, G., Arentoft, T., et al. 2002, MNRAS, 333, L35
  • Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
  • Ballot et al. (2010) Ballot, J., Lignières, F., Reese, D. R., & Rieutord, M. 2010, A&A, 518, A30
  • Borucki et al. (2010) Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977
  • Bouabid et al. (2013) Bouabid, M.-P., Dupret, M.-A., Salmon, S., et al. 2013, MNRAS, 429, 2500
  • Breger (1990) Breger, M. 1990, Delta Scuti Star Newsletter, 2, 13
  • Cox & Giuli (1968) Cox, J. P. & Giuli, R. T. 1968, Principles of stellar structure (New York: Gordon and Breach)
  • Dupret et al. (2005) Dupret, M.-A., Grigahcène, A., Garrido, R., Gabriel, M., & Scuflaire, R. 2005, A&A, 435, 927
  • Guzik et al. (2000) Guzik, J. A., Kaye, A. B., Bradley, P. A., Cox, A. N., & Neuforge, C. 2000, ApJ, 542, L57
  • Hambleton et al. (2013) Hambleton, K. M., Kurtz, D. W., Prša, A., et al. 2013, MNRAS, 434, 925
  • Herwig (2000) Herwig, F. 2000, A&A, 360, 952
  • Iglesias & Rogers (1993) Iglesias, C. A. & Rogers, F. J. 1993, ApJ, 412, 752
  • Iglesias & Rogers (1996) Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943
  • Keen et al. (2015) Keen, M. A., Bedding, T. R., Murphy, S. J., et al. 2015, MNRAS, 454, 1792
  • Kurtz et al. (2014) Kurtz, D. W., Saio, H., Takata, M., et al. 2014, MNRAS, 444, 102
  • Ledoux (1951) Ledoux, P. 1951, ApJ, 114, 373
  • Maceroni et al. (2014) Maceroni, C., Lehmann, H., da Silva, R., et al. 2014, A&A, 563, A59
  • Maceroni et al. (2009) Maceroni, C., Montalbán, J., Michel, E., et al. 2009, A&A, 508, 1375
  • Miglio et al. (2008) Miglio, A., Montalbán, J., Noels, A., & Eggenberger, P. 2008, MNRAS, 386, 1487
  • Moravveji et al. (2015) Moravveji, E., Aerts, C., Papics, P. I., Andres Triana, S., & Vandoren, B. 2015, A&A, 580, A27
  • Moravveji et al. (2016) Moravveji, E., Townsend, R. H. D., Aerts, C., & Mathis, S. 2016, ApJ, 823, 130
  • Pápics et al. (2015) Pápics, P. I., Tkachenko, A., Aerts, C., et al. 2015, ApJ, 803, L25
  • Paxton et al. (2011) Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3
  • Paxton et al. (2013) Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4
  • Paxton et al. (2015) Paxton, B., Marchant, P., Schwab, J., et al. 2015, ApJS, 220, 15
  • Perryman et al. (2001) Perryman, M. A. C., de Boer, K. S., Gilmore, G., et al. 2001, A&A, 369, 339
  • Saio et al. (2015) Saio, H., Kurtz, D. W., Takata, M., et al. 2015, MNRAS, 447, 3264
  • Schmid et al. (2015) Schmid, V. S., Tkachenko, A., Aerts, C., et al. 2015, A&A, 584, A35
  • Southworth (2012) Southworth, J. 2012, in Orbital Couples: Pas de Deux in the Solar System and the Milky Way, ed. F. Arenou & D. Hestroffer, 51–58
  • Stellingwerf (1979) Stellingwerf, R. F. 1979, ApJ, 227, 935
  • Tassoul (1980) Tassoul, M. 1980, ApJS, 43, 469
  • Townsend (2003) Townsend, R. H. D. 2003, MNRAS, 343, 125
  • Townsend & Teitler (2013) Townsend, R. H. D. & Teitler, S. A. 2013, MNRAS, 435, 3406
  • Van Reeth et al. (2015a) Van Reeth, T., Tkachenko, A., Aerts, C., et al. 2015a, A&A, 574, A17
  • Van Reeth et al. (2015b) Van Reeth, T., Tkachenko, A., Aerts, C., et al. 2015b, ApJS, 218, 27
  • Welsh et al. (2011) Welsh, W. F., Orosz, J. A., Aerts, C., et al. 2011, ApJS, 197, 4
  • Willems & Aerts (2002) Willems, B. & Aerts, C. 2002, A&A, 384, 441

Appendix A MESA input file

! inlist for KIC10080943 grid calculations
! evolve main sequence model starting from ZAMS
! Parameters input via run_star_extras.f
! not in inlist
! Input via run_star_extras.f:
! new_Z, new_Y, initial_mass, min_D_mix, and
! step_overshoot_f_above_burn_h or
! overshoot_f_above_burn_h


  ! history and profile columns
    history_columns_file = ’hist.list’
    profile_columns_file = ’prof.list’

!!! begin with a pre-main sequence model
    create_pre_main_sequence_model = .false.

!!! start from ZAMS:
    relax_Z = .true.
    change_Z = .true.
    relax_initial_Z = .true.
    change_initial_Z = .true.
    ! new_Z set via run_star_extras.f

    relax_Y = .true.
    change_Y = .true.
    relax_initial_Y = .true.
    change_initial_Y = .true.
    ! new_Y set via run_star_extras.f

!!! Composition and opacities (Asplund09+OPAL):
    kappa_file_prefix = ’a09’
    kappa_lowT_prefix = ’lowT_fa05_a09p’
    initial_zfracs = 6

    change_lnPgas_flag = .true.
    change_initial_lnPgas_flag = .true.
    new_lnPgas_flag = .true.

    change_net = .true.
    new_net_name = ’’
    change_initial_net = .true.

  ! display on-screen plots
    pgstar_flag = .false.

/ !end of star_job namelist


!!! Starting specifications
    ! initial_mass set via run_star_extras.f

!!! Atmosphere
    which_atm_option = ’photosphere_tables’

!!! Mixing
  ! Diffusive mixing
  ! Minimal mixing in radiative zone
    set_min_D_mix = .true.
    ! min_D_mix set via run_star_extras.f
    remove_small_D_limit = 0

  ! Convection MLT
    mixing_length_alpha = 1.8

    use_Ledoux_criterion = .false.

  ! Over- and undershooting
    D_mix_ov_limit = 5d-2

    overshoot_f0_above_burn_h = 0.001

  ! Exponential overshooting and coarse grid
    ! overshoot_f_above_burn_h set
    ! via run_star_extras.f

  ! Step function overshooting grid
    ! step_overshoot_f_above_burn_h set
    ! via run_star_extras.f

    step_overshoot_D = 0
    step_overshoot_D0_coeff = 1

    num_cells_for_smooth_brunt_B = 0
    num_cells_for_smooth_gradL_composition_term = 0

!!! Time step control
    varcontrol_target = 2.5d-5

    max_age = 2d10
    max_years_for_timestep = 1d7

  ! Stop condition
  ! stop when the center mass fraction of h1
  ! drops below this limit
    xa_central_lower_limit_species(1) = ’h1’
    xa_central_lower_limit(1) = 1d-3
    when_to_stop_rtol = 1d-3
    when_to_stop_atol = 1d-3

!!! Mesh grid
    max_allowed_nz = 10000
    mesh_delta_coeff = 0.5
    cubic_interpolation_in_Z = .true.

    xa_function_species(1) = ’he4’
    xa_function_weight(1) = 100
    xa_function_param(1) = 1d-2

    xa_function_species(2) = ’he3’
    xa_function_weight(2) = 100
    xa_function_param(2) = 1d-5

    mesh_dlogX_dlogP_extra(:) = 0.1         ! resol coeff for chemical gradients
    mesh_dlogX_dlogP_full_on(:) = 2
    mesh_dlogX_dlogP_full_off(:) = 1

    mesh_logX_species(1) = ’he4’            ! taking into account abundance of He4
    mesh_logX_min_for_extra(1) = -6         ! for abundances larger than ~
    mesh_logX_species(2) = ’n14’            ! taking into account abundance of N14
    mesh_logX_min_for_extra(2) = -6
    mesh_logX_species(3) = ’c12’
    mesh_logX_min_for_extra(3) = -6

    xtra_coef_czb_full_on = 1.0             ! use this coef if center_he4 below this
    xtra_coef_czb_full_off = 1.0            ! i.e. always on MS

    xtra_coef_a_l_nb_czb = 0.1
    xtra_dist_a_l_nb_czb = 10.0

    xtra_coef_b_l_nb_czb = 0.1
    xtra_dist_b_l_nb_czb = 3.0

    okay_to_remesh = .true.

!!! Input and Output
    max_num_profile_models = 500
    profile_interval = 100
    history_interval = 1
    terminal_cnt = 10
    write_header_frequency = 5

  ! GYRE output
    write_pulse_info_with_profile = .true.
    pulse_info_format = ’GYRE’
    add_atmosphere_to_pulse_info = .true.

/ ! end of controls namelist
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description