Assume that $\mathbf{J} = \{ J \}_{ij}$ are centered independent standard Gaussian with variance $1/n$ random variables for $i, j = 1, \dots, n$. Why is the function $\|\mathbf{J}\|_{\infty}$ then $1/n$-Lipschitz w.r.t to the Euclidean norm, that is, $$ \left| \sup_{\|u\| = 1} \langle u, \mathbf{J}u \rangle - \sup_{\|v\| = 1} \langle v, \mathbf{J}v \rangle \right| \leq \frac{1}{n} \|u - v\|_2 . $$ Thus, it has a sub-Gaussian tail which is $$ P(\|\mathbf{J}\|_{\infty} - E[\|\mathbf{J}\|_{\infty}] \geq x) \leq e^{-n x^2 / 2}. $$
The original statement is as follows.
Moreover, it is easy to see that $$ \|\mathbf{J}\|_\infty = \sup_{\|u\| = 1} \langle u, \mathbf{J} u \rangle $$ has sub-Gaussian tail. Indeed, is is a Lipschitz function of the Gaussian entries of $\mathbf{J}$ (with respect to the Euclidean norm) with constant bounded by $N^{-1/2}$ so that by Herbst argument (see [24, 25]), we have the concentration inequality $$ P(\|\mathbf{J}\|_\infty \geq E[\|\mathbf{J}\|_\infty] + x) \leq e^{-\frac{1}{2} N x^2} . $$