AbstractWe derive a shape derivative formula for the family of principal Dirichlet eigenvalues $$\lambda _s(\Omega )$$
λ
s
(
Ω
)
of the fractional Laplacian $$(-\Delta )^s$$
(
-
Δ
)
s
associated with bounded open sets $$\Omega \subset \mathbb {R}^N$$
Ω
⊂
R
N
of class $$C^{1,1}$$
C
1
,
1
. This extends, with a help of a new approach, a result in Dalibard and Gérard-Varet (Calc. Var. 19(4):976–1013, 2013) which was restricted to the case $$s=\frac{1}{2}$$
s
=
1
2
. As an application, we consider the maximization problem for $$\lambda _s(\Omega )$$
λ
s
(
Ω
)
among annular-shaped domains of fixed volume of the type $$B\setminus \overline{B}'$$
B
\
B
¯
′
, where B is a fixed ball and $$B'$$
B
′
is ball whose position is varied within B. We prove that $$\lambda _s(B\setminus \overline{B}')$$
λ
s
(
B
\
B
¯
′
)
is maximal when the two balls are concentric. Our approach also allows to derive similar results for the fractional torsional rigidity. More generally, we will characterize one-sided shape derivatives for best constants of a family of subcritical fractional Sobolev embeddings.