For two different 2D materials, their unit cells are mostly incommensurate (that is, they have different values for the a lattice parameter). Therefore, we seek to construct the bilayers such that they are approximately commensurate; that is, for each of the two monolayers, we increase the number of monolayer images for each monolayer until the difference in the lattice constant of each monolayer supercell is < 2%. We wrote a Java program that automates the search for approximately commensurate supercells where the number of images of either monolayer does not exceed 5. For simplicity, the twist angle adopted (angle by which one monolayer lattice rotates with respect to the other) is 0. The number of bilayers we choose from this set is 267. These bilayers will be used for the prediction of the interlayer distance.