The density derived from NNMD was 2.441 g/cm 3 , consistent with the experimental result (2.457 g/cm ) with an error of -0.64%. ), consistent with the value reported by the previous XAFS 127 study of ISG-1 (2.09 Å). The Zr 5 polyhedra, as shown in Fig. right, showed a square pyramidal structure rather than the bipyramidal shape expected on the basis of the hybrid orbital of the five-coordinated structure. This may be related to the presence of a lone pair in the axial direction. Two types of edge-sharing bonds of Zr and Si have been found in zircon and zirconia crystals. Edge-sharing rings composed of a combination of Al and Si are found in aluminosilicate crystals. The formation of small rings containing edge-sharing bonds appeared to be due to the contribution of Al and Zr coordination structures, analogous to these reported crystalline structures. The protocol proposed here to construct NNP using pre-training is applicable to other amorphous multi-component systems. This could expand our understanding the structures of rare elements in glasses, as well as particular structures such as surfaces and defects. The NNP and datasets used in this study are publicly available in Zenodo repository. This is in contrast to commercial NNMD software using black-box NNP such as the Matlantis package. Therefore, all researchers can reproduce the simulations at no cost. With NNMD and this open data repository, further pioneering studies involving previously impossible structural modeling of novel glasses can be expected.