# Sokhotski–Plemelj theorem

The Sokhotski–Plemelj theorem (Polish spelling is Sochocki) is a theorem in complex analysis, which helps in evaluating certain integrals. The real-line version of it (see below) is often used in physics, although rarely referred to by name. The theorem is named after Julian Sochocki, who proved it in 1868, and Josip Plemelj, who rediscovered it as a main ingredient of his solution of the Riemann–Hilbert problem in 1908.

## Statement of the theorem

Let C be a smooth closed simple curve in the plane, and $\varphi$ an analytic function on C. Note that the Cauchy-type integral
    $\phi (z)={\frac {1}{2\pi i}}\int _{C}{\frac {\varphi (\zeta )\,d\zeta }{\zeta -z}},$
cannot be evaluated for any z on the curve C. However, on the interior and exterior of the curve, the integral produces analytic functions, which will be denoted $\phi _{i}$ inside C and $\phi _{e}$ outside. The Sokhotski–Plemelj formulas relate the limiting boundary values of these two analytic functions at a point z on C and the Cauchy principal value ${\mathcal {P}}$ of the integral:
    $\lim _{w\to z}\phi _{i}(w)={\frac {1}{2\pi i}}{\mathcal {P}}\int _{C}{\frac {\varphi (\zeta )\,d\zeta }{\zeta -z}}+{\frac {1}{2}}\varphi (z),$
    $\lim _{w\to z}\phi _{e}(w)={\frac {1}{2\pi i}}{\mathcal {P}}\int _{C}{\frac {\varphi (\zeta )\,d\zeta }{\zeta -z}}-{\frac {1}{2}}\varphi (z).$
Subsequent generalizations relax the smoothness requirements on curve C and the function φ.

## Version for the real line

Especially important is the version for integrals over the real line.
    $\lim _{\varepsilon \to 0^{+}}{\frac {1}{x\pm i\varepsilon }}=\mp i\pi \delta (x)+{\mathcal {P}}{{\Big{\frac {1}{x}}{\Big )}}.$
where $\delta (x)$ is the Dirac delta function where ${\mathcal {P}}$ denotes the Cauchy principal value. One may take the difference of these two equalities to obtain
    $\lim _{\varepsilon \to 0^{+}}\left[{\frac {1}{x-i\varepsilon }}-{\frac {1}{x+i\varepsilon }}\right]=2\pi i\delta (x).$
These formulae should be interpreted as integral equalities, as follows: Let f be a complex-valued function which is defined and continuous on the real line, and let a and b be real constants with $a<0<b$. Then
    $\lim _{\varepsilon \to 0^{+}}\int _{a}^{b}{\frac {f(x)}{x\pm i\varepsilon }}\,dx=\mp i\pi f(0)+{\mathcal {P}}\int _{a}^{b}{\frac {f(x)}{x}}\,dx$
and
    $\lim _{\varepsilon \to 0^{+}}\int _{a}^{b}\left[{\frac {f(x)}{x-i\varepsilon }}-{\frac {f(x)}{x+i\varepsilon }}\right]\,dx=2\pi if(0)$
Note that this version makes no use of analyticity.

## Proof of the real version

A simple proof is as follows.
    $\lim _{\varepsilon \to 0^{+}}\int _{a}^{b}{\frac {f(x)}{x\pm i\varepsilon }}\,dx=\mp i\pi \lim _{\varepsilon \to 0^{+}}\int _{a}^{b}{\frac {\varepsilon }{\pi (x^{2}+\varepsilon ^{2})}}f(x)\,dx+\lim _{\varepsilon \to 0^{+}}\int _{a}^{b}{\frac {x^{2}}{x^{2}+\varepsilon ^{2}}}\,{\frac {f(x)}{x}}\,dx.$
For the first term, ε⁄π(x$^2$ + ε$^2$) is a nascent delta function, and therefore approaches a Dirac delta function in the limit. Therefore, the first term equals ∓iπ f(0).
For the second term, the factor x$^2$⁄(x$^2$ + ε$^2$) approaches 1 for |x| ≫ ε, approaches 0 for |x| ≪ ε, and is exactly symmetric about 0. Therefore, in the limit, it turns the integral into a Cauchy principal value integral.

## Physics application

In quantum mechanics and quantum field theory, one often has to evaluate integrals of the form
    $\int _{-\infty }^{\infty }dE\,\int _{0}^{\infty }dt\,f(E)\exp(-iEt)$
where E is some energy and t is time. This expression, as written, is undefined (since the time integral does not converge), so it is typically modified by adding a negative real term to -iEt in the exponential, and then taking that to zero, i.e.:
    $\lim _{\varepsilon \to 0^{+}}\int _{-\infty }^{\infty }dE\,\int _{0}^{\infty }dt\,f(E)\exp(-iEt-\varepsilon t)=-i\lim _{\varepsilon \to 0^{+}}\int _{-\infty }^{\infty }{\frac {f(E)}{E-i\varepsilon }}\,dE=\pi f(0)-i{\mathcal {P}}\int _{-\infty }^{\infty }{\frac {f(E)}{E}}\,dE,$
where the latter step uses the real version of the theorem.

### Heitler function

In theoretical quantum optics, the derivation of a master equation in Lindblad form often requires the following integral function, which is a direct consequence of the Sokhotski–Plemelj theorem and is often called the Heitler-function:
    $\int _{0}^{\infty }d\tau \,\exp(-i(\omega \pm \nu )\tau )=\pi \delta (\omega \pm \nu )-i{\mathcal {P}}{\Big{\frac {1}{\omega \pm \nu }}{\Big )}$