Massive star evolution is dominated by key physical processes such as mass loss, convection and rotation, yet these effects are poorly constrained, even on the main sequence. We utilise a detached, eclipsing binary HD166734 as a testbed for single star evolution to calibrate new MESA stellar evolution grids. We introduce a novel method of comparing theoretical models with observations in the ’Mass-Luminosity Plane’, as an equivalent to the HRD (see Higgins & Vink, 2018). We reproduce stellar parameters and abundances of HD166734 with enhanced overshooting (0.5), mass loss and rotational mixing. When comparing the constraints of our testbed to the systematic grid of models we find that a higher value of 0.5 (rather than 0.1) results in a solution which is more likely to evolve to a neutron star than a black hole, due to a lower value of the compactness parameter.
Keywords. stars: mass loss, stars: evolution, stars: rotation, convection
Massive star evolution revealed in the Mass-Luminosity plane] Massive star evolution revealed in the Mass-Luminosity plane Erin R. Higgins & Jorick S. Vink] Erin R. Higgins & Jorick S. Vink
The lives of massive stars (M 8 ) are influenced by the key processes acting on their structure, including stellar winds, convection and rotation. The degree at which these processes effect evolutionary paths is dictated by the initial mass, metallicity and multiplicity of the star. In the mass range of 8-30 rotation is the dominant effect, whereas above 30 the evolution is heavily dominated by mass loss via stellar winds, e.g. Maeder & Meynet (2000). The main-sequence (MS) lifetime is dependent on the extension of the convective core by overshooting due to extra mixing of hydrogen (H) in the core, though this processes remains largely unresolved it is key in determining the final phases of evolution (e.g. final masses which form neutron stars or black holes) since the MS accounts for 90% of O star lifetimes, (Claret & Torres, 2017). Martins & Palacios (2013) present an overview of current evolutionary codes e.g. Ekström et al. (2012), Brott et al. (2011b) as well as their implementations of various physical processes. The differing treatment of rotation, mass loss and convective overshooting has led to a diversity of evolutionary masses and MS-lifetimes. Weidner & Vink (2010) discuss the effects this may have on the mass determination of O stars and further discrepancies when compared to spectroscopic masses (i.e. the ’mass discrepancy’ problem). In this paper we aim to compare new and existing prescriptions with one code, utilising the highly flexible MESA code e.g. Paxton et al. 2011. We investigate the evolution of a high mass, detached, eclipsing binary HD166734, modelled as a testbed for single star evolution with constraints on and . Dynamical masses of 39.5 and 33.5 , for the primary and secondary respectively, have allowed exploration of the dominant processes for the entire O star mass range, since in this mass regime the effects of rotation, mass loss and overshooting all play as role as they interact and overlap. We develop a novel method of analysing the evolution of our models alongside observations of HD166734 in the Mass-Luminosity plane, and present a calibrated grid of rotating and non-rotating models with a sample of Galactic O stars.
The stellar evolution code MESA has been utilised in this study in developing a set of theoretical models for calibration of physical processes such as convective overshooting, mass loss and rotation. The convective core boundary is defined by the Ledoux criterion with step-overshooting applied as the extension of the core by a factor of the pressure scale height H. Mass-loss rates are employed by the Vink et al. (2001) prescription, accounting for the bi-stability jump at 21kK and metallicity dependancies. We adopt the default metallicity of Z=0.02 in MESA to provide comparisons with a galactic observations. Rotation is applied in a fully diffusive approach with Eddington-Sweet circulation and dynamical and secular shear instabilities accounted for.
Mahy et al. (2017) provide a rare opportunity in the analysis the non-interacting binary HD166734, since the dynamical masses of both stars are in excellent agreement with their spectroscopic masses, thus allowing calibration of evolutionary masses through our models. Table 2 highlights the main stellar parameters of HD166734, including exact positions in the HRD and surface nitrogen abundances. We apply an equal-age assumption for the binary such that they have evolved from the same initial stage, allowing for constraints of the MS-width and rotation rates.
|[K]||32000 1000||30500 1000|
|log(L/ )||5.840 0.092||5.732 0.104|
|[ ]||39.5 5.4||33.5 4.6|
|[ ]||37.7 29.2||31.8 26.6|
|v sin i [ km s]||95 10||98 10|
3 Mixing and Mass loss
In order to minimise interacting processes, we firstly calculated non-rotating models which solely employ mass loss and overshooting as mixing processes, with the aim of reproducing the 3:1 ratio of surface nitrogen abundances between the primary and secondary, and HRD positions simultaneously. We explore a range of factors of the Vink et al. (2001) mass-loss prescription (0.1-10 times), alongside a variety of overshooting = 0.1-0.8, and find that all variations of this parameter space result in either negligible surface nitrogen enrichment, or a factor of 10 enrichment, representative of CN-equilibrium. We conclude that a large value of coupled with an increased mass-loss rate results in stripping of the stellar envelope exposing CNO material at the core boundary. Thus another mechanism of mixing is required to reproduce any intermediate surface enrichment, we hence explore the addition of rotation to our models. Figure 1 demonstrates the surface nitrogen abundances of both non-rotating (red) and rotating (blue) models for a range of and , we note that intermediate enrichment only occurs for rotating models. However, when analysing these models we were unable to simultaneously reproduce the stellar parameters in a HRD for both primary and secondary due to the dependance of rotation on mass loss. We found that only by removing the default implementation of rotationally-induced mass loss could we disentangle each physical process to reach a solution in the HRD.
4 The Mass-Luminosity Plane
Observations of HD166734 suggest high initial masses representative of the observed luminosities, yet present challenges when comparing this to the much lower dynamical masses. We find that an extreme factor of mass-loss rate would be required to reproduce both the observed luminosities and dynamical masses, though in this case factors 2 and above result in a significant drop in luminosity such that the observed luminosity cannot be reached.
We hence developed a method of simultaneously reproducing the mass and luminosity of our observations via a novel tool termed here as the ’mass-luminosity plane’ (M-L plane). Figure 2 illustrates an evolutionary model beginning the ZAMS at the red dot and evolving along the vector with time, or temperature as in the HRD. We find that at a given evolutionary stage or observed temperature, we may calculate the length of this vector with respect to initial mass and luminosity setting the initial position, and an observed temperature setting the final position or current evolutionary stage. We thus can use this vector length to constrain physical processes in our theoretical models to reproduce observations provided by the testbed HD166734. We observe that the length (from initial - final positions) can only be further extended by rotation or overshooting, while the gradient of the vector relies on the multiplication factor of the mass-loss rate.
We exclude extreme factors of the Vink et al. (2001) mass-loss rate since a factor of 0.1 results in an almost vertical vector which cannot reproduce dynamical masses, as well as factors beyond 1.5 times the Vink et al. (2001) prescription since this results in a much too shallow gradient to reproduce observed luminosities, leaving an accepted parameter range of 0.5-1.5 time the Vink et al. (2001) prescription. When constraining our models in the M-L plane with observations of HD166734 we found that enhanced overshooting was required to reproduce stellar parameters simultaneously (=0.3 0.1 and 0.5 0.1 for the primary and secondary respectively), since initial rotation rates of 250 km s and 120 km s for the primary and secondary respectively were also calibrated via surface nitrogen enrichments, and extra mixing was still required.
5 Grid analysis
We consolidate our results for HD166734 by comparing our constraints to a galactic sample of 30 O stars from Markova et al. (2018). We calculated a systematic grid of models for initial masses 8-60 , initial rotation rates of 0-500 km s and for two extreme values of 0.1 and 0.5, see figure 3 for example. We identify the main calibrations from this study as extended overshooting to 0.5, and reduced mass loss with rotation by excluding rotationally-induced mass loss. We note that comparing models of 0.1 and 0.5, we find that the luminosity cut-off for red supergiant evolution is much higher ( log L/ 6.0) for 0.1, than for 0.5 ( log L/ 5.5-5.8), regardless of rotation rate or mass-loss factor. We also note that enlarging the overshooting region may have consequences for the compactness parameter () in the final phases of our models. We find a lower value of with an enhanced overshooting of 0.5 (when compared to 0.1), thus making it easier to create neutron stars and more challenging to evolve to black holes.
6 Discussion and Conclusion
We present a novel method of comparing evolutionary tracks with observations in the Mass-Luminosity plane (Higgins & Vink, 2018), demonstrated by our constraints of the testbed, detached binary HD166734. We overlay these calibrations with a systematic grid of models alongside a galactic sample of O stars, which reveals dependancies of and on the final stages of evolution. We emphasise the necessity of rotational mixing in stellar models for reproducing observed intermediate surface enrichments.
We conclude that extra mixing by overshooting is required to simultaneously reproduce the observed luminosities and dynamical masses of HD166734 in the M-L plane, while we exclude extreme factors of the mass-loss rate that lie beyond 0.5-1.5 times the Vink et al. (2001) rates. Finally, we omit the application of rotationally-induced mass loss since stellar parameters of HD166734 cannot be simultaneously reproduced with this theory.
- Brott et al. (2011b) Brott, I., de Mink, S. E., Cantiello, M., Langer, N., de Koter, A., Evans, C. J., Hunter, I., Trundle, C., and Vink, J. S. 2011, Astronomy & Astrophysics, 530, A115.
- (Claret & Torres, 2017) Claret, A., & Torres, G. 2017, Astrophysical Journal, 849(1), 18.
- Ekström et al. (2012) Ekström, S., Georgy, C., Eggenberger, P., Meynet, G., Mowlavi, N., Wyttenbach, A., Granada, A., Decressin, T., Hirschi, R., Frischknecht, U. and Charbonnel, C. 2012, Astronomy & Astrophysics, 537, A146.
- Maeder & Meynet (2000) Maeder, A., & Meynet, G. 2000, Annual Review of Astronomy & Astrophysics, 38(1), 143-190.
- Mahy et al. (2017) Mahy, L., Damerdji, Y., Gosset, E., Nitschelm, C., Eenens, P., Sana, H. and Klotz, A. 2017, Astronomy & Astrophysics, 607, A96.
- Markova et al. (2018) Markova, N., Puls, J., & Langer, N. 2018, Astronomy & Astrophysics, 613, A12.
- Martins & Palacios (2013) Martins, F., & Palacios, A. 2013, Astronomy & Astrophysics, 560, A16.
- Paxton et al. 2011 Paxton, B., Bildsten, L., Dotter, A., Herwig, F., Lesaffre, P. and Timmes, F. 2011, Astrophysical Journal Supplement Series, 192(1), 3.
- Vink et al. (2001) Vink, J.S., de Koter, A. and Lamers, H.J.G.L.M. 2001, Astronomy & Astrophysics, 369(2), 574-588.
- Weidner & Vink (2010) Weidner, C., & Vink, J. S. 2010, Astronomy & Astrophysics, 524, A98.