Figure (d) illustrates the force MAE distributions predicted for default (0NN), 1NN, and 1+2NN models, each trained on 100, 300, and 500 structures, respectively. The force MAE distributions for graphene with a pyridinic-N defect display a single Gaussian-like peak across all dataset sizes. For 100 training points, the default model yields a force MAE distribution peaking around 0.40 kcal/(mol•Å) with RMSE of 0.59 kcal/(mol•Å). In contrast, the 1NN and 1+2NN models produce noticeably narrower distributions that peak near 0.38 and 0.30 kcal/(mol•Å), with corresponding RMSE values of 0.57 and 0.50 kcal/(mol•Å), respectively.