Issue 
A&A
Volume 601, May 2017



Article Number  A127  
Number of page(s)  7  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201629825  
Published online  19 May 2017 
Asteroseismology of the exoplanethost Ftype star 94 Ceti: Impact of atomic diffusion on the stellar parameters
^{1} Laboratoire Univers et Particules de Montpellier (LUPM), UMR 5299, Université de Montpellier, CNRS, Place Eugène Bataillon, 34095 Montpellier Cedex 5 France
email: morgan.deal@umontpellier.fr
^{2} Université de Toulouse, UPSOMP, IRAP, 31400 Toulouse, France
^{3} CNRS, IRAP, 14 avenue Edouard Belin, 31400 Toulouse, France
^{4} Programa Explora, CONICYT, Ministerio de Educación, Moneda 1375, Santiago, Chile
Received: 1 October 2016
Accepted: 8 February 2017
Context. A precision of the order of one percent is needed on the parameters of exoplanethosts stars to correctly characterise the planets themselves. This can be achieved by asteroseismology. It is important in this context to test the influence of introducing atomic diffusion with radiative accelerations in the models upon the derived parameters. In this paper, we begin this study with the case of the star 94 Ceti A.
Aims. We aim to perform a complete asteroseismic analysis of the exoplanethost Ftype star 94 Ceti A, from the first radialvelocity observations with the High Accuracy Radial velocity Planet Searcher (HARPS) up to the final computed best models. We also aim to test the influence of atomic diffusion, including radiative accelerations, upon the computed frequencies and on the determined stellar parameters. 94 Ceti A is hot enough to suffer these effects. We also aim to test the effect of including a complete atmosphere in the stellar models.
Methods. The radial velocity observations were performed with HARPS in 2007. The low degree modes were derived and identified using classical methods and compared with the results obtained from stellar models computed with the Toulouse Geneva Evolution Code (TGEC).
Results. We obtained precise parameters for the star 94 Ceti A. We showed that including atomic diffusion with radiative accelerations can modify the age by a few percent, whereas adding a complete atmosphere does not change the results by more than one percent.
Conclusions. Atomic diffusion including radiative accelerations should be taken into account in all the computations of stellar models with masses larger than 1.3 M_{⊙} in order to determine accurate parameters for observed stars.
Key words: asteroseismology / planetary systems / diffusion / instabilities / stars: evolution
© ESO, 2017
1. Introduction
Asteroseismic studies of G and Ftype stars were strongly boosted by space observations, especially with the Kepler satellite. Big data have been obtained, leading to global studies referred to as “ensemble asteroseismology”. However, it is still important to study particular stars for deeper insights, including detailed internal physical structure. Radial velocity observations of the exoplanethost star 94 Ceti obtained with the High Accuracy Radial velocity Planet Searcher (HARPS) allowed derivation of seismic frequencies precise enough to test the importance of new stellar physics on the derived parameters. This star is on the main sequence, hot enough to suffer the effects of atomic diffusion including radiative accelerations on heavy elements in an important way. This leads to element accumulation in specific internal layers and to extra mixing (Deal et al. 2016). We thought it would be interesting to test, at least in a preliminary way, this new physics on an Ftype star through asteroseismology. We also wanted to study the impact of introducing a complete atmosphere instead of the usual grey limit.
94 Ceti is a binary system. An exoplanet has been found in orbit around 94 Ceti A, at a distance of 1.42 AU. The star 94 Ceti A is classified as an F8V star, with a visual magnitude of V = 5.08. Its parallax, π = 44.29 ± 0.28 mas (Hipparcos catalogue, van Leeuwen 2007), puts it at a distance of d = 22.38 pc. Its companion is an M dwarf that is 6.4 mag fainter at a distance of approximately 100 AU (Hale 1994).
In the present paper, we first obtain the parameters of the star 94 Ceti A (HD 19994, HR 962, HIP 14954, GJ 128) by using classical standard evolutionary models and fitting these models with the observations gathered with the spectrograph HARPS. Then we obtain parameters with more sophisticated models including detailed radiative accelerations, as computed in the Toulouse Geneva Evolution Code (TGEC), and discuss them. Finally, we add complete atmospheres to characterise the impact of more sophisticated treatments of the stellar outer layers.
2. The observational data
2.1. Spectroscopic observations
Several spectroscopic studies have been performed for 94 Ceti A. Some of them suggest that it has a nearly solar metallicity, such as Edvardsson et al. (1993) who, in their study of the chemical evolution of the galactic disk, gave it a metallicity of 0.09. Eight years later, Smith et al. (2001) gave a similar metallic value for this star in their work exclusively focused on 94 Ceti A and its abundance distribution. However most studies mainly focused on the exoplanet hoststar status of 94 Ceti and led to a stellar overmetallicity. This is the case for Santos et al. (2001, 2003, 2004), who present 94 Ceti amongst their general studies of statistical properties, metalrich nature, and spectroscopic [Fe/H] for exoplanethost stars. In our study, we use the metallicity values of their 2004 work, which presents a summary and a revision of their previous results and provides separate values for [Fe/H] derived from the highresolution spectra obtained with four different instruments: CORALIE, FEROS, UVES and UES.
Valenti & Fischer (2005) also included 94 Ceti in their spectroscopic study of cool stars. Their results are in agreement with the results of Santos et al., meaning that this star shows an overmetallicity with respect to the Sun. They also gave estimates of the age, mass, and radius of the star. Maldonado et al. (2012) spectroscopically analysed 94 Ceti as part of their survey of the metallicity of solarlike stars that have debris disks and planets. Their results are also in agreement with the previous mentioned analyses. A summary of the spectroscopic studies previously mentioned is presented in Table 1.
2.2. Interferometric observations
Direct measurements for the radius of 94 Ceti A were performed with interferometric techniques by van Belle & von Braun (2009). Using the Palomar Testbed Interferometer (PTI) and the CHARA array, they measured linear radii and effective temperatures of twelve stars with known planetary companions. For 94 Ceti A, a radius of R/R_{⊙} = 1.898 ± 0.070 was obtained, which provides an important constraint for stellar modelling. The associated effective temperature derived was T_{eff} = 6109 ± 111 K.
Summary of previous spectroscopic studies of 94 Ceti.
2.3. The exoplanet data
A Jupiterlike companion to 94 Ceti A was first announced by Queloz et al. (2000). A more precise orbital solution was later provided by Mayor et al. (2004). Discovered by using the radial velocity technique, this planet has a minimum mass of 1.68 M_{J} and orbits its parent star at a distance of 1.42 AU with a period of approximately 535.7 days. A summary of the properties of HD 19994b is given in Table 2.
Summary of HD 19994b properties (Mayor et al. 2004).
3. Seismic observations and analysis
3.1. Mode identification
We observed the star 94 Ceti A with the HARPS spectrograph mounted on the ESO 3.6 m telescope, at La Silla Observatory, Chile. We carried out the observations during seven consecutive nights in November 2007. The global Fourier transform of the observations is given in Fig. 1.
The Fourier transform was first analysed with the Period 04 code (Lenz & Breger 2005) to extract the most important peaks. However, due to aliasing it may happen that high amplitude peaks are combinations of aliases instead of real modes. For this reason we also manually investigated the Fourier transform, using the window function to pick out the real peaks and eliminate the surrounding aliases. We then used stellar models to identify the modes corresponding to the detected frequencies. This method was somewhat successful, allowing us to recognise twentysix modes in a convincing way. The frequencies of these identified modes are given in Table 3 with an uncertainty evaluated as 1 μHz. The average large separation is Δν = 64.1 μHz and the average small separation is δν_{0,2} = 4.0 μHz.
Identified Pmode frequencies (in μHz) for 94 Ceti.
Fig. 1 Global Fourier transform of 94 Ceti data. 
3.2. Rotational splitting
Given the projected rotational velocity vsini, the rotation period P_{rot}, and the linear radius R_{star}, the inclination angle i of the rotation axis on the line of sight can be determined. In the case of 94 Ceti, P_{rot} = 12.2 days (Mayor et al. (2004), from R_{HK} activity indicator) and R/R_{⊙} = 1.898 ± 0.070 from interferometric observations (van Belle & von Braun 2009). We obtain sini ~ 1.02 which means that the star is seen equatoron. Gizon & Solanki (2003) showed that in this case, for l = 1 modes, the central m = 0 component should vanish. With a rotation period of P_{rot} = 12.2 days, the shift induced by rotational splitting is ~0.95 μHz. Thus, in the 94 Ceti A case, we expect to see the m = ± 1 components of l = 1 modes separated by ~1.9 μHz. For l = 2 modes, Gizon & Solanki (2003) predict the visibility of the three components, such that m = 0 and m = ± 2. Also, the frequency shift between the visible components should be ~1.9 μHz. However, because the predicted amplitudes for the m = 0 components are lower than those of the m = ± 2 components, there may be cases where only m = ± 2 components are visible. The present observations are not precise enough to test these rotational effects, because the uncertainties on the frequency determinations are of order 1 μHz. However, this may be the reason for some of the shifts obtained between the observed and computed frequencies as seen in the échelle diagram (Fig. 3).
Because the star is seen almost equatoron, the mass of the planet HD 19994b should be close to its derived minimum value of 1.68 M_{J}. This planet may also be transiting in front of the star. To our knowledge, this sort of transit has not yet been looked for observationally.
4. Asteroseismic study
4.1. Stellar models
We used the ToulouseGeneva Evolution Code (TGEC) to compute stellar models that fitted the seismic observations of 94 Ceti A. This code performs complete computations of atomic diffusion for twentyone species, including radiative accelerations for all heavy elements. The twentyone species include twelve elements and their main isotopes: H, ^{3}He, ^{4}He, ^{6}Li, ^{7}Li, ^{9}Be, ^{10}B, ^{12}C, ^{13}C, ^{14}N, ^{15}N, ^{16}O, ^{17}O,^{18}O, ^{20}Ne, ^{22}Ne, ^{24}Mg, ^{25}Mg, ^{26}Mg, ^{40}Ca, and ^{56}Fe (Théado et al. 2012). The diffusion coefficients used in the code are those derived by Paquette et al. (1986). The Rosseland opacities are recalculated inside the model, at each time step and at every mesh point, using OPCD v3.3 and data from Seaton (2005), to take into account the local chemical composition. In this way, the stellar structure is consistently computed all along the evolutionary tracks, as well as the individual radiative accelerations of C, N, O, Ne, Mg, Ca, and Fe. This is performed by using the improved semianalytical prescription proposed by Alecian & Leblanc (2004). A more detailed discussion of these computations is provided in Deal et al. (2016).
The equation of state used in the code is the OPAL2001 equation (Rogers & Nayfonov 2002). The nuclear reaction rates are from the NACRE compilation (Angulo 1999). The mixing length formalism is used for the convective zones with a mixing length parameter of 1.8, calibrated on the Sun (Théado et al. 2009).
4.2. Determination of the best model
The stellar oscillation modes were derived for each model using the PULSE code (Brassard & Charpinet 2008). We first computed models including diffusion, without the effects of radiative accelerations, with masses ranging from 1.36 to 1.50 M_{⊙}. The initial helium mass fraction Y_{init} was chosen either as “solar” (Y_{init} = 0.271, Grevesse & Noels 1993) or “galactic” (Y_{init} = Y_{G}), the latter of which means that it increases with Z following the chemical evolution of galaxies, as derived by Izotov & Thuan 2004, 2010. The initial heavyelement massfraction Z_{init} was chosen so that the surface metallicity at the present age lies inside the error bars of the observed values. This lead to Z_{init} = 0.0245 for the low metallicity case ([Fe/H] = 0.09 ± 0.10) and to Z_{init} = 0.0305 for the high one ([Fe/H] = 0.21 ± 0.08). The larger metallicities derived by Santos et al. (2004) with CORALIE and FEROS were not considered. Evolutionary tracks computed for these four different initial compositions are presented in Fig. 2.
For each track, we determined the model with the same average large separation as the observed one and then we derived the best of all these models for each initial composition. We did this by fitting the échelle diagrams with the observed one, using χ^{2} minimisations. This selection gave four potential candidates for the best model, one for each initial chemical composition. The first one has a mass of 1.40 M_{⊙} with Y_{init} = Y_{G} and Z_{init} = 0.0245 at 2.32 Gyr (bottom right panel of Fig. 2, model 1), the second one has a mass of 1.44 M_{⊙} with Y_{init} = Y_{G} and Z_{init} = 0.0305 at 2.38 Gyr (bottom left panel of Fig. 2, model 2), the third one has a mass of 1.42 M_{⊙} with Y_{init} = 0.271 and Z_{init} = 0.0245 at 2.79 Gyr (upper right panel of Fig. 2, model 3) and the last one has a mass of 1.44 M_{⊙} with Y_{init} = 0.271 and Z_{init} = 0.0305 at 2.97 Gyr (upper left panel of Fig. 2, model 4). Comparing the χ^{2} values, we find that models 1 and 2 best fitted the observational values.
An interesting result of this study is the constraint on log g. For all initial compositions used for these model the log g value ranges between 4 and 4.1. This rules out the spectroscopic log g and T_{eff} determinations of Santos et al. (2004) obtained with the instruments CORALIE and FEROS, and the work of Valenti & Fischer (2005) and Maldonado et al. (2012). On the other hand, the results obtained by Santos et al. (2004) with UVES and UES are consistent with the seismic determination of log g .
Comparing the seismic results with the spectroscopic observations (Fig. 2), we find that the best model is model 2, with a mass of 1.44 M_{⊙}, because model 1 lies outside all spectroscopic boxes (the 1.40 M_{⊙} model on the bottom right panel of Fig. 2). The échelle diagram corresponding to this model is presented in Fig. 3. Its parameters are given in Table 4.
Fig. 2 log g vs. T_{eff} diagrams for models with four different initial chemical compositions. The masses range between 1.36 and 1.5 M_{⊙}. The error boxes come from the spectroscopic studies detailed in Sect. 2.1 and Table 1. The crosses correspond to Edvardsson et al. (1993), the pluses correspond to Smith et al. (2001), the right triangles correspond to Santos et al. (2004) (CORALIE), the triangles correspond to Santos et al. (2004) (FEROS), the dots correspond to Santos et al. (2004) (UVES), the squares correspond to Santos et al. (2004) (UES), the stars correspond to Valenti & Fischer (2005) and the left triangles correspond to Maldonado et al. (2012). The blue squares represent the models with the same mean large separation as the one determined from the seismic observations. The heavy line segments correspond to models consistent with the interferometric radius. 
Properties of 94 Ceti from this work.
Fig. 3 Echelle diagram of the 1.44 M_{⊙} at 2.38 Gyr. The frequencies computed for the models with Y_{init} = 0.297 and Z_{init} = 0.0305 are represented by blue dots (l = 0), green triangles (l = 1), red squares (l = 2), and black diamonds (l = 3). The observed frequencies are represented by crosses (error bars: 1 μHz). 
5. Impact of radiative accelerations and fingering convection on the model parameters
The radiative accelerations are larger than the gravity for stars with masses larger than 1.3 M_{⊙} at solar metallicity (Michaud et al. 1976). The best model for 94 Ceti has a mass of 1.44 M_{⊙}. Even though the initial metallicity is larger than solar, the effect of radiative accelerations may be important and has to be investigated. As shown by Richer et al. (2000) and Richard et al. (2001), and more recently by Théado et al. (2009) and Deal et al. (2016), radiative accelerations that are larger than the gravity and decrease towards the surface may lead to element accumulation in specific stellar layers. These accumulations may also trigger fingering convection if they lead to inversions of μgradients (Théado et al. 2009; Deal et al. 2016).
Our final aim is to evaluate the differences obtained for the bestmodel parameters when we include or not the radiative accelerations and the fingering convection in the computations. As a first step, we characterised the impact of radiative accelerations and fingering convection on the oscillation frequencies along evolutionary tracks. We then chose to compute the difference in ages of our best models, using the same initial conditions of mass and chemical composition.
5.1. Impact of radiative accelerations
We studied the impact of including radiative accelerations in the stellar models all along the evolutionary tracks upon the age determinations of the star. We used as initial conditions those of the best model derived in Sect. 4.2. The mean large frequency separations computed in the models were used as a proxy.
We first computed a 1.44 M_{⊙} evolutionary track, including radiative accelerations all over the main sequence. The best model that we obtained in this case has an age of 2.48 Gyr. This is roughly 4% larger than the age determined without including radiative accelerations. The most important elements in this case are iron and calcium. As we see in Fig. 4, the radiative acceleration on calcium is larger than the gravity at the bottom of the surface convective zone. This produces an accumulation of calcium and increases the surface abundance (lower panel of Fig. 5).
Furthermore, although the radiative acceleration on iron is lower than the gravity in the radiative zone (Fig. 4), it is large enough to produce a gradient in the diffusion velocity and to slow down the gravitational settling of iron between −2 < log (ΔM/M_{∗}) < −1. This leads to an iron accumulation at the bottom of the surface convective zone (lower panel of Fig. 5).
Fig. 4 Radiative acceleration on calcium (green dashed lines) and iron (blue solid lines) for the model of 1.44 M_{⊙} at 2.48 Gyr with Y_{init} = 0.297 and Z_{init} = 0.0305. The dark grey area represents the surface convective zone. 
Fig. 5 Abundance profiles of calcium (green dashed lines) and iron (blue solid lines) for the model of 1.44 M_{⊙} with Y_{init} = 0.297 and Z_{init} = 0.0305 with gravitational settling only (upper panel) and including radiative acceleration (lower panel). The dark grey areas represent the surface convective zone. 
Fig. 6 Opacity profiles for the model of 1.44 M_{⊙} with Y_{init} = 0.297 and Z_{init} = 0.0305 with gravitational settling only (blue solid curve) and including radiative acceleration (green dashed curve). The light grey area represents the surface convective zone of the model when radiative accelerations are not included. The dark grey area represents the extension of this convective zone obtained when radiative accelerations are included in the computations. 
The abundance profiles are very different when radiative accelerations are taken into account. This modifies the mass of the surface convective zone from log (ΔM/M_{∗}) = −2.56 when taking into account only gravitational settling to log (ΔM/M_{∗}) = −2.47 when including radiative accelerations. This change of the mass of the convective zone by roughly 20% induces modifications of the oscillation frequencies that partly explains the difference in age obtained between the two models. Another reason for the difference in age is the modification of the opacity profile. Figure 6 shows the opacity profiles of the two models. The modification of the abundance profiles changes the opacity in the surface convective zone and at its bottom. This is due to the accumulation of calcium in the surface convective zone and the accumulation of iron around log (ΔM/M_{∗}) ~ −2. These differences in the opacities are only of a few percent, but still have an impact on the mean density of the model and modify the age of the model for which the mean large frequency separation is consistent with the observations.
5.2. Induced fingering convection
Théado et al. (2009) and Deal et al. (2016) show that the accumulation of heavy elements in Atype stars due to the effect of radiative accelerations may produce an inverse mean molecularweight gradient and trigger fingering (thermohaline) convection. This instability deeply modifies the structure of the star and should have an impact on the oscillation frequencies. In this section we study the impact of the fingering convection on the age determination for 94 Ceti A in the same way as in Sect. 5.1. Fingering convection may occur in stars every time a local inversion of the mean molecularweight gradient appears in the presence of a stable temperature gradient. For example, this may happen in the case of heavy element accumulation or in the case of accretion of metal rich planetary matter (Vauclair 2004; Garaud 2011; Deal et al. 2013, 2015).
Fingering convection is characterised by the socalled density ratio R_{0}, which is the ratio between the thermal and μgradients: (1)The instability can only develop if this ratio is higher than one and lower than the Lewis number, which is the ratio of the thermal to the molecular diffusivities. In this case, a heavy blob of fluid falls down inside the star and continues to fall because it exchanges heat with the surroundings more quickly than it exchanges particles. If R_{0} is smaller than one the region is dynamically convective (Ledoux criterium), and if it is larger than the Lewis number the region is stable.
Various analytical treatments of fingering convection in stars were undertaken in the past, leading to mixing coefficients that differed by orders of magnitude (Ulrich 1972; Kippenhahn et al. 1980). More recently, 2D and 3D numerical simulations were performed that converged on coefficients of similar orders (Denissenkov 2010; Traxler et al. 2011).
Here we test the occurrence of fingering convection by using the recent prescription given by Brown et al. (2013) that has been confirmed by the 3D simulations of Zemskova et al. (2014). Whereas fingering convection may occur later on the evolutionary track, we found that in our best model for 94 Ceti the radiative layers are not yet subject to this instability.
The Ftype star 94 Ceti lies in the range of effective temperatures for which radiative accelerations on individual elements begin to have important effects on the stellar internal structure. In a general way, the effects are larger for larger effective temperatures and lower metallicities. Many stars observed in the future by space missions such as TESS or PLATO will be affected by these sorts of effects. The characterisation of exoplanets needs very precise stellar parameters such as the mass, the radii and the age. With the precision of these new instruments, it will be necessary to always perform complete computations of atomic diffusion, including radiative accelerations as well as the induced mixing processes, for comparison with the observations.
6. Effect of modifying the atmospheric boundary conditions
When computing acoustic frequencies in stellar models obtained with the usual evolution codes, the waves are assumed completely reflected at some layer in the atmosphere. This is a crude approximation that becomes worse and worse for larger frequencies close to the cutoff value. This has been widely discussed in the literature, and a recipe has been provided by Kjeldsen et al. (2008) to try correcting the effect.
For the present paper, we chose to compute one complete model including a realistic atmosphere and compare the results with those of the same model computed with the usual atmospheric approximations. This is not trivial; we discuss below how we proceeded.
In stellar interiors, the medium is optically thick and the radiative transfer is computed using the diffusion approximation. The outermost layers of a star being an optically thin medium, their structure requires the detailed computation of the radiative transfer. The approach we chose was to merge the internal structure of a TGEC model with a detailed ATLAS9 (Kurucz 1993) model atmosphere, to show the difference in the seismic properties obtained.
6.1. Computation of the complete model
In this section, we present the method used to merge a stellar interior with a realistic stellar atmosphere (HuiBonHoa 2013 contains further details). For a given TGEC model, we computed a model atmosphere with the same surface chemical composition and fundamental parameters. We merged the atmosphere and the interior at the layer where the temperature at the bottom of the atmosphere matches that of the interior. Once the atmosphere was anchored to the internal structure, slight adjustments were made to ensure continuity for the pressure. The density was then recomputed accordingly, assuming a perfect gas law. In the overlapping regions, the transition for each physical quantity from the interior to the atmospheric value was computed by using a third degree polynomial weighting.
The mixing length parameter α_{MLT} is usually slightly larger in the internal structure than in the atmosphere, where observational constraints can lead to values as low as 0.5 (van’t VeerMenneret & Megessier 1996). In our case, we checked that modifying the atmospheric mixing length parameter led to negligible effects on the computed frequencies.
6.2. Effect on the mean large frequency separation
We built a complete stellar model as described in Sect. 6.1, using the best TGEC model found in Sect. 5.2. This model has a mass of 1.44 M_{⊙}, an initial composition Y_{init} = 0.297 and Z_{init} = 0.0305, an effective temperature T_{eff} = 6142 K, and a gravity log g = 4.0. The average large separation found for the model without complete atmosphere is 64.1 μHz. Computing the frequencies and separations in the new model including the complete atmosphere leads to a slightly smaller average large separation of 63.8 μHz. This difference, smaller than 1 μHz, is mainly due to the more realistic location of the outer turning point. Instead of being artificially set at the outermost layer of the internal structure model, it now occurs at the place where the atmosphere becomes radiative, that is to say where the BruntVäisälä frequency rises again. The échelle diagram obtained in this case is very similar to the one obtained from the original model. We thus conclude that our results described in the previous sections above are robust.
7. Conclusion
The observations of 94 Ceti A made with HARPS allowed precise parameters for this star to be obtained, resulting in the exclusion of several spectroscopic error boxes. Our best model lies inside the error boxes provided by Santos et al. (2004) that were obtained with UVES and UES observations. On the other hand, those obtained with FEROS and CORALIE give too large a gravity and too high an effective temperature.
We have tested the effect on the final results of including atomic diffusion with radiative accelerations, as well as induced fingering convection, in the models. The star 94 Ceti A is hot enough to have important radiative effects, in spite of its high metallicity reducing them. We found that the age of the final best model is a few percent larger than the age obtained without introducing radiative accelerations in the computations. This difference should be larger for hotter stars and for stars with lower metallicities. We note that this age difference may still be larger than the one derived in the present paper. We computed radiative accelerations by using SVP tables (Alecian & Leblanc 2004) in which some elements are not yet included. Nickel, for example, may be important in that respect (Turcotte et al. 1998). Introducing a complete ATLAS atmosphere for the computation of the best model does not change the results in a significant way (less than one percent).
Considering the precision obtained with present and future space missions, our results show that complete atomic diffusion should be introduced in the models for stars with masses larger than 1.4 solar masses.
Acknowledgments
We thank the referee for his/her important remarks on the first version of this paper. The observations were performed with the HARPS spectrograph, ESO programme number 080.D0408. We also thank the French Programme National de Physique Stellaire (PNPS) of INSU/CNRS for financial support. This work was performed using HPC resources from GENCICINES/IDRIS and from LUPM.
References
 Alecian, G., & Leblanc, F. 2004, in The AStar Puzzle, eds. J. Zverko, J. Ziznovsky, S. J. Adelman, & W. W. Weiss, IAU Symp., 224, 587 [Google Scholar]
 Angulo, C. 1999, AIP Conf. Ser., 495, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Brassard, P., & Charpinet, S. 2008, Ap&SS, 316, 107 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, J. M., Garaud, P., & Stellmach, S. 2013, ApJ, 768, 34 [NASA ADS] [CrossRef] [Google Scholar]
 Deal, M., Deheuvels, S., Vauclair, G., Vauclair, S., & Wachlin, F. C. 2013, A&A, 557, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deal, M., Richard, O., & Vauclair, S. 2015, A&A, 584, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deal, M., Richard, O., & Vauclair, S. 2016, A&A, 589, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Denissenkov, P. A. 2010, ApJ, 723, 563 [NASA ADS] [CrossRef] [Google Scholar]
 Edvardsson, B., Andersen, J., Gustafsson, B., et al. 1993, A&A, 275, 101 [NASA ADS] [Google Scholar]
 Garaud, P. 2011, ApJ, 728, L30 [Google Scholar]
 Gizon, L., & Solanki, S. K. 2003, ApJ, 589, 1009 [Google Scholar]
 Grevesse, N., & Noels, A. 1993, in Origin and Evolution of the Elements, eds. N. Prantzos, E. VangioniFlam, & M. Casse, 15 [Google Scholar]
 Hale, A. 1994, AJ, 107, 306 [NASA ADS] [CrossRef] [Google Scholar]
 HuiBonHoa, A. 2013, in EAS Publ. Ser. 63, eds. G. Alecian, Y. Lebreton, O. Richard, & G. Vauclair, 35 [Google Scholar]
 Izotov, Y. I., & Thuan, T. X. 2004, ApJ, 602, 200 [CrossRef] [Google Scholar]
 Izotov, Y. I., & Thuan, T. X. 2010, ApJ, 710, L67 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
 Kippenhahn, R., Ruschenplatt, G., & Thomas, H.C. 1980, A&A, 91, 175 [NASA ADS] [Google Scholar]
 Kjeldsen, H., Bedding, T. R., & ChristensenDalsgaard, J. 2008, ApJ, 683, L175 [NASA ADS] [CrossRef] [Google Scholar]
 Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km s grid. Kurucz CDROM 13 (Cambridge, Mass.: Smithsonian Astrophysical Observatory) [Google Scholar]
 Lenz, P., & Breger, M. 2005, Commun. Asteroseismol., 146, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Maldonado, J., Eiroa, C., Villaver, E., Montesinos, B., & Mora, A. 2012, A&A, 541, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mayor, M., Udry, S., Naef, D., et al. 2004, A&A, 415, 391 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Michaud, G., Charland, Y., Vauclair, S., & Vauclair, G. 1976, ApJ, 210, 447 [Google Scholar]
 Paquette, C., Pelletier, C., Fontaine, G., & Michaud, G. 1986, ApJS, 61, 177 [NASA ADS] [CrossRef] [Google Scholar]
 Queloz, D., Mayor, M., Naef, D., et al. 2000, in From Extrasolar Planets to Cosmology: The VLT Opening Symp., eds. J. Bergeron, & A. Renzini, 548 [Google Scholar]
 Richard, O., Michaud, G., & Richer, J. 2001, ApJ, 558, 377 [NASA ADS] [CrossRef] [Google Scholar]
 Richer, J., Michaud, G., & Turcotte, S. 2000, ApJ, 529, 338 [NASA ADS] [CrossRef] [Google Scholar]
 Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064 [Google Scholar]
 Santos, N. C., Israelian, G., & Mayor, M. 2001, A&A, 373, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, N. C., Israelian, G., Mayor, M., Rebolo, R., & Udry, S. 2003, A&A, 398, 363 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, N. C., Israelian, G., & Mayor, M. 2004, A&A, 415, 1153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Seaton, M. J. 2005, MNRAS, 362, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, V. V., Cunha, K., & Lazzaro, D. 2001, AJ, 121, 3207 [NASA ADS] [CrossRef] [Google Scholar]
 Théado, S., Vauclair, S., Alecian, G., & LeBlanc, F. 2009, ApJ, 704, 1262 [NASA ADS] [CrossRef] [Google Scholar]
 Théado, S., Alecian, G., LeBlanc, F., & Vauclair, S. 2012, A&A, 546, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Traxler, A., Garaud, P., & Stellmach, S. 2011, ApJ, 728, L29 [Google Scholar]
 Turcotte, S., Richer, J., Michaud, G., Iglesias, C. A., & Rogers, F. J. 1998, ApJ, 504, 539 [NASA ADS] [CrossRef] [Google Scholar]
 Ulrich, R. K. 1972, ApJ, 172, 165 [NASA ADS] [CrossRef] [Google Scholar]
 Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 van Belle, G. T., & von Braun, K. 2009, ApJ, 694, 1085 [NASA ADS] [CrossRef] [Google Scholar]
 van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 van’t VeerMenneret, C., & Megessier, C. 1996, A&A, 309, 879 [Google Scholar]
 Vauclair, S. 2004, ApJ, 605, 874 [Google Scholar]
 Zemskova, V., Garaud, P., Deal, M., & Vauclair, S. 2014, ApJ, 795, 118 [Google Scholar]
All Tables
All Figures
Fig. 1 Global Fourier transform of 94 Ceti data. 

In the text 
Fig. 2 log g vs. T_{eff} diagrams for models with four different initial chemical compositions. The masses range between 1.36 and 1.5 M_{⊙}. The error boxes come from the spectroscopic studies detailed in Sect. 2.1 and Table 1. The crosses correspond to Edvardsson et al. (1993), the pluses correspond to Smith et al. (2001), the right triangles correspond to Santos et al. (2004) (CORALIE), the triangles correspond to Santos et al. (2004) (FEROS), the dots correspond to Santos et al. (2004) (UVES), the squares correspond to Santos et al. (2004) (UES), the stars correspond to Valenti & Fischer (2005) and the left triangles correspond to Maldonado et al. (2012). The blue squares represent the models with the same mean large separation as the one determined from the seismic observations. The heavy line segments correspond to models consistent with the interferometric radius. 

In the text 
Fig. 3 Echelle diagram of the 1.44 M_{⊙} at 2.38 Gyr. The frequencies computed for the models with Y_{init} = 0.297 and Z_{init} = 0.0305 are represented by blue dots (l = 0), green triangles (l = 1), red squares (l = 2), and black diamonds (l = 3). The observed frequencies are represented by crosses (error bars: 1 μHz). 

In the text 
Fig. 4 Radiative acceleration on calcium (green dashed lines) and iron (blue solid lines) for the model of 1.44 M_{⊙} at 2.48 Gyr with Y_{init} = 0.297 and Z_{init} = 0.0305. The dark grey area represents the surface convective zone. 

In the text 
Fig. 5 Abundance profiles of calcium (green dashed lines) and iron (blue solid lines) for the model of 1.44 M_{⊙} with Y_{init} = 0.297 and Z_{init} = 0.0305 with gravitational settling only (upper panel) and including radiative acceleration (lower panel). The dark grey areas represent the surface convective zone. 

In the text 
Fig. 6 Opacity profiles for the model of 1.44 M_{⊙} with Y_{init} = 0.297 and Z_{init} = 0.0305 with gravitational settling only (blue solid curve) and including radiative acceleration (green dashed curve). The light grey area represents the surface convective zone of the model when radiative accelerations are not included. The dark grey area represents the extension of this convective zone obtained when radiative accelerations are included in the computations. 

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.