One natural way to study the nonbonding interaction is through energy decomposition analysis (EDA) methods, such as Symmetry-Adapted Perturbation Theory (SAPT). In SAPT, the total nonbonding energy can be decomposed into different physically meaningful components including: electrostatic, induction/polarization, dispersion, and exchange repulsion etc. In the long-range region where the electron density overlap between two molecules is negligible, most of these components can be further distributed and expanded using population analysis and multipole expansion techniques. This procedure eventually leads to a series of well-converged force field terms with simple analytical expressions, the parameters of which can be directly computed at monomer level. Using these techniques, we developed accurate and transferrable intermolecular force fields for a wide range of systems including: small molecular fluids (such as CO2 and N2), ionic liquids, and bulk materials such as metal organic frameworks. Meanwhile, in previous studies, the short-and medium-range interactions were usually fitted at dimer level using predetermined functional forms (i.e., Lennard-Jones, Born-Mayer, or Slater type functions etc.). While these functions are physically motivated, they are derived with heavy approximations and can show significant errors in strong interacting systems. In short-range, the non-negligible electron density overlaps cause the failure of multipole expansion and give rise to strong charge penetration and charge transfer effects. Complicated by these effects, the short-range interactions can be quite anisotropic and many-body in nature, and its exact functional form remains unknown.