An important issue that our models do not address is the role of fluoride anions and organic structure-directing agents (OSDAs). Earlier works , which consider only a limited amount of configurations, have shown that the role of fluoride/OSDA is significant and that their inclusion in the model commonly leads to an increased preference for germanium occupation of the D4R units and to the lowering of the asymmetry in Ge occupations among the D4R units present in the unit cell. Gramatikov et al. also showed that the stability of the configurations depends heavily on fluoride location and OSDA orientation, which could not be extensively sampled using the costly periodic DFT calculations. In addition, one can expect that the size and the chemical composition of the OSDA as well as the zeolite topology would play a sizable role as well. Hence, it appears that, until a more general ML potential, similar to the examples proposed recently , capable of comprehensive sampling of configuration space that also encompasses fluoride/OSDAs is available, it is difficult to draw reliable, quantitative predictions regarding the germanium distributions for any of the (D4R-containing) zeolite frameworks. Also, during zeolite synthesis, besides the thermodynamic stability, kinetic factors may play a role, which are not captured by Monte Carlo simulations focusing on low-energy structures. Despite these limitations, the qualitative agreement of our calculations with some of the experimental observations (see Section 3.5 below) suggest that even by employing only the bare-framework models, one can capture an important aspect of germanium distributions, which may turn out to be decisive for some of the topologies.