SeismicBased on data obtained from the CoRoT (Convection, Rotation and planetary Transits) space mission, developed by the French Space agency CNES in collaboration with the Science Programs of ESA, Austria, Belgium, Brazil, Germany and Spain. and spectroscopicBased on data obtained using the Télescope Bernard Lyot at Observatoire du Pic du Midi, CNRS and Université Paul Sabatier, France. characterization of the solar-like pulsating CoRoT target HD 49385

Seismicthanks: Based on data obtained from the CoRoT (Convection, Rotation and planetary Transits) space mission, developed by the French Space agency CNES in collaboration with the Science Programs of ESA, Austria, Belgium, Brazil, Germany and Spain. and spectroscopicthanks: Based on data obtained using the Télescope Bernard Lyot at Observatoire du Pic du Midi, CNRS and Université Paul Sabatier, France. characterization of the solar-like pulsating CoRoT target HD 49385

S. Deheuvels 1LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France
   H. Bruntt 1LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France
   E. Michel 1LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France
   C. Barban 1LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France
   G. Verner 2Astronomy Unit, Queen Mary, University of London Mile End Road, London E1 4NS, UK5 2    C. Régulo 3Instituto de Astrofísica de Canarias, 38205, La Laguna, Tenerife, Spain 34Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain 4    B. Mosser 1LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France
   S. Mathur 5Indian Institute of Astrophysics, Koramangala, Bangalore 560034, India 5    P. Gaulme 6Institut d’Astrophysique Spatiale, UMR8617, Université Paris XI, Bâtiment 121, 91405 Orsay Cedex, France 6    R. A. Garcia 7Laboratoire AIM, CEA/DSM-CNRS-Université Paris Diderot; CEA, IRFU, SAp, centre de Saclay, 91191, Gif-sur-Yvette, France 78GEPI, Observatoire de Paris, CNRS, Université Paris Diderot; 5 place Jules Janssen, 92190 Meudon, France 8    P. Boumier 6Institut d’Astrophysique Spatiale, UMR8617, Université Paris XI, Bâtiment 121, 91405 Orsay Cedex, France 6    T. Appourchaux 6Institut d’Astrophysique Spatiale, UMR8617, Université Paris XI, Bâtiment 121, 91405 Orsay Cedex, France 6    R. Samadi 1LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France
   C. Catala 1LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France
   F. Baudin 6Institut d’Astrophysique Spatiale, UMR8617, Université Paris XI, Bâtiment 121, 91405 Orsay Cedex, France 6    A. Baglin 1LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France
   M. Auvergne 1LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France
   I. W. Roxburgh 1LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France 2Astronomy Unit, Queen Mary, University of London Mile End Road, London E1 4NS, UK5 2
   F. Pérez Hernández 3Instituto de Astrofísica de Canarias, 38205, La Laguna, Tenerife, Spain 34Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain 4
Submitted …
Key Words.:
Methods: data analysis – Methods: statistical – Methods: observational – Stars: oscillations – Stars: individual: HD 49385
offprints: S. Deheuvels

Context:The star HD 49385 is the first G-type solar-like pulsator observed in the seismology field of the space telescope CoRoT. The satellite collected 137 days of high-precision photometric data on this star, confirming that it presents solar-like oscillations. HD 49385 was also observed in spectroscopy with the NARVAL spectrograph in January 2009.

Aims:Our goal is to characterize HD 49385 using both spectroscopic and seismic data.

Methods:The fundamental stellar parameters of HD 49385 are derived with the semi-automatic software VWA, and the projected rotational velocity is estimated by fitting synthetic profiles to isolated lines in the observed spectrum. A maximum likelihood estimation is used to determine the parameters of the observed p modes. We perform a global fit, in which modes are fitted simultaneously over nine radial orders, with degrees ranging from to (36 individual modes).

Results:Precise estimates of the atmospheric parameters (, [M/H], ) and of the  of HD 49385 are obtained. The seismic analysis of the star leads to a clear identification of the modes for degrees . Around the maximum of the signal (Hz), some peaks are found significant and compatible with the expected characteristics of modes. Our fit yields robust estimates of the frequencies, linewidths and amplitudes of the modes. We find amplitudes of ppm for radial modes at the maximum of the signal. The lifetimes of the modes range from one day (at high frequency) to a bit more than two days (at low frequency). Significant peaks are found outside the identified ridges and are fitted. They are attributed to mixed modes.


1 Introduction

In the Sun, oscillations are excited by the turbulent motions in the outer part of the external convective envelope and are further propagated into the interior of the star. The study of these oscillations has yielded constraints on the inner structure of the Sun, allowing us to estimate the sound speed and density profiles (2003ApJ...591..432B, 2001ApJ...555L..69T), the position of the base of the convective zone (1991ApJ...378..413C), and the rotation profile (2003ARA&A..41..599T, 2008A&A...484..517M). However, the very low amplitude of these oscillations (a few ppm in photometry) makes it very challenging to detect and analyze them in other stars than the Sun.

Achieving a better understanding of the interiors of solar-like pulsations is one of the main objectives of the space mission CoRoT (Convection, Rotation and planetary Transits). CoRoT is a space telescope performing high-precision photometry over quasi-uninterrupted long observing runs (Baglin06). Solar-like oscillations have already been studied in several other stars with CoRoT data. The analyses of these stars have encountered difficulties identifying the degrees of the modes, either because of a too low signal-to-noise ratio (HD 175726: 2009A&A...506...33M, HD 181906: 2009A&A...506...41G), or because of a too short lifetime of the modes, inducing large mode linewidths (HD 49933: 2008A&A...488..705A, HD 181420: HD181420 ).

The star HD 49385 is the first G-type solar-like pulsator observed in the seismology field of CoRoT. It is cooler than the solar-like pulsators previously analysed with CoRoT data, and probably more evolved (at the end of the Main Sequence or shortly after it). The choice of HD 49385 as a CoRoT target has motivated us to lead spectroscopic observations, performed with the NARVAL spectrograph at the Pic du Midi Observatory. The fundamental parameters of HD 49385 are derived from these observations, as described in Sect. 2. The photometric observations with CoRoT are presented in Sect. 3. Section 4 presents the study of the low-frequency part of the power spectrum, in search of a signature of the stellar rotation. The extraction of p-mode parameters is described in Sect. 5, and Sect. 6 is dedicated to conclusions.

2 Stellar parameters of HD 49385

The solar-like pulsator HD 49385 is a G0-type star with an apparent magnitude of (uvby catalog, 1998A&AS..129..431H). To determine the fundamental parameters of HD 49385 we analysed a high-quality spectrum from the NARVAL spectrograph mounted on the 2-m Bernard Lyot Telescope at the Pic du Midi Observatory. We stacked two spectra collected on 2009 January 10, totalizing an exposure time of 6 600 s. The spectrum was normalized by identifying continuum windows in a synthetic spectrum and fitting a low-order spline through these points. This was done order-by-order and we made sure the line depths agreed for the overlapping part of adjacent echelle orders. We measured a mean signal-to-noise (S/N) ratio of in the continuum for several line-free regions in the range  Å for a data sampling of data points per resolution element ().

2.1 Temperature, and metalicity

We used the semi-automatic VWA software (bruntt09) to fit synthetic profiles for more than 600 lines in the range from to  Å. We used a differential approach, meaning all abundances are measured relative to exactly the same lines in the solar spectral atlas from This differential approach is described in more detail by bruntt08. As part of the analysis, we determined the atmospheric model parameters , , and the microturbulence. This was done by requiring that abundances determined from Fe i lines do not correlate with the equivalent width (EW) or the excitation potential (EP). Furthermore we required that the same mean abundances are measured from the Fe i and Fe ii lines (”ionization balance”). We estimated the uncertainty on the model parameters by adjusting them until the correlations of Fe i with EW or EP became significant or the ionization balance deviated (see bruntt08 for details). We stress that the uncertainties are strictly internal errors because the underlying assumption is that the model atmosphere represents the star. Any systematic error in the temperature profile or departures from local thermal equilibrium (LTE) will likely affect the results. We therefore quadratically added 50 K and 0.05 dex to the uncertainty on  and . The spectroscopic parameters are listed in Table 1.

As an additional check of the surface gravity we fitted the wings of the pressure sensitive Ca lines at and  Å. This was done by first adjusting the van der Waals constants to fit the solar spectrum for the canonical value . The fitted values for the two lines for HD 49385 are and , which agrees very well with the Fe i/Fe ii ionization balance. The low value we obtain for the of HD 49385 indicates that it certainly is an evolved object.

  km s
Table 1: Spectroscopic parameters of HD 49385.

The parameters agree well with the calibration of the Strömgren indices. From TEMPLOGG (rogers95, we get  K, , and . TEMPLOGG determines a significant interstellar reddening . If we assume zero interstellar reddening, we get from the , using the 2005ApJ...626..465R calibration,  K and  K from . The uncertainties include the calibration uncertainty and the photometric uncertainty on the indices. Finally, 2004A&A...427..933K used line-depth ratios to determine effective temperatures of stars and found  K (internal error). They also determined  and Fe/H, but did not give uncertainties on these two parameters: and . In summary, our result from VWA agrees well with the photometric calibrations and the line-depth ratio method for .

We found a high abundance of lithium in HD 49385 of (relative to the Sun). This value appears to be typical for a dwarf star with solar metalicity (e.g., boesgaard05). We used the line list from ghezzi09, but did not include the relatively weak molecular CN bands; hence the Li abundance may be slightly overestimated. The Li feature around 6707.8 is known to be strongly affected by NLTE effects, but when using 3D instead of 1D atmosphere models these two effects appear to cancel to first order (asplund05). In other words, because we are using 1D LTE models, the imposed error should be relatively small. Furthermore, because we are calculating abundances relative to the same Li line feature as in the Sun, our estimate of [Li] should be robust as the atmospheric parameters of the two stars are quite similar.

The abundance pattern we determine is shown in Fig. 1 and listed in Table 2. The overall metalicity is determined as the mean of the metals with at least 10 lines used in the analysis (Si, Ca, Ti, V, Cr, Fe, Co, Ni) giving . This range is marked by the horizontal bar in Fig. 1.

Figure 1: Relative abundances of 22 elements measured in HD 49385. Open and solid circles are mean values for neutral and singly ionized lines, respectively. The yellow horizontal bar marks the mean metalicity with 1- uncertainty range.

2.2 Projected rotational velocity

Because we have a relatively high-resolution spectrum with very high S/N we were able to estimate the projected rotational velocity (). We fitted synthetic profiles to ten isolated lines in the range  Å. We calculated the synthetic profiles with the SYNTH software 1996A&AS..118..595V, thus taking into account weak blends. We considered three broadening mechanisms which were convolved with the synthetic profiles: (1) the instrumental resolution, (2) macroturbulence () and (3) . The instrumental power of resolution of NARVAL is , but to check this value we fitted six narrow telluric lines around 6290 Å. We found that a Gaussian function with a FWHM of 6290 Å/65 000 = 0.097 Å fitted these lines very well and we used the same value of for the stellar lines. For  and  we calculated a grid of convolved profiles for 0–8  km s in steps of 0.5  km s. To identify the best combination of  and  we compared the synthetic profile and the observed spectrum by computing the  value.

In Fig. 2 we show an example of the  contour for the Fe i line at 6156.1 Å. The plots for all ten lines are very similar and show the same strong correlation between  and . We marked three points on the contour in Fig. 2 with filled circles, and the synthetic profiles corresponding to these grid points are shown in Fig. 3. The observed spectrum is shown with open circles. From our analysis we can put a firm upper limit on  of 5  km s, corresponding to   km s. The values that provide the best fit based on the analysis of the  contours for all ten lines are   km s and   km s.

Figure 2: Example of the contour surface used to determine and for the Fe i line  Å. The three circles mark the minimum of the surface and the 1- uncertainty range. The corresponding three synthetic spectra are compared to the observed line in Fig. 3.
Figure 3: The upper panel shows the fit of three synthetic profiles to the observed line (marked by open circles). Each profile was convolved with different combinations of and as given in the top panel (cf. Fig. 2). The lower panel shows the difference between the observed and synthetic profiles, indicating a slightly better fit for the profile represented by a dashed line (the value for each fit is given in the lower panel).
El. El.
C i 9 V i 12
N i 1 Cr i 18
O i 4 Cr ii 7
Na i 4 Fe i 336
Mg i 4 Fe ii 27
Mg ii 2 Co i 12
Al i 4 Ni i 81
Si i 38 Zn i 2
Si ii 2 Sr i 1
S i 2 Sr ii 1
Ca i 13 Y ii 7
Ca ii 2 Ba ii 1
Sc ii 6 Ce ii 1
Ti i 30 Li i 1
Ti ii 9
Table 2: Abundances of 22 elements in HD 49385 relative to the Sun. The third column lists the number of lines used to determine the mean abundance listed in the second column. The quoted uncertainties are intrinsic RMS errors on the mean value.

2.3 Luminosity

Using the Hipparcos parallax of the star, mas (2007A&A...474..653V), we could derive the absolute visual magnitude of the object . The bolometric correction was determined by interpolating in the grid provided by 1998A&A...333..231B for appropriate values of temperature and for HD 49385. We obtained , which yields a luminosity . The values of and provided an estimate of the radius .

3 CoRoT observations

The star HD 49385 was one of the targets of the second Long Run in the asteroseismology field of the space mission CoRoT. 136.9 days of photometric data have been collected from October 2007 to March 2008, with a duty cycle of 88.2%. After being corrected from known instrumental effects as described in, the lightcurve consists of a series of equally-spaced measurements (32 s) in the heliocentric frame (also known as N2 data). Small gaps remain in the lightcurve (a few minutes long) mainly due to the passage of the CoRoT satellite in the South Atlantic Anomaly (SAA). A larger gap (3.5 days) was caused by an unexpected reset of the DPU (Digital Processing Unit), which was probably the consequence of the impact of an energetic particle, because this event occured while the satellite was crossing the SAA. These gaps were linearly interpolated, based on the measurements collected 1 000 s on either sides of the gap. A low-frequency trend due to the aging of the CCD was also corrected by removing a linear gain variation of day (see 2009A&A...506..411A for more details). The resulting lightcurve is shown in Fig. 4.

Figure 5 shows the power spectrum of the lightcurve, computed with the Fast Fourier Transform algorithm (FFT). It was normalized so that the integrated power spectrum from 0 to twice the Nyquist frequency corresponds to the variance of the time series. At low frequency the background rises due to stellar granulation and possibly to stellar activity (see discussion in Sect. 4). The signature of perturbations due to the orbit of the satellite (Hz) and the day (Hz) remains in the power spectrum. These perturbations were expected (see 2009A&A...506..411A) and are of relatively low amplitudes (a few tens of ppmHz for the strongest ones at low frequency). They appear as combinations of the harmonics of both frequencies. For instance, a peak appears at 150.1 Hz (corresponding to ). Finally, a broad excess of power, centered around 1 mHz, is clearly seen in the power spectrum and corresponds to the acoustic modes of HD 49385, which are studied in Sect. 5.

Figure 4: Lightcurve of 137 days of CoRoT observations on HD 49385, detrended using instrumental information.
Figure 5: Power spectrum of 137 days of CoRoT data for HD 49385. The peaks indicated by arrows are from left to right: the 14-Hz-peak discussed in Sect. 4, a peak at Hz, and a peak at Hz. The grey area shows the excess of power due to the p-mode oscillations.

4 Search for a signature of stellar rotation in the Fourier domain

The measurement of the rotation period from the lightcurve analysis was found to be useful in previous analyses of p-mode parameters (see 2008A&A...488..705A, HD181420, 2009A&A...506...33M). It indeed provides a direct measurement of the surface rotation period and hence a first estimate of the rotational splitting value.

Contrary to other solar-like pulsating CoRoT targets, no immediate signature of the activity can be seen in the lightcurve (see Fig. 4). Only a long period trend stands out. The same trend was also found in several stars observed on the same CCD during the same specific observation run. Its instrumental or environmental origin is thus established and is currently under study.

We first identified a peak around Hz (period of about 29 days), the significance of which is hard to assess because it presents a relatively low amplitude in a frequency domain where the intrumental noise can dominate the power spectrum. We performed a time-frequency analysis with Morlet wavelets (as described in morlet) which showed that this peak presents a certain stability in time. But because the periods we were testing correspond to a non-negligible fraction of the total observation period, the cone of influence (1998BAMS...79...61T) prevented us from drawing any conclusion.

A prominent peak can also be detected at a frequency of Hz (period of about 0.8 day) with a signal-to-noise ratio of more than 10. This peak is however unlikely to be directly related to rotation. Indeed, let us assume ad absurdum that the rotation frequency is Hz. Our estimate of the radius implies a surface velocity of more than 120  km s. To obtain the measured , we would need a value of the inclination angle lower than , which means that the star is seen pole-on. In this case it is extremely difficult to derive the signature of the activity from the light curve (Mosser09) and no peak related to the rotation period can be detected in the power spectrum, which contradicts our assumption. In the power spectrum of HD 49933, which is a target of the seismo-field for the same run, we also found a small peak around Hz, suggesting this peak might be instrumental. However, the frequencies do not match exactly, and the peak observed for HD 49933 does not appear to be significant. The origin of this peak remains unclear. But it cannot be related to rotation.

In summary, we found no clear signature of the stellar rotation of HD 49385. The low-frequency part of the power spectrum can therefore not be used to obtain a first estimate of the rotational splitting. We may however infer from the absence of a clear activity signal either that the star has a low activity, which cannot be detected even with the CoRoT sensitivity, or that the inclination of the star is close to pole-on.

5 Analysis of p mode oscillations

Figure 6 shows the power spectrum smoothed with a 10-Hz boxcar. The comb-like structure of this excess of power is typical of p-mode oscillations. This section presents the analysis of these modes.

Figure 6: Smoothed power spectrum of HD 49385 (with a 10 Hz boxcar) in the frequency range of the oscillations. The green curve represents the fitted background (granulation and white noise components). The blue dashed curve corresponds to a Gaussian fit of the raw power spectrum giving Hz (dotted line).

5.1 Frequency of maximum power and large spacing

To estimate the frequency of the maximum of signal, we fitted a Gaussian profile on the raw power spectrum in the frequency range of the oscillations and above a fitted background (fit described in Sect. 5.3). The result is shown in Fig. 6. We obtained as the frequency of maximum power Hz.

In Fig. 6 a comb-like structure is clearly visible between 600 and Hz. The autocorrelation of the power spectrum provides a first estimate of the mean value of the large spacing Hz.

5.2 Mode identification

We used the value of the large spacing given in the previous section to build the échelle diagram shown in Fig. 7. It is obtained by piling onto one another sections of 56 Hz of the power spectrum. Three clear ridges appear in the échelle diagram: two neighboring ones, referred to as ridges A and B in Fig. 7, and a third one (ridge C).

Figure 7: Échelle diagram of HD 49385 computed with a large spacing of Hz. The power spectrum is smoothed with a 1-Hz boxcar and normalized by the background (estimated in Sect. 5.3.1). Three clear ridges appear, labeled here as ridges A, B, and C. The white dashed lines correspond to the frequencies obtained from a fit to the data presented in Sect. 5.3, and summarized in Tables 3 and 4.

The pattern we observe in Fig. 7 is the one we expect to observe for low-degree high-radial-order modes in solar-like pulsators: two main ridges separated by about half a large spacing corresponding to and modes, and alongside the ridge a fainter one corresponding to modes. It then seems reasonable to identify the ridges as follows: neighboring ridges A and B correspond to and modes, respectively, and ridge C to modes. We note that 2009A&A...508..877M and bedding10 obtained the same identification for the degrees of the ridges using different approaches. The following analysis is based on this identification. However, since the previous analyses of CoRoT solar-like targets encountered ambiguity indentifying the modes (see HD181420, 2009A&A...506...41G), we decided to also consider the possibility of the alternate identification (for which ridge C corresponds to modes), which is discussed in Sect. 5.3.4.

Figure 7 also shows fainter peaks outside the identified ridges. We here study their significance by testing the H hypothesis in the frequency range of the oscillations. Because we searched for short-lived modes, we binned the power spectrum over bins, as prescribed in app04. The size of the boxcar ( bins Hz) for the binning was chosen to approximately match the expected linewidth for the p modes in solar-like pulsators. In this case, the noise is distributed as a with degrees of freedom. We then applied the H hypothesis, i.e. we supposed that the observations are due to pure noise. For each peak we computed the probability that pure noise is able to reach the height of the peak over the whole studied window. Low values of indicate a poor compatibility between the observations and noise, and the significance of the peaks is defined as . We show in Fig. 8 all the peaks with a significance higher than 95%.

We notice four peaks which appear with a significance above 98%, and lie outside the identified ridges. They are labeled as , , and on Fig. 8. The peak , clearly outside the pattern expected from asymptotic theory, could correspond to a mixed mode, given that HD 49385 is an evolved object. The peak might be an mode, assuming that the curvature of ridge B changes at low frequency. This hypothesis seems to be corroborated by the fact that such a change exists in ridge C. However, a more precise inspection shows that including in ridge B generates an abrupt step in the large separation, whereas the large separation varies quite smoothly. The peak could also be an mixed mode in avoided crossing. Indeed, at low frequency, the and ridges have different curvatures. This can be seen in Fig. 7 and will be confirmed by the large separation profiles derived from the fitted mode frequencies (see Sect. 5.3.2). It was shown in deheuvels10 that a low-degree avoided crossing creates a characteristic distortion in the ridge of degree . If we assume that is an mixed mode in avoided crossing, the expected curvature of the ridge is very similar to that of the observed ridge. We cannot establish this identification more firmly at this stage though. Because of the uncertainty regarding the identification of the mode, we preferred to consider and fit it individually (see Sect. 5.3). The peaks and are located at about 7 Hz on the left of ridge C. We studied the possibility for these peaks to correspond to mixed modes. In this case, they would have to be of different degree . Indeed, g modes of same degree and increasing radial order have frequencies increasingly close to each other. Two mixed modes of same degree spaced by would imply many more mixed modes in the ridge of degree at lower frequency. It is clear from ridges and that this is not the case. And if and were mixed modes of different degree, we would expect them to have different amplitudes, which does not seem to be the case. Their position in the échelle diagram on the left of ridge C suggests that they could also be modes. Several reasons led us to favor this identification. First, these two modes appear vertically aligned in the échelle diagram and around the maximum of the oscillations, i.e. at the place where it is most likely to observe modes if they are present. Then, and might be part of a fainter ridge (referred to as ridge D), because other peaks, which have a lower significance but seem to follow the alignment, appear in Fig. 7. Finally, the distance between ridges C and D is also consistent with this identification. Indeed, the asymptotic theory gives an expression for the average distance between the ridges


where the expression of is found in coursCD. By definition, the quantity corresponds to the distance between the and ridges, and corresponds to the distance between the and ridges. Figure 7 shows that we have Hz. We therefore expect modes to be approximately 8 Hz on the left of modes, which approximately matches the position of ridge D.

The significance is an efficient tool to select peaks, but it is not a guarantee for the detection of a signal, as stressed in app09. To confirm the existence of a signal, one needs to compute the posterior probability . The details of the calculation along with the assumptions we make are given in Appendix B. For and , assumed to be mixed modes, we found posterior probabilities of 0.001% and 8.7%, respectively. These peaks are therefore probably due to signal. The larger posterior probability found for is caused by the fact that the expected profile of mixed modes was built with very conservative priors . For and , assumed to correspond to modes, we obtained a very low posterior probability (1.4% and 1.1%, respectively), indicating that these peaks are highly likely due to signal, and a good compatibility between their observed profiles and the expected profile of modes for HD 49385. This confirms our identification for and . It is the first time that modes can be detected in a solar-like pulsator (other than the Sun). We explain in Sect. 5.3.2 why HD 49385 was a particularly favorable target to detect modes.

We remark that ridge D could also be seen as the component of the rotationally splitted mode. This hypothesis is studied in 5.3.4.

Figure 8: Échelle diagram of the peaks with a significance greater than 95% in the power spectrum smoothed over bins (Hz). Four peaks appear, which do not belong to the identified ridges A, B, or C, and are labeled as .

5.3 Extraction of p mode parameters

5.3.1 Method

Solar-like oscillations are the result of stochastic excitations of p modes in the convective envelope. They respond to the equation of a forced, damped and randomly excitated harmonic oscillator (see 1986ssds.proc..105D). Therefore their signature in the power spectrum can be modeled as Lorentzian profiles with a linewidth proportional to the inverse of the mode lifetimes, added to the background and perturbated by a multiplicative noise following a distribution with two degrees of freedom. We used a maximum likelihood estimation (MLE) to determine the most probable parameters of the observed p modes (see 1990ApJ...364..699A). This method was already applied to determine p-mode parameters for the Sun (Toutain92, Appourchaux98) as well as for CoRoT solar-like targets (2008A&A...488..705A, HD181420, 2009A&A...506...41G). As described in 1990ApJ...364..699A, the likelihood is defined as the product of the probability density of each frequency bin of the power spectrum


where is the observed power spectrum, the model described above, and the parameters of the Lorentzian profiles used in the model. The set of parameters that maximizes the likelihood corresponds to the most probable parameters. Maximizing is the same as minimizing the negative logarithmic likelihood , defined as


As mentioned above, perturbations remain in the power spectrum at frequencies corresponding to combinations of the harmonics of the orbital frequency and those of the day. Some of these perturbations can be seen in the frequency domain of the oscillations. We therefore gave a null weight to the corresponding bins so that our estimates of the p-mode parameters are not influenced by those peaks.

We estimated the background with a power law to describe the granulation component, as prescribed by 1985ESASP.235..199H. Prior to the fit of the p-mode parameters, we fitted a profile of the type


on the power spectrum, excluding the frequency domain of the p-mode oscillations to avoid any bias, with and as free parameters. We then used this background to fit the mode parameters.

The analysis of HD 49933, which was the first solar-like CoRoT target, showed that a global fit of all the modes gave more robust estimates of the mode parameters (see 2008A&A...488..705A). We therefore fitted modes simultaneously on a frequency range corresponding to nine overtones in the power spectrum. The assumptions made for the fit are the same as in 2008A&A...488..705A. We briefly recall them here. The splitting was assumed to be constant over the frequency range of the oscillations, and the components of non-radial modes were supposed to be symmetrical with respect to the component. We sliced the power spectrum in -wide intervals, in which the linewidth is assumed to be the same for all modes, and where only the height of the mode () is fitted. The heights of non-radial modes are determined from those of the mode. Assuming an equal repartition of the energy in the modes of different degrees , and taking the value of the limb darkening for this star into account, we obtain

where . Inside non-radial multiplets, the heights are given by the intensity visibilities (see 2003ApJ...589.1009G, 2006MNRAS.369.1281B).

5.3.2 Mode parameters estimates

Five teams worked independently on the extraction of the p-mode parameters and applied the method we described in Sect. 5.3.1 to fit the four ridges associated to modes. As was done in previous analyses of this type, we present in the following section one set of reference results, and we use results from the other teams to cross-check them and comment the robustness of the solution. The reference results of the fit are given in Tables 3 and 4.

(Hz) agr.() (ppmHz) agr.() (Hz) agr.() (ppm) agr.()
0 2.23
0 1.94
0 3.03
0 4.09
0 4.78
0 5.18
0 3.90
0 2.63
0 1.36
Table 3: Reference results for the parameters of the radial p modes. The heights given here are those of the modes (the heights of non-radial modes are computed using height ratios, see Sect. 5.3.1). The agreement (agr.) between the teams is given for each fitted parameter. Filled circles indicate that at least four of the five teams agreed on the value within 1 . indicates that more than one team disagrees with the value, and that the five teams agree within ().
Mode frequencies

There is a very good agreement between the different teams on the fitted mode eigenfrequencies. Only two frequencies in the ridge differ by more than 1 . Indeed, this ridge is too faint to offer reliable estimates of the mode frequencies away from the maximum of the signal. An échelle diagram of the fitted modes is given in Fig. 10. At the maximum of the signal, the obtained precision on the mode frequencies ranges from about 0.2 to about 0.6 Hz for modes and increases at the edges of the frequency domain of the oscillations.

At Hz, the mode is found merged together with the mode. This could be because the signal-to-noise ratio decreases when reaching the edges of the fitted domain. Indeed, with a lower signal-to-noise ratio, the analysis of previous CoRoT targets yielded several modes merged with the neighboring mode (see 2008A&A...488..705A, HD181420). However, the mode might also have been shifted because of an mixed mode in the neighborhood.

Figure 9 shows the large separation profiles derived from the fitted frequencies of the and modes. It appears, as mentioned in Sect. 5.2, that the and ridges have a different curvature at low frequency, which might be caused by a low-frequency avoided crossing.

Figure 9: Large separation profiles derived from the fitted frequencies of the (squares and blue line) and (triangles and red line) modes. We overplotted the 1- error bars. The cross corresponds to the large separation we would obtain, if were considered as an mode.
Figure 10: Échelle diagram of the fitted modes. The colors correspond to the following degrees: blue , red , green and purple . The surface of the disks is proportional to the heights of the modes. 3- error bars are overplotted.
Mode linewidths

Figure 11 shows the linewidths of the fitted modes. As was observed in the Sun, the linewidths of the modes increase with frequency and are almost constant around the maximum of the signal. They are larger than those of the Sun by a factor 2. The lifetimes of the modes can be deduced from the linewidth by the relation . For HD 49385 they range from about one day to two days. By comparison the mode lifetimes are shorter in other CoRoT solar-like pulsators, e.g. for HD 181420 (about six hours for the shortest ones, and 12 hours around the maximum of the signal).

It appears that one of the fitted intervals (Hz for the mode) has a linewidth somewhat larger than what could be expected considering the linewidths of neighboring orders (see Fig. 11). Indeed, the mode of this interval (Hz) is larger than the other modes in the neighborhood.

Figure 11: Linewidths of the fitted modes with 1- error bars (squares and full line). The linewidths for the alternate identification are also given (diamonds and dashed line).
Mode heights and amplitudes

As was mentioned before, we fixed the ratios of heights beween modes of different degrees in our fit. Because HD 49385 is the first CoRoT target for which the and ridges are clearly distinct in the échelle diagram, we had the opportunity to check the chosen height ratios, which are those usually adopted in this type of analysis. We performed a new fit releasing the height ratios and found

For modes, the theoretical value of the ratio lies within the 1- error bars of the fitted value. For and modes, the theoretical ratios are at 2 and 3 of the fitted values, respectively. At this stage, this difference is too small to be interpreted as an inconsistency with the theoretical height ratios. It would be interesting however to investigate this matter in future analyses of this type on different objects. We noticed that with free height ratios the results for mode frequencies, linewidths and amplitudes remained within the 1- error bars of the parameters obtained with fixed height ratios.

From the heights and linewidths of the modes we derived the amplitudes (in ppm), given by . To estimate the error bars on this quantity, we took into account the anticorrelation which exists between the determination of the height and the linewidth. The mode amplitudes are represented for modes as a function of frequency in Fig. 12. They range from 3.6 ppm to 5.4 ppm. However, these estimates correspond to instrumental values. Taking into account the instrumental response functions of CoRoT (as described in michel09), they can be converted into intrinsic bolometric amplitudes per radial mode by multiplying them by a factor 1.032. The obtained amplitudes are consistent with the prediction of the maximum amplitudes using scaling laws (, see Appendix A). The star HD 49385 is then quite specific among the other solar-like pulsators observed with CoRoT for which amplitudes were found significantly lower than the predicted ones (see 2008A&A...488..705A, HD181420, 2009A&A...506...41G, 2008Sci...322..558M).

The signal-to-noise (SNR) ratio of each mode computed from the fitted height and the background estimated in Sect. 5.3.1 is given in Table 3. The SNR is found to be very small for the fitted modes (about 0.25 at the maximum of the signal). However, when releasing the height ratios, we find . The SNR of the modes should be closer to 1 at the maximum of the signal.

Figure 12: Amplitudes of the fitted modes with 1- error bars (squares and full line). The amplitudes for the alternate identification also appear (diamonds and dashed line).
(Hz) agr.()
1 3.26
1 2.83
1 4.44
1 6.00
1 7.03
1 7.62
1 5.76
1 3.89
1 2.01
2 1.11
2 0.97
2 1.51
2 2.04
2 2.38
2 2.58
2 1.95
2 1.31
2 0.68
3 0.11
3 0.09
3 0.15
3 0.20
3 0.23
3 0.25
3 0.19
3 0.13
3 0.07
Table 4: Reference results for the parameters of the non-radial p modes (symbols are the same as in Table 3).
Splitting and inclination angle estimates

Most of the fits we performed converge toward a null value of the splitting and consequently an undefined value of the inclination angle. This can correspond to two different cases:

  • The star could be seen almost pole-on, i.e. with an inclination angle close to . In that case the components of the non-radial modes are not visible, and we therefore cannot estimate the splitting. However, the measured of the object and the estimate of the stellar radius allow us to give constraints on the splitting through the relation


    For a low inclination angle we expect higher values of the rotational splitting (e.g. Hz for lower than ).

  • We could also explain this result by a low value of the splitting compared to the linewidth of the modes, so that the fits are unable to separate the components of the non-radial modes. Equation 5 shows that a small rotational splitting corresponds to a high inclination angle.

This fitting procedure does not significantly favor any of these two hypotheses. However, they lead to almost identical p mode parameters (all the parameters agree within 1- error bars).

In both cases, the effect of the rotational splitting on the observed mode profiles is weak. Indeed, either the angle is small, and the components of the multiplets have very small amplitudes compared to the component, or the splitting is small, and the components are mixed together. We note that this favors the detection of lower degree modes and therefore justifies why modes seem to be detected in HD 49385, while they were not in previously analyzed CoRoT targets.

5.3.3 Possible mixed modes parameters estimates

The position of HD 49385 in the HR diagram suggests that it either reaches the end of the main sequence or lies shortly after it. The star HD 49385 is therefore evolved and its acoustic spectrum may contain mixed modes, i.e. modes which have a g-mode behavior in the deep interior and a p-mode behavior below the surface (see 1974A&A....36..107S). The mixed modes are known to experience the so-called avoided crossing phenomenon, during which their frequency does not follow the asymptotic approximation and they can be located out of the ridges (unno89). It was already noticed (e.g. 1991A&A...248L..11D) that the detection of these modes could yield a very valuable constraint on the size of the convective core.

As was noticed in Sect. 5.1, two peaks are detected in the power spectrum with a confidence level above 95% and are not located within any of the observed ridges in the échelle diagram (namely and ). They could be the signature of mixed modes. We note that if the peak , found at Hz is an mixed mode in avoided crossing, it could explain why the p mode at Hz seems not to follow the curvature of the ridge. We fitted these peaks as Lorentzian profiles with free frequencies, heights and linewidths. Because the mixed modes have a g-mode character near the center, their inertia is larger, and we expect them to have a longer lifetime and therefore a smaller linewidth than regular modes. We obtained the results given in Table 5. We notice that the fitted frequencies have small error bars, but the other parameters are ill-determined.

(Hz) (ppmHz) (Hz) (ppm)
Table 5: Results of a Lorentzian fit performed on the peaks identified as potential mixed modes.

5.3.4 Other mode identification scenarii

(Hz) agr.() (ppmHz) agr.() (Hz) agr.() (ppm) agr.()
0 2.76
0 2.85
0 4.28
0 5.26
0 5.47
0 6.88
0 4.95
0 3.38
0 1.60
Table 6: Reference results for the parameters of the radial p modes for the alternate identification (symbols are the same as in Table 3).
Can ridge D be the component of the modes?

Still considering the identification preferred in this paper (i.e. ridge C corresponding to modes), ridge D could also possibly be seen as the of the ridge, as already mentioned in Sect. 5.1. When not including the modes in the fit, we indeed find solutions with a splitting of about 7 Hz, corresponding to the distance between ridges C and D.

However, we ran models representative of HD 49385 and found that with this rotation frequency the components of non-radial multiplets are no longer symmetric with respect to the component. To treat this case correctly, we therefore need to take this asymmetry into account in our fits.

For this purpose, we studied this asymmetry on evolutionary models reproducing the and that we obtained for HD 49385 in Sect. 2.1 and with different angular velocities. The models were computed with the evolutionary code CESAM2k (1997A&AS..124..597M) and the eigenfrequencies were derived from them with the oscillation code FILOU (2008Ap&SS.316..155S). They showed that the effects of asymmetry should be non-negligible even for values of the splitting as low as a few Hz. We describe the asymmetry of the components of a mode of degree and order by the shift defined as


These shifts are represented in Fig. 13 for a model representative of HD 49385 for eigenmodes in the frequency range of the detected oscillations () and with a rotational splitting of 3 Hz. We observe that the shifts are expected to vary approximately linearly with respect to frequency. We therefore model the shift as


where corresponds to the order where the signal is maximal (). For modes in non-radial multiplets we now fit eigenfrequencies of the type


This adds two free parameters to the fit per fitted component.

The obtained results for the eigenfrequencies, heights, linewidths and amplitudes of modes of degrees all lie within the 1- error bars of those given in Table 3. The fit converges toward a splitting of Hz and an angle of . As expected, ridge D is fitted as the component of the modes. The obtained rotational shift is such that Hz per order , and Hz. However, using Eq. 5 and the fitted values of the splitting and the angle, we obtain  km s. This result is not compatible with the firm upper limit of 5  km s established for  in the spectroscopic analysis (see Sect. 2.2).

Figure 13: Frequency shift (defined in Eq. 6) of the components of non-radial multiplets for an evolutionary model representative of HD 49385 in the frequency range of the detected oscillations and with a splitting of 3 Hz. Crosses stand for modes and plus signs for modes. A linear regression of the variations of the shifts with respect to the frequency is overplotted (dotted lines).
Can ridge C correspond to modes?

As was said in Sect. 5.1, for the sake of completeness and in the light of previous analyses of p-mode oscillations in solar-like pulsators, we decided to also investigate the alternate identification. In this case, ridges C and D correspond to and modes, respectively. Ridge A could not possibly be an ridge, because it has much too strong amplitudes. Ridges A and B would then be the components of rotationally splitted modes.

We fitted the modes exactly the same way it was done before. The results are given in Tables 6 and 7, and an échelle diagram of the fitted modes is shown in Fig. 14. The linewidths, heights and amplitudes are very comparable with those obtained in the first identification. We obtain a splitting of Hz and an angle of , which is consistent with the idea that ridges A and B are indeed identified as modes because they are separated by about twice the fitted splitting. Using Eq. 5 as before, these fitted values of the splitting and angle give an estimate of  km s, in complete disagreement with the measured  km s.

The alternate identifications we considered lead to major disagreements with the measured  and were therefore rejected. This confirms the identification we adopted in Sect. 5.2 to perform the global fit of the p-mode parameters.

Figure 14: Échelle diagram for the fitted modes of the alternate identification. Same symbols as in Fig. 10.
(Hz) agr.()
1 4.22
1 4.39
1 6.56
1 8.05
1 8.36
1 10.49
1 7.53
1 5.15
1 2.44
2 1.36
2 1.43
2 2.13
2 2.62
2 2.73
2 3.43
2 2.47
2 1.69
2 0.80
Table 7: Reference results for the parameters of the non-radial p modes for the alternate identification (symbols are the same as in Table 3).

6 Conclusion

The star HD 49385 was characterized from both CoRoT seismic observations and NARVAL spectroscopic data. The atmospheric parameters of the star were derived and the 137-day-long photometric time series was analyzed.

A clear series of peaks associated to p modes was detected in the power spectrum around 1 mHz. Up to now, the CoRoT solar-like pulsators presented some ambiguities in the identification of ridges corresponding to different degrees . Here three very clear ridges appear in the échelle diagram, which were readily identified as , 1 and 2 modes. The ridge appears clearly distinct from the ridge in the échelle diagram. Furthermore, two peaks, part of a fainter ridge, were found to be significant and compatible with the characteristics we expect for modes in HD 49385. This probably constitutes the first photometric detection of modes in a solar-like pulsating star (other than the Sun).

We performed a global fit over nine radial orders with degrees ranging from to modes (simultaneous fit of 36 individual modes). We obtained precise estimates of the mode frequencies (uncertainty of about Hz for radial modes at the maximum of the signal). The obtained value of the maximum amplitude for modes ( ppm) is consistent with the estimate deduced from samadi07, unlike the other CoRoT solar-like pulsators (e.g. 2008A&A...488..705A).

We found no evidence of a rotational splitting of the modes. This can be explained either by a small inclination angle or by a low rotational velocity (inducing a small rotational splitting compared to the linewidth of the modes). In passing we stress that even for very slow rotators (rotational splitting as low as a few Hz) the components of non-radial multiplets are not expected to be symmetrical with respect to the component. We proposed a simple way of treating this assymetry in future analyses of solar-like pulsators.

The p modes of HD 49385 were found to have lifetimes ranging from about one day to two days, i.e. somewhat shorter than the mode lifetimes in the Sun, but significantly larger than those of the previously observed CoRoT solar-like targets. This explains why the spectrum of HD 49385 is much clearer than for previous CoRoT pulsators, for which the large mode linewidths made it harder to separate the ridge from the ridge. The results obtained for HD 49385 confirm that the linewidths of the modes in G-type pulsators are smaller than those of F-type pulsators, making their analysis easier.

The very high quality of the spectrum also enabled us to detect significant peaks outside the identified ridges. These peaks were found to be compatible with mixed modes, whose presence can be expected in the spectrum of evolved objects such as HD 49385. The existence of mixed modes in avoided crossing can explain some specific behaviors we observe in the low-frequency eigenmodes. In particular we found that the ridge is distorted compared to the ridge at low frequency. deheuvels10 showed that this type of distortion could be associated with a low-degree mixed mode in avoided crossing. The identification of the mode as an mixed mode would be consistent with the observed pattern. Other features like the mode found to overlap the mode (around 855 Hz) might also result from avoided crossing phenomena. This needs to be further investigated in the seismic interpretation of this star.

This work was supported by the Centre National d’Etudes Spatiales (CNES). We thank M.-A. Dupret for providing estimates of the mode linewidths for HD 49385, computed with the non-adiabatic pulstion code MAD. I. W. Roxburgh and G. Verner thank the UK Science and Technology Facilities Council for support under grant PP/E001793/1


Appendix A Estimate of the expected amplitudes for the oscillations in HD 49385

Scaling laws give an estimate of the amplitude of modes at the maximum of the signal. 1995A&A...293...87K established the following relation between the luminosity amplitude of the oscillation and the velocity amplitude :


We assume that the frequency at the maximum of the signal scales as the acoustic cutoff frequency , i.e.


samadi07 obtained


Inserting Eq. 10 and 11 into Eq. 9 we obtain


With the spectroscopic measure of the temperature K and our estimate Hz we obtain the expected amplitude for modes at the signal maximum ppm.

Appendix B Computation of the posterior probability

By applying the Bayes theorem and assuming an equiprobability of H and H, app09 obtained the following expression for the posterior probability of the H hypothesis given the observed data:


The expression of for a spectrum binned over bins is given in app04


where is the Gamma function. To solve Eq. 13, we need to define the alternate hypothesis H and to make prior assumptions on the expected signal.

Among the four studied peaks, two are assumed to be modes around the maximum of the signal ( and ). To obtain an expression for , we need an estimate of the expected height and linewidth of modes.

An estimate of the expected mode linewidths of HD 49385 was computed with the non-adiabatic pulsation code MAD combined with a 1D model reproducing the effective temperature, the luminosity and the of the object in the same way as done in app09. We obtain around the maximum of the signal a linewidth between 1 Hz and 2.5 Hz. We adopt a uniform prior for the linewidth, taking into account an uncertainty factor in the theoretical model: we assume that the maximum linewidth is twice as large as that given from the models ([Hz).

The height of the modes can be deduced from the amplitude of modes through the relation


for an angle of inclination . The visibility factors are computed from 2003ApJ...589.1009G. The expected amplitude of the modes around the maximum of the signal is computed in Appendix A. We impose a uniform prior for the amplitude of the modes ( [], where and are the theoretical values of at ) and for the angle of inclination ( []).

Based on app09 we obtain the following expression for the probability of the observed data given the signal


where is the normalization factor. and are given in app04 and are functions of the mode height and linewidth .

Inserting Eq. 14 and 16 into Eq. 13, we obtain for and a posterior probability of 1.4% and 1.1%, respectively.

Two peaks are assumed to be the signature of mixed modes ( and ). For this type of modes, a part of the energy is located in the g-mode cavity and the inertia is larger. We therefore expect a smaller linewidth for these modes. But we can have no clue of their degree . Because both their amplitudes and their linewidths greatly depend on , we use very conservative priors: uniform priors on the amplitude of the mode ( [] ppm) and on the linewidth ([] Hz).

We obtain posterior probabilities of 0.001% and 8.7% for and , respectively.

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