For the optimization and density of states calculations, the Monkhorst-Pack grids of 9×9×9 and 27×27×27 were utilized, respectively, to accurately sample the first Brillouin zone. This ensured reliable integration over the Brillouin zone for the calculations. Furthermore, to account for the weak vdW interactions, vdW correction methods were incorporated into the calculations.