As an alternative to BH&H-LYP, one might consider using newer meta-GGA functionals. We consider the semilocal SCAN functional, for which it has been claimed that many-body SIE is small in neutral water clusters. For F -(H 2 O) 15 , however, SCAN exhibits divergent behavior in MBE(n) calculations (see Fig. ). The same is true for the hybrid SCAN0 functional, 130 which uses α hfx = 0.25 (Fig. ). We also tried the ωB97X-V functional that sets α hfx = 0.167 for short-range exchange and is a very good all-around density functional, but it also exhibits serious oscillations for F -(H 2 O) 15 as shown in Fig. . Each of these meta-GGAs is inferior to BH&H-LYP for the SIE-dominated "SIE4×4" subset of GMTKN55, 131 which suggests that these meta-GGAs exhibit larger all-around SIE as compared to BH&H-LYP.