Exact results for Wilson loops in superconformal Chern-Simons theories with matter

We use localization techniques to compute the expectation values of supersymmetric Wilson loops in Chern-Simons theories with matter. We find the path-integral reduces to a non-Gaussian matrix model. The Wilson loops we consider preserve a single complex supersymmetry, and exist in any N = 2 theory, though the localization requires superconformal symmetry. We present explicit results for the cases of pure Chern-Simons theory with gauge group U(N), showing agreement with the known results, and ABJM, showing agreement with perturbative calculations. Our method applies to other theories, such as Gaiotto-Witten theories, BLG, and their variants.