On the nature of the core of α Centauri A: the impact of the metallicity mixture

Forward asteroseismic modelling plays an important role towards a complete understanding of the physics taking place in deep stellar interiors. With a dynamical mass in the range over which models develop convective cores while in the main sequence, the solar-like oscillator α Centauri A presents itself as an interesting case study. We address the impact of varying the metallicity mixture on the determination of the energy transport process at work in the core of α Centauri A. We ﬁnd that (cid:38) 70 % of models reproducing the revised dynamical mass of α Centauri A have convective cores, regardless of the metallicity mixture adopted. This is consistent with the ﬁndings of Nsamba et al., where nuclear reaction rates were varied instead. Given these results, we propose that α Centauri A be adopted in the calibration of stellar model parameters when modelling solar-like stars with convective cores.


INTRODUCTION
Most parameters used in stellar modelling are calibrated based on the Sun, e.g., the mixing length parameter, the helium-to-heavy-element ratio, chemical abundances etc. This approach is well-suited to the modelling of stars with similar properties to the Sun, i.e., for solar-type stars with a mass below 1.1 M . The quest for a more massive star with well-known properties and interior structure is of the utmost importance, as such star could become a potential model calibrator for solar-like stars having convection as the main energy transport process in their cores.
The bright, nearby binary α Centauri is amongst the best characterized star systems, with a plethora of available high-precision observations, e.g., parallaxes, angular diameters, interferometric radii, metallicities, effective temperatures, luminosities, and oscillation frequencies (Söderhjelm 1999;Pourbaix et al. 2002;Bouchy and Carrier 2002;Bedding et al. 2005;Bazot et al. 2007;de Meulenaer et al. 2010;Pourbaix and Boffin 2016;Kervella et al. 2016Kervella et al. , 2017. α Centauri A is of particular interest, since its dynamical mass is in the range (1.1 -1.15 M ; Aerts et al. 2010) over which models constructed at solar metallicity are expected to develop convective cores while in the main sequence. This has given rise to studies that aimed at establishing the nature of its core and at exploring the physics that affect core properties.
Forward stellar modelling, when coupled with observational constraints from asteroseismology, constitutes a valuable tool in the exploration of the physics of stellar interiors (e.g., Miglio and Montalbán 2005;Mathur et al. 2012;Metcalfe et al. 2012Metcalfe et al. , 2014Lebreton and Goupil 2014;Silva Aguirre et al. 2015;Nsamba et al. 2018b). Bazot et al. (2016) complemented asteroseismic and spectroscopic observables with the dynamical mass determined by Pourbaix et al. (2002) to tightly constrain stellar models of α Centauri A. They explored the impact of varying the nuclear reaction rates, overshoot, and diffusion of hydrogen and helium, having found a noticeable change in the number of best-fit models 1 with convective cores when nuclear reaction rates were varied. Following the revision of the dynamical mass of α Centauri A (Pourbaix and Boffin 2016), Nsamba et al. (2018a) (hereafter Paper I) carried out a detailed modelling of this star again allowing the nuclear reaction rates to vary, and found about 70% of best-fit models to have convective cores. More recently, Joyce and Chaboyer (2018) suggested that, if α Centauri A has a convective core, then it would be necessary to modify standard physical prescriptions (e.g., enhancing diffusion) in order to correctly model the star. Amongst the different model physics explored in Paper I and in Bazot et al. (2016), the impact of the metallicity mixture on the core properties of α Centauri A has, however, not been investigated.
In this work we investigate the impact of the metallicity mixture on the inferred nature of the core of α Centauri A. The paper is organized as follows: Section 2 describes the stellar model grids and sets of observables used in the optimization process. A discussion of the results is presented in Sect. 3. Section 4 presents the conclusions.

MODEL GRIDS AND OBSERVATIONAL CONSTRAINTS
To explore the impact of varying the metallicity mixture on the nature of the core of α Centauri A, we set up two grids (A and B) with the same model physics except for the metallicity mixture (see Table 1). The stellar evolution code MESA (Modules for Experiments in Stellar Astrophysics; Paxton et al. 2011Paxton et al. , 2013Paxton et al. , 2015Paxton et al. , 2018 version 9793 was used to generate the grids. We set the metallicity mixture in Grid A according to Grevesse et al. (1998) with a solar surface heavy element mass fraction Z surface, = 0.016, while Grid B uses the metallicity mixture from Asplund et al. (2009) with Z surface, = 0.0134. The main motivation for considering these two mixtures goes back to the theoretically predicted sound speed profiles for solar models constructed with the different composition mixtures. Delahaye and Pinsonneault (2006) reported that solar models using the Grevesse et al. (1998) and Asplund et al. (2009) mixtures yield a sound speed profile close to that of the real Sun, as opposed to models that use the Asplund et al. (2005) mixture. We recall that the metallicity, [Fe/H], is defined as: where X surface is the surface hydrogen mass fraction. Notice that varying the metallicity mixture requires setting the corresponding appropriate opacities. We used opacities from OPAL tables (Iglesias and Rogers 1996) at high temperatures, whereas at low temperatures tables from Ferguson et al. (2005) were used instead, for the respective metallicity mixtures. We employed the Joint Institute for Nuclear Astrophysics Reaction Library (JINA REACLIB; Cyburt et al. 2010) in both grids. The 14 N(p, γ) 15 O and 12 C(α, γ) 16 O rates were described according to Imbriani et al. (2005) and Kunz et al. (2002), respectively. We note that the 12 C(α, γ) 16 O reaction rate is less relevant for stars on the main-sequence phase but is vital in more evolved stars, i.e., stars at the core helium-burning evolution stage. Table 1 lists the macrophysics used in either grid. We note that core overshoot was included as described by Herwig (2000) for models identified to have developed convective cores. Atomic diffusion was included in all our models according to Thoul et al. (1994). The latter is known to be a vital chemical transport process in low-mass stars (e.g., Silva Aguirre et al. 2017; Nsamba et al. 2018b). The mixing length theory, as described by Böhm-Vitense (1958), was used to describe convection. We also implemented the Grey-Eddington atmosphere to integrate the atmospheric structure from the photosphere to an optical depth of 10 −4 . The helium mass fraction (Y ) was estimated using the relation where ∆Y /∆Z is the helium-to-heavy-element relation (set to 2.0; Serenelli and Basu 2010) and Y 0 is the big bang nucleosynthesis value (set to 0.2484; Cyburt et al. 2003).
Evolutionary tracks were evolved from the zero-age main sequence to the end of the subgiant evolution phase. The terminal criterion affecting the tracks is similar to that implemented in Paper I. The parameter space of the model grids is as follows: M ∈ [1.0, 1.2] M in steps of 0.01 M ; mixing length parameter, α mlt ∈ [1.3, 2.5] in steps of 0.1; overshoot parameter, f ov ∈ [0, 0.03] in steps of 0.005; and Z ∈ [0.023, 0.039] in steps of 0.005. Each model grid contains about 156,000 models. The corresponding adiabatic oscillation frequencies of each model, for spherical degrees l = 0, 1, 2, and 3, were determined using GYRE (Townsend and Teitler, 2013). The surface effects were accounted for using the combined-term surface correction method described by Ball and Gizon (2014). This surface correction method has been reported to yield the least internal systematic uncertainties among the different available corrections (e.g., Nsamba et al. 2018b;Compton et al. 2018;Jørgensen et al. 2019). with the interferometric radius from Pourbaix and Boffin (2016). Run II adopts T eff and interferometric radius values from Kervella et al. (2016). We further note that Pourbaix et al. (2002) derived a dynamical mass of 1.105 ± 0.0070 M . This dynamical mass was then revised by Pourbaix and Boffin (2016), who obtained a value of 1.133 ± 0.0050 M .
Finally, the same asteroseismic constraints as in Paper I were adopted. The Bayesian code AIMS (Asteroseismic Inference on a Massive Scale; Lund and Reese 2018; Rendle et al. 2019), a software for fitting stellar pulsation data, was used as our optimization tool. Stellar parameters and their associated uncertainties were taken as the mean and standard deviation of the resulting posterior probability distribution functions (PDFs), as output by AIMS. Table 3 presents the stellar parameters determined using the model grids described in Table 1 and the sets of observables in Table 2. Our results show that we are able to reproduce the dynamical masses of Pourbaix and Boffin (2016) (Run I) and Pourbaix et al. (2002) (Run II) within 1σ (see Fig. 1 and Table 3). We note that the observed luminosity (i.e., 1.521 ± 0.015 L ; Kervella et al. 2017) of α Centauri A was not included among the sets of observables as shown in Table 2. This is constrained via the combination of the interferometric radius and effective temperature. Our derived luminosity values are in agreement with the observed values (see Table 4). Run II and Run I luminosity values from both grids are consistent within 1σ and 2σ, respectively. The slight increase in luminosity values obtained in Run I is attributed to the larger interferometric radius used (see Table 2).  (2016) and Pourbaix et al. (2002) are shown in orange and grey, respectively. A clear difference can be seen in the right panel of Fig. 1 (Run II) between the stellar mass posterior PDFs obtained using the two grids. Since varying the metallicity mixture has been shown to have a   Nsamba et al. 2018b), this feature can instead be explained by the different core properties of the best-fit models. Table 3 shows that the number of best-fit models with convective cores changes from 46% to 79% when the metallicity mixture is changed from that of Grevesse et al. (1998) to that of Asplund et al. (2009). This happens since the dynamical mass of Pourbaix et al. (2002) lies within a range in which the onset of the CNO (carbon-nitrogen-oxygen) cycle, and thus core convection, is highly sensitive to the adopted metallicity mixture.

DISCUSSION
A different scenario is found when considering models that reproduce the dynamical mass of Pourbaix and Boffin (2016), with the stellar mass posterior PDFs showing excellent agreement (Run I; see left panel of Fig. 1). The percentage of models with convective cores is now consistent (i.e., 70%) irrespective of the model grid adopted. We note that best-fit models are on average higher in mass compared to Run II and most have already developed convective cores, with any variation in the metallicity mixture generating no significant difference on their core properties.
In Table 3, we show the upper limits of the convective-core mass (M c ) of our best-fit models. From all runs, we find the core radius to have an upper limit of 0.11 R . Bazot et al. (2012) derived an upper limit for the radius and the mass of a possible convective core in α Centauri A to be 0.059 R and 0.035 M , respectively. These limits were derived while taking into account the small frequency separation (δν) in the optimisation, as this parameter can provide a direct estimation of the convective core characteristics. Furthermore, when exploring the contribution of the different model physics to the nature of the core of α Centauri A, Bazot et al. (2016), in their table 3, report the core radius of their best-fit models to vary between 0.026 and 0.084 R , which is consistent with our findings.
The top left panel of Fig. 2 shows that the derived stellar ages are in excellent agreement irrespective of the grid and observational constraints used. Furthermore, these ages are consistent with literature values (Kim 1999;Yıldız 2007;Bazot et al. 2016;Nsamba et al. 2018a;Joyce and Chaboyer 2018). Table 3 and the top right panel of Fig. 2 show that the α mlt estimated based on either Run I (solid lines) or Run II (dashed lines) are consistent within 1σ. The values of α mlt across runs are however different, this being mainly due to the different radius constraints used (see Table 2). We note that the interferometric radius measurements used in each run indirectly constrain the model mass. Moreover, α mlt is known to have a significant degree of correlation with the stellar mass and effective temperature (Pinheiro and Fernandes, 2013).
A clear contrast can be seen in the bottom panels of Fig. 2 between the best-fit models obtained using the two grids. As expected, the grid based on the metallicity mixture from Asplund et al. (2009) (Grid B) leads to best-fit models with a lower Z compared to those based on the mixture from Grevesse et al. (1998) (Grid A). A similar feature can be seen for the surface helium mass fraction, Y surf .
The model properties that influence the onset of the CNO cycle (and associated convective core) include the adopted physics, metallicity, and mass. As mentioned in Sect. 2, both grids contain the same physics apart from the metallicity mixture. We note that models with high metallicity have a higher chance of developing convective cores. This is because a high metallicity leads to an increase in opacity, which in turn reduces the efficiency of radiative energy transport. This ultimately results in an increase in core temperature which favours the onset of the CNO cycle. Similarly, models with a higher mass have higher core temperatures, hence higher chances of developing a convective core. The top left panel of Fig. 3 shows that best-fit models with high mass develop large and massive convective cores. This could be explained by their high overshoot parameter values as shown in the bottom panel of Fig 3. Best-fit models with f ov = 0 have smaller core masses and radii. The top right panel of Fig. 3 shows no clear trend regarding the contribution from the initial metal mass fraction. We note that consistent results as those in Fig. 3 are obtained when different runs are used.
It is interesting to assess the dominant model property that facilitates the occurrence of convective cores for the best-fit models in either run. Despite the high metallicity of best-fit models from Grid A (Run II) (see bottom left panel of Fig. 2), the majority of these models have masses 1.1 M (see right panel of Fig. 1) resulting into 46% of models with convective cores. Grid B (Run II) contains most of the best-fit models with low metallicity but with masses 1.1 M (see right panel of Fig. 1), leading to 79% of models with convective cores. Hence, for Run II (both grids), model mass is the dominant model property responsible for the onset of the CNO cycle. For Grid A (Run I), it is challenging to determine the dominant model property that yields convective cores. This is because most of the best-fit models have masses 1.1 M and high metallicity (see left panel of Fig. 1 and bottom left panel of Fig. 2). However, for Grid B (Run I), the majority of best-fit models have low metallicity but masses 1.1 M , with 79% of models having convective cores. Therefore, also in this case model mass is the dominant model property contributing to the onset of the CNO cycle. Further, our results show that the mass range over which models constructed at different metallicities are expected to develop convective cores is 1.05 -1.15 M .

CONCLUSIONS
In Paper I we assessed the impact on the nature of the core of α Centauri A of varying the nuclear reaction rates, which showed that 70% of best-fit models reproducing the revised dynamical mass of Pourbaix and Boffin (2016) have convective cores. In this article, we expanded on the previous work by exploring the impact of varying the metallicity mixture (and corresponding opacities). Our findings show that 70% of best-fit models reproducing the revised dynamical mass have convective cores.
In sum, the percentage of best-fit models with convective cores remains above 70% when imposing the most up-to-date set of observational constraints. This happens irrespective of the adopted metallicity mixture and nuclear reaction rates. Therefore, we propose that α Centauri A be adopted in the calibration of stellar model parameters when modelling solar-like stars with small convective cores.
During the "standard" solar calibration process, the initial metal mass fraction (Z), initial helium mass fraction (Y ), and the mixing length parameter (α mlt ) are varied until a model is attained that satisfies the observed oscillation frequencies, effective temperature, metallicity, luminosity, and radius at the current solar age. The same model physics and solar calibrated parameters are then used to create grids for modelling other stars. Unlike the case of the Sun, there is no model-independent age for α Centauri A, but we do have a precise dynamical mass, interferometric radius, effective temperature, metallicity, and luminosity. In addition, we currently have ground-based seismic data, with the quality of those data expected to improve following the star's planned observations by space-based missions such as NASA's Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015;Campante et al. 2016) and ESA's PLAnetary Transits and Oscillations of stars (PLATO; Rauer et al. 2014). This will improve the precision of observed oscillation frequencies and is also expected to increase the number of oscillation frequencies for all observable spherical degrees (l). This will support a more comprehensive asteroseismic analysis than the one presented here and in Paper I. Therefore, with all these sets of observables, it will be possible to carry out a calibration procedure similar to the "standard" solar calibration routine briefly described above, without having the age among the constraints. It will also be feasible to provide effective constraints on some aspects of the physics, namely convection (mixing length, overshoot, surface effects), diffusion and opacities. The potential for constraining reaction rates is also a possibility when two or all stellar components of this triple system have seismic data of high precision available.