# Riesz–Thorin theorem

In mathematics, the Riesz–Thorin theorem, often referred to as the Riesz–Thorin interpolation theorem or the Riesz–Thorin convexity theorem, is a result about interpolation of operators. It is named after Marcel Riesz and his student G. Olof Thorin.
This theorem bounds the norms of linear maps acting between L$^p$ spaces. Its usefulness stems from the fact that some of these spaces have rather simpler structure than others. Usually that refers to L$^2$ which is a Hilbert space, or to L$^1$ and L$^∞$. Therefore one may prove theorems about the more complicated cases by proving them in two simple cases and then using the Riesz–Thorin theorem to pass from the simple cases to the complicated cases. The Marcinkiewicz theorem is similar but applies also to a class of non-linear maps.

## Motivation

First we need the following definition:
    Definition. Let p$_0$, p$_1$ be two numbers such that 1 ≤ p$_0$ < p$_1$ ≤ ∞. Then for 0 < θ < 1 define p$_θ$ by:  ⁠1/p$_θ$⁠ = ⁠1 − θ/p$_0$⁠ + ⁠θ/p$_1$⁠.
By splitting up the function  f  in L$^{pθ}$ as the product | f | = | f |$^{1−θ}$ | f |$^θ$ and applying Hölder's inequality to its p$_θ$ power, we obtain the following result, foundational in the study of L$^p$-spaces:
Proposition (log-convexity of L$^p$-norms) — Each  f  ∈ L$^{p0}$ ∩ L$^{p1}$ satisfies:
This result, whose name derives from the convexity of the map 1⁄p ↦ log || f ||$_p$ on [0, ∞], implies that L$^{p0}$ ∩ L$^{p1}$ ⊂ L$^{pθ}$.
On the other hand, if we take the layer-cake decomposition  f  =  f 1$_{{|f|>1}}$ +  f 1$_{{|f|≤1}}$, then we see that  f 1$_{{|f|>1}}$ ∈ L$^{p0}$ and  f 1$_{{|f|≤1}}$ ∈ L$^{p1}$, whence we obtain the following result:
Proposition — Each  f  in L$^{pθ}$ can be written as a sum:  f  = g + h, where g ∈ L$^{p0}$ and h ∈ L$^{p1}$.
In particular, the above result implies that L$^{pθ}$ is included in L$^{p0}$ + L$^{p1}$, the sumset of L$^{p0}$ and L$^{p1}$ in the space of all measurable functions. Therefore, we have the following chain of inclusions:
Corollary — L$^{p0}$ ∩ L$^{p1}$ ⊂ L$^{pθ}$ ⊂ L$^{p0}$ + L$^{p1}$.
In practice, we often encounter operators defined on the sumset L$^{p0}$ + L$^{p1}$. For example, the Riemann–Lebesgue lemma shows that the Fourier transform maps L$^1$(R$^d$) boundedly into L$^∞$(R$^d$), and Plancherel's theorem shows that the Fourier transform maps L$^2$(R$^d$) boundedly into itself, hence the Fourier transform ${\mathcal {F}}$ extends to (L$^1$ + L$^2$) (R$^d$) by setting ${\mathcal {F}}(f_{1}+f_{2})={\mathcal {F}}_{L^{1}}(f_{1})+{\mathcal {F}}_{L^{2}}(f_{2})$ for all  f$_1$  ∈ L$^1$(R$^d$) and  f$_2$  ∈ L$^2$(R$^d$). It is therefore natural to investigate the behavior of such operators on the intermediate subspaces L$^{pθ}$.
To this end, we go back to our example and note that the Fourier transform on the sumset L$^1$ + L$^2$ was obtained by taking the sum of two instantiations of the same operator, namely ${\mathcal {F}}_{L^{1}}:L^{1}(\mathbf {R} ^{d})\to L^{\infty }(\mathbf {R} ^{d}),$ ${\mathcal {F}}_{L^{2}}:L^{2}(\mathbf {R} ^{d})\to L^{2}(\mathbf {R} ^{d}).$
These really are the same operator, in the sense that they agree on the subspace (L$^1$ ∩ L$^2$) (R$^d$). Since the intersection contains simple functions, it is dense in both L$^1$(R$^d$) and L$^2$(R$^d$). Densely defined continuous operators admit unique extensions, and so we are justified in considering ${\mathcal {F}}_{L^{1}}$ and ${\mathcal {F}}_{L^{2}}$ to be the same.
Therefore, the problem of studying operators on the sumset L$^{p0}$ + L$^{p1}$ essentially reduces to the study of operators that map two natural domain spaces, L$^{p0}$ and L$^{p1}$, boundedly to two target spaces: L$^{q0}$ and L$^{q1}$, respectively. Since such operators map the sumset space L$^{p0}$ + L$^{p1}$ to L$^{q0}$ + L$^{q1}$, it is natural to expect that these operators map the intermediate space L$^{pθ}$ to the corresponding intermediate space L$^{qθ}$.

## Statement of the theorem

There are several ways to state the Riesz–Thorin interpolation theorem; to be consistent with the notations in the previous section, we shall use the sumset formulation.
Riesz–Thorin interpolation theorem — Let (Ω$_1$, Σ$_1$, μ$_1$) and (Ω$_2$, Σ$_2$, μ$_2$) be σ-finite measure spaces. Suppose 1 ≤ p$_0$ , q$_0$ , p$_1$ , q$_1$ ≤ ∞, and let T : L$^{p0}$(μ$_1$) + L$^{p1}$(μ$_1$) → L$^{q0}$(μ$_2$) + L$^{q1}$(μ$_2$) be a linear operator that boundedly maps L$^{p0}$(μ$_1$) into L$^{q0}$(μ$_2$) and L$^{p1}$(μ$_1$) into L$^{q1}$(μ$_2$). For 0 < θ < 1, let p$_θ$, q$_θ$ be defined as above. Then T boundedly maps L$^{pθ}$(μ$_1$) into L$^{qθ}$(μ$_2$) and satisfies the operator norm estimate
In other words, if T is simultaneously of type (p$_0$, q$_0$) and of type (p$_1$, q$_1$), then T is of type (p$_θ$, q$_θ$) for all 0 < θ < 1. In this manner, the interpolation theorem lends itself to a pictorial description. Indeed, the Riesz diagram of T is the collection of all points (⁠1/p⁠, ⁠1/q⁠) in the unit square [0, 1] × [0, 1] such that T is of type (p, q). The interpolation theorem states that the Riesz diagram of T is a convex set: given two points in the Riesz diagram, the line segment that connects them will also be in the diagram.
The interpolation theorem was originally stated and proved by Marcel Riesz in 1927. The 1927 paper establishes the theorem only for the lower triangle of the Riesz diagram, viz., with the restriction that p$_0$ ≤ q$_0$ and p$_1$ ≤ q$_1$. Olof Thorin extended the interpolation theorem to the entire square, removing the lower-triangle restriction. The proof of Thorin was originally published in 1938 and was subsequently expanded upon in his 1948 thesis.

## Proof

We will first prove the result for simple functions and eventually show how the argument can be extended by density to all measurable functions.

### Simple functions

By symmetry, let us assume p$_0$<p$_1${\textstyle p_{0}<p_{1}} (the case p$_0$=p$_1${\textstyle p_{0}=p_{1}} trivially follows from (1)). Let f{\textstyle f} be a simple function, that is $f=\sum _{j=1}^{m}a_{j}\mathbf {1} _{A_{j}}$ for some finite m∈N{\textstyle m\in \mathbb {N} }, a$_j$=|a$_j$|e$^{iαj}$∈C{\textstyle a_{j}=\left\vert a_{j}\right\vert \mathrm {e} ^{\mathrm {i} \alpha _{j}}\in \mathbb {C} } and A$_j$∈Σ$_1${\textstyle A_{j}\in \Sigma _{1}}, j=1, 2, …, m{\textstyle j=1, 2, \dots , m}. Similarly, let g{\textstyle g} denote a simple function Ω$_2$→C{\textstyle \Omega _{2}\to \mathbb {C} }, namely $g=\sum _{k=1}^{n}b_{k}\mathbf {1} _{B_{k}}$ for some finite n∈N{\textstyle n\in \mathbb {N} }, b$_k$=|b$_k$|e$^{iβk}$∈C{\textstyle b_{k}=\left\vert b_{k}\right\vert \mathrm {e} ^{\mathrm {i} \beta _{k}}\in \mathbb {C} } and B$_k$∈Σ$_2${\textstyle B_{k}\in \Sigma _{2}}, k=1, 2, …, n{\textstyle k=1, 2, \dots , n}.
Note that, since we are assuming Ω$_1${\textstyle \Omega _{1}} and Ω$_2${\textstyle \Omega _{2}} to be σ{\textstyle \sigma }-finite metric spaces, f∈L$^r$(μ$_1$){\textstyle f\in L^{r}(\mu _{1})} and g∈L$^r$(μ$_2$){\textstyle g\in L^{r}(\mu _{2})} for all r∈[1, ∞]{\textstyle r\in [1, \infty ]}. Then, by proper normalization, we can assume ‖f‖$_{pθ}$=1{\textstyle \lVert f\rVert _{p_{\theta }}=1} and ‖g‖$_{qθ′}$=1{\textstyle \lVert g\rVert _{q_{\theta }'}=1}, with q$_θ$$^′$=q$_θ$(q$_θ$−1)$^{−1}${\textstyle q_{\theta }'=q_{\theta }(q_{\theta }-1)^{-1}} and with p$_θ${\textstyle p_{\theta }}, q$_θ${\textstyle q_{\theta }} as defined by the theorem statement.
Next, we define the two complex functions ${\begin{aligned}u:\mathbb {C} &\to \mathbb {C} &v:\mathbb {C} &\to \mathbb {C} \\z&\mapsto u(z)={\frac {1-z}{p_{0}}}+{\frac {z}{p_{1}}}&z&\mapsto v(z)={\frac {1-z}{q_{0}}}+{\frac {z}{q_{1}}}.\end{aligned}}$ Note that, for z=θ{\textstyle z=\theta }, u(θ)=p$_θ$$^{−1}${\textstyle u(\theta )=p_{\theta }^{-1}} and v(θ)=q$_θ$$^{−1}${\textstyle v(\theta )=q_{\theta }^{-1}}. We then extend f{\textstyle f} and g{\textstyle g} to depend on a complex parameter z{\textstyle z} as follows: ${\begin{aligned}f_{z}&=\sum _{j=1}^{m}\left\vert a_{j}\right\vert ^{\frac {u(z)}{u(\theta )}}\mathrm {e} ^{\mathrm {i} \alpha _{j}}\mathbf {1} _{A_{j}}\\g_{z}&=\sum _{k=1}^{n}\left\vert b_{k}\right\vert ^{\frac {1-v(z)}{1-v(\theta )}}\mathrm {e} ^{\mathrm {i} \beta _{k}}\mathbf {1} _{B_{k}}\end{aligned}}$ so that f$_θ$=f{\textstyle f_{\theta }=f} and g$_θ$=g{\textstyle g_{\theta }=g}. Here, we are implicitly excluding the case q$_0$=q$_1$=1{\textstyle q_{0}=q_{1}=1}, which yields v≡1{\textstyle v\equiv 1}: In that case, one can simply take g$_z$=g{\textstyle g_{z}=g}, independently of z{\textstyle z}, and the following argument will only require minor adaptations.
Let us now introduce the function $\Phi (z)=\int _{\Omega _{2}}(Tf_{z})g_{z}\,\mathrm {d} \mu _{2}=\sum _{j=1}^{m}\sum _{k=1}^{n}\left\vert a_{j}\right\vert ^{\frac {u(z)}{u(\theta )}}\left\vert b_{k}\right\vert ^{\frac {1-v(z)}{1-v(\theta )}}\gamma _{j,k}$ where γ$_{j, k}$=e$^{i(αj+βk)}$∫$_{Ω2}$(T1$_{Aj}$)1$_{Bk}$dμ$_2${\textstyle \gamma _{j, k}=\mathrm {e} ^{\mathrm {i} (\alpha _{j}+\beta _{k})}\int _{\Omega _{2}}(T\mathbf {1} _{A_{j}})\mathbf {1} _{B_{k}}\, \mathrm {d} \mu _{2}} are constants independent of z{\textstyle z}. We readily see that Φ(z){\textstyle \Phi (z)} is an entire function, bounded on the strip 0≤Re⁡z≤1{\textstyle 0\leq \operatorname {\mathbb {R} e} z\leq 1}. Then, in order to prove (2), we only need to show that
for all f$_z${\textstyle f_{z}} and g$_z${\textstyle g_{z}} as constructed above. Indeed, if (3) holds true, by Hadamard three-lines theorem, $\left\vert \Phi (\theta +\mathrm {i} 0)\right\vert ={\biggl \vert }\int _{\Omega _{2}}(Tf)g\,\mathrm {d} \mu _{2}{\biggr \vert }\leq \|T\|_{L^{p_{0}}\to L^{q_{0}}}^{1-\theta }\|T\|_{L^{p_{1}}\to L^{q_{1}}}^{\theta }$ for all f{\textstyle f} and g{\textstyle g}. This means, by fixing f{\textstyle f}, that $\sup _{g}{\biggl \vert }\int _{\Omega _{2}}(Tf)g\,\mathrm {d} \mu _{2}{\biggr \vert }\leq \|T\|_{L^{p_{0}}\to L^{q_{0}}}^{1-\theta }\|T\|_{L^{p_{1}}\to L^{q_{1}}}^{\theta }$ where the supremum is taken with respect to all g{\textstyle g} simple functions with ‖g‖$_{qθ′}$=1{\textstyle \lVert g\rVert _{q_{\theta }'}=1}. The left-hand side can be rewritten by means of the following lemma.
Lemma — Let 1≤p, p$^′$≤∞{\textstyle 1\leq p, p'\leq \infty } be conjugate exponents and let f{\textstyle f} be a function in L$^p$(μ$_1$){\textstyle L^{p}(\mu _{1})}. Then $\lVert f\rVert _{p}=\sup {\biggl |}\int _{\Omega _{1}}fg\,\mathrm {d} \mu _{1}{\biggr |}$ where the supremum is taken over all simple functions g{\textstyle g} in L$^{p′}$(μ$_1$){\textstyle L^{p'}(\mu _{1})} such that ‖g‖$_{p′}$≤1{\textstyle \lVert g\rVert _{p'}\leq 1}.
In our case, the lemma above implies $\lVert Tf\rVert _{q_{\theta }}\leq \|T\|_{L^{p_{0}}\to L^{q_{0}}}^{1-\theta }\|T\|_{L^{p_{1}}\to L^{q_{1}}}^{\theta }$ for all simple function f{\textstyle f} with ‖f‖$_{pθ}$=1{\textstyle \lVert f\rVert _{p_{\theta }}=1}. Equivalently, for a generic simple function, $\lVert Tf\rVert _{q_{\theta }}\leq \|T\|_{L^{p_{0}}\to L^{q_{0}}}^{1-\theta }\|T\|_{L^{p_{1}}\to L^{q_{1}}}^{\theta }\lVert f\rVert _{p_{\theta }}.$

#### Proof of (3)

Let us now prove that our claim (3) is indeed certain. The sequence (A$_j$)$_{j=1}$$^m${\textstyle (A_{j})_{j=1}^{m}} consists of disjoint subsets in Σ$_1${\textstyle \Sigma _{1}} and, thus, each ξ∈Ω$_1${\textstyle \xi \in \Omega _{1}} belongs to (at most) one of them, say A$_{ȷ^}${\textstyle A_{\hat {\jmath }}}. Then, for z=iy{\textstyle z=\mathrm {i} y}, ${\begin{aligned}\left\vert f_{\mathrm {i} y}(\xi )\right\vert &=\left\vert \left\vert a_{\hat {\jmath }}\right\vert ^{\frac {u(\mathrm {i} y)}{u(\theta )}}\right\vert \\&=\left\vert \exp {\biggl\log \left\vert a_{\hat {\jmath }}\right\vert {\frac {p_{\theta }}{p_{0}}}{\biggr )}\exp {\biggl-\mathrm {i} y\log \left\vert a_{\hat {\jmath }}\right\vert p_{\theta }{\biggl{\frac {1}{p_{0}}}-{\frac {1}{p_{1}}}{\biggr )}{\biggr )}\right\vert \\&=\left\vert a_{\hat {\jmath }}\right\vert ^{\frac {p_{\theta }}{p_{0}}}\\&=\left\vert f(\xi )\right\vert ^{\frac {p_{\theta }}{p_{0}}}\end{aligned}}$ which implies that ‖f$_{iy}$‖$_{p0}$≤‖f‖$_{pθ}$$^{pθ/p0}${\textstyle \lVert f_{\mathrm {i} y}\rVert _{p_{0}}\leq \lVert f\rVert _{p_{\theta }}^{\frac {p_{\theta }}{p_{0}}}}. With a parallel argument, each ζ∈Ω$_2${\textstyle \zeta \in \Omega _{2}} belongs to (at most) one of the sets supporting g{\textstyle g}, say B$_{k^}${\textstyle B_{\hat {k}}}, and $\left\vert g_{\mathrm {i} y}(\zeta )\right\vert =\left\vert b_{\hat {k}}\right\vert ^{\frac {1-1/q_{0}}{1-1/q_{\theta }}}=\left\vert g(\zeta )\right\vert ^{\frac {1-1/q_{0}}{1-1/q_{\theta }}}=\left\vert g(\zeta )\right\vert ^{\frac {q_{\theta }'}{q_{0}'}}\implies \lVert g_{\mathrm {i} y}\rVert _{q_{0}'}\leq \lVert g\rVert _{q_{\theta }'}^{\frac {q_{\theta }'}{q_{0}'}}.$
We can now bound Φ(iy){\textstyle \Phi (\mathrm {i} y)}: By applying Hölder’s inequality with conjugate exponents q$_0${\textstyle q_{0}} and q$_0$$^′${\textstyle q_{0}'}, we have ${\begin{aligned}\left\vert \Phi (\mathrm {i} y)\right\vert &\leq \lVert Tf_{\mathrm {i} y}\rVert _{q_{0}}\lVert g_{\mathrm {i} y}\rVert _{q_{0}'}\\&\leq \|T\|_{L^{p_{0}}\to L^{q_{0}}}\lVert f_{\mathrm {i} y}\rVert _{p_{0}}\lVert g_{\mathrm {i} y}\rVert _{q_{0}'}\\&=\|T\|_{L^{p_{0}}\to L^{q_{0}}}\lVert f\rVert _{p_{\theta }}^{\frac {p_{\theta }}{p_{0}}}\lVert g\rVert _{q_{\theta }'}^{\frac {q_{\theta }'}{q_{0}'}}\\&=\|T\|_{L^{p_{0}}\to L^{q_{0}}}.\end{aligned}}$
We can repeat the same process for z=1+iy{\textstyle z=1+\mathrm {i} y} to obtain |f$_{1+iy}$(ξ)|=|f(ξ)|$^{pθ/p1}${\textstyle \left\vert f_{1+\mathrm {i} y}(\xi )\right\vert =\left\vert f(\xi )\right\vert ^{p_{\theta }/p_{1}}}, |g$_{1+iy}$(ζ)|=|g(ζ)|$^{qθ′/q1′}${\textstyle \left\vert g_{1+\mathrm {i} y}(\zeta )\right\vert =\left\vert g(\zeta )\right\vert ^{q_{\theta }'/q_{1}'}} and, finally, $\left\vert \Phi (1+\mathrm {i} y)\right\vert \leq \|T\|_{L^{p_{1}}\to L^{q_{1}}}\lVert f_{1+\mathrm {i} y}\rVert _{p_{1}}\lVert g_{1+\mathrm {i} y}\rVert _{q_{1}'}=\|T\|_{L^{p_{1}}\to L^{q_{1}}}.$

### Extension to all measurable functions in L$^{pθ}$

So far, we have proven that
when f{\textstyle f} is a simple function. As already mentioned, the inequality holds true for all f∈L$^{pθ}$(Ω$_1$){\textstyle f\in L^{p_{\theta }}(\Omega _{1})} by the density of simple functions in L$^{pθ}$(Ω$_1$){\textstyle L^{p_{\theta }}(\Omega _{1})}.
Formally, let f∈L$^{pθ}$(Ω$_1$){\textstyle f\in L^{p_{\theta }}(\Omega _{1})} and let (f$_n$)$_n${\textstyle (f_{n})_{n}} be a sequence of simple functions such that |f$_n$|≤|f|{\textstyle \left\vert f_{n}\right\vert \leq \left\vert f\right\vert }, for all n{\textstyle n}, and f$_n$→f{\textstyle f_{n}\to f} pointwise. Let E={x∈Ω$_1$:|f(x)|>1}{\textstyle E=\{x\in \Omega _{1}:\left\vert f(x)\right\vert >1\}} and define g=f1$_E${\textstyle g=f\mathbf {1} _{E}}, g$_n$=f$_n$1$_E${\textstyle g_{n}=f_{n}\mathbf {1} _{E}}, h=f−g=f1$_{Ec}${\textstyle h=f-g=f\mathbf {1} _{E^{\mathrm {c} }}} and h$_n$=f$_n$−g$_n${\textstyle h_{n}=f_{n}-g_{n}}. Note that, since we are assuming p$_0$≤p$_θ$≤p$_1${\textstyle p_{0}\leq p_{\theta }\leq p_{1}}, ${\begin{aligned}\lVert f\rVert _{p_{\theta }}^{p_{\theta }}&=\int _{\Omega _{1}}\left\vert f\right\vert ^{p_{\theta }}\,\mathrm {d} \mu _{1}\geq \int _{\Omega _{1}}\left\vert f\right\vert ^{p_{\theta }}\mathbf {1} _{E}\,\mathrm {d} \mu _{1}\geq \int _{\Omega _{1}}\left\vert f\mathbf {1} _{E}\right\vert ^{p_{0}}\,\mathrm {d} \mu _{1}=\int _{\Omega _{1}}\left\vert g\right\vert ^{p_{0}}\,\mathrm {d} \mu _{1}=\lVert g\rVert _{p_{0}}^{p_{0}}\\\lVert f\rVert _{p_{\theta }}^{p_{\theta }}&=\int _{\Omega _{1}}\left\vert f\right\vert ^{p_{\theta }}\,\mathrm {d} \mu _{1}\geq \int _{\Omega _{1}}\left\vert f\right\vert ^{p_{\theta }}\mathbf {1} _{E^{\mathrm {c} }}\,\mathrm {d} \mu _{1}\geq \int _{\Omega _{1}}\left\vert f\mathbf {1} _{E^{\mathrm {c} }}\right\vert ^{p_{1}}\,\mathrm {d} \mu _{1}=\int _{\Omega _{1}}\left\vert h\right\vert ^{p_{1}}\,\mathrm {d} \mu _{1}=\lVert h\rVert _{p_{1}}^{p_{1}}\end{aligned}}$ and, equivalently, g∈L$^{p0}$(Ω$_1$){\textstyle g\in L^{p_{0}}(\Omega _{1})} and h∈L$^{p1}$(Ω$_1$){\textstyle h\in L^{p_{1}}(\Omega _{1})}.
Let us see what happens in the limit for n→∞{\textstyle n\to \infty }. Since |f$_n$|≤|f|{\textstyle \left\vert f_{n}\right\vert \leq \left\vert f\right\vert }, |g$_n$|≤|g|{\textstyle \left\vert g_{n}\right\vert \leq \left\vert g\right\vert } and |h$_n$|≤|h|{\textstyle \left\vert h_{n}\right\vert \leq \left\vert h\right\vert }, by the dominated convergence theorem one readily has ${\begin{aligned}\lVert f_{n}\rVert _{p_{\theta }}&\to \lVert f\rVert _{p_{\theta }}&\lVert g_{n}\rVert _{p_{0}}&\to \lVert g\rVert _{p_{0}}&\lVert h_{n}\rVert _{p_{1}}&\to \lVert h\rVert _{p_{1}}.\end{aligned}}$ Similarly, |f−f$_n$|≤2|f|{\textstyle \left\vert f-f_{n}\right\vert \leq 2\left\vert f\right\vert }, |g−g$_n$|≤2|g|{\textstyle \left\vert g-g_{n}\right\vert \leq 2\left\vert g\right\vert } and |h−h$_n$|≤2|h|{\textstyle \left\vert h-h_{n}\right\vert \leq 2\left\vert h\right\vert } imply ${\begin{aligned}\lVert f-f_{n}\rVert _{p_{\theta }}&\to 0&\lVert g-g_{n}\rVert _{p_{0}}&\to 0&\lVert h-h_{n}\rVert _{p_{1}}&\to 0\end{aligned}}$ and, by the linearity of T{\textstyle T} as an operator of types (p$_0$, q$_0$){\textstyle (p_{0}, q_{0})} and (p$_1$, q$_1$){\textstyle (p_{1}, q_{1})} (we have not proven yet that it is of type (p$_θ$, q$_θ$){\textstyle (p_{\theta }, q_{\theta })} for a generic f{\textstyle f}) ${\begin{aligned}\lVert Tg-Tg_{n}\rVert _{p_{0}}&\leq \|T\|_{L^{p_{0}}\to L^{q_{0}}}\lVert g-g_{n}\rVert _{p_{0}}\to 0&\lVert Th-Th_{n}\rVert _{p_{1}}&\leq \|T\|_{L^{p_{1}}\to L^{q_{1}}}\lVert h-h_{n}\rVert _{p_{1}}\to 0.\end{aligned}}$
It is now easy to prove that Tg$_n$→Tg{\textstyle Tg_{n}\to Tg} and Th$_n$→Th{\textstyle Th_{n}\to Th} in measure: For any ϵ>0{\textstyle \epsilon >0}, Chebyshev’s inequality yields $\mu _{2}(y\in \Omega _{2}:\left\vert Tg-Tg_{n}\right\vert >\epsilon )\leq {\frac {\lVert Tg-Tg_{n}\rVert _{q_{0}}^{q_{0}}}{\epsilon ^{q_{0}}}}$ and similarly for Th−Th$_n${\textstyle Th-Th_{n}}. Then, Tg$_n$→Tg{\textstyle Tg_{n}\to Tg} and Th$_n$→Th{\textstyle Th_{n}\to Th} a.e. for some subsequence and, in turn, Tf$_n$→Tf{\textstyle Tf_{n}\to Tf} a.e. Then, by Fatou’s lemma and recalling that (4) holds true for simple functions, $\lVert Tf\rVert _{q_{\theta }}\leq \liminf _{n\to \infty }\lVert Tf_{n}\rVert _{q_{\theta }}\leq \|T\|_{L^{p_{\theta }}\to L^{q_{\theta }}}\liminf _{n\to \infty }\lVert f_{n}\rVert _{p_{\theta }}=\|T\|_{L^{p_{\theta }}\to L^{q_{\theta }}}\lVert f\rVert _{p_{\theta }}.$

## Interpolation of analytic families of operators

The proof outline presented in the above section readily generalizes to the case in which the operator T is allowed to vary analytically. In fact, an analogous proof can be carried out to establish a bound on the entire function $\varphi (z)=\int (T_{z}f_{z})g_{z}\,d\mu _{2},$ from which we obtain the following theorem of Elias Stein, published in his 1956 thesis:
Stein interpolation theorem — Let (Ω$_1$, Σ$_1$, μ$_1$) and (Ω$_2$, Σ$_2$, μ$_2$) be σ-finite measure spaces. Suppose 1 ≤ p$_0$ , p$_1$ ≤ ∞, 1 ≤ q$_0$ , q$_1$ ≤ ∞, and define:
    S = {z ∈ C : 0 < Re(z) < 1},
    S = {z ∈ C : 0 ≤ Re(z) ≤ 1}.
We take a collection of linear operators {T$_z$ : z ∈ S} on the space of simple functions in L$^1$(μ$_1$) into the space of all μ$_2$-measurable functions on Ω$_2$. We assume the following further properties on this collection of linear operators:
* The mapping $z\mapsto \int (T_{z}f)g\,d\mu _{2}$ is continuous on S and holomorphic on S for all simple functions  f  and g.
* For some constant k < π, the operators satisfy the uniform bound: $\sup _{z\in S}e^{-k|{\text{Im}}(z)|}\log \left|\int (T_{z}f)g\,d\mu _{2}\right|<\infty$
* T$_z$ maps L$^{p0}$(μ$_1$) boundedly to L$^{q0}$(μ$_2$) whenever Re(z) = 0.
* T$_z$ maps L$^{p1}$(μ$_1$) boundedly to L$^{q1}$(μ$_2$) whenever Re(z) = 1.
* The operator norms satisfy the uniform bound $\sup _{{\text{Re}}(z)=0,1}e^{-k|{\text{Im}}(z)|}\log \left\|T_{z}\right\|<\infty$ for some constant k < π.
Then, for each 0 < θ < 1, the operator T$_θ$ maps L$^{pθ}$(μ$_1$) boundedly into L$^{qθ}$(μ$_2$).
The theory of real Hardy spaces and the space of bounded mean oscillations permits us to wield the Stein interpolation theorem argument in dealing with operators on the Hardy space H$^1$(R$^d$) and the space BMO of bounded mean oscillations; this is a result of Charles Fefferman and Elias Stein.

## Applications
### Hausdorff–Young inequality

It has been shown in the first section that the Fourier transform ${\mathcal {F}}$ maps L$^1$(R$^d$) boundedly into L$^∞$(R$^d$) and L$^2$(R$^d$) into itself. A similar argument shows that the Fourier series operator, which transforms periodic functions  f  : T → C into functions ${\hat {f}}:\mathbf {Z} \to \mathbf {C}$ whose values are the Fourier coefficients ${\hat {f}}(n)={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(x)e^{-inx}\,dx,$ maps L$^1$(T) boundedly into ℓ$^∞$(Z) and L$^2$(T) into ℓ$^2$(Z). The Riesz–Thorin interpolation theorem now implies the following: ${\begin{aligned}\left\|{\mathcal {F}}f\right\|_{L^{q}(\mathbf {R} ^{d})}&\leq \|f\|_{L^{p}(\mathbf {R} ^{d})}\\\left\|{\hat {f}}\right\|_{\ell ^{q}(\mathbf {Z} )}&\leq \|f\|_{L^{p}(\mathbf {T} )}\end{aligned}}$ where 1 ≤ p ≤ 2 and ⁠1/p⁠ + ⁠1/q⁠ = 1. This is the Hausdorff–Young inequality.
The Hausdorff–Young inequality can also be established for the Fourier transform on locally compact Abelian groups. The norm estimate of 1 is not optimal. See the main article for references.

### Convolution operators

Let  f  be a fixed integrable function and let T be the operator of convolution with  f , i.e., for each function g we have Tg =  f  ∗ g.
It follows from Fubini's theorem that T is bounded from L$^1$ to L$^1$ and it is trivial that it is bounded from L$^∞$ to L$^∞$ (both bounds are by || f ||$_1$). Therefore the Riesz–Thorin theorem gives $\|f*g\|_{p}\leq \|f\|_{1}\|g\|_{p}.$
We take this inequality and switch the role of the operator and the operand, or in other words, we think of S as the operator of convolution with g, and get that S is bounded from L$^1$ to L$^p$. Further, since g is in L$^p$ we get, in view of Hölder's inequality, that S is bounded from L$^q$ to L$^∞$, where again ⁠1/p⁠ + ⁠1/q⁠ = 1. So interpolating we get $\|f*g\|_{s}\leq \|f\|_{r}\|g\|_{p}$ where the connection between p, r and s is ${\frac {1}{r}}+{\frac {1}{p}}=1+{\frac {1}{s}}.$

### The Hilbert transform

The Hilbert transform of  f  : R → C is given by ${\mathcal {H}}f(x)={\frac {1}{\pi }}\,\mathrm {p.v.} \int _{-\infty }^{\infty }{\frac {f(x-t)}{t}}\,dt=\left({\frac {1}{\pi }}\,\mathrm {p.v.} {\frac {1}{t}}\ast f\right)(x),$ where p.v. indicates the Cauchy principal value of the integral. The Hilbert transform is a Fourier multiplier operator with a particularly simple multiplier: ${\widehat {{\mathcal {H}}f}}(\xi )=-i\,\operatorname {sgn}(\xi ){\hat {f}}(\xi ).$
It follows from the Plancherel theorem that the Hilbert transform maps L$^2$(R) boundedly into itself.
Nevertheless, the Hilbert transform is not bounded on L$^1$(R) or L$^∞$(R), and so we cannot use the Riesz–Thorin interpolation theorem directly. To see why we do not have these endpoint bounds, it suffices to compute the Hilbert transform of the simple functions 1$_{(−1,1)}$(x) and 1$_{(0,1)}$(x) − 1$_{(0,1)}$(−x). We can show, however, that $({\mathcal {H}}f)^{2}=f^{2}+2{\mathcal {H}}(f{\mathcal {H}}f)$ for all Schwartz functions  f  : R → C, and this identity can be used in conjunction with the Cauchy–Schwarz inequality to show that the Hilbert transform maps L$^{2n}$(R$^d$) boundedly into itself for all n ≥ 2. Interpolation now establishes the bound $\|{\mathcal {H}}f\|_{p}\leq A_{p}\|f\|_{p}$ for all 2 ≤ p < ∞, and the self-adjointness of the Hilbert transform can be used to carry over these bounds to the 1 < p ≤ 2 case.

## Comparison with the real interpolation method

While the Riesz–Thorin interpolation theorem and its variants are powerful tools that yield a clean estimate on the interpolated operator norms, they suffer from numerous defects: some minor, some more severe. Note first that the complex-analytic nature of the proof of the Riesz–Thorin interpolation theorem forces the scalar field to be C. For extended-real-valued functions, this restriction can be bypassed by redefining the function to be finite everywhere—possible, as every integrable function must be finite almost everywhere. A more serious disadvantage is that, in practice, many operators, such as the Hardy–Littlewood maximal operator and the Calderón–Zygmund operators, do not have good endpoint estimates. In the case of the Hilbert transform in the previous section, we were able to bypass this problem by explicitly computing the norm estimates at several midway points. This is cumbersome and is often not possible in more general scenarios. Since many such operators satisfy the weak-type estimates $\mu \left(\{x:Tf(x)>\alpha \}\right)\leq \left({\frac {C_{p,q}\|f\|_{p}}{\alpha }}\right)^{q},$ real interpolation theorems such as the Marcinkiewicz interpolation theorem are better-suited for them. Furthermore, a good number of important operators, such as the Hardy-Littlewood maximal operator, are only sublinear. This is not a hindrance to applying real interpolation methods, but complex interpolation methods are ill-equipped to handle non-linear operators. On the other hand, real interpolation methods, compared to complex interpolation methods, tend to produce worse estimates on the intermediate operator norms and do not behave as well off the diagonal in the Riesz diagram. The off-diagonal versions of the Marcinkiewicz interpolation theorem require the formalism of Lorentz spaces and do not necessarily produce norm estimates on the L$^p$-spaces.

## Mityagin's theorem

B. Mityagin extended the Riesz–Thorin theorem; this extension is formulated here in the special case of spaces of sequences with unconditional bases (cf. below).
Assume: $\|A\|_{\ell _{1}\to \ell _{1}},\|A\|_{\ell _{\infty }\to \ell _{\infty }}\leq M.$
Then $\|A\|_{X\to X}\leq M$
for any unconditional Banach space of sequences X, that is, for any $(x_{i})\in X$ and any $(\varepsilon _{i})\in \{-1,1\}^{\infty }$, $\|(\varepsilon _{i}x_{i})\|_{X}=\|(x_{i})\|_{X}$.
The proof is based on the Krein–Milman theorem.