Partial regularity of optimal transport with Coulomb cost

We prove that for two-marginal optimal transport with Coulomb cost, the optimal map is a $C^{1,\alpha }$ diffeomorphism outside a closed set of Lebesgue measure zero provided the marginals are $\alpha$-Hölder continuous and bounded away from zero and infinity. Excluding a set of measure zero is necessary as optimal maps for the Coulomb cost have long been known to exhibit jump singularities across codimension $1$ surfaces (even for smooth marginals on convex domains).