A smooth sphere-to-cube transition is experimentally, computationally and theoretically studied in plasmonic Au nanoparticles, including retardation effects. Localized surface plasmon-polariton resonances were described with precision, discriminating among the influences of shape statistics, particle polydispersity, electrochemistry of excess (surface) charges. Sphere, cube and semicubes in between all show well-defined secular electrostatic eigenvalues, producing a wealthy of topological modes afterwards quenched by charge relaxation processes. The way both eigenvalues and plasmon wavelength vary as a function of a shape descriptor, parametrizing the transition, is explained by a minimal model based on the key concepts of crystal (or ligand) field theory (CFT), bringing for the first time to an electromagnetic analog of crystal field splitting. For any orbital angular momentum, eigenvalues evolve as in a Tanabe-Sugano correlation diagram, relying on the symmetry set by particle topology and a charge defect between cube and sphere. Expressions for non-retarded and retarded plasmon wavelengths are given and succeffully applied to both experimental UV-Vis and numerically simulated values. The CFT analogy can be promising to delve into the role of shape in nanoplasmonics and nanophotonics.Propagating and localized surface plasmon-polaritons (LSP) are coherent collective excitations by which photons couple to quasifree metal electrons. As light may be confined to a smaller scale than the photon wavelength, photonic and electronic characteristic lengths can be tuned into a nanoscale device. 1-3 Plasmonics arises from a complex interplay of electronic and geometric properties, where size and shape, 4 chemical composition, 5 surface charge 6 and dielectric environment 7 turn out to dramatically affect LSP resonance and absorption efficiency in plasmon-induced phenomena such as hot carriers injection, 8,9 radiative and resonance energy transfers. 10 However, while the relevance of topology was well established in electronic structure theory, 11 a little is known about purely topological shape effects in nanoplasmonics. The notion of shape, standing in between geometry and topology, is used therein mostly in connection with the (differential) geometry of surfaces, e.g. SP in curved media, scattering and radiation at bends and interfaces. [12][13][14][15] Topology optimization of nanostructures 16 can be regarded as well as an unrestricted shape optimization of a geometric functional. A primary shape effect in LSP resonance may be detected by the so-called figure of merit, the ratio between enhanced local and incident fields, whose formal dependence on the real and imaginary parts of the complex dielectric function is changing with the particle shape. 17 The linear optical response then may be derived from Drude's model or some semiclassical extension of it, 18,19