Abstract We developed new parameterizations of local regularized finite-range pseudopotentials up to next-to-next-to-next-to-leading order (N³LO), used as generators of nuclear density functionals. When supplemented with zero-range spin–orbit and density-dependent terms, they provide a correct single-reference description of binding energies and radii of spherical and deformed nuclei. We compared the obtained results to experimental data and discussed benchmarks against the standard well-established Gogny D1S functional.