The hybrid B-type pulsator \gamma Pegasi

The hybrid B-type pulsator Pegasi: mode identification and complex seismic modelling

P. Walczak, J. Daszyńska-Daszkiewicz, A.A. Pamyatnykh, T. Zdravkov
Instytut Astronomiczny, Uniwersytet Wrocławski, ul. Kopernika 11, 51-622 Wrocław, Poland
Nicolaus Copernicus Astronomical Center, Bartycka 18, 00-716, Warsaw, Poland
E-mail: walczak@astro.uni.wroc.plE-mail: daszynska@astro.uni.wroc.plE-mail:
Accepted … Received …; in original form …

We present interpretation of the oscillation spectrum of the early B-type star Pegasi, in which both low order p/g and high-order g-modes are observed. Using amplitudes and phases of the photometric and radial velocity variations, we identify/constrain the mode degree, , for all 14 detected frequencies. Seismic models fitting two pulsational frequencies corresponding to the modes , p and , g were constructed. This set of models contains those which reproduce also the empirical values of the complex nonadiabatic parameter associated to these two mode frequencies. Unfortunately, there are no models reproducing the values of for both frequencies simultaneously, regardless of model atmospheres, opacity data, chemical mixture as well as opacity enhancement in the and Deep Opacity Bumps. Most probably, some modifications of the opacities in stellar interiors are still required.

stars: early-type – stars: oscillations – stars: individual: Peg – atomic data: opacities.
pagerange: The hybrid B-type pulsator Pegasi: mode identification and complex seismic modellingLABEL:lastpagepubyear: 2002

1 Introduction

The simultaneous presence of both acoustic (p) and internal gravity (g) modes in a star, allows potentially for probing almost the entire interior. Such hybrid oscillations have been detected in a few main sequence stars of early B spectral type: Eridani (Handler et al., 2004; Aerts et al., 2004; Jerzykiewicz et al., 2005), 12 Lacertae (Handler et al., 2006; Desmet et al., 2009), Pegasi (Handler et al., 2009), HD50230 (Degroote et al., 2012), HD43317 (Pápics et al., 2012) and many candidates have been found (Pigulski& Pojmański 2008, Degroote et al. 2009, Balona et al. 2011). It was only possible thanks to multisite campaigns and space observations from the satellite missions like MOST (Walker et al., 2003), CoRoT (Baglin et al., 2006) and Kepler (Koch et al., 2010). The existence of both low order p/g modes as well as high-order g-modes makes the early B-type pulsators the attractive targets for asteroseismic studies.

The basic goal of asteroseismology is to find pulsational models with eigenfrequency reproducing the observational values. A more advanced seismic model should additionally account for instability and properties of oscillation modes. Daszyńska-Daszkiewicz, Dziembowski & Pamyatnykh (2003, 2005) introduced a new asteroseismic tool associated with each pulsational frequency. It is a ratio of the amplitude of the bolometric flux perturbations to the radial displacement, called the nonadiabatic -parameter. The theoretical values of are obtained from stellar pulsation computations and their empirical counterparts are derived from multicolour photometry and radial velocity data simultaneously with the mode degree, . The parallel fitting of pulsational frequencies and corresponding values of the -parameter was termed complex asteroseimology by Daszyńska-Daszkiewicz & Walczak (2009). The hybrid pulsators are of particular interest for such an in-depth modelling because the -parameter depends on pulsational frequency and the shape of eigenfunctions. Therefore its behavior is very different for the different frequency ranges. In the case of low order p/g modes (the high frequency range), the values of are independent of the mode degree and slowly vary with the frequency. On the contrary, for high order g-modes (the low frequency range), the values of are strongly -dependent and change rapidly with the frequency.

Complex seismic modeling has been already applied to the Cep star Oph (Daszyńska-Daszkiewicz & Walczak, 2009) and to the two Cep/SPB stars: Eri (Daszyńska-Daszkiewicz & Walczak, 2010) and 12 Lac (Daszyńska-Daszkiewicz, Szewczuk & Walczak 2013). In the case of Oph, we got a strong preference for the OPAL (Iglesias & Rogers, 1996) opacity tables. From the analysis of Eri a contradictory result was obtained: the Cep-type modes indicated the OPAL data whereas the SPB-type modes preferred the OP (Seaton, 2005) opacities. For 12 Lac, somewhat better agreement was found with the OP data. It can indicate, that the used opacity tables are not properly defined in the whole range of temperature and various physical conditions. Up to now, the sample of analyzed stars is too small to draw more general conclusions.

In this paper, we present mode identification and complex seismic modeling of Pegasi. The analysis of the space based observations from the MOST satellite (Handler et al., 2009) and ground based photometry and spectroscopy (Handler, 2009) led to the discovery of 14 pulsational frequencies with 8 typical for the Cephei stars and 6 typical for the Slowly Pulsating B-type stars (SPB). Because the star is a very slow rotator, here, we neglect all effects of rotation on pulsations.

In Section 2, we give a short description of the star. Section 3 contains results on mode identification for all detected pulsational frequencies using two approaches. Results of our seismic modeling are presented in Section 4. Conclusions are summarised in Section 5.

2 The hybrid pulsator: Pegasi

Peg (HR 39, HD 886) is a pulsating star of B2IV spectral type with the brightness of V=2.83 mag. The radial velocity variations were discovered by Burns (1911) but a very small amplitude made the variability uncertain. Series of spectrograms obtained by McNamara (1953) showed, that this star is a variable with the period of about and the amplitude of the radial velocity variations of 3.5 km/s. McNamara (1953) attributed Peg to the Cephei class of variable stars on the basis of its short period and spectral type. The follow-up spectroscopic and photometric observations confirmed changes of the radial velocity (McNamara, 1955, 1956; Sandberg & McNamara, 1960) and led to the discovery of the light variations in the yellow filter with the range of 0.015 mag (Williams, 1954). Jerzykiewicz (1970) determined the following light ranges in the UBV passbands: , and mag, respectively. Sareyan, Valtier & Le Contel (1975) determined period more precisely, d, and confirmed its constancy within 0.06 s per century. Smith & McCall (1978) and Cugier, Dziembowski & Pamyatnykh (1994) showed, that the variability of Peg can be explained by the radial pulsation.

In the catalogue of the galactic Cephei stars by Stankov & Handler (2005), Peg is the lowest mass variable, close to the Slowly Pulsating B stars (SPB) instability strip. Peg had been thought to be one of the monoperiodic Cep star until Chapellier et al. (2006) reported three additional pulsational frequencies; one in the Cephei range and two in the SPB frequency domain. Chapellier et al. (2006) claimed also that the star is a spectroscopic binary with an orbital period of 370.5 d. However, the analysis of data collected by the MOST satellite and from ground-based photometric and spectroscopic observations (Handler et al., 2009; Handler, 2009), showed that Peg is a single star and the hypothetical orbital variations can be explained by the high-order g-mode pulsation. Moreover, Handler et al. (2009) confirmed the frequencies of Chapellier et al. (2006) and discovered ten new ones. Consequently, we know that the light variations of Peg are caused by at least 14 pulsational frequencies: 8 of them are of the Cep type and 6 are of the SPB type. Thus, Peg is the hybrid pulsator in which low-order pressure/mixed modes and high-order gravity modes are excited simultaneously.

The rotational velocity of Peg was determined by Telting et al. (2006) who derived from the SiIII lines. Similar result was given by Pandey et al. (2011) and Nieva & Przybilla (2012), who got km/s and km/s, respectively. Handler et al. (2009) assuming that two pulsational frequencies are components of the , g triplet, determined km/s from the rotational splitting.

Determinations of the heavy elements abundance are ambiguous. Daszyńska-Daszkiewicz & Niemczura (2005) obtained the metal abundance of [m/H] (equivalent to the metallicity parameter of ) from the IUE ultraviolet spectra. From the optical spectra, Morel et al. (2006) and Pandey et al. (2011) derived and , respectively. The recent determinations from the optical and ultraviolet IUE spectra give (Wu et al., 2011) and (Nieva & Przybilla, 2012), while the analysis of the ultraviolet HST spectra indicates (Koleva & Vazdekis, 2012).

In Fig. 1 we show the observational error box of Peg in the HR diagram. We included the most recent determination of effective temperature by Wu et al. (2011), Koleva & Vazdekis (2012) and Nieva & Przybilla (2012). The luminosity was calculated from the Hipparcos parallax, mas (van Leeuwen, 2007), and the bolometric correction from Flower (1996). The total error box is as follows: and .

We depicted also the evolutionary tracks from the zero-age main sequence (ZAMS) to the terminal-age main sequence (TAMS) for masses , 8.5, 9.0 and 9.5, the hydrogen abundance of , metallicity of , two values of the overshooting parameter, and , the initial equatorial rotational velocity of km/s and the element mixture by Asplund et al. (2009), hereafter AGSS09. The tracks were computed by means of the Warsaw-New Jersey evolutionary code (e.g Pamyatnykh et al., 1998) adopting the OPAL opacities (Iglesias & Rogers, 1996). Lines labeled as and and models marked with a diamond and asterisk will be discussed later on.

Figure 1: The observational error boxes of Peg in the HR diagram. The evolutionary tracks were computed for metallicity of and two values of the overshooting parameter, . Lines of the constant period (0.15175 d) for the fundamental () and first overtone () radial mode are also drawn. Points labeled as a diamond and asterisk are discussed in Section 4.
c/d [mmag] [rad] [mmag] [rad] [mmag] [rad] [km/s] [rad]
=6.58974(2) 12.26(13) 2.669(11) 6.74(10) 2.628(15) 6.074(94) 2.607(15) 3.3582(69) 4.490(2)
=0.63551(10) 2.44(13) 0.870(57) 1.55(10) 0.734(67) 1.269(92) 0.947(76) 0.4992(70) 4.735(14)
=0.68241(7) 2.21(13) 5.607(61) 1.64(10) 5.681(63) 1.500(94) 5.644(63) 0.7364(72) 3.402(10)
=0.73940(10) 1.77(13) 3.542(75) 1.18(10) 3.604(89) 1.166(95) 3.782(83) 0.5214(71) 1.468(14)
=6.01616(14) 1.58(14) 4.592(87) 1.12(10) 4.664(92) 0.828(97) 4.54(12) 0.3612(85) 0.148(24)
=0.88550(7) 1.54(14) 3.733(87) 1.10(10) 4.014(94) 0.896(94) 4.00(11) 0.7236(68) 0.873(10)
=6.9776(5) 0.48(13) 3.08(28) 0.25(10) 3.27(42) 0.328(95) 3.36(29) 0.0962(72) 4.930(75)
=0.91442(11) 0.94(13) 4.45(15) 0.60(10) 4.80(17) 0.417(95) 4.95(23) 0.4652(68) 1.191(16)
=6.5150(8) 0.18(13) 3.46(73) 0.15(10) 2.44(68) 0.258(94) 2.64(37) 0.0628(70) 4.15(11)
=8.1861(8) 0.26(13) 4.92(50) 0.26(10) 4.89(40) 0.143(94) 5.57(66) 0.0639(69) 1.60(11)
=0.8352(3) 0.82(13) 1.30(16) 0.60(10) 0.97(17) 0.700(93) 1.43(14) 0.1785(67) 4.743(41)
=6.0273(5) 0.41(14) 1.55(34) 0.35(10) 0.53(30) 0.389(95) 0.57(25) 0.1104(82) 1.823(75)
=9.1092(12) 0.38(13) 3.61(34) 0.26(10) 3.38(39) 0.279(94) 3.24(34) 0.0408(68) 4.51(17)
=8.552(2) 0.20(13) 0.06(67) 0.37(10) 0.59(28) 0.222(94) 5.54(42) 0.0260(69) 1.70(27)
Table 1: Pulsational observables of Peg. In the first column, we give the pulsational frequencies of Peg. In the following columns there are values of the Strömgren , , amplitudes and phases. In the last two columns there is the radial velocity amplitude and corresponding phase.

3 Identification of oscillation modes

To identify the mode degree, , we made use of the light variations in the Strömgren passbands (Handler, 2009) and the radial velocity changes (Handler et al., 2009). In Table 1, we give amplitudes and phases of the light and radial velocity variations determined by fitting 14 frequencies detected in the MOST data. Therefore, our values of the amplitudes and phases can differe from those determined by Handler (2009).

We performed identification of the degree, , by applying two methods. In the first case, we compared theoretical and observational values of the amplitude ratios and phase differences between the available passbands and relied on the theoretical values of the -parameter. In the second approach, we used amplitudes and phases themselves and the -parameter was determined from the observations together with the mode degree, (Daszyńska-Daszkiewicz, Dziembowski & Pamyatnykh, 2003, 2005). In the case of B-type pulsators, the latter method demands the radial velocity measurements to get a unique identification of .

In both cases we need input from model atmospheres. Here, we present results obtained with the Kurucz (2004) models and the microturbulent velocity of km/s. Coefficients of the non-linear limb darkening law were adopted from Claret (2000). Identification of does not change if other value of or the non-LTE model atmospheres (Lanz & Hubeny, 2007) were used. The theoretical values of the -parameter are calculated with the nonadiabatic pulsational code of Dziembowski (1977).

In Table 2, we give the most probable values of the mode degrees for the Peg frequencies from the two approaches. The dominant mode, , is certainly radial. Also the identification of the , and is unique; they are modes with , 1 and 2, respectively. The frequencies , and can be either dipoles () or quadruples (). We got also two possibilities for the frequencies and . They can be or 2 and or 5 modes, respectively. Identification of the remaining frequencies are ambiguous. In the case of , , and we can say only that their mode degrees should be different from 4 and 6, and in the case of only degrees higher than 3 are possible. It is interesting that is close to the first overtone radial mode in many Peg models, but the degree was excluded by our photometric identification.

frequency phot. phot.+
c/d theoret. empir.
=6.58974 =0 =0 0
=0.63551 =1,2,4 =2,1,3,5 2,1
=0.68241 =1 =1 1
=0.73940 =1,2 =1,2 1,2
=6.01616 =1 =1,3 1
=0.88550 =2 =2 2
=6.9776 =?
=0.91442 4 =3,2,5 3,2
=6.5150 =?
=8.1861 =?
=0.8352 =1,2 =1,2 1,2
=6.0273 =? =2,5 2,5
=9.1092 =?
Table 2: Identification of from the two methods. In the first column, we give the pulsational frequencies of Peg, the second and third columns contain the mode degrees, , identified from photometric observables using the theoretical values of and from photometric observables and radial velocity data using the empirical values of , respectively. In the last column, we put values of consistent for both approaches.

In the HR diagram (Fig. 1), we plot lines of a constant period (0.15175 d) corresponding to for the fundamental and first overtone radial modes. As we can see the line is included only marginally in the error box. Nevertheless, two possibilities have to be considered. The discrimination of the radial order, , can be done by a comparison of the empirical and theoretical values of the -parameter (Daszyńska-Daszkiewicz, Dziembowski & Pamyatnykh 2005). In Fig. 2, we plotted this comparison on the complex plane assuming that the dominant frequency is the fundamental (the top panel) or the first overtone mode (the bottom panel). The theoretical values of were calculated for three different metallicities, , 0.015, and 0.020, overshooting parameter , hydrogen abundance of and the OPAL data. In all cases, we included models that are inside the observational error box of Peg. In the upper panel, we showed the line of the constant instabillity parameter , defined as (Stellingwerf 1978):


where is the work integral and is the stellar radius. Models located to the right of the line are unstable, whereas in the lower panel all models are stable. As we can see, agreement between the empirical and theoretical values of the -parameter can be achieved only if is the radial fundamental mode, . Moreover, constraints on the metal abundance were obtained, viz , 0.014). The conclusion about the radial order of did not change if the OP opacities were used. In this case, we got the metallicity in the range of . Moreover, identification of the radial order is independent of the value of the core overshooting parameter, . For the higher values of , the allowed metallicity range is shifted to the smaller value of .

Figure 2: Comparison of empirical and theoretical values of for the dominant frequency, , on the complex plane (, ) for models from inside the total error box. The OPAL data and AGSS09 mixture were adopted. The top and bottom panels correspond to the hypothesis of the fundamental and first overtone radial mode, respectively. In the upper panel, models to the right of the line are stable; in the lower panel all models are stable.

4 Complex asteroseismology

Basic seismic modeling consists in fitting pulsational frequency taking into account instability condition. This approach can be extended by a requirement of reproducing also the empirical values of the nonadiabatic parameter , which corresponds to each mode frequency.

In all computations, we used the Warsaw-New Jersey stellar evolution code (e.g. Pamyatnykh et al., 1998) and nonadiabatic pulsational code of Dziembowski (1977). Seismic models were calculated with both the OPAL and OP data. We assumed the initial rotational velocity of 3 km/s and two chemical mixtures: AGSS09 and that determined for Peg (Nieva & Przybilla, 2012), hereafter the Peg mixture. Due to the excitation problem of some modes, we tested also the opacity enhancement near the -bump and Deep Opacity Bump (DOB).

4.1 Fitting the centroid frequencies of Peg

Figure 3: The overshooting parameter, , as a function of metallicity, , for seismic models of Peg found from the fitting of the frequency (the , p mode) and the frequency (the , g mode) for the hydrogen abundance and OPAL opacities. We plot also the lines of constant mass, effective temperature and luminosity. Indicated vales of and are for the center and 1 error of the total error box. Models below the lines are unstable.

In the first step we searched for models fitting the value of the dominant frequency, , identified as the radial fundamental mode in previous section. Next, from this set of models we selected only those, which fit also the observational value of the frequency . In the whole set of models, we found that if is the radial fundamental mode, then can be only the dipole g mode. This identification of is unique, because the difference between frequencies of consecutive dipole modes is of the order of 1 c/d. Identification of the radial order is also independent of the azimuthal number, , since the rotational splitting is of the order of 0.005 c/d. The frequencies and were chosen because they are well identified low-order pressure/mixed modes. High-order gravity modes have a very dense frequency spectrum and fixing the radial order is problematic. The accuracy of the fitting and is of the order of c/d, which is equivalent to the observational errors.

These seismic models of Peg, calculated with the hydrogen abundance of , the AGSS09 chemical composition and OPAL opacities, are shown in Fig. 3 on the plane. We depicted the lines of constant mass, effective temperature and luminosity, as well as the instability borders for the radial mode, (thick solid line) and for the dipole mode, (thick dashed line). We indicated masses from 6.8 up to 9.8 and the values of and corresponding to the center and edges of the total error box. Models that are located below the lines are unstable. Almost all models between the lines and are inside the observational error box of Peg.

As we can see, for a given mass, models with higher metallicity require smaller values of the overshooting parameter, effective temperature and luminosity. The same results were obtained for Oph (Daszyńska-Daszkiewicz & Walczak 2009) and 12 Lac (Daszyńska-Daszkiewicz, Szewczuk & Walczak 2013). Lines of constant effective temperature and luminosity are nearly parallel to each other and more steep than lines of constant mass. Having such a huge number of seismic models fitting two frequencies of Peg, we determined an approximate relation between their parameters:


The very important point of the seismic analysis is mode instability. As can be seen in Fig. 3, all models with have the modes , p and , g stable. It is caused by a small abundance of heavy elements, which are crucial for exciting pulsations in the B-type stars. A lot of unstable models appear for higher value of metallicity, , and overshooting less than about for the radial mode and about for the dipole mode . Models with high overshooting are stable because they have too low masses and fall outside the instability region of the Cep stars (Pamyatnykh, 1999). For more thorough discussion of the mode instability analysis see Zdravkov et al. (2013).

Figure 4: The same as in Fig. 3, but we added lines of seismic models fitting the high-order g-mode frequencies which were identified as or modes. These are: , , , and . The line of models fitting the p-mode frequency (, p) is plotted as . The grey area indicates models located inside the total error box of Peg. In the case of two values of were considered: and .
Figure 5: The same as in Fig. 4, but models fitting frequencies identified as high-order g-modes with were added. These are: , , , and .

In Fig. 4 and 5, we plotted the same as in Fig. 3, but we marked also lines of models fitting some of seven frequencies: , , , , , and . In most cases, we assumed that these modes are axisymmetric () because there is no determination of . For we considered two possibilities, and , since there are some indications that the mode is prograde (Pandey et al., 2011). We considered simple rotational splitting in order to calculate the centroid frequency, that is , where is the Ledoux constant and is the rotational frequency.

Models located on each line fit simultaneously three frequencies. Fig. 4 includes frequencies which were identified as or 3, whereas Fig. 5 - those which have the degree . These seven frequencies were fitted with an accuracy of - c/d, equivalent to the observational errors.

In the case of high-order g-modes (, , , , , ), we usually had to include more than one radial order. In a considered range of metallicity and overshooting parameters, the frequency is the , g mode (in both cases, for and ) and the frequency is the , g or g mode. Identification of the mode degree of , , and is ambiguous and we had to consider more possibilities. For we got the modes , g and , g, g, g, for - the modes , g and , g. The frequency can be a mode with or . In the former case, it has to be the g or g mode and in the later - the g or g mode. If the frequency is a dipole, then it has to be g, but if it is a quadruple, it is g or g. Because in many models the pressure mode , p has a frequency close to , we plotted also a model line reproducing this frequency.

Models fitting three frequencies require higher overshooting, effective temperature, luminosity and mass if metallicity decrease. Unfortunately, high-order g-modes are usually stable. The modes (, g) and (, g) are exited if the metallicity parameter ; (, g), (, g) and (, g) if and (, g) and (, g) if .

Moreover, we were able to find a lot of models fitting four frequencies and two models fitting five frequencies. A model with , , and fits the modes (, p), (, g), (, g, ), (, g) and (, g). We will call it MODEL1 and it is marked with a diamond in Fig. 1. The parameters of the other model fitting five modes are , , and . The fitted modes are: (, p), (, g), (, g), (, g) and (, p). It is worth to note that these two models are located inside the observational error box of Peg. However, neither of these two models can reproduce the frequency , what can indicate that it is not an axisymmetric mode, as we have assumed.

In Fig. 6, we present the detailed comparison of the observational frequency spectrum and their possible theoretical counterparts for MODEL1 described above. We considered modes with from 0 to 3. Each theoretical frequency peak is labeled with the radial order, . The well identified frequencies have the best theoretical counterparts for mode degrees generally consistent with the photometric identification. The exception is . Despite the fact that it was identified as a dipole or quadruple mode, the best match is for the mode. A comparison with the theoretical spectrum may suggest, that this is the prograde mode (). The frequency can be the prograde mode, , m=-1, g as assumed by Handler et al. (2009). The frequencies and seems to be the , g and p modes, respectively. The frequencies and can be consecutive dipole modes, i.e., p and p. We can see also, that the frequency is consistent with the first overtone radial mode, but, as we mentioned in Section 3, it does not agree with our identification. The other close mode to is , g (not shown in Fig. 6), which is not inconsistent with photometric identification.

Figure 6: Comparison of the observed frequencies of Peg (upper panel) with the theoretical counterparts (bottom panel) corresponding to MODEL1 of Peg which has the following parameters: , , and . The height of lines in the upper panel corresponds to the amplitude in the Strömgren y passband while in the bottom panel to the mode degree, . The radial order is given at each theoretical frequency peak. The values of shown in the upper panel were derived in Section 3.

Let us now check the instability conditions for MODEL1. In the left panel of Fig. 7, we plotted the instability parameter, , as a function of frequency for modes with . For a comparison, in the right panel of Fig. 7, we plotted the same for a model with very similar parameters but computed with the OP data. Short, vertical lines represent the observed frequencies of Peg. Modes with are excited in the model. As we can see, with the OPAL opacities only frequencies in the range from about 5 up to 7 c/d are unstable. In the domain of the high order g-modes only high-degree modes () are excited with frequencies c/d. The frequencies in the range 8-9.5 c/d and 0.6-0.9 c/d are stable. Slightly better situation appears with the OP data. The model computed with these opacity data excite high-order g-modes with . But neither with the OP nor OPAL date we were able to excite modes with frequencies higher than about 7 c/d. It is tempting to explain this unsatisfactory results assuming higher opacities in the driving zones, but, as we will discuss it in the Section 4.3, it causes other problems.

Figure 7: Run of the instability parameter, , with frequency for MODEL1 (the left panel) and a model with very similar parameters but computed with the OP tables (the right panel). The vertical lines corresponds to the observed oscillation spectrum of Peg.

4.2 Fitting the nonadiabatic parameter,

As we have mentioned in the Introduction, an adequate seismic model should reproduce not only pulsational frequencies but also other characteristic parameters. Such parameter, which is determinable from observations, is the -parameter associated with each frequency. The theoretical values of depend mostly on properties of the driving zone.

To derive the empirical values of the -parameter, we used the LTE models of stellar atmospheres with the microturbulent velocity of km/s. The results do not differ qualitatively from those computed with the non-LTE atmosphere models (Lanz & Hubeny, 2007) as well as with higher microturbulent velocity km/s.

In Fig. 8, we put seismic models fitting, within the observational errors, the empirical values of the -parameter (hatched areas) for the radial fundamental mode (labeled as )) and for the dipole g mode (labeled as )). The area is larger than because of larger observational errors in the photometric amplitudes and phases of the frequency . The left and right panels of Fig. 8 correspond to computations with the OPAL and OP tables, respectively. As we can see, the requirement of fitting the -parameter for the mode reduce significantly the allowed range of stellar parameters of Peg. With the OPAL tables and the hydrogen abundance , only models with and can be considered. The OPAL models fitting the -parameter for the mode are inside the observational error box of Peg, while all models fitting the -parameter for the mode are outside the box. In the case of the OP models, the areas and part of are inside the error box.

Figure 8: The same as in Fig. 3 but we marked models fitting the empirical values of the nonadiabatic -parameter (hatched areas) of the radial fundamental mode (labeled as )) and of the dipole mode g (labeled as )). In the left panel we used the OPAL data and in the right panel the OP data.

Results for the OPAL and OP tables are similar and the biggest difference is in the position of the instability borders. For the OP models, these lines are shifted to the lower values of metallicity, , and the overshooting parameter, . Furthermore, for a given value of and , models with the OP tables have slightly higher mass whereas effective temperature and luminosity are almost the same. Moreover, the areas indicating models fitting the -parameter appear at the lower overshooting and metallicity. Unfortunately, there is no model fitting the -parameter for both and , simultaneously; neither with the OPAL nor OP data.

In the case of the most high-order g-modes, the empirical and theoretical values of the real part of agreed and the imaginary part differed significantly. This occurs both with the OPAL and OP opacities. Only for a mode we were able to find some seismic models which fit its empirical value of (both with the OPAL and OP data). These models, obtained with the OP data, have and , and are located inside the observational error box of Peg. However, they do not overlap neither with the area nor . The OPAL models reproducing the -parameter of are outside the allowed range of parameters and have and .

In the next step, we checked the effects of the hydrogen abundance and the heavy elements composition. In comparison with the results computed with the standard hydrogen abundance, , models with for a given value of and have larger masses (about ), smaller effective temperatures () and smaller luminosities (). Also models fitting the -parameters required much smaller value of the overshooting parameter. Still, there were no models fitting the -parameter for both and frequency simultaneously.

Then, we changed the chemical mixture. In Fig. 9, we show seismic models calculated with the chemical composition of Peg as determined by Nieva & Przybilla (2012). The Peg mixture has smaller abundance of iron-group elements than the Sun. These elements are extremely important for exciting pulsation and because of their deficiency, the instability region is smaller. Nevertheless, there are still a lot of unstable models fitting two pulsational frequencies and the -parameter corresponding to or . In this case seismic models fitting the values of the -parameter require slightly less efficient core overshooting.

Because the formal errors of can be underestimated, in Fig. 9, we marked additionally models fitting the -parameters within 2, (shaded-in, dark areas around the hatched regions) and 3 (shaded-in, bright areas around the hatched regions), where is the empirical error of . Now, the majority of models reproducing for fit also, within 3, the -parameter corresponding to .

Figure 9: The same as in the left panel of Fig. 8 but for the chemical composition of Peg as determined by Nieva & Przybilla (2012). The dark shaded-in and bright shaded-in areas around the hatched regions indicate models fitting the -parameters within 2 and 3, respectively. The asterisk marks model chosen for a comparison of the values of .

It is important to add, that in all considered cases, we were able to find also a model fitting, within 1 error, the - parameter for and the real part of for . The disagreement remains for the imaginary part of . This is clearly visible in Fig. 10, where we compare the theoretical and empirical values of in a function of the frequency. The left and right panels show the real () and imaginary () part, respectively. The parameters of the chosen seismic model marked with an asterisk in Fig. 1 and 9 are: , , , , metallicity , , the OPAL data and Peg mixture. We called it MODEL2. Modes with the degrees were considered.

We marked the empirical values of the -parameters for , and six high-order g-modes: ( and 1), (), ( and 1), (), ( and 3), ( and 2). The value of the imaginary part of the empirical -parameter for () is out of the scale. The empirical values of corresponding to different degrees have different symbols. We can conclude that only in the case of the real part of an agreement between the theoretical and empirical values of is quite good. The same result was obtained by Daszyńska-Daszkiewicz, Dziembowski & Pamyatnykh (2005) and Daszyńska-Daszkiewicz & Walczak (2010) for Eri.

The values of the -parameter for p modes are almost independent of the mode degree, , and change slowly with the mode frequency, while for high-order g-modes -parameter depend strongly on and changes rapidly with . The fact, that the empirical values of the real part of , especially for high-order g-modes, are located almost perfectly along the lines of theoretical counterparts are very encouraging for further studies and it shows a great potential of the -parameter, particularly for the SPB-type modes.

The empirical values of the nonadiabatic -parameter, plotted in Fig. 10, are listed in Table 3. In the fifth column we give also the values of the intrinsic amplitude, , multiplied by spherical harmonic, , where is the inclination angle. In the case of the radial mode, , we have the exact value of , whereas for other modes we can get only their lower limit. As we can see, the intrinsic amplitude of the dominant radial mode is very small and amounts to 0.26 per cent of the stellar radius and it is about three times larger than for other modes. The last column contains the values of the discriminant, , which measures goodness of the fit between the theoretical and observational values of the photometric amplitudes and phases.

6.58974 0 -9.43(11) 1.46(11) 0.002697(12) 4.50
0.63551 1 5.03(42) 5.54(41) 0.000584(11) 1.70
2 26.9(1.5) 20.8(1.5) 0.000376(6) 1.40
0.68241 1 3.81(15) 3.08(15) 0.000910(7) 0.54
0.7394 1 4.10(29) 2.84(29) 0.000684(10) 1.15
2 23.1(1.3) 10.4(1.3) 0.000451(8) 1.52
6.01616 1 -7.94(41) 1.99(40) 0.000381(7) 0.62
0.8855 2 8.87(37) 9.67(37) 0.000734(4) 0.37
0.91442 2 5.05(78) 7.54(77) 0.000485(6) 0.74
3 3.4(4.3) 56.7(4.2) 0.000717(7) 0.41
0.8352 1 2.31(92) 6.57(90) 0.000255(12) 1.50
2 16.3(3.7) 23.5(3.6) 0.000172(9) 1.77
Table 3: The empirical values of the real and imaginary part of , and the absolute value of the intrinsic mode amplitude for 8 frequencies of Peg. Columns from left to right are: pulsational frequency, mode degree, , empirical values of the real, , and imaginary, , part of the nonadiabatic parameter, the intrinsic amplitude multiplied by spherical harmonic, , and the discriminator, , describing the goodness of the fit. These quantities were derived for the parameters of MODEL2, marked with asterisk in Fig. 1 and 9.
Figure 10: The real (left panel) and imaginary (right panel) part of the nonadiabatic parameter as a function of pulsational frequency for MODEL2 of Peg marked as asterisk in Fig. 1 and Fig. 9. The empirical values of were put for the eight well identified frequencies; circles correspond to , triangles to , squares to and stars to .

4.3 Effect of the opacity enhancement

In the previous section, we have tried to find seismic models of Peg reproducing two frequencies and corresponding values of the -parameter. Using the standard opacity data and 1 error in we did not succeed in finding a model fitting the values of for and simultaneously. Because of this, we decided to change artificially the OPAL opacity tables. Two cases were examined.

In the first case we tested the 50% opacity enhancement near the -bump, i.e., at . The lines of constant effective temperature, luminosity and mass on the plane were nearly unchanged in comparison with the standard models, but opacity enhancement had a huge impact on the instability borders. Almost all our models had the modes , p and , g unstable for metallicity as low as . Unfortunately, with the modified opacities no model reproduces the empirical values of , neither for nor for .

As we have already mentioned, the theoretical values of the real part of calculated with the standard opacities (OP, OPAL) agreed very well with their empirical counterparts for almost all modes. The problem was to fit the imaginary part of . The change of the opacity near the -bump significantly altered the real part of , whereas the imaginary part was almost unaffected.

In the second case, opacities near the Deep Opacity Bump (DOB), occurring in the temperature range of , were increased by 20%. This modification did not change our models significantly. Also the instability regions for the modes and were nearly unaffected. This is caused by the fact that this bump is located deep inside the star and contains a small amount of mass. Therefore, this opacity modification did not substantially change the star structure. The only noticeable effect was the small reduction of the overshooting parameter of models fitting the -parameters for both and modes. However, there is still no seismic model fitting the -parameter for the modes and simultaneously. The effect of the opacity enhancement on the frequencies and their instabilities will be discussed in details in our next paper (Zdravkov et al., 2013).

5 Conclusions

The aim of this paper was to give a more detailed interpretation of the oscillation spectrum of the hybrid pulsator Pegasi. We began with identification of the mode degree, , for the 14 frequencies: 8 of the Cep type and 6 of the SPB type. Based on the two approaches, we were able to determine unambiguously four frequencies. For the other five frequencies two possible values of were obtained. In the case of the remaining five frequencies only some constraints were derived.

Then, we tried to construct seismic models which fit the two low order p/g-mode frequencies and their corresponding values of the nonadiabatic complex parameter . We chose the frequencies and , which were identified as the radial fundamental mode and dipole g mode, respectively. The problem we encountered was that there was no seismic model reproducing the -parameter for these two frequencies simultaneously. This inconsistence can be caused either by the underestimated errors or/and indicate that some additional effects should be included in pulsation modeling. One of the reasons could be inadequacies in the opacity data. Zdravkov & Pamyatnykh (2009) have suggested increasing opacity by 20-50% around the Z-bump and DOB (Deep Opacity Bump) to explain the observed frequency range of Peg. However, our studies showed that these artificially increased opacities spoiled even more the agreement between the empirical and theoretical values of the -parameter. With the modified opacities there is no model fitting the empirical values of , neither for nor .

Although we did not fully succeeded in constructing complex seismic models of Peg, we have showed directions and problems that need to be solved. There are the two main messages from this paper. The first one is a need for more accurate data on multi-colour time-series photometry and radial velocity data to better identify the observed modes and determine the empirical values of the -parameter. The second one is that if opacities are to ,,blame” for these disagreements and problems with mode instability, the improvement in computations of this microphysics data should be done in a more sophisticated way.


We gratefully thank Gerald Handler for kindly providing data on photometric and radial velocity variations. PW’s work was supported by the Human Capital Programme grant financed by the European Social Fund. AAP and TZ acknowledge partial financial support from the Polish NCN grants 2011/01/B/ST9/05448 and 2011/01/M/ST9/05914.


  • Aerts et al. (2004) Aerts C., de Cat P., Handler G., et al.: 2004, MNRAS 347, 463
  • Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A.J., Scott, P.: 2009, ARA&A 47, 481
  • Baglin et al. (2006) Baglin, A., Michel, E., Auvergne, M., The COROT Team 2006, in Proc. SOHO 18/GONG 2006/HELAS I, Beyond the spherical Sun, ESA SP, Sheffield, 624
  • Balona et al. (2011) Balona, L.A., Pigulski, A., Cat, P., et al.: 2011, MNRAS, 413, 2403
  • Burns (1911) Burns, K.: 1911, LicOB 6, 141
  • Chapellier et al. (2006) Chapellier, E., Le Contel, D., Le Contel, J.M., Mathias, P., Valtier, J.-C.: 2006, A&A 448, 697
  • Claret (2000) Claret, A., 2000, A&A, 363, 1081
  • Cugier, Dziembowski & Pamyatnykh (1994) Cugier, H., Dziembowski, W.A., Pamyatnykh, A.A.: 1994, A&A, 291, 143
  • Daszyńska-Daszkiewicz, Dziembowski & Pamyatnykh (2003) Daszyńska-Daszkiewicz J., Dziembowski, W.A., Pamyatnykh, A.A.: 2003, A&A 407, 999
  • Daszyńska-Daszkiewicz, Dziembowski & Pamyatnykh (2005) Daszyńska-Daszkiewicz J., Dziembowski, W.A., Pamyatnykh, A.A.: 2005, A&A 441, 641
  • Daszyńska-Daszkiewicz & Niemczura (2005) Daszyńska-Daszkiewicz, J., Niemczura, E.: 2005, A&A, 433, 431
  • Daszyńska-Daszkiewicz & Walczak (2009) Daszyńska-Daszkiewicz, J., Walczak, P.: 2009, MNRAS 398, 1961
  • Daszyńska-Daszkiewicz & Walczak (2010) Daszyńska-Daszkiewicz, J., Walczak, P.: 2010, MNRAS 403, 496
  • Daszyńska-Daszkiewicz, Szewczuk & Walczak (2013) Daszyńska-Daszkiewicz, J., Szewczuk, W., Walczak, P.: 2013, in preparation
  • Degroote et al. (2009) Degroote, P., Briquet, M., Catala, C., et al.: 2009, A&A, 506, 111
  • Degroote et al. (2012) Degroote, P., Aerts, C., Michel, E., et al.: 2012, A&A, 542, 88
  • Desmet et al. (2009) Desmet, M., Briquet, M., Thoul, A. et al.: 2009, MNRAS 396, 1460
  • Dziembowski (1977) Dziembowski, W.A.: 1977, AA 27, 203
  • Flower (1996) Flower P. J., ApJ 469, 355
  • Handler (2009) Handler, G.: 2009, MNRAS 398, 1339
  • Handler et al. (2004) Handler, G., Shobbrook, R. R., Jerzykiewicz, M. et al.: 2004, MNRAS 347, 454
  • Handler et al. (2006) Handler, G., Jerzykiewicz, M., Rodríguez, E. et al.: 2006, MNRAS 365, 327
  • Handler et al. (2009) Handler, G., Matthews, J.M., Eaton, J.A. et al.: 2009, ApJ 698, 56
  • Huang & Gies (2008) Huang, W., Gies, D.R.: 2008, ApJ 683, 1045
  • Iglesias & Rogers (1996) Iglesias, C.A., Rogers, F.J.: 1996, ApJ 464, 943
  • Jerzykiewicz (1970) Jerzykiewicz, M.: 1970, AcA 20, 93
  • Jerzykiewicz et al. (2005) Jerzykiewicz, M., Handler, G., Shobbrook, R.R., Pigulski, A., Medupe, R., Mokgwetsi, T. Tlhagwane, P., Rodrííguez, E.: 2005, MNRAS 360, 619
  • Koch et al. (2010) Koch, D.G., Borucki, W.J., Basri, G.: 2010, ApJ 713, 79
  • Koleva & Vazdekis (2012) Koleva, M., Vazdekis, A.: 2012, A&A 538, 143
  • Kurucz (2004) Kurucz R.L.: 2004,
  • Lanz & Hubeny (2007) Lanz T., Hubeny I., 2007, ApJS, 169, 83
  • McNamara (1953) McNamara, D.H.: 1953, PASP 65, 144
  • McNamara (1955) McNamara, D.H.: 1955, ApJ 122, 95
  • McNamara (1956) McNamara, D.H.: 1956, PASP 68, 158
  • Morel et al. (2006) Morel, T., Butler, K., Aerts, C., Neiner, C., Briquet, M.: 2006, A&A 457, 651
  • Nieva & Przybilla (2012) Nieva, M.-F., Przybilla, N.: 2012A&A, 539, 143
  • Pamyatnykh (1999) Pamyatnykh, A.A.: 1999, AcA, 49, 119
  • Pamyatnykh et al. (1998) Pamyatnykh, A.A., Dziembowski, W.A., Handler, G., Pikall, H.: 1998, A&A, 333, 141
  • Pandey et al. (2011) Pandey, C.P., Morel, T., Briquet, M., Jayakumar, K., Bisht, S., Sanwal, B.B.: 2011, CoAst 162, 21
  • Pápics et al. (2012) Pápics, P.I., Briquet, M., Baglin, A., et al.: 2012, A&A, 542, 55
  • Pigulski & Pojmański (2008) Pigulski, A., Pojmański, G.: 2008, A&A, 477, 917
  • Seaton (2005) Seaton, M.J.: 2005, MNRAS 362, 1
  • Sandberg & McNamara (1960) Sandberg, H.E., McNamara, D.H.: 1955, PASP 72, 508
  • Sareyan, Valtier & Le Contel (1975) Sareyan, J.-P., Valtier, J.-C., Le Contel, J. M.: 1975, A&A 44, 215
  • Smith & McCall (1978) Smith, M. L., McCall, M. L.: 1978, ApJ 221, 861
  • Stankov & Handler (2005) Stankov, A., Handler, G.: 2005, yCat, 21580193
  • Stellingwerf (1978) Stellingwerf, R. F.: 1978, AJ 83, 1184
  • Telting et al. (2006) Telting, J.H., Schrijvers, C., Ilyin, I.V., Uytterhoeven, K., De Rider, J., Aerts, C., Henrichs, H.F.: 2006, A&A 452, 945
  • Valtier et al. (1985) Valtier, J.-C., Chapellier, E., Morel, P.-J., Le Contel, J.M.: 1985, IBVS 2843
  • van Leeuwen (2007) van Leeuwen, F.: 2007, A&A 474, 653
  • Walker et al. (2003) Walker, G., Matthews, J., Kuschnig, R. at al.: 2003, PASP 115, 1023
  • Williams (1954) Williams, A.D.: 1954, PASP 66, 25
  • Wu et al. (2011) Wu, Y., Singh, H.P., Prugniel, P., Gupta, R., Koleva, M.: 2011, A&A 525, 71
  • Zdravkov & Pamyatnykh (2009) Zdravkov, T., Pamyatnykh, A.A.: 2009, in Stellar Pulsation: Challanges for Theory and Observation, AIP Conf. Proc., eds. J. A. Guzik & P. A. Bradley, Vol. 1170, 388
  • Zdravkov et al. (2013) Zdravkov, T., Pamyatnykh, A.A., Daszyńska-Daszkiewicz, J., Walczak, P.: 2013, in preparation
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