Stellar evolution codes, as all other numerical tools, need to be verified. One of the standard stellar objects that allow stringent tests of stellar evolution theory and models, are detached eclipsing binaries. We have used 19 such objects to test our stellar evolution code, in order to see whether standard methods and assumptions suffice to reproduce the observed global properties. In this paper we concentrate on three effects that contain a specific uncertainty: atomic diffusion as used for standard solar model calculations, overshooting from convective regions, and a simple model for the effect of stellar spots on stellar radius, which is one of the possible solutions for the radius problem of M dwarfs. We find that in general old systems need diffusion to allow for, or at least improve, an acceptable fit, and that systems with convective cores indeed need overshooting. Only one system (AI Phe) requires the absence of it for a successful fit. To match stellar radii for very low-mass stars, the spot model proved to be an effective approach, but depending on model details, requires a high percentage of the surface being covered by spots. We briefly discuss improvements needed to further reduce the freedom in modelling and to allow an even more restrictive test by using these objects.

The extent of mixed regions around convective zones is one of the biggest uncertainties in stellar evolution. One-dimensional overshooting descriptions introduce a free parameter (fov) that is, in general, not well constrained from observations. Especially in small central convective regions, the value is highly uncertain due to its tight connection to the pressure scale height. Long-term multi-dimensional hydrodynamic simulations can be used to study the size of the overshooting region as well as the involved mixing processes. Here we show how one can calibrate an overshooting parameter by performing two-dimensional Maestro simulations of zero-age-main-sequence stars ranging from 1.3 to 3.5 M⊙. The simulations cover the convective cores of the stars and a large fraction of the surrounding radiative envelope. We follow the convective flow for at least 20 convective turnover times, while the longest simulation covers 430 turnover time scales. This allows us to study how the mixing as well as the convective boundary itself evolve with time, and how the resulting entrainment can be interpreted in terms of overshooting parameters. We find that increasing the overshooting parameter fov beyond a certain value in the initial model of our simulations changes the mixing behaviour completely. This result can be used to put limits on the overshooting parameter. We find 0.010 <  fov <  0.017 to be in good agreement with our simulations of a 3.5 M⊙ mass star. We also identify a diffusive mixing component due to internal gravity waves that is active throughout the convectively stable layer, but it is most likely overestimated in our simulations. Furthermore, applying our calibration method to simulations of less massive stars suggests a need for a mass-dependent overshooting description where the mixing in terms of the pressure scale height is reduced for small convective cores.

Context. TZ Fornacis (TZ For) is an evolved detached binary system that is difficult to model and interpret, but very useful for testing stellar evolution theory and physics. Aims. We aim to search for solutions that are self-consistent and to determine the necessary stellar physics input. We also check solutions found previously for their internal consistency and for reproducibility. Methods. We use both a single and a binary stellar evolution code, and take into account all known system properties. We determine the physical stellar parameters by imposing that the models match the known radii for identical stellar ages. The evolution has to be consistent with a binary system in classical Roche geometry. Results. We obtained two different solutions to model TZ For successfully. Both depend on avoiding a long evolution on the first giant branch and imply a sufficiently large convective core on the main sequence. TZ For can be modelled consistently as a detached binary system by invoking either a substantial amount of core overshooting or a tidally enhanced wind mass loss along the red giant branch. An evolution with Roche-lobe overflow can definitely be excluded. Conclusions. A comparison of our results with previous studies also reveals that in addition to uncertainties associated with the input physics, the modelling of overshooting by different algorithms can have a strong impact.