Sharp L² estimate of Schr\"odinger maximal function in higher dimensions

We show that, for $n\geq 3$, $\lim_{t \to 0} e^{it\Delta}f(x) = f(x)$ holds almost everywhere for all $f \in H^s (\mathbb{R}^n)$ provided that $s>\frac{n}{2(n+1)}$. Due to a counterexample by Bourgain, up to the endpoint, this result is sharp and fully resolves a problem raised by Carleson. Our main theorem is a fractal $L^2$ restriction estimate, which also gives improved results on the size of divergence set of Schr\"odinger solutions, the Falconer distance set problem and the spherical average Fourier decay rates of fractal measures. The key ingredients of the proof include multilinear Kakeya estimates, decoupling and induction on scales.