A real random variable is a sub-Gaussian random variable with exponent when its moment-generating function satisfies
Differentiating at from both sides gives . Exponential integrability near zero justifies this differentiation. Conversely, if and , convexity gives
The last inequality follows by comparing power series using . The case is immediate. Thus a bounded random variable is sub-Gaussian exactly when it is centered; the bound is one possible exponent.
For , the Markov inequality applied to gives . Choose and repeat with . The union bound yields
If , the same exponential estimate with arbitrarily large shows almost surely, so all subsequent inequalities are trivial.
Integrating this tail probability estimate gives, for every real ,
We need . For integer this follows inductively from and . To cover real as well, set . Differentiating the Gamma function integral gives for a Gamma distribution with shape and rate . The Jensen inequality bounds this by . Hence and . Taking the th root now proves
Here is the Lp norm.
For the final assertion we prove the required Littlewood interpolation inequality, rather than assuming it. Apply Holder inequality with conjugate exponents and :
Raising to the power gives . The preceding Lp norm estimate with says . Therefore