In the present work we detail how the many-body potential energy landscape of elemental carbon, as described by interatomic potentials, can be explored by utilising the nested sampling algorithm, allowing the calculation of their pressure-temperature phase diagram up to high pressures. We present a comparison of four interatomic potential models: Tersoff, EDIP, GAP-20 and its recently updated version, GAP-20U. Our evaluation is focused on their macroscopic properties, particularly on their melting transition and on identifying thermodynamically
stable solid structures up to at least 100 GPa. The studied models all form graphite structures upon freezing at lower pressure, then the diamond structure as the pressure increases. In the cases of the Tersoff, EDIP and GAP-20U potentials we were able to locate the transition between these phases. We placed particular focus on the state-of-the-art machine-learned models, the GAP-20 and GAP-20U, and calculated their phase diagrams up to 1 and 0.1 TPa, respectively, to evaluate their predictive capabilities well outside of the models’ fitting range. The phase diagrams of the GAP models showed remarkably good agreement with the experimental phase diagram up to 200 GPa. However, we found that the accuracy of the models’ descriptions of graphite was adversely affected by the relative stability of phases with incorrect layer spacing. By adding a suitable selection of graphite structures to the GAP training set and re-training the model, we have derived an improved model — referred to as the GAP-20U+gr — that suppresses erroneous local minima in the graphitic energy landscape. At extreme high pressure nested sampling identified two novel stable solid structures in the GAP-20 model, a strained diamond structure above 270 GPa and a strained hexagonal-close-packed structure above 890 GPa. However, the stability of these two phases were not confirmed by DFT calculations, highlighting potential routes to further improve the GAP model.