# Cauchy–Riemann equations

(Redirected from Cauchy-Riemann equations)

In the field of complex analysis in mathematics, the Cauchy–Riemann equations, named after Augustin Cauchy and Bernhard Riemann, consist of a system of two partial differential equations which, together with certain continuity and differentiability criteria, form a necessary and sufficient condition for a complex function to be complex differentiable, that is holomorphic. This system of equations first appeared in the work of Jean le Rond d'Alembert (d'Alembert 1752). Later, Leonhard Euler connected this system to the analytic functions (Euler 1797). Cauchy (1814) then used these equations to construct his theory of functions. Riemann's dissertation (Riemann 1851) on the theory of functions appeared in 1851.

The Cauchy–Riemann equations on a pair of real-valued functions of two real variables u(x,y) and v(x,y) are the two equations:

 (1a) $\dfrac{ \partial u }{ \partial x } = \dfrac{ \partial v }{ \partial y }$ (1b) $\dfrac{ \partial u }{ \partial y } = -\dfrac{ \partial v }{ \partial x }$

Typically u and v are taken to be the real and imaginary parts respectively of a complex-valued function of a single complex variable z = x + iy, f(x + iy) = u(x,y) + iv(x,y). Suppose that u and v are real-differentiable at a point in an open subset of C (C is the set of complex numbers), which can be considered as functions from R2 to R. This implies that the partial derivatives of u and v exist (although they need not be continuous) and we can approximate small variations of f linearly. Then f = u + iv is complex-differentiable at that point if and only if the partial derivatives of u and v satisfy the Cauchy–Riemann equations (1a) and (1b) at that point. The sole existence of partial derivatives satisfying the Cauchy–Riemann equations is not enough to ensure complex differentiability at that point. It is necessary that u and v be real differentiable, which is a stronger condition than the existence of the partial derivatives, but it is not necessary that these partial derivatives be continuous.

Holomorphy is the property of a complex function of being differentiable at every point of an open and connected subset of C (this is called a domain in C). Consequently, we can assert that a complex function f, whose real and imaginary parts u and v are real-differentiable functions, is holomorphic if and only if, equations (1a) and (1b) are satisfied throughout the domain we are dealing with.

The reason why Euler and some other authors relate the Cauchy–Riemann equations with analyticity is that a major theorem in complex analysis says that holomorphic functions are analytic and vice versa. This means that, in complex analysis, a function that is complex-differentiable in a whole domain (holomorphic) is the same as an analytic function. This is not true for real differentiable functions.

## Interpretation and reformulation

The equations are one way of looking at the condition on a function to be differentiable in the sense of complex analysis: in other words they encapsulate the notion of function of a complex variable by means of conventional differential calculus. In the theory there are several other major ways of looking at this notion, and the translation of the condition into other language is often needed.

### Conformal mappings

First, the Cauchy–Riemann equations may be written in complex form

(2)    ${ i \dfrac{ \partial f }{ \partial x } } = \dfrac{ \partial f }{ \partial y } .$

In this form, the equations correspond structurally to the condition that the Jacobian matrix is of the form

$\begin{pmatrix} a & -b \\ b & \;\; a \end{pmatrix},$

where $\scriptstyle a=\partial u/\partial x=\partial v/\partial y$ and $\scriptstyle b=\partial v/\partial x=-\partial u/\partial y$. A matrix of this form is the matrix representation of a complex number. Geometrically, such a matrix is always the composition of a rotation with a scaling, and in particular preserves angles. The Jacobian of a function f(z) takes infinitesimal line segments at the intersection of two curves in z and rotates them to the corresponding segments in f(z). Consequently, a function satisfying the Cauchy–Riemann equations, with a nonzero derivative, preserves the angle between curves in the plane. That is, the Cauchy–Riemann equations are the conditions for a function to be conformal.

### Complex differentiability

Suppose that

$f(z) = u(z) + i \cdot v(z)$

is a function of a complex number z. Then the complex derivative of f at a point z0 is defined by

$\lim_{\underset{h\in\mathbf{C}}{h\to 0}} \frac{f(z_0+h)-f(z_0)}{h} = f'(z_0)$

provided this limit exists.

If this limit exists, then it may be computed by taking the limit as h → 0 along the real axis or imaginary axis; in either case it should give the same result. Approaching along the real axis, one finds

$\lim_{\underset{h\in\mathbf{R}}{h\to 0}} \frac{f(z_0+h)-f(z_0)}{h} = \frac{\partial f}{\partial x}(z_0).$

On the other hand, approaching along the imaginary axis,

$\lim_{\underset{h\in \mathbf{R}}{h\to 0}} \frac{f(z_0+ih)-f(z_0)}{ih} =\frac{1}{i}\frac{\partial f}{\partial y}(z_0).$

The equality of the derivative of f taken along the two axes is

$i\frac{\partial f}{\partial x}(z_0)=\frac{\partial f}{\partial y}(z_0),$

which are the Cauchy–Riemann equations (2) at the point z0.

Conversely, if f : C → C is a function which is differentiable when regarded as a function on R2, then f is complex differentiable if and only if the Cauchy–Riemann equations hold. In other words, if u and v are real-differentiable functions of two real variables, obviously u + iv is a (complex-valued) real-differentiable function, but u + iv is complex-differentiable if and only if the Cauchy–Riemann equations hold.

Indeed, following Rudin (1966), suppose f is a complex function defined in an open set Ω ⊂ C. Then, writing z = x + iy for every z ∈ Ω, one can also regard Ω as an open subset of R2, and f as a function of two real variables x and y, which maps Ω ⊂ R2 to C. We consider the Cauchy–Riemann equations at z = z0. So assume f is differentiable at z0, as a function of two real variables from Ω to C. This is equivalent to the existence of the following linear approximation

$f(z_0 + \Delta z) - f(z_0) = f_x \Delta x + f_y \Delta y + \eta(\Delta z) \Delta z \,$

where z = x + iy and η(Δz) → 0 as Δz0. Since $\Delta z + \Delta \bar{z}= 2 \Delta x$ and $\Delta z - \Delta \bar{z}=2i \Delta y$, the above can be re-written as

$\Delta f(z_0) = \frac{f_x - if_y}{2}\Delta z + \frac{f_x + if_y}{2}\Delta \bar{z} + \eta(\Delta z) \Delta z\,$

Defining the two Wirtinger derivatives as

$\frac{\partial}{\partial z} = \frac{1}{2} \Bigl( \frac{\partial}{\partial x} - i \frac{\partial}{\partial y} \Bigr), \;\;\; \frac{\partial}{\partial\bar{z}}= \frac{1}{2} \Bigl( \frac{\partial}{\partial x} + i \frac{\partial}{\partial y} \Bigr),$

in the limit $\Delta z \rightarrow 0,\Delta \bar{z} \rightarrow 0$ the above equality can be written as

$\left.\frac{df}{dz}\right|_{z=z_0} = \left.\frac{\partial f}{\partial z}\right |_{z=z_0} + \left.\frac{\partial f}{\partial\bar{z}}\right |_{z=z_0} \cdot \frac{\bar{dz}}{dz} + \eta(\Delta z), \;\;\;\;(\Delta z \neq 0).$

For real values of z, we have $\bar{dz}/dz = 1$ and for purely imaginary z we have $\bar{dz}/dz = -1$. Similarly, when approaching z0 from different directions in the complex plane, the value of $\bar{dz}/dz$ is different. But since for complex differentiability the derivative should be the same, approaching from any direction, hence f is complex differentiable at z0 if and only if $(\partial f/\partial\bar{z}) = 0$ at $z=z_0$. But this is exactly the Cauchy–Riemann equations, thus f is differentiable at z0 if and only if the Cauchy–Riemann equations hold at z0.

### Independence of the complex conjugate

The above proof suggests another interpretation of the Cauchy–Riemann equations. The complex conjugate of z, denoted $\bar{z}$, is defined by

$\overline{x + iy} := x - iy$

for real x and y. The Cauchy–Riemann equations can then be written as a single equation

(3)    $\dfrac{\partial f}{\partial\bar{z}} = 0$

by using the Wirtinger derivative with respect to the conjugate variable. In this form, the Cauchy–Riemann equations can be interpreted as the statement that f is independent of the variable $\bar{z}$. As such, we can view analytic functions as true functions of one complex variable as opposed to complex functions of two real variables.

### Physical interpretation

Contour plot of a pair u and v satisfying the Cauchy-Riemann equations. Streamlines (v=const, red) are perpendicular to equipotentials (u=const, blue). The point (0,0) is a stationary point of the potential flow, with six streamlines meeting, and six equipotentials also meeting and bisecting the angles formed by the streamlines.

A standard physical interpretation of the Cauchy–Riemann equations going back to Riemann's work on function theory (see Klein 1893) is that u represents a velocity potential of an incompressible steady fluid flow in the plane, and v is its stream function. Suppose that the pair of (twice continuously differentiable) functions $u,v$ satisfies the Cauchy–Riemann equations. We will take u to be a velocity potential, meaning that we imagine a flow of fluid in the plane such that the velocity vector of the fluid at each point of the plane is equal to the gradient of u, defined by

$\nabla u = \frac{\partial u}{\partial x}\mathbf i + \frac{\partial u}{\partial y}\mathbf j$

By differentiating the Cauchy–Riemann equations a second time, one shows that u solves Laplace's equation:

$\frac{\partial^2u}{\partial x^2} + \frac{\partial^2u}{\partial y^2} = 0.$

That is, u is a harmonic function. This means that the divergence of the gradient is zero, and so the fluid is incompressible.

The function v also satisfies the Laplace equation, by a similar analysis. Also, the Cauchy–Riemann equations imply that the dot product $\nabla u\cdot\nabla v = 0$. This implies that the gradient of u must point along the the $v=const$ curves; so these are the streamlines of the flow. The $u=const$ curves are the equipotential curves of the flow.

A holomorphic function can therefore be visualized by plotting the two families of level curves $u=const$ and $v=const$. Near points where the gradient of u (or, equivalently, v) is not zero, these families form an orthogonal family of curves. At the points where $\nabla u=0$, the stationary points of the flow, the equipotential curves of $u=const$ intersect. The streamlines also intersect at the same point, bisecting the angles formed by the equipotential curves.

### Harmonic vector field

Another interpretation of the Cauchy–Riemann equations can be found in Pólya & Szegő (1978). Suppose that u and v satisfy the Cauchy–Riemann equations in an open subset of R2, and consider the vector field

$\bar{f} = \begin{bmatrix}u\\ -v\end{bmatrix}$

regarded as a (real) two-component vector. Then the second Cauchy–Riemann equation (1b) asserts that $\bar{f}$ is irrotational (its curl is 0):

$\frac{\partial (-v)}{\partial x} - \frac{\partial u}{\partial y} = 0.$

The first Cauchy–Riemann equation (1a) asserts that the vector field is solenoidal (or divergence-free):

$\frac{\partial u}{\partial x} + \frac{\partial (-v)}{\partial y}=0.$

Owing respectively to Green's theorem and the divergence theorem, such a field is necessarily a conservative one, and it is free from sources or sinks, having net flux equal to zero through any open domain without holes. (These two observations combine as real and imaginary parts in Cauchy's integral theorem.) In fluid dynamics, such a vector field is a potential flow (Chanson 2007). In magnetostatics, such vector fields model static magnetic fields on a region of the plane containing no current. In electrostatics, they model static electric fields in a region of the plane containing no electric charge.

This interpretation can equivalently be restated in the language of differential forms. The pair u,v satisfy the Cauchy–Riemann equations if and only if the one-form $v\,dx + u\, dy$ is both closed and coclosed (a harmonic differential form).

### Preservation of complex structure

Another formulation of the Cauchy–Riemann equations involves the complex structure in the plane, given by

$J=\begin{bmatrix}0&1\\-1&0\end{bmatrix}.$

This is a complex structure in the sense that the square of J is the negative of the 2×2 identity matrix: $J^2=-I$. As above, if u(x,y),v(x,y) are two functions in the plane, put

$f(x,y) = \begin{bmatrix}u(x,y)\\v(x,y)\end{bmatrix}.$

The Jacobian matrix of f is the matrix of partial derivatives

$Df(x,y) = \begin{bmatrix}\frac{\partial u}{\partial x} & \frac{\partial u}{\partial y}\\\frac{\partial v}{\partial x} & \frac{\partial v}{\partial y}\end{bmatrix}$

Then the pair of functions u, v satisfies the Cauchy–Riemann equations if and only if the 2×2 matrix Df commutes with J (Kobayashi & Nomizu 1969, Proposition IX.2.2)

### Other representations

Other representations of the Cauchy–Riemann equations occasionally arise in other coordinate systems. If (1a) and (1b) hold for a differentiable pair of functions u and v, then so do

$\frac{\partial u}{\partial n} = \frac{\partial v}{\partial s},\quad \frac{\partial v}{\partial n} = -\frac{\partial u}{\partial s}$

for any coordinate system (n(x, y), s(x, y)) such that the pair (∇n, ∇s) is orthonormal and positively oriented. As a consequence, in particular, in the system of coordinates given by the polar representation z = re, the equations then take the form

${ \partial u \over \partial r } = {1 \over r}{ \partial v \over \partial \theta},\quad{ \partial v \over \partial r } = -{1 \over r}{ \partial u \over \partial \theta}.$

Combining these into one equation for f gives

${\partial f \over \partial r} = {1 \over i r}{\partial f \over \partial \theta}.$

The inhomogeneous Cauchy–Riemann equations consist of the two equations for a pair of unknown functions u(x,y) and v(x,y) of two real variables

$\frac{\partial u}{\partial x}-\frac{\partial v}{\partial y} = \alpha(x,y)$
$\frac{\partial u}{\partial y}+\frac{\partial v}{\partial x} = \beta(x,y)$

for some given functions α(x,y) and β(x,y) defined in an open subset of R2. These equations are usually combined into a single equation

$\frac{\partial f}{\partial\bar{z}} = \varphi(z,\bar{z})$

where f = u + iv and φ = (α + iβ)/2.

If φ is Ck, then the inhomogeneous equation is explicitly solvable in any bounded domain D, provided φ is continuous on the closure of D. Indeed, by the Cauchy integral formula,

$f(\zeta,\bar{\zeta}) = \frac{1}{2\pi i}\iint_D \varphi(z,\bar{z})\frac{dz\wedge d\bar{z}}{z-\zeta}$

for all ζ ∈ D.

## Generalizations

### Goursat's theorem and its generalizations

Suppose that f = u + iv is a complex-valued function which is differentiable as a function f : R2R2. Then Goursat's theorem asserts that f is analytic in an open complex domain Ω if and only if it satisfies the Cauchy–Riemann equation in the domain (Rudin 1966, Theorem 11.2). In particular, continuous differentiability of f need not be assumed (Dieudonné 1969, §9.10, Ex. 1).

The hypotheses of Goursat's theorem can be weakened significantly. If f = u + iv is continuous in an open set Ω and the partial derivatives of f with respect to x and y exist in Ω, and satisfies the Cauchy–Riemann equations throughout Ω, then f is holomorphic (and thus analytic). This result is the Looman–Menchoff theorem.

The hypothesis that f obey the Cauchy–Riemann equations throughout the domain Ω is essential. It is possible to construct a continuous function satisfying the Cauchy–Riemann equations at a point, but which is not analytic at the point (e.g., f(z) = z5 / |z|4). Similarly, some additional assumption is needed besides the Cauchy–Riemann equations (such as continuity), as the following example illustrates (Looman 1923, p. 107)

$f(z) = \begin{cases} \exp(-z^{-4})&\mathrm{if\ }z\not=0\\ 0&\mathrm{if\ }z=0 \end{cases}$

which satisfies the Cauchy–Riemann equations everywhere, but fails to be continuous at z = 0.

Nevertheless, if a function satisfies the Cauchy–Riemann equations in an open set in a weak sense, then the function is analytic. More precisely (Gray & Morris 1978, Theorem 9):

• If f(z) is locally integrable in an open domain Ω ⊂ C, and satisfies the Cauchy–Riemann equations weakly, then f agrees almost everywhere with an analytic function in Ω.

This is in fact a special case of a more general result on the regularity of solutions of hypoelliptic partial differential equations.

### Several variables

There are Cauchy–Riemann equations, appropriately generalized, in the theory of several complex variables. They form a significant overdetermined system of PDEs. As often formulated, the d-bar operator

$\bar{\partial}$

annihilates holomorphic functions. This generalizes most directly the formulation

${\partial f \over \partial \bar z} = 0,$

where

${\partial f \over \partial \bar z} = {1 \over 2}\left({\partial f \over \partial x} + i{\partial f \over \partial y}\right).$

### Bäcklund transform

Viewed as conjugate harmonic functions, the Cauchy–Riemann equations are a simple example of a Bäcklund transform. More complicated, generally non-linear Bäcklund transforms, such as in the sine-Gordon equation, are of great interest in the theory of solitons and integrable systems.

### Definition in Clifford algebra

In Clifford algebra the complex number $z = x+iy$ is represented as $z \equiv x + I y$ where $I \equiv \sigma_1 \sigma_2$. The fundamental derivative operator in Clifford algebra of Complex numbers is defined as $\nabla \equiv \sigma_1 \partial_x + \sigma_2\partial_y$. The function $f=u + I v$ is considered analytic ⇔ $\nabla f = 0$, which can be calculated in following way:

$0=\nabla f= ( \sigma_1 \partial_x + \sigma_2 \partial_y )(u + \sigma_1 \sigma_2 v)=\sigma_1 \partial_x u + \underbrace{\sigma_1 \sigma_1 \sigma_2}_{=\sigma_2} \partial_x v + \sigma_2 \partial_y u + \underbrace{\sigma_2 \sigma_1 \sigma_2}_{=-\sigma_1} \partial_y v =0$

Grouping by $\sigma_1$ and $\sigma_2$:

$\nabla f = \sigma_1 ( \partial_x u - \partial_y v) + \sigma_2 ( \partial_x v + \partial_y u) = 0$$\begin{cases} \partial_x u - \partial_y v = 0\\ \partial_x v + \partial_y u = 0 \end{cases}$

$\begin{cases} \dfrac{ \partial u }{ \partial x } = \dfrac{ \partial v }{ \partial y }\\ \dfrac{ \partial u }{ \partial y } = -\dfrac{ \partial v }{ \partial x } \end{cases}$

### Conformal mappings in higher dimensions

Let Ω be an open set in the Euclidean space Rn. The equation for an orientation-preserving mapping $f:\Omega\to\mathbf{R}^n$ to be a conformal mapping (that is, angle-preserving) is that

$Df^TDf = (\det(Df))^{2/n}I$

where Df is the Jacobian matrix, with transpose $Df^T$, and I denotes the identity matrix (Iwaniec & Martin 2001, p. 32). For n = 2, this system is equivalent to the standard Cauchy–Riemann equations of complex variables, and the solutions are holomorphic functions. In dimension n > 2, this is still sometimes called the Cauchy–Riemann system, and Liouville's theorem implies, under suitable smoothness assumptions, that any such mapping is a Möbius transformation.