To generate the final training set, we performed the MD simulation at 500 K and selected a single structure from each trajectory. To obtain the target energy and force values, all of the selected structures were evaluated with the B3LYP-D3BJ/def2-SV(P) reference method described in Section Reference method. An equivalent procedure was used to generate the