“…Aside from providing a better understanding of the geometry of the state space, the coupling method for Brownian motions is a great tool for many analysis results involving the harmonic functions and the heat semi-group such as Harnack, Poincaré, Sobolev or Wasserstein inequalities (see [23,30,15,14] for some examples). This method has been studied these last decades in the case of Riemannian manifolds (see for example [21,27]). The case of subRiemannian manifolds is a current topic of interest and have been investigated on the Heisenberg group in [9,8,21,20,19,4,11,5,25], on SU(2) [12,13,25] and on SL(2, R) [13,25].…”