Laurent Series

What kind of power series can represent a function like e1/ze^{1/z}, which blows up at the origin? A Taylor series cannot do it. Taylor series only capture functions that are actually analytic at the center point. And e1/ze^{1/z} is not even defined at z=0z=0. Yet e1/ze^{1/z} is holomorphic everywhere else. We would like an expansion valid on the punctured neighborhood 0<∣z∣<∞0<|z|<\infty. The resolution is to allow negative powers of zz alongside the usual nonnegative ones. This produces a two-sided series ∑n=−∞∞an(z−z0)n\sum_{n=-\infty}^{\infty} a_n(z-z_0)^n. This is the Laurent series, named for Pierre Alphonse Laurent, who introduced it in 1843. It is the single most important computational and conceptual tool for understanding the local behavior of holomorphic functions. Specifically, it illuminates their behavior near points where they fail to be analytic. The series splits naturally into a principal part (the negative-power terms, encoding the singular behavior) and an analytic part (the nonnegative-power terms, behaving just like an ordinary Taylor series). The coefficient a−1a_{-1} of the term (z−z0)−1(z-z_0)^{-1} turns out to carry enormous weight. It is the residue, and its extraction from a Laurent expansion is the seed of the residue calculus. The big picture is this: while Taylor series expand functions on disks, Laurent series expand functions on annuli. These are regions of the form r<∣z−z0∣<Rr<|z-z_0|<R. Every function holomorphic on such an annulus has a unique Laurent expansion converging there. Conversely, every convergent two-sided power series defines a holomorphic function on its annulus of convergence. This article proves that existence-and-uniqueness theorem in full. It works two asymmetric examples showing how expansions change annulus-to-annulus for a single function. It also shows precisely why each hypothesis is indispensable: holomorphy, the correct domain, and its punctured/annular shape. The article closes by setting up the classification of isolated singularities and the residue theorem. Both of these are read directly off the Laurent coefficients.

Key takeaways
  • A Laurent series is a two-sided series ∑n=−∞∞an(z−z0)n\sum_{n=-\infty}^\infty a_n(z-z_0)^n, splitting into a principal part (n<0n<0) and an analytic part (n≥0n\ge 0).
  • Every function holomorphic on an annulus r<∣z−z0∣<Rr<|z-z_0|<R has a unique Laurent expansion converging there.
  • Coefficients are given by contour integrals an=12πi∮f(w)(w−z0)n+1 dwa_n = \frac{1}{2\pi i}\oint \frac{f(w)}{(w-z_0)^{n+1}}\,dw over any circle in the annulus.
  • The same function can have different Laurent series on different annuli around the same point — the annulus is part of the data.
  • The coefficient a−1a_{-1} is the residue, the key quantity in the future Residue Theorem.
  • Key idea: decompose ff into a piece holomorphic inside the annulus and a piece holomorphic outside, then Taylor-expand each.

Where this lives today

The Laurent expansion is the computational backbone of the residue theorem, which reduces contour integrals to reading off a single coefficient (a−1a_{-1}) from a local expansion — this is how physicists and engineers evaluate real integrals like ∫−∞∞dx1+x4\int_{-\infty}^\infty \frac{dx}{1+x^4} via complex contours. It underlies the classification of isolated singularities (removable, pole, essential) purely by counting negative-power terms, which in turn controls the behavior of meromorphic functions used throughout analytic number theory (e.g. the Laurent expansion of the Riemann zeta function at its pole s=1s=1) and in the study of generating functions in combinatorics. In physics, Laurent expansions around poles of propagators and Green's functions extract physically meaningful quantities (masses, coupling residues) in quantum field theory and statistical mechanics. In engineering, the z-transform of a discrete-time signal is literally a Laurent series in z−1z^{-1}, and questions of stability and causality of digital filters translate directly into questions about the annulus of convergence of that series. In numerical analysis, Laurent/Fourier-type expansions on annuli underlie fast algorithms for evaluating functions on circles and for constructing conformal maps.

Ask the tutor: How does the number of nonzero negative-power Laurent coefficients determine whether a singularity is removable, a pole, or essential?

DefinitionLaurent Series

Let z0∈Cz_0 \in \mathbb{C} and let 0≤r<R≤∞0 \le r < R \le \infty. A Laurent series centered at z0z_0 is a formal bi-infinite series

∑n=−∞∞an(z−z0)n,an∈C,\sum_{n=-\infty}^{\infty} a_n (z-z_0)^n, \qquad a_n \in \mathbb{C},

understood as the sum of two ordinary series,

∑n=0∞an(z−z0)n(the analytic part)and∑n=1∞a−n(z−z0)−n(the principal part).\sum_{n=0}^{\infty} a_n(z-z_0)^n \quad \text{(the analytic part)} \qquad \text{and} \qquad \sum_{n=1}^{\infty} a_{-n}(z-z_0)^{-n} \quad \text{(the principal part).}

The Laurent series is said to converge at a point zz if both of these one-sided series converge there (as ordinary series), and its sum is the sum of the two limits. The domain of convergence of a Laurent series is an open annulus

A(z0;r,R)={ z∈C:r<∣z−z0∣<R },A(z_0; r, R) = \{\, z \in \mathbb{C} : r < |z - z_0| < R \,\},

where rr is the radius of convergence of the principal part (viewed as a power series in 1/(z−z0)1/(z-z_0)) and RR is the radius of convergence of the analytic part. Given a function ff holomorphic on such an annulus, we say ff admits a Laurent expansion on A(z0;r,R)A(z_0;r,R) if there exist coefficients ana_n (n∈Zn \in \mathbb{Z}) such that

f(z)=∑n=−∞∞an(z−z0)nf(z) = \sum_{n=-\infty}^{\infty} a_n (z-z_0)^n

for all z∈A(z0;r,R)z \in A(z_0;r,R), with both one-sided pieces converging.

Intuition

Think of a Taylor series as a lens that only sees the interior of a disk: it is built from data at a single point z0z_0 and it necessarily converges out to the nearest singularity. A Laurent series instead is built for an annular neighborhood, a ring-shaped region that can have singularities both at its center and beyond its outer edge. The negative-power terms are exactly what is needed to describe behavior as z→z0z \to z_0, when z0z_0 itself is a hole in the domain, while the nonnegative-power terms describe ordinary analytic behavior filling out toward the outer radius RR. A very useful mental model: the analytic part ∑n≥0an(z−z0)n\sum_{n\ge 0} a_n(z-z_0)^n converges on the whole disk ∣z−z0∣<R|z-z_0|<R and represents a function holomorphic there; the principal part ∑n≥1a−n(z−z0)−n\sum_{n\ge1} a_{-n}(z-z_0)^{-n}, upon substituting w=1/(z−z0)w = 1/(z-z_0), becomes an ordinary power series ∑n≥1a−nwn\sum_{n\ge 1} a_{-n}w^n in ww, convergent for ∣w∣<1/r|w| < 1/r, i.e. for ∣z−z0∣>r|z-z_0|>r. So the principal part converges on the exterior of a disk of radius rr. A Laurent series therefore converges exactly on the overlap of an "inside" region and an "outside" region — the annulus r<∣z−z0∣<Rr<|z-z_0|<R — and the function it represents is genuinely built from two different analytic pieces glued along that ring.

Figure
The annulus of convergence: inner radius from the principal part, outer radius from the analytic part
z₀ r R annulus r < |z−z₀| < R analytic part: converges |z−z₀|<R principal part: converges |z−z₀|>r

The two-radii picture in coordinates

Concretely, if f(z)=∑n≥0(z−z0)nf(z) = \sum_{n\ge 0}(z-z_0)^n has radius of convergence 11 and g(w)=∑n≥1(13)nwng(w)=\sum_{n\ge1}\left(\tfrac{1}{3}\right)^n w^n has radius of convergence 33, then setting w=1/(z−z0)w=1/(z-z_0) turns gg into a series converging for ∣z−z0∣>1/3|z-z_0|>1/3. The combined Laurent series f+gf+g converges exactly on the annulus 1/3<∣z−z0∣<11/3 < |z-z_0| < 1 — a concrete, checkable instance of the general picture above.

TheoremLaurent's Theorem

Let z0∈Cz_0\in\mathbb C, let 0≤r<R≤∞0\le r<R\le\infty, and let ff be holomorphic on the open annulus

A={z∈C:r<∣z−z0∣<R}.A = \{z\in\mathbb C: r<|z-z_0|<R\}.

Then:

(Existence.) There exist complex numbers ana_n, n∈Zn\in\mathbb Z, such that

f(z)=∑n=−∞∞an(z−z0)nf(z) = \sum_{n=-\infty}^{\infty} a_n(z-z_0)^n

for every z∈Az\in A, with the series converging absolutely and locally uniformly on AA (in the sense of the two one-sided series defined above). Moreover, for every ρ\rho with r<ρ<Rr<\rho<R,

an=12πi∮∣w−z0∣=ρf(w)(w−z0)n+1 dw,n∈Z,a_n = \frac{1}{2\pi i}\oint_{|w-z_0|=\rho} \frac{f(w)}{(w-z_0)^{n+1}}\,dw, \qquad n\in\mathbb Z,

where the circle is traversed once counterclockwise, and this value of ana_n is independent of the choice of ρ∈(r,R)\rho\in(r,R).

(Uniqueness.) The coefficients ana_n are uniquely determined by ff and AA: if f(z)=∑nbn(z−z0)nf(z)=\sum_n b_n(z-z_0)^n is any other Laurent expansion of ff valid throughout AA, then bn=anb_n=a_n for all nn.

Proof

We use the Cauchy Integral Formula for an annulus (obtained from the ordinary Cauchy theory by cutting the annulus with two radial slits), together with the geometric series expansion, exactly as in the derivation of Taylor's theorem for a disk. We give full details.

Step 1: A Cauchy formula on the annulus. Fix z∈Az\in A and choose ρ1,ρ2\rho_1,\rho_2 with r<ρ1<∣z−z0∣<ρ2<Rr<\rho_1<|z-z_0|<\rho_2<R. Let C1C_1 and C2C_2 denote the circles ∣w−z0∣=ρ1|w-z_0|=\rho_1 and ∣w−z0∣=ρ2|w-z_0|=\rho_2, both oriented counterclockwise. We claim

f(z)=12πi∮C2f(w)w−z dw  −  12πi∮C1f(w)w−z dw.(∗)f(z) = \frac{1}{2\pi i}\oint_{C_2}\frac{f(w)}{w-z}\,dw \;-\; \frac{1}{2\pi i}\oint_{C_1}\frac{f(w)}{w-z}\,dw. \tag{$\ast$}

To see this, form the closed region KK between C1C_1 and C2C_2 (the closed annulus ρ1≤∣w−z0∣≤ρ2\rho_1\le|w-z_0|\le\rho_2), which is contained in AA since r<ρ1r<\rho_1 and ρ2<R\rho_2<R, so ff is holomorphic on an open set containing KK. Cut KK along two disjoint radial segments joining C1C_1 to C2C_2 (avoiding the point zz), splitting KK into two simply-connected closed regions K1,K2K_1,K_2, each bounded by a simple closed contour consisting of an arc of C2C_2, a radial segment, an arc of C1C_1, and the other radial segment. The point zz lies in the interior of exactly one of K1,K2K_1,K_2; say z∈int⁡(K1)z\in \operatorname{int}(K_1). Since f(w)/(w−z)f(w)/(w-z) is holomorphic on (a neighborhood of) K2K_2 (which does not contain zz), Cauchy's theorem gives that the integral of f(w)/(w−z)f(w)/(w-z) around ∂K2\partial K_2 vanishes. Since f(w)/(w−z)f(w)/(w-z) is holomorphic on K1K_1 except at the interior point zz, the Cauchy Integral Formula applied to K1K_1 gives

f(z)=12πi∮∂K1f(w)w−z dw.f(z) = \frac{1}{2\pi i}\oint_{\partial K_1} \frac{f(w)}{w-z}\,dw.

Adding the (zero) integral over ∂K2\partial K_2 to this, and observing that the two radial segments are traversed once in each direction across ∂K1\partial K_1 and ∂K2\partial K_2 combined (so their contributions cancel exactly), the sum ∂K1+∂K2\partial K_1 + \partial K_2 reduces to C2C_2 traversed counterclockwise minus C1C_1 traversed counterclockwise. This proves (∗)(\ast).

Step 2: Expand each integral as a power series. Outer integral. For w∈C2w\in C_2 we have ∣w−z0∣=ρ2>∣z−z0∣|w-z_0| = \rho_2 > |z-z_0|, so

1w−z=1(w−z0)−(z−z0)=1w−z0⋅11−z−z0w−z0=∑n=0∞(z−z0)n(w−z0)n+1,\frac{1}{w-z} = \frac{1}{(w-z_0)-(z-z_0)} = \frac{1}{w-z_0}\cdot\frac{1}{1-\frac{z-z_0}{w-z_0}} = \sum_{n=0}^{\infty} \frac{(z-z_0)^n}{(w-z_0)^{n+1}},

the geometric series converging since ∣z−z0w−z0∣=∣z−z0∣ρ2<1\left|\frac{z-z_0}{w-z_0}\right| = \frac{|z-z_0|}{\rho_2} < 1, and in fact uniformly in w∈C2w\in C_2 for zz fixed (the ratio is a fixed number <1<1). Since ff is bounded on the compact set C2C_2, the series

f(w)w−z=∑n=0∞f(w)(z−z0)n(w−z0)−n−1\frac{f(w)}{w-z} = \sum_{n=0}^\infty f(w)(z-z_0)^n (w-z_0)^{-n-1}

converges uniformly in w∈C2w\in C_2, so we may integrate term by term:

12πi∮C2f(w)w−z dw=∑n=0∞an(z−z0)n,an:=12πi∮C2f(w)(w−z0)n+1 dw,n≥0.\frac{1}{2\pi i}\oint_{C_2}\frac{f(w)}{w-z}\,dw = \sum_{n=0}^{\infty} a_n (z-z_0)^n, \qquad a_n := \frac{1}{2\pi i}\oint_{C_2} \frac{f(w)}{(w-z_0)^{n+1}}\,dw, \quad n\ge 0.

Inner integral. For w∈C1w\in C_1 we have ∣w−z0∣=ρ1<∣z−z0∣|w-z_0|=\rho_1 < |z-z_0|, so instead expand around the other variable:

−1w−z=1z−w=1(z−z0)−(w−z0)=1z−z0⋅11−w−z0z−z0=∑m=0∞(w−z0)m(z−z0)m+1,-\frac{1}{w-z} = \frac{1}{z-w} = \frac{1}{(z-z_0)-(w-z_0)} = \frac{1}{z-z_0}\cdot \frac{1}{1-\frac{w-z_0}{z-z_0}} = \sum_{m=0}^\infty \frac{(w-z_0)^m}{(z-z_0)^{m+1}},

convergent since ∣w−z0z−z0∣=ρ1∣z−z0∣<1\left|\frac{w-z_0}{z-z_0}\right| = \frac{\rho_1}{|z-z_0|}<1, again uniformly for w∈C1w \in C_1. As before f(w)f(w) is bounded on C1C_1, so term-by-term integration is valid:

−12πi∮C1f(w)w−z dw=∑m=0∞(12πi∮C1f(w)(w−z0)m dw)(z−z0)−m−1.-\frac{1}{2\pi i}\oint_{C_1}\frac{f(w)}{w-z}\,dw = \sum_{m=0}^{\infty}\left(\frac{1}{2\pi i}\oint_{C_1} f(w)(w-z_0)^m\,dw\right)(z-z_0)^{-m-1}.

Reindexing with n=−m−1n=-m-1 (so nn runs over −1,−2,−3,…-1,-2,-3,\dots as m=0,1,2,…m=0,1,2,\dots), and writing

an:=12πi∮C1f(w)(w−z0)n+1 dw,n≤−1,a_n := \frac{1}{2\pi i}\oint_{C_1} \frac{f(w)}{(w-z_0)^{n+1}}\,dw, \qquad n\le -1,

(noting (w−z0)m=(w−z0)−n−1(w-z_0)^m = (w-z_0)^{-n-1} matches f(w)/(w−z0)n+1f(w)/(w-z_0)^{n+1} with n=−m−1n=-m-1), this sum equals ∑n=−∞−1an(z−z0)n\sum_{n=-\infty}^{-1} a_n(z-z_0)^n. Combining Step 1 and Step 2:

f(z)=∑n=0∞an(z−z0)n+∑n=−∞−1an(z−z0)n=∑n=−∞∞an(z−z0)n,f(z) = \sum_{n=0}^\infty a_n(z-z_0)^n + \sum_{n=-\infty}^{-1} a_n(z-z_0)^n = \sum_{n=-\infty}^{\infty} a_n(z-z_0)^n,

with both one-sided series converging (absolutely, by construction from the convergent geometric series comparison; local uniformity in zz on compact subsets of AA follows by the standard argument bounding the ratio ∣z−z0∣/ρ2|z-z_0|/\rho_2 or ρ1/∣z−z0∣\rho_1/|z-z_0| uniformly away from 11 on compact sets).

Step 3: Independence of ρ\rho and the single-circle formula. The coefficient ana_n for n≥0n\ge 0 was defined using C2={∣w−z0∣=ρ2}C_2 = \{|w-z_0|=\rho_2\}, and for n≤−1n\le -1 using C1={∣w−z0∣=ρ1}C_1=\{|w-z_0|=\rho_1\}; but for any fixed n∈Zn\in\mathbb Z, the function w↦f(w)(w−z0)−n−1w\mapsto f(w)(w-z_0)^{-n-1} is holomorphic throughout the annulus AA (no poles are introduced by the fixed integer power, since w≠z0w\ne z_0 throughout AA). Hence by Cauchy's theorem applied to the annular region between any two circles ∣w−z0∣=ρ|w-z_0|=\rho and ∣w−z0∣=ρ′|w-z_0|=\rho' with r<ρ,ρ′<Rr<\rho,\rho'<R (same slit-and-cancel argument as Step 1, but now with no singularity of the integrand inside), the integral

∮∣w−z0∣=ρf(w)(w−z0)n+1 dw\oint_{|w-z_0|=\rho} \frac{f(w)}{(w-z_0)^{n+1}}\,dw

is the same for every ρ∈(r,R)\rho\in(r,R). This proves both that ana_n can be computed on any circle CρC_\rho in the annulus (not just the particular C1C_1 or C2C_2 used to define it), and gives the stated formula.

Step 4: Uniqueness. Suppose f(z)=∑n=−∞∞bn(z−z0)nf(z) = \sum_{n=-\infty}^\infty b_n(z-z_0)^n for all z∈Az\in A, with both one-sided pieces converging, in particular converging uniformly on every circle ∣z−z0∣=ρ|z - z_0|=\rho, r<ρ<Rr<\rho<R (uniform convergence on compact subsets of the disk of convergence for the analytic part, and of the exterior-disk of convergence for the principal part, are standard facts about ordinary power series, inherited by Laurent series by definition). Fix ρ∈(r,R)\rho\in(r,R) and fix k∈Zk\in\mathbb Z. Divide by (z−z0)k+1(z-z_0)^{k+1} and integrate term-by-term over Cρ={∣z−z0∣=ρ}C_\rho=\{|z-z_0|=\rho\}, justified by the uniform convergence:

∮Cρf(z)(z−z0)k+1 dz=∑n=−∞∞bn∮Cρ(z−z0)n−k−1 dz.\oint_{C_\rho} \frac{f(z)}{(z-z_0)^{k+1}}\,dz = \sum_{n=-\infty}^{\infty} b_n \oint_{C_\rho} (z-z_0)^{n-k-1}\,dz.

Parametrize z−z0=ρeiθz-z_0 = \rho e^{i\theta}, θ∈[0,2π]\theta\in[0,2\pi], so dz=iρeiθdθdz = i\rho e^{i\theta}d\theta and

∮Cρ(z−z0)n−k−1 dz=∫02πρn−k−1ei(n−k−1)θ⋅iρeiθ dθ=iρn−k∫02πei(n−k)θ dθ.\oint_{C_\rho}(z-z_0)^{n-k-1}\,dz = \int_0^{2\pi} \rho^{n-k-1}e^{i(n-k-1)\theta}\cdot i\rho e^{i\theta}\,d\theta = i\rho^{n-k}\int_0^{2\pi} e^{i(n-k)\theta}\,d\theta.

The last integral equals 2π2\pi if n=kn=k and 00 if n≠kn\ne k (integral of a nonzero-frequency complex exponential over a full period vanishes). Hence the right side collapses to 2πi bk2\pi i\, b_k, giving

bk=12πi∮Cρf(z)(z−z0)k+1 dz=ak,b_k = \frac{1}{2\pi i}\oint_{C_\rho} \frac{f(z)}{(z-z_0)^{k+1}}\,dz = a_k,

by Step 3. Since k∈Zk\in\mathbb Z was arbitrary, bn=anb_n=a_n for all nn, proving uniqueness. ■\blacksquare

∎

Worked examples

Example 1: One function, three different annuli. Let

f(z)=1(z−1)(z−3),z0=0.f(z) = \frac{1}{(z-1)(z-3)}, \qquad z_0 = 0.

This function is holomorphic on C∖{1,3}\mathbb{C}\setminus\{1,3\}, so it has three natural regions of interest centered at 00: the disk ∣z∣<1|z|<1, the annulus 1<∣z∣<31<|z|<3, and the exterior ∣z∣>3|z|>3. Partial fractions give

f(z)=12(1z−3−1z−1).f(z) = \frac{1}{2}\left(\frac{1}{z-3} - \frac{1}{z-1}\right).

On ∣z∣<1|z|<1: both ∣z∣<1|z|<1 and ∣z∣<3|z|<3, so expand both terms as geometric series in zz:

−1z−1=11−z=∑n=0∞zn,1z−3=−13⋅11−z/3=−∑n=0∞zn3n+1.\frac{-1}{z-1} = \frac{1}{1-z} = \sum_{n=0}^\infty z^n, \qquad \frac{1}{z-3} = -\frac{1}{3}\cdot\frac{1}{1-z/3} = -\sum_{n=0}^\infty \frac{z^n}{3^{n+1}}.

Hence, valid for ∣z∣<1|z|<1:

f(z)=12∑n=0∞(1−13n+1)zn=13+49z+1327z2+⋯f(z) = \frac{1}{2}\sum_{n=0}^\infty \left(1 - \frac{1}{3^{n+1}}\right) z^n = \frac{1}{3} + \frac{4}{9}z + \frac{13}{27}z^2+\cdots

This is an ordinary Taylor series (pure analytic part, no negative powers), as expected since ff is actually analytic at 00 with nearest singularity at z=1z=1. On 1<∣z∣<31<|z|<3: here ∣z∣>1|z|>1 so we must expand 1/(z−1)1/(z-1) in powers of 1/z1/z, while ∣z∣<3|z|<3 still allows the geometric expansion of 1/(z−3)1/(z-3) in zz:

−1z−1=−1z⋅11−1/z=−∑n=0∞1zn+1,(∣z∣>1)\frac{-1}{z-1} = \frac{-1}{z}\cdot\frac{1}{1-1/z} = -\sum_{n=0}^\infty \frac{1}{z^{n+1}}, \qquad (|z|>1)
1z−3=−13∑n=0∞zn3n,(∣z∣<3).\frac{1}{z-3} = -\frac13\sum_{n=0}^\infty \frac{z^n}{3^n}, \qquad (|z|<3).

So, valid for 1<∣z∣<31<|z|<3:

f(z)=−12∑n=1∞z−n  −  16∑n=0∞zn3n=⋯−12z2−12z−16−z18−⋯f(z) = -\frac{1}{2}\sum_{n=1}^{\infty} z^{-n} \;-\; \frac{1}{6}\sum_{n=0}^{\infty} \frac{z^n}{3^n} = \cdots -\frac{1}{2z^2} - \frac{1}{2z} - \frac{1}{6} - \frac{z}{18} - \cdots

Now there genuinely is a principal part: a−1=−1/2a_{-1} = -1/2, a−2=−1/2a_{-2}=-1/2, a−n=−1/2a_{-n} = -1/2 for all n≥1n\ge1 in fact (check: −12z−n-\frac12 z^{-n} has coefficient −1/2-1/2 for every n≥1n\ge1), and analytic part coefficients an=−16⋅3na_n = -\frac{1}{6\cdot 3^n} for n≥0n\ge 0. On ∣z∣>3|z|>3: now both 1/z1/z and 3/z3/z have modulus <1<1, so expand both terms in powers of 1/z1/z:

−1z−1=−∑n=0∞z−n−1,1z−3=1z⋅11−3/z=∑n=0∞3nz−n−1.\frac{-1}{z-1} = -\sum_{n=0}^\infty z^{-n-1}, \qquad \frac{1}{z-3} = \frac1z\cdot\frac{1}{1-3/z} = \sum_{n=0}^\infty 3^n z^{-n-1}.

So, valid for ∣z∣>3|z|>3:

f(z)=12∑n=0∞(3n−1)z−n−1=0z+1z2+4z3+13z4+⋯f(z) = \frac{1}{2}\sum_{n=0}^\infty (3^n - 1)z^{-n-1} = \frac{0}{z} + \frac{1}{z^2}+\frac{4}{z^3}+\frac{13}{z^4}+\cdots

Here the analytic part is identically zero — every coefficient ana_n with n≥0n\ge 0 vanishes — because f(z)→0f(z)\to 0 like 1/z21/z^2 as z→∞z\to\infty, so there is no room for a nonzero constant, linear, etc. term. This single computation is the clearest possible illustration of why the annulus is part of the data: the same function ff has three genuinely different Laurent expansions centered at the same point z0=0z_0=0, one for each of the three natural annuli determined by the singularities at 11 and 33.

Example 2: An essential singularity, and extracting the residue. Let g(z)=z2e1/zg(z) = z^2 e^{1/z}, holomorphic on C∖{0}\mathbb{C}\setminus\{0\}, and expand around z0=0z_0=0 on the annulus 0<∣z∣<∞0<|z|<\infty. Using the known Taylor series of the exponential, ew=∑k=0∞wk/k!e^{w} = \sum_{k=0}^\infty w^k/k! for all w∈Cw\in\mathbb C, substitute w=1/zw=1/z:

e1/z=∑k=0∞1k! zk=1+1z+12z2+16z3+⋯ ,z≠0,e^{1/z} = \sum_{k=0}^\infty \frac{1}{k!\,z^k} = 1 + \frac1z + \frac{1}{2z^2}+\frac{1}{6z^3}+\cdots, \qquad z\ne 0,

convergent for all z≠0z\ne 0 (i.e. r=0r=0, R=∞R=\infty) since the exponential series converges for every w∈Cw\in\mathbb C. Multiplying by z2z^2 termwise (justified since we are multiplying a convergent series by a fixed polynomial, which preserves convergence and simply shifts/scales each term):

g(z)=z2+z+12+16z+124z2+1120z3+⋯=∑n=−∞2anzn,g(z) = z^2 + z + \frac{1}{2} + \frac{1}{6z} + \frac{1}{24z^2} + \frac{1}{120 z^3}+\cdots = \sum_{n=-\infty}^{2} a_n z^n,

where explicitly a2=1a_2=1, a1=1a_1=1, a0=12a_0=\frac12, and for n≤−1n\le -1, an=1(2−n)!a_n = \dfrac{1}{(2-n)!} (check n=−1n=-1: 13!=16\frac{1}{3!}=\frac16 ✓; n=−2n=-2: 14!=124\frac{1}{4!}=\frac{1}{24} ✓). By Laurent's theorem this expansion is unique on 0<∣z∣<∞0<|z|<\infty, so we may read off the residue directly:

a−1=Res⁡z=0g(z)=13!=16.a_{-1} = \operatorname{Res}_{z=0} g(z) = \frac{1}{3!} = \frac{1}{6}.

No contour integral needed to be computed by hand — the uniqueness clause guarantees that any correct method of producing a two-sided series representing gg on this punctured plane yields the same a−1a_{-1}, which is exactly the value the defining contour integral 12πi∮∣z∣=1g(z) dz\frac{1}{2\pi i}\oint_{|z|=1} g(z)\,dz would produce. This is precisely the mechanism the future Residue Theorem exploits.

CounterexampleWhy the domain must be an honest annulus (holomorphy alone is not enough)

Laurent's theorem requires ff to be holomorphic on a full annulus r<∣z−z0∣<Rr<|z-z_0|<R — a connected ring with no gaps. If the domain is instead only an arc-shaped or slit region, or if ff fails to be holomorphic somewhere inside the claimed annulus, uniqueness and even existence of a global two-sided expansion can fail. Concretely, let f(z)=log⁡zf(z)=\log z (any fixed branch, say the principal branch, holomorphic on C∖(−∞,0]\mathbb{C}\setminus(-\infty,0]), and try to "Laurent expand" around z0=0z_0=0 on the punctured disk 0<∣z∣<10<|z|<1. This set is an annulus (with r=0r=0), so at first glance Laurent's theorem should apply — but ff is not holomorphic on all of it: the principal branch of log⁡z\log z is discontinuous (and not holomorphic) along the negative real axis (−1,0)(-1,0), which lies inside 0<∣z∣<10<|z|<1. The hypothesis "ff holomorphic on the whole open annulus" genuinely fails. Consequently there is no Laurent series ∑anzn\sum a_n z^n converging to log⁡z\log z on all of 0<∣z∣<10<|z|<1. Indeed, if such a series existed, the coefficient formula from the proof would force

an=12πi∮∣z∣=ρlog⁡zzn+1 dza_n = \frac{1}{2\pi i}\oint_{|z|=\rho} \frac{\log z}{z^{n+1}}\,dz

for any fixed radius ρ∈(0,1)\rho\in(0,1) avoiding the branch cut — but the principal branch of log⁡z\log z is not even single-valued/continuous all the way around the full circle ∣z∣=ρ|z|=\rho (it jumps by 2πi2\pi i crossing the negative real axis), so the contour integral around the full circle is not even well-defined for this branch. The failure is not a technicality: log⁡z\log z genuinely has no single-valued holomorphic extension to the full punctured disk (it has a branch point, not an isolated singularity in the Laurent sense), and no two-sided power series in zz can represent a multivalued function. This shows the hypothesis "ff is holomorphic on the entire open annulus" cannot be weakened to "holomorphic on the annulus minus a slit" — the conclusion (a globally convergent two-sided series) genuinely breaks.

Scope, edge cases & common traps

SituationWhat happens
r=0r=0Annulus is a punctured disk 0<∣z−z0∣<R0<|z-z_0|<R; this is the setting for isolated singularities.
R=∞R=\inftyAnnulus is the exterior of a disk (possibly all of C∖{z0}\mathbb{C}\setminus\{z_0\} if also r=0r=0); used to expand "at infinity".
Principal part ≡0\equiv 0Laurent series reduces to the ordinary Taylor series; happens exactly when ff extends holomorphically across ∣z−z0∣=r|z-z_0|=r down to z0z_0 (removable case).
Finitely many negative-power termsCorresponds to a pole; see Singularities and Classification.
Infinitely many nonzero negative-power termsCorresponds to an essential singularity (as in e1/ze^{1/z} above).
Same ff, different annulus around same z0z_0Genuinely different coefficients ana_n — the annulus is part of the data, not a formality (Example 1).
Domain not a full annulus (slit, branch cuts)Theorem does not apply; no global two-sided series need exist (Counterexample).

The mathematician behind it

Origins · 1 minPierre Alphonse Laurent and how this was discovered
PA
Pierre Alphonse Laurent
1813–1854 · French military engineer and mathematician, worked largely outside the Parisian academic mainstream

Laurent trained at the École Polytechnique and worked for years as a military engineer, notably on harbor fortification projects at Le Havre, publishing mathematics on the side. His theorem on two-sided series expansions was submitted to the Paris Academy of Sciences in 1843, essentially completing the theory of power-series representations that Cauchy had begun for the analytic (Taylor) case, by showing what happens when the point of expansion is itself excluded from the domain. Laurent's paper arrived too late for a prize competition deadline judged partly by Cauchy himself, and its full recognition took time, but the result now bears his name universally across complex analysis.

History
- 1821–1825 — Cauchy develops the foundations of complex integration and the Cauchy Integral Theorem/Formula, the machinery Laurent's proof directly builds on.
- 1841 — Karl Weierstrass, independently and roughly contemporaneously, develops closely related power-series ideas for complex functions, though his work on this point circulated later.
- 1843 — Laurent submits his memoir to the Paris Academy establishing the two-sided series expansion on an annulus.
- 1863 — Laurent's memoir is published posthumously (Laurent died in 1854), cementing the name "Laurent series."
- late 1800s onward — The theorem becomes a standard chapter of complex analysis, feeding directly into the classification of isolated singularities and the residue calculus.
Is the Laurent series of a function unique?

Yes, but only relative to a fixed annulus. Laurent's theorem states that if ff is holomorphic on a specific annulus r<∣z−z0∣<Rr<|z-z_0|<R, then the coefficients ana_n representing ff on that annulus are unique, and the proof shows this by extracting each ana_n as a contour integral, a quantity determined entirely by ff and the circle, not by the algebraic method used to find the series. However, if ff is holomorphic on a larger punctured domain containing several natural annuli around z0z_0 (separated by other singularities), the function typically has a different Laurent series on each such annulus, as Example 1 in this article demonstrates concretely.

How is a Laurent series different from a Fourier series?

They are close relatives. If f(z)=∑anznf(z)=\sum a_n z^n converges on an annulus containing the unit circle, restricting to z=eiθz=e^{i\theta} gives f(eiθ)=∑aneinθf(e^{i\theta}) = \sum a_n e^{in\theta}, exactly a complex Fourier series in θ\theta, with the same coefficient formula (a contour integral over ∣z∣=1|z|=1 becomes an integral over θ∈[0,2π]\theta\in[0,2\pi]). The difference is that a Laurent series additionally requires the coefficients to extend to a function holomorphic off the circle in both directions, whereas a general Fourier series only needs L2L^2 or continuity conditions on the boundary circle itself, with no holomorphy requirement inward or outward.

What is the "principal part" used for, concretely?

The principal part ∑n≥1a−n(z−z0)−n\sum_{n\ge1} a_{-n}(z-z_0)^{-n} is exactly the part of ff that is singular at z0z_0; subtracting it from ff leaves a function that extends holomorphically across z0z_0. Counting how many nonzero terms the principal part has (zero, finitely many, or infinitely many) is precisely how isolated singularities are classified as removable, poles, or essential, a classification developed immediately from this expansion. The coefficient a−1a_{-1} specifically, called the residue, further governs the value of contour integrals encircling z0z_0, independent of all the other coefficients.

Can a Laurent series have infinitely many negative-power terms with nonzero coefficients?

Yes — this is exactly what happens at an essential singularity, as with e1/ze^{1/z} expanded around z0=0z_0=0, where every coefficient a−n=1/n!a_{-n}=1/n! for n≥0n\ge0 is nonzero. In this case the principal part is itself an infinite series (in 1/(z−z0)1/(z-z_0)) with infinite radius of convergence in that variable, meaning it converges for every z≠z0z\ne z_0. This is qualitatively different from a pole, where only finitely many negative-power coefficients are nonzero and the principal part is a genuine (finite) rational expression in 1/(z−z0)1/(z-z_0).

Does every convergent two-sided series automatically come from a holomorphic function?

Yes, in the following precise sense: if ∑an(z−z0)n\sum a_n(z-z_0)^n converges (meaning both one-sided pieces converge) throughout some annulus r<∣z−z0∣<Rr<|z-z_0|<R, then it defines a function holomorphic on that entire annulus, because each one-sided piece is an ordinary (Taylor-type) power series representing a holomorphic function on its disk of convergence, and a sum of two holomorphic functions is holomorphic. So Laurent series are exactly the two-sided power series that arise from — and fully characterize — holomorphic functions on annuli; existence and uniqueness together make this a genuine equivalence of data.

What goes wrong if the "annulus" I use isn't simply an annulus, e.g. it has a slit removed?

Then Laurent's theorem simply does not apply, because its hypothesis is holomorphy on the entire open annulus. As the counterexample with log⁡z\log z shows, a function can be perfectly holomorphic on an annulus-minus-a-slit and yet have no single-valued two-sided power series representing it on the full annulus, because going around the annulus picks up a discontinuity (a branch point) that no power series — single-valued by construction — can reproduce. Recognizing branch points versus honest isolated singularities is exactly the first thing to check before attempting a Laurent expansion.

How do I actually compute Laurent coefficients in practice, without doing contour integrals?

In virtually all textbook and research situations, you never evaluate the defining contour integral directly; instead you use known Taylor series (geometric series, exponential, sine/cosine, binomial series) and algebraic manipulation — partial fractions, substitution, term-by-term multiplication — exactly as in the two worked examples here. The uniqueness clause of Laurent's theorem is what licenses this shortcut: however you produce a convergent two-sided series equal to ff on the specified annulus, it must be the Laurent series, matching what the contour integral would have given.

Before this

Where this leads