We present a measurement of the angular bispectrum of the millimeter-wave sky in observing bands centered at roughly 95, 150, and 220 GHz, on angular scales of 1' 2 of three-band South Pole Telescope data, and we use a multi-frequency fitting procedure to separate the bispectrum of the tSZ effect from the extragalactic source contribution. We simultaneously detect the bispectrum of the tSZ effect at >10s, the unclustered component of the extragalactic source bispectrum at >5s in each frequency band, and the bispectrum due to the clustering of DSFGs—i.e., the clustered cosmic infrared background (CIB) bispectrum—at >5s. This is the first reported detection of the clustered CIB bispectrum. We use the measured tSZ bispectrum amplitude, compared to model predictions, to constrain the normalization of the matter power spectrum to be s8 = 0.787 ± 0.031 and to predict the amplitude of the tSZ power spectrum at l = 3000. This prediction improves our ability to separate the thermal and kinematic contributions to the total SZ power spectrum. The addition of bispectrum data improves our constraint on the tSZ power spectrum amplitude by a factor of two compared to power spectrum measurements alone and demonstrates a preference for a nonzero kinematic SZ (kSZ) power spectrum, with a derived constraint on the kSZ amplitude at l = 3000 of A kSZ = 2.9 ± 1.6 µK2, or A kSZ = 2.6 ± 1.8 µK2 if the default A kSZ > 0 prior is removed.
