Seismic and spectroscopic characterization of the solarlike pulsating CoRoT target HD 49385^{}^{,}
S. Deheuvels^{1}  H. Bruntt^{1}  E. Michel^{1}  C. Barban^{1}  G. Verner^{2}  C. Régulo^{3,4}  B. Mosser^{1}  S. Mathur^{5}  P. Gaulme^{6}  R. A. Garcia^{7,8}  P. Boumier^{6}  T. Appourchaux^{6}  R. Samadi^{1}  C. Catala^{1}  F. Baudin^{6}  A. Baglin^{1}  M. Auvergne^{1}  I. W. Roxburgh^{1,2}  F. Pérez Hernández^{3,4}
1  LESIA, UMR8109, Université Pierre et Marie
Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon
Cedex, France
2  Astronomy Unit, Queen Mary, University of London Mile End Road, London E1 4NS, UK
3  Instituto de Astrofísica de Canarias, 38205, La Laguna, Tenerife, Spain
4  Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain
5  Indian Institute of Astrophysics, Koramangala, Bangalore 560034, India
6  Institut d'Astrophysique Spatiale, UMR8617, Université Paris XI, Bâtiment 121, 91405 Orsay Cedex, France
7  Laboratoire AIM, CEA/DSMCNRSUniversité Paris Diderot; CEA, IRFU, SAp, centre de Saclay, 91191, GifsurYvette, France
8  GEPI, Observatoire de Paris, CNRS, Université Paris Diderot; 5 place Jules Janssen, 92190 Meudon, France
Received 16 October 2009 / Accepted 8 March 2010
Abstract
Context. The star HD 49385 is the first Gtype
solarlike pulsator observed in the seismology field of the space
telescope CoRoT. The satellite collected 137 days of highprecision
photometric data on this star, confirming that it presents solarlike
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 semiautomatic 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.
Key words: methods: data analysis  methods: statistical  methods: observational  stars: oscillations  stars: individual: HD 49385
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 (Basu et al. 2003; TurckChièze et al. 2001), the position of the base of the convective zone (ChristensenDalsgaard et al. 1991), and the rotation profile (Thompson et al. 2003; Mathur et al. 2008). 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 solarlike pulsations is one of the main objectives of the space mission CoRoT (Convection, Rotation and planetary Transits). CoRoT is a space telescope performing highprecision photometry over quasiuninterrupted long observing runs (Baglin et al. 2006). Solarlike 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 signaltonoise ratio (HD 175726: Mosser et al. 2009b; HD 181906: García et al. 2009), or because of a too short lifetime of the modes, inducing large mode linewidths (HD 49933: Appourchaux et al. 2008; HD 181420: Barban et al. 2009).
The star HD 49385 is the first Gtype solarlike pulsator observed in the seismology field of CoRoT. It is cooler than the solarlike 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 lowfrequency part of the power spectrum, in search of a signature of the stellar rotation. The extraction of pmode parameters is described in Sect. 5, and Sect. 6 is dedicated to conclusions.
2 Stellar parameters of HD 49385
The solarlike pulsator HD 49385 is a G0type star with an apparent magnitude of (uvby catalog, Hauck & Mermilliod 1998). To determine the fundamental parameters of HD 49385 we analysed a highquality spectrum from the NARVAL spectrograph mounted on the 2m Bernard Lyot Telescope at the Pic du Midi Observatory. We stacked two spectra collected on 2009 January 10, totalizing an exposure time of 6600 s. The spectrum was normalized by identifying continuum windows in a synthetic spectrum and fitting a loworder spline through these points. This was done orderbyorder and we made sure the line depths agreed for the overlapping part of adjacent echelle orders. We measured a mean signaltonoise (S/N) ratio of 600 in the continuum for several linefree regions in the range 50007000 Å for a data sampling of 2.7 data points per resolution element (R=65 000).
2.1 Temperature, and metalicity
We used the semiautomatic VWA software (Bruntt 2009) to fit synthetic profiles for more than 600 lines in the range from 4135 to 8545 Å. We used a differential approach, meaning all abundances are measured relative to exactly the same lines in the solar spectral atlas from Kurucz et al. (1984). This differential approach is described in more detail by Bruntt et al. (2008). 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 Bruntt et al. 2008, 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 6122.2 and 6162.2 Å. 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.
Table 1: Spectroscopic parameters of HD 49385.
The parameters agree well with the calibration of the Strömgren indices. From TEMPLOGG (Rogers 1995; Kupka & Bruntt 2001) we get K, , and . TEMPLOGG determines a significant interstellar reddening . If we assume zero interstellar reddening, we get from the VK, using the Ramírez & Meléndez (2005) calibration, K and K from by. The uncertainties include the calibration uncertainty and the photometric uncertainty on the indices. Finally, Kovtyukh et al. (2004) used linedepth 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 linedepth 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., Boesgaard et al. 2005). We used the line list from Ghezzi et al. (2009), 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 (Asplund 2005). 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. 

Open with DEXTER 
2.2 Projected rotational velocity
Because we have a relatively highresolution 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 46006800 Å. We calculated the synthetic profiles with the SYNTH software Valenti & Piskunov (1996), 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 R=65 000, 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 R for the stellar lines. For and we calculated a grid of convolved profiles for 08 km s^{1} in steps of 0.5 km s^{1}. 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^{1}, corresponding to km s^{1}. The values that provide the best fit based on the analysis of the contours for all ten lines are km s^{1} and km s^{1}.
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. 

Open with DEXTER 
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). 

Open with DEXTER 
Table 2: Abundances of 22 elements in HD 49385 relative to the Sun.
2.3 Luminosity
Using the Hipparcos parallax of the star, mas (van Leeuwen 2007), we could derive the absolute visual magnitude of the object . The bolometric correction was determined by interpolating in the grid provided by Bessell et al. (1998) for appropriate values of temperature and for HD 49385. We obtained , which yields a luminosity . The values of and L 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 Samadi et al. (2007a), the lightcurve consists of a series of equallyspaced 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 1000 s on either sides of the gap. A lowfrequency trend due to the aging of the CCD was also corrected by removing a linear gain variation of 5.32 10 ^{5} day^{1} (see Auvergne et al. 2009 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 Auvergne et al. 2009) 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. 

Open with DEXTER 
Figure 5: Power spectrum of 137 days of CoRoT data for HD 49385. The peaks indicated by arrows are from left to right: the 14Hzpeak discussed in Sect. 4, a peak at Hz, and a peak at Hz. The grey area shows the excess of power due to the pmode oscillations. 

Open with DEXTER 
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 pmode parameters (see Appourchaux et al. 2008; Barban et al. 2009; Mosser et al. 2009b). It indeed provides a direct measurement of the surface rotation period and hence a first estimate of the rotational splitting value.
Contrary to other solarlike 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 timefrequency analysis with Morlet wavelets (as described in Mathur et al. 2010) which showed that this peak presents a certain stability in time. But because the periods we were testing correspond to a nonnegligible fraction of the total observation period, the cone of influence (Torrence & Compo 1998) 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 signaltonoise 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^{1}. To obtain the measured , we would need a value of the inclination angle lower than , which means that the star is seen poleon. In this case it is extremely difficult to derive the signature of the activity from the light curve (Mosser et al. 2009a) 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 seismofield 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 lowfrequency 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 poleon.
5 Analysis of p mode oscillations
Figure 6 shows the power spectrum smoothed with a 10Hz boxcar. The comblike structure of this excess of power is typical of pmode 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). 

Open with DEXTER 
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 comblike 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 1Hz 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. 

Open with DEXTER 
The pattern we observe in Fig. 7 is the one we expect to observe for lowdegree highradialorder modes in solarlike 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 Mosser & Appourchaux (2009) and Bedding & Kjeldsen (2010) 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 solarlike targets encountered ambiguity indentifying the modes (see Barban et al. 2009, García et al. 2009), 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 8: Échelle diagram of the peaks with a significance greater than 95% in the power spectrum smoothed over n=15 bins (Hz). Four peaks appear, which do not belong to the identified ridges A, B, or C, and are labeled as . 

Open with DEXTER 
Figure 7 also shows fainter peaks outside the identified ridges. We here study their significance by testing the H_{0} hypothesis in the frequency range of the oscillations. Because we searched for shortlived modes, we binned the power spectrum over n bins, as prescribed in Appourchaux (2004). The size of the boxcar (n=15 bins Hz) for the binning was chosen to approximately match the expected linewidth for the p modes in solarlike pulsators. In this case, the noise is distributed as a with 2n degrees of freedom. We then applied the H_{0} hypothesis, i.e. we supposed that the observations are due to pure noise. For each peak we computed the probability p that pure noise is able to reach the height of the peak over the whole studied window. Low values of p indicate a poor compatibility between the observations and noise, and the significance of the peaks is defined as 1p. 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 Deheuvels & Michel (2009) that a lowdegree
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
(1) 
where the expression of D_{0} is found in ChristensenDalsgaard (2003). 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 Appourchaux et al. (2009). To confirm the existence of a signal, one needs to compute the posterior probability P(H_{0}X=x). 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 solarlike 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 m=+1 component of the rotationally splitted mode. This hypothesis is studied in Sect. 5.3.4.
5.3 Extraction of p mode parameters
5.3.1 Method
Solarlike 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 Duvall & Harvey 1986).
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 Anderson et al. 1990). This method was already
applied to determine pmode parameters for the Sun (Toutain & Froehlich 1992; Appourchaux et al. 1998) as well as for CoRoT solarlike targets
(Appourchaux et al. 2008; Barban et al. 2009; García et al. 2009).
As described in Anderson et al. (1990), the likelihood L is defined as the product of the probability density of each frequency bin
of the power spectrum
(2) 
where S is the observed power spectrum, M 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 L is the same as minimizing the negative logarithmic likelihood , defined as
(3) 
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 pmode parameters are not influenced by those peaks.
We estimated the background with a power law to describe the granulation component, as prescribed by Harvey (1985). Prior to the fit of the pmode parameters,
we fitted a profile of the type
(4) 
on the power spectrum, excluding the frequency domain of the pmode 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 solarlike CoRoT
target, showed that a global fit of all the modes gave more robust
estimates of the mode parameters
(see Appourchaux et al. 2008). 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 Appourchaux et al. (2008). We briefly recall them here.
The splitting was assumed to be constant over the frequency range of the oscillations, and the
components of nonradial modes were supposed to be symmetrical
with respect to the m=0 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 nonradial 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 nonradial multiplets, the heights are given by the intensity visibilities (see Gizon & Solanki 2003; Ballot et al. 2006).
5.3.2 Mode parameters estimates
Five teams worked independently on the extraction of the pmode 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 crosscheck them and comment the robustness of the solution. The reference results of the fit are given in Tables 3 and 4.
Table 3: Reference results for the parameters of the radial p modes.
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. 

Open with DEXTER 
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 signaltonoise ratio decreases when reaching the edges of the fitted domain. Indeed, with a lower signaltonoise ratio, the analysis of previous CoRoT targets yielded several modes merged with the neighboring mode (see Appourchaux et al. 2008; Barban et al. 2009). 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 lowfrequency avoided crossing.
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. 

Open with DEXTER 
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 solarlike 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). 

Open with DEXTER 
Mode heights and amplitudes
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 Michel et al. 2009), 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 solarlike pulsators observed with CoRoT for which amplitudes were found significantly lower than the predicted ones (see Appourchaux et al. 2008; Barban et al. 2009; García et al. 2009; Michel et al. 2008).
The signaltonoise (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). 

Open with DEXTER 
Table 4: Reference results for the parameters of the nonradial 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 poleon, i.e. with an inclination angle close to .
In that case the
components of the nonradial
modes are not visible, and we therefore cannot estimate the splitting. However, the measured
of the object and the estimate of the stellar radius R 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 i 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 nonradial modes. Equation 5 shows that a small rotational splitting corresponds to a high inclination angle.
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 m=0 component, or the splitting is small, and the m 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 parameter 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 gmode behavior in the deep interior and a pmode behavior below the surface (see Scuflaire 1974). The mixed modes are known to experience the socalled avoided crossing phenomenon, during which their frequency does not follow the asymptotic approximation and they can be located out of the ridges (Unno et al. 1989). It was already noticed (e.g. Dziembowski & Pamyatnykh 1991) 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 gmode 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 illdetermined.
Table 5: Results of a Lorentzian fit performed on the peaks identified as potential mixed modes.
5.3.4 Other mode identification scenarii
Table 6: Reference results for the parameters of the radial p modes for the alternate identification (symbols are the same as in Table 3).
Figure 13: Frequency shift (defined in Eq. (6)) of the components of nonradial 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). 

Open with DEXTER 
Can ridge D be the m=+1 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 m=+1 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 nonradial multiplets are no longer symmetric with respect to the m=0 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 L that we obtained for HD 49385 in Sect. 2.1 and with different angular velocities.
The models were computed with the evolutionary code CESAM2k (Morel 1997)
and the eigenfrequencies were derived from them with the oscillation code FILOU (Suárez & Goupil 2008).
They showed that the effects of asymmetry should be nonnegligible 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 n 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
(7) 
where n_{0} corresponds to the order where the signal is maximal (n_{0}=17). For modes in nonradial multiplets we now fit eigenfrequencies of the type
(8) 
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 m=+1 component of the modes. The obtained rotational shift is such that Hz per order n, and Hz. However, using Eq. (5) and the fitted values of the splitting and the angle, we obtain km s^{1}. This result is not compatible with the firm upper limit of 5 km s^{1} established for in the spectroscopic analysis (see Sect. 2.2).
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 pmode oscillations in solarlike 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^{1}, in complete
disagreement with the measured
km s^{1}.
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 pmode parameters.
Figure 14: Échelle diagram for the fitted modes of the alternate identification. Same symbols as in Fig. 10. 

Open with DEXTER 
Table 7: Reference results for the parameters of the nonradial 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 137daylong 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 solarlike 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 solarlike 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 Samadi et al. (2007b), unlike the other CoRoT solarlike pulsators (e.g. Appourchaux et al. 2008).
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 nonradial multiplets are not expected to be symmetrical with respect to the m=0 component. We proposed a simple way of treating this assymetry in future analyses of solarlike 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 solarlike 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 Gtype pulsators are smaller than those of Ftype 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 lowfrequency eigenmodes. In particular we found that the ridge is distorted compared to the ridge at low frequency. Deheuvels & Michel (2009) showed that this type of distortion could be associated with a lowdegree 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.
AcknowledgementsThis 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 nonadiabatic 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. Kjeldsen & Bedding (1995) established the following relation between the luminosity amplitude A 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.
Samadi et al. (2007b) obtained
Inserting Eqs. (A.2) and (A.3) into Eq. (A.1) we obtain
(A.4) 
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_{0} and H_{1}, Appourchaux et al. (2009) obtained the following expression for the posterior probability
of the H_{0} hypothesis given the observed data:
The expression of P(xH_{0}) for a spectrum binned over n bins is given in Appourchaux (2004)
where is the Gamma function. To solve Eq. (B.1), we need to define the alternate hypothesis H_{1} 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 P(xH_{1}), 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 nonadiabatic 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 Appourchaux et al. (2009). 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 ( [0,5] Hz).
The height of the
modes can be deduced from the amplitude of
modes through the relation
(B.3) 
for an angle of inclination i. The visibility factors are computed from Gizon & Solanki (2003). 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 Appourchaux et al. (2009) we obtain the following expression for the probability of the observed data given the signal
where is the normalization factor. and are given in Appourchaux (2004) and are functions of the mode height H and linewidth .
Inserting Eqs. (B.2) and (B.4) into Eq. (B.1), 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 gmode 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 ( [0,5] Hz).
We obtain posterior probabilities of 0.001% and 8.7% for and , respectively.
References
 Anderson, E. R., Duvall, Jr., T. L., & Jefferies, S. M. 1990, ApJ, 364, 699 [NASA ADS] [CrossRef] [Google Scholar]
 Appourchaux, T. 2004, A&A, 428, 1039 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Appourchaux, T., Gizon, L., & RabelloSoares, M.C. 1998, A&AS, 132, 107 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Appourchaux, T., Michel, E., Auvergne, M., et al. 2008, A&A, 488, 705 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Appourchaux, T., Samadi, R., & Dupret, M. A. 2009, A&A, 506, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Asplund, M. 2005, ARA&A, 43, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Baglin, A., Auvergne, M., Barge, P., et al. 2006, in ESA SP 1306, ed. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, 33 [Google Scholar]
 Ballot, J., García, R. A., & Lambert, P. 2006, MNRAS, 369, 1281 [NASA ADS] [CrossRef] [Google Scholar]
 Barban, C., Deheuvels, S., Baudin, F., et al. 2009, A&A, 506, 51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Basu, S., ChristensenDalsgaard, J., Howe, R. et al. 2003, ApJ, 591, 432 [NASA ADS] [CrossRef] [Google Scholar]
 Bedding, T. R., & Kjeldsen, H. 2010, Commun. Asteroseismol., 161, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Bessell, M. S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 [NASA ADS] [Google Scholar]
 Boesgaard, A. M., Stephens, A., & Deliyannis, C. P. 2005, ApJ, 633, 398 [NASA ADS] [CrossRef] [Google Scholar]
 Bruntt, H. 2009, A&A, 506, 235 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bruntt, H., De Cat, P., & Aerts, C. 2008, A&A, 478, 487 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 ChristensenDalsgaard, J. 2003, Lect. Not. Stellar Oscillat. [Google Scholar]
 ChristensenDalsgaard, J., Gough, D. O., & Thompson, M. J. 1991, ApJ, 378, 413 [NASA ADS] [CrossRef] [Google Scholar]
 Deheuvels, S., & Michel, E. 2009, Ap&SS, 241 [Google Scholar]
 Duvall, Jr., T. L., & Harvey, J. W. 1986, in NATO ASIC Proc. 169: Seismology of the Sun and the Distant Stars, ed. D. O. Gough, 105 [Google Scholar]
 Dziembowski, W. A., & Pamyatnykh, A. A. 1991, A&A, 248, L11 [NASA ADS] [Google Scholar]
 García, R. A., Régulo, C., Samadi, R., et al. 2009, A&A, 506, 41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ghezzi, L., Cunha, K., Smith, V. V., et al. 2009, ApJ, 698, 451 [NASA ADS] [CrossRef] [Google Scholar]
 Gizon, L., & Solanki, S. K. 2003, ApJ, 589, 1009 [NASA ADS] [CrossRef] [Google Scholar]
 Harvey, J. 1985, in Future Missions in Solar, Heliospheric & Space Plasma Physics, ed. E. Rolfe & B. Battrick, ESA SP, 235, 199 [Google Scholar]
 Hauck, B., & Mermilliod, M. 1998, A&AS, 129, 431 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Kjeldsen, H., & Bedding, T. R. 1995, A&A, 293, 87 [NASA ADS] [Google Scholar]
 Kovtyukh, V. V., Soubiran, C., & Belik, S. I. 2004, A&A, 427, 933 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kupka, F., & Bruntt, H. 2001, Using TEMPLOGG for determining stellar parameters of MONS targets, ed. C. Sterken (C), 39 [Google Scholar]
 Kurucz, R. L., Furenlid, I., Brault, J., & Testerman, L. 1984, Solar flux atlas from 296 to 1300 nm, ed. R. L. Kurucz, I. Furenlid, J. Brault, & L. Testerman [Google Scholar]
 Mathur, S., EffDarwich, A., García, R. A., & TurckChièze, S. 2008, A&A, 484, 517 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathur, S., García, R. A., Régulo, C., et al. 2010, A&A, 511, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Michel, E., Baglin, A., Auvergne, M., et al. 2008, Science, 322, 558 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Michel, E., Samadi, R., Baudin, F., et al. 2009, A&A, 495, 979 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Morel, P. 1997, A&AS, 124, 597 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Mosser, B., & Appourchaux, T. 2009, A&A, 508, 877 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Baudin, F., Lanza, A. F., et al. 2009a, A&A, 506, 245 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Michel, E., Appourchaux, T., et al. 2009b, A&A, 506, 33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ramírez, I., & Meléndez, J. 2005, ApJ, 626, 465 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Rogers, N. Y. 1995, Commun. Asteroseismol., 78, 1 [Google Scholar]
 Samadi, R., Fialho, F., Costa, J. E. S., et al. 2007a, [astroph/0703354] [Google Scholar]
 Samadi, R., Georgobiani, D., Trampedach, R., et al. 2007b, A&A, 463, 297 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Scuflaire, R. 1974, A&A, 36, 107 [NASA ADS] [Google Scholar]
 Suárez, J. C., & Goupil, M. J. 2008, Ap&SS, 316, 155 [NASA ADS] [CrossRef] [Google Scholar]
 Thompson, M. J., ChristensenDalsgaard, J., Miesch, M. S., & Toomre, J. 2003, ARA&A, 41, 599 [NASA ADS] [CrossRef] [Google Scholar]
 Torrence, C., & Compo, G. P. 1998, Bull. Am. Meteorol. Soc., 79, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Toutain, T., & Froehlich, C. 1992, A&A, 257, 287 [NASA ADS] [Google Scholar]
 TurckChièze, S., Couvidat, S., Kosovichev, A. G., et al. 2001, ApJ, 555, L69 [NASA ADS] [CrossRef] [Google Scholar]
 Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial oscillations of stars (University of Tokyo Press), 2nd ed. [Google Scholar]
 Valenti, J. A., & Piskunov, N. 1996, A&AS, 118, 595 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Footnotes
 ... HD 49385^{}
 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.
 ...^{}
 Based on data obtained using the Télescope Bernard Lyot at Observatoire du Pic du Midi, CNRS and Université Paul Sabatier, France.
 ... CoRoT^{}
 The CoRoT space mission, launched on 2006 December 27, was delopped an d is operated by the CNES with participation of the Science Programs of ESA; ESA's RSSD, Austria, Belgium, Brazil, Germany and Spain.
All Tables
Table 1: Spectroscopic parameters of HD 49385.
Table 2: Abundances of 22 elements in HD 49385 relative to the Sun.
Table 3: Reference results for the parameters of the radial p modes.
Table 4: Reference results for the parameters of the nonradial p modes (symbols are the same as in Table 3).
Table 5: Results of a Lorentzian fit performed on the peaks identified as potential mixed modes.
Table 6: Reference results for the parameters of the radial p modes for the alternate identification (symbols are the same as in Table 3).
Table 7: Reference results for the parameters of the nonradial p modes for the alternate identification (symbols are the same as in Table 3).
All Figures
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. 

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

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

Open with DEXTER  
In the text 
Figure 4: Lightcurve of 137 days of CoRoT observations on HD 49385, detrended using instrumental information. 

Open with DEXTER  
In the text 
Figure 5: Power spectrum of 137 days of CoRoT data for HD 49385. The peaks indicated by arrows are from left to right: the 14Hzpeak discussed in Sect. 4, a peak at Hz, and a peak at Hz. The grey area shows the excess of power due to the pmode oscillations. 

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

Open with DEXTER  
In the text 
Figure 7: Échelle diagram of HD 49385 computed with a large spacing of Hz. The power spectrum is smoothed with a 1Hz 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. 

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

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

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

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

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

Open with DEXTER  
In the text 
Figure 13: Frequency shift (defined in Eq. (6)) of the components of nonradial 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). 

Open with DEXTER  
In the text 
Figure 14: Échelle diagram for the fitted modes of the alternate identification. Same symbols as in Fig. 10. 

Open with DEXTER  
In the text 
Copyright ESO 2010