The most straightforward approach involves the application of two-body interaction potentials such as Lennard-Jones in coarse-grained models. Models employing this strategy include, for example, the Kremer-Grest (KG) polymer model 1 and the Martini model for polymer physics and biological systems, respectively. Despite their strengths in handling complex systems, these models face significant computational demands compared to other mesoscopic alternatives. The Lennard-Jones potential used in these models, which includes a strong repulsive component (e.g., 1/r ), restricts the time steps to maintain numerical stability, and the more intricate potential energy surface compared to other methods may lead to kinetic trapping.