Asteroseismic modelling of the two Ftype hybrid pulsators KIC 10080943A and KIC 10080943B
^{1} Institute for Astronomy, KU Leuven Celestijnenlaan 200D 3001 Leuven Belgium
email: valentina.schmid@ster.kuleuven.be
^{2} Department of Astrophysics/IMAPP, Radboud University Nijmegen, PO Box 9010, 6500 GL Nijmegen, The Netherlands
Received: 31 March 2016
Accepted: 22 May 2016
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. Highorder gravity modes are sensitive to the conditions near the convective core and therefore allow for a determination of parameters describing interior physics, especially the convectivecore overshooting parameter. KIC 10080943 is a binary system that contains two gravity and pressuremode 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 gmode 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 HertzsprungRussell 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 gmode period spacing. We find that the amount of core overshooting and diffusive mixing can be well constrained by the equalage 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 firstorder 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 firstorder perturbation approximation ignores.
Key words: stars: individual: KIC 10080943 / binaries: spectroscopic / stars: oscillations / stars: evolution
© ESO, 2016
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 highorder gmodes of a nonrotating, 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 (1.5 ≤ M (M_{⊙}) ≤ 2.2) 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 gmodes are thought to be excited by the convective flux blocking mechanism (Guzik et al. 2000; Dupret et al. 2005). The heatdriven pmodes 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 gmode 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. Doublelined 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 doublelined 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 onedimensional (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
Fig. 1 Observed gmodes in frequency (top) and period (bottom) as a function of degree ℓ are shown for the primary in red and the secondary in blue. The m values are distinguished by the different symbols. The yellow stars denote the ℓ = 2 modes, which could not be assigned to either component. 

Open with DEXTER 
Fig. 2 Observed pmode 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 m values are distinguished by the different symbols. 

Open with DEXTER 
KIC 10080943 is an eccentric, doublelined binary system, containing two γ Dor/δ Sct hybrid pulsators (for a detailed description of γ Dor and δ Sct type pulsators, see Aerts et al. 2010, chap. 2). Keen et al. (2015) published the analysis of the period spacings of the gmodes 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 highresolution followup spectroscopy covering the orbit. They were able to assign most period spacing series and highamplitude pmodes to either the primary or the secondary. This was achieved by analysing the phase modulations of the pmodes that arise from the binary motion. Modes of the primary have a time delay of 36.2 ± 4.8 s, while those of the secondary are delayed by 44.0 ± 3.5 s. Combination frequencies between the pmodes 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 ℓ = 2 modes could not be assigned to either star, as it does not show combination frequencies or rotational splitting. In Fig. 2 the pmodes 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 HertzsprungRussell diagram (HRD) inside a 1σ error box in Fig. 3. It is obvious that the parameters lack the precision typical for doublelined 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.
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 highorder gmodes. For both components, the mean period spacing value is rather low at ΔP< 3000 s for stars near 2 M_{⊙}. Indeed, comparing with the sample of γ Dor stars in Van Reeth et al. (2015b), the measured ΔP values of both components hint at masses <2 M_{⊙}. 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 hydrogenburning core has left behind a chemical composition gradient ∇_{μ} leading to a peak in the BruntVäisälä frequency N and therefore a decrease in the asymptotic period spacing (e.g. Miglio et al. 2008).
In addition to the gmodes, two pmodes were identified as candidate radial modes and can help distinguish between possible models. To estimate their radial overtone, we calculate the pulsation constant Q (Breger 1990) as (1)or, expressed in observable parameters (2)For the highest amplitude pmode (f_{s,2} = 13.947586 d^{1}) and for the parameters of the secondary, we find Q = 0.0332 ± 0.0045, which is indicative of the radial fundamental mode. The pmode with frequency f_{s,1} = 15.68333 d^{1} belongs to the primary and has Q = 0.018 ± 0.001, which is consistent with 3 ≤ n ≤ 5 (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 α_{MLT} = 1.8. All equilibrium models along the evolutionary tracks are nonrotating, as both components are slow rotators (rotational period near the core P_{c,1} ≈ 7 d and P_{c,2} ≈ 11 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 BruntVäisälä frequency N 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 N and hence the periods and period spacings of highorder gmodes. This can be illustrated by the asymptotic period spacing, which is defined by Tassoul (1980) as (3)where , x = r/R and x_{0} is the convectivecore boundary. To probe the strength of mixing beyond the convective regions, we included convectivecore overshooting and diffusive mixing as free parameters into our models. Convectivecore overshooting mainly increases the size of the convective core, thus decreasing the integral ^{∫}  N  /x dx. It is implemented in MESA using two different definitions that both depend on the MLT diffusion coefficient (D_{conv}) and an adjustable parameter. In the exponential description, the region immediately outside the core has a diffusion coefficient as determined by (4)where H_{P} is the local pressure scale height, z is the distance from the convectivecore boundary, and f_{ov} is a free parameter (Herwig 2000). This definition results in an exponentially decaying D_{ov}. Alternatively overshooting can also be defined as a step function, where the extent of the overshooting layer is d_{ov} = α_{ov}H_{P} with the adjustable parameter α_{ov}. Typically α_{ov} ≃ 10f_{ov} 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 D_{mix} of 1 cm^{2} s^{1}. 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 X_{c} at a given age when only the mixing strength is increased.
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 firstorder perturbative approach Δf = mβ_{n,ℓ}Ω (Ledoux 1951), where β_{n,ℓ} = 1−C_{n,ℓ} and C_{n,ℓ} is the Ledoux constant. We use β_{n,ℓ} as calculated by GYRE and the observed frequency splitting Δf 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) 0.14 < 2Ω_{c}/f_{co}< 0.27, where Ω_{c} is the rotation rate near the core and f_{co} are the gmode frequencies in the corotating 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 gmodes in γ Dor stars. Bouabid et al. (2013) studied the influence of the Coriolis force on highorder gmodes, using the TAR, and compared their computations to predictions by the firstorder 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
Fig. 3 Top panel: HRD with the observed positions and 1σ 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 ΔP as a function of central hydrogen mass fraction X_{c}. The observed mean period spacing of the primary (secondary) is given as the horizontal solid (dashed) line and the 1σ areas are shown as the grey patch. 

Open with DEXTER 
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 M, metallicity Z, convectivecore overshooting f_{ov} in the exponential description, and diffusive mixing D_{mix} (the extent of the grid is given in Table 2), while the initial hydrogen mass fraction was fixed at the solar value X_{i} = 0.71.
At every time step of each track, we calculated the asymptotic period spacing ΔP_{c} (Eq. (3)). This parameter, as well as the fundamental model parameters log (L/L_{⊙}), log (R/R_{⊙}), log T_{eff}, and X_{c} were then interpolated on an equidistant age scale with steps of Δa = 10^{7} yr. 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 zeroage main sequence (ZAMS; Z_{1} = Z_{2}) and which fulfil the observed mass ratio M_{2}/M_{1} = q = 0.96 ± 0.01 within 2σ. We calculated the as (5)where ΔP_{o,1} and ΔP_{o,2} refer to the mean of the observed period spacing of the zonal dipole modes for the primary and the secondary, respectively. For about 4.1% of all 19 456 valid combinations and at the same age, meaning that the asymptotic period spacing of the model is within 2σ of the observations. When we compare the model tracks to the observations it becomes obvious that those with high diffusion (log D_{mix} = 3) are too luminous and too old. Thus, we can limit the number of good models even further, when only considering those whose log T_{eff} and log (R/R_{⊙}) lie within 2σ of the observed values, leaving 60 combinations of primary and secondary models fulfilling all observational constraints derived from the gmode 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 M_{1} = 1.9 M_{⊙}; the best models have M_{1} = 1.8−1.85 M_{⊙} and M_{2} = 1.7−1.75 M_{⊙}) and low metallicity, with the majority having Z = 0.012. While the core overshooting is limited to small values f_{ov} = 0.0−0.1, diffusion with values log D_{mix} = 1−2 are possible. As already predicted by the mean value and structure of the observed period spacing, both stars are old with a ≈ 1 Gyr and low central hydrogen mass fractions X_{c,1} = 0.05−0.15 and X_{c,2} = 0.15−0.3. 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 exponentialovershoot and stepovershoot 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 χ^{2}. 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 χ^{2} approach, the best models are old at around 1.1 Gyr and have a low X_{c} (see Figs. 4 and 5).
Fig. 4 Chisquared distribution of age for all coeval models in the two denser grids, which fall within 2σ of the HRD position of both stars and have equal composition. 

Open with DEXTER 
Fig. 5 Chisquared distribution of central hydrogen mass fraction X_{c} 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. 

Open with DEXTER 
So far, we have only taken into account the mean of the observed gmode 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 χ^{2}. 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 Ω_{c} near the core for both components. We derive Ω_{c,1} = 0.13963 ± 0.00007 d^{1} and Ω_{c,2} = 0.09034 ± 0.00002 d^{1} and a negligible difference between the models with exponential and those with step overshooting. For these rotation rates and the corotating dipolemode frequencies we find 0.198 < 2Ω_{c,1}/f_{co}< 0.264 and 0.144 < 2Ω_{c,2}/f_{co}< 0.237. 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 pmodes. 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 ℓ = 2 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 ΔP_{c}(ℓ = 2) / ΔP_{c}(ℓ = 1) = 0.577. We use the derived Ω_{c,1} to shift the observed quadrupole sectoral modes to their zonal counterparts, and find ΔP_{o}(ℓ = 2) = 1664 ± 21 s and ΔP_{o}(ℓ = 2) / ΔP_{o,1}(ℓ = 1) = 0.591 ± 0.009. For the secondary and Ω_{c,2} we find ΔP_{o}(ℓ = 2) = 1840 ± 14 s and ΔP_{o}(ℓ = 2) / ΔP_{o,2}(ℓ = 1) = 0.633 ± 0.009. These results support that the ℓ = 2 series originates in the primary.
Fig. 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 k are indicated near selected modes. Dipole modes of radial order −31 ≤ k ≤ −20 are shown for the primary and −37 ≤ k ≤ −25 for the secondary, while the radial orders of the quadrupole modes are −28 ≤ k ≤ −35 and −29 ≤ k ≤ −36, respectively. 

Open with DEXTER 
Fig. 7 Comparison of the theoretical zonal pmodes (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. 

Open with DEXTER 
To test the accuracy of the rotation rates derived with the Ledoux splitting, i.e. Ω_{c,1} = 0.13963 and Ω_{c,2} = 0.09034, we recompute all modes using the TAR. At the rotation rates that were found from the firstorder 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 χ^{2} for the asymptotic period spacing, which is only valid for zonal modes in a nonrotating 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 ΔP_{c,1} is about 64 s higher and ΔP_{c,2} 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 Figs. 9 and 6, which mainly results from adopting younger models.
Fig. 8 Difference between firstorder perturbation approximation and TAR for the rotation rates Ω_{c,1} = 0.13963 d^{1} and Ω_{c,2} = 0.09034 d^{1}. 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. 

Open with DEXTER 
Fig. 9 Same as Fig. 6, but using the TAR instead of the Ledoux splitting to estimate the influence of rotation. 

Open with DEXTER 
4. Discussion
4.1. On the mixing parameters of the coeval models
During the analysis of all three grids of coeval models based on ΔP_{c}, T_{eff}, and R, 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 nuclearburning core. By directly affecting the stellar structure they also alter the BruntVäisälä frequency N and thus ΔP_{c}. 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 X_{c}. Yet, the evolution of ΔP_{c} 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 f_{ov} for the primary is needed to slow the evolution of ΔP_{c} down. Additionally, it extends the evolutionary track to reach the observed T_{eff,1}, when it would otherwise turn off before the observed error box. For the secondary, additional mixing is mainly required to adjust the evolution of ΔP_{c} with that of the primary. As models with f_{ov} = 0.0 and log (D_{mix}) = 0.0 already cross the observational error box on the HRD, increasing f_{ov} 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 coreovershooting 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.
Fig. 10 Asymptotic period spacing ΔP as a function of central hydrogen mass fraction X_{c} for different masses and with different values of f_{ov} (top panel) and D_{mix} (bottom panel). The red points denote ages of 0.6 Gyr. 

Open with DEXTER 
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 X_{c} as more and more hydrogen is burned into helium. In Fig. 11 we show how the period spacing evolves for different values of X_{c} for a model of M = 1.81 M_{⊙}. 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 log (D_{mix}). Hence, the period spacing morphology points to stars with a higher X_{c} than we found with the grid search, explained in Sect. 3. Indeed, we find a better fit with models where X_{c} ≈ 0.35 and 0.5 ≤ log (D_{mix}) ≤ 0.75, 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 X_{c} the stars move closer to the ZAMS and outside the error box defined by T_{eff} and R from binary modelling. In even stronger disagreement, however, are the asymptotic period spacings, as they are several σ above the observed values.
Fig. 11 Period spacing for different values of X_{c} along two evolutionary tracks incorporating exponential coreovershooting (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. 

Open with DEXTER 
Fig. 12 Period spacing for different values of log (D_{mix}) along two evolutionary tracks, one with exponential coreovershooting (left panel) and one with step overshooting (right panel). The radial orders of the first and last mode in the period spacing is given. 

Open with DEXTER 
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 ΔP_{c}. For a given X_{c}, 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 ΔP_{c} as much as possible, we use a secondary mass of M_{2}−3σ_{M2} = 1.6 M_{⊙} and determine the primary mass via the mass ratio as M_{1} = M_{2}/q = 1.67 M_{⊙}. We choose a metallicity of Z = 0.01, which is still in agreement with the slightly subsolar 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 Ω_{c,1} = 0.13969 ± 0.00003 d^{1} and Ω_{c,2} = 0.09047 ± 0.00004 d^{1}. We show the result for this test in Fig. 13, where the computed g and pmodes 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 2σ of the observed mean values (Fig. 13a). Some room for improvement exists for the global parameters of the primary model, as can be seen in Fig. 14. The T_{eff} and log (L/L_{⊙}) of the model still lie within 3σ of the observed values, but its radius R = 2.12 R_{⊙} lies outside this range.
Fig. 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. 

Open with DEXTER 
Fig. 14 Same as Fig. 3 but for Model 6, indicated by stars and its evolutionary track as the black solid line. Here the 2σ areas for T_{eff,1,2} and log (L_{1,2}/L_{⊙}) are also shown as dashed boxes. 

Open with DEXTER 
4.3. Varying the He mass fraction
In the model grids discussed above, we have assumed an initial hydrogen mass fraction X_{i} = 0.71. Varying X_{i} and thereby adapting the initial helium mass fraction Y_{i} influences the positions of the tracks on the HRD and the asymptotic period spacing. Increasing Y_{i}, while decreasing X_{i} causes a star of the same mass and Z to become more luminous. Furthermore, for a given X_{c}, ΔP_{c} 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 X_{i} = 0.71. We computed two additional small grids of models for X_{i} = 0.7 and X_{i} = 0.72, using exponential overshooting. The best fitting model with X_{i} = 0.72 has a primary mass of M_{1} = 1.83 M_{⊙}, which is slightly higher than the mass of Model 2. However, this increase in mass is accompanied with a lower metallicity of Z = 0.0115. The secondary mass is the same as for Model 2, while the overshooting and diffusion parameters only change slightly (f_{ov,1} = 0.01, log D_{mix,1} = 0, f_{ov,2} = 0.005, and log D_{mix,2} = 1.75).
4.4. The pmodes
The observed pmodes provide weaker constraints for both components than the gmodes 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 pmode f_{s,2} fits the first overtone mode better than the radial fundamental. For the primary, the situation is somewhat more ambiguous. For the singlet f_{s,1} and the triplet centre frequency f = 19.51028 d^{1}, neither ℓ nor k can be assigned without doubt. As for the gmodes, it is not possible to distinguish between exponential and step overshooting using the pmodes.
Schmid et al. (2015) used the rotational splittings observed in the pmodes and the gmodes to provide an estimate of the coretosurface 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 C_{n,ℓ} and derive the envelope rotation rates based on the pmode splittings. Since the splittings detected by Schmid et al. (2015) show a significant asymmetry, which cannot be reproduced by the firstorder Ledoux approximation, we use the average Δf per multiplet.
For the triplet associated with the secondary (see Fig. 2) we find that β_{n,ℓ} = (1−C_{n,ℓ}) is close to 1, which is the expected value for pure pmodes. This yields an envelope rotation rate of Ω_{e,2} = 0.122 ± 0.002 d^{1} for the secondary. Thus, we observe differential rotation in the radial direction with a slower core than envelope, given the ratio of surfacetocore rotation rate Ω_{e,2}/ Ω_{c,2} = 1.35 ± 0.02. This value is equal for Models 2 and 3.
For the primary, the situation is slightly more complicated. Using β_{n,ℓ} for the modes of radial orders n = 1 and n = 4 for the dipole modes and n = 1 for the quadrupole mode, we derive rotation rates in the range 0.137 ≤ Ω_{e,1} (d^{1}) ≤ 0.173. The envelope of the primary can thus be either faster or slower than its core, depending on the pmode multiplet. When we examine the β_{n,ℓ} values of the individual modes more closely, we find that some have values below 0.9. For the mode near f = 19.51 d^{1} we find β_{4,1} ≈ 0.8 and a rotation rate of 0.16−0.17 d^{1}, while β_{3,1} ≈ 0.97 and hence Ω_{e,1} ≈ 0.135 d^{1}. This is because the n = 4 mode has a mixed mode character and not only probes the envelope (like pure pmodes) but also has a high amplitude near the convective core (like pure gmodes). This is also the case for the quadrupole mode, where we find envelope rotation rates of 0.14−0.15 d^{1}. To obtain a better estimate of the surfacetocore rotation rate, we only use the rotation rates derived from the pure pmode triplet near f = 12.76 d^{1}. For this triplet Ω_{e,1} ranges from 0.131 ± 0.003 d^{1} to 0.137 ± 0.003 d^{1}, using either modes of n = 1 or n = 2 of both Models 2 and 3. Thus, we find almost uniform rotation with a slightly faster core than envelope, with Ω_{e,1}/ Ω_{c,1} = 0.94 ± 0.2 up to 0.98 ± 0.2.
Moreover, we use the pmodes 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 pmodes of Model 6 have frequencies that are too high and therefore do not agree with the observations. Aside from the binary parameters, the pmode 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 Ftype g and pmode hybrid pulsators that reside in the binary system KIC 10080943 (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 2σ of the binary solution, together with asymptotic period spacing values within 2σ of the observed mean period spacing. With this approach, we provide the first such consistent seismic modelling of a binary Ftype gmode pulsator. We find that the stars have an age of ≃1.1 Gyr and have a low central hydrogen mass fraction of X_{c,1} ≃ 0.1 and X_{c,2} ≃ 0.2. This is a consequence of the observed mean period spacing values, which are below 3000 s and, thus, rather low for 2 M_{⊙} stars. In fact, our best models have masses 2σ 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 Ω_{c,1} = 0.13963 ± 0.00007 d^{1} and Ω_{c,2} = 0.09034 ± 0.00002 d^{1}. When the traditional approximation of rotation is used at these rotation rates instead of the firstorder 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 convectivecore 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 M_{1} = 1.82 M_{⊙}, M_{2} = 1.76 M_{⊙}, Z = 0.0125, f_{ov,1} = 0.008, f_{ov,2} = 0.005, log (D_{mix,1}) = 0.25, and log (D_{mix,2}) = 1.75; f_{ov} is expressed in local pressure scale height H_{P} and D_{mix} has a unit of cm^{2} s^{1}. For the step overshooting grid, the best model is similar at M_{1} = 1.81 M_{⊙}, M_{2} = 1.76 M_{⊙}, Z = 0.0125, α_{ov,1} = 0.11, α_{ov,2} = 0.05, log (D_{mix,1}) = 0.25, and log (D_{mix,2}) = 1.5.
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 X_{c} 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 pmode 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 + 3σ level. Especially, the primary radius of the model is 4σ 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 gmodes, whose probing power is situated near the core. We can estimate the tidal influence by the tidegenerating potential, which depends on the small parameter ε_{T} = (R_{1}/a)^{3}(M_{2}/M_{1}) (e.g. Willems & Aerts 2002) and find ε_{T} ≈ 3.7 × 10^{4} for KIC 10080943. The influence of the centrifugal force due to the rotation on the pulsations can be estimated by (Ω /f_{co})^{2}, 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 highprecision 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 ℓ = 1gmode, which has a frequency of exactly 15f_{orb}, and in the rotational splitting of the pmodes (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 convectivecore 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 pulsationconvection 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.
Acknowledgments
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.
References
 Aerts, C., Handler, G., Arentoft, T., et al. 2002, MNRAS, 333, L35 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Aerts, C., ChristensenDalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology, 1st edn. (Springer) [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Ballot, J., Lignières, F., Reese, D. R., & Rieutord, M. 2010, A&A, 518, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Bouabid, M.P., Dupret, M.A., Salmon, S., et al. 2013, MNRAS, 429, 2500 [NASA ADS] [CrossRef] [Google Scholar]
 Breger, M. 1990, Delta Scuti Star Newsletter, 2, 13 [NASA ADS] [Google Scholar]
 Cox, J. P., & Giuli, R. T. 1968, Principles of stellar structure(New York: Gordon and Breach) [Google Scholar]
 Dupret, M.A., Grigahcène, A., Garrido, R., Gabriel, M., & Scuflaire, R. 2005, A&A, 435, 927 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Guzik, J. A., Kaye, A. B., Bradley, P. A., Cox, A. N., & Neuforge, C. 2000, ApJ, 542, L57 [NASA ADS] [CrossRef] [Google Scholar]
 Hambleton, K. M., Kurtz, D. W., Prša, A., et al. 2013, MNRAS, 434, 925 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Herwig, F. 2000, A&A, 360, 952 [NASA ADS] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1993, ApJ, 412, 752 [NASA ADS] [CrossRef] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Keen, M. A., Bedding, T. R., Murphy, S. J., et al. 2015, MNRAS, 454, 1792 [NASA ADS] [CrossRef] [Google Scholar]
 Kurtz, D. W., Saio, H., Takata, M., et al. 2014, MNRAS, 444, 102 [NASA ADS] [CrossRef] [Google Scholar]
 Ledoux, P. 1951, ApJ, 114, 373 [NASA ADS] [CrossRef] [Google Scholar]
 Maceroni, C., Montalbán, J., Michel, E., et al. 2009, A&A, 508, 1375 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Maceroni, C., Lehmann, H., da Silva, R., et al. 2014, A&A, 563, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Miglio, A., Montalbán, J., Noels, A., & Eggenberger, P. 2008, MNRAS, 386, 1487 [NASA ADS] [CrossRef] [Google Scholar]
 Moravveji, E., Aerts, C., Papics, P. I., Andres Triana, S., & Vandoren, B. 2015, A&A, 580, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Moravveji, E., Townsend, R. H. D., Aerts, C., & Mathis, S. 2016, ApJ, 823, 130 [NASA ADS] [CrossRef] [Google Scholar]
 Pápics, P. I., Tkachenko, A., Aerts, C., et al. 2015, ApJ, 803, L25 [NASA ADS] [CrossRef] [Google Scholar]
 Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Paxton, B., Marchant, P., Schwab, J., et al. 2015, ApJS, 220, 15 [NASA ADS] [CrossRef] [Google Scholar]
 Perryman, M. A. C., de Boer, K. S., Gilmore, G., et al. 2001, A&A, 369, 339 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Saio, H., Kurtz, D. W., Takata, M., et al. 2015, MNRAS, 447, 3264 [NASA ADS] [CrossRef] [Google Scholar]
 Schmid, V. S., Tkachenko, A., Aerts, C., et al. 2015, A&A, 584, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Southworth, J. 2012, in Orbital Couples: Pas de Deux in the Solar System and the Milky Way, Workshop held October 2011, eds. F. Arenou, & D. Hestroffer, 51 [Google Scholar]
 Stellingwerf, R. F. 1979, ApJ, 227, 935 [NASA ADS] [CrossRef] [Google Scholar]
 Tassoul, M. 1980, ApJS, 43, 469 [NASA ADS] [CrossRef] [Google Scholar]
 Townsend, R. H. D. 2003, MNRAS, 343, 125 [NASA ADS] [CrossRef] [Google Scholar]
 Townsend, R. H. D., & Teitler, S. A. 2013, MNRAS, 435, 3406 [NASA ADS] [CrossRef] [Google Scholar]
 Van Reeth, T., Tkachenko, A., Aerts, C., et al. 2015a, A&A, 574, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Van Reeth, T., Tkachenko, A., Aerts, C., et al. 2015b, ApJS, 218, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Welsh, W. F., Orosz, J. A., Aerts, C., et al. 2011, ApJS, 197, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Willems, B., & Aerts, C. 2002, A&A, 384, 441 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
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 &star_job ! history and profile columns history_columns_file = ’hist.list’ profile_columns_file = ’prof.list’ !!! begin with a premain 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 = ’pp_cno_extras_o18_ne22.net’ change_initial_net = .true. ! display onscreen plots pgstar_flag = .false. / !end of star_job namelist &controls !!! 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 = 5d2 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.5d5 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) = 1d3 when_to_stop_rtol = 1d3 when_to_stop_atol = 1d3 !!! 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) = 1d2 xa_function_species(2) = ’he3’ xa_function_weight(2) = 100 xa_function_param(2) = 1d5 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
All Tables
Fundamental and seismic parameters of the primary and secondary component of KIC 10080943.
All Figures
Fig. 1 Observed gmodes in frequency (top) and period (bottom) as a function of degree ℓ are shown for the primary in red and the secondary in blue. The m values are distinguished by the different symbols. The yellow stars denote the ℓ = 2 modes, which could not be assigned to either component. 

Open with DEXTER  
In the text 
Fig. 2 Observed pmode 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 m values are distinguished by the different symbols. 

Open with DEXTER  
In the text 
Fig. 3 Top panel: HRD with the observed positions and 1σ 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 ΔP as a function of central hydrogen mass fraction X_{c}. The observed mean period spacing of the primary (secondary) is given as the horizontal solid (dashed) line and the 1σ areas are shown as the grey patch. 

Open with DEXTER  
In the text 
Fig. 4 Chisquared distribution of age for all coeval models in the two denser grids, which fall within 2σ of the HRD position of both stars and have equal composition. 

Open with DEXTER  
In the text 
Fig. 5 Chisquared distribution of central hydrogen mass fraction X_{c} 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. 

Open with DEXTER  
In the text 
Fig. 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 k are indicated near selected modes. Dipole modes of radial order −31 ≤ k ≤ −20 are shown for the primary and −37 ≤ k ≤ −25 for the secondary, while the radial orders of the quadrupole modes are −28 ≤ k ≤ −35 and −29 ≤ k ≤ −36, respectively. 

Open with DEXTER  
In the text 
Fig. 7 Comparison of the theoretical zonal pmodes (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. 

Open with DEXTER  
In the text 
Fig. 8 Difference between firstorder perturbation approximation and TAR for the rotation rates Ω_{c,1} = 0.13963 d^{1} and Ω_{c,2} = 0.09034 d^{1}. 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. 

Open with DEXTER  
In the text 
Fig. 9 Same as Fig. 6, but using the TAR instead of the Ledoux splitting to estimate the influence of rotation. 

Open with DEXTER  
In the text 
Fig. 10 Asymptotic period spacing ΔP as a function of central hydrogen mass fraction X_{c} for different masses and with different values of f_{ov} (top panel) and D_{mix} (bottom panel). The red points denote ages of 0.6 Gyr. 

Open with DEXTER  
In the text 
Fig. 11 Period spacing for different values of X_{c} along two evolutionary tracks incorporating exponential coreovershooting (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. 

Open with DEXTER  
In the text 
Fig. 12 Period spacing for different values of log (D_{mix}) along two evolutionary tracks, one with exponential coreovershooting (left panel) and one with step overshooting (right panel). The radial orders of the first and last mode in the period spacing is given. 

Open with DEXTER  
In the text 
Fig. 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. 

Open with DEXTER  
In the text 
Fig. 14 Same as Fig. 3 but for Model 6, indicated by stars and its evolutionary track as the black solid line. Here the 2σ areas for T_{eff,1,2} and log (L_{1,2}/L_{⊙}) are also shown as dashed boxes. 

Open with DEXTER  
In the text 