Cauchy principal value

In mathematics, the Cauchy principal value, named after Augustin-Louis Cauchy, is a method for assigning values to certain improper integrals which would otherwise be undefined. In this method, a singularity on an integral interval is avoided by limiting the integral interval to the non singular domain.

Depending on the type of singularity in the integrand f, the Cauchy principal value is defined according to the following rules:

$${\displaystyle \lim _{\;\varepsilon \to 0^{+}\;}\,\,\left[\,\int _{a}^{b-\varepsilon }f(x)\,\mathrm {d} x~+~\int _{b+\varepsilon }^{c}f(x)\,\mathrm {d} x\,\right]}$$ with ${\displaystyle a<b<c}$ and where b is the difficult point, at which the behavior of the function f is such that $${\displaystyle \int _{a}^{b}f(x)\,\mathrm {d} x=\pm \infty \quad }$$ for any ${\displaystyle a<b}$ and $${\displaystyle \int _{b}^{c}f(x)\,\mathrm {d} x=\mp \infty \quad }$$ for any ${\displaystyle c>b}$. (See plus or minus for the precise use of notations ± and ∓.)

For a singularity at infinity (${\displaystyle \infty }$)

$${\displaystyle \lim _{a\to \infty }\,\int _{-a}^{a}f(x)\,\mathrm {d} x}$$ where $${\displaystyle \int _{-\infty }^{0}f(x)\,\mathrm {d} x=\pm \infty }$$ and $${\displaystyle \int _{0}^{\infty }f(x)\,\mathrm {d} x=\mp \infty .}$$

In some cases it is necessary to deal simultaneously with singularities both at a finite number b and at infinity. This is usually done by a limit of the form $${\displaystyle \lim _{\;\eta \to 0^{+}}\,\lim _{\;\varepsilon \to 0^{+}}\,\left[\,\int _{b-{\frac {1}{\eta }}}^{b-\varepsilon }f(x)\,\mathrm {d} x\,~+~\int _{b+\varepsilon }^{b+{\frac {1}{\eta }}}f(x)\,\mathrm {d} x\,\right].}$$ In those cases where the integral may be split into two independent, finite limits, $${\displaystyle \lim _{\;\varepsilon \to 0^{+}\;}\,\left|\,\int _{a}^{b-\varepsilon }f(x)\,\mathrm {d} x\,\right|\;<\;\infty }$$ and $${\displaystyle \lim _{\;\eta \to 0^{+}}\;\left|\,\int _{b+\eta }^{c}f(x)\,\mathrm {d} x\,\right|\;<\;\infty ,}$$ then the function is integrable in the ordinary sense. The result of the procedure for principal value is the same as the ordinary integral; since it no longer matches the definition, it is technically not a "principal value". The Cauchy principal value can also be defined in terms of contour integrals of a complex-valued function ${\displaystyle f(z):z=x+i\,y\;,}$ with ${\displaystyle x,y\in \mathbb {R} \;,}$ with a pole on a contour C. Define ${\displaystyle C(\varepsilon )}$ to be that same contour, where the portion inside the disk of radius ε around the pole has been removed. Provided the function ${\displaystyle f(z)}$ is integrable over ${\displaystyle C(\varepsilon )}$ no matter how small ε becomes, then the Cauchy principal value is the limit:^[1] $${\displaystyle \operatorname {p.\!v.} \int _{C}f(z)\,\mathrm {d} z=\lim _{\varepsilon \to 0^{+}}\int _{C(\varepsilon )}f(z)\,\mathrm {d} z.}$$ In the case of Lebesgue-integrable functions, that is, functions which are integrable in absolute value, these definitions coincide with the standard definition of the integral. If the function ${\displaystyle f(z)}$ is meromorphic, the Sokhotski–Plemelj theorem relates the principal value of the integral over C with the mean-value of the integrals with the contour displaced slightly above and below, so that the residue theorem can be applied to those integrals. Principal value integrals play a central role in the discussion of Hilbert transforms.^[2]

Let ${\displaystyle {C_{c}^{\infty }}(\mathbb {R} )}$ be the set of bump functions, i.e., the space of smooth functions with compact support on the real line ${\displaystyle \mathbb {R} }$. Then the map $${\displaystyle \operatorname {p.\!v.} \left({\frac {1}{x}}\right)\,:\,{C_{c}^{\infty }}(\mathbb {R} )\to \mathbb {C} }$$ defined via the Cauchy principal value as $${\displaystyle \left[\operatorname {p.\!v.} \left({\frac {1}{x}}\right)\right](u)=\lim _{\varepsilon \to 0^{+}}\int _{\mathbb {R} \setminus [-\varepsilon ,\varepsilon ]}{\frac {u(x)}{x}}\,\mathrm {d} x=\lim _{\varepsilon \to 0^{+}}\int _{\varepsilon }^{+\infty }{\frac {u(x)-u(-x)}{x}}\,\mathrm {d} x\quad {\text{for }}u\in {C_{c}^{\infty }}(\mathbb {R} )}$$ is a distribution. The map itself may sometimes be called the principal value (hence the notation p.v.). This distribution appears, for example, in the Fourier transform of the sign function and the Heaviside step function.

To prove the existence of the limit $${\displaystyle \lim _{\varepsilon \to 0^{+}}\int _{\varepsilon }^{+\infty }{\frac {u(x)-u(-x)}{x}}\,\mathrm {d} x}$$ for a Schwartz function ${\displaystyle u(x)}$, first observe that ${\displaystyle {\frac {u(x)-u(-x)}{x}}}$ is continuous on ${\displaystyle [0,\infty ),}$ as $${\displaystyle \lim _{\,x\searrow 0\,}\;{\Bigl [}u(x)-u(-x){\Bigr ]}~=~0~}$$ and hence $${\displaystyle \lim _{x\searrow 0}\,{\frac {u(x)-u(-x)}{x}}~=~\lim _{\,x\searrow 0\,}\,{\frac {u'(x)+u'(-x)}{1}}~=~2u'(0)~,}$$ since ${\displaystyle u'(x)}$ is continuous and L'Hopital's rule applies.

Therefore, ${\displaystyle \int _{0}^{1}\,{\frac {u(x)-u(-x)}{x}}\,\mathrm {d} x}$ exists and by applying the mean value theorem to ${\displaystyle u(x)-u(-x),}$ we get:

${\displaystyle \left|\,\int _{0}^{1}\,{\frac {u(x)-u(-x)}{x}}\,\mathrm {d} x\,\right|\;\leq \;\int _{0}^{1}{\frac {{\bigl |}u(x)-u(-x){\bigr |}}{x}}\,\mathrm {d} x\;\leq \;\int _{0}^{1}\,{\frac {\,2x\,}{x}}\,\sup _{x\in \mathbb {R} }\,{\Bigl |}u'(x){\Bigr |}\,\mathrm {d} x\;\leq \;2\,\sup _{x\in \mathbb {R} }\,{\Bigl |}u'(x){\Bigr |}~.}$

And furthermore:

${\displaystyle \left|\,\int _{1}^{\infty }{\frac {\;u(x)-u(-x)\;}{x}}\,\mathrm {d} x\,\right|\;\leq \;2\,\sup _{x\in \mathbb {R} }\,{\Bigl |}x\cdot u(x){\Bigr |}~\cdot \;\int _{1}^{\infty }{\frac {\mathrm {d} x}{\,x^{2}\,}}\;=\;2\,\sup _{x\in \mathbb {R} }\,{\Bigl |}x\cdot u(x){\Bigr |}~,}$

we note that the map $${\displaystyle \operatorname {p.v.} \;\left({\frac {1}{\,x\,}}\right)\,:\,{C_{c}^{\infty }}(\mathbb {R} )\to \mathbb {C} }$$ is bounded by the usual seminorms for Schwartz functions ${\displaystyle u}$. Therefore, this map defines, as it is obviously linear, a continuous functional on the Schwartz space and therefore a tempered distribution.

Note that the proof needs ${\displaystyle u}$ merely to be continuously differentiable in a neighbourhood of 0 and ${\displaystyle x\,u}$ to be bounded towards infinity. The principal value therefore is defined on even weaker assumptions such as ${\displaystyle u}$ integrable with compact support and differentiable at 0.

The principal value is the inverse distribution of the function ${\displaystyle x}$ and is almost the only distribution with this property: $${\displaystyle xf=1\quad \Leftrightarrow \quad \exists K:\;\;f=\operatorname {p.\!v.} \left({\frac {1}{x}}\right)+K\delta ,}$$ where ${\displaystyle K}$ is a constant and ${\displaystyle \delta }$ the Dirac distribution.

In a broader sense, the principal value can be defined for a wide class of singular integral kernels on the Euclidean space ${\displaystyle \mathbb {R} ^{n}}$. If ${\displaystyle K}$ has an isolated singularity at the origin, but is an otherwise "nice" function, then the principal-value distribution is defined on compactly supported smooth functions by $${\displaystyle [\operatorname {p.\!v.} (K)](f)=\lim _{\varepsilon \to 0}\int _{\mathbb {R} ^{n}\setminus B_{\varepsilon }(0)}f(x)K(x)\,\mathrm {d} x.}$$ Such a limit may not be well defined, or, being well-defined, it may not necessarily define a distribution. It is, however, well-defined if ${\displaystyle K}$ is a continuous homogeneous function of degree ${\displaystyle -n}$ whose integral over any sphere centered at the origin vanishes. This is the case, for instance, with the Riesz transforms.

Consider the values of two limits: $${\displaystyle \lim _{a\to 0+}\left(\int _{-1}^{-a}{\frac {\mathrm {d} x}{x}}+\int _{a}^{1}{\frac {\mathrm {d} x}{x}}\right)=0,}$$

This is the Cauchy principal value of the otherwise ill-defined expression $${\displaystyle \int _{-1}^{1}{\frac {\mathrm {d} x}{x}},{\text{ (which gives }}{-\infty }+\infty {\text{)}}.}$$

Also: $${\displaystyle \lim _{a\to 0+}\left(\int _{-1}^{-2a}{\frac {\mathrm {d} x}{x}}+\int _{a}^{1}{\frac {\mathrm {d} x}{x}}\right)=\ln 2.}$$

Similarly, we have $${\displaystyle \lim _{a\to \infty }\int _{-a}^{a}{\frac {2x\,\mathrm {d} x}{x^{2}+1}}=0,}$$

This is the principal value of the otherwise ill-defined expression $${\displaystyle \int _{-\infty }^{\infty }{\frac {2x\,\mathrm {d} x}{x^{2}+1}}{\text{ (which gives }}{-\infty }+\infty {\text{)}}.}$$ but $${\displaystyle \lim _{a\to \infty }\int _{-2a}^{a}{\frac {2x\,\mathrm {d} x}{x^{2}+1}}=-\ln 4.}$$

Different authors use different notations for the Cauchy principal value of a function ${\displaystyle f}$, among others: $${\displaystyle PV\int f(x)\,\mathrm {d} x,}$$ $${\displaystyle \mathrm {p.v.} \int f(x)\,\mathrm {d} x,}$$ $${\displaystyle \int _{L}^{*}f(z)\,\mathrm {d} z,}$$ $${\displaystyle -\!\!\!\!\!\!\int f(x)\,\mathrm {d} x,}$$ as well as ${\displaystyle P,}$ P.V., ${\displaystyle {\mathcal {P}},}$ ${\displaystyle P_{v},}$ ${\displaystyle (CPV),}$ ${\displaystyle {\mathcal {C}},}$ and V.P.

• Hadamard finite part integral
• Hilbert transform
• Sokhotski–Plemelj theorem