Current massive star evolution grids are not able to simultaneously reproduce the empirical upper luminosity limit of red supergiants, the Humphrey-Davidson (HD) limit at high and low metallicity. In this study, we provide a better understanding of what drives massive star evolution to blue and red supergiant phases, with the ultimate aim of reproducing the HD limit at varied metallicities (Z). For solar, LMC, and SMC Z, we develop eight grids of MESA models for the mass range 20-60M to probe the effect of semiconvection and overshooting. We compare rotating and non-rotating models with efficient (alpha_semi = 100) and inefficient semi-convection (alpha_semi = 0.1), with high and low core overshooting (alpha_ov of 0.1 or 0.5). The red and blue supergiant evolutionary phases are investigated by comparing the fraction of core He-burning lifetimes spent in each phase. We find that the extension of the convective core by overshooting alpha_ov = 0.5 has an effect on the post-MS evolution which can disable semiconvection leading to more RSGs, but a lack of BSGs. We therefore implement alpha_ov = 0.1 which switches on semiconvective mixing, though for standard alpha_semi = 1, would result in an HD limit which is higher than observed at low Z. Therefore, we need to implement very efficient semiconvection of alpha_semi = 100 which reproduces the HD limit at log L ~ 5.5 for the Magellanic Clouds while simultaneously reproducing the Galactic HD limit of log L ~ 5.8 naturally. The effect of semiconvection is not active at high Z due to the depletion of the envelope structure by strong mass loss such that semiconvective regions could not form. Z-dependent mass loss plays an indirect, yet decisive role in setting the HD limit as a function of Z. For a combination of efficient semiconvection and low overshooting with standard Z-dependent mass loss, we find a natural HD limit at all metallicities.

We aim to constrain massive star evolution models using the unique testbed eclipsing binary HD166734 with new grids of MESA stellar evolution models, adopting calibrated prescriptions of overshooting, mass loss, and rotation. We introduce a novel tool: the "mass-luminosity plane" or "M-L plane", as an equivalent to the traditional HR diagram, utilising it to reproduce the testbed binary HD166734 with newly calibrated MESA stellar evolution models for single stars. We can only reproduce the Galactic binary system with an enhanced amount of core overshooting (alpha = 0.5), mass loss, and rotational mixing. We can utilise the gradient in the M-L plane to constrain the amount of mass loss to 0.5 - 1.5 times the standard Vink et al. 2001 prescriptions, and we can exclude extreme reduction or multiplication factors. The extent of the vectors in the M-L plane leads us to conclude that the amount of core overshooting is larger than is normally adopted in contemporary massive star evolution models. We furthermore conclude that rotational mixing is mandatory to get the nitrogen abundance ratios between the primary and secondary components to be correct (3:1) in our testbed binary system. Our calibrated grid of models, alongside our new M-L plane approach, present the possibility of a widened main sequence due to an increased demand for core overshooting. The increased amount of core overshooting is not only needed to explain the extended main sequence, but the enhanced overshooting is also needed to explain the location of the upper-luminosity limit of the red supergiants. Finally, the increased amount of core overshooting has -- via the compactness parameter -- implications for supernova explodibility.

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 (alpha=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 alpha = 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.