We discuss a variational problem defined on couples of functions that are constrained totake values into the 2-dimensional unit sphere. The energy functional contains, besidesstandard Dirichlet energies, a non-local interaction term that depends on the distancebetween the gradients of the two functions. Different gradients are preferred or penalizedin this model, in dependence of the sign of the interaction term. In this paper we studythe lower semicontinuity and the coercivity of the energy and we find an explicitrepresentation formula for the relaxed energy. Moreover, we discuss the behavior of theenergy in the case when we consider multifunctions with two leaves rather than couples offunctions.