Free Access
Issue
A&A
Volume 534, October 2011
Article Number L3
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201117349
Published online 29 September 2011

© ESO, 2011

1. Introduction

The F5V solar-like pulsator HD 49933 (HR 2530, HIP 32851) is one of the primary targets of the CoRoT space mission (Baglin et al. 2006). The discovery of its oscillation was made from the ground by observing highly resolved time-series spectra (Mosser et al. 2005). Since the launch of the spacecraft in December 2006, two runs of 60 and 137 days have been performed leading to a rich p-mode spectrum of 51 frequencies corresponding to the  = 0,1,2 eigenmodes (Appourchaux et al. 2008; Benomar et al. 2009; Kallinger et al. 2010).

thumbnail Fig. 1

Monochromatic (730 nm) center-to-limb emergent intensity at a representative time of the simulation for various cosines of the viewing angle: μ = 1.0 (disk center) down to 0.1 (towards the limb). The horizontal sizes of the simulation domain are 21   000  ×  21   000 km.

Open with DEXTER

The interpretation of these eigenmodes to derive constraints on the stellar interior depends on the accuracy of the determination of the fundamental parameters of the star (Teff,log    g, [Fe/H] ). All chemical abundance analyses of HD 49933 converge to that of a slightly metal-poor star with  [Fe/H]  between −0.3 and −0.5 (e.g. Gillon & Magain 2006; Bruntt 2009). The determination of its effective temperature and gravity have been revisited several times. As reported in Bruntt (2009), the Teff obtained by different methods ranges from 6450 to 6780 K, which is by far too uncertain to tightly constrain the stellar evolution models in the HR diagram. Similarly, log    g is also not very well constrained with values ranging from 4.0 to 4.3. Several authors have attempted to determine Teff and log    g by fitting the Hβ and Mg I lines, respectively (Ryabchikova et al. 2009; Kallinger et al. 2010). Their solutions are most consistent with a cooler temperature (≈ 6500 K) and a lower gravity (≈ 4.0).

Since spectrometry and photometry lead to stellar parameters with a large uncertainty, in this work we use a different approach by determining its diameter using interferometry. Combined with seismology, it provides a very good constraint of the stellar parameter space (Creevey et al. 2007). This approach was proposed for the first time by Kervella et al. (2003) for the binary star αCen AB and later with similar success for other asteroseismic targets (Cunha et al. 2007; Teixeira et al. 2009; Bruntt et al. 2009; Mazumdar et al. 2009; Bazot et al. 2011).

In this letter we present the first campaign of interferometric measurements of HD 49933 obtained using the VEGA instrument (Mourard et al. 2009, 2011) at the CHARA Array (ten Brummelaar et al. 2005). We used state-of-art 3D hydrodynamical simulations of the surface of the star to calculate the limb darkening and to compute the visibility curves prior to extracting the angular diameter. The derived stellar radius is then used to constrain the other fundamental stellar parameters by fitting the oscillation frequency separations.

2. CHARA/VEGA observations

Observations of HD 49933 were performed on October 16, 2010 using the VEGA instrument of the CHARA Array and the instrument CLIMB (Sturmann et al. 2010) for 3T group delay tracking. The telescopes E1, E2, and W2 were used giving access to ground baselines of 66, 156, and 221 m. The seeing was stable with a value of r0 between 10 and 13 cm. For the absolute calibration of the squared visibilities, we observed two different calibrators, C1 = HD 55185 and C2 = HD 46487. The observing sequence was C1-T-C2-T-C2-T-C1, each block being almost 20 mn. The data were processed using the standard V2 procedure of the VEGA instrument as described in Mourard et al. (2011). To validate the absolute calibration, we first considered C1 as a target and thus estimated its uniform disk (UD) diameter using C2 as the calibrator. A spectral band of Δλ = 30 nm centered around λ = 735 nm was used for the data processing. The equivalent UD diameter of C2 is taken as θHD   46487 = 0.180  ±  0.013 mas according to the SearchCal tool1 of the JMMC (Bonneau et al. 2006). Fitting a UD model to the three squared visibilities provides an estimate of the angular diameter of HD 55185 of θHD   55185 = 0.474  ±  0.014 mas. Spectral bandwidth and baseline smearing are included but only introduce a bias less than 0.1% on the diameter determination. Surface brightness relations, as provided by SearchCal, provide an independent estimate of the angular diameter. We found 0.478 mas and 0.476 mas from the (V − R) and (B − V) relations, respectively. Our result agrees very well with these photometric estimations. We finally use C1 and C2 to calibrate the measurements on the target. To avoid a strong telluric feature, data are processed using a spectral band of Δλ = 20 nm centered around λ = 740 nm. The final error bars in the target measurements take into account the measurement noise, the dispersion in the individual measurements, the diameter of the calibrators, and their errors.

3. Determination of the diameter

To derive a reliable diameter of HD 49933, we base our analysis on a 3D hydrodynamical simulation of the surface of the star. Allende Prieto et al. (2002) demonstrated that for Procyon, a star similar to HD 49933, the 3D limb darkening is significantly different from the corresponding 1D hydrostatic and homogeneous one. The differences were found to be larger in the visible (the VEGA domain) than the infrared. We adopt the same methodology as used for both α Cen B (Bigot et al. 2006) and Procyon (Aufdenberg et al. 2005). A hydrodynamical simulation of the surface is used to obtain 3D time-dependent limb-darkened intensities, instead of the simplified 1D homogeneous or UD models, to compute the visibility curves. The angular diameter is found by minimizing the differences between the synthetic and observed visibility curves.

3.1. The 3D simulation of the surface

We use the state-of-the-art radiative hydrodynamical code (STAGGER CODE, Nordlund & Galsgaard2) to simulate the surface convection and stratification of HD 49933. In a local box, the code solves the full set of conservative hydrodynamical equations coupled to an accurate treatment of the radiative transfer. The code is based on a sixth order explicit finite difference scheme. The equations are solved on a staggered mesh where the thermodynamical variables are cell centered, while the fluxes are shifted to the cell edge. The domain of simulation contains the entropy minimum located at the surface and is extended deep enough to have a flat entropy profile at the bottom (adiabatic regime). The code uses periodic boundary conditions horizontally and open boundaries vertically. At the bottom of the simulation, the inflows have constant entropy and pressure. The outflows are not constrained and are free to pass through the boundary. We used a realistic equation-of-state that accounts for ionization, recombination, and dissociation (Mihalas et al. 1988) and continuous + line opacities (Gustafsson et al. 2008). Radiative transfer is crucial since it drives convection through entropy losses at the surface (Stein & Nordlund 1998) and is solved using the Feautrier’s scheme along several inclined rays (one vertical, eight inclined) through each grid point. The wavelength dependence of the radiative transfer is taken into account using a binning scheme in which the monochromatic lines are collected into 12 bins. The numerical resolution used for the present simulation is 2403. The geometrical sizes are 21  ×  21 Mm horizontally and 5 Mm vertically. The horizontal sizes of the domain are defined to contain a sufficient number of granules at each time-step and the vertical one to be deep enough to ensure an adiabatic regime at the bottom. The stellar parameters that define our 3D model are Teff = 6690  ±  35 K, log    g = 4.21 and a scaled solar chemical composition (Asplund et al. 2009) down to −0.5 dex. The uncertainty in Teff represents the fluctuations with time around the mean value. The simulation was run for a couple of stellar hours to get the full hydrodynamical relaxation. The convection in such an F star is more efficient than for the Sun: the rms vertical velocity is 3.9 km s-1, which is about twice the solar value. The granulation pattern is shown in Fig. 1.

3.2. Limb darkening and visibility curves

thumbnail Fig. 2

Our best-fit model of the observed squared visibilities (black dots) with the calculated one (full line) with a reduced χ2 = 0.47. The angular diameter is θLD = 0.445 mas.

Open with DEXTER

Hereafter, we use the methodology presented in Bigot et al. (2006). The snapshots are taken after relaxation and cover a sequence of several convective turnover times. We extract snapshots every five minutes over a sequence of one hour. For each snapshot, we calculate the radiative transfer in the selected wavelength bandwidth of the VEGA instrument Δλ =  [730−750]  nm to get the monochromatic limb-darkened intensities for each time step, each position at the stellar surface, and each cosine μ of the angle made between the vertical and the line of sight.

The normalized fringe visibilities are obtained using the Cittert-Zernike theorem (1)where B is the projected baseline (in meters), J0 is the zero order of the Bessel function, and  ⟨ Iλ(μ) ⟩  is the time- and horizontal-average of the 3D intensities. We integrate the monochromatic visibilities over the spectral domain Δλ of the instrument weighted by the spectral transmission T(λ) (2)The function T(λ) is derived from the direct flux measurement of the target and checked on several bright stars.

3.3. The χ2 fit of the data

The angular diameter is determined by a Levenberg-Marquardt least-squares fit to the observed squared visibilities with the model. The result is shown in Fig. 2. The value derived from the 3D limb darkening is θLD = 0.445  ±  0.012 mas. The reduced χ2 of the fit is 0.47. Using the revised Hipparcos parallax π = 33.69  ±  0.42 mas (van Leeuwen 2007), we obtain a linear radius of R = 1.42  ±  0.04   R. The corresponding radius derived from 1D hydrostatic limb-darkening (Claret 2000) is R = 1.45  ±  0.04   R, which is significantly larger by 1σ because of the more pronounced limb darkening in 1D. We note that the derived radius is not very sensitive to the choice of stellar parameters (those selected in Sect. 3.1). We indeed tried several values of Teff (= 6500,6750   K) and log    g (=4.0,4.5) and found very small differences in terms of derived radius (≤ 0.002   R). The derived angular diameter θLD agrees with the prediction of 0.452  ±  0.007 mas obtained by the surface brightness relation (Thévenin et al. 2006). The relative uncertainty in the angular diameter (2.7%) is larger by an order of magnitude than the smallest uncertainty ever achieved for an interferometric radius, (0.2%) derived in the case of the α Centauri A (Kervella et al. 2003). However, we emphasize that the target in the present work is 20 times smaller. This uncertainty dominates the total uncertainty of about 3%, including the error in the parallax. Nonetheless, the relative uncertainty in the radius is sufficiently accurate to tightly constrain the mass using asteroseismic frequencies (Creevey et al. 2007).

4. Asteroseismic fits of acoustic frequencies

Table 1

Our stellar evolution model for HD 49933.

thumbnail Fig. 3

Large (left panel) and small (right panel) separations for  = 0 (blue) and  = 1 (red) calculated with the stellar model of Table 1. The observed data (bullets) are from Benomar et al. (2009). For each frequency, we took the maximum error bar at 1σ.

Open with DEXTER

We estimate the mass and the age of HD 49933 using asteroseismic data. The aim is to find a stellar evolutionary model whose radius equals the interferometric one, whose Teff, log    g, and  [Fe/H]  match the spectrometric ones (Bruntt 2009), and whose eigenmode frequency separations match the CoRoT data (Benomar et al. 2009). We compute a grid of models of mass between 1.7 and 2.1 M, and search for the one that most closely matches the observed small and large frequency separations. The stellar model evolves until the radius equals the interferometric one. The initial metallicity is adjusted to match the spectroscopic  [Fe/H] . The stellar evolutionary models are computed using the CESAM2k code (Morel 1997; Morel & Lebreton 2008). In addition, we adopt Chemical diffusion (Burgers 1969) combined with extra-mixing (Morel & Thévenin 2002) and the Canuto & Mazzitelli (1991) formalism of convection. The rotation of HD 49933 has a minor impact on the modeling (Piau et al. 2009) and is ignored. For each stellar model, we compute eigenfrequencies between 1100 and 2500 μHz for different Legendre polynomial degrees  = 0,1,2 by solving3 the standard linearized equations of continuity and momentum for adiabatic oscillations (see e.g. Unno et al. 1989). Since these modes pulsate in high overtones (n), we added the shift in frequency caused by the centrifugal distortion of the star (Dziembowski & Goode 1992).

We fit both the large and small frequency separations that we define as follows (3)Since the radius is given, the large separations are mainly sensitive to the mass of the star. The small separations are more sensitive to the deep interior and are therefore a good indicator of the stellar age (e.g. Roxburgh & Vorontsov 2003).

The solution that minimizes the differences between the observed and computed separations is given in Table 1 and the corresponding frequency separations are shown in Fig. 3. We note that a core overshoot of 0.35 is needed to fit the small separations. We get a better fit of the small separations by including the centrifugal effects on the oscillations. The fits to both small and large separations are very satisfactory and better than in Benomar et al. (2010). To roughly estimate the errors in the mass and Teff, we use the observed errors in the average large separation (~M1/2R−3/2), in the luminosity L (), and in the interferometric radius. Our stellar evolution model agrees well with one of the solutions of Creevey & Bazot (2011), who used a local minimization method based on radial modes. The mass (1.20  ±  0.08   M) is close to that of other works (Appourchaux et al. 2008; Benomar et al. 2010; Kallinger et al. 2010). The luminosity log    L/L = 0.55 matches the one observed (Michel et al. 2008). The derived atmospheric parameters (Teff = 6640  ±  100 K, log    g  =  4.21  ±  0.14,  [Fe/H]  = −0.38) are in good agreement with the values obtained by Strömgren photometry (Bruntt 2009). Our Teff is larger than the value (≈ 6500 K) derived from Hβ (Ryabchikova et al. 2009; Kallinger et al. 2010). However, we emphasize that this strong line is not an excellent temperature indicator since it spreads over a wide wavelength range.

5. Conclusion

We have reported our long-baseline interferometric observations of the CoRoT target HD 49933 using the VEGA instrument at the CHARA Array. Combined with 3D hydrodynamical simulations of the surface, we derived its radius, R = 1.42  ±  0.04   R, which agrees with the radius derived from the surface brightness relation. The use of this interferometric radius reduces the stellar parameter space and helped us to more tightly constrain the position of the star in the HR diagram. The fit of the p-mode frequencies led us to determine the other stellar parameters. The derived atmospheric ones (Teff = 6640  ±  100 K, log    g  =  4.21  ±  0.14,  [Fe/H]  = −0.38) agree well with photometric studies.


3

Code kindly provided by W. Dziembowski.

Acknowledgments

L.B. thanks A. Chiavassa, R. Collet, and R. Trampedach who provided the latest EOS and opacities for the 3D modeling, and J. Provost for valuable discussions about the fitting of frequencies and the stellar evolution modeling. Computations have been done on the “Mesocentre SIGAMM” machine, hosted by the Observatoire de la Côte d’Azur. The CHARA Array is operated with support from the National Science Foundation through grant AST-0908253, from Georgia State University, the W. M. Keck Foundation, and the NASA Exoplanet Science Institute.

References

All Tables

Table 1

Our stellar evolution model for HD 49933.

All Figures

thumbnail Fig. 1

Monochromatic (730 nm) center-to-limb emergent intensity at a representative time of the simulation for various cosines of the viewing angle: μ = 1.0 (disk center) down to 0.1 (towards the limb). The horizontal sizes of the simulation domain are 21   000  ×  21   000 km.

Open with DEXTER
In the text
thumbnail Fig. 2

Our best-fit model of the observed squared visibilities (black dots) with the calculated one (full line) with a reduced χ2 = 0.47. The angular diameter is θLD = 0.445 mas.

Open with DEXTER
In the text
thumbnail Fig. 3

Large (left panel) and small (right panel) separations for  = 0 (blue) and  = 1 (red) calculated with the stellar model of Table 1. The observed data (bullets) are from Benomar et al. (2009). For each frequency, we took the maximum error bar at 1σ.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.