Linking loops in ABJM and refined theory

We consider the link average of the half-BPS Wilson loop operators in N = 6 superconformal Chern-Simons-matter theory, which is called ABJM theory. We show that this loop average is reduced to a (super)matrix integral by the localization method, in a similar way to the bosonic U(N) Chern-Simons theory. Using this matrix integral, we compute the two- and three-link averages with an operator formalism inspired by a three-dimensional topological field theory. We obtain a factorization of the link average, and the Verlinde formula in a sector of supergroup representations. We also propose a refined version of ABJM theory, and compute some refined link averages.