In contrast, grid-based errors need not be strictly stabilizing so the error pattern is more scattered in cases where the errors are dominated by grid artifacts rather than SIE. This is the case for SCAN-based MBE(5) calculations using SG-2 (Fig. ). Only upon saturating the grid does SIE take over as the dominant source of error, leading to oscillatory results for MBE(n) with meta-GGA functionals, as in Fig. . However, fluctuations of almost 50 kcal/mol persist between MBE(4) and MBE(5) results, even when the high-quality EML(99,590) grid is used. In previous work, we tested several strategies to mitigate accumulation of SIE, including hybrid functionals, dielectric boundary conditions, and density-corrected (DC-)DFT. Although these procedures reduced the errors in DFT-MBE(n) calculations, none was sufficient to restore convergence, at least not for standard hybrid functionals with 20-25% exact exchange. (Functionals with 50% exact exchange did afford convergent results, indicating that SIE plays a prominent role in divergent behavior for semilocal functionals and for hybrid functionals such as B3LYP, PBE0, and SCAN0.) Because the calculations in Ref. 16 used the SG-2 quadrature grid for meta-GGA functionals, SIE-driven errors were intertwined with MBE-induced grid errors. We next revisit these calculations using the unpruned EML(99,590) grid. MBE(n) errors for F -(H 2 O) 15 clusters are plotted in Fig. for several SCAN-based calculations, all using the EML(99,590) grid. For reference, the SCAN-MBE(n) results in Fig. are the same as those in Fig. , and they are compared in Fig. alongside MBE(n) calculations using SCAN + PCM, SCAN0, and DC-SCAN, all with the same high-quality grid.