In the medium-and short-range, Pauli repulsion interaction emerges as a consequence of electron density overlap. This term is further complicated by charge penetration and charge transfer effects, which are difficult to fit using predetermined analytical functions. Conventionally, we fit the short-range interactions to dimer calculations, assuming an isotropic pairwise additive form, such as Lennard-Jones terms completely neglect the anisotropicity and the many-body effects, they do provide a reasonable zeroth order approximation that prevents the atoms from entering unphysical regions. Empirically, adding these terms may significantly improve the stability of the ML model in numerical simulations. Therefore, we further partition the short-range interaction into the isotropic pairwise and the anisotropic many-body components: isotropic Born-Mayer function is first utilized to fit the short-range energy, and the residual of the fitting is considered to be the anisotropic many-body part, which is then tackled using ML techniques, , , , exp( )