
1 Introduction
In this paper we study a parabolic equation associated to the (normalized) infinity fractional Laplacian
operator. We recall that the local version of the game had been introduced by Peres et al. in 2009
([31]) where it is shown that the standard infinity Laplace equation is solved by the value function for
a two-players random turn “tug–of–war” game. The game is as follows: a token is initially placed at a
position x0∈Ωand every turn a fair coin is tossed to choose which of the players plays. This player
moves the token to any point in the ball of radius ε > 0around the current position. If, eventually,
iterating this process, the token reaches a point xe∈∂Ω, the players are awarded (or penalized) f(xe)
(payoff function). For a PDE overview of the infinity Laplacian operator and its role as an absolute
minimizer for the L∞norm of the gradient, see [24,25].
In 2012 Bjorland, Caffarelli and Figalli ([5]) introduced equations involving the so-called infinity
fractional Laplacian as a model for a nonlocal version of the “tug-of-war” game. Following their
explantation, instead of flipping a coin at every step, every player chooses a direction and it is an
s-stable Lévy process that chooses both the active player and the distance to travel. The infinity
fractional Laplacian, with symbol ∆s
∞, is a nonlinear integro-differential operator, the original definition
is given in Lemma 2.1 below. However, for the purpose of this paper, we also consider the alternative
equivalent definition introduced in [5] (see also [14]) given by
∆s
∞φ(x) := Cssup
|y|=1
inf
|˜y|=1 ˆ∞
0
(φ(x+ηy) + φ(x−η˜y)−2φ(x)) dη
η1+2swhere s∈(1/2,1). (1.1)
The constant is usually taken as Cs= (4ssΓ(1
2+s))/(π1
2Γ(1 −s)) but the value is irrelevant for our
results. In their paper [5] the authors study two stationary problems involving the infinity fractional
Laplacian posed in bounded space domains, namely, a Dirichlet problem and a double-obstacle problem.
Here, we consider the evolution problem
(∂tu(x, t)=∆s
∞u(x, t), x ∈Rn, t > 0,
u(x, 0) = u0(x), x ∈Rn,
(1.2)
(1.3)
with s∈(1/2,1) and n≥2. When n= 1 the operator −∆s
∞is just the usual linear fractional Laplacian
operator (−∆)sof order s, and equation (1.2) is just the well-known fractional heat equation [6,21].
See also a detailed study of that equation using PDE techniques in [18,3,7,36]. Note that for n≥2
the operator is nonlinear so a new theory is needed. A non-normalized version is introduced in [10]
along with a well-posedness theory for the corresponding equations of the type (1.2)–(1.3). However,
the two problems are not equivalent nor closely related.
Here we develop an existence theory of suitable viscosity solutions for the parabolic problem (1.2)–
(1.3), based on approximation with monotone schemes. We show that the obtained class of solutions
enjoys a number of good properties. As in the elliptic case [5], we lack a uniqueness result in the
context of viscosity solutions. However, we are able to prove an important comparison theorem relat-
ing two types of solutions, classical and viscosity ones, see Theorem 2.6. As a counterpart, we also
obtain uniqueness and comparison of classical solutions. Moreover, we show that for smooth, radially
symmetric functions and nonincreasing along the radius in Rnwith n≥2, the operator −∆s
∞reduces
to the classical fractional Laplacian (−∆)sin dimension n= 1 (Theorem 6.1). A similar example
regarding nondecreasing one-dimensional profiles can be found in Lemma 6.3. In this way we may con-
struct a large class of classical solutions that make the comparison theorem relevant (Theorems 2.10
and 2.12). Note that no similar reduction applies in general, even in the radial case (see Subsection
6.2 for a counterexample).
Using the developed tools, we study the asymptotic behavior of the constructed solutions, and obtain
a global Harnack type principle, see Theorem 2.13.
2