A matrix model for the latitude Wilson loop in ABJM theory

In ABJ(M) theory, we propose a matrix model for the exact evaluation of BPS Wilson loops on a latitude circular contour, so providing a new weak-strong interpolation tool. Intriguingly, the matrix model turns out to be a particular case of that computing torus knot invariants in $U(N_1|N_2)$ Chern-Simons theory. At weak coupling we check our proposal against a three-loop computation, performed for generic framing, winding number and representation. The matrix model is amenable of a Fermi gas formulation, which we use to systematically compute the strong coupling and genus expansions. For the fermionic Wilson loop the leading planar behavior agrees with a previous string theory prediction. For the bosonic operator our result provides a clue for finding the corresponding string dual configuration. Our matrix model is consistent with recent proposals for computing Bremsstrahlung functions exactly in terms of latitude Wilson loops. As a by-product, we extend the conjecture for the exact $B^{\theta}_{1/6}$ Bremsstrahlung function to generic representations and test it with a four-loop perturbative computation. Finally, we propose an exact prediction for $B_{1/2}$ at unequal gauge group ranks.