arxiv: v1 [astro-ph.sr] 19 May 2017

Transcription

1 Astronomy& Astrophysics manuscript no JNp c ESO 2017 May 22, 2017 Photometric and spectroscopic variability of the B5IIIe star HD L. Andrade 1, E. Janot-Pacheco 2, M. Emilio 1, 3, Y. Frémat 4, C. Neiner 5, E. Poretti 6, P. Mathias 7, M. Rainer 6, J.C. Suárez 8, 9, K. Uytterhoeven 10, M. Briquet 5, 11, P.D. Diago 12, J. Fabregat 13 and J. Gutiérrez-Soto 12 arxiv: v1 [astro-ph.sr] 19 May Observatório Astronômico/DEGEO, Universidade Estadual de Ponta Grossa, Av. Carlos Cavalcanti, 4748, , Ponta Grossa, PR, Brazil, laerteandrade@uepg.br 2 Universidade de São Paulo/IAG-USP, rua do Matão, Cidade Universitária, , São Paulo, SP, Brazil 3 Institute for Astronomy, University of Hawaii, Honolulu, HI 96822, USA 4 Royal Observatory of Belgium, 3 avenue circulaire, B1180 Brussels, Belgium 5 LESIA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, UPMC Univ. Paris 06, Univ. Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, Meudon, France 6 INAF Osservatorio Astronomico di Brera, via E. Bianchi 46, Merate (LC), Italy 7 Institut de Recherche en Astrophysique et Planétologie, CNRS, 14 avenue Edouard Belin, Université de Toulouse, UPS-OMP, IRAP, Toulouse, France 8 Universidad de Granada, Dept. Física Teórica y del Cosmos, 18071, Granada, Spain 9 Instituto de Astrofísica de Andalucía (CSIC), Glorieta de la Astronomía s/n, Granada, Spain 10 Instituto de Astrofísica de Canarias, E La Laguna, Tenerife, Spain 11 Institut d Astrophysique et de Géophysique, Université de Liège, Quartier Agora, Allée du 6 août 19C, 4000 Liège, Belgium 12 Depto. Didáctica de la Matemática, Universidad de Valencia, Avda. Tarongers, 4, Valencia, Spain 13 Observatorio Astronómico, Universidad de Valencia, Burjassot, Spain Received: June 23, 2016; accepted: March 24, 2017 ABSTRACT We analyzed the star HD , one of the relatively bright Be stars observed in the seismo field of the CoRoT satellite, in order to determine its physical and pulsation characteristics. Classical Be stars are main-sequence objects of mainly B-type, whose spectra show, or had shown at some epoch, Balmer lines in emission and an infrared excess. Both characteristics are attributed to an equatorially concentrated circumstellar disk fed by non-periodic mass-loss episodes (outbursts). Be stars often show nonradial pulsation gravity modes and, as more recently discovered, stochastically excited oscillations. Applying the CLEANEST algorithm to the highcadence and highly photometrically precise measurements of the HD light curve led us to perform an unprecedented detailed analysis of its nonradial pulsations. Tens of frequencies have been detected in the object compatible with nonradial g-modes. Additional high-resolution ground-based spectroscopic observations were obtained at La Silla (HARPS) and Haute Provence (SOPHIE) observatories during the month preceding CoRoT observations. Additional information was obtained from low-resolution spectra from the BeSS database. From spectral line fitting we determined physical parameters of the star, which is seen equator-on (i=90 ). We also found in the ground data the same frequencies as in CoRoT data. Additionally, we analyzed the circumstellar activity through the traditional method of violet to red emission Hα line variation. A quintuplet was identified at approximately c d 1 (12.88 µhz) with a separation of c d 1 that can be attributed to a pulsation degreel 2. The light curve shows six small- to mediumscale outbursts during the CoRoT observations. The intensity of the main frequencies varies after each outburst, suggesting a possible correlation between the nonradial pulsations regime and the feeding of the envelope. Key words. Stars: early-type Stars: emission-line, Be Stars: individual: HD Stars: oscillations Stars: rotation 1. Introduction Classical Be stars are main-sequence objects of mainly B-type, whose spectrum show, or has shown at some epoch, Balmer lines in emission and an infrared excess. Both characteristics are attributed to an equatorially concentrated circumstellar disk produced by sporadic mass ejections occurring during light brightening episodes, called outbursts. The ejection episodes can be explained by the rapid rotation of these objects, but The CoRoT space mission was developed and operated by the French space agency CNES, with participation of ESA s RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. This work is partially based on observations made with the 3.6-m telescope at La Silla Observatory under the ESO Large Programme LP185.D since most Be stars rotate below their critical rotation velocities (Frémat et al. 2005; Cranmer 2005; Yudin 2001) an additional explanation is needed. Because nonradial pulsations (NRP) in the Be starµcen (Rivinius et al. 1998) correlated with mass ejection episodes, NRP could provide the additional energy required to trigger the mass ejection episodes, converting a high Vsin i B star into a Be star. Huat et al. (2009) found for the hybrid pulsating Be HD 49330, from CoRoT and ground-based spectroscopic data, an increase in amplitude for g-modes and a decrease for p-modes just before an outburst and a reverse behavior after the outburst. Neiner et al. (2012) discovered stochasticallyexcited gravito-inertial modes in the Be star, HD 51452, enhanced by rapid rotation. They found that the mode amplitudes were linked to the presence of very-small-scale outbursts. This Article number, page 1 of 7

2 A&A proofs: manuscript no JNp and Huat et al. (2009) led them to propose a scenario to explain Be outbursts by the transport of angular momentum by waves (Neiner et al. 2013; Neiner & Mathis 2014). Despite the empirical evidence, the physical process responsible for the outbursts and mass loss in Be stars is still poorly understood. Furthermore, the effect of an outburst on NRP mode production by means of change in stellar structure is worth investigating (see Rivinius et al. 2013, for a comprehensive review on Be stars). Dziembowski et al. (1993) showed that high-order g-modes driven by the κ mechanism are unstable in the region of the HR diagram occupied by Be stars. The predicted radial velocity amplitude caused by anl=3g-mode is 20 km s 1 for a 4 M star for pulsation periods in the range d. For more massive and hotter stars such as HD , unstable modes are predicted with periods of d. Townsend (2005) found that mixed NRP modes (hybrid between Rossby and Poincaré modes) are unstable for B stars. Inertial-frame typical periods for Be stars occur in the range d. Pulsation analysis of a star requires the measurement of frequencies with both high photometric precision and high timefrequency resolution in addition to being free of one-day aliases, typical of ground-based observations. Only space missions provide both conditions. Stellar seismology made a great leap forward thanks to MOST, CoRoT, and Kepler satellites. The present paper reports the analysis of the Be star HD The observations with CoRoT and ground-based spectroscopy are described in Sect. 2. The star parameters, its light intensity behavior, and circumstellar variations are presented in Sect.3. We discuss the pulsation analysis in Sect. 4 and present our conclusions in Sect. 5. Fig. 1. CoRoT light curve of the Be star HD in the satellite reference frame (CoRoT Julian Date, CJD = HJD d). The colored regions (A) through (F) are divided by the occurrences of a medium and minor outbursts. 2. Observations 2.1. CoRoT photometric observations The main light curve of HD was obtained with the seismology core program of the CoRoT (Convection, Rotation and planetary Transits) space observatory (Baglin 2016), which studied bright pulsating stars (6<V< 9) with cadences of 32 and 512 seconds. For detailed information on the CoRoT satellite, see Auvergne et al. (2009) and Samadi et al. (2007). Data on HD were taken during d of the LRc06 long run from July 08th 2010 to September 24th 2010 with a cadence of 32 seconds. The light curve is shown in Fig. 1. It contains photometric observations. Five small-scale outbursts were observed in the light curve during the CoroT observations. For the purpose of the pulsation analysis we then divided the light curve into six regions as shown in Fig. 1 (see section 4.1.2) Ground-based spectroscopic observations Spectroscopic observations of HD were obtained with HARPS (3.6m telescope, ESO, La Silla, Chile) and SOPHIE (1.93m telescope, Haute Provence Observatory, OHP, France) spectrographs as a part of the ground-based campaign on the CoRoT targets (see Table 1). We also used low-resolution observations from the BeSS database (Neiner et al. 2011) from about the same epoch. HARPS was used in the high-efficiency EGGS mode R=80,000 with a fast readout (Mayor et al. 2003). From June 14 to July 2, 2010, 41 spectra were collected. The average exposure time was 1200s and the signal-to-noise ratio (S/N) ranges from 97 to 236. HARPS spectra are available in the SISMA Fig. 2. Average HARPS spectrum of HD centered around the Hα line (top panel) and the Hβ line (bottom panel). database (Rainer et al. 2016). Strong fringes in the reduced data were carefully corrected, as detailed in Neiner et al. (2012). Figure 2 shows the average HARPS spectrum around two of the Balmer lines. The Hα line profile depicts the status of emission activity during observations. SOPHIE was used in the high efficiency mode (HE) yielding a resolving power of R=40,000. Twelve spectra were obtained in late June, 2010 with an average exposure time of about 900s. The resulting S/N ratio varies from 90 to 213. Usual bias, flatfield and wavelength calibration was performed using the SOPHIE pipeline. We also retrieved reduced spectra in the Hα line region obtained by amateur astronomers in June and July 2010 from the BeSS database. The data were obtained with telescopes ranging from 21 to 30 cm in diameter, and a resolving power of R 15,000. Article number, page 2 of 7

3 Andrade et al.: Photometric and spectroscopic variability of the B5IIIe star HD Table 1. Spectroscopy observations of HD obtained close to the CoRoT photometric observations DATE Instrument # of Spectra Covered Time[hours] HARPS HARPS HARPS HARPS HARPS SOPHIE SOPHIE SOPHIE SOPHIE SOPHIE SOPHIE SOPHIE HARPS HARPS F/Fc Si II Si II Hδ C II Hγ Mg II F/Fc Si II Mg II Wavelength (Å) Fig. 3. Average SOPHIE blue spectrum of HD (black points). Red line is our best fit used to derive the stellar parameters (Table 2) and green line is the solution obtained previously with only one FEROS spectrum (Frémat et al. 2006). We highlight the sharp bottom shape of Balmer shell lines. Main spectral lines used to derive the apparent parameters are identified Wavelength (Å) Fig. 4. The same as Fig. 3 for some spectral regions. Predicted photospheric lines (red and green curves) are top flattened as expected from line transitions formed at the stellar equator. In the observations, they are dominated by the shell features. 3. The Be star HD Stellar parameters HD is a classical V= 7.65 mag Be star. First values of the fundamental parameters of the star were given in Frémat et al. (2006). This determination was based on one FEROS spectrum with S/N(5500 Å) 130. Apparent parameters, not considering gravitational darkening, were derived by fitting synthetic spectra with plane-parallel model atmospheres (1D spectra) to the FEROS data. A table of corrections was then used to correct the apparent parameters values for gravitational darkening and estimate the fundamental parameters of the non-rotating stellar counterpart as well as the "true" Vsin i. Currently, 53 good quality spectra obtained with the SOPHIE and HARPS spectrographs are available. We decided to rederive the apparent parameters by directly comparing 2D spectra that account for first-order rapidrotation effects (stellar flattening and gravitational darkening) to observations made of a combination of the best available spectra. In the procedure, the spectra were computed with the FASTROT routine, using the same atomic data and atmosphere models as in Frémat et al. (2006). The presence of quite stable, deep, narrow shell-like absorptions superimposed to the broader photospheric lines over the period of observations makes any fundamental parameter deter- mination very difficult. We therefore started using 1D spectra to fit the wings of the helium and hydrogen lines present in theλ wavelength range. Since theλ4026 Hei line did not seem to be affected by shell lines, we added the complete profile in the subsequent fits and rederived the apparent parameters and Vsin i which were found to be equivalent to those given in Frémat et al. (2006), as can be seen in the green lines in Figs. 3 and 4. In the next step, a grid of 2D spectra was computed for Ω/Ω c = 0.90, 0.95, and 0.99, i=70-90 o ( i=5 o, T eff = K ( T eff = 100 K), log(g)= ( log g=0.1). Lower inclinations i and rotation rates were not considered because they were not providing enough broadening to account for the extended helium line wings. Each synthetic spectrum was compared individually to the observations and the closest match was obtained for the parameters given in Table 2. The best spectrum is shown in red in Figs. 3 and 4 and compared to the 1D spectrum (green line). As the figures show, most shell lines are attributable to Feii and appear top flattened in the photospheric 2D synthetic spectrum, as we expect for transitions formed at the equator in the cooler parts of the star s surface. The determined stellar parameters are compatible with a rapidly rotating B5 III star. Article number, page 3 of 7

4 A&A proofs: manuscript no JNp Table 2. Stellar parameters of HD Superscript "nrcp" indicates the derivation for a non-rotating counterpart of the star (Frémat et al. 2005). Subscript "true" indicates the apparent Vsin i value corrected for gravitational darkening effects. Parameter Value T nrcp eff (K) log g nrcp 3.60 i 90 Vsin i true (km s 1 ) 345 Ω/Ω c 0.99 f rot (cd 1 ) R equ (R ) 8.18 R equ /R pole 1.46 T pole eff (K) T equ eff (K) 7632 log g pole 3.68 log g equ 2.30 V crit (km s 1 ) 356 V/R V+R HARPS SOPHIE BeSS HARPS SOPHIE BeSS We find that i = 90, that is, the star is seen equator-on. MoreoverΩ/Ω c is 0.99, and thus the star rotates close to critical rotation Circumstellar variations and occurrence of minor outbursts Be stars undergo aperiodic mass ejection episodes that feed their circumstellar disks. These so-called outbursts produce brightness variations and changes in the emission line profiles. Be stars are known to present short and long-term variations of the violet (V) and red (R) components of emission lines (V/R variations). The long-term variations (5 to 10 years) can be attributed to one-armed oscillations in the Be circumstellar envelope (Okazaki 1997). Short-term variations were first thought to be caused by rotating circumstellar structures but space-based observations give increasing evidence that they are caused by pulsations (Rivinius et al. 2013). There are also cases where the V/R variability was shown to be phase-locked to orbital period in binaries (Zharikov et al. 2013). Close examination of the HD CoRoT light curve shows the occurrence of a medium-sized outburst around CJD (HJD )= 3876 and four minor outbursts close to CJD = 3847, 3853, 3894 and 3905 (see Fig. 1). Some of these outbursts could actually be instrumental jumps in the light curve. Such jumps have been observed in other (non-be) stars. The only way to be sure that they are outbursts is to check light curves of other stars observed at the same time in the same CCD. If they also show a jump at the same time, then it is an instrumental jump. Only HD showed a jump in our verification, proving that it is in fact a real outburst. The relative intensity of pulsation frequencies found in the star change between the outbursts (see section 4.1.2). We measured the relative intensity of the violet to red components of the Hα emission line from HARPS and SO- PHIE observations, namely V/R = [V V c ]/[R R c ] (Carrier & Burki 2003). Three measurements from the BeSS database (Neiner et al. 2011) were also included (see Fig. 5). The V/R variations indicate two epochs of sudden change in V/R, centered around CJD 3818 and In both cases there is a definite increase in the V/R ratio and in the V+R relation, which can be attributed to minor outbursts. The time interval between CJD = HJD (d) Fig. 5. Violet over red (V/R) emission ratio (top panel) and sum of the V and R emission (bottom panel) measured in the emission peaks of the Hα line observed with HARPS, SOPHIE, and additional spectra from the BeSS database. The hashed regions indicate the occurrence of minor outbursts. these episodes is compatible with the outburst occurrence (10 to 20 d) seen in the CoRoT observations that began the following month. 4. Pulsation analysis Fourier analysis of photometric and spectroscopic data was performed with the CLEANEST algorithm (Foster 1995, 1996). It removes most of the alias patterns in time-series analysis as discussed previously in Emilio et al. (2010) Time-frequency analysis of CoRoT data Time-Frequency analysis of the whole light curve We assume a precision ν on the frequency determination equal to the inverse of the duration of the observing run, that is, ν=1/77.56 d=0.013 c d 1. Periods shorter than about one minute (ν=1440 c d 1 ) and longer than the total duration of the observing run of d (ν=0.013 c d 1 ) could not be detected. We applied the CLEANEST algorithm to look for pulsation frequencies on the HD light curve. In Fig. 6 we show the CLEANEST spectrum of the reconstructed time series. Table 3 shows a list of the main frequencies found. Frequencies that correspond to the star rotation period, c d 1 (compatible with the stellar parameters in Table 2 within our frequency precision) and a half of this period appear in the CLEANEST power spectrum. The ensemble of frequencies ( f 1 = c d 1, f 2 = c d 1, f 3 = c d 1, f 4 = c d 1 and f 5 = c d 1 ) seems to be a g-mode quintuplet. These frequencies are above the significance limit of 99.9%. They are equidistant with f c d 1. The distance between f 2 and f 3 is c d 1, Article number, page 4 of 7

5 Andrade et al.: Photometric and spectroscopic variability of the B5IIIe star HD Power Frequency (c/d) Fig. 6. Power spectrum of the frequencies extracted from the full CoRoT data (Fig. 1) using CLEANEST. Inset: Zoom of the region containing the most relevant powerful frequencies ( c/d). Fig. 7. Periodogram of the results of the time-frequency analysis of the CoRoT light curve (Fig. 1) of HD showing the variation in amplitude, including disappearance and reappearance, of the main frequencies during the CoRoT run. Table 3. The most powerful frequencies detected in CoRoT data with CLEANEST Frequency Remarks c d 1 µhz Power f rot / f 2 / f 4 / f rot f f f f f * f * f Fourier analysis between outbursts The five small-scale outbursts observed in the light curve suggest naturally a pulsation analysis in each time interval between them to search for potential correlations between the NRP regime and the outbursts, as observed in Huat et al. (2009). A moving window frequency analysis was made to investigate the frequency multiplet apparently found (Figure 7). We divide the photometric light curve into six regions (from A to F in Fig. 1) that is, between CJD=3847, 3853, 3876, 3894, and For each region, a CLEANEST spectrum was made (cf. Fig. 8). Although there is considerable loss of resolution because of the smaller time spans, some clear power variations can be seen. The dominant multiplet appears strong mainly in regions A, C, E, and F. NRP power reaches a minimum in region D, after the medium-sized outburst observed at the end of region C. The strength variation in time of the multiplet can also be seen in Fig. 7. It is possible that some of the pulsation energy was used to produce the outburst in HD Rivinius et al. (1998) and Huat et al. (2009) found correlations between variability in NRP intensity and the Be phenomenon for two other Be stars. More extensive and systematic work is needed to clearly establish such a correlation Ground-based data Fourier analysis of spectral line profiles compatible with c d 1 within the precision in our frequency determination from the light curve (0.013 c d 1 ). Deviations from equidistancy are expected for multiplets of fast rotating stars (Provost et al. 1981). Differential rotation can also alter the frequencies of the modes. First harmonics of some of these frequencies can be seen at about 2.2 c d 1. There is also a group of less regularly spaced, relatively strong frequencies around c d 1. As will be seen in the following section, each one of the quintuplet frequencies found in the whole time-frequency analysis appear with different power at different times of the photometric observations. We also applied the CLEANEST method to the time series of the spectroscopic data. Because of the scarcity of spectra, there was less information on NRP from line profile variations to be extracted (Neiner et al. 2012; Jankov et al. 2000). Nevertheless, peaks around the multiplet frequencies and their harmonics appear clearly. We fitted, by least squares, a Lorentzian function in the region of the multiplet, around 1.0 c d 1. The average half width at half-maximum (HWHM) found was ν = 0.026±0.009 c d 1 that we adopt as the frequency resolution of the spectroscopic data. Variations in those data can be visually inspected. We searched for phase variations in the Hβ spectral series formed by HARPS and SOPHIE data for the main frequencies found with CoRoT (Table 3). Article number, page 5 of 7

6 A&A proofs: manuscript no JNp Fig. 8. CLEANEST power spectrum calculated for each piece (A) through (F) in Fig. 1, similarly colored, in order to evaluate the effect of outbursts. The dashed lines indicate four of the most powerful frequencies in the full light curve (Fig. 6). Coherent phase variations were found for some of the quintuplet frequencies. The results are presented in Fig. 9. As pointed out above, nonradial pulsations can contribute to trigger the Be phenomenon but Rivinius et al. (2013) found that for a majority of Be stars, an ejection velocity of about 50 km s 1 is needed to form a Keplerian disk, even for very fast rotating stars. In order to look for possible velocity excursions caused by NRP, Fourier transform analysis has been applied to the Hei lines λ 4026, 4144 and 4471 in each high-resolution spectrum of HD The first minimum of the Fourier transform of line profiles can be reduced to velocity units so that it indicates the projected rotation velocity of the star (Gray 1976). However, the method works for symmetric, unblended lines, which is not the case of HD (see Figs. 3 and 4). Therefore, the Vsin i will be underestimated, as can be seen in Fig. 10, which in fact shows real line width variations. These variations can be interpreted as being caused by a changing NRP velocity field super-imposed to the supposedly constant rotation velocity field. Variations in the line width first minimum of HD occur nightly as well as during the same night with an amplitude of 15 km s 1. For a star seen almost equator-on, this is near the true amplitude of the velocity field and is not enough to trigger mass ejection episodes. The same kind of results were obtained for the Be starζ Oph by Jankov et al. (2000). Indeed, such a velocity amplitude is compatible to predictions by Dziembowski et al. (1993), who showed that g-modes of harmonic degreel 3 are unstable for early to late B stars. They predicted velocity amplitudes caused by NRP between 10 and 30 km s 1 for pulsation periods around 1 day. Thus, the velocity excursions detected in HD are probably insufficient to fully explain mass ejection episodes NRP and line-profile variations Telting & Schrijvers (1997) showed that pulsation quantum numbersl and m can be estimated in a series of spectra from the blue-to-red phase variations of the observed pulsation frequency and its first harmonics, respectively. The method is especially valid for rapid rotators as is the case of HD Again, our spectroscopic observations are limited, yielding a resolution in frequency of the order of just c d 1. This resolution hinders a precise analysis of the CoRoT frequencies. Nevertheless, variations in spectroscopic data can be inspected visually. We thus searched for phase variations in the Hβ spectral series formed by HARPS and SOPHIE data for the main frequencies found with CoRoT (Table 3). In Fig. 9 we plotted the amplitudes and phases of the variations of the multiplet frequencies as a function of the position in the Hβ Balmer line profile. They were calculated by CLEAN- EST along the mean SOPHIE + HARPS observations. Two cycles of pulsations were plotted shifted by multiples of π for clarity. For each frequency in Fig. 9 we read off the blue-to-red phase differences within the region where the variational amplitudes are sufficiently strong. The phase differences may be underestimated because of the low variational power in the line wings. Following Telting & Schrijvers (1997), the blue-to-red phase variations of frequencies that showed clear variations along the Hβ line were used to estimate the value of the NRP degreel. The value obtained for frequencies f 3 and f 4 yieldsl=1.7±0.2. A harmonic degreel=2 if rotationally split would produce a quintuplet (see Section 4.1.1) corresponding to a NRP order m with values -2, -1, 0,+1 and+2. Some harmonic frequencies of the quintuplet seem to be present around 2.2 c d 1. In principle, pulsation orders m could be determined from the phase variations across the Hβ line (Telting & Schrijvers 1997). Unfortunately, the low intensity of the harmonics prevent their phases to be estimated. 5. Conclusions We have detected twenty periodicities in the CoRoT light curve of HD compatible with nonradial pulsation g-modes (see Table 3). Indeed, they are expected to be present in a B5e star. Frequencies f 1 to f 5 are identified as a quintuplet and some of their harmonics may be present. High-resolution spectroscopic data were obtained during the month preceding the CoRoT observations; a detailed analysis of spectral data points out to a quite rapidly rotating B5 III star, seen equator-on. The phase variations of fundamental tones along the Hβ line are compatible with an azimuthal numberl=2±1. Four relatively small and a medium outbursts were observed in the light curve during the almost 80d time span of CoRoT observations. The frequency amplitudes vary neatly between outbursts. They pass through a minimum after the medium-sized outburst and this could be interpreted as pulsation power being temporarily transferred and consumed for pumping outbursts. It is also possible that the pulsations are less visible at the surface because they have been destabilized by the outbursts (see Neiner et al. 2013). Further detailed and extensive observations are still needed to establish a firm correlation between NRP regime and mass loss from Be stars. Acknowledgements. We wish to thank the CoRoT team for the acquisition and reduction of the CoRoT data. This work was supported in Brazil by FAPESP (project 2008/57866), CNPq (projects / and /2013-7) and CAPES (PNPD scholarship). EP acknowledges financial support from PRIN INAF MR acknowledges financial support from the FP7-Space project SpaceInn. JCS acknowledges funding support from Spanish public funds for research under project ESP C5-5-R (MINECO/FEDER), and under Research Fellowship program Ramón y Cajal (MINECO/FEDER). This work has made use of the BeSS database, operated at LESIA, Observatoire de Meudon, France: Article number, page 6 of 7

7 Andrade et al.: Photometric and spectroscopic variability of the B5IIIe star HD Fig. 9. Power (top panel) and phase (middle panel) variation of quintuplet frequencies along the mean HARPS and SOPHIE Hβ line profile (Fig. 2, bottom panel). Fig. 10. Position of the first (red) and second (blue) minima of Heiλ 4471 Fourier-transform profiles of HARPS and SOPHIE spectra. References Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411 Baglin, A., ed. 2016, The CoRoT Legacy Book: The adventure of the ultra high precision photometry from space, by the CoRot Team (EDP sciences) Carrier, F. & Burki, G. 2003, A&A, 401, 271 Cranmer, S. R. 2005, ApJ, 634, 585 Dziembowski, W. A., Moskalik, P., & Pamyatnykh, A. A. 1993, MNRAS, 265, 588 Emilio, M., Andrade, L., Janot-Pacheco, E., et al. 2010, A&A, 522, A43 Foster, G. 1995, AJ, 109, 1889 Foster, G. 1996, AJ, 111, 541 Frémat, Y., Neiner, C., Hubert, A.-M., et al. 2006, A&A, 451, 1053 Frémat, Y., Zorec, J., Hubert, A.-M., & Floquet, M. 2005, A&A, 440, 305 Gray, D. 1976, The observation and analysis of stellar photospheres, Wiley- Interscience publication (Wiley) Huat, A.-L., Hubert, A.-M., Baudin, F., et al. 2009, A&A, 506, 95 Jankov, S., Janot-Pacheco, E., & Leister, N. V. 2000, ApJ, 540, 535 Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20 Neiner, C., de Batz, B., Cochard, F., et al. 2011, AJ, 142, 149 Neiner, C., Floquet, M., Samadi, R., et al. 2012, A&A, 546, A47 Neiner, C. & Mathis, S. 2014, in IAU Symposium, Vol. 301, Precision Asteroseismology, ed. J. A. Guzik, W. J. Chaplin, G. Handler, & A. Pigulski, Neiner, C., Mathis, S., Saio, H., & Lee, U. 2013, in Astronomical Society of the Pacific Conference Series, Vol. 479, Progress in Physics of the Sun and Stars: A New Era in Helio- and Asteroseismology, ed. H. Shibahashi & A. E. Lynas-Gray, 319 Okazaki, A. T. 1997, A&A, 318, 548 Provost, J., Berthomieu, G., & Rocca, A. 1981, A&A, 94, 126 Rainer, M., Poretti, E., Mistò, A., et al. 2016, AJ, 152, 207 Rivinius, T., Baade, D., Stefl, S., et al. 1998, A&A, 333, 125 Rivinius, T., Carciofi, A. C., & Martayan, C. 2013, A&A Rev., 21, 69 Samadi, R., Fialho, F., Costa, J. E. S., et al. 2007, ArXiv Astrophysics e-prints [astro-ph/ ] Telting, J. H. & Schrijvers, C. 1997, A&A, 317, 723 Townsend, R. H. D. 2005, MNRAS, 364, 573 Yudin, R. V. 2001, A&A, 368, 912 Zharikov, S. V., Miroshnichenko, A. S., Pollmann, E., et al. 2013, A&A, 560, A30 Article number, page 7 of 7