Chapter 1. Vector Analysis

1.1. Introduction

This chapter will give an introduction to vector analysis. A vector, or displacement, has a direction and a magnitude. The following vector notations will be used in this course:

1. : the vector A

2. : the magnitude of vector A

The opposite of a vector is a vector with the same magnitude as but pointing in a direction opposite to that of . The opposite of vector is written as .
There are four different vector operations that we will be using in this course: vector addition, vector multiplication by a scalar, the dot product of two vectors, and the cross product of two vectors. We will start Chapter 1 by discussing these four vector operations in some detail.

1.1.1. Vector Addition

Two vectors and can be added. The result of this operation is a new vector (see Figure 1.1a). Using vector notation we can write vector addition as follows:



Vector addition is commutative. This means that the order in which two vectors are added does not affect the result (see Figure 1.1). The commutative properties of vector addition can be written as



To subtract vector from vector is equivalent to adding the opposite of vector to . In other words:




Figure 1.1. Vector Addition.

1.1.2. Multiplication by a scalar

A vector can be multiplied by a scalar. If the scalar a is a positive number (see Figure 1.2a) than the result of the multiplication of the vector by a is a new vector with a magnitude equal to and a direction equal to the direction of . If the scalar a is a negative number (see Figure 1.2b) than the result of the multiplication of the vector by a is a new vector with a magnitude equal to and a direction opposite to the direction of .

Figure 1.2. Vector multiplication.
Scalar multiplication is distributive. This means that the result of the multiplication of the vector sum of vector and by a scalar a is equal to the vector sum of and :



1.1.3. Dot product of two vectors


Figure 1.3. The scalar product of two vectors.
The dot product, also called the scalar product, is defined as



where is the angle between vectors and (see Figure 1.3). The dot product is commutative which means that the order of the vectors and does not effect the result of the dot product. In other words:



The dot product is also distributive:



The dot product will be frequently used to determine whether two vectors and are perpendicular. When and are perpendicular the angle is equal to 90°. The definition of the dot product shows that in this case the dot product between and is equal to zero.

1.1.4. Cross products of two vectors

The cross product of two vectors and is a third vector . The vector is perpendicular to both and , and has a length equal to



where is the smallest angle between and (see Figure 1.4). The direction of the vector can be determined using the right-hand rule. By applying the right-hand rule it can be shown easily that the vector product is not commutative. This implies that the order of the vectors and is important, and reversing the order will change the result of the vector product:



The vector product is distributive which requires that




Figure 1.4. The vector product of and .

Example: Problem 1.2
Is the cross product associative? If so, prove it; if not, provide a counter example.

If the cross product of two vectors is associative then the following relation must hold:



Consider the special case in which = and is perpendicular to and . In this case the cross product between and is equal to the null vector. The left-hand side of the equation is thus equal to the null vector. The cross product of and is a new vector, perpendicular to both and , and with a length equal to . The cross product between and is a new vector. Since and are perpendicular, the magnitude of the cross product of and is equal to which is non-zero. We thus conclude that is not equal to which shows that the cross product is not associative.

1.2. Vector components

A vector in three dimensions can be identified uniquely by specifying three coordinates. In a Cartesian coordinate frame three base vectors , , and are defined, parallel to the x, y, and z axes, respectively (see Figure 1.5). Each of these base vectors has unit length and is perpendicular to the other two base vectors. Any vector is uniquely defined by specifying its components along the x, y, and z axes. A vector is defined in terms of the three Cartesian coordinates Ax, Ay, and Az. Using the three base vectors, one can reconstruct the vector :




Figure 1.5. The Cartesian coordinate frame.
The vector operations discussed in Section 1.1 can be easily expressed in terms of vector components. We will now discuss each of these four vector operations.

1.2.1. Vector addition

The components of the vector sum of vectors and is equal to the sum of their like components. If a vector is the vector sum of the vectors and then the components of are equal to







In these equations, Ax, Ay, and Az are the vector components of vector , and Bx, By, and Bz are the vector components of vector .

1.2.2. Multiplication by a scalar

The components of a vector , which is the result of a scalar multiplication of vector with a scalar a, are equal to the components of multiplied by a:







1.2.3. Scalar product

The scalar product of vectors and is equal to the sum of the products of their like components:



This relation can be derived using the commutative properties of the scalar product:



In the last step of this derivation we have used the fact that the vectors , , and are perpendicular to each other and have unit length. The scalar products of these unit vectors are easy to evaluate and are equal to





1.2.4. Vector Product

The components of the vector , which is the vector product of vectors and , can be calculated using the distributive properties of the vector product:



In this derivation we have used the fact that the vectors , , and are perpendicular to each other and have unit length. The vector products of these unit vectors are easy to evaluate and are equal to









The expression of the vector product of and in terms of the components of these vectors can be neatly rewritten as a determinant:



Example: Problem 1.4
Use the vector product to find the components of the unit vector perpendicular to the plane shown in Figure 1.6.

Figure 1.6. Problem 1.4
The orientation of a plane in three dimensions is completely determined by specifying two vectors that are parallel to this plane. Two possible vectors are the vector , connecting (1,0,0) and (0,2,0), and the vector , connecting (1,0,0) and (0,0,3). The vector product of and is a vector which is perpendicular to both and . The vector will therefore be pointing in the same direction as . Using the definition of the vector product in terms of the components of vectors and we can calculate :



The length of the vector is equal to



The unit vector is parallel to and has a length equal to 1. Therefore



1.3. Vector Transformation

The three coordinates required to specify the direction and length of a vector are not uniquely defined. The same vector will in general have different coordinates in different coordinate systems. There is an infinite number of ways to choose a coordinate system, although the choice of the coordinate axes is usually influenced by the symmetry of the problem. Figure 1.7 illustrates two possible choices of the y and z axes. Coordinate system S is defined by the x, y, and z axes. Coordinate system S' is defined by the x', y', and z' axes. Coordinate system S' can be obtained by rotating coordinate system S through an angle around the x axis.

Figure 1.7. Coordinate Transformation.
The angle between the vector and the y axis is equal to . The y and z coordinates of the vector are therefore equal to





The angle between the vector and the y' axis is equal to . The y' and z' coordinates of the vector are therefore equal to





Using simple trigonometric relations we can obtain the component of the vector in S' in terms of the components of the vector in S:





These relations show, not unexpectedly, that the coordinates of the vector in coordinate system S are related to the coordinates of vector in coordinate system S'. This relation can be rewritten in matrix notation as



The rotation of the coordinate system S around the x axis leaves the x axis unchanged. The x coordinate of any vector in S will therefore be equal to the x' coordinate of this vector in S'. The three dimensional version of this coordinate transformation is therefore given by



The rotation of the coordinate system S around the x axis will not change the length of the vector . This can be verified by calculating the length of the vector in coordinate system S. This length is equal to



The vector in coordinate system S' is specified by the following coordinates:



The length of in coordinate system S' is defined in terms of its coordinates in S' as



Using the relation between the coordinates of in S' and the coordinates of in S we obtain



In general the coordinate transformation describing a rotation around an arbitrary axis can be written as



or, more compactly,



where A1 = Ax, A2 = Ay, and A3 = Az.

1.4. Differential calculus

The derivative df/dx of a function f(x) tells us how rapidly the function varies when the argument x is changed by a tiny amount dx:



A position dependent scalar function in three dimensions f(x, y, z) will in general be a function of three variables. The variation of f(x, y, z) between two closely spaced points depends not only on the distance between the two points, but also on their orientation. In this case, the theorem of partial derivatives can be used to calculate the change in f(x, y, z):



This equation can be rewritten as:



The first term in this expression is called the gradient of f, written as , and is defined as



The second term is called the infinitesimal displacement vector which is equal to



The change in the scalar function f can thus be written as



From the definition of the gradient it is clear that it is a vector. The direction of points in the direction of maximum increase of the function f. This follows immediately from the expression of df in terms of :



The change in f, df, will be maximum when = 0°. In this case, and are parallel. The magnitude of gives the slope of the function f along the direction of maximum change. If the gradient of f vanishes at a certain point then the function f has a maximum, a minimum, a saddle point, or a shoulder at that point.

Example: Problem 1.13
Let be the vector from some fixed point (x0, y0, z0) to the point (x, y, z), and let r be its length. Show that

a)

b)

c) What is he general formula for ?

a) The vector is equal to



The scalar function r2 is equal to the square of the length of the vector and will be a function of x, y, and z:



The gradient of r2 can be obtained using the following partial derivatives:







The gradient of r2 is equal to



where is the unit vector in the direction of .

b) The scalar function 1/r can be written as



The gradient of 1/r can be obtained using the following partial derivatives:







The gradient of 1/r is thus equal to



c) The most general form of can be guessed by comparing the results obtained in part a) and b). These answers suggest that the general form of is given by



The operator has the formal appearance of a vector and can be written as



The operator is also called the vector operator. If behaves like a vector than we expect that behaves like a vector, and consequently rotates like a vector.

Example: Problem 1.14
Suppose that f is a function of two variables (y and z) only. Show that the gradient transforms as a vector under a rotation about the x axis.

Consider a coordinate system S with x, y, and z axes. Coordinate system S is related to a coordinate system S' via a rotation by an angle about the x axis. The y and z coordinates of a vector in S are related to the y' and z' coordinates in S':



The gradient of f in the (y, z) frame is equal to



The gradient of f in the (y', z') frame is equal to



If transforms like a vector than the components of must be related to the components of in the following manner:





Using standard differential algebra we obtain the following relations for ∂f/∂y' and ∂f/∂z':





To evaluate ∂y/∂y', ∂y/∂z', ∂z/∂y', and ∂z/∂z' we must express y and z in terms of y' and z'. This can be achieved by manipulating the following two relations:





After a little algebra we obtain for y and z





Using these two relations we now can calculate the various partial derivatives of y and z:



Using these partial derivatives we can now express in terms of :





These last two equations clearly show that rotates like a vector.

If mimics a vector it may be used in vector multiplication. Three different types of vector multiplication can be carried out using :

1. scalar multiplication with scalar function f, , also called the gradient of f.

2. dot (scalar) product with vector function , , also called the divergence of .

3. cross (vector) product with vector function , , also called the curl of .

We will now discuss in detail the divergence and curl of .

1.4.1. The Divergence

The divergence of a vector is defined as



and is a scalar. The geometrical interpretation of the divergence of a vector function can be obtained by considering three special cases. First consider the vector function defined as



This vector function is shown at a few points in the x-y plane in Figure 1.8a. The divergence of is equal to



Now consider a vector function defined such that the vectors point toward the origin (see Figure 1.8b):



For this vector function the divergence is equal to




Figure 1.8. Various vector functions discussed in the text.
Finally consider the case in which the vector function is a constant vector of unit length along the y axis (see Figure 1.8c), independent of position:



For this vector function the divergence is equal to



We conclude that the divergence of a vector function evaluated at a particular point P is a measure of how much the vector function spreads out:

1. > 0: vector function is spreading out.

2. < 0: vector function is narrowing in.

3. = 0: vector function is not spreading out or narrowing in.

One vector field that we will be using frequently is the electric field . The three cases just discussed are relevant in electrostatics:

1. the divergence of generated by a positive point charge is positive (see Figure 1.9a).

2. the divergence of generated by a negative point charge is negative (see Figure 1.9b).

3. the divergence of generated by an infinitely large parallel-plate capacitor is zero (see Figure 1.9c).

Figure 1.9. a) Electric field generated by a positive point charge. b) Electric field generated by a negative point charge. c) Electric field generated by an infinitely large parallel-plate capacitor.

Example: Problem 1.17
In two dimensions, show that the divergence transforms as a scalar under rotations. Hint: use the method of Problem 1.14 to calculate the derivatives.

If the divergence of a vector function transforms as a scalar under rotation then the divergence of in a coordinate system S(x, y, z) must be identical to the divergence of in a coordinate system S'(x', y', z'). This requires that



Consider a rotation about the x axis. In this case, x = x' and vx = vx', and consequently



The y' and z' components of vector are related to the y and z components of vector in the following manner:





The partial derivative of vy' with respect to y' is equal to



Using the expressions for ∂y/∂y' and ∂z/∂y' obtained in problem 1.14 we can rewrite this expression as



In a similar manner the partial derivative of vz' with respect to z' can be obtained:



Combining these last two equations, we obtain



and therefore, the divergence of a vector function translates like a scalar under rotation.

1.4.2. The curl

The curl of vector function is equal to



The curl of a vector function evaluated at a certain point P is a measure of how much the vector function curls around this point. For the vector functions shown in Figure 1.8 the curl is zero.

Example: The magnetic field
Consider a wire of radius R carrying a current I along the z axis. The magnetic field produced by this wire depends only on the distance to the center of the wire. Its strength is equal to



and its direction is tangent to the circle centered on the wire (see Figure 1.10).

Figure 1.10. Magnetic vector field.
The magnetic field is a vector field, and the vector function is equal to



In order to calculate the curl of at r > R we have to evaluate ∂By/∂x and ∂Bx/∂y at r > R:





The curl of at r > R is thus equal to



In order to calculate the curl of at r < R we have to evaluate ∂By/∂x and ∂Bx/∂y at r < R:





The curl of at r < R is thus equal to



These calculations suggest that the curl of evaluated at a certain point is proportional to the current density at that point: outside the wire (r > R) the current density is 0, while inside the wire (r < R) the current density is I/(πR2).

1.4.3. Product Rules

Without proof I mention the following product rules involving that will be used frequently in this course:

1. : f and g are scalar functions

2. : and are vector functions

3. : f is a scalar function and is a vector function

4. : and are vector functions

5. : f is a scalar function and is a vector function

6. : and are vector functions

1.4.4. Second derivatives

By applying twice we can construct five species of second derivatives:

1. The divergence of a gradient:



where T is a scalar function of x, y, and z. is also called the Laplacian of T.

2. The curl of a gradient (is always zero):



This can be shown easily:



3. The gradient of a divergence:



4. The divergence of a curl (is always zero):



This can be shown easily:



5. The curl of a curl of a vector function can be expressed in terms of the Laplacian and the gradient of the divergence of :



1.5. Integral calculus

Consider a one-dimensional function df(x)/dx to be integrated between x = a and x = b. This integral can be obtained using the fundamental theorem of calculus, which states that



In three dimensions this fundamental theorem is replaced by the fundamental theorem for gradients, which states that



In a one-dimensional world there is just a single path from a to b. In a three-dimensional world there are an infinite number of ways to move from a to b. Nevertheless, the line integral of depends only on the function values at the end points and not on the path chosen. Thus we conclude:

1. The line integral is independent of the path chosen between a and b.

2. The line integral around any closed loop is zero.

Example 1.6
Although the line integral of the gradient of a vector function is independent of the path, the same is not true for arbitrary vector functions. Calculate

from a = (1, 1, 0) to b = (2, 2, 0)

for the vector function . Do it first by path (1), shown in Figure 1.11, and then by path (2).

Figure 1.11. Example 1.6
Consider first path (1). The first part of this path is parallel to the x axis (y = 1) and along this segment the scalar product between and is equal to



The line integral along this segment is equal to



The second part of path (1) is parallel to the y axis (x = 2) and along this segment the scalar product between and is equal to



The line integral along this segment is equal to



The line integral of along path (1) is equal to the sum of the line integrals along the two segments, and is thus equal to 1 + 10 = 11.
Consider now path (2). Along this path the x and y coordinates are equal (x = y and dx = dy). The displacement vector along this path is equal to



The scalar product between and along path (2) is equal to



The line integral along segment (2) is equal to



Comparing the line integral of for path (1) and the line integral of for path (2), we conclude that for the vector function the line integral is path dependent.

In a three-dimensional world we can have, besides line integrals, surface integrals and volume integrals. Various fundamental theorems can simplify the calculation of surface and volume integrals. For example, the fundamental theorem of divergences states that the integral of a divergence of a vector function over a volume is equal to the surface integral of the vector function over the surface that bounds the volume:



The right-hand side of this equation is called the flux of through the surface that bounds the volume. The vector is a vector whose magnitude is equal to the area of an infinitesimal surface element and whose direction is perpendicular to the surface, pointing outwards. The left-hand side of this equation represents a source term (for example, integrated over a volume is proportional to the total charge present in that volume).

Example: Problem 1.32
Test the divergence theorem for the function . Take as your volume the cube shown in Figure 1.12, with sides of length 2.

The divergence of is equal to



The volume integral of over the volume of the cube is equal to




Figure 1.12. Problem 1.32.
To evaluate the surface integral of across the surface of the cube, we have to consider each of the six faces separately. Start with considering the face of the cube in the x-y plane (with z = 0). The vector of this surface is equal to . The scalar product of and is equal to



In the last step we have used the fact that z = 0 on this face. For the same reason the scalar product of and is equal to zero on the face in the x-z plane and on the face in the y-z plane. On the face of the cube parallel to the x-y plane and with z = 2 the scalar product of and is equal to



The surface integral of over this face is equal to



On the face the cube parallel to the x-z plane and with y = 2 the scalar product of and is equal to



The surface integral of over this face is equal to



On the face of the cube parallel to the y-z plane and with x = 2 the scalar product of and is equal to



The surface integral of over this face is equal to



The surface integral across the surface of the cube is equal to the sum of the surface integrals across each of the 6 faces of the cube, and is thus equal to 24 + 16 + 8 = 48. This is equal to the volume integral of over the volume of the cube.

The fundamental theorem for curls, also known as Stokes' theorem, states that



Here, is an infinitesimal surface element. It is a vector whose magnitude is equal to the area of the surface element and whose direction is perpendicular to the surface. The vector is tangential to the boundary of the surface. The orientation of the surface vector and the direction of integration of the boundary should be consistent with he right-hand rule. The following two corollaries follow from Stokes' theorem:

1. depends only on the boundary line, not on the particular surface used.

2. around any closed surface.

Example: Problem 1.33
Test Stokes' theorem for the vector function , using the triangular shaded area shown in Figure 1.13.

The curl of is equal to




Fig. 1.13. Problem 1.33.
The surface vector is perpendicular to the surface and is equal to . The scalar product between and is equal to



The surface integral of is equal to



To evaluate the line integral around the boundary of the surface we have to evaluate the line integral along each of the three sides of the triangle. The direction of evaluation of the integral must be consistent with the chosen direction of . The right-hand rule requires that the line integral is evaluated counter clockwise. First consider the segment between (0, 0, 0) and (0, 2, 0). Along this segment . The vector product between and is equal to



since z = 0 along this segment. The line integral along this segment is therefore equal to zero. Now consider the line segment between (0, 2, 0) and (0, 0, 2). Along this segment and . The vector product between and is equal to



since x = 0 along this segment. Along this segment z = 2 - y and thus



The line integral of along this segment is thus equal to



Note: the limits of the integral are chosen such that the line integral is evaluated in a counter clockwise direction. The third segment of the boundary to be considered connects (0, 0, 2) and (0, 0, 0). The vector along this segment is equal to . The vector product between and is equal to



since x = 0 along this segment. The line integral of is therefore equal to zero. The line integral along the boundary of the surface is equal to the sum of the line integral of along each of the three segments. Thus



which is equal to the surface integral of .

1.6. Curvilinear coordinates

The Cartesian coordinate system is a coordinate system that is often used in calculations involving systems with no apparent symmetry. To describe systems that have spherical or cylindrical symmetry it is often more convenient to use spherical coordinates or cylindrical coordinates, respectively. These two coordinate systems will be discussed in this section.

1.6.1. Spherical coordinates

Spherical coordinates are always used when the system under consideration has spherical symmetry. The location of a point P (see Figure 1.14) is completely determined by specifying the following three coordinates:

1. r : the distance between the origin and P.

2. : the polar angle which is the angle between the vector P and the z-axis (see Figure 1.14).

3. : the azimuthal angle which is the angle between the projection of the vector to P in the x-y plane and the x axis.

Figure 1.14. Spherical coordinates.
In general, any vector can be expressed in terms of these three coordinates:



In contrast to the unit vectors , , and in a Cartesian coordinate system, the unit vectors , , and in a spherical coordinate system are not constant; they change direction as P moves around.
Consider a point P (see Figure 1.14) which is defined by the three spherical coordinates (r, , and ). The corresponding Cartesian coordinates are







The unit vectors in the spherical coordinate system can be calculated as follows:

1. : The direction of this unit vector points from the origin to point P:



2. : A change d in the polar angle will cause a change in the direction of the vector , without changing its length:



The unit vector is defined as



3. : A change d in the polar angle will cause a change in the direction of the vector without changing its length:



The unit vector is defined as



The expressions for , , and clearly show that the direction of these unit vectors depends on the point being described. In contrast, in a Cartesian coordinate system the direction of the three unit vectors is fixed, and independent of the point being described.
An infinitesimal element of length in the direction is simply dr:



The length of an infinitesimal element of length in the direction (as a result of a change in the polar angle of d) is equal to



The length of an infinitesimal element of length in the direction (as a result of a change in the azimuthal angle d) is equal to



The most general infinitesimal displacement is thus equal to



and this expression is used to evaluate line integrals in spherical coordinates. The infinitesimal volume element , in spherical coordinates, is the product of , , and :



There is no general expression for an infinitesimal surface element since it depends entirely on the orientation of the surface. For example, consider the surface of a sphere of radius r. On this surface r is constant, and the orientation of the surface is perpendicular to . In this case is equal to



Most calculations in spherical coordinates involve the application of vector derivatives. I will therefore summarize here, without derivation, the vector derivatives in spherical coordinates:

1. the gradient of a scalar function T:



2. the divergence of a vector function :



3. the curl of a vector function :



4. the Laplacian of a scalar function T:



1.6.2. Cylindrical coordinates

Cylindrical coordinates are always used when the system under consideration has cylindrical symmetry. The z axis is defined such that it coincides with the center axis of the cylinder. The location of a point P (see Figure 1.15) is completely determined by specifying the following three coordinates:

1. r : the distance between P and the z axis (see Figure 1.15).

2. : the azimuthal angle which is the angle between the projection of the vector to P in the x-y plane and the x axis (see Figure 1.15).

3. z : the z coordinate of point P (see Figure 1.15).

In general, any vector can be expressed in terms of these three coordinates:



However, the unit vectors and are not constant; they change direction as P moves around.

Figure 1.15. Cylindrical coordinates.
Consider a point P (see Figure 1.15) which is defined by three spherical coordinates (r, , and z). The corresponding Cartesian coordinates are







An infinitesimal change in r, , or z will produce the following infinitesimal displacements:







The infinitesimal displacement vector, to be used when evaluating line integrals in spherical coordinates, is equal to



The infinitesimal volume element , to be used in volume integration, is equal to



The infinitesimal area element , to be used in surface integration, is not uniquely determined and depends on the orientation of the surface. For example, if we carry out a surface integration over the surface of a cylinder of radius r (fixed), the infinitesimal surface element to be used is equal to



Most calculations in cylindrical coordinates involve the application of vector derivatives. I will therefore summarize here, without derivation, the vector derivatives in cylindrical coordinates:

1. The gradient of a scalar function T:



2. The divergence of a vector function :



3. The curl of a vector function :



4. The Laplacian of a scalar function T:



Example: Problem 1.42
a) Find the divergence of the function



b) Test the divergence theorem for this function, using the quarter-cylinder (radius 2, height 5) shown in Figure 1.16.
c) Find the curl of .

Figure 1.16. Problem 1.42.
a) To determine the divergence of we have to calculate the following partial derivatives:







Using the definition of the divergence of in terms of cylindrical coordinates we find that



b) The volume integral of the divergence of is equal to



To calculate the surface integral of we have to consider 5 different surfaces:

1. The surface parallel to the x-y plane at z = 0. The infinitesimal surface element for this surface is equal to



(Note: the direction of is perpendicular to the surface and pointing outwards). The scalar product between and is equal to



since z = 0 on this surface. The surface integral of across this surface is therefore equal to zero.

2. The surface parallel to the x-y plane at z = 5. The infinitesimal surface element for this surface is equal to



The scalar product between and is equal to



since z = 5 on this surface. The surface integral of across this surface is equal to



3. The surface in the x-z plane with = 0. The infinitesimal surface element for this surface is equal to



The scalar product between and is equal to



since = 0 on this surface. The surface integral of across this surface is therefore equal to zero.

4. The surface in the y-z plane with = π/2. The infinitesimal surface element for this surface is equal to



The scalar product between and is equal to



since = π/2 on this surface. The surface integral of across this surface is therefore equal to zero.

5. The cylinder wall (r = 2). The infinitesimal surface element for this surface is equal to



The scalar product between and is equal to



The surface integral of across this surface is equal to



The surface integral of across the whole surface is equal to the sum of the surface integral of across each of these five surfaces:



which is equal to the volume integral of .

1.7. The Dirac delta function

Consider the vector function :



Consider the surface integral of across the surface of a sphere of radius R. The surface element vector for this surface is equal to



The scalar product of and is equal to



The surface integral of across the surface of this sphere is equal to



and is independent of R. Applying the divergence theorem we conclude that



for every sphere, centered on r = 0, and independent of the radius R. This suggests that the only contribution to the volume integral of comes from a single point at r = 0. The divergence of is zero at every point except at r = 0:



At r = 0, the divergence of is undefined since 0/0 is undefined.
The function is called the Dirac delta function and will be used frequently in this course. The Dirac delta function has the following properties:

for any volume that includes r = 0

for any volume that does not include r = 0

The Dirac delta function is thus defined as



In Problem 1.13 we showed that the gradient of 1/r is equal to . This relation can be used to rewrite the expression for the Dirac delta function as



In this section we will discuss the use of the one- and three-dimensional Dirac delta functions.

1.7.1. The one-dimensional Dirac delta function

The one-dimensional Dirac delta function is defined as



and its integral is equal to



Any integral of between x = a and x = b will be equal to 1 if a < 0 < b.
Every time we will be using the Dirac delta function it will be used as part of the integrand of an integral. For example:



Example: Problem 1.43
Evaluate the following integrals:

a)

b)

c)

d)

a) The Dirac delta function will be 0 for all x except x = 3. The limits of the integral include x = 3. Thus



b) The Dirac delta function will be 0 for all x except x = π. The limits of the integral include x = π. Thus



c) The Dirac delta function will be 0 for all x except x = - 1. The limits of the integral do not include x = - 1. Thus



d) The Dirac delta function will be 0 for all x except x = - 2. The limits of the integral include x = - 2. Thus



Note: Most mistakes made in the evaluation of integrals containing the Dirac delta function occur when the argument of the Dirac delta function has a form other than (x + a). This is best illustrated by looking at the following example:



1.7.2. The three-dimensional Dirac delta function

The three-dimensional Dirac delta function is zero everywhere except at = 0 (x = 0, y = 0, z = 0). Every time we will be using the three-dimensional Dirac delta function it will be used as part of the integrand of an integral. For example:



if the integration volume V includes = 0.

Example: Problem 1.47a
Evaluate the following integral:



The Dirac delta function is zero everywhere except at . The volume integral is thus equal to



1.8. Helmholtz Theorem

Consider being told that the divergence of a vector function is the scalar function D and that the curl of the vector function is the vector function . Is this sufficient information to determine the vector function uniquely?
The vector function satisfies the following relations:





The divergence of must be zero since the divergence of the curl of any vector function is zero (). Consider the following solution:



where





If is a correct solution than the divergence of must be equal to D:



Here we have used one of the properties of second derivatives, which states that for any vector function .
If is a correct solution than the curl of must be equal to :



Here we have used one of the properties of second derivatives which states that for any scalar function U. The Laplacian of is equal to



The divergence of is equal to



In this derivation we have used the product rule for divergences (Griffiths page 21), the divergence theorem and the fact that the divergence of is equal to zero. The surface integral of /(r - r') will be equal to zero if goes to zero faster than 1/r2 when r → ∞. If this is the case, the divergence of is equal to zero, and consequently



So far we have verified that we can obtain a vector function if its divergence and curl are given. However, the vector function found is not unique. Consider a vector function whose curl and divergence vanish. Since the curl is distributive we can express the curl of the new vector function in terms of the curl of and the curl of :



Since the divergence is distributive we can express the divergence of the new vector function in terms of the divergence of and the divergence of :



However, there is no function that has zero divergence and zero curl and goes to zero at infinity. Since it is expected that the electric and magnetic fields go to zero at large distances we will exclude this type of contributions to by requiring that goes to zero at large distances. With this requirement, is uniquely defined if its curl and divergence are given. This conclusion is know as the Helmholtz theorem:

If the divergence and the curl of a vector function are specified, and if they both go to zero faster than 1/r2 as r → ∞ and if goes to zero as r → ∞ then is given uniquely by



where





1.9. Scalar and vector potentials

Consider a vector function which is the gradient of a scalar function : . The scalar function is also called the scalar potential of the field . The vector field , defined by the scalar function , has the following properties:

1. for any closed loop.

2. is independent of path, for any given set of end points.

3. everywhere.

The vector field generated by the scalar function is called a curl-less field. In Example 1.6 we showed that for the vector function the line integral is path dependent. The curl of this vector function is equal to and consequently this vector function can not be written as the gradient of a scalar function U. Therefore, the line integral of this vector function is path dependent, as was demonstrated in Example 1.6.
Now consider a vector field which is the curl of a vector function : . The vector function is called the vector potential for the field . The vector field defined by the vector potential W has the following properties:

1. for any closed surface.

2. is independent of the surface chosen, and depends only on the boundary line.

3. everywhere.

The vector field generated by the vector potential is called a divergence-less field.
The scalar potential U and the vector potential are not unique. Any constant can be added to U without effecting its gradient, since the gradient of a constant is equal to zero. Any gradient of a scalar function can be added to without effecting its curl, since the curl of a gradient is equal to zero.