Mass-loss implementation and temperature evolution of very massive stars
Very massive stars (VMS) dominate the physics of young clusters due to their ionising radiation and extreme stellar winds. It is these winds that determine their lifepaths until expiration. Observations in the Arches cluster show that VMS all have similar temperatures. The VLT-Flames Tarantula survey analysed VMS in the 30 Dor region of the LMC also finding a narrow range of temperatures, albeit at higher values - likely a metallicity effect. Using MESA, we study the main-sequence evolution of VMS with a new mass-loss recipe that switches from optically-thin O-star winds to optically-thick Wolf-Rayet type winds through the model-independent transition mass-loss rate of Vink & Gräfener. We examine the temperature evolution of VMS with mass loss that scales with the luminosity-over-mass (L/M) ratio and the Eddington parameter ($Γ_{\rm e}$), assessing the relevance of the surface hydrogen (H) abundance which sets the number of free electrons. We present grids of VMS models at Galactic and LMC metallicity and compare our temperature predictions with empirical results. Models with a steep $Γ_{\rm e}$-dependence evolve horizontally in the Hertzsprung-Russel (HR) diagram at nearly constant luminosities, requiring a delicate and unlikely balance between envelope inflation and enhanced mass loss over the entire VMS mass range. By contrast, models with a steep L/M-dependent mass loss are shown to evolve vertically in the HR-diagram at nearly constant Teff, naturally reproducing the narrow range of observed temperatures, as well as the correct trend with metallicity. This distinct behavior of a steeply dropping luminosity is a self-regulatory mechanism that keeps temperatures constant during evolution in the HR-diagram.