While our findings support that a three-body potential is enough to reproduce previously discovered phases, higher-body interactions have not been investigated and their role in NP assembly thus remains unclear. However, the fact that we could reproduce all known phases of polymer-grafted NPs using a three-body potential implies the role of higher-body interactions may not be significant. Nonetheless, the accuracy of the developed potential could be further improved by including higher-body interactions. However, this would require computation of high dimensional PMFs, along 𝑛(𝑛 -1)/2 distinct interparticle distances for each 𝑛-body contribution, which can quickly become prohibitive with increasing 𝑛. In such cases, deep neural networks may offer a better alternative for fitting such higher-body contributions.