Holomorphic Functions and Cauchy-Riemann Equations

⚡ Key Takeaways

A function f:U→Cf:U\to\mathbb{C} is holomorphic (complex-differentiable on an open set) if and only if it satisfies the Cauchy–Riemann equations ux=vyu_x = v_y and uy=−vxu_y = -v_x together with differentiability of its real and imaginary parts. Key result in words: complex differentiability is enormously stronger than real differentiability — it forces analyticity (local power-series expansion), infinite differentiability, and harmonicity of the real/imaginary parts. Where it applies: complex analysis, conformal mapping, potential theory, fluid flow, electrostatics, signal processing. Who/when: d'Alembert (1752), Euler, Cauchy (1814–1825), Riemann (1851).

💡 The Core Insight

In real calculus, "differentiable" only constrains behavior along the single real line. In the complex plane, differentiability must hold along infinitely many directions of approach simultaneously, and all those directional limits must agree on one common value f′(z0)f'(z_0). That coherence condition, written out in terms of u=Re⁡fu=\operatorname{Re}f and v=Im⁡fv=\operatorname{Im}f, is the Cauchy–Riemann equations. Geometrically, f′(z0)f'(z_0) acts as multiplication by a single complex number — a rotation-plus-scaling — so holomorphic maps are locally conformal (angle-preserving) wherever f′≠0f'\neq 0.

Formal Statement

DefinitionComplex Differentiability and Holomorphy

Let U⊆CU\subseteq\mathbb{C} be open and f:U→Cf:U\to\mathbb{C}. We say ff is complex-differentiable at z0∈Uz_0\in U if the limit

1
f′(z0)=lim⁡h→0f(z0+h)−f(z0)h,h∈C∖{0}\quad f'(z_0)=\lim_{h\to 0}\frac{f(z_0+h)-f(z_0)}{h},\qquad h\in\mathbb{C}\setminus\{0\}

exists in C\mathbb{C}. The function ff is holomorphic on UU if it is complex-differentiable at every point of UU. A function holomorphic on all of C\mathbb{C} is called entire. Symbols: z0=x0+iy0z_0=x_0+iy_0 is the base point; h=Δx+iΔyh=\Delta x+i\Delta y is a complex increment; f′(z0)∈Cf'(z_0)\in\mathbb{C} is the complex derivative; UU is an open domain (every point has a neighborhood inside UU).

Derivation hint: This is the direct transplant of the difference-quotient definition into C\mathbb{C}; the crucial novelty is that division by the complex number hh requires hh to approach 00 from all planar directions, not just ±\pm along a line.

DefinitionThe Cauchy–Riemann Equations

Write f(x+iy)=u(x,y)+i v(x,y)f(x+iy)=u(x,y)+i\,v(x,y) with u,v:U→Ru,v:U\to\mathbb{R} real-valued. The Cauchy–Riemann (CR) equations are

2
∂u∂x=∂v∂y,∂u∂y=−∂v∂x\quad \frac{\partial u}{\partial x}=\frac{\partial v}{\partial y},\qquad \frac{\partial u}{\partial y}=-\frac{\partial v}{\partial x}

equivalently written ux=vyu_x=v_y and uy=−vxu_y=-v_x. Symbols: ux=∂u/∂xu_x=\partial u/\partial x etc. denote first-order partial derivatives. These two real PDEs couple uu and vv.

Derivation hint: Compute the limit [1] first along the real direction h=Δx→0h=\Delta x\to 0 and then along the imaginary direction h=i Δy→0h=i\,\Delta y\to 0, and equate the two expressions for f′(z0)f'(z_0).

TheoremCauchy–Riemann Characterization of Holomorphy

Let f=u+ivf=u+iv be defined on an open set UU. Then:

(a) If ff is complex-differentiable at z0z_0, then the partials ux,uy,vx,vyu_x,u_y,v_x,v_y exist at z0z_0 and satisfy the CR equations [2] there, with

3
f′(z0)=ux+i vx=vy−i uy.\quad f'(z_0)=u_x+i\,v_x=v_y-i\,u_y.

(b) Conversely, if u,vu,v are real-differentiable at z0z_0 (in particular if ux,uy,vx,vyu_x,u_y,v_x,v_y exist and are continuous near z0z_0) and satisfy the CR equations at z0z_0, then ff is complex-differentiable at z0z_0.

Derivation hint: (a) follows from directional limits as in [2]; (b) uses the first-order Taylor (total-differential) expansion of uu and vv and substitutes the CR relations to factor the increment as f′(z0) h+o(∣h∣)f'(z_0)\,h + o(|h|).

TheoremHolomorphic ⟹ Harmonic

If f=u+ivf=u+iv is holomorphic on UU and (as is automatic) C∞C^\infty, then uu and vv are harmonic:

4
Δu=uxx+uyy=0,Δv=vxx+vyy=0,\quad \Delta u = u_{xx}+u_{yy}=0,\qquad \Delta v = v_{xx}+v_{yy}=0,

and vv is the harmonic conjugate of uu.

Derivation hint: Differentiate ux=vyu_x=v_y by xx and uy=−vxu_y=-v_x by yy, then add; equality of mixed partials vxy=vyxv_{xy}=v_{yx} cancels the cross terms.

Why It Holds / Derivation

ProofDeriving the CR Equations from Complex Differentiability

Assume f′(z0)f'(z_0) exists. We exploit that the limit [1] is independent of the path of approach.

Step 1 — Approach along the real axis. Let h=Δx∈Rh=\Delta x\in\mathbb{R}, Δx→0\Delta x\to 0. Then

f(z0+Δx)−f(z0)Δx=u(x0+Δx,y0)−u(x0,y0)Δx+i v(x0+Δx,y0)−v(x0,y0)Δx.\frac{f(z_0+\Delta x)-f(z_0)}{\Delta x}=\frac{u(x_0+\Delta x,y_0)-u(x_0,y_0)}{\Delta x}+i\,\frac{v(x_0+\Delta x,y_0)-v(x_0,y_0)}{\Delta x}.

Taking Δx→0\Delta x\to 0 gives

5
f′(z0)=ux(z0)+i vx(z0).\quad f'(z_0)=u_x(z_0)+i\,v_x(z_0).

Derivation hint: A real increment leaves the imaginary coordinate fixed, isolating the xx-partials of uu and vv.

Step 2 — Approach along the imaginary axis. Let h=i Δyh=i\,\Delta y, Δy∈R\Delta y\in\mathbb{R}, Δy→0\Delta y\to 0. Since 1/(i Δy)=−i/Δy1/(i\,\Delta y)=-i/\Delta y,

f(z0+iΔy)−f(z0)i Δy=−i[u(x0,y0+Δy)−u(x0,y0)Δy+i v(x0,y0+Δy)−v(x0,y0)Δy].\frac{f(z_0+i\Delta y)-f(z_0)}{i\,\Delta y}=-i\left[\frac{u(x_0,y_0+\Delta y)-u(x_0,y_0)}{\Delta y}+i\,\frac{v(x_0,y_0+\Delta y)-v(x_0,y_0)}{\Delta y}\right].

Taking Δy→0\Delta y\to 0,

6
f′(z0)=−i uy(z0)+vy(z0)=vy(z0)−i uy(z0).\quad f'(z_0)=-i\,u_y(z_0)+v_y(z_0)=v_y(z_0)-i\,u_y(z_0).

Derivation hint: The factor −i-i from 1/i1/i swaps and signs the real/imaginary contributions of the yy-partials.

Step 3 — Equate real and imaginary parts. Because the limit is unique, [5] and [6] must agree:

ux+i vx=vy−i uy.u_x+i\,v_x = v_y - i\,u_y.

Matching real parts: ux=vyu_x=v_y. Matching imaginary parts: vx=−uyv_x=-u_y, i.e. uy=−vxu_y=-v_x. These are exactly the CR equations [2]. ■\blacksquare

∎
ProofConverse: CR + Real-Differentiability ⟹ Complex-Differentiability

Assume u,vu,v are real-differentiable at z0z_0 with continuous partials satisfying [2]. By the total differential,

u(z0+h)−u(z0)=ux Δx+uy Δy+ε1(h),v(z0+h)−v(z0)=vx Δx+vy Δy+ε2(h),u(z_0+h)-u(z_0)=u_x\,\Delta x+u_y\,\Delta y+\varepsilon_1(h),\quad v(z_0+h)-v(z_0)=v_x\,\Delta x+v_y\,\Delta y+\varepsilon_2(h),

where εj(h)=o(∣h∣)\varepsilon_j(h)=o(|h|). Then, with h=Δx+iΔyh=\Delta x+i\Delta y:

Step 1 — Combine.

f(z0+h)−f(z0)=(uxΔx+uyΔy)+i(vxΔx+vyΔy)+o(∣h∣).f(z_0+h)-f(z_0)=(u_x\Delta x+u_y\Delta y)+i(v_x\Delta x+v_y\Delta y)+o(|h|).

Step 2 — Substitute CR (uy=−vx, vy=uxu_y=-v_x,\ v_y=u_x).

=(uxΔx−vxΔy)+i(vxΔx+uxΔy)+o(∣h∣)=(ux+ivx)(Δx+iΔy)+o(∣h∣).=(u_x\Delta x - v_x\Delta y)+i(v_x\Delta x+u_x\Delta y)+o(|h|)=(u_x+i v_x)(\Delta x+i\Delta y)+o(|h|).

Step 3 — Factor. The right side is (ux+ivx) h+o(∣h∣)(u_x+i v_x)\,h+o(|h|), so

7
f(z0+h)−f(z0)h=(ux+ivx)+o(∣h∣)h→h→0ux+ivx.\quad \frac{f(z_0+h)-f(z_0)}{h}=(u_x+i v_x)+\frac{o(|h|)}{h}\xrightarrow[h\to0]{}u_x+i v_x.

Thus f′(z0)=ux+ivxf'(z_0)=u_x+iv_x exists. ■\blacksquare Derivation hint: The whole trick is that CR is precisely the algebraic condition under which the real-linear map (Δx,Δy)↦(du,dv)(\Delta x,\Delta y)\mapsto(\text{du},\text{dv}) becomes complex-linear, i.e. equals multiplication by a single complex number.

∎
⚠️ Real-differentiability is needed, not just existence of partials

The mere existence of the four partial derivatives satisfying CR at a single point is not enough for complex differentiability; one needs full real-differentiability (e.g. continuous partials). There exist functions with all partials existing and CR holding at the origin that are not complex-differentiable there. Always verify differentiability of u,vu,v, not just that the partials exist.

mathematician · Bernhard Riemann | 1826–1866

Though the equations were known to d'Alembert (1752) in fluid dynamics and developed computationally by Euler and Cauchy, it was Bernhard Riemann in his revolutionary 1851 doctoral dissertation who placed them at the geometric heart of function theory. Riemann defined a holomorphic function essentially by the CR system, emphasizing conformality, harmonicity, and the deep link to potential theory, and introduced the notion of a Riemann surface to handle multivaluedness. His geometric, almost physical viewpoint — functions as conformal flows — complemented Cauchy's analytic, integral-based approach and Weierstrass's power-series approach, and together these three perspectives constitute modern complex analysis.

Building Intuition

🧠 Intuition

Think of ff as a map deforming the plane. Its real Jacobian at z0z_0 is

J=(uxuyvxvy).J=\begin{pmatrix}u_x & u_y\\ v_x & v_y\end{pmatrix}.

A general real-differentiable map can stretch differently in different directions and shear. But the CR equations force ux=vyu_x=v_y and uy=−vxu_y=-v_x, so

J=(a−bba),a=ux, b=vx,J=\begin{pmatrix}a & -b\\ b & a\end{pmatrix},\qquad a=u_x,\ b=v_x,

which is exactly the matrix of multiplication by the complex number a+iba+ib — a uniform rotation by arg⁡(a+ib)\arg(a+ib) and scaling by ∣a+ib∣=∣f′(z0)∣|a+ib|=|f'(z_0)|. That's why holomorphic maps preserve angles (are conformal) wherever f′≠0f'\neq 0: there's no shear, no anisotropic stretch, just spin-and-scale. The "infinitely many directions agree" requirement is geometrically the statement "the local map is a similarity transformation."

Figure
Conformality: a holomorphic map sends the orthogonal grid through z0z_0 to an orthogonal grid through f(z0)f(z_0), rotated by arg⁡f′\arg f' and scaled by ∣f′∣|f'|.
z₀ z-plane (domain) f, rotate+scale f(z₀) w-plane (image): angles preserved

Worked Examples

ExampleExample 1 — f(z)=z2f(z)=z^2 is entire; verify CR and find f′f'

Given: f(z)=z2f(z)=z^2. Write z=x+iyz=x+iy.

Step 1 — Separate into u,vu,v.

z2=(x+iy)2=x2−y2+2ixy ⇒ u=x2−y2,v=2xy.z^2=(x+iy)^2=x^2-y^2+2ixy\ \Rightarrow\ u=x^2-y^2,\quad v=2xy.

Step 2 — Compute partials.

ux=2x,uy=−2y,vx=2y,vy=2x.u_x=2x,\quad u_y=-2y,\quad v_x=2y,\quad v_y=2x.

Step 3 — Check CR equations [2].

ux=2x=vy ✓,uy=−2y=−vx=−(2y) ✓.u_x=2x=v_y\ \checkmark,\qquad u_y=-2y=-v_x=-(2y)\ \checkmark.

Partials are polynomials, hence continuous everywhere, so ff is holomorphic on all of C\mathbb{C} (entire).

Step 4 — Derivative via [3].

f′(z)=ux+ivx=2x+i 2y=2(x+iy)=2z.f'(z)=u_x+iv_x=2x+i\,2y=2(x+iy)=2z.
 f(z)=z2 is entire, and f′(z)=2z \boxed{\,f(z)=z^2\text{ is entire, and } f'(z)=2z\,}
ExampleExample 2 — f(z)=zˉf(z)=\bar z is nowhere holomorphic

Given: f(z)=zˉ=x−iyf(z)=\bar z=x-iy.

Step 1 — Identify u,vu,v. u=x, v=−y.u=x,\ v=-y.

Step 2 — Partials.

ux=1,uy=0,vx=0,vy=−1.u_x=1,\quad u_y=0,\quad v_x=0,\quad v_y=-1.

Step 3 — Test CR equations [2].

ux=1,vy=−1 ⇒ ux≠vy.u_x=1,\quad v_y=-1\ \Rightarrow\ u_x\neq v_y.

CR fails at every point, since 1≠−11\neq-1 identically.

Step 4 — Confirm via the limit [1]. Along h=Δxh=\Delta x: quotient =Δx‾/Δx=1=\overline{\Delta x}/\Delta x=1. Along h=iΔyh=i\Delta y: quotient =iΔy‾/(iΔy)=(−iΔy)/(iΔy)=−1=\overline{i\Delta y}/(i\Delta y)=(-i\Delta y)/(i\Delta y)=-1. The two directional limits (11 and −1-1) disagree, so f′f' exists nowhere.

 f(z)=zˉ is real-smooth but nowhere complex-differentiable \boxed{\,f(z)=\bar z\text{ is real-smooth but nowhere complex-differentiable}\,}

Common Misconceptions

⚠️ CR equations alone are insufficient

Satisfying the CR equations at a point does not guarantee complex differentiability there. You additionally need u,vu,v to be (real-)differentiable, which is guaranteed if the partials exist and are continuous in a neighborhood. CR at an isolated point, without surrounding regularity, can hold while f′f' fails to exist.

CounterexampleCR holds at a point yet ff is not differentiable there

Define

f(z)={zˉ 2z,z≠0,0,z=0.f(z)=\begin{cases}\dfrac{\bar z^{\,2}}{z}, & z\neq 0,\\[4pt] 0,& z=0.\end{cases}

At the origin the difference quotients along the axes both give 00, so the partials of uu and vv exist at 00 and one can verify the CR equations hold at the origin. Yet approaching along z=reiθz=re^{i\theta} gives

f(z)−f(0)z−0=zˉ2/zz=zˉ2z2=e−4iθ,\frac{f(z)-f(0)}{z-0}=\frac{\bar z^{2}/z}{z}=\frac{\bar z^{2}}{z^{2}}=e^{-4i\theta},

which depends on θ\theta. Hence the limit [1] does not exist and ff is not complex-differentiable at 00, despite CR holding there. The missing ingredient is real-differentiability of u,vu,v at the origin.

⚠️ Holomorphic is not the same as "smooth" in the real sense

A function can be infinitely real-differentiable (e.g. f(z)=∣z∣2=zzˉf(z)=|z|^2=z\bar z) yet fail to be holomorphic anywhere except possibly isolated points. ∣z∣2|z|^2 has u=x2+y2, v=0u=x^2+y^2,\ v=0, so ux=2x, vy=0u_x=2x,\ v_y=0; CR forces x=y=0x=y=0, meaning ff is complex-differentiable only at the origin — and even there it is not holomorphic, since holomorphy requires differentiability on an open set.

Connections

connection

The CR equations are the gateway to the rest of complex analysis. Once holomorphy is established, Cauchy's Integral Theorem shows contour integrals of ff over closed curves vanish, and Cauchy's Integral Formula reproduces ff inside a contour from its boundary values, instantly yielding analyticity via Power Series and Analyticity. The harmonicity result [4] ties the subject to Harmonic Functions and Potential Theory and Laplace's equation. Conformality links holomorphy to Conformal Mapping and the Riemann Mapping Theorem. The "∂/∂zˉ=0\partial/\partial\bar z=0" reformulation connects to Wirtinger Derivatives, and the rigidity of holomorphic functions underlies Liouville's Theorem and the Maximum Modulus Principle.

💡 The Core Insight

Introducing the operators

8
∂∂z=12 ⁣(∂∂x−i∂∂y),∂∂zˉ=12 ⁣(∂∂x+i∂∂y),\quad \frac{\partial}{\partial z}=\tfrac12\!\left(\frac{\partial}{\partial x}-i\frac{\partial}{\partial y}\right),\qquad \frac{\partial}{\partial\bar z}=\tfrac12\!\left(\frac{\partial}{\partial x}+i\frac{\partial}{\partial y}\right),

the entire CR system collapses to the single elegant equation

9
∂f∂zˉ=0.\quad \frac{\partial f}{\partial\bar z}=0.

Derivation hint for [8]–[9]: Apply ∂/∂zˉ\partial/\partial\bar z to f=u+ivf=u+iv, expand, and observe that the real and imaginary parts of the result vanish exactly when ux=vyu_x=v_y and uy=−vxu_y=-v_x. Holomorphic functions are precisely those "independent of zˉ\bar z."

Real-World Applications

🔬 Where this appears

Historical Development

⏳ History
- 1752 — Jean le Rond d'Alembert writes down what are essentially the CR equations while studying two-dimensional fluid flow, relating velocity components.
- 1777 — Leonhard Euler uses the same relations in evaluating real integrals via complex substitution, developing computational facility with u,vu,v pairs.
- 1814–1825 — Augustin-Louis Cauchy systematically develops the equations within his integral-based theory of complex functions, proving the foundational integral theorem; the equations appear as a necessary condition for path-independence of integrals.
- 1841 — Cauchy formalizes the role of the equations and contour integration in his published memoirs, establishing the residue calculus.
- 1851 — Bernhard Riemann, in his doctoral thesis Grundlagen für eine allgemeine Theorie der Functionen einer veränderlichen complexen Größe, makes the equations the geometric definition of a holomorphic function, links them to conformality and Riemann surfaces.
- 1857 — Riemann applies the theory to Abelian functions and develops potential-theoretic methods (the Dirichlet principle).
- 1860s–1880s — Karl Weierstrass builds an alternative, rigorous foundation via convergent power series, proving the equivalence of differentiability and analyticity and clarifying logical gaps in the Dirichlet principle.
- 1900–1902 — Édouard Goursat proves that mere complex differentiability (without assuming continuity of f′f') suffices for Cauchy's theorem (the Cauchy–Goursat theorem), removing earlier regularity hypotheses.
- 20th century — The Wirtinger calculus and the operator form ∂f/∂zˉ=0\partial f/\partial\bar z=0 unify the equations with the theory of several complex variables and PDE.

Subtleties & Limits

open-problem

Several deep refinements and boundaries surround the basic theory. The Looman–Menchoff theorem (1923–1936) shows that if ff is continuous on UU and its partials exist and satisfy CR everywhere — without assuming continuity of the partials — then ff is holomorphic; this sharpens exactly how weak the hypotheses can be, and removing continuity of ff breaks the conclusion. In several complex variables, the system ∂f/∂zˉj=0\partial f/\partial\bar z_j=0 becomes an overdetermined PDE (the ∂ˉ\bar\partial-problem), whose solvability and regularity (Hörmander's L2L^2 estimates, the Cauchy–Riemann complex on CR manifolds) remain an active research area with subtle obstructions absent in one variable. The interplay between holomorphy, harmonicity, and the inhomogeneous equation ∂ˉu=f\bar\partial u=f underlies modern complex geometry, and questions about boundary behavior, the Dirichlet problem on rough domains, and quasiconformal generalizations (where CR is relaxed to a Beltrami equation) continue to drive analysis.

NF
Note for the applied reader

For physicists: whenever you see a 2D harmonic field with an associated "conjugate" field (potential/stream function, equipotential/field line), you are implicitly using a holomorphic function, and the CR equations are the statement that those two fields are orthogonal and co-determined. The complex derivative f′f' gives the local field strength and orientation in a single complex number — a major computational economy.

See Also

Cauchy's Integral Formula · Conformal Mapping · Harmonic Functions and Potential Theory · Wirtinger Derivatives · Liouville's Theorem · Maximum Modulus Principle · Riemann Mapping Theorem · Power Series and Analyticity