Ch0 Total
Ch0 Total
Ch0 Total
Complex Numbers
A complex number z is an ordered pair ( x, y ) of real
numbers x and y, written
z = ( x, y ) .
x = Re z ,
y = Im z.
By definition, two complex numbers are equal if and only
if their real parts are equal and their imaginary parts are
equal.
In particular, (0,1) is called the imaginary unit and is
denoted by i ,
i = (0,1).
The following abbreviations are frequently used:
( x,0) = x,
(0, y ) = iy.
0- 1
Chapter 0
x1 x2 + y1 y2
x2 y1 x1 y2
.
i
+
2
2
2
2
x2 + y2
x2 + y2
0- 2
Chapter 0
Complex Plane
Complex numbers can be represented by points in a
Cartesian coordinate system, whose x-axis and y-axis are
respectively called the real and imaginary axes.
0- 3
Chapter 0
Addition
of
complex
numbers
Subtraction of complex
numbers
z = x iy.
Useful relations:
1
Re z = x = ( z + z ),
2
Im z = y =
1
( z z ).
2i
0- 4
Chapter 0
z = r (cos + i sin ).
PHY 4602 Methods in Theoretical Physics II
0- 5
Chapter 0
counterclockwise sense.
For z = 0 this angle is undefined.
For given z 0 it is determined only up to integer
multiples of 2 .
0- 6
Chapter 0
Triangle inequality
For any complex numbers we have the important triangle
inequality,
z1 + z2 z1 + z2 .
Proof:
The triangle inequality is a direct consequence of an
important theorem in geometry:
The
sum
of
the
z1 + z2 + + zn
0- 7
Chapter 0
z1 r1
z1
z1
=
.
z2
z2
z1
) = arg z1 arg z2 .
z2
De Moivres formula
z n = r n (cos n + i sin n ).
(cos + i sin ) n = cos n + i sin n .
0- 8
Chapter 0
Roots
If z = wn (n = 1, 2, ) , each of these possible values of w
is called an n-th root of z, and written as
w = n z.
+ 2 k
+ 2 k
z = n r cos
+ i sin
n
n
k = 0,1, n 1.
z is n-valued.
Example:
n
2k
2k
1 = cos
,
+ i sin
n
n
k = 0,1, n 1.
Some terminologies
Circle:
A circle of radius and centre a in the complex plane
is given by the solution of the equation:
za =
0- 9
Chapter 0
za < .
Closed circular disk
A closed circular disk is the interior of a circle, i.e.
z a < , plus the circle itself, given by the equation:
za .
Neighbourhood
0- 10
Chapter 0
1 < z a < 2 .
1 z a 2
Half-planes.
0- 11
Chapter 0
Open sets:
A set S is called open if every point of S has a
neighbourhood consisting entirely of points that
belong to S .
Domains:
An open set S is called connected if any two of its
points can be joined by a broken line of finitely many
straight-line segments all of whose points belong to
S . An open connected set is called a domain.
Complement:
The complement of a set S in the complex plane is
the set of all points of the complex plane that do not
belong to S .
Boundary points:
A boundary point of a set S is a point every
neighbourhood of which contains both points that
belong to S and points that do not belong to S .
PHY 4602 Methods in Theoretical Physics II
0- 12
Chapter 0
Complex Functions
Let S be a set of complex numbers. A complex function f
defined on S is a rule assigning to every z = x + iy in S a
complex number w, called the value of the function f at z
and denoted as:
w = f ( z ).
z varies in S, a set of complex numbers, and is hence called
a complex variable. The set S is called the domain of
definition of f. The set consisting of all values of function
f is called the range of f.
= x 2 + 2ixy y 2 + 3 x + 3iy
PHY 4602 Methods in Theoretical Physics II
0- 13
Chapter 0
u = Re f ( z ) = x 2 y 2 + 3 x , and
= Im f ( z ) = 2 xy + 3 y .
Consider w = f ( z ) = 2iz + 6 z .
f ( z ) = 2i ( x + iy ) + 6( x iy )
u ( x, y ) = 6 x 2 y , and
( x, y ) = 2 x 6 y .
Limit
Consider a function f ( z ) defined in a neighbourhood of
z0 (expect perhaps at z0 itself). If the value of f is close to
l for all z close to z0 , f ( z ) is said to have the limit l as z
approaches z0 , written
lim f ( z ) = l ,
z z0
Note:
z may approach z0 from any direction in the complex
plane and the way
0- 14
Chapter 0
z-plane
w-plane
z z0
Continuity
A function f ( z ) is said to be continuous at z = z0 if f ( z0 )
is defined and
lim f ( z ) = f ( z0 ).
z z0
0- 15
Chapter 0
Derivative
Definition:
The derivative of a complex function f at a point z0 is
written f ( z0 ) and is defined by
f ( z0 ) = lim
z z0
f ( z ) f ( z0 )
f ( z0 + z ) f ( z0 )
= lim
z 0
z
z z0
0- 16
Chapter 0
Differentiation rules:
(cf ) = cf ,
( f + g ) = f + g ,
( fg ) = f g + fg ,
f f g fg
.
=
2
g
g
They are the same as in real calculus.
Example:
Despite that z is obviously a differentiable function, its
complex conjugate, z , is not differentiable.
Proof:
f ( z + z ) f ( z ) ( z + z ) z z x iy
=
=
=
.
z
z
z x + iy
If y = 0 , the limit is +1.
If x = 0 , it is 1.
The limit is path-dependent and z is not differentiable.
PHY 4602 Methods in Theoretical Physics II
0- 17
Chapter 0
Analytic Functions
Definition of Analytic Functions
A function f ( z ) is said to be analytic in a domain D (or,
equivalently, holomorphic in D) if f ( z ) is defined and
differentiable at all points of D. The function f ( z ) is said
to be analytic at a point z = z0 in D if f ( z ) is analytic in a
neighbourhood of z0 .
Analytic functions are functions that are differentiable in
some domain.
Examples:
The nonnegative integer powers 1, z , z 2 , and also
the polynomials,
f ( z ) = c0 + c1 z + c2 z 2 + + cn z n
f ( z) =
g ( z)
h( z )
( a rational function),
0- 18
Chapter 0
Cauchy-Riemann Equations
Let f ( z ) = u ( x, y ) + iv( x, y ) be defined and continuous in
some
neighbourhood
of
point
z = x + iy
and
u y = vx
Proof:
By definition,
f ( z + z ) f ( z )
z 0
z
u ( x + x, y + y ) u ( x, y )
= lim
z 0
x + iy
v( x + x, y + y ) v( x, y )
.
+i lim
z 0
x + iy
f ( z ) = lim
0- 19
Chapter 0
u ( x + x, y ) u ( x, y )
x 0
x
v( x + x, y ) v( x, y )
.
+i lim
x 0
x
f ( z ) = lim
f ( z ) = u x + ivx .
z = iy (i.e. along a vertical path).
u ( x, y + y ) u ( x, y )
y 0
i y
v( x, y + y ) v( x, y )
.
+i lim
y 0
iy
f ( z ) = lim
f ( z ) = iu y + v y .
The two results should be the same, resulting in the
Cauchy-Riemann equations:
u x = v y ,
u y = vx
0- 20
Chapter 0
Examples:
Consider f ( z ) = z 2 = ( x + iy ) 2 = u + iv
u = x 2 y 2
v = 2 xy
Since f ( z ) = z 2 is analytic for all z, its CauchyRiemann equations are given by:
ux = 2 x = vy ;
and vx = 2 y = u y
Consider f ( z ) = z = x iy .
We have:
u = x, v = y .
u x = 1 v y = 1.
Conclusion: f ( z ) = z is not analytic.
domain
D,
then
the
complex
function
f ( z ) = u ( x, y ) + iv( x, y ) is analytic in D.
PHY 4602 Methods in Theoretical Physics II
0- 21
Chapter 0
Question:
Is f ( z ) = z 3 an analytic function?
Solution:
u = x 3 3 xy 2 and v = 3 x 2 y y 3 .
ux = 3x 2 3 y 2 ,
v y = 3x 2 3 y 2
u y = 6 xy,
vx = 6 xy.
u y = vx
are satisfied and hence f ( z ) = z 3 is analytic for every z.
Classwork:
Express z in terms of the polar coordinates:
z = r (cos + i sin )
and consider an analytic function:
f ( z ) = u (r , ) + iv(r , ) .
Show that the Cauchy-Riemann equations in the polar
coordinates are:
1
u
=
r r v ,
v = 1 u .
r
r
PHY 4602 Methods in Theoretical Physics II
0- 22
Chapter 0
Laplaces Equation
If f ( z ) = u ( x, y ) + iv( x, y ) is analytic in a domain D, then u
and v satisfy Laplaces equation
2 u = u xx + u yy = 0,
2 v = vxx + v yy = 0,
respectively, in D and have continuous second partial
derivatives.
Proof:
Differentiate the Cauchy-Riemann equations once:
=
u
x x x v y ,
u xx = v yx
u yy = vxy .
u y = ( vx )
y
y
One then obtains:
2 u = u xx + u yy = 0 .
The other Laplaces equation,
2 v = vxx + v yy = 0,
0- 23
Chapter 0
u
=
y x y v y ,
u xy = v yy
u yx = vxx .
u = ( v )
x
x y x
Harmonic Functions
Solutions of Laplaces equation having continuous second
partial derivatives are called harmonic functions and
relevant theory is called potential theory.
Note:
Real and imaginary parts of an analytic function are
harmonic functions.
If a pair of harmonic functions u and v satisfy the
Cauchy-Riemann equations:
u x = v y and u y = vx ,
in a domain D, they are the real and imaginary parts
of an analytic functions f in D. v is then said to be a
conjugate harmonic function of u in D.
PHY 4602 Methods in Theoretical Physics II
0- 24
Chapter 0
Example:
Show that u = x 2 y 2 y is a harmonic function and find
a conjugate harmonic function v of u.
Solution:
It is straightforward to show that 2u = 0 and hence
u ( x, y ) is indeed a harmonic function.
To find a conjugate harmonic function of u, we write
down the Cauchy-Riemann equations:
v y = u x = 2 x,
vx = u y = 2 y + 1.
vx = 2 y +
dh
dh
=1
= 2 y +1
dx
dx
h( x ) = x + c .
v = 2 xy + x + c
f ( z ) = u + iv = x 2 y 2 y + i (2 xy + x + c) = z 2 + iz + ic.
0- 25
Chapter 0
Note:
A conjugate of a given harmonic function is uniquely
determined up to an arbitrary real additive constant.
Conformal Mapping
Mapping
A given function f
f ( z ) = u ( x, y ) + iv( x, y )
assigns to each point z in its domain of definition D the
corresponding point
w = f ( z)
in the w-plane. We say that f defines a mapping of D into
the w-plane.
For each point z in D , the corresponding point w = f ( z )
is called the image of z with respect to the function f.
0- 26
Chapter 0
horizontal
lines
y = k = constant
are
0- 27
Chapter 0
Conformal Mapping
A conformal mapping is a mapping which preserves
angles between any two oriented curves (i.e. the angle
between the oriented tangents at the intersection point)
both in magnitude and in sense.
the mapping w = f ( z )
z -plane
Theorem of conformality
The mapping defined by an analytic function f ( z ) is
conformal, except at critical points where the derivative
f ( z ) is zero.
Example:
The
mapping
w = z n , n = 1, 2,3,
conformal, except at z = 0 .
0- 28
is
Chapter 0
Proof:
the mapping w = f ( z )
z -plane
1 = lim arg [ z1 (t ) z0 ]
t t0
2 = lim arg [ z2 (t ) z0 ]
t t0
0- 29
Chapter 0
z1 (t ) z0
= lim arg
t t0
z
t
z
(
)
2
0
( z1 (t ) z0 ) /(t t0 )
dz1 / dt
= lim arg
= arg
t t0
(
z
(
t
)
z
)
/(
t
t
)
dz
/
d
t
0
0
2
1* = lim arg [ f ( z1 (t )) f ( z0 ) ]
t t0
2* = lim arg [ f ( z2 (t )) f ( z0 ) ]
t t0
* = 1* *
2
f ( z1 (t )) f ( z0 )
= lim arg
t t0
f
(
z
(
t
))
f
(
z
)
2
0
( f ( z1 (t )) f ( z0 )) /(t t0 )
= lim arg
t t0
(
f
(
z
(
t
))
f
(
z
))
/(
t
t
)
2
0
0
df ( z1 ) / dt
= arg
d
f
(
z
)
/
dt
2
0- 30
Chapter 0
df ( z1 ) / dt
df
z
dt
(
)
/
2
* = arg
df
dz
= arg
df
dz
dz1
dt
dz2
dt
If f ( z ) 0 , then we have:
dz1
dt
*
= arg
dz
2
dt
0- 31
Chapter 0
Exponential Function
The complex exponential function is defined as:
e z = e x (cos y + i sin y ).
Properties:
e z is an entire function, that is, analytic for all z.
(e z ) = e z .
e z1 + z2 = e z1 e z2
Eulers formula:
ez 0
Periodicity:
e z + 2 i = e z
Parity:
e z + i = e z
0- 32
Chapter 0
Trigonometric Functions
For real variable we have:
eix = cos x + i sin x,
e ix = cos x i sin x.
Hence,
1
cos x = (eix + e ix ),
2
sin x =
1 ix ix
(e e ).
2i
1
,
cos z
1
csc z =
.
sin z
sec z =
0- 33
Chapter 0
Hyperbolic Functions
The hyperbolic cosine and sine functions for complex
variables are defined by:
1 z
(e + e z ),
2
1
sinh z = (e z e z ).
2
cosh z =
Similarly,
sinh z
,
cosh z
1
sech z =
,
cosh z
tanh z =
cosh z
,
sinh z
1
csc hz =
.
sinh z
coth z =
Differentiation:
(cosh z ) = sinh z ,
(sinh z ) = cosh z ,
0- 34
Chapter 0
Logarithm Function
The inverse function of the exponential function, e w = z , is
defined as the natural logarithm ln z , i.e.
w = ln z
Properties:
ln(e z ) = z
a = eln a
a z = e z ln a
ln( z1 z2 ) = ln z1 + ln z2
ln( z1 / z2 ) = ln z1 ln z2
(ln z ) =
1
z
ln( rei ) = ln r + i
The
complex
natural
logarithm ln z ( z 0) is
ln z = w i 2n
0- 35
Chapter 0
Principal value
Since
ln z = ln(| z | ei ) = ln | z | +i ,
where = arg z , one can define the principal value of ln z
by:
Ln z = ln z + iArg z .
Therefore,
ln z = Ln z 2n i .
Examples:
ln1 = 0, 2 i, 4 i,
Ln 1 = 0
ln i = i / 2, 3 i / 2,5 i / 2,
Ln i = i / 2
ln 4i = 1.386294 + i / 2 2n i
Ln 4i = 1.386294 + i / 2
ln(4i ) = 1.386294 i / 2 2n i
Ln (4i ) = 1.386294 i / 2
PHY 4602 Methods in Theoretical Physics II
0- 36
Chapter0
Ln z = Ln r + i ,
r > 0,
< .
(ii)
0- 37
Chapter 0
Domain D
0- 38
Chapter 0
Proof:
We write Ln z = u ( x, y ) + iv( x, y ) with
u = Ln r ,
v =,
and
u 1 v
=
, and
r r
1 u
v
=
,
r
r
u 1
= ,
r r
u
and v
1 v 1
= ,
r r
obviously
v
1 u
= 0,
= 0,
r
r
satisfy
the
Cauchy-Riemann
Note:
ln z is a multi-valued function in the entire z -plane.
Ln z is a single-valued function in the entire z -plane.
Ln z is an analytic single-valued function in domain
D.
0- 39
Chapter 0
Example:
Ln z is called the principal branch of the logarithmic
Branch cut
A line used to create a domain of analyticity is called a
branch line or branch cut.
0- 40
Chapter 0
Examples:
The principal value of the logarithmic function, Ln z ,
defined by:
Ln z = Ln r + i ,
r > 0,
< .
Domain D
The semi-infinite
line y = 0 , x 0
is the branch cut.
Another branch of ln z :
f ( z ) = Ln r + i ,
where 3 / 2 < / 2 .
is analytic in domain D1 .
The
infinite
semiline
x = 0 , y 0 is
0- 41
Chapter 0
Branch point
When we are trying to establish branches of a multi-valued
function, we find that all the possible branch cuts share a
common point. We call that point a branch point of the
multi-valued function.
Examples:
In the case of ln z the origin z = 0 is a branch point.
branch point z = 3 + 4i
PHY 4602 Methods in Theoretical Physics II
0- 42
Chapter 0
Power functions
In general, we can define z c , where c can be any real or
complex number, by:
z c = exp ( c ln z ) .
Since the logarithm of z is multi-valued, depending on the
value of c, z c may have more than one numerical value. If
we employ the principal value of the logarithm in the
above definition, we obtain the principal value of z c .
Examples:
Square-root function z1/ 2
1/ 2 )[ Ln z + i (2 k )]
z1/ 2 = e(
1/ 2 ) ln z
= e(
= e(
1/ 2 ) Ln z
= e(
1/ 2 ) Ln z + ik
[cos(k ) + i sin(k )]
= ( 1) e(
k
1/ 2 ) Ln z
0- 43
Chapter 0
( z z0 )
=e
c ln ( z z0 )
( z z0 ) .
c
line y = 0 , x 1.
r1 = z 1 ,
1 = arg ( z 1) ,
0 1 < 2 .
( z 1)
1/ 3
1/ 3
= r1ei ( + 2 k )
= 3 r1 e ( 1
i / 3 + 2 k / 3 )
0- 44
Chapter0
f ( z )dz.
z (t ) = x(t ) + iy (t ) ,
where t is a real parameter, and x(t ) and y (t ) are two real-
valued functions of t.
0- 45
Chapter 0
The curve C
S n = f ( m )zm
m =1
where
zm = zm zm1.
0- 46
Chapter 0
( m = 1, 2, , n )
f ( z )dz
, or by
f ( z ) dz
if C is a closed path.
0- 47
Chapter 0
[k f ( z) + k
C
1 1
f ( z ) ] dz = k1 f1 ( z )dz + k2 f 2 ( z )dz.
2 2
Directionality
z0
z0
f ( z )dz = f ( z )dz.
Z
Additivity
C2
0- 48
Chapter 0
Examples:
1. Integral of 1/ z around the unit circle C in the counterclockwise direction.
z (t ) = cos t + i sin t = eit
( 0 t 2 ) .
2
2
dz
it
it
=
e
(
ie
dt
)
=
i
C z 0
0 dt = 2 i.
2. Consider
dz
C z = 2 i.
( z z0 ) m dz .
( z z0 ) m = m eimt ,
( z z0 ) dz = i
m +1
= i
m +1
( 0 t 2 ) .
dz = i eit dt
ei ( m +1) t dt
2 i, m = 1
C ( z z0 ) dz = 0, m 1
m
0- 49
Chapter 0
(0 t 1)
(0 t 1)
z (t ) = 1,
f ( z (t )) = x(t ) = t
C2 : z (t ) = 1 + it ,
(0 t 2)
z (t ) = i,
f ( z (t )) = x(t ) = 1 .
(a)
(b)
Re zdz = t (1 + 2i )dt =
*
0
1
1
(1 + 2i ) = + i.
2
2
1
C2
1
+ 2i.
2
0- 50
Chapter 0
f ( z )dz ML ,
0- 51
Chapter 0
An open annulus:
An open circle:
za <
is
connected domain.
0- 52
Chapter 0
f ( z ) dz = 0.
Note:
f ( z ) is of course analytic in the domain enclosed by C.
Corollary:
If f ( z ) is analytic in the domain enclosed by a contour C,
f ( z ) dz = 0.
0- 53
Chapter 0
Examples:
cos z dz = 0,
e z dz = 0,
z n dz = 0.
sec z dz = 0
dz
=0
2
z +4
(i)
(ii)
z 1 dz = 2 i
0- 54
Chapter 0
f ( z )dz =
f ( z ) = u ( x, y ) + iv( x, y ) ,
and
Ia =
(udx vdy ),
Ib =
(udy + vdx).
v u
=
udx
v
dy
(
)
C
R x y dxdy
u v
C (vdx + udy) = R x y dxdy
0- 55
Chapter 0
v u
=
x y
u v
=
x y
v u
=
udx
v
dy
(
)
C
R x y dxdy = 0
u v
C (vdx + udy) = R x y dxdy = 0
f ( z )dz =
Theorem
If f ( z ) is analytic in a simply connected domain D, then
the integral of f ( z ) is independent of path in D.
0- 56
Chapter 0
Proof:
Consider two directed paths C1 and C2 in D from z1 to z2 .
Reverse the direction of C2 to obtain the path C2* .
C2
fdz = * fdz.
Apply
Cauchys
integral
C2
C1
fdz + * fdz = 0,
C2
Therefore,
C1
C1
fdz = * fdz.
C2
f ( z )dz = f ( z ) dz
C2
C1
f ( z )dz = f ( z )dz =
C2
z2
z1
f ( z ) dz
0- 57
Chapter 0
Deformation of path
If
path
deforms
C1
continuously to path C2
without crossing any point
at which
f ( z ) is non-
analytic, then
C1
f ( z )dz = f ( z )dz
C2
z1
z0
f ( z )dz = F ( z1 ) F ( z0 ) .
0- 58
Chapter 0
Examples:
1+ i
1 3
z dz = z
3
2
1+ i
0
1
2 2
3
= (1 + i ) = + i
3
3 3
i
i
= 2sin( i ) = 23.097i
Ln
(
)
=
) = i .
i z
2
2
i
Proof:
Consider the line integral from z0 to z :
z
F ( z ) = f ( z * )dz *
z0
z +z
z0
f ( z * )dz *
z
F ( z + z ) F ( z ) 1 z + z
*
*
*
*
(
)
(
)
f
z
dz
f
z
dz
=
z0
z
z z0
1 z +z
*
*
(
)
.
=
f
z
dz
z
z
0- 59
D.
Chapter 0
f ( z )z ]
[
F ( z + z ) F ( z )
= lim
= f ( z ).
F ( z ) = lim
z 0
0
z
z
Therefore, F ( z ) is an indefinite integral of f ( z ) , written
F ( z ) = f ( z )dz.
There are many possible indefinite integrals for f ( z ) ,
corresponding to different choices of the starting point:
f ( z * )dz * = Fj ( z ) .
zj
Fj ( z ) Fk ( z ) = f ( z )dz f ( z * ) dz *
*
zj
zk
zk
zj
= f ( z * )dz * +
f ( z * )dz * =
zk
zj
f ( z * )dz *
Hence, if by definition
z
z1
f ( z )dz = G ( z )
F ( z1 ) = +C
G ( z ) = F ( z ) F ( z1 )
z2
z1
f ( z )dz = F ( z2 ) F ( z1 )
0- 60
Chapter 0
f ( z)
dz = 2 if ( z0 ) , or
z z0
1
f ( z0 ) =
2 i
f ( z)
dz
z z0
Proof:
f ( z)
dz
dz = f ( z0 )
+
C
z z0
z z0
f ( z0 )
f ( z ) f ( z0 )
dz.
z z0
dz
= f ( z0 ) 2 i .
z z0
f ( z ) f ( z0 )
dz = 0
z z0
0- 61
Chapter 0
Arguments:
(i)
Since f ( z ) is analytic at z = z0 ,
lim
z z0
f ( z ) f ( z0 )
= f ( z0 ) .
z z0
Hence
f ( z ) f ( z0 )
is bounded if z z0
z z0
is
(iii)
f ( z ) f ( z0 )
dz =
z z0
z z0 =
f ( z ) f ( z0 )
dz .
z z0
0- 62
Chapter 0
(iv)
z z0 =
f ( z ) f ( z0 )
dz < M (2 ) < M (2 )
z z0
(v)
z z0 =
f ( z ) f ( z0 )
dz = 0 =
z z0
f ( z ) f ( z0 )
dz
z z0
f ( z)
= 2 if ( z0 )
z z0
Examples:
ez
dz = 2 ie z
z2
z3 6
dz =
2z i
z =2
= 2 ie 2 46.4268i.
z3 3
dz
1
z 2i
1
2
C
= 2 i 12 z 3 3
z =i / 2
6 i.
C : z = 3/ 2 in counter-clockwise direction.
Use method of partial fractions to expand ( z 2 1) 1 .
PHY 4602 Methods in Theoretical Physics II
0- 63
Chapter 0
1
1 1
1
=
.
2
z 1 2 z 1 z +1
tan z
1 tan z
tan z
dz
dz
dz
=
C
C
2 z 1
z +1
z 1
2 i
=
tan1 tan( 1)] = 2 i tan1 9.785i.
[
2
1
2 i
C1
f ( z)
1
dz +
2 i
z z0
C2
f ( z)
dz ,
z z0
where the outer integral (over C1 ) is taken along counterclockwise direction and the inner integral (over C2 ) is
taken along clockwise direction.
PHY 4602 Methods in Theoretical Physics II
0- 64
Chapter 0
Proof:
By adding two equal but
opposite curves, C3 and C4 ,
we
can
join
the
two
f ( z0 ) =
1
2 i
1
2 i
1
=
2 i
f ( z)
dz
z z0
C1 + C2 + C3 + C4
C1
f ( z)
dz
z z0
f ( z)
1
dz +
z z0
2 i
C2
f ( z)
dz.
z z0
0- 65
Chapter 0
f ( z0 ) =
1
2 i
f ( z)
C ( z z0 )2 dz,
f ( z0 ) =
2!
2 i
f ( z)
C ( z z0 )3 dz.
f ( n ) ( z0 ) =
n!
2 i
f ( z)
C ( z z0 )n +1 dz.
0- 66
Chapter 0
Proof:
f ( z0 ) = lim
z 0
f ( z0 + z ) f ( z0 )
.
z
f ( z0 + z ) f ( z0 )
z
f ( z)
1
dz
=
2 iz C z ( z0 + z )
=
1
2 i
f ( z)
dz
z z0
f ( z)
C ( z z0 z )( z z0 ) dz
1
2 i
f ( z)
C ( z z0 )2 dz
0- 67
Chapter 0
Examples:
Let C be a contour enclosing the point z = i
z 4 3z 2 + 6
2 i 4
2
dz
=
(
z
3
z
+ 6)
3
2!
( z + i)
= i 12 z 2 6
Find
z = i
z = i
= 18 i.
cos z
C ( z 1)3 ( z 5)2 dz , where C is the circle
z 4 = 2.
z 1 is nonzero both inside and on C.
cos z
( z 1)3 dz
64sin 5 48cos 5
d cos z
=
2
i
dz ( z 1)3 z = 5
46
0- 68
Chapter 0
f ( z )dz = 0
f ( z0 ) =
1
2
f (z0 + rei )d ,
0- 69
Chapter0
Complex Sequences
A complex sequence z1 , z2 , is a sequence of complex
numbers.
Convergent sequence
A complex convergent sequence z1 , z2 , is a sequence of
complex numbers that has a finite limit c, written
lim zn = c
n
In other words, as n , zn c.
Rigorous definition:
For every > 0 , no matter how small it is, we can always
find a sufficiently large N > 0 such that
zn c <
Chapter0
Examples:
The
sequence
{i / n} = {i, 1/ 2, i / 3,1/ 4, }
n
is
zn = xn + iyn ,
where n = 1, 2, ,
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-71
Chapter0
Complex Series
Given a complex sequence z1 , z2 , , zm , , we can form
another sequence of sums,
s1 = z1 ,
s2 = z1 + z2 ,
s3 = z1 + z2 + z3 ,
sn = z1 + z2 + + zn .
sn is called the n-th partial sum of the infinite series and
z1 , z2 , are called the terms of the series.
Convergent series
A (complex) convergent series is a series whose sequence
of partial sums converges,
lim sn = s.
n
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-72
Chapter0
Theorem
A complex series with zm = xm + iym converges and has the
sum s = u + iv if and only if
x1 + x2 + converges and has the sum u ; and
y1 + y2 + converges and has the sum v.
lim zm 0 ,
Absolute convergence
A series z1 + z2 + is called absolutely convergent if the
series of the absolute values of the terms
z
m =1
= z1 + z2 +
is convergent.
Theorem:
If a series is absolutely convergent, it is convergent.
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-73
Chapter0
Conditional convergence
If the series z1 + z2 + converges while the series
z1 + z2 + diverges, then the series z1 + z2 + is called
conditionally convergent.
Example:
1 1 1
The series 1 + + converges, and is a
2 3 4
conditionally
convergent
z1 + z2 + = 1 +
series
since
the
series
1 1 1
+ + + diverges.
2 3 4
zn bn
for n = 1, 2, ,
Chapter0
= 1 + q + q 2 +
m =0
1 q n +1
1
q n +1
sn = q =
=
.
1 q
1 q 1 q
m=0
n
Ratio test 1
If a series z1 + z2 + with zn 0 (n = 1, 2, ) has the
property that for every n > N ,
Chapter0
Ratio test 2
If a series z1 + z2 + with zn 0 (n = 1, 2, ) is such that
lim
n
zn+1
= L , then:
zn
Example:
Consider the series
(100 + 75i ) n
1
=
1
+
(100
+
75
i
)
+
(100 + 75i ) 2 +
n!
2!
n =0
n +1
zn +1
=0 .
zn
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-76
Chapter0
Let
z n = i / 23 n / 2
, n = 0, 2, 4,
zn = 1/ 2(3n 1) / 2
, n = 1,3,5,
1 i 1
i
1
z0 + z1 + z2 + z3 + = i + + + + +
+
2 8 16 64 128
Root test 1
If a (complex) series z1 + z2 + is such that for every n
greater than some N ,
n
zn q < 1
zn 1
Chapter0
Root test 2
If a series z1 + z2 + is such that lim n zn = L , then:
n
Power series
A power series in powers of z z0 is a series of the form
a (z z )
n=0
= a0 + a1 ( z z0 ) + a2 ( z z0 ) 2 +
n
2
a
z
=
a
+
a
z
+
a
z
+
0
0
1
2
n =0
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-78
Chapter0
Examples:
Consider the geometric series
= 1 + z + z 2 +
n =0
z n +1
= z.
n
z
zn
z 2 z3
= 1 + z + + +
2! 3!
n =0 n !
is absolutely convergent for every z . By the ratio test,
z
z n+1 /(n + 1)!
=
0
n
z / n!
n +1
as
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-79
Chapter0
Theorem
Consider a power series of the form:
n
2
a
(
z
z
)
=
a
+
a
(
z
z
)
+
a
(
z
z
)
+
0
0
0
1
0
2
0
n=0
z z0 < z1 z0 .
c) If it diverges at a z = z2 , it diverges for every z farther
away from z0 than z2 .
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-80
Chapter0
n
2
a
(
z
z
)
=
a
+
a
(
z
z
)
+
a
(
z
z
)
+
0
0
0
1
0
2
0
n=0
z z0 < R
that is centered at center z0 and with radius R, and, in
addition, includes all the points at which the series
converges. The circle
z z0 = R
Chapter0
a (z z )
n=0
= a0 + a1 ( z z0 ) + a2 ( z z0 ) 2 +
R=
an
1
=
lim
.
*
n
L
an+1
(Cauchy-Hadamard formula).
Chapter0
Method 2
Suppose that the sequence of
an , ( n = 1, 2, ) ,
R=
1
1
.
=
n
L lim an
n
R = lim
(2n)!
n
(
3
)
z
i
2
n = 0 ( n !)
an
(2n)! (2n + 2)!
= lim
an +1 n (n !) 2 ((n + 1)!) 2
(n + 1) 2
1
= lim
= .
n (2 n + 2)(2n + 1)
4
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-83
Chapter0
f ( z ) = an ( z z0 ) n = a0 + a1 ( z z0 ) + a2 ( z z0 ) 2 +
n=0
Example:
For z < 1, the power series 1 + z + z 2 +
converges and is
equal to (1 z ) 1 .
Theorem
Continuity of the sum of a power series
The function f ( z ) defined by
f ( z ) = an ( z z0 ) n = a0 + a1 ( z z0 ) + a2 ( z z0 ) 2 +
n=0
z = z0 .
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-84
Chapter0
Theorem
Uniqueness of representation
an ( z z0 )
n
b
(
z
z
)
n
0
and
n=0
n=0
Corollary:
If a function f ( z ) has a power series representation with
any center z0 , this representation is unique.
Example:
The function (1 z ) 1 is represented by the power series
1+ z + z2 +
= (1 z ) 1
(1 z ) 1 = [1 ( z ) ]
1
=
1
z z 2
1 +
+
+
1 1
Chapter0
k
a
(
z
z
)
= a0 + a1 ( z z0 ) + = f ( z )
k
0
k =0
b (z z )
k =0
= b0 + b1 ( z z0 ) + = g ( z )
k
(
a
b
)(
z
z
)
k k
0
k =0
k
(
a
b
)(
z
z
)
= f ( z) g ( z)
k k
0
k =0
( a
k =0
bk )( z z0 ) k = f ( z ) g ( z )
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-86
Chapter0
Termwise multiplication:
Termwise multiplication of the two power series means
the multiplication of each term of the first series by each
term of the second series and the collection of like powers
of z . That is,
s ( z ) = ak ( z z0 ) k bm ( z z0 ) m
k =0 m=0
= (a0 bn + a1bn 1 + + an b0 )( z z0 ) n .
n=0
k
s ( z ) = ak ( z z0 ) bk ( z z0 ) k
k =0
k =0
.
= f ( z) g ( z)
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-87
Chapter0
Termwise differentiation:
Termwise differentiation of a power series,
d
n
n 1
a
z
z
n
a
z
z
(
)
(
)
n
n
0
0
dz
n =1
n =1
= a1 + 2a2 ( z z0 ) + 3a3 ( z z0 ) 2 +
gives the derived series of the original one. It is
permissible in the sense that
d
d
n
n
(
)
(
)
a
z
z
a
z
z
n
0
0
n
dz n =1
n =1 dz
= nan ( z z0 ) n 1
n =1
n an
(n + 1) an +1
a
n
lim n
n ( n + 1) n a
n +1
= lim
= lim
an
an +1
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-88
Chapter0
Termwise integration:
Termwise integration of the series
a0 + a1 ( z z0 ) + a2 ( z z0 ) 2 + yields a new series:
an
a1
n +1
2
z
z
F
z
a
z
z
z
z
(
(
)
=
(
)
+
(
)
+ ,
0
0
0
0
2
n=0 n + 1
which has the same radius of convergence as the original
series, and F ( z ) = f ( z ).
a (z z )
n=0
= f ( z ) , with a radius of
df
= nan ( z z0 ) n 1 , etc.
dz n =1
Chapter0
f ( z ) = an ( z z0 ) n
n =0
where an =
1 (n)
f ( z0 )
n!
(n)
n!
( z0 ) =
2 i
f ( z* )
*
dz
C ( z* z0 )n +1 .
The counterclockwise
integration is taken along
a simple closed curve C
that encloses z0 and f ( z )
is analytic in the interior
of C.
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-90
Chapter0
f ( z* )
*
dz
C ( z* z0 )n +1
f ( z* )
*
dz
C ( z* z0 )n +1 ( z* z )
Note:
A Taylor series is also called a Maclaurin series if the
center z0 = 0 .
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-91
Chapter0
Proof:
From theory of contour
integration, we know that:
1
f ( z) =
2 i
f ( z* ) *
dz .
z* z
Note:
z z0
< 1 for all z in the
*
z z0
interior of the circle C.
1 r n +1
+r =
,
1 r
n
1
:
we rewrite the quotient *
z z
1
z* z
=
1
=
*
z z0 ( z z0 )
1
( z * z0 )(1
z z0
)
*
z z0
n
n +1
z z0
z z0
1
1 z z0
1 + *
= *
+ + *
+ *
*
z z 0 z z0
z
z
z
z
z
z
0
0
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-92
Chapter0
Hence,
1
f ( z) =
2 i
+
f ( z * ) * z z0
dz +
2 i
z * z0
( z z0 )
2 i
f ( z* )
*
dz
C ( z* z0 )2 +
f (z )
*
dz
C ( z* z0 )n +1 + Rn ( z )
where
( z z0 ) n +1
Rn ( z ) =
2 i
f ( z* )
*
dz
n
+
*
1
*
C ( z z0 ) ( z z ) .
z z0
( z z0 ) 2
f ( z ) = f ( z0 ) +
f ( z0 ) +
f ( z0 ) +
1!
2!
( z z0 ) n ( n )
+
f ( z0 ) + Rn ( z )
n!
by evaluating the integrals with the theory of contour
integral. Furthermore, it can be shown that
lim Rn ( z ) = 0.
n
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-93
Chapter0
Remarks:
The limit lim Rn ( z ) = 0 can be obtained as follows:
n
( z z0 ) n +1
Rn ( z ) =
2 i
f ( z* )
*
dz
C ( z* z0 )n +1 ( z* z )
( z z0 ) n +1
2 z z0
*
n +1
,
ML
f ( z* )
where M is the maximum value of *
on the
z z
Since
z z0
z z0
*
n +1
n +1
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-94
Chapter0
f ( z* )
M
*
dz
C ( z* z0 )n +1
rn
1
= z n = 1 + z + z 2 +
1 z n=0
( z < 1).
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-95
Chapter0
Exponential function
zn
z2
e = = 1 + z + + .
2!
n =0 n!
z
z 2n
z2 z4
cos z = (1)
= 1 + +
(2n)!
2! 4!
n =0
n
z (2 n +1)
z3 z5
sin z = (1)
= z + +
(2n + 1)!
3! 5!
n =0
z 2n
z2 z4
cosh z =
= 1 + +
2! 4!
n = 0 (2n )!
z (2 n +1)
z3 z5
sinh z =
= z + +
3! 5!
n = 0 (2n + 1)!
Logarithm
z 2 z3
Ln (1 + z ) = z + +
2 3
( z < 1).
_____________________________________________________________________
PHY 4602 Methods in Theoretical Physics II
0-96
Chapter0
Theorem
Laurents theorem
bn
n
n =1 ( z z0 )
= an ( z z0 ) +
n
n=0
= a0 + a1 ( z z0 ) + a2 ( z z0 ) 2 +
+
b1
b2
+
+
2
z z0 ( z z0 )
0- 97
Chapter 0
1
f ( z* )
*
an =
dz
,
n +1
*
v
C
2 i ( z z0 )
n = 0,1, 2,3,"
1
*
n 1
*
*
bn =
z
z
f
z
dz
(
)
(
)
,
0
C
2 i v
n = 1, 2,3,"
Note:
The domain of validity of this expansion can be
0- 98
Chapter 0
a (z z )
n =
1
f ( z* )
*
an =
dz
C ( z* z0 )n +1
2 i v
with
(n = 0, 1, 2, ).
Proof:
(a). Apply the Cauchys integral formula for multiply
connected domains to f ( z ) .
If f ( z ) is analytic on C1 and
C2 and in domain D bounded
Cauchys
integral
formula
reads:
PHY 4602 Methods in Theoretical Physics II
0- 99
Chapter 0
f ( z) =
1
f ( z) *
1
f ( z) *
dz
+
dz ,
*
*
v
C
C
2 i 1 z z
2 i 2 z z
where the outer integral (over C1 ) is taken along counterclockwise direction and the inner integral (over C2 ) is
taken along clockwise direction.
In this case,
f ( z ) = g ( z ) + h( z )
1
f ( z* ) *
=
C1 z* z dz
2 i v
1
f ( z* ) *
dz
*
v
C
2
2 i
z z
Note: Both the integrals along C1 and C2 are taken in the
counter-clockwise direction.
1
f ( z* ) *
n
g ( z) =
dz
=
a
z
z
(
)
n
0
C1 ( z* z )
2 i v
n=0
where
1
f ( z* )
*
an =
dz
.
n
+
*
1
v
C
1
2 i
( z z0 )
0- 100
Chapter 0
C
2 i
z z
We note that
z z0
< 1 for all z in
*
z z0
the annulus.
Here C is the outer-circle C1 .
1
z* z
=
1
=
*
z z0 ( z z0 )
1
( z * z0 )(1
z z0
)
*
z z0
n
n +1
z z0
z z0
1
1 z z0
1 + *
= *
+ + *
*
+ *
z z 0 z z0
z
z
z
z
z
z
0
0
Therefore,
0- 101
Chapter 0
f ( z * ) * z z0
1
g ( z) =
dz +
*
v
C
2 i
2 i
z z0
+
( z z0 )
2 i
f ( z* )
*
dz
v C ( z* z0 )2 +
f (z )
*
dz
v C ( z* z0 )n +1 + Rn ( z )
where
( z z0 ) n +1
Rn ( z ) =
2 i
f ( z* )
*
dz
v C ( z* z0 )n +1 ( z* z ) .
f ( z* ) *
1
n
g ( z) =
dz
=
a
z
z
with
(
)
n
0
*
v
C
1
2 i
( z z)
n=0
f ( z* )
1
*
an =
dz
C1 ( z* z0 )n +1 .
2 i v
(c). For the integral along the inner circle we have
z * z0
1
< 1. Therefore, we rewrite *
as:
z z0
z z
1
1
=
=
*
*
z z z z0 ( z z 0 )
1
z z0
( z z0 )(1
)
z z0
*
0- 102
Chapter 0
z z0
1 z z 0 z z 0
1
=
+
1 +
+
*
z z0 z z0
z
z
z z z z0
0
1 z * z0
z z * z z0
n +1
f ( z* ) *
1
h( z ) =
dz
*
v
C
2 i 2 z z
1 1
*
*
=
f
z
dz
(
)
C2
2 i z z0 v
+
1
( z z0 ) 2
C2
( z * z0 ) f ( z * )dz *
1
+ +
( z z0 ) n +1
Since
+
f
z
dz
(
z
z
)
(
)
R
( z)
n
0
v C2
z * z0
< 1 , we can also show that the
z z0
remainder
1
Rn* ( z ) =
2 i ( z z0 ) n +1
( z * z0 ) n +1
*
*
(
)
f
z
dz
v C2 z z*
0- 103
Chapter 0
Examples:
Find the Laurent series of z 5 sin z with center z = 0 .
1 (1) n 2 n +1
z sin z = 5
z
z n = 0 (2n + 1)!
1
1
1
1 2
= 4 2+
z +
6 z 120 5040
z
5
z 2 e1/ z = z 2 1 +
+
+
2
1! z 2! z
1 1
1
= z2 + z + +
+
+
2
2 3! z 4! z
PHY 4602 Methods in Theoretical Physics II
( z > 0)
0- 104
Chapter 0
z > 1.
1
(a)
= zn
1 z n =0
(b)
if z < 1.
If z > 1, we have
1
1
=
1 z z (1 z 1 )
1 1
1 1
= n = 2
z n=0 z
z z
(b) z > 1.
(a) If 0 < z < 1, we have:
0- 105
Chapter 0
1
1
1 n
= 3
= 3 z
3
4
z z
z (1 z ) z n = 0
1 1 1
= 3 + 2 + + 1 + z
z
z
z
(b) If z > 1
1
1
=
z 3 z 4 z 4 (1 z 1 )
1 1
1 1
= 4 n = 4 5
z n=0 z
z
z
2 z + 3
The function f ( z ) = 2
is singular at z = 1, 2 .
z 3z + 2
Find its Laurent series with center z = 0 if (a) 1 > z ;
(b) 2 > z > 1; and (c) z > 2 .
Note:
f ( z) =
1
1
.
z 1 z 2
1
= zn
1 z n =0
if z < 1.
1
1 1
= 2
1 z
z z
if z > 1.
0- 106
Chapter 0
1
1
1 n
=
=
z if z < 2 .
z 2 2(1 1 z ) n = 0 2n +1
2
1
1
2n
=
= n +1 if z > 2 .
2
z2
n=0 z
z (1 )
z
Region I: z < 1
Region II: 1 < z < 2
Region III: z > 2
1
3 5
9
f ( z ) = 1 + n+1 z n = + z + z 2 + .
2
2 4
8
n =0
(b) For 1 < z < 2 ,
0- 107
Chapter 0
1 n 1
f ( z ) = n +1 z n+1
n =0 2
n =0 z
1 1
1
1 1
= + z + z 2 + 2 .
2 4
8
z z
(c) For z > 2 ,
f ( z ) = (2n + 1)
n =0
1
z n+1
2 3 5 9
= 2 3 4 .
z z
z
z
Singular functions
Singularity
A function f ( z ) is singular or has a singularity at a point
z = z0 if f ( z ) is not analytic (or not even defined) at
z = z0 , but every neighbourhood of z = z0 contains points
at which f ( z ) is analytic.
Isolated singularity
z = z0 is called an isolated singularity of f ( z ) if there
Examples:
0- 108
Chapter 0
(2n + 1)
,
2
n = 0,1, 2,".
2
1
, n = 0,1, 2,". The
tan( ) is singular at z =
(2n + 1)
z
singularities
densely
populate
around
z=0 .
bn
,
n
n =1 ( z z0 )
f ( z ) = an ( z z0 ) +
n
n=0
0- 109
Chapter 0
Example:
1/ z
1
1
1
1
=
+
+
+
n
2
z 2! z
n=0 n! z
0- 110
Chapter 0
0- 111
Chapter 0
f ( z ) = an ( z z0 ) n + an +1 ( z z0 ) n +1 +
= ( z z0 ) n an + an +1 ( z z0 ) + an + 2 ( z z0 ) 2 +
where an 0 .
Theorem
Let f ( z ) be analytic at z = z0 and has a zero of n-th order
at z = z0 . Then 1/ f ( z ) has a pole of n-th order at z = z0 .
0- 112
Chapter 0
f ( z ) = an ( z z0 ) n +
n =0
b1
b2
+
+
2
z z0 ( z z0 )
C
2 i
Since we can obtain Laurent series (and hence b1 ) by
other methods, we can in fact use the above formula
to evaluate the integral, namely:
f ( z )dz = 2 ib1 .
Examples:
sin z
Integrate f ( z ) = 4 counter-clockwise around the
z
contour C : z = 1.
Expand the function about z = 0 :
0- 113
Chapter 0
sin z 1
1
z z3
+ +
f ( z) = 4 = 3
z
z 3! z 5! 7!
b1 = 1/ 3!.
sin z
i
dz
2
ib
.
=
1
v C z 4
3
Integrate
circle C: z = 1/ 2 .
f ( z ) is singular at z = 0 and z = 1 . We have to
seek the Laurent series that converges for 0 < z < 1.
From previous calculation we have:
1
1 1 1
=
+ 2 + + 1 +
3
4
3
z
z z
z
z
dz
= 2 i Res [ f ( z = 0)] = 2 i
3
4
z z
Evaluation of residues
PHY 4602 Methods in Theoretical Physics II
0- 114
Chapter 0
b1
+ a0 + a1 ( z z0 ) + a2 ( z z0 ) 2 + ,
z z0
we have
( z z0 ) f ( z ) = b1 + ( z z0 ) [ a0 + a1 ( z z0 ) + ] ,
and hence
Res f ( z ) = b1 = lim( z z0 ) f ( z ).
z = z0
z z0
Example:
Re s
z =i
9z + i
9z + i
=
lim
z
i
(
)
z ( z + i )( z i )
z ( z 2 + 1) z i
9z + i
10i
=
= 2 = 5i.
z
z
i
(
+
)
z =i
p( z )
, where p ( z0 ) 0 and q ( z ) has a simple
q( z )
0- 115
Chapter 0
Re s f ( z )
z = z0
= lim ( z z0 )
z z0
p( z )
q( z )
( z z0 ) p ( z )
z z0 ( z z ) q ( z ) + ( z z ) q ( z ) / 2 +"
]
0 [
0
0
0
= lim
Re s f ( z ) = Re s
z = z0
z = z0
p ( z ) p ( z0 )
=
.
q ( z ) q ( z0 )
Example:
9z + i 9z + i
10i
9z + i
=
=
=
= 5i.
Re s 2
3
z =i
z ( z + 1) ( z + z ) z =i 3 z + 1 z =i 2
bm
bm 1
b2
+
+
+
( z z0 ) m ( z z0 ) m 1
( z z0 ) 2
b1
+ a0 + a1 ( z z0 ) +
z z0
0- 116
Chapter 0
( z z0 ) m f ( z )
g ( z)
= bm + bm 1 ( z z0 ) + + b2 ( z z0 ) m 2
+b1 ( z z0 ) m 1 + a0 ( z z0 ) m + a1 ( z z0 ) m +1 +
This is just the Taylor series of g ( z ) . Therefore,
b1 =
1
g ( m 1) ( z0 )
(m 1)!
d m 1
1
m
Re s f ( z ) =
lim
( z z0 ) f ( z ) .
z = z0
( m 1)! z z0 dz m 1
f ( z) =
Example:
50 z
( z + 4)( z 1) 2
d
d 50 z
( z 1) 2 f ( z ) = lim
= 8.
z 1 dz
z 1 dz
z+4
Re s f ( z ) = lim
z =1
Residue integration
Residue theorem
Let f ( z ) be analytic inside a simple closed curve C and
on C, except for finitely many singular points z1 , z2 , , zk
inside C. Then the integral of f ( z ) taken counter-
0- 117
Chapter 0
f ( z )dz = 2 i Res f ( z ).
j =1
z=zj
Proof:
Let
Cn (n = 1, 2,.., k )
be a tiny circle
enclosing
the
singularity zn .
f ( z )dz + v
f ( z )dz + v f ( z )dz + + v
C1
C2
Ck
f ( z )dz = 0,
0- 118
Chapter 0
f ( z )dz = v
f ( z )dz + v f ( z )dz + + v
C1
C2
Ck
f ( z )dz
Cj
f ( z )dz = 2 i Res f ( z ),
z=zj
k
f ( z )dz = 2 i Res f ( z ).
j =1
z=zj
Examples:
Evaluate the integral
4 3z
v C z 2 z dz where C is a simple
Res
z =0
4 3z 4 3z
=
= 4,
z ( z 1) z 1 z =0
Res
z =1
4 3z 4 3z
=
= 1.
z ( z 1) z z =1
0- 119
Chapter 0
tan z
dz around the circle
2
z 1
C: z = 3/ 2 .
The circle C encircles two singularities of the
integrand, namely, z = 1.
tan z
v C z 2 1dz
tan z
tan z
= 2 i Res 2
+ Res 2
z =1 z 1 z = 1 z 1
tan z
= 2 i
2z
tan z
+
2z
z =1
z = 1
= 2 i tan1.
Evaluate the integral
ze z
/z
v C z 4 16 + ze dz ,
where C is the ellipse 9 x 2 + y 2 = 9 in the counterclockwise direction.
PHY 4602 Methods in Theoretical Physics II
0- 120
Chapter 0
ze
2
3
= z 1 + +
+
+
z 2! z 2 3! z 3
2 1
= z + +
+
2 z
ze z
/z
2
1
+
=
)i
ze
dz
4
v C z 4 16
0- 121
Chapter0
I = F (cos ,sin )d .
0
Since
1 i
1
1
(e + e i ) = ( z + ) ,
2
2
z
1
1
1
sin = (ei e i ) = ( z ) ,
2i
2i
z
cos =
dz / d = iei = iz ,
d = dz / iz ,
I=
dz
f ( z)
iz
0- 121
Chapter 0
Example:
I =
d
2 cos
dz / iz
=
1
1
2 (z + )
2
z
dz
i
( z 2 2 2 z + 1)
2
2
dz
.
=
C
i
( z 2 1)( z 2 + 1)
C
There are two simple poles in the integrand. The first one,
z1 = 2 + 1, lies outside the unit circle C. The second one,
z2 = 2 1, is inside C with the residue given by:
Res
z = z2
1
1
=
( z 2 1)( z 2 + 1) z 2 + 1 z =
1
= .
2
2 1
1
C ( z 2 1)( z 2 + 1) = 2 i 2
dz
I = ( 2 / i ) [ 2 i ( 1/ 2) ] = 2
0- 122
Chapter 0
f ( x)dx = lim
a a
If both lim
f ( x) dx is defined as:
b
a a
b 0
R R
principal value of
f ( x)dx = lim
R R
f ( x)dx.
f ( x)dx .
integral
f ( z )dz
0- 123
Chapter 0
f ( z )dz = f ( z )dz +
S
R
R
f ( x)dx = 2 i Re s f ( z )
k
z
Hence,
f ( z )dz <
k
k
R
=
.
2
R
R
f ( x)dx = 2 i Res f ( z ) .
0- 124
Chapter 0
Example:
Evaluate the improper integral I =
dx
.
4
1+ x
Solution:
By symmetry we have I =
dx
1 dx
=
.
4
4
2 1+ x
1+ x
z2 = e3 i / 4 ,
z3 = e 3 i / 4 ,
z 4 = e i / 4 .
1
1 3 i / 4
1 i / 4
1
= 3
= e
= e ,
Res f ( z ) =
4
z = z1
4
(1 + z ) z = z1 4 z z = z1 4
1
1 9 i / 4 1 i / 4
1
=
=
= e
Res f ( z ) =
.
e
4
3
z = z2
4
(1 + z ) z = z2 4 z z = z2 4
dx
2 i i / 4 i / 4
=
e e
(
)
1 + x 4
4
2 i
=
2i sin = sin =
.
4
4
4
2
1 dx
dx
=
=
.
4
4
1+ x
2 1+ x
2 2
0- 125
Chapter 0
f ( x) cos x dx
and
f ( x) sin x dx .
Jordans Lemma.
If as R , f ( z ) 0 uniformly in for 0 < < 2 ,
lim
R SU
f ( z )ei z dz = 0
lim f ( z )e i z dz = 0
R S L
f ( z ) 1,
for = 0
, for < < 2
0- 126
Chapter 0
Applications:
If as R , f ( z ) 0 uniformly in for 0 < < 2 ,
then
By closing the contour in the upper-half plane,
where
f ( z )e
i z
upper
upper-half plane.
By closing the contour in the lower-half plane,
where
f ( x)e i x dx = 2 i Res f ( z )e i z ,
Res f ( z )e
lower
lower
i z
lower-half plane.
0- 127
Chapter 0
i z
Res f ( z )e ,
f ( x) cos xdx = 2 Im upper
i z
Res f ( z )e ,
f ( x) sin xdx = 2 Re upper
f ( x) cos x dx = 2 Im Res f ( z )e i z ,
lower
i z
f ( x) sin x dx = 2 Re Res f ( z )e .
lower
Proof:
Let S be a large
semi-circle with a
radius R. Consider
the
contour
consisting of the
semicircle S and
the diameter of S.
PHY 4602 Methods in Theoretical Physics II
0- 128
Chapter 0
f ( z )e dz = f ( z )e dz +
i z
i z
f ( x)ei z dz
= 2 i Re s f ( z )ei z
On the large semi-circle S ( z = R > R0 ) , we have:
f ( z ) M ( R) ,
where M ( R ) is the maximum value of f ( z ) on S and
lim M ( R ) = 0 . Besides, for / 2 0 we have:
/ 2 0 . Hence,
f ( z )ei z dz
f ( z )ei z dz
M ( R) e
R sin
2M ( R)
/2
/2
Rd = 2 M ( R )
/2
e R sin Rd
exp( 2 R / )Rd
exp( 2 R / ) Rd =
[1 exp( R)]
2
0- 129
Chapter 0
[1 exp( R)]
f ( z )e dz M ( R )
i z
Since lim M ( R ) = 0 ,
R
i z
,
f
(
x
)
cos
xdx
2
Im
Res
f
(
z
)
e
=
upper
i z
Res f ( z )e ,
f ( x) sin xdx = 2 Re upper
i z
f ( x) cos x dx = 2 Im Res f ( z )e ,
lower
f ( x) sin x dx = 2 Re Res f ( z )e i z .
lower
0- 130
Chapter 0
Example:
Consider a case with f ( x) = (k 2 + x 2 ) 1 .
Assume for the moment the parameter s > 0 .
eisz
eisx
s 2
2
k 2 + x 2 dx = 2 i Re
z = ik
+
k
z
eisz
= 2 i
= e ks .
2 z z =ik k
cos sx
ks
=
dx
e ,
k 2 + x 2
k
sin sx
k 2 + x 2 dx = 0
cos s x
k 2 + x2
dx =
k s
sin s x
cos sx
dx ,
k 2 + x 2
sin sx
k 2 + x 2 dx = 0 = k 2 + x 2 dx
sin sx
k 2 + x 2 dx = 0
0- 131
Chapter 0
B
A
f ( x)dx
B
A
f ( x)dx = lim
0
a
A
f ( x)dx + lim
0 a +
f ( x)dx
lim f ( x)dx + f ( x) dx P f ( x) dx
a +
A
0
A
Example:
The improper integral
dx
1 dx
dx
1
1
+ lim 3 = lim 2
1 x3 = lim
0 1 x 3
0 x
2 0 x
1
1
1
+ lim 2
0 x
0- 132
Chapter 0
Theorem
If f ( z ) has a simple pole at z = a on the real axis, then
lim f ( z )dz = i Res f ( z ) ,
r 0
z =a
C2
Proof:
From the Laurent expansion we have:
f ( z) =
b1
+ g ( z) ,
za
C2
f ( z )dz =
b1
i
ire
d + g ( z )dz.
i
C2
re
C2
z =a
0- 133
Chapter 0
Application:
If f ( z ) has several simple poles on real axis and the
integral
along
an
infinitely
large
semicircle
R S
Proof:
lim f ( z )dz = 0
R S
lim f ( z )dz
r 0
C2
= i Res f ( z )
z =a
ar
a+r
r 0
0- 134
Chapter 0
f ( z )dz + f ( z ) dz +
a+r
C2
f ( z ) dz + f ( z )dz
S
= 2 i Res f ( z )
upper
f ( x)dx = 2 i Res f ( z ) + i
UHP
Res f ( z ).
x axis
Examples:
dx
.
( x 1)( x 2)( x 2 + 1)
Evaluate P
z = 1,
1
1
Res f ( z ) =
= 2,
2
z =1
+
(
2)(
1)
z
z
z =1
z = 2,
1
1
=
Res f ( z ) =
,
2
z =2
( z 1)( z + 1) z =2 5
z = i,
1
Res f ( z ) =
z =i
(
1)(
2)(
)
z
z
z
i
z =i
1
3i
,
=
=
6 + 2i
20
dx
=
.
( x 1)( x 2)( x 2 + 1)
10
0- 135
Chapter 0
sin x
dx.
I =
Note:
The integral
sin z
dz does not vanish as R and
z
eiz
dz does vanish as R . Therefore,
z
eiz
eix
P
dx = i Res
x
z =0
z
= i.
cos x + i sin x
eix
P
dx = P
dx = i.
x
cos x
dx = 0,
sin x
dx = .
0- 136
Chapter 0
dx f ( x) ,
0 < 2
C1
0- 137
Chapter 0
dx f ( x) =
1
dz f ( z )Ln z.
C
2 i 1
1
dz f ( z )Ln z
+
C
C
1
2
2 i
= Res [ f ( z )Ln z ]
dx f ( x) =
dx f ( x) = Res f ( zk ) Ln zk .
k
dx
Ln zk
1
=
n
n 1
1+ x
k nz k
zk = n th roots of ( -1) .
0- 138
Chapter 0
dx
Ln zk
1
=
k nz n 1
1 + x3
k
i
5 i / 3
i /3
= 2 i / 3 + 2 i + 10 i / 3
3e
3e
3e
i 1 1
3
5 1
3
= i
+1+ + i
3 3 2
2
3 2
2
2 3
=
.
9
Evaluation of integrals with branch cuts
Consider integrals of the form
x 1 R ( x)dx ,
0- 139
Chapter 0
0 < 2 ,
z 1 = r 1 = x 1
if 0+ .
z 1 R ( z )dz
= e 2 i ( 1) x 1 R ( x)dx
0
+ x 1 R ( x)dx
0
= 2 i
inside C
1 e
2 i ( 1)
Res z 1 R( z )
x
0
2i sin 1
R ( x) dx =
0 x R ( x)dx.
i
e
0- 140
Chapter 0
Therefore,
e i
R ( x)dx =
sin
inside C
Res z 1 R ( z )
(1) 1
Res z 1 R ( z ) .
=
sin inside C
Example:
Evaluate the integral
I =
x 1
dx
1+ x
z -1
x 1
(1) 1
dx =
Res
1+ x
sin
1 + z z =1
.
sin
0- 141
Chapter 0
Transformation of integrals
An integral can be transformed to a more convenient form
by using Cauchys theorem to deform the contour. As an
example, we consider the integrals
F1 = dx cos ax 2 ,
0
F2 = dx sinax 2 ,
0
F = F1 + iF2
= dx e
iax 2
F = dz e
iaz 2
C1
= dz e
iaz 2
C2
+ dz e
C3
iaz 2
angular
1/ R 2
range
near = 0 .
Since
iaz 2
i / 4
=e
dr e
ar 2
1+ i 1
,
2 2 a
dx cos ax = dx sin ax 2 = / 8a .
2
0- 142
Chapter 0
m= p
f ( m) =
1
2 i
cot z f ( z )dz
Res [ cot( z ) f ( z )].
poles of f ( z )
inside C
cos z f ( z )
(sin z )
= f ( m)
z =m
cot z f ( z )dz
= 2 i
Res [ cot z f ( z )]
poles of f ( z )
inside C
0- 143
Chapter 0
2x
.
2
2 2
x
n
+
n =1
Solution:
(a) We consider the value of a related series:
N
m= N
2x
x 2 + m 2 2
m= N
2x
1
=
x 2 + m 2 2 2 i
cot zf ( z )dz
Res [ cot( z ) f ( z )],
C
poles of f ( z )
inside C
0- 144
Chapter 0
cot zf ( z )dz
1
2
cot z
2x
dz .
2
2 2
x +z
cot zf ( z )dz 0.
2x
2 x( cot z )
Res
=
2
2 2
x 2 + z 2 2
m = x + m
z = ix /
2 x cot (ix / ) cot (ix / )
=
+
2ix /
2ix /
= 2i cot ix = 2 coth x.
2x
2
+
= 2coth x
2
2 2
x
m =1 x + m
1
2x
coth x = 2
,
2 2
x m=1 x + m
0- 145
Chapter 0
cos z
1/ 2
cos 2 x + sinh 2 y
= 2
cot z =
sinh y
cot z =
2
1 + sinh y
2
= tanh y 1.
cos x + sinh y
lim cot z = lim 2
y
y sin x + sinh 2 y
= 1.
1)
f ( n)
0- 146
Chapter 0
Hence,
1
has first-order poles at
sin z
z = 0, 1, 2, n, .
( 1)
1
1
=
=
=
.
(sin z ) z = n cos n
n
Theorem
Let f ( z ) be a meromorphic function and let C be a
contour that encloses the zeros of sin z , located at
z = p, p + 1, , n . If we assume that the poles of f ( z ) and
(1)m f (m) =
m= p
1
2 i
C sin z f ( z )dz
f
z
Res
(
)
sin z
.
poles of f ( z )
inside C
0- 147
Chapter 0
dz
Re s f ( zk )
f ( z)
2 i
sin z
sin z k
zk
= 2 i
( 1)
f (n).
( sin z )
0 , we have:
( 1)
n =
f ( n) =
zk
Re s f ( zk )
sin zk
0- 148
Chapter 0
Example:
( 1)
n =
In this case, f ( z ) =
1
.
n
1
and there is only one simple pole
z
located at zk = . Besides,
Res f ( zk ) = 1 .
From the above theorem, we have:
1
=
, or
( 1)
n
sin
n =
1
1 2
=
sin
( 1)
n
n =1
0- 149
Chapter 0
The surfaces
= constant
are called the equipotential surfaces. At each point P the
gradient of is perpendicular to the surface = constant
through P ; that is, the electrical force has the direction
perpendicular to the equipotential surface.
If the potential is independent of one of the space
coordinates, say z , then the Laplaces equation becomes
2 2
= 2 + 2 = 0.
x
y
2
0- 151
Chapter 0
Examples:
Potential between parallel plates
( x)
=
1
1
( 2 1 ) x + ( 2 + 1 ) .
2
2
on
the
plates,
respectively located at x = 1.
r 2 rr + r r + = 0
Since = 0 , we then have:
r + = 0.
0- 152
Chapter 0
1
= ,
r
d ln
1
= ,
dr
r
ln = ln r + C ,
a
,
r
= a ln r + b.
r 2 rr + r r + = 0
We will make use of the theory of harmonic
functions to find the potential.
PHY 4602 Methods in Theoretical Physics II
0- 153
Chapter 0
1
1
+
( 2 1 ) ( 2 1 ) ,
where
= arctan
y
.
x
0- 154
Chapter 0
Complex Potential
Let ( x, y ) be harmonic in some domain D and ( x, y )
a conjugate harmonic of in D . The complex function
F ( z ) = ( x , y ) + i ( x, y )
is an analytic function of z = x + iy . This function F is
called the complex potential corresponding to the real
potential .
Theorem:
The curves = constant intersect the equipotential lines
Proof:
Assume that the curves = A and = B intersect at a
point z0 = x0 + iy0 . The gradient vectors of these two
functions are respectively given by:
= i
+j
x
y
and
= i
.
+j
x
y
0- 155
Chapter 0
i
i
j
i = i
+j
+
x
y
+
=
x x y y
From Cauchy-Riemann relations
=
=
and
,
x
y
y
x
i =
The
gradient
=0 .
+
y y y y
vectors
and
hence
the
curves
are
Note:
The gradient vectors vanish if F ( z ) = 0 and the proof
sketched above becomes invalid.
With the complex potential F = + i we can
simultaneously
handle
equipotential
lines
0- 156
Chapter 0
Examples:
Consider the complex potential
F ( z ) = az + b = ax + b + iay.
The corresponding real potential is ( x, y ) = ax + b ,
which gives the potential in between two parallel
plates. The lines of force are given by:
= ay = constant
= a ln r + b = a ln z + b .
The family of lines of force are given by:
= aArg z = constant
= constant .
0- 157
Chapter 0
is harmonic in D.
0- 158
Chapter 0
Proof:
Let * (u , v) be a harmonic conjugate of * (u , v) . The
complex function
F * ( w) * (u , v) + i * (u , v)
is then an analytic function of w. Since the composite of
analytic functions is also analytic, the composite function
F ( z ) F * ( f ( z ))
is an analytic function in D. Consequently, the real part of
F ( z ) , given by
Re F ( z ) ( x , y )
Re F * ( f ( z )) = * (u ( x, y ), v( x, y ))
is a harmonic function.
0- 159
Chapter 0
w=
az + b
cz + d
where
a b
.
c d
Note:
If
a b
= = k (i.e. ad bc = 0 ), the mapping becomes
c d
w=
kcz + kd
=k ,
cz + d
Identity
w= z
Translation
w = z+b
Rotation
w = az with a = 1
Linear transformation
w = az + b
0- 160
Chapter 0
w = 1/ z
Fixed points:
Fixed points of a mapping are points that are mapped onto
themselves. That is,
w = f ( z) = z .
For linear fractional transformation w =
az + b
, the fixed
cz + d
z=
az + b
.
cz + d
cz 2 ( a d ) z b = 0 .
This equation has at most two solutions except when
c = 0 = b and a = d .
PHY 4602 Methods in Theoretical Physics II
0- 161
Chapter 0
Theorem
A linear fractional transformation, not the identity, has at
most two fixed points. If a linear fractional transformation
is known to have three or more fixed points, it must be the
identity mapping w = z such that c = 0 = b and a = d .
f ( zi ) = g ( zi ) for i = 1, 2,3 ,
then these two mappings are identical.
Theorem
Three given distinct points z1 , z2 , z3 can always be mapped
onto three prescribed distinct points w1 , w2 , w3 by one, and
only one, linear fractional transformation w = f ( z ) . This
mapping is given implicitly by the equation
w w1 w2 w3 z z1 z2 z3
.
w w3 w2 w1 z z3 z2 z1
PHY 4602 Methods in Theoretical Physics II
0- 162
Chapter 0
Proof:
(a) It is obvious that the mapping
w w1 w2 w3 z z1 z2 z3
w w3 w2 w1 z z3 z2 z1
is a linear fractional mapping from z-plane to the wplane.
(b) It can be easily shown that
F ( w1 ) = 0 = G ( z1 ) ,
F ( w2 ) = 1 = G ( z2 ) ,
F ( w3 ) = = G ( z3 ) ,
where
F ( w) =
w w1 w2 w3
, and
w w3 w2 w1
G( z) =
z z1 z2 z3
.
z z3 z2 z1
i = 1, 2,3 .
(d) By the uniqueness theorem there is only one linear
fractional mapping that can map three distinct zi onto
wi and hence the theorem is proved.
0- 163
Chapter 0
Theorem
A bilinear transformation always maps circles into
circles.
Note: Here circles include straight lines as a limiting
case.
0- 164
Chapter 0
Standard mappings:
Mapping of a half-plane onto a disk
If we let z1 = 1 , z2 = 0 , z3 = 1 ; w1 = 1 , w2 = i ,
w3 = 1, then from
w w1 w2 w3 z z1 z2 z3
w w3 w2 w1 z z3 z2 z1
we have:
w ( 1) i 1
z ( 1) 0 1
w 1 i ( 1)
z 1 0 ( 1)
z i
w=
.
iz + 1
Mapping
z i
w=
iz + 1
0- 165
Chapter 0
respectively.
w w1 w2 w3 z z1 z2 z3
,
w w3 w2 w1 z z3 z2 z1
Since
we have:
w (0) i z ( 1) i 1
w i (0)
z 1 i (1)
w=
z +1
.
z 1
z
w
= 1,
= 0 = , z = , etc
0- 166
Chapter 0
If c =
1
2z 1
, then w =
.
2
z2
Examples:
Potential between noncoaxial cylinders
Find the potential between the cylinders C1 : z = 1
(potential U1 = 0 ) and C2 : z 2 / 5 = 2 / 5 (potential
U 2 = 110 volts).
PHY 4602 Methods in Theoretical Physics II
0- 167
Chapter 0
z b
, (b is a real number)
bz 1
z = 4 / 5 w = r0 .
r0 =
0b
= b,
0 1
r0 =
4 / 5 r0
4/5 b
.
=
4b / 5 1 4 r0 / 5 1
and b = 1/ 2 .
The mapping is given by
w = f ( z) =
2z 1
.
z2
0- 168
Chapter 0
0- 169
Chapter 0
We map the unit disk onto the right half of the w-plane by
using the linear fractional transformation
w = f ( z) =
1+ z
.
1 z
0- 170
Chapter 0
6000
, where = Arg w .
F * ( w) =
6000i
Ln w .
6000i
Ln
1+ z
,
1 z
and
( x, y ) = Re F ( z )
1+ z
=
.
Im Ln
1 z
6000
1+ z
=
Arg
1 z
6000
0- 171
Chapter 0
= 0.
x
y
dxdy = 0 ,
y
x
C (V1dx + V2 dy ) =
D
( x, y )
( a ,b )
(V1dx + V2 dy )
is independent of path in D .
0- 172
Chapter 0
dx +
dy ,
x
y
we have:
V1 = / x ,
V2 = / y .
0- 173
Chapter 0
,
V1 =
x
V2 =
.
y
F ( z ) = x i x = x + i y = V1 + iV2 = V .
0- 174
Chapter 0
Applications:
Flow around a 90 corner
The complex potential
F ( z ) = z 2 = x 2 y 2 + 2ixy
describes an incompressible irrotational flow around a 90
corner.
Equipotential lines = x 2 y 2 = const (Hyperbolas)
Streamlines
= 2 xy = const
Velocity
V = ( z 2 ) = 2 z = 2( x iy ) ,
(Hyperbolas)
V1 = 2 x , V2 = 2 y .
0- 175
Chapter 0
1
z
( x, y ) = r sin = const .
r
0- 176