Metal nanoparticles (NPs) are ubiquitous in many fields, from nanotechnology to heterogeneous catalysis, with properties differing from those of single-crystal surfaces and bulks. A key aspect is the size-dependent evolution of NP properties toward the bulk limit, including the adoption of different NP shapes, which may bias the NP stability based on the NP size. Herein, the stability of different Pd_n NPs (n = 10–1504 atoms) considering a myriad of shapes is investigated by first-principles energy optimisation, leading to the determination that icosahedron shapes are the most stable up to a size of ca. 4 nm. In NPs larger than that size, truncated octahedron shapes become more stable, yet a presence of larger {001} facets than the Wulff construction is forecasted due to their increased stability, compared with (001) single-crystal surfaces, and the lower stability of {111} facets, compared with (111) single-crystal surfaces. The NP cohesive energy breakdown in terms of coordination numbers is found to be an excellent quantitative tool of the stability assessment, with mean absolute errors of solely 0.01 eV·atom⁻¹, while a geometry breakdown allows only for a qualitative stability screening.