Phonon-phonon interactions play a crucial role in determining lattice thermal conductivity in crystals. In first principles-based calculations, these interactions are generally taken into account up to the three-phonon scattering order. However, recent studies have highlighted the significant role of higher-order phonon-phonon interactions in accurately determining lattice thermal conductivity, even for graphene. The previous tests for the GAP model were based on static calculations of interatomic forces. To further evaluate the accuracy of the generated potentials in representing physical properties, we performed VACF-based calculations using molecular dynamics trajectories generated by the LAMMPS code. We conducted NVE simulations for all the considered crystals following equilibration NVT simulations using the potentials developed in this study.