Issue 
A&A
Volume 534, October 2011



Article Number  L3  
Number of page(s)  4  
Section  Letters  
DOI  https://doi.org/10.1051/00046361/201117349  
Published online  29 September 2011 
Letter to the Editor
The diameter of the CoRoT target HD 49933
Combining the 3D limb darkening, asteroseismology, and interferometry
^{1}
Université NiceSophia Antipolis, Observatoire de la Côte d’Azur, CNRS UMR
6202,
BP 4229,
06304
Nice Cedex,
France
email: lbigot@oca.eu
^{2}
Université NiceSophia Antipolis, Observatoire de la Côte d’Azur,
CNRS UMR 6525, BP
4229, 06304
Nice Cedex,
France
^{3}
UCBL/CNRS CRAL, 9 avenue Charles André, 69561
Saint Genis Laval Cedex,
France
^{4}
UJFGrenoble 1/CNRSINSU, Institut de Planétologie et
d’Astrophysique de Grenoble, UMR 5274, 38041
Grenoble,
France
^{5}
CHARA and Department of Physics and Astronomy, Georgia State
University, PO Box
4106, Atlanta,
GA
303024106,
USA
^{6}
LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Université
Paris Diderot, 5 place Jules
Janssen, 92195
Meudon,
France
^{7}
Georgia State University, PO Box 3969, Atlanta
GA
303023969,
USA
^{8}
CHARA Array, Mount Wilson Observatory,
91023
Mount Wilson
CA,
USA
^{9}
National Optical Astronomy Observatory,
PO Box 26732, Tucson, AZ
85726,
USA
Received: 26 May 2011
Accepted: 22 August 2011
Context. The interpretation of stellar pulsations in terms of internal structure depends on the knowledge of the fundamental stellar parameters. Longbase interferometers permit us to determine very accurate stellar radii, which are independent constraints for stellar models that help us to locate the star in the HR diagram.
Aims. Using a direct interferometric determination of the angular diameter and advanced threedimensional (3D) modeling, we derive the radius of the CoRoT target HD 49933 and reduce the global stellar parameter space compatible with seismic data.
Methods. The VEGA/CHARA spectrointerferometer is used to measure the angular diameter of the star. A 3D radiative hydrodynamical simulation of the surface is performed to compute the limb darkening and derive a reliable diameter from visibility curves. The other fundamental stellar parameters (mass, age, and T_{eff}) are found by fitting the large and small pmode frequency separations using a stellar evolution model that includes microscopic diffusion.
Results. We obtain a limbdarkened angular diameter of θ_{LD} = 0.445 ± 0.012 mas. With the Hipparcos parallax, we obtain a radius of R = 1.42 ± 0.04 R_{⊙}. The corresponding stellar evolution model that fits both large and small frequency separations has a mass of 1.20 ± 0.08 M_{⊙} and an age of 2.7 Gy. The atmospheric parameters are T_{eff} = 6640 ± 100 K, log g = 4.21 ± 0.14, and [Fe/H] = −0.38.
Key words: stars: oscillations / hydrodynamics / stars: fundamental parameters / techniques: interferometric / asteroseismology / stars: individual: HD 49933
© ESO, 2011
1. Introduction
The F5V solarlike 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 timeseries 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 pmode spectrum of 51 frequencies corresponding to the ℓ = 0,1,2 eigenmodes (Appourchaux et al. 2008; Benomar et al. 2009; Kallinger et al. 2010).
Fig. 1
Monochromatic (730 nm) centertolimb 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. 
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 (T_{eff},log g, [Fe/H] ). All chemical abundance analyses of HD 49933 converge to that of a slightly metalpoor 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 T_{eff} 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 T_{eff} 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 stateofart 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 C1TC2TC2TC1, each block being almost 20 mn. The data were processed using the standard V^{2} 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 tool^{1} 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 timedependent limbdarkened 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 stateoftheart radiative hydrodynamical code (STAGGER CODE, Nordlund & Galsgaard^{2}) 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 equationofstate 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 240^{3}. 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 timestep 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 T_{eff} = 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 T_{eff} 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
Fig. 2
Our bestfit 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. 
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 limbdarkened 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 CittertZernike theorem (1)where B is the projected baseline (in meters), J_{0} is the zero order of the Bessel function, and ⟨ I_{λ}(μ) ⟩ is the time and horizontalaverage 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 LevenbergMarquardt leastsquares 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 limbdarkening (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 T_{eff} (= 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
Our stellar evolution model for HD 49933.
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σ. 
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 T_{eff}, 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 extramixing (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 solving^{3} 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 T_{eff}, we use the observed errors in the average large separation (~M^{1/2}R^{−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 (T_{eff} = 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 T_{eff} 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 longbaseline 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 pmode frequencies led us to determine the other stellar parameters. The derived atmospheric ones (T_{eff} = 6640 ± 100 K, log g = 4.21 ± 0.14, [Fe/H] = −0.38) agree well with photometric studies.
Available at http://www.jmmc.fr/searchcal
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 AST0908253, from Georgia State University, the W. M. Keck Foundation, and the NASA Exoplanet Science Institute.
References
 Allende Prieto, C., Asplund, M., García López, R. J., & Lambert, D. L. 2002, ApJ, 567, 544 [NASA ADS] [CrossRef] [Google Scholar]
 Appourchaux, T., Michel, E., Auvergne, M., et al. 2008, A&A, 488, 705 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Aufdenberg, J. P., Ludwig, H., & Kervella, P. 2005, ApJ, 633, 424 [NASA ADS] [CrossRef] [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]
 Bazot, M., Ireland, M. J., Huber, D., et al. 2011, A&A, 526, L4 [Google Scholar]
 Benomar, O., Baudin, F., Campante, T. L., et al. 2009, A&A, 507, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Benomar, O., Baudin, F., Marques, J. P., et al. 2010, Astron. Nachr., 331, 956 [NASA ADS] [CrossRef] [Google Scholar]
 Bigot, L., Kervella, P., Thévenin, F., & Ségransan, D. 2006, A&A, 446, 635 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bonneau, D., Clausse, J., Delfosse, X., et al. 2006, A&A, 456, 789 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bruntt, H. 2009, A&A, 506, 235 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bruntt, H., Kurtz, D. W., Cunha, M. S., et al. 2009, MNRAS, 396, 1189 [NASA ADS] [CrossRef] [Google Scholar]
 Burgers, J. M. 1969, Flow Equations for Composite Gases (New York: Academic Press) [Google Scholar]
 Canuto, V. M., & Mazzitelli, I. 1991, ApJ, 370, 295 [NASA ADS] [CrossRef] [Google Scholar]
 Claret, A. 2000, A&A, 363, 1081 [NASA ADS] [Google Scholar]
 Creevey, O. L., & Bazot, M. 2011, J. Phys. Conf. Ser., 271, 012038 [NASA ADS] [CrossRef] [Google Scholar]
 Creevey, O. L., Monteiro, M. J. P. F. G., Metcalfe, T. S., et al. 2007, ApJ, 659, 616 [NASA ADS] [CrossRef] [Google Scholar]
 Cunha, M. S., Aerts, C., ChristensenDalsgaard, J., et al. 2007, A&AR, 14, 217 [Google Scholar]
 Dziembowski, W. A., & Goode, P. R. 1992, ApJ, 394, 670 [NASA ADS] [CrossRef] [Google Scholar]
 Gillon, M., & Magain, P. 2006, A&A, 448, 341 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kallinger, T., Gruberbauer, M., Guenther, D. B., Fossati, L., & Weiss, W. W. 2010, A&A, 510, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kervella, P., Thévenin, F., Ségransan, D., et al. 2003, A&A, 404, 1087 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mazumdar, A., Mérand, A., Demarque, P., et al. 2009, A&A, 503, 521 [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]
 Mihalas, D., Dappen, W., & Hummer, D. G. 1988, ApJ, 331, 815 [NASA ADS] [CrossRef] [Google Scholar]
 Morel, P. 1997, A&AS, 124, 597 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Morel, P., & Lebreton, Y. 2008, Ap&SS, 316, 61 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Morel, P., & Thévenin, F. 2002, A&A, 390, 611 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Bouchy, F., Catala, C., et al. 2005, A&A, 431, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mourard, D., Clausse, J. M., Marcotto, A., et al. 2009, A&A, 508, 1073 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mourard, D., Bério, P., Perraut, K., et al. 2011, A&A, 531, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Piau, L., TurckChièze, S., Duez, V., & Stein, R. F. 2009, A&A, 506, 175 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Roxburgh, I. W., & Vorontsov, S. V. 2003, A&A, 411, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ryabchikova, T., Fossati, L., & Shulyak, D. 2009, A&A, 506, 203 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stein, R. F., & Nordlund, A. 1998, ApJ, 499, 914 [Google Scholar]
 Sturmann, J., Ten Brummelaar, T., Sturmann, L., & McAlister, H. A. 2010, in SPIE Conf. Ser. 7734, ed. W. Danchi, C. William, F. Delplancke, & J. K. Rajagopal [Google Scholar]
 Teixeira, T. C., Kjeldsen, H., Bedding, T. R., et al. 2009, A&A, 494, 237 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 tenBrummelaar, T. A., McAlister, H. A., Ridgway, S. T., et al. 2005, ApJ, 628, 453 [NASA ADS] [CrossRef] [Google Scholar]
 Thévenin, F., Bigot, L., Kervella, P., et al. 2006, Mem. Soc. Astron. It., 77, 411 [Google Scholar]
 Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial oscillations of stars, 2nd ed. (Tokyo: University of Tokyo Press) [Google Scholar]
 van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
All Figures
Fig. 1
Monochromatic (730 nm) centertolimb 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. 

In the text 
Fig. 2
Our bestfit 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. 

In the text 
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σ. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.