“…For this, we will have to handle not only the usual nonlocal character of such fractional operators, but also the difficulties given by the corresponding nonlinear behavior. As a consequence, we can make use neither of the powerful framework provided by the Caffarelli-Silvestre s-harmonic extension ( [4]) nor of various tools as, e. g., the sharp 3-commutators estimates introduced in [5] to deduce the regularity of weak fractional harmonic maps, the strong barriers and density estimates in [26,28,29], the commutator and energy estimates in [25,27], and so on. Indeed, the aforementioned tools seem not to be trivially adaptable to a nonlinear framework; also, increasing difficulties are due to the non-Hilbertian structure of the involved fractional Sobolev spaces W s,p when p is different than 2.…”