  For the singular integral definition of the fractional Laplacian, we consider
an adaptive finite element method steered by two-level error indicators. For
this algorithm, we show linear convergence in two and three space dimensions as
well as convergence of the algorithm with optimal algebraic rates in 2D, when
newest vertex bisection is employed for mesh refinement.
