A Treatise on Electricity and Magnetism/Part I/Chapter IV
CHAPTER IV.
GENERAL THEOREMS.
95.] In the preceding chapter we have calculated the potential function and investigated its properties on the hypothesis that there is a direct action at a distance between electrified bodies, which is the resultant of the direct actions between the various electrified parts of the bodies.
If we call this the direct method of investigation, the inverse method will consist in assuming that the potential is a function characterised by properties the same as those which we have already established, and investigating the form of the function.
In the direct method the potential is calculated from the distribution of electricity by a process of integration, and is found to satisfy certain partial differential equations. In the inverse method the partial differential equations are supposed given, and we have to find the potential and the distribution of electricity.
It is only in problems in which the distribution of electricity is given that the direct method can be used. When we have to find the distribution on a conductor we must make use of the inverse method.
We have now to shew that the inverse method leads in every case to a determinate result, and to establish certain general theorems deduced from Poisson's partial differential equation

The mathematical ideas expressed by this equation are of a different kind from those expressed by the equation
. 
In the differential equation we express that the values of the second derivatives of V in the neighbourhood of any point, and the density at that point are related to each other in a certain manner, and no relation is expressed between the value of at that point and the value of at any point at a sensible distance from it.
In the second expression, on the other hand, the distance between the point at which exists from the point at which exists is denoted by , and is distinctly recognised in the expression to be integrated.
The integral, therefore, is the appropriate mathematical expression for a theory of action between particles at a distance, whereas the differential equation is the appropriate expression for a theory of action exerted between contiguous parts of a medium.
We have seen that the result of the integration satisfies the differential equation. We have now to shew that it is the only solution of that equation fulfilling certain conditions.
We shall in this way not only establish the mathematical equi valence of the two expressions, but prepare our minds to pass from the theory of direct action at a distance to that of action between contiguous parts of a medium.
96.] The potential function , considered as derived by integration from a known distribution of electricity either in the substance of bodies with the volumedensity or on certain surfaces with the surfacedensity and being everywhere finite, has been shewn to have the following characteristics:—
(1) is finite and continuous throughout all space.
(2) vanishes at the infinite distance from the electrified system.
(3) The first derivatives of are finite throughout all space, and continuous except at the electrified surfaces.
(4) At every point of space, except on the electrified surfaces, the equation of Poisson

is satisfied. We shall refer to this equation as the General Characteristic equation.
At every point where there is no electrification this equation becomes the equation of Laplace,
. 
(5) At any point of an electrified surface at which the surfacedensity is , the first derivative of , taken with respect to the normal to the surface, changes its value abruptly at the surface, so that
, 
where and are the normals on either side of the surface, and and are the corresponding potentials. We shall refer to this equation as the Superficial Characteristic equation.
(6) If denote the potential at a point whose distance from any fixed point in a finite electrical system is , then the product , when increases indefinitely, is ultimately equal to , the total charge in the finite system.
97.] Lemma. Let V be any continuous function of and let be functions of subject to the general solenoidal condition
,  (1) 
where these functions are continuous, and to the superficial solenoidal condition
,  (2) 
at any surface at which these functions become discontinuous, being the directioncosines of the normal to the surface, and and the values of the functions on opposite sides of the surface, then the triple integral
 (3) 
vanishes when the integration is extended over a space bounded by surfaces at which either is constant, or
,  (4) 
being the directioncosines of the surface.
Before proceeding to prove this theorem analytically we may observe, that if be taken to represent the components of the velocity of a homogeneous incompressible fluid of density unity, and if be taken to represent the potential at any point of space of forces acting on the fluid, then the general and superficial equations of continuity ((1) and (2)) indicate that every part of the space is, and continues to be, full of the fluid, and equation (4) is the condition to be fulfilled at a surface through which the fluid does not pass.
The integral represents the work done by the fluid against the forces acting on it in unit of time.
Now, since the forces which act on the fluid are derived from the potential function , the work which they do is subject to the law of conservation of energy, and the work done on the whole fluid within a certain space may be found if we know the potential at the points where each line of flow enters the space and where it issues from it. The excess of the second of these potentials over the first, multiplied by the quantity of fluid which is transmitted along each line of flow, will give the work done by that portion of the fluid, and the sum of all such products will give the whole work.
Now, if the space be bounded by a surface for which , a constant quantity, the potential will be the same at the place where any line of flow enters the space and where it issues from it, so that in this case no work will be done by the forces on the fluid within the space, and .
Secondly, if the space be bounded in whole or in part by a surface satisfying equation (4), no fluid will enter or leave the space through this surface, so that no part of the value of can depend on this part of the surface.
The quantity is therefore zero for a space bounded externally by the closed surface , and it remains zero though any part of this space be cut off from the rest by surfaces fulfilling the condition (4).
The analytical expression of the process by which we deduce the work done in the interior of the space from that which takes place at the bounding surface is contained in the following method of integration by parts.
Taking the first term of the integral ,
, 
where 
+ &c.; 
and where , , &c. are the values of and at the points whose coordinates are , , &c., , , &c. being the values of where the ordinate cuts the bounding surface or surfaces, arranged in descending order of magnitude.
Adding the two other terms of the integral , we find
; 
the integration of the first term being extended over the bounding surface, and that of the second throughout the entire space.
For all spaces within which are continuous, the second term vanishes in virtue of equation (1). If for any surface within the space are discontinuous but subject to equation (2), we find for the part of depending on this surface,
, 
; 
where the suffixes and , applied to any symbol, indicate to which of the two spaces separated by the surface the symbol belongs.
Now, since is continuous, we have at every point of the surface,
we have also 
; 
but since the normals are drawn in opposite directions, we have
so that the total value of M, so far as it depends on the surface of discontinuity, is
. 
The quantity under the integral sign vanishes at every point in virtue of the superficial solenoidal condition or characteristic (2).
Hence, in determining the value of , we have only to consider the surfaceintegral over the actual bounding surface of the space considered, or
. 
Case 1. If V is constant over the whole surface and equal to ,
. 
The part of this expression under the sign of double integration represents the surfaceintegral of the flux whose components are and by Art. 21 this surfaceintegral is zero for the closed surface in virtue of the general and superficial solenoidal conditions (1) and (2).
Hence for a space bounded by a single equipotential surface.
If the space is bounded externally by the surface V = C, and internally by the surfaces , , &c., then the total value of for the space so bounded will be
where is the value of the integral for the whole space within the surface , and &c. are the values of the integral for the spaces within the internal surfaces. But we have seen that , &c. are each of them zero, so that the integral is zero also for the periphractic region between the surfaces.
Case 2. If is zero over any part of the bounding surface, that part of the surface can contribute nothing to the value of , because the quantity under the integral sign is everywhere zero. Hence will remain zero if a surface fulfilling this con dition is substituted for any part of the bounding surface, provided that the remainder of the surface is all at the same potential.
98.] We are now prepared to prove a theorem which we owe to Sir William Thomson ^{[1]}.
As we shall require this theorem in various parts of our subject, I shall put it in a form capable of the necessary modifications.
Let be any functions of (we may call them the components of a flux) subject only to the condition
 (5) 
where has given values within a certain space. This is the general characteristic of .
Let us also suppose that at certain surfaces (S) , and are discontinuous, but satisfy the condition
;  (6) 
where are the directioncosines of the normal to the surface, the values of on the positive side of the surface, and those on the negative side, and a quantity given for every point of the surface. This condition is the superficial characteristic of .
;  (7) 
and the superficial characteristic at the surfaces
,
 (8) 
being a quantity which may be positive or zero but not negative, given at every point of space.
Finally, let represent the triple integral
,  (9) 
extended over a space bounded by surfaces, for each of which either
or 
,  (10) 
where the value of is given at every point of the surface; then, if be supposed to vary in any manner, subject to the above conditions, the value of will be a unique minimum, when
.  (11) 
If we put for the general values of
;  (12) 
then, by substituting these values in equations (5) and (7), we find that satisfy the general solenoidal condition
We also find, by equations (6) and (8), that at the surfaces of discontinuity the values of and satisfy the superficial solenoidal condition
The quantities , therefore, satisfy at every point the solenoidal conditions as stated in the preceding lemma.
We may now express in terms of and ,
The last term of may be written , where is the quantity considered in the lemma, and which we proved to be zero when the space is bounded by surfaces, each of which is either equipotential or satisfies the condition of equation (10), which may be written
is therefore reduced to the sum of the first and second terms.
In each of these terms the quantity under the sign of integration consists of the sum of three squares, and is therefore essentially positive or zero. Hence the result of integration can only be positive or zero.
Let us suppose the function known, and let us find what values of will make a minimum.
If we assume that at every point , , and , these values fulfil the solenoidal conditions, and the second term of is zero, and is then a minimum as regards the variation of .
For if any of these quantities had at any point values differing from zero, the second term of would have a positive value, and would be greater than in the case which we have assumed.
But if , , and , then
Hence these values of make a minimum.
But the values of , as expressed in equations (12), are perfectly general, and include all values of these quantities con sistent with the conditions of the theorem. Hence, no other values of can make a minimum.
Again, is a quantity essentially positive, and therefore is always capable of a minimum value by the variation of . Hence the values of which make a minimum must have a real existence. It does not follow that our mathematical methods are sufficiently powerful to determine them.
Corollary I. If and are given at every point of space, and if we write with the condition (1)
, 
then can be found without ambiguity from these four equations.
Corollary II. The general characteristic equation

where a finite quantity of single value whose first derivatives are finite and continuous except at the surface , and at that surface fulfil the superficial characteristic

can be satisfied by one value of , and by one only, in the following cases.
Case 1. When the equations apply to the space within any closed surface at every point of which .
For we have proved that in this case have real and unique values which determine the first derivatives of , and hence, if different values of exist, they can only differ by a constant. But at the surface is given equal to , and therefore is determinate throughout the space.
As a particular case, let us suppose a space within which bounded by a closed surface at which . The characteristic equations are satisfied by making for every point within the space, and therefore is the only solution of the equations.
Case 2. When the equations apply to the space within any closed surface at every point of which is given.
For if in this case the characteristic equations could be satisfied by two different values of , say and , put , then subtracting the characteristic equation in from that in , we find a characteristic equation in . At the closed surface because at the surface , and within the surface the density is zero because . Hence, by Case 1, throughout the enclosed space, and therefore throughout this space. Case 3. When the equations apply to a space bounded by a closed surface consisting of two parts, in one of which is given at every point, and in the other
, 
where is given at every point.
For if there are two values of ,let represent, as before, their difference, then we shall have the equation fulfilled within a closed surface consisting of two parts, in one of which , and in the other
; 
and since satisfies the equation it is the only solution, and therefore there is but one value of possible.
Note.—The function in this theorem is restricted to one value at each point of space. If multiple values are admitted, then, if the space considered is a cyclic space, the equations may be satisfied by values of containing terms with multiple values. Examples of this will occur in Electromagnetism.
99.] To apply this theorem to determine the distribution of electricity in an electrified system, we must make throughout the space occupied by air, and throughout the space occupied by conductors. If any part of the space is occupied by dielectrics whose inductive capacity differs from that of air, we must make K in that part of the space equal to the specific inductive capacity.
The value of , determined so as to fulfil these conditions, will be the only possible value of the potential in the given system.
Green's Theorem shews that the quantity , when it has its minimum value corresponding to a given distribution of electricity, represents the potential energy of that distribution of electricity. See Art. 100, equation (11).
In the form in which is expressed as the result of integration over every part of the field, it indicates that the energy due to the electrification of the bodies in the field may be considered as the result of the summation of a certain quantity which exists in every part of the field where electrical force is in action, whether elec trification be present or not in that part of the field.
The mathematical method, therefore, in which , the symbol of electrical energy, is made an object of study, instead of , the symbol of electricity itself, corresponds to the method of physical speculation, in which we look for the seat of electrical action in every part of the field, instead of confining our attention to the electrified bodies.
The fact that attains a minimum value when the components of the electric force are expressed in terms of the first derivatives of a potential, shews that, if it were possible for the electric force to be distributed in any other manner, a mechanical force would be brought into play tending to bring the distribution of force into its actual state. The actual state of the electric field is therefore a state of stable equilibrium, considered with reference to all variations of that state consistent with the actual distribution of free electricity.
100.] The following remarkable theorem was given by George Green in his essay ‘On the Application of Mathematics to Electricity and Magnetism.’
I have made use of the coefficient , introduced by Thomson, to give greater generality to the statement, and we shall find as we proceed that the theorem may be modified so as to apply to the most general constitution of crystallized media.
We shall suppose that and are two functions of , which, with their first derivatives, are finite and continuous within the space bounded by the closed surface .
We shall also put for conciseness
,  (1) 
and 
,  (2) 
where is a real quantity, given for each point of space, which may be positive or zero but not negative. The quantities and correspond to volumedensities in the theory of potentials, but in this investigation they are to be considered simply as abbreviations for the functions of and to which they are here equated.
In the same way we may put
,  (3) 
and 
,  (4) 
where are the directioncosines of the normal drawn inwards from the surface . The quantities and correspond to superficial densities, but at present we must consider them as defined by the above equations.
Green's Theorem is obtained by integrating by parts the expression
 (5) 
throughout the space within the surface .
If we consider as a component of a force whose potential is , and as a component of a flux, the expression will give the work done by the force on the flux.
If we apply the method of integration by parts, we find

;  (6) 
or 
.  (7) 
In precisely the same manner by exchanging and , we should find
 (8) 
The statement of Green's Theorem is that these three expressions for are identical, or that

 (9) 
There are cases in which the resultant force at any point of a certain region fulfils the ordinary condition of having a potential, while the potential itself is a manyvalued function of the coordinates. For instance, if
we find , a manyvalued function of and , the values of forming an arithmetical series whose common difference is , and in order to define which of these is to be taken in any particular case we must make some restriction as to the line along which we are to integrate the force from the point where to the required point.
In this case the region in which the condition of having a potential is fulfilled is the cyclic region surrounding the axis of z, this axis being a line in which the forces are infinite and therefore not itself included in the region.
The part of the infinite plane of for which is positive may be taken as a diaphragm of this cyclic region. If we begin at a point close to the positive side of this diaphragm, and integrate along a line which is restricted from passing through the diaphragm, the lineintegral will be restricted to that value of which is positive but less than .
Let us now suppose that the region bounded by the closed surface in Green s Theorem is a cyclic region of any number of cycles, and that the function is a manyvalued function having any number of cyclic constants.
The quantities , , and will have definite values at all points within , so that the volumeintegral
has a definite value, and have also definite values, so that if is a single valued function, the expression
has also a definite value.
The expression involving has no definite value as it stands, for is a manyvalued function, and any expression containing it is manyvalued unless some rule be given whereby we are directed to select one of the many values of V at each point of the region.
To make the value of definite in a region of cycles, we must conceive diaphragms or surfaces, each of which completely shuts one of the channels of communication between the parts of the cyclic region. Each of these diaphragms reduces the number of cycles by unity, and when n of them are drawn the region is still a connected region but acyclic, so that we can pass from any one point to any other without cutting a surface, but only by reconcileable paths.
Let be the first of these diaphragms, and let the lineintegral of the force for a line drawn in the acyclic space from a point on the positive side of this surface to the contiguous point on the negative side be , then is the first cyclic constant.
Let the other diaphragms, and their corresponding cyclic constants, be distinguished by suffixes from 1 to n, then, since the region is rendered acyclic by these diaphragms, we may apply to it the theorem in its original form.
We thus obtain for the complete expression for the first member of the equation
.
The addition of these terms to the expression of Green's Theorem, in the case of manyvalued functions, was first shewn to be necessary by Helmholtz^{[2]}, and was first applied to the theorem by Thomson^{[3]}.
Physical Interpretation of Green’s Theorem.
The expressions and denote the quantities of electricity existing on an element of the surface S and in an element of volume respectively. We may therefore write for either of these quantities the symbol e, denoting a quantity of electricity. We shall then express Green's Theorem as follows—
;
where we have two systems of electrified bodies, e standing in succession for e_{1}, e_{2}, &c., any portions of the electrification of the first system, and V denoting the potential at any point due to all these portions, while e' stands in succession for e_{1}', e_{2}' , &c., portions of the second system, and V' denotes the potential at any point due to the second system.
Hence Ve' denotes the product of a quantity of electricity at a point belonging to the second system into the potential at that point due to the first system, and denotes the sum of all such quantities, or in other words, represents that part of the energy of the whole electrified system which is due to the action of the second system on the first.
In the same way represents that part of the energy of the whole system which is due to the action of the first system on the second.
If we define V as , where r is the distance of the quantity e of electricity from the given point, then the equality between these two values of M may be obtained as follows, without Green Theorem—
This mode of regarding the question belongs to what we have called the direct method, in which we begin by considering certain portions of electricity, placed at certain points of space, and acting on one another in a way depending on the distances between these points, no account being taken of any intervening medium, or of any action supposed to take place in the intervening space.
Green's Theorem, on the other hand, belongs essentially to what we have called the inverse method. The potential is not supposed to arise from the electrification by a process of summation, but the electrification is supposed to be deduced from a perfectly arbitrary function called the potential by a process of differentiation.
In the direct method, the equation is a simple extension of the law that when any force acts directly between two bodies, action and reaction are equal and opposite.
In the inverse method the two quantities are not proved directly to be equal, but each is proved equal to a third quantity, a triple integral which we must endeavour to interpret.
If we write R for the resultant electromotive force due to the potential V, and l, m, n for the directioncosines of R, then, by Art. 71,
If we also write R for the force due to the second system, and l, m, n for its directioncosines,
and the quantity M may be written
 (10) 
where
being the angle between the directions of R and .
Now if K is what we have called the coefficient of electric inductive capacity, then KR will be the electric displacement due to the electromotive force R, and the product will represent the work done by the force on account of the displacement caused by the force , or in other words, the amount of intrinsic energy in that part of the field due to the mutual action of and .
We therefore conclude that the physical interpretation of Green's theorem is as follows :
If the energy which is known to exist in an electrified system is due to actions which take place in all parts of the field, and not to direct action at a distance between the electrified bodies, then that part of the intrinsic energy of any part of the field upon which the mutual action of two electrified systems depends is per unit of volume.
The energy of an electrified system due to its action on itself is, by Art. 85,
which is by Green's theorem, putting U = V,
 (11) 
and this is the unique minimum value of the integral considered in Thomson's theorem.
101.] Let a closed surface S be maintained at potential zero. Let P and Q be two points on the positive side of the surface S (we may suppose either the inside or the outside positive), and let a small body charged with unit of electricity be placed at P; the potential at the point Q will consist of two parts, of which one is due to the direct action of the electricity on P, while the other is due to the action of the electricity induced on S by P. The latter part of the potential is called Green's Function, and is denoted by G_{pq}.
This quantity is a function of the positions of the two points P and Q, the form of which depends on that of the surface S. It has been determined in the case in which S is a sphere, and in a very few other cases. It denotes the potential at Q due to the electricity induced on S by unit of electricity at P.
The actual potential at any point due to the electricity at and on is
where denotes the distance between and .
At the surface and at all points on the negative side of , the potential is zero, therefore
(1) 
where the suffix indicates that a point on the surface is taken instead of .
Let denote the surfacedensity induced by at a point of the surface , then, since is the potential at due to the superficial distribution,
(2) 
where is an element of the surface at , and the integration is to be extended over the whole surface .
But if unit of electricity had been placed at , we should have had by equation (1),
(3) 
(4) 
where is the density induced by on an element at , and is the distance between and . Substituting this value of in the expression for , we find
(5) 
Since this expression is not altered by changing into and into ,we find that
(6) 
a result which we have already shewn to be necessary in Art. 88, but which we now see to be deducible from the mathematical process by which Green’s function may be calculated.
If we assume any distribution of electricity whatever, and place in the field a point charged with unit of electricity, and if the surface of potential zero completely separates the point from the assumed distribution, then if we take this surface for the surface , and the point for , Green’s function, for any point on the same side of the surface as , will be the potential of the assumed distribution on the other side of the surface. In this way we may construct any number of cases in which Green’s function can be found for a particular position of . To find the form of the function when the form of the surface is given and the position of is arbitrary, is a problem of far greater difficulty, though, as we have proved, it is mathematically possible.
Let us suppose the problem solved, and that the point is taken within the surface. Then for all external points the potential of the superficial distribution is equal and opposite to that of . The superficial distribution is therefore centrobaric^{[4]}, and its action on all external points is the same as that of a unit of negative electricity placed at .
102.] Let a region be completely bounded by a number of surfaces , &c., and let be a quantity, positive or zero but not negative, given at every point of this region. Let be a function subject to the conditions that its values at the surfaces , &c. are the constant quantities , &c., and that at the surface
(1) 
where is a normal to the surface . Then the integral
(2) 
taken over the whole region, has a unique minimum when satisfies the equation
(3) 
throughout the region, as well as the original conditions.
We have already shewn that a function exists which fulfils the conditions (1) and (3), and that it is determinate in value. We have next to shew that of all functions fulfilling the surfaceconditions it makes a minimum.
Let be the function which satisfies (1) and (3), and let
(4) 
be a function which satisfies (1).
It follows from this that at the surfaces , &c. .
The value of becomes
(5) 
Let us confine our attention to the last of these three groups of terms, merely observing that the other groups are essentially positive. By Green’s theorem
(6) 
the first integral of the second member being extended over the surface of the region and the second throughout the enclosed space. But on the surfaces &c. , so that these contribute nothing to the surfaceintegral.
Again, on the surface , , so that this surface contributes nothing to the integral. Hence the surfaceintegral is zero.
The quantity within brackets in the volumeintegral also disappears by equation (3), so that the volumeintegral is also zero. Hence is reduced to
(7) 
Both these quantities are essentially positive, and therefore the minimum value of is when
(8) 
or when is a constant. But at the surfaces , &c. . Hence everywhere, and gives the unique minimum value of .
The quantity in its minimum form can be expressed by means of Green’s theorem in terms of , &c., the potentials of , and , &c., the charges of these surfaces,
(9) 
or, making use of the coefficients of capacity and induction as defined in Article 87,
(10) 
The accurate determination of the coefficients is in general difficult, involving the solution of the general equation of statical electricity, but we make use of the theorem we have proved to determine a superior limit to the value of any of these coefficients.
To determine a superior limit to the coefficient of capacity , make , and , &c. each equal to zero, and then take any function which shall have the value 1 at , and the value 0 at the other surfaces.
From this trial value of calculate by direct integration, and let the value thus found be . We know that is not less than the absolute minimum value , which in this case is .
Hence
(11) 
If we happen to have chosen the right value of the function , then , but if the function we have chosen differs slightly from the true form, then, since is a minimum, will still be a close approximation to the true value.
We may also determine a superior limit to the coefficients of potential defined in Article 86 by means of the minimum value of the quantity in Article 98, expressed in terms of .
By Thomson’s theorem, if within a certain region bounded by the surfaces &c. the quantities are subject to the condition
(12) 
and if
(13) 
be given all over the surface, where are the directioncosines of the normal, then the integral
(14) 
is an absolute and unique minimum when
(15) 
When the minimum is attained is evidently the same quantity which we had before.
If therefore we can find any form for which satisfies the condition (12) and at the same time makes
&c.;  (16) 
and if be the value of calculated by (14) from these values of , then is not less than
(17) 
If we take the case in which one of the surfaces, say , surrounds the rest at an infinite distance, we have the ordinary case of conductors in an infinite region; and if we make , and for all the other surfaces, we have at infinity, and is not greater than .
In the very important case in which the electrical action is entirely between two conducting surfaces and , of which completely surrounds and is kept at potential zero, we have and .
Hence in this case we have
(18) 
and we had before
(19) 
so that we conclude that the true value of , the capacity of the internal conductor, lies between these values.
This method of finding superior and inferior limits to the values of these coefficients was suggested by a memoir 'On the Theory of Resonance,' by the Hon. J. W. Strutt, Phil. Trans., 1871. See Art. 308.
 ↑ Cambridge and Dublin Mathematical Journal, February, 1848.
 ↑ ‚Ueber Integrate der Hydrodynamischen Gleichungen welche den Wirbelbewegungen entsprechen‘, Crelle, 1858. Translated by Tait in Phil. Mag., 1867.
 ↑ ‚On Vortex Motion‘, Trans. R. S. Edin., xxv. part i. p. 241 (1868).
 ↑ Thomson and Tait’s Natural Philosophy, § 526.