Gradient Estimate of subelliptic harmonic maps with potential

In this paper, we investigate subelliptic harmonic maps with potential from noncompact complete sub-Riemannian manifolds corresponding to totally geodesic Riemannian foliations. Under some suitable conditions, we give the gradient estimates of these maps and establish a Liouville type result.