Chapter III
Ideal Linear Oscillators

135. The main problem of the theory of heat radiation is to determine the energy distribution in the normal spectrum of black radiation, or, what amounts to the same thing, to find the function which has been left undetermined in the general expression of Wien’s displacement law (119), the function which connects the entropy of a certain radiation with its energy. The purpose of this chapter is to develop some preliminary theorems leading to this solution. Now since, as we have seen in Sec. 48, the normal energy distribution in a diathermanous medium cannot be established unless the medium exchanges radiation with an emitting and absorbing substance, it will be necessary for the treatment of this problem to consider more closely the processes which cause the creation and the destruction of heat rays, that is, the processes of emission and absorption. In view of the complexity of these processes and the difficulty of acquiring knowledge of any definite details regarding them, it would indeed be quite hopeless to expect to gain any certain results in this way, if it were not possible to use as a reliable guide in this obscure region the law of Kirchhoff derived in Sec. 51. This law states that a vacuum completely enclosed by reflecting walls, in which any emitting and absorbing bodies are scattered in any arrangement whatever, assumes in the course of time the stationary state of black radiation, which is completely determined by one parameter only, namely, the temperature, and in particular does not depend on the number, the nature, and the arrangement of the material bodies present. Hence, for the investigation of the properties of the state of black radiation the nature of the bodies which are assumed to be in the vacuum is perfectly immaterial. In fact, it does not even matter whether such bodies really exist somewhere in nature, provided their existence and their properties are consistent with the laws of thermodynamics and electrodynamics. If, for any special arbitrary assumption regarding the nature and arrangement of emitting and absorbing systems, we can find a state of radiation in the surrounding vacuum which is distinguished by absolute stability, this state can be no other than that of black radiation.

Since, according to this law, we are free to choose any system whatever, we now select from all possible emitting and absorbing systems the simplest conceivable one, namely, one consisting of a large number N of similar stationary oscillators, each consisting of two poles, charged with equal quantities of electricity of opposite sign, which may move relatively to each other on a fixed straight line, the axis of the oscillator.

It is true that it would be more general and in closer accord with the conditions in nature to assume the vibrations to be those of an oscillator consisting of two poles, each of which has three degrees of freedom of motion instead of one, i.e., to assume the vibrations as taking place in space instead of in a straight line only. Nevertheless we may, according to the fundamental principle stated above, restrict ourselves from the beginning to the treatment of one single component, without fear of any essential loss of generality of the conclusions we have in view.

It might, however, be questioned as a matter of principle, whether it is really permissible to think of the centers of mass of the oscillators as stationary, since, according to the kinetic theory of gases, all material particles which are contained in substances of finite temperature and free to move possess a certain finite mean kinetic energy of translatory motion. This objection, however, may also be removed by the consideration that the velocity is not fixed by the kinetic energy alone. We need only think of an oscillator as being loaded, say at its positive pole, with a comparatively large inert mass, which is perfectly neutral electrodynamically, in order to decrease its velocity for a given kinetic energy below any preassigned value whatever. Of course this consideration remains valid also, if, as is now frequently done, all inertia is reduced to electrodynamic action. For this action is at any rate of a kind quite different from the one to be considered in the following, and hence cannot influence it.

Let the state of such an oscillator be completely determined by its moment f(t), that is, by the product of the electric charge of the pole situated on the positive side of the axis and the pole distance, and by the derivative of f with respect to the time or

df(t) dt = (t). (204)

Let the energy of the oscillator be of the following simple form:

U = 1 2Kf2 + 1 2L2, (205)

where K and L denote positive constants, which depend on the nature of the oscillator in some way that need not be discussed at this point.

If during its vibration an oscillator neither absorbed nor emitted any energy, its energy of vibration, U, would remain constant, and we would have:

dU = Kfdf + Ld = 0, (205a)

or, on account of (204),

Kf(t) + Lf̈(t) = 0. (206)

The general solution of this differential equation is found to be a purely periodical vibration:

f = C cos(2πνt θ) (207)

where C and θ denote the integration constants and ν the number of vibrations per unit time:

ν = 1 2πK L . (208)

136. If now the assumed system of oscillators is in a space traversed by heat rays, the energy of vibration, U, of an oscillator will not in general remain constant, but will be always changing by absorption and emission of energy. Without, for the present, considering in detail the laws to which these processes are subject, let us consider any one arbitrarily given thermodynamic state of the oscillators and calculate its entropy, irrespective of the surrounding field of radiation. In doing this we proceed entirely according to the principle advanced in the two preceding chapters, allowing, however, at every stage for the conditions caused by the peculiarities of the case in question.

The first question is: What determines the thermodynamic state of the system considered? For this purpose, according to Sec. 124, the numbers N1N2, N3, of the oscillators, which lie in the region elements 12, 3, of the “state space” must be given. The state space of an oscillator contains those coordinates which determine the microscopic state of an oscillator. In the case in question these are only two in number, namely, the moment f and the rate at which it varies, , or instead of the latter the quantity

ψ = L, (209)

which is of the dimensions of an impulse. The region element of the state plane is, according to the hypothesis of quanta (Sec. 126), the double integral

dfdψ = h. (210)

The quantity h is the same for all region elements. A priori, it might, however, depend also on the nature of the system considered, for example, on the frequency of the oscillators. The following simple consideration, however, leads to the assumption that h is a universal constant. We know from the generalized displacement law of Wien (equation (119)) that in the universal function, which gives the entropy radiation as dependent on the energy radiation, there must appear a universal constant of the dimension c3u ν3 and this is of the dimension of a quantity of action47 ( erg sec.). Now, according to (210), the quantity h has precisely this dimension, on which account we may denote it as “element of action” or “quantity element of action.” Hence, unless a second constant also enters, h cannot depend on any other physical quantities.

137. The principal difference, compared with the calculations for an ideal gas in the preceding chapter, lies in the fact that we do not now assume the distribution densities w1w2, w3 of the oscillators among the separate region elements to vary but little from region to region as was assumed in Sec. 129. Accordingly the w’s are not small, but finite proper fractions, and the summation over the region elements cannot be written as an integration.

In the first place, as regards the shape of the region elements, the fact that in the case of undisturbed vibrations of an oscillator the phase is always changing, whereas the amplitude remains constant, leads to the conclusion that, for the macroscopic state of the oscillators, the amplitudes only, not the phases, must be considered, or in other words the region elements in the fψ plane are bounded by the curves C =  const., that is, by ellipses, since from (207) and (209)

f C2 + ψ 2πνLC2 = 1. (211)

The semi-axes of such an ellipse are:

a = C andb = 2πνLC. (212)

Accordingly the region elements 12, 3, n are the concentric, similar, and similarly situated elliptic rings, which are determined by the increasing values of C:

0,C1,C2,C3,Cn1,Cn. (213)

The nth region element is that which is bounded by the ellipses C = Cn1 and C = Cn. The first region element is the full ellipse C1. All these rings have the same area h, which is found by subtracting the area of the full ellipse Cn1 from that of the full ellipse Cn; hence

h = (anbn an1bn1)π
or, according to (212),
h = (Cn2 C n12)2π2νL,
where n = 1, 2, 3,.

From the additional fact that C0 = 0, it follows that:

Cn2 = nh 2π2νL. (214)

Thus the semi-axes of the bounding ellipses are in the ratio of the square roots of the integral numbers.

138. The thermodynamic state of the system of oscillators is fixed by the fact that the values of the distribution densities w1w2, w3, of the oscillators among the separate region elements are given. Within a region element the distribution of the oscillators is according to the law of elemental chaos (Sec. 122), i.e., it is approximately uniform.

These data suffice for calculating the entropy S as well as the energy E of the system in the given state, the former quantity directly from (173), the latter by the aid of (205). It must be kept in mind in the calculation that, since the energy varies appreciably within a region element, the energy En of all those oscillators which lie in the nth region element is to be found by an integration. Then the whole energy E of the system is:

E = E1 + E2 + En + . (215)

En may be calculated with the help of the law that within every region element the oscillators are uniformly distributed. If the nth region element contains, all told, Nn oscillators, there are per unit area Nn h  oscillators and hence Nn h df dψ per element of area. Hence we have:

En = Nn h Udfdψ.
In performing the integration, instead of f and ψ we take C and φ, as new variables, and since according to (211),

f = C cos φψ = 2πνLC sin φ (216)

we get:

En = 2πνLNn h UCdCdφ
to be integrated with respect to φ from 0 to 2π and with respect to C from Cn1 to Cn. If we substitute from (205)(209) and (216)

U = 1 2KC2, (217)

we obtain by integration

En = π2 2 νLKNn h (Cn4 C n14)
and from (214) and (208):

En = Nn(n 1 2)hν = Nwn(n 1 2)hν, (218)

that is, the mean energy of an oscillator in the nth region element is (n 1 2)hν. This is exactly the arithmetic mean of the energies (n 1)hν and nhν which correspond to the two ellipses C = Cn1 and C = Cn bounding the region, as may be seen from (217), if the values of Cn1 and Cn are therein substituted from (214).

The total energy E is, according to (215),

E = Nhν n=1n=(n 1 2)wn. (219)

139. Let us now consider the state of thermodynamic equilibrium of the oscillators. According to the second principle of thermodynamics, the entropy S is in that case a maximum for a given energy E. Hence we assume E in (219) as given. Then from (179) we have for the state of equilibrium:

δS = 0 = 1(log w n + 1)δwn,
where according to (167) and (219)
1δw n = 0 and 1(n 1 2)δwn = 0.
From these relations we find:
log wn + βn +  const. = 0
or

wn = αγn. (220)

The values of the constants α and γ follow from equations (167) and (219):

α = 2Nhν 2E Nhνγ = 2E Nhν 2E + Nhν. (221)

Since wn is essentially positive it follows that equilibrium is not possible in the system of oscillators considered unless the total energy E has a greater value than Nhν 2 , that is unless the mean energy of the oscillators is at least hν 2 . This, according to (218), is the mean energy of the oscillators lying in the first region element. In fact, in this extreme case all N oscillators lie in the first region element, the region of smallest energy; within this element they are arranged uniformly.

The entropy S of the system, which is in thermodynamic equilibrium, is found by combining (173) with (220) and (221)

S = kN{ E Nhν + 1 2 log E Nhν + 1 2 E Nhν 1 2 log E Nhν 1 2 }.  (222)

140. The connection between energy and entropy just obtained allows furthermore a certain conclusion as regards the temperature. For from the equation of the second principle of thermodynamics, dS = dE T and from differentiation of (222) with respect to E it follows that

E = N hν 2 1 + ehν kT 1 ehν kT = Nhν 1 2 + 1 ehν kT 1 . (223)

Hence, for the zero point of the absolute temperature E becomes, not 0, but N hν 2 . This is the extreme case discussed in the preceding paragraph, which just allows thermodynamic equilibrium to exist. That the oscillators are said to perform vibrations even at the temperature zero, the mean energy of which is as large as hν 2 and hence may become quite large for rapid vibrations, may at first sight seem strange. It seems to me, however, that certain facts point to the existence, inside the atoms, of vibrations independent of the temperature and supplied with appreciable energy, which need only a small suitable excitation to become evident externally. For example, the velocity, sometimes very large, of secondary cathode rays produced by Roentgen rays, and that of electrons liberated by photoelectric effect are independent of the temperature of the metal and of the intensity of the exciting radiation. Moreover the radioactive energies are also independent of the temperature. It is also well known that the close connection between the inertia of matter and its energy as postulated by the relativity principle leads to the assumption of very appreciable quantities of intra-atomic energy even at the zero of absolute temperature.

For the extreme case, T = , we find from (223) that

E = NkT, (224)

i.e., the energy is proportional to the temperature and independent of the size of the quantum of action, h, and of the nature of the oscillators. It is of interest to compare this value of the energy of vibration E of the system of oscillators, which holds at high temperatures, with the kinetic energy L of the molecular motion of an ideal monatomic gas at the same temperature as calculated in (200). From the comparison it follows that

E = 2 3L. (225)

This simple relation is caused by the fact that for high temperatures the contents of the hypothesis of quanta coincide with those of the classical statistical mechanics. Then the absolute magnitude of the region element, G or h respectively, becomes physically unimportant (compare Sec. 125) and we have the simple law of equipartition of the energy among all variables in question (see below Sec. 169). The factor 2 3 in equation (225) is due to the fact that the kinetic energy of a moving molecule depends on three variables (ξηζ,) and the energy of a vibrating oscillator on only two (fψ).

The heat capacity of the system of oscillators in question is, from (223),

dE dT = Nk hν kT2 ehν kT (ehν kT 1)2. (226)

It vanishes for T = 0 and becomes equal to Nk for T = . A. Einstein48 has made an important application of this equation to the heat capacity of solid bodies, but a closer discussion of this would be beyond the scope of the investigations to be made in this book.

For the constants α and γ in the expression (220) for the distribution density w we find from (221):

α = ehν kT 1γ = ehν kT (227)

and finally for the entropy S of our system as a function of temperature:

S = kN hν kT ehν kT 1 log 1 ehν kT . (228)