Example of the impact of SFR burstiness on the MZR for a system with M_dm = 3 × 10¹⁰ M_⊙ at z = 6 evolved with our minimal physical model. Different lines indicate different delay times for the SN feedback, as is indicated in the legend. Note that the model results have been resampled in 0.2 Myr steps and we do not consider the evolutionary phases when the galaxies have a negligible amount of gas (M_g ≤ 10⁻³ M_dm), which allows us to have a well-defined Z_g even in periods when no gas and metals are present in the galaxy (see Fig. 2). For the effect of a coarser time sampling, see Appendix B, in particular Fig. B.1.