Zeros, Poles, and Isolated Singularities

Almost every function a physicist actually computes with is analytic except at a handful of points — and those points carry the physics. The Feynman propagator blows up when a particle goes on shell; a driven oscillator's response diverges at resonance; a partition function's singularities mark phase transitions; a scattering amplitude's poles are bound states and its branch cuts are thresholds. Complex analysis turns this folklore into a theorem: near an isolated singularity — a point z0z_0 where ff fails to be analytic but is analytic on some punctured disc 0<∣z−z0∣<R0<|z-z_0|<R — the function admits a convergent Laurent series, and the whole local character of ff is read off from that series' principal part, the sum of its negative-power terms. Three mutually exclusive cases exhaust the possibilities, and they are decided by counting negative powers. None: the singularity is removable and ff extends analytically across z0z_0. Finitely many, the deepest being (z−z0)−m(z-z_0)^{-m}: a pole of order mm, a clean polynomial-rate blow-up. Infinitely many: an essential singularity, where the function is so wild that in every neighbourhood it comes arbitrarily close to every complex value. Dual to poles are zeros, classified the same way by where the positive part of the series begins, and the two are tied by reciprocation with an order arithmetic so simple it becomes bookkeeping. This article builds the classification rigorously from the Laurent series, proves the three cornerstone theorems — Riemann, the pole–growth equivalence, and Casorati–Weierstrass — and sets up the residue calculus that turns all of this into computable integrals.

Key takeaways
  • An isolated singularity is a point z0z_0 where ff is non-analytic but analytic on a punctured disc 0<∣z−z0∣<R0<|z-z_0|<R, so it has a convergent Laurent series there.
  • Counting the negative powers of that Laurent series classifies the point: none ⇒ removable, a highest power (z−z0)−m(z-z_0)^{-m} ⇒ pole of order mm, infinitely many ⇒ essential.
  • A zero of order mm factors as f(z)=(z−z0)mg(z)f(z)=(z-z_0)^m g(z) with g(z0)≠0g(z_0)\neq 0; a pole of order mm factors as f(z)=h(z)/(z−z0)mf(z)=h(z)/(z-z_0)^m with h(z0)≠0h(z_0)\neq 0.
  • Poles and zeros are reciprocal — ff has a pole of order mm iff 1/f1/f has a zero of order mm — and orders add: ord(fg)=ord(f)+ord(g)\mathrm{ord}(fg)=\mathrm{ord}(f)+\mathrm{ord}(g).
  • The type is fixed by the limit of ∣f∣|f|: bounded ⇒ removable (Riemann), →∞\to\infty ⇒ pole, no limit ⇒ essential (Casorati–Weierstrass says ff approaches every value).
  • Key idea: The Laurent series in the punctured disc is the master tool — count the negative-power terms of the principal part; zero means removable, a highest term (z−z0)−m(z-z_0)^{-m} means a pole of order mm, infinitely many means essential, and zeros are the same idea read from where the positive part begins.
The core insight

The Laurent series in the punctured disc is the master tool: count the negative-power terms in the principal part. Zero negative powers ⇒ removable; a highest negative power (z−z0)−m(z-z_0)^{-m} ⇒ pole of order mm; infinitely many ⇒ essential. Zeros are the same idea applied to where the positive part begins.

Where this lives today

The pole–zero anatomy is the working language of modern spectral physics. In quantum field theory a stable particle is literally a simple pole of the propagator on the real axis of the complex energy plane, its position fixing the mass and its residue normalizing the field; an unstable particle is a pole pushed off the real axis onto the second Riemann sheet, and the size of the imaginary part sets the decay width. The Z boson, for instance, shows up in electron–positron collision data at LEP and the LHC as a Breit–Wigner resonance whose complex pole sits at mass mZ≈91.19 GeVm_Z\approx 91.19\ \mathrm{GeV} with width ΓZ≈2.495 GeV\Gamma_Z\approx 2.495\ \mathrm{GeV} — reading that pole is how the particle's lifetime is measured. The same grammar runs through condensed matter (Green's-function poles are quasiparticles, their residues the spectral weight ZZ), through control engineering (the Nyquist criterion is the argument principle counting right-half-plane poles for stability), and through the analytic SS-matrix program, where bound states, resonances, and thresholds are catalogued as poles and branch points of scattering amplitudes. Essential singularities, far from being pathological, encode non-perturbative physics: the e−1/g2e^{-1/g^2} dependence of instanton and tunnelling amplitudes on a coupling gg is exactly the essential singularity at g=0g=0 that makes perturbation theory an asymptotic, non-convergent series.

The physical picture

Think of an analytic function as, locally, a power series. A zero is a series that starts late: the lowest surviving term is (z−z0)m(z-z_0)^m, so ff presses flat against zero to mm-th order. A pole is the mirror image — the series starts before the constant term, led by (z−z0)−m(z-z_0)^{-m}, and runs off to infinity at a definite polynomial rate. The dividing line between tame and wild behaviour is a single yes/no question: do infinitely many negative powers survive? A pole diverges but can be tamed — multiply by (z−z0)m(z-z_0)^m and it becomes analytic. An essential singularity resists every finite power: no matter how many times you multiply by (z−z0)(z-z_0), infinitely many negative powers remain, and the function sweeps through the whole complex plane in every neighbourhood of z0z_0. The first figure makes this trichotomy visible by watching the modulus ∣f(z)∣|f(z)| as zz approaches z0z_0 — bounded, diverging, or neither — which is not just a heuristic but exactly the content of the three classification theorems.

Figure
The trichotomy of isolated singularities read from the modulus of f as z → z₀: a removable point stays bounded (sin z / z → 1), a pole of order m diverges like 1/rᵐ → ∞, and an essential singularity has no limit at all — e^(1/z) sweeps every value in every neighbourhood
One question decides everything: what does |f(z)| do as z → z₀? removable sin z / z → 1 limit finite z₀ pole of order m |f| ≈ |a₋ₘ| / rᵐ → ∞ z₀ essential no limit — every value z₀ Bounded ⇒ removable (Riemann) · diverges cleanly ⇒ pole · neither ⇒ essential (Casorati–Weierstrass)

A useful physical translation: poles are resonances — a denominator vanishing — and the order counts the degree of degeneracy or the multiplicity of the coincident modes. Essential singularities, by contrast, are the mark of genuinely non-perturbative behaviour: e−1/z2e^{-1/z^2} near 00 is the textbook function that is infinitely differentiable on the real axis yet has no convergent Taylor series, the reason asymptotic (divergent) series pervade perturbation theory. The bookkeeping that ties zeros and poles together — orders adding under products, subtracting under quotients — is not just algebra: integrated around a contour it becomes the argument principle, a machine that literally counts zeros minus poles. The second figure shows that machine on a concrete meromorphic function.

Figure
Pole–zero map of f(z) = (z−i)² / (z (z+2)³): a double zero at z = i weighs Z = 2 against a simple pole at 0 and a triple pole at −2 weighing P = 4, so the argument principle gives (1/2πi) ∮ f′/f dz = Z − P = 2 − 4 = −2 for any contour enclosing all three
Orders add up: the argument principle counts zeros minus poles inside a contour Re z Im z contour C −2 i zero, order 2 z = i pole, order 1 z = 0 pole, order 3 z = −2 The count Z = 2 (double zero at i) P = 1 + 3 = 4 (poles) (1/2πi) ∮ f′/f dz = Z − P = −2

The core equations

Everything begins with the Laurent series. On the punctured disc where ff is analytic, it is represented by a two-sided power series whose coefficients are contour integrals; the negative-index coefficients are precisely the data that classify the singularity.

f(z)=∑n=−∞∞an(z−z0)n,an=12πi∮Cf(ζ)(ζ−z0)n+1 dζ(1)f(z)=\sum_{n=-\infty}^{\infty} a_n (z-z_0)^n,\qquad a_n=\frac{1}{2\pi i}\oint_{C}\frac{f(\zeta)}{(\zeta-z_0)^{n+1}}\,\mathrm{d}\zeta \qquad (1)

where: z0z_0 — the isolated singular point; ana_n — the Laurent coefficient; CC — a positively oriented circle lying in 0<∣z−z0∣<R0<|z-z_0|<R. Reading: the coefficients are unique, and the sum splits into an analytic part (n≥0n\ge 0) plus the principal part (n<0n<0). The principal part is the whole classifier. Isolating the negative powers and asking how many are nonzero is all three definitions at once.

∑n=1∞a−n(z−z0)−n⏟principal part(2)\underbrace{\sum_{n=1}^{\infty} a_{-n}(z-z_0)^{-n}}_{\text{principal part}} \qquad (2)

where: a−na_{-n} — the coefficients of the negative powers. Reading: all a−n=0a_{-n}=0 ⇒ removable; a highest nonzero a−ma_{-m} ⇒ pole of order mm; infinitely many nonzero ⇒ essential. A zero is classified by where the positive part starts. Vanishing to order mm is the same as factoring out (z−z0)m(z-z_0)^m from an analytic function that survives at z0z_0.

f(z)=(z−z0)m g(z),g(z0)≠0, g analytic(3)f(z)=(z-z_0)^m\,g(z),\qquad g(z_0)\neq 0,\ g\ \text{analytic} \qquad (3)

where: mm — the order of the zero, a positive integer; gg — the analytic cofactor, nonzero at z0z_0. Reading: equivalently f(z0)=⋯=f(m−1)(z0)=0f(z_0)=\cdots=f^{(m-1)}(z_0)=0 but f(m)(z0)≠0f^{(m)}(z_0)\neq 0. The pole is the exact reciprocal statement: a finite principal part ending at order mm is the same as an analytic numerator over (z−z0)m(z-z_0)^m.

f(z)=h(z)(z−z0)m,h(z0)≠0, h analytic(4)f(z)=\frac{h(z)}{(z-z_0)^m},\qquad h(z_0)\neq 0,\ h\ \text{analytic} \qquad (4)

where: hh — analytic and nonzero at z0z_0; mm — the order of the pole. Reading: h(z0)=a−mh(z_0)=a_{-m}, the leading Laurent coefficient. The two structures are locked together by taking reciprocals, which is why zeros and poles are never studied apart.

f has a pole of order m at z0  ⟺  1f has a zero of order m at z0(5)f\ \text{has a pole of order } m \text{ at } z_0 \iff \frac{1}{f}\ \text{has a zero of order } m \text{ at } z_0 \qquad (5)

Because orders live in the exponent, they simply add and subtract when functions are multiplied and divided — this is the divisor calculus of a meromorphic function.

ordz0(fg)=ordz0(f)+ordz0(g),ordz0(f/g)=ordz0(f)−ordz0(g)(6)\mathrm{ord}_{z_0}(fg)=\mathrm{ord}_{z_0}(f)+\mathrm{ord}_{z_0}(g),\qquad \mathrm{ord}_{z_0}(f/g)=\mathrm{ord}_{z_0}(f)-\mathrm{ord}_{z_0}(g) \qquad (6)

where: ordz0(f)\mathrm{ord}_{z_0}(f) — the signed integer order: +m+m at a zero, −m-m at a pole, 00 at an ordinary analytic point where f≠0f\neq 0. The essential case has its own signature theorem: near such a point the image of every punctured neighbourhood is dense in the whole plane.

f({0<∣z−z0∣<δ})‾=Cfor every δ>0(7)\overline{f\big(\{0<|z-z_0|<\delta\}\big)}=\mathbb{C}\quad\text{for every } \delta>0 \qquad (7)

where: δ\delta — the radius of an arbitrarily small punctured disc. Reading: no matter how close you zoom, ff still comes arbitrarily near every complex number. Finally, the single Laurent coefficient that survives contour integration is the residue, and at a pole of order mm it is extracted by differentiation — the formula that powers the entire residue theorem.

Res⁡z=z0f=a−1=lim⁡z→z01(m−1)!dm−1dzm−1[(z−z0)mf(z)](8)\operatorname*{Res}_{z=z_0} f = a_{-1} = \lim_{z\to z_0}\frac{1}{(m-1)!}\frac{\mathrm{d}^{m-1}}{\mathrm{d}z^{m-1}}\Big[(z-z_0)^m f(z)\Big] \qquad (8)

where: a−1a_{-1} — the coefficient of (z−z0)−1(z-z_0)^{-1}. Reading: for a simple pole (m=1m=1) this is just lim⁡z→z0(z−z0)f(z)\lim_{z\to z_0}(z-z_0)f(z).

Derivation ARiemann's removable singularity theorem

The goal is to prove the least obvious of the three characterizations: a bounded isolated singularity is always removable. Concretely, if ff is analytic on the punctured disc 0<∣z−z0∣<R0<|z-z_0|<R and merely stays bounded there, ∣f(z)∣≤M|f(z)|\le M, then not only does the limit at z0z_0 exist — every coefficient of the principal part vanishes identically, so ff was secretly analytic all along and one value at z0z_0 repairs it. The proof is a single sharp estimate applied to the Laurent coefficients and then a limit.

Step 1. Write the negative coefficients from the Laurent formula (1) on a circle of radius ρ<R\rho<R: a−n=12πi∮∣ζ−z0∣=ρf(ζ)(ζ−z0) n−1 dζa_{-n}=\frac{1}{2\pi i}\oint_{|\zeta-z_0|=\rho} f(\zeta)(\zeta-z_0)^{\,n-1}\,\mathrm{d}\zeta. Why: this expresses each principal-part coefficient purely in terms of values of ff on a circle we are free to shrink toward z0z_0.

Step 2. Apply the ML (estimation) inequality — integrand bounded by Mρ n−1M\rho^{\,n-1}, contour length 2πρ2\pi\rho: ∣a−n∣≤12π⋅2πρ⋅Mρ n−1=Mρ n|a_{-n}|\le \frac{1}{2\pi}\cdot 2\pi\rho\cdot M\rho^{\,n-1}=M\rho^{\,n}. Why: boundedness is the only hypothesis, and here is where it enters — it caps the integrand uniformly on the circle regardless of how small ρ\rho is.

Step 3. Let the circle shrink. For every n≥1n\ge 1 the exponent is positive, so ∣a−n∣≤Mρ n→0|a_{-n}|\le M\rho^{\,n}\to 0 as ρ→0\rho\to 0; since a−na_{-n} does not depend on ρ\rho, it must be exactly 00. Why: a fixed number bounded by something that goes to zero is zero — the entire principal part is annihilated.

Step 4. What remains is f(z)=∑n≥0an(z−z0)nf(z)=\sum_{n\ge 0}a_n(z-z_0)^n, an ordinary convergent power series; defining f(z0):=a0f(z_0):=a_0 makes ff analytic on the full disc. Why: a function equal to a convergent power series on a punctured disc extends to its centre, so the singularity was removable — a hole in the domain, not in the function. *The key check is reproducing the ML bound ∣a−n∣≤Mρ n|a_{-n}|\le M\rho^{\,n} and taking ρ→0\rho\to 0.*

Go deeper · Derivation B · 1 minPole order equals growth rate, and the reciprocal zero
Derivation BPole order equals growth rate, and the reciprocal zero

The goal is to weld together three descriptions of a pole that a physicist uses interchangeably: the algebraic one (a Laurent series whose deepest term is (z−z0)−m(z-z_0)^{-m}), the analytic one (the operational limit lim⁡z→z0(z−z0)mf=a−m≠0\lim_{z\to z_0}(z-z_0)^m f=a_{-m}\neq 0 that lets you read off an order by hand), and the dual one (1/f1/f has a genuine zero of the same order mm). Establishing these as one statement is what makes order arithmetic (6) and pole–zero duality (5) legitimate rather than mnemonic.

Step 1. Start from the definition of a pole of order mm: the Laurent series begins at n=−mn=-m, so f(z)=a−m(z−z0)m+⋯+a−1z−z0+∑n≥0an(z−z0)nf(z)=\dfrac{a_{-m}}{(z-z_0)^m}+\cdots+\dfrac{a_{-1}}{z-z_0}+\sum_{n\ge 0}a_n(z-z_0)^n with a−m≠0a_{-m}\neq 0. Why: this is the precise meaning of "finitely many negative powers, deepest one of order mm", and it is the only assumption we make.

Step 2. Factor (z−z0)−m(z-z_0)^{-m} out of the whole series: f(z)=(z−z0)−m[a−m+a−m+1(z−z0)+⋯]=(z−z0)−mh(z)f(z)=(z-z_0)^{-m}\big[a_{-m}+a_{-m+1}(z-z_0)+\cdots\big]=(z-z_0)^{-m}h(z). Why: the bracketed series has no negative powers and converges, so it defines an analytic hh with h(z0)=a−m≠0h(z_0)=a_{-m}\neq 0 — exactly form (4).

Step 3. Multiply by (z−z0)m(z-z_0)^m and take the limit: lim⁡z→z0(z−z0)mf(z)=h(z0)=a−m≠0\lim_{z\to z_0}(z-z_0)^m f(z)=h(z_0)=a_{-m}\neq 0. Why: this is the practical test — if (z−z0)kf(z-z_0)^k f tends to a finite nonzero limit for some integer kk, then kk is the pole order, no series expansion required.

Step 4. Take the reciprocal: 1f(z)=(z−z0)m1h(z)\dfrac{1}{f(z)}=(z-z_0)^m\dfrac{1}{h(z)}. Since h(z0)≠0h(z_0)\neq 0, 1/h1/h is analytic and nonzero at z0z_0, so 1/f1/f vanishes to exactly order mm. Why: this is pole–zero duality (5) proved outright, and it also confirms that ∣f∣→∞|f|\to\infty at a pole because 1/f→0.1/f\to 0.

Go deeper · Derivation C · 1 minThe essential singularity of e^(1/z) and Casorati–Weierstrass
Derivation CThe essential singularity of e^(1/z) and Casorati–Weierstrass

The goal is to exhibit the essential case in the flesh and to see why it is qualitatively unlike a pole. We show that e1/ze^{1/z} has infinitely many negative Laurent powers at 00, that its one-sided real limits flatly contradict each other so no value in the extended plane can be assigned, and — the punchline — that it hits (almost) every target value infinitely often as z→0z\to 0, an explicit instance of Casorati–Weierstrass sharpened to Picard.

Step 1. Expand the exponential and substitute w=1/zw=1/z: ew=∑n=0∞wnn!e^{w}=\sum_{n=0}^{\infty}\dfrac{w^n}{n!} converges for all ww, so e1/z=∑n=0∞1n! z−n=1+1z+12! z2+13! z3+⋯e^{1/z}=\sum_{n=0}^{\infty}\dfrac{1}{n!}\,z^{-n}=1+\dfrac1z+\dfrac{1}{2!\,z^2}+\dfrac{1}{3!\,z^3}+\cdots for every z≠0z\neq 0. Why: this is the Laurent series about 00 directly, no contour integrals needed, because the exponential series is entire.

Step 2. Read the principal part: a−n=1/n!≠0a_{-n}=1/n!\neq 0 for all n≥0n\ge 0. Why: infinitely many negative coefficients are nonzero — by definition (2), an essential singularity, and by contrast with a pole, no finite power of zz clears them all.

Step 3. Probe the modulus along the real axis: as z→0+z\to 0^+, e1/z→+∞e^{1/z}\to +\infty (e.g. e1/0.25=54.6e^{1/0.25}=54.6), while as z→0−z\to 0^-, e1/z→0+e^{1/z}\to 0^+ (e.g. e−1/0.25=0.0183e^{-1/0.25}=0.0183). Why: the two one-sided limits disagree, so no single value — finite or ∞\infty — works, ruling out both removable and pole and leaving only essential.

Step 4. Solve e1/z=ce^{1/z}=c for any c≠0c\neq 0: 1z=log⁡c+2πik\tfrac1z=\log c+2\pi i k, so zk=1log⁡c+2πik→0z_k=\dfrac{1}{\log c+2\pi i k}\to 0 as k→∞k\to\infty. Why: every nonzero target is attained infinitely often at points accumulating on 00 — a concrete Casorati–Weierstrass, and in fact Picard's great theorem, with c=0c=0 the single omitted value. *The three concrete checks are: the Laurent coefficients a−n=1/n!a_{-n}=1/n!, the one-sided numeric limits e4=54.6e^{4}=54.6 and e−4=0.0183e^{-4}=0.0183, and the root formula zk=1/(log⁡c+2πik)z_k=1/(\log c+2\pi i k).*

Limiting cases

The three types, plus the point at infinity, are distinguished entirely by the behaviour of ff as z→z0z\to z_0, which is the operational face of the Laurent classification:

CaseResultMeaning
Removablelim⁡z→z0f(z)\lim_{z\to z_0} f(z) exists and is finitePrincipal part vanishes; ff extends analytically. Example: sin⁡zz\dfrac{\sin z}{z} at 00, limit 11.
Pole of order mmlim⁡z→z0(z−z0)mf(z)=a−m≠0\lim_{z\to z_0}(z-z_0)^m f(z)=a_{-m}\neq 0, and ∣f(z)∣→∞|f(z)|\to\inftyFinite principal part ending at (z−z0)−m(z-z_0)^{-m}. Example: 1(z−1)2\dfrac{1}{(z-1)^2}, a double pole at z=1z=1.
Essentiallim⁡z→z0f(z)\lim_{z\to z_0} f(z) does not exist in C^\widehat{\mathbb{C}}Infinitely many negative powers; by Casorati–Weierstrass ff nears every value. Example: e1/ze^{1/z} at 00.
At infinityclassify g(w)=f(1/w)g(w)=f(1/w) at w=0w=0A degree-nn polynomial has a pole of order nn at ∞\infty; eze^z has an essential singularity at ∞\infty.

Worked examples

1 · Classifying and residue-computing a pole (fully worked)

Classify the singularity of f(z)=cos⁡zz3f(z)=\dfrac{\cos z}{z^3} at z=0z=0, and find its residue two independent ways.

Step 1. Expand the numerator and divide, which reads the Laurent series off directly: cos⁡z=1−z22+z424−⋯\cos z = 1-\dfrac{z^2}{2}+\dfrac{z^4}{24}-\cdots, so f(z)=1z3−12 1z+z24−⋯f(z)=\dfrac{1}{z^3}-\dfrac{1}{2}\,\dfrac{1}{z}+\dfrac{z}{24}-\cdots. The deepest negative power is z−3z^{-3} with nonzero coefficient, so this is a pole of order 3.

Step 2. Read the residue as the coefficient of z−1z^{-1}: the only z−1z^{-1} term is −12z−1-\tfrac12 z^{-1}, so Res⁡z=0f=−12\operatorname{Res}_{z=0}f=\mathbf{-\tfrac12}. Note the z−2z^{-2} coefficient is 00 — a pole of order 3 need not have every intermediate term.

Step 3. Cross-check with the order-mm residue formula (8) using m=3m=3: Res⁡=12!lim⁡z→0d2dz2[z3f(z)]=12lim⁡z→0d2dz2cos⁡z=12lim⁡z→0(−cos⁡z)=−12\operatorname{Res}=\dfrac{1}{2!}\lim_{z\to0}\dfrac{\mathrm{d}^2}{\mathrm{d}z^2}\big[z^3 f(z)\big]=\dfrac12\lim_{z\to0}\dfrac{\mathrm{d}^2}{\mathrm{d}z^2}\cos z=\dfrac12\lim_{z\to0}(-\cos z)=\mathbf{-\tfrac12}. The two routes agree exactly. Pole of order 3, residue = −1/2.

2 · Reading orders through a product and quotient

Find and classify every isolated singularity of g(z)=z2+4(z2+1)(ez−1)g(z)=\dfrac{z^2+4}{(z^2+1)(e^z-1)} near the origin, and compute the residue at z=0z=0.

Step 1. Locate candidates where the denominator vanishes: ez−1=0e^z-1=0 at z=0z=0 (and z=2πikz=2\pi i k), and z2+1=0z^2+1=0 at z=±iz=\pm i. Near the origin the relevant points are z=0z=0 and z=iz=i.

Step 2. Use order arithmetic (6) at z=0z=0. The numerator z2+4→4≠0z^2+4\to 4\neq 0 (order 00); z2+1→1≠0z^2+1\to 1\neq 0 (order 00); ez−1e^z-1 has a simple zero (since ez−1=z+z2/2+⋯e^z-1=z+z^2/2+\cdots, order +1+1). So ord0(g)=0−0−1=−1\mathrm{ord}_0(g)=0-0-1=-1: a simple pole.

Step 3. At z=iz=i: z2+1=(z−i)(z+i)z^2+1=(z-i)(z+i) has a simple zero (order +1+1); z2+4→3≠0z^2+4\to 3\neq 0; ez−1→ei−1≠0e^z-1\to e^i-1\neq 0. Hence ordi(g)=0−1−0=−1\mathrm{ord}_i(g)=0-1-0=-1: another simple pole.

Step 4. Residue at the simple pole z=0z=0 via (8) with m=1m=1: Res⁡z=0g=lim⁡z→0z g(z)=z2+4z2+1⋅lim⁡z→0zez−1=41⋅1=4\operatorname{Res}_{z=0}g=\lim_{z\to0} z\,g(z)=\dfrac{z^2+4}{z^2+1}\cdot\lim_{z\to0}\dfrac{z}{e^z-1}=\dfrac{4}{1}\cdot 1=\mathbf{4}, using the standard limit z/(ez−1)→1z/(e^z-1)\to 1. Two simple poles (at 0 and i); Res(g, 0) = 4.

Physics of the portrait

The classification of singularities is where nineteenth-century analysis became rigorous. The names attached to its theorems — Laurent, Riemann, Casorati, Weierstrass, Picard — trace the passage from a formal calculus of series to a precise anatomy of how functions can fail to be analytic, and no one did more to enforce that rigour than Weierstrass.

Origins · 1 minKarl Weierstrass and how this was discovered
KW
Karl Weierstrass
1815–1897 · the arithmetization of analysis and the theory of essential singularities

Karl Weierstrass came to mathematics late and by an unusual road — for fourteen years he was a schoolteacher in remote Prussian towns, teaching gymnastics and penmanship alongside mathematics while writing research he could barely publish, until an 1854 paper on Abelian functions brought him a sudden honorary doctorate and a chair in Berlin. There he built the modern foundations of analysis brick by rigorous brick, insisting that limits, continuity, and convergence rest on explicit inequalities rather than geometric intuition — the famous ε\varepsilon–δ\delta discipline is his legacy, and his construction of a continuous but nowhere-differentiable function warned a generation that intuition alone is treacherous. In complex analysis he championed the power series as the primary object, defining an analytic function by its local expansions and analytic continuation rather than by Cauchy's integrals or Riemann's geometry. It was in this program that he sharpened the notion of an essential singularity: the theorem that a function comes arbitrarily close to every value near such a point, which he proved in 1876, is now shared with Felice Casorati (who published a version in 1868) and was soon deepened by Émile Picard. Weierstrass's students — Schwarz, Frobenius, Mittag-Leffler, and Sofia Kovalevskaya, whom he tutored privately when the university barred her — carried his standard of rigour into the twentieth century, and the classification of poles and essential singularities in this article is essentially his framework.

History
1843
— Laurent presents the two-sided series expansion valid in an annulus
1851
— Riemann's dissertation founds geometric function theory and the removable-singularity viewpoint
1868
— Casorati states the density theorem for essential singularities
1876
— Weierstrass publishes his theory of essential singularities and analytic continuation
1879
— Picard proves the great theorem: near an essential singularity every value but one is taken infinitely often

Going further

Texts. Stone & Goldbart, Mathematics for Physics, §17–18 (Function Theory) · Arfken, Weber & Harris, Mathematical Methods for Physicists, 7th ed., Ch. 11 · Hassani, Mathematical Physics, Ch. 10–11 · Byron & Fuller, Mathematics of Classical and Quantum Physics, Ch. 6 · Riley, Hobson & Bence, Mathematical Methods for Physics and Engineering, 3rd ed., Ch. 24 · Ahlfors, Complex Analysis, Ch. 4–5 · Conway, Functions of One Complex Variable, Ch. V. Key numerical anchors in this article: the Riemann ML estimate ∣a−n∣≤Mρ n|a_{-n}|\le M\rho^{\,n}; the order-mm residue formula (8); the Laurent coefficients a−n=1/n!a_{-n}=1/n! for e1/ze^{1/z}; order arithmetic (6); and the results Res⁡(cos⁡z/z3,0)=−12\operatorname{Res}(\cos z/z^3,0)=-\tfrac12, Res⁡(g,0)=4\operatorname{Res}(g,0)=4, sin⁡z/z→1\sin z/z\to 1, (cos⁡z−1)/z2→−12(\cos z-1)/z^2\to -\tfrac12, e4=54.6e^{4}=54.6, e−4=0.0183e^{-4}=0.0183. The particle-physics values mZ=91.19 GeVm_Z=91.19\ \mathrm{GeV} and ΓZ=2.495 GeV\Gamma_Z=2.495\ \mathrm{GeV} are standard PDG values.

Prerequisites. Functions of a Complex Variable and Limits · Analytic Functions and the Cauchy–Riemann Equations · Taylor Series of Analytic Functions · Laurent Series. Leads to. Residues and the Residue Theorem · Consequences: Liouville, Maximum Modulus, Morera · Cauchy's Integral Formula and Derivatives.

Ask the tutor: Hand me any function and a point, and I will walk you through counting its Laurent principal part to classify the singularity and read off the residue.

What actually makes a singularity "isolated", and why does it matter so much?

A singularity at z0z_0 is isolated when there is some punctured disc 0<∣z−z0∣<R0<|z-z_0|<R on which ff is analytic — the bad point sits alone, with a clean analytic ring around it. This matters because that ring is exactly the annulus on which a Laurent series is guaranteed to converge, and the Laurent series is the entire classification engine. Without isolation there is no Laurent expansion, so the removable/pole/essential trichotomy simply does not apply. Branch points (like 00 for log⁡z\log z) and accumulation points of poles (like 00 for 1/sin⁡(1/z)1/\sin(1/z)) are the standard non-isolated singularities: they are genuinely singular, but of a different species that needs branch cuts or other tools rather than a principal-part count.

How do I tell a pole from an essential singularity in practice?

Fastest test: look at ∣f(z)∣|f(z)| as z→z0z\to z_0. If it tends to +∞+\infty from every direction, you have a pole; find the order by the smallest integer mm with (z−z0)mf(z-z_0)^m f tending to a finite nonzero limit. If ∣f∣|f| does something inconsistent — blowing up along one approach and staying bounded or going to zero along another — it is essential. The cleanest confirmation is to write the Laurent series and count: a finite number of negative powers is a pole (the deepest one gives the order), infinitely many is essential. The exponential-of-a-reciprocal pattern, like e1/ze^{1/z} or sin⁡(1/z)\sin(1/z), is the classic tell for an essential singularity.

Why are zeros and poles always discussed together?

Because reciprocation swaps them exactly: ff has a pole of order mm at z0z_0 if and only if 1/f1/f has a zero of order mm there, and vice versa. This duality means every fact about zeros has a mirror fact about poles for free, and it lets you compute one by studying the other — often it is easier to find where a clean analytic function vanishes than to analyze where its reciprocal blows up. The unifying device is the signed order ordz0\mathrm{ord}_{z_0}: positive for zeros, negative for poles, additive under multiplication. Physically, a transfer function's zeros (where output vanishes) and poles (resonances) are read from the same diagram for exactly this reason.

Is a removable singularity a real singularity at all?

Not really — it is an artifact of how the function was written, a hole in the domain rather than in the function. When ff is bounded near z0z_0, Riemann's theorem guarantees that a single well-chosen value at z0z_0 makes ff fully analytic, so the "singularity" was only ever a 0/00/0 bookkeeping gap. The canonical case is sin⁡z/z\sin z / z at 00: define it to be 11 there and it becomes the perfectly smooth entire function sinc⁡\operatorname{sinc}. Removable singularities matter precisely because they are removable — they let you cancel apparent divergences, and recognizing one prevents you from mistaking a superficial 0/00/0 for genuine singular behaviour.

Where do poles actually show up in physics, concretely?

Everywhere a denominator can vanish. A driven damped oscillator's response 1/(ω02−ω2−iγω)1/(\omega_0^2-\omega^2-i\gamma\omega) has poles at the complex resonance frequencies, and the pole locations are the natural modes. In quantum field theory the propagator 1/(p2−m2+iϵ)1/(p^2-m^2+i\epsilon) has a pole whose position is the particle's mass and whose residue normalizes the field; unstable particles are poles off the real axis, with the imaginary part fixing the decay width. In circuit theory and control, the poles of a transfer function determine stability and ringing. And inverse Laplace transforms reduce to summing residues at the poles of the transformed function — the whole method of solving linear ODEs by transforms is pole bookkeeping.

What is Casorati–Weierstrass really saying, and how does Picard sharpen it?

Casorati–Weierstrass says that near an essential singularity the function is maximally erratic: in every punctured neighbourhood, no matter how tiny, the image of ff is dense in the whole complex plane — you can get arbitrarily close to any complex number you name. Picard's great theorem strengthens "arbitrarily close" to "actually equal, infinitely often": in every neighbourhood ff attains every complex value infinitely many times, with at most a single exceptional value omitted. For e1/ze^{1/z} near 00 that lone exception is 00 itself (an exponential is never zero), while every other value cc is hit at the infinitely many points zk=1/(log⁡c+2πik)z_k=1/(\log c+2\pi i k). This is the precise sense in which an essential singularity is not just "worse than a pole" but categorically wilder.

Before this

Where this leads