The extragalactic background light at far-infrared wavelengths 1-3 originates from opticallyfaint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year 4 . Due to the relatively poor spatial resolution of farinfrared telescopes 5, 6 , the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations 7-10 . A previous attempt 11 at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model 12 . Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 µm. From this excess, we find that submillimetre galaxies are located in dark matter halos with a minimum mass of log[M min /M ⊙ ] = 11.5 +0.7 −0.2 at 350 µm. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe 13 , and is lower than that predicted by semi-analytical models for galaxy formation 14 .Despite recent successes in attributing most of the extragalactic background light at submillimetre wavelengths to known galaxy populations through stacking analyses 15-17 , we have not individually detected the faint galaxies that are responsible for more than 85% of the total extragalactic intensity at these wavelengths 18 . The faint star-forming galaxies are expected to trace the large-scale structure of the Universe, especially in models where galaxy formation and evolution is closely connected to dark matter halos. While not individually detected in low resolution observations, the clustering of galaxies is expected to leave a distinct signature in the total intensity variations at sub-millimetre wavelengths. The amplitude of the power spectrum of intensity vari-2 ations as a function of the angular scale provides details on the redshift distribution and the dark matter halo mass scale of dusty, star-forming galaxies in the universe 7 .For this analysis, we used data from the Herschel Multi-tiered Extra-galactic survey (HerMES 18 ), taken with the Spectral and Photometric Imaging Receiver (SPIRE 19 ) onboard the Herschel Space Observatory 20 , during the Science Demonstration Phase (SDP) of Herschel. The data are composed of a wide 218 ′ by 218 ′ area in the Lockman Hole complemented by a narrow, but very deep (30 repeated scans), map of the Great Observatories Origins Deep Survey (GOODS) North field covering 30 ′ by 30 ′ . These fields have been very well studied at other wavelengths and they are known to have a low Galactic dust density, making it easier to distinguish the extragalactic component we wish to study. The observing time to complete each of the two fields was about 13.5 hours, observing simultaneously at 250, 350, and 500 µm.To limit the influence of a few ...