2025MNRAS.539..347V

The Herschel-SPIRE Dark Field is the deepest field produced by the SPIRE instrument pushing down below the galaxy confusion limit in each of the 250, 350, 500 $\mu$m bands. Standard source extraction techniques inevitably fail because of this, and we must turn to statistical methods. Here, we present a P(D) - probability of deflection - analysis of a 12$^{\prime }$ diameter region of uniform coverage at the centre of the Herschel-SPIRE Dark Field. Comparing the distribution of pixel fluxes from our observations to the distributions predicted by current literature models, we find that none of the most recent models can accurately recreate our observations. Using a P(D) analysis, we produce a fitted differential source count spline with a bump in the source counts at faint flux densities, followed by a turnover at fainter fluxes, required to fit the observations. This indicates a possible missing component from the current literature models that could be interpreted perhaps as a new population of galaxies, or a missing aspect of galaxy evolution. Taking our best-fitting results, we also calculate the contribution to the cosmic infrared background (CIB) in each of the bands, which all agree with the Planck CIB measurements in this field.