Derivation of Bose’s Entropy Spectral Density from the Multiplicity of Energy Eigenvalues

The modern textbook analysis of the thermal state of photons inside a three-dimensional reflective cavity is based on the three quantum numbers that characterize photon’s energy eigenvalues coming out when the boundary conditions are imposed. The crucial passage from the quantum numbers to the continuous frequency is operated by introducing a three-dimensional continuous version of the three discrete quantum numbers, which leads to the energy spectral density and to the entropy spectral density. This standard analysis obscures the role of the multiplicity of energy eigenvalues associated to the same eigenfrequency. In this paper we review the past derivations of Bose’s entropy spectral density and present a new analysis of energy spectral density and entropy spectral density based on the multiplicity of energy eigenvalues. Our analysis explicitly defines the eigenfrequency distribution of energy and entropy and uses it as a starting point for the passage from the discrete eigenfrequencies to the continuous frequency.


Introduction
The most important step after the fundamental paper [1] of Planck on the energy spectral density of the blackbody radiation was moved by Bose in [2], whose result was immediately extended by Einstein in [3] to massive particles, leading to the Bose-Einstein distribution of energy.Entropy plays a key role both in Planck's derivation and in Bose's derivation of the energy spectral density of the radiation.Today, the entropy calculated by Bose is recognized to be the equal to the thermodynamic entropy of the blackbody radiation.
Bose's derivation is based on two ingredients.One is the geometric probability distribution of the number of energy quanta that occupy one quantum state, the other is the so-called density of states, that is, the infinitesimal number of quantum states in an infinitesimal interval of frequency.Bose's calculation of the density of states is based on the quantization of the phase space into cells of volume h 3 , where h is Planck's constant.This approach, where one cell is a quantum state, is today surpassed by the modern quantum approach, where the quantum states used for the calculation of the density of states are the quantum eigenstates obtained by imposing the boundary conditions determined by the spatial boundaries of the body.
In the case of massive particles, the relationship between multiplicity of energy eigenvalues and density of states is well known and widely discussed, see for instance [4] and 1.4 of [5].However, in the case of photons, today's textbook treatment completely skips the multiplicity of energy eigenvalues in the derivation of the density of states and, with it, it skips the discrete energy spectrum at the eigenfrequencies corresponding to the energy eigenvalues.
In this paper we fill this gap.First we derive the discrete spectrum of energy and entropy at system's eigenfrequencies, then we pass from the discrete to the continuous frequency, finding in this way Planck's energy distribution and Bose's entropy distribution.
The outline of the paper is as follows.In Section 2, we discuss the relationship between the multiplicity of energy eigenvalues and the density of states in the continuous frequency domain for photonic systems.Section 3 introduces the discrete spectrum of energy and entropy and shows that the passage to the continuous frequency leads to Planck's energy spectral density and Bose's entropy spectral density.Section 4 compares our derivation to other derivations.Section 5 discusses the difference between the entropy of the thermal state and that of a set of harmonic oscillators at the coherent state with the same energy distribution, showing that the difference between the two is really small.Finally, in Section 6 we draw the conclusions.

Multiplicity of Energy Eigenvalues and Density of States
Consider photons inside a cubic reflective cavity of side L. By imposing the boundary conditions, one finds by standard arguments that the energy eigenvalues are where the triple of natural numbers (i, j, k) identifies the eigenstate of the photon in the considered polarization, ν i,j,k is the eigenfrequency (i, j, k), c is the speed of the light and h is the universal constant that Planck determined in [1] to be equal to 6.55 • 10 −34 J• s, today h = 6.62 • 10 −34 J• s, by imposing compatibility between his theory and the experimental evidence.Let r i,j,k be the random occupancy number of eigenstate (i, j, k), i.e., the random number of photons whose eigenstate is (i, j, k).The energy and the frequency of the photon depend not on the specific triple (i, j, k) but only on the sum of their squares, that is, only on the following integer n: Let M n be the set of triples (i, j, k) that are distinct roots of (1) and let m n be the multiplicity, or degeneracy, of ϵ n , that is the number of triples in M n .For instance, The first values of m n are ), m 63 = 0.By Legendre's three squares theorem, all the natural numbers excepting those of the form n = 4 a (8b + 7), (a, b) ∈ N 2 , can be represented as the sum of three squares, therefore m 4 a (8b+7) = 0. Also, for convenience in the following we will put m 0 = 0.The relationship between the multiplicity of energy eigenvalues and the density of states is well known and widely discussed in the context of massive particles, see, e.g., [4].These arguments are based on the observation that the number of the points (eigenstates) with non-negative integer coordinates (i, j, k) that satisfy the inequality is equal to the number of points contained in the positive octant of a sphere of radius √ N. By quantizing the three-dimensional space into cubic cells of the unit side, hence of unit volume, we recognize that the number of points (eigenstates) that satisfy the above inequality is approximated to the volume of the positive octant of the sphere of radius √ N: Actually, since the right hand side of (3) is just the mentioned volume.By uniformly sampling the argument of the integral with unit step, we have the approximation which suggests that we should regard the argument of the sum in the right hand side as a smooth but still discrete version of the true multiplicity.For N → ∞, the relative error in the approximation (3) becomes vanishingly small, leading to The above limit is deeply discussed in 1.4 of [5] in the context of the analysis of multiplicity of energy eigenvalues of massive particles.
The arguments discussed till here can be used both in the analysis of multiplicity of energy eigevalues of systems of massive particles and of systems of photons.The difference between the two cases comes out when we consider Equation ( 2), which is specific for photons and leads us to define the continuous frequency ν as where V is the eigenfrequency that maps on the surface of the sphere.Looking at the integral in (3) as at a sum between n = 0 and n = N, we regard the integration variable s as the continuous version of index n, therefore, after the change of variables from s to ν in the integral (3), we obtain The fraction appearing inside the integral is called the density of states per polarization, which, multiplied by dν, is interpreted as the infinitesimal number of quantum states between frequency ν and frequency ν + dν per polarization.
The above density of states holds also for cavities with other geometries than the cube, provided that the volume of the cavity is equal to the volume L 3 of the cubic cavity, see [4] for the case of massive particles.For instance, when the cavity is a rectangular parallelepiped with sides L x , L y , L z , the eigenfrequecies are with V cavity = L x L y L z and If L x , L y , L z , are non-commensurable between them, then the multiplicity of all the eigenfrequencies is one and the number of eigenfrequencies such that is equal to the number of points with integer coordinates (i, j, k) contained in the first octant of the ellipsoid with semi-axes The mentioned number of points is approximated to the volume of the first octant of the ellipsoid, which, expressed in spherical coordinates is, where, in the last equality, we substitute ν = ρV.

Energy Spectral Density and Entropy Spectral Density of the Photonic Thermal State
Consider the thermal state and let r n be the random number of photons that occupy a quantum eigenstate whose three quantum numbers belong to M n .By standard maximization of Shannon entropy with constrained mean value, one finds that, whichever is the specific triple (i, j, k) ∈ M n , the probability distribution of the random occupancy number is geometric: where β = (k B T) −1 is the inverse temperature, k B = 1.38 • 10 −23 J/K is the Boltzmann constant, T is the temperature in Kelvin degrees and ϵ is the energy of the photon.The mean value of the random variable r n is the Bose statistics Probability distribution, variance and entropy can be expressed as a function of the mean value: where the entropy is expressed in k B units and the first term in the last line is β times the expectation of the sum of the energies of the random number of photons that populate the quantum state at hand.It is easy to see that which holds with equality for µ r n → 0.
Assuming independency between the occupancy numbers, the joint probability distribution of vector r of the occupancy numbers is Due to independency of the individual occupancy numbers, the total entropy of the considered polarization is the sum of the entropies of the occupancy numbers of the individual eigenstates where the first term in the second line is β times the energy distribution in the eigenfrequencies.Figure 1 reports the cumulative sum with the true multiplicity g n = m n and with the smooth multiplicity g n = π √ n/4, see the right hand side of (4).Defining the entropy spectral density as the continuous version of the entropy distribution with smooth multiplicity S r n π √ n/4, hence as the product between the density of states and the entropy of the geometric distribution expressed as a function of the continuous frequency, that is where, from (8), and integrating it, we arrive at Bose's entropy per polarization.Note that, due to an error, the sign of the exponent in the exponential in the second-last equation of [2] is flipped.
where the first addend inside the integral is β times Planck's energy spectral density per polarization.The calculus of the integral is standard.Specifically, it is based on

Comparison with Other Derivations
Bose's approach is semi-classical.He computes the infinitesimal number of quantum states that a photon with frequency between ν and ν + dν can occupy as the infinitesimal number of cells of volume h 3 contained in the region of the six-dimensional phase space of the photon where the sum of photon's three squared momenta is between (hν/c) 2 and (h(ν + dν)/c) 2 .The volume of the mentioned region in the three spatial coordinates is L 3 , while in the domain of the three momenta (p x , p y , p z ) the infinitesimal volume is that of a spherical shell between radia hν/c and h(ν + dν)/c, therefore the sought infinitesimal number of cells is where the first equality is based on the change of variables from Cartesian coordinates to spherical coordinates with radial coordinate hu/c: and ( 14) is just the argument of the integral (7).Pathria, in chapter 6 of [5], completely skips the discrete domain and basically follows the same approach of Bose, excepting that Pathria refers to (14) as to the Rayleigh expression for the (infinitesimal, we add) number of normal modes of vibration.
The standard modern approach of, among others [6][7][8][9][10], starts from the three quantum numbers (i, j, k) and derives the energy spectral density (the entropy spectral density can be derived in a completely similar manner) by considering a continuous triple (x, y, z) in place of (i, j, k): In the above sequence of passages, the last equality is obtained by changing the integration variables from Cartesian coordinates to spherical coordinates, that is, putting where the radial coordinate 2Lν/c is the continuous version of our √ n.Since now only non-negative values of (x, y, z) are allowed, the factor 1/8 in front of the integral restricts the domain of the integral to the positive octant of the three-dimensional space spanned by (x, y, z).A very similar derivation can be found in 52 of [11].An exception is Feynman, who derives the energy spectrum in an inherently continuous manner that is neither based on the discretization of phase space nor on the imposition of the boundary conditions, see 41-2 and 41-3 of [12].However, Feynman discusses only the expected energy, not the entire probability distribution, so he does not derive the entropy.The approach of Feynman seems to better fit the dense nature of the energy spectrum obtained from the measurements, while the approach based on the imposition of the boundary conditions leads to a spectral distribution of energy and the successive passage to the energy spectral density is justified only for L → ∞.In the case of box of finite size, even if the spectral distribution obtained by imposing the boundary conditions is not dense and the relative frequency spacing between the spectral lines with lower frequency is wide, it happens that the spectral density obtained after the passage from discrete spectrum to dense spectrum fits the experimental evidence, and this is the only reason that, to our opinion, can be invoked to support the approach based on the boundary conditions and on the discrete spectrum.
Exactly as it happens with Bose's derivation, the use of the triple integral (16) obscures the role of the multiplicity of energy eigenvalues, because the continuous variables (x, y, z) lead, through the successive change of continuous variables, to the continuous frequency, skipping at all the discrete frequency and, with it, the multiplicity of energy eigenvalues and the distribution of energy and entropy in the domain of the discrete frequency.
To further substantiate our claim of novelty, we point out that the density of states of massive particles is derived in many papers from the multiplicity of energy eigenvalues but the density of states of photons differs from the density of states of massive particles because, as it is apparent from (6) and as already discussed in one of the previous sections, the continuous version of the discrete energy of the photon, hence the continuous version of its discrete frequency, is proportional to √ s, while the energy of the massive particle is proportional to s.To our best knowledge, the passage from the discrete index n to the continuous frequency ν through (3) and through the successive change of variables ( 6) is proposed here for the first time.This passage allows us to derive energy spectral density and entropy spectral density from energy distribution and entropy distribution in the domain of the eigenfrequencies, hence from the multiplicity of energy eigenvalues.

Comparison with the Poisson Distribution of the Occupancy Numbers
The photonic occupancy number of the coherent state follows the Poisson distribution.One could figure out a collection of independent coherent states at the frequencies {ν n }.If the expected occupancy number is small at each one of the frequencies, then it is reasonable to expect that the resulting radiation is virtually non-coherent and that the statistics of the photonic occupancy numbers approaches that of the thermal state.
In the following we compare entropies of different probability distributions with the same energy spectral distribution.To this aim, we take for the probability P i,j,k that one photon occupies state (i, j, k) where The expected energy of the n-th eigenfrequency is that is the one already introduced in Equation ( 10).
It is worth observing that the Poisson distribution, besides being the distribution that characterizes the coherent state, is also in strict connection with the multinomial distribution, recently proposed in [13] for the occupancy numbers of the canonical state of the ideal gas.According to [13], the joint probability distribution of the occupancy numbers in the canonical ensemble is multinomial: where R is the total number of photons, which, in the canonical ensemble approach, is assumed to be fixed and known.However, in the case of photons inside a cavity, R is random, so we take a weighted average of the multinomial distributions with weights equal to the distribution of N: where, again, (20) is understood.The actual probability distribution of R is that of the random variable obtained from the sum of geometrically distributed random variables, see [14].Hower, this distribution is untractable.Assuming that the expectation of R is large, say, greater than 30, and since R is the sum of a large number of discrete random variables, we can approximate its distribution to a Poisson distribution with expected value µ R of the Poissonian random variable, see also [15] P Inserting the above distribution and (20) in (21), we see that the distribution of the occupancy numbers is the product of Poisson distributions: where, in the last equality, we substitute (19).In practice, the product of geometric distributions that characterize the thermal state becomes here the product of Poisson distributions with the same mean values.The entropy of the individual Poisson distribution with parameter µ r is equal to where E{•} is the expectation over the Poisson distribution of the function inside the curly brackets: When the expected number of photons in a quantum state tends to zero, the expectation in (23) tends to zero, leading to where the first inequality is consequence of the maximum entropy property of the geometric distribution, while the rightmost term is obtained by neglecting the non-negative expectation in the entropy of the Poisson distribution (23) and is equal to the right-hand side of (9). Figure 2 reports the entropy distribution in the domain of the eigenfrequency calculated with the smooth discrete multiplicity π √ n/4 for the geometric distribution, the Poisson distribution, and for the rightmost term of (25).
Using again the density of states and the rightmost term of (25), and using

Conclusions
In the paper we have reviewed the past derivations of Bose's entropy spectral density and presented a new one that is explicitly based on the multiplicity of the energy eigenvalues resulting from the imposition of the boundary conditions on the spatial boundaries of the blackbody.Among the approached mentioned in the paper, the author's favor goes to the one of Feynman, which seems to better reproduce the physics of the blackbody.This approach, discussed by Feynman only for energy, not for entropy, is not based on the eigenstates of the photons inside the box and leads directly to the energy spectral density, without any passage from the discrete domain of the eigenfrequencies to the domain of the continuous frequency.The derivation of the entropy spectral density from Feynman's approach is left to future investigation.
Funding: This research received no external funding.