_{1}

Electromagnetic (EM) radiation is both wave and heat. Waves are characterized by spectral distribution, spatial distribution, time coherence, spatial coherence, energy flux, and polarization. Heat, namely energy transferred from a hot body to a cold one, is characterized by its energy and entropy, and the ratio between them is the temperature. Here we calculate the entropy and temperature of a single radiation mode from the wave properties of the radiation. Using the Heisenberg uncertainty principle and Planck law , we calculate, from the optical properties of the radiation, the number of modes and their occupation number. Then we calculate the entropy and temperature of a single-mode EM radiation. It is shown that the entropy of a single-mode varies from zero for low occupation number, namely in the quantum limit to one Boltzmann constant for high occupation number, namely in the classical limit. The temperature varies from zero Kelvin in the quantum limit to infinity at high energies in the classical limit. This analysis is consistent with Fourier optics and statistical mechanics namely: Stephan-Boltzmann Law, Wein Law, Zipf law, IT File’s entropy, and the canonical distribution.

Thermodynamics is the science of energy flow. Most energy flow in nature is electromagnetic radiation. The electromagnetic radiation has different names for different wavelengths: The longest wavelengths (longer than one centimeter) are called radiowaves, in the millimetric range, microwaves in the microns range, infrared (that we feel as “heat”) in the range between around 0.8 to 0.3 micron, visible waves or light. In the less than microns, ultraviolet or UV, the shortest wavelength is the X-ray and γ rays that can ionize molecules. The energy of the electromagnetic radiation is moving at the speed of light [

Entropy is a measure of the uncertainty of a system. If a system can be found in W distinguishable configurations (microstates), its entropy is the logarithm of W multiplied by the Boltzmann constant k_{B}. Therefore, if we assume that each microstate has identical probability as the others, entropy is the measure of lack of knowledge in what microstate the system is. Nevertheless, entropy interpreted as uncertainty mostly in information theory, and in sparse physical systems, like gases, it is described erroneously as a disorder. The concept of uncertainty in physics is dedicated to the Heisenberg uncertainty principle that formulates the limit of the accuracy of simultaneous measurements of energy and time or momentum and space, i.e. if we want to measure the energy with very high accuracy, we need to give up the time resolution, etc. In fact, the Heisenberg uncertainty principle is a direct outcome of Fourier wave analysis [

δ v δ t ≥ 1 / 2 π and δ k x δ x ≥ 1 / 2 π (1)

where δ v is the spectral width, and δ t is the coherence time, δ k x is the width of the wavenumber distribution and δ x is the spatial coherence. This uncertainty bound 1/2π, depends on the profile of the beam, i.e. in a Gaussian beam the limit is 1/4π. These relations were known about 100 years before Heisenberg. Basically, Heisenberg following Planck invoked the Planck constant to connect between energy and frequency and momentum and wave number, namely ϵ = h v and p x = h k x where ϵ is the energy and p x is the momentum in the direction of the propagation axis. Therefore, for a wave described by Ψ = A e − i 2 π ( v t + k x x ) , Equation (1) yields the famous uncertainty principle,

δ ϵ δ t ≥ h / 2 π and δ p x δ x ≥ h / 2 π (2)

Here we address the question, what is the contribution of the Heisenberg uncertainty principle to the entropy, which is also the uncertainty of the energy emission of blackbody radiation.

A blackbody is any materialistic body in thermal equilibrium and therefore has a well-defined temperature T. The materialistic particles reach thermal equilibrium by exchanging electromagnetic radiation, namely photons, with each other. The blackbody’s volume V contains many radiation modes, N, in all the frequencies and directions that are possible in its geometric boundaries. The blackbody is linked to a thermal bath of temperature T such that its temperature remains constant despite the electromagnetic radiation that it emits. Till Planck’s publication in 1901 [_{B}T energy, where k_{B} is the Boltzmann constant. This value is the average energy of a classical oscillator in a bath at temperature T (i.e. each oscillation degree of freedom in ideal gas has internal energy k_{B}T). Therefore, to calculate the energy emitted by blackbody radiation we have to count the number of modes and multiply it by the average energy. The amount of heat Q emitted by the blackbody depends on the measuring time Δ t . Equation (1) yields that δ t ≥ 1 / 2 π δ v . However, we can increase the measuring time without changing the spectrum as long as we wish, namely Δ t > δ t = 1 / 2 π δ v . The ratio

Δ t δ t = N l , which is the number of the coherence length in the pulse, is called the number of longitudinal

modes.

The number of the possible oscillators per unit volume V, in a frequency interval v + δ v , is derived from the number of radiation modes N calculated from Maxwell equation [

N = 8 π v 2 δ v V c 3 (3)

The fraction of the modes contained within a solid angle Ω will be,

N ( Ω ) = N 4 π = 2 v 2 c 3 V δ v (4)

Now assume that a beam of duration Δ t emerges from a square on the blackbody surface having a cross-section a 2 . Its volume is V = a 2 c Δ t and the number of modes that participate in its emission is,

N ( Ω ) = 2 v 2 c 2 δ v a 2 Δ t = 2 a 2 λ 2 δ v Δ t = 2 N l ( a 2 λ 2 ) (5)

Now we go back to Equation (2) δ p x δ x ≥ h / 2 π and divide both sides of Equation (2) by p x . Since the wave vector k x = 2 π λ , where λ is the wavelength, we obtain that ( δ p x p x ) δ x ≥ λ 2 π , δ p x p x = δ Ω x is

the diffraction limit angular divergence in the x direction and δ x is the size of the spot of the focused beam in the x direction, therefore, δ x > λ / 2 π and the minimum area resolution of the spot size is

δ x δ y > λ 2 / ( 2 π ) 2 . With analogy to the definition of longitudinal modes, the ratio between the spot size resolution a 2 / ( 2 π ) 2 and the diffraction limit resolution, namely a 2 λ 2 is called the number of spatial

modes N s . We see that a spatially coherent beam, in which N s x = 1 , cannot be resolved better than λ / 2 π . Equation (5) now can be written as,

N = 2 N l N s (6)

The number 2 stands for the two independent polarizations, N l is the number of the coherence length of a beam and N s is our ability to focus the beam as compared to the diffraction limit spot size [

N is a function of the temporal coherence and the spatial coherence and the polarization. The heat emitted from a blackbody is thus,

Q = N k B T (7)

The entropy reduction of a Blackbody due to energy emission is,

S = Q T = N k B (8)

Equations (7-8) are the expressions of the heat and entropy of N classical modes of electromagnetic

radiation. A single classical radiation mode has an energy q = Q N = k B T and an entropy s = S N = k B . This

well-known result is sometimes overlooked and it is assumed by some, that a single-mode pulse has zero entropy. The reason for this somewhat common mistake is that the expression of the entropy in statistical mechanics is proportional to the logarithm of the number of microstates. Since a single oscillator may be considered as one microstate, and since the logarithm of one is zero, one may conclude that a solitary event, well defined in space and time, has zero entropy. However, single-mode radiation, because of the Heisenberg uncertainty principle, has a built-in uncertainty (entropy) in the oscillation itself, namely the potential and the kinetic energy and the spatial resolution, therefore even a single classical oscillator has k_{B} entropy. Moreover, if a single oscillator would have zero entropy since entropy is extensive, it would lead to the conclusion that any combination of single oscillators has zero entropy, including the blackbody radiation that is a plurality of oscillators in different modes. This conclusion is, of course, a violation of the second law. A simple proof that electromagnetic radiation should carry entropy can be demonstrated by building a heat machine which consists of a hot bath that is a blackbody at temperature T_{H} and a cold bath that is also a blackbody with a temperature T_{C}. We replace the cylinder with a piston that contains gas by an ideal fast shutter that can transmit or reflect the electromagnetic emitted radiation. Now we open the shutter for a very short time such that the energy of the electromagnetic radiation leaves the hot blackbody but not yet be absorbed by the cold blackbody. In the short time interval that the electromagnetic radiation is between the two black bodies, the entropy of the hot blackbody decreases but the entropy of the cold blackbody is not yet increased, therefore, if electromagnetic radiation does not carry, entropy, the second law is violated. The only explanation for this “paradox” is that the electromagnetic radiation must carry entropy equal to that removed from the hot blackbody.

Equation (7) is known as Rayleigh-Jeans approximation. The problem in Equation (7) is that when λ is very small, the number of the spatial modes increases to infinity and the blackbody should then emit an infinite amount of energy which is, of course, not possible. This problem was known as the “UV catastrophe”. In 1901 Planck replaced the frequency-independent average energy of a mode q = k B T in Equation (2) with the famous frequency-dependent formula,

q = n h v = h v exp ( h v k B T ) − 1 (9)

and added a new universal constant h to physics.

The famous Planck’s assumption is that energy is not continuous and can be divided into quanta hv that are a linear function of their frequencies. Therefore, a spatial radiation mode having a very small wavelength has a very high frequency and therefore its quantum, the photon, may have energy that is higher than k_{B}T. Therefore, according to classical mechanics, it cannot emit energy.

In quantum mechanics, the emission of photons with energy higher than the average energy k_{B}T is possible. The reason for this strange phenomenon is the uncertainty expressed by the entropy. The second law is about the propensity of entropy to grow. In a closed system, it will reach its maximum value where it cannot grow anymore. This point is the equilibrium point. Planck was the first to calculate the distribution of photons among the radiation modes that maximized the Boltzmann entropy S = k B ln W under the constraint of conservation of energy [

W = ( N + B − 1 ) ! ( N − 1 ) ! B ! , applying Stirling’s approximation to calculate ln W he obtained that,

S ≅ k B N [ ( n + 1 ) ln ( n + 1 ) − n ln n ] (10)

Here Planck replaced the total number of particles B by n = B / N which is the average number of particles in a box. To maximize the entropy, he applied the conservation of energy constraint

Q ( n ) = N q = N n h v , where, q is the energy per mode, and therefore, to find the equilibrium distribution, he used Lagrange multiplier technique, namely to find the maximum of the function F = S ( n ) − β Q (n)

with respect to n, where β is a Lagrange multiplier, namely ∂ F ∂ n = 0 , where,

F ≅ k B N [ ( n + 1 ) ln ( n + 1 ) − n ln n − β ( q − n h v ) ] (11)

Since it can be shown that β = 1 k B T , he found out his famous result for n that maximizes the entropy namely,

n ( v , T ) = 1 e h v k B T − 1 (12)

Now the energy carried by a radiation mode is not k B T but n ( v , T ) h v as given by Equation (9).

In the classical limit, the energy is continuous, the photon’s energy is infinitely small, and therefore the number of photons is infinitely large. We rewrite Equation (12) as

h v = k B T ln ( 1 + 1 n ) (13)

Since the energy of the mode is q = n h v , in the classical limit where n → ∞

q = n h v = k B T lim n → ∞ ln ( 1 + 1 n ) n = k B T

And the entropy of a classical mode is s = k B , namely, one Boltzmann constant as was assumed in the Rayleigh-Jeans formula.

It is worth noting that both the entropy and the frequency of a classical oscillator are independent of its energy. The energy-independent frequency of the classical oscillators is used for time measurements i.e. pendulum clocks and watches. The energy-independent entropy of classical oscillators is used for digital communication [

In the classical limit where the energy of the photon is very small, e h v k B T ≈ 1 + h v k B T , therefore, from Equation (12) we can see that n ≈ k B T h v . Namely, the number of photons is the average energy of a mode

divided by the energy of the photon. That means that the higher the energy of photons the higher their frequency and smaller their number, n. Namely, there are more poor-energy photons than high-energy photons.

To find the exact ratio between the frequency v and the number of photons n in N modes, we normalize the photon energy of Equation (13), by dividing it in the total sum of the energies over the modes. Namely,

∑ n = 1 N h v ( n , T ) = k B T ∑ n = 1 N ln ( 1 + 1 n ) = k B T ln ( N + 1 )

The normalized distribution is given by ε ( n , T ) = h v ( n , T ) / ∑ n = 1 N h v ( n , T ) or:

ε ( n , N ) = ln ( 1 + 1 n ) ln ( N + 1 ) = log N + 1 ( 1 + 1 n ) (14)

ε ( n ) , which is a parameter-free long-tail distribution, called the Planck-Benford distribution [7 - 9], is the relative energy of the modes having n photons as compared to the energy of all the N modes. It is seen that the higher the photon’s energy, the smaller the number of photons having this energy.

The reason for the surprising result, that the Planck-Benford distribution has no physical constants, is that in the classical limit the radiation modes are pure harmonic oscillators and therefore have the universal statistics of harmonic oscillators.

This outcome that in the energy normalization process all “physics” disappeared, namely: h, k_{B} and T are canceled out, makes the Planck-Benford distribution more universal than the Maxwell-Boltzmann distribution. The Planck-Benford law is also the Maximum Entropy distribution of identical balls in N distinguishable boxes, where the number of balls is greater than the number of boxes [

Here we show that Planck-Benford law is the empirical Zipf law. Zipf law is a probability function found in many distributions in nature i.e. n people in N cities or the redundancy n of words in a text of N different words. Zipf law was found, empirically, to be, φ ( n , N ) = 1 / n s H N , where H_{N} is the harmonic number and s is some parameter (usually close to one). To show that Zipf law is Planck-Benford law with s = 1 , we assume that n is a continuous function, therefore, we can use the Riemann sum to find the

area between two sequential integers n and n + 1 , namely, ∫ n n + 1 d n ′ n ′ = ln ( n + 1 ) − ln ( n ) = ln ( 1 + 1 n ) . Similarly, the Harmonic number defined by H N = ∑ n = 1 N 1 n becomes,

H N + 1 = ∫ 1 N + 1 d n ′ n ′ = ln ( N + 1 ) .

Note that the integration is until N + 1 because we start to count the N ranks from n = 1 and not from zero (from zero to one we are in the quantum regime). Now we can write Planck-Benford law as,

ε ( n , N ) = ln ( 1 + 1 n ) ln ( N + 1 ) = 1 n ∑ n = 1 N 1 n = 1 n H N (15)

This is the Zipf law with slope 1.

In the quantum limit, the photon energy has approximately 1/n times more energy than the average energy k_{B}T. One may expect that the probability to emit photons will be also 1/n. However, from Equation (13) we see that in the quantum limit where n < 1 , the frequency of a photon increases linearly with the temperature of the blackbody and logarithmically with the number of photons n. For very small n’s,

1 + 1 n ≈ 1 n and h v ( n , T ) = k B T ln ( 1 + 1 n ) ≈ − k B T ln n . This equation usually is written as n i ≈ exp ( − h v i k B T ) . When it is normalized, by dividing it by the probability partition function, Z P ( n , T ) = ∑ i n i , we obtain

the famous exponential “canonical energy distribution”.

p i = exp ( − h v i k B T ) ∑ i exp ( − h v i k B T ) (16)

The canonical distribution derivation deserves some remarks. Zipf law was obtained by energy normalization. However, the energy is inversely proportional to the number of photons in a mode that is proportional to probability. The canonical distribution was obtained directly by probability normalization

in which n was changed to p. The energy emission of a single-mode in this limit is, ε i = h v i exp ( − h v i k B T ) .

Stephan-Boltzmann law enables us to calculate the energy flux of a blackbody as a function of its temperature [

q ( v , T ) = n h v = h v e h v k B T − 1 = k B T h v k B T e h v k B T − 1 (17)

And the entropy of a single mode of frequency v is,

s ( v , T ) = q ( v , T ) T = k B h v k B T e h v k B T − 1 (18)

Let’s designate u = h v k B T , therefore, the expression for the energy of a single-mode is q ( v , T ) = k B T u e u − 1 , and its entropy is s ( v , T ) = k B u e u − 1 . See

To calculate the energy emission of a blackbody we have to integrate Equation (17) over all the temporal modes. Let us start with the emission of a single spatial mode. The energy flux I ( v , T ) is q ( v , T ) , is divided by the area of the spot size of a single spatial mode λ 2 on the surface of the blackbody and multiplied by the two polarizations. Thus the energy flux is,

I E ( T ) = ∫ 0 ∞ ( 2 q ( v , T ) λ 2 ) d v (19)

Or,

I E ( T ) = 2 k B ∫ 0 ∞ T u λ 2 ( e u − 1 ) d v

Since d u = h k B T d v and 1 λ 2 = v 2 c 2 , we obtain,

I E ( T ) = 2 k B ∫ 0 ∞ T u v 2 c 2 ( e u − 1 ) d v = 2 k B 4 ∫ 0 ∞ T 4 u 3 c 2 h 3 ( e u − 1 ) d u

The finite integral ∫ 0 ∞ u 3 e u − 1 d u = ζ ( 4 ) = π 4 15 where ζ is Reimann zeta function, therefore, the energy flux of a single spatial mode is,

I E ( T ) = 2 ( π k B T ) 4 15 c 2 h 3 = σ T 4 (20)

where I ( T ) has units of energy per unit area per unit time.

The entropy flux is,

I s ( T ) = σ T 3 (21)

One should remember that the area of a single-mode λ 2 that is part of the integration, vary from zero to the significant part of the blackbody surface. This point deserves careful consideration. The

function f ( u ) = u e u − 1 = s k B is monotonically decreasing function as is seen in

The reason for this paradox is that quantum modes have very short wavelengths and the energy flux

increases by the factor 1 / λ 2 as was derived in Equation (19). The function that obtained f ′ = u 3 e u − 1 has a maximum.

Stephan-Boltzmann law was derived in the Ph. D. thesis of Boltzmann under the supervision of Stephan. The derivation was done by classical thermodynamics argumentations and therefore the value of the Stephan-Boltzmann constant was not calculated. The constant that was calculated later was an additional validation of Planck’s theory.

Another quantity that derives from Planck’s law is the maximum frequency of a blackbody’s emission as a function of its temperature [

in 1901. Here it is derived from the Planck law. We start with q ( v , T ) = k B T u e u − 1 remembering that u = h v / k B T , the energy per mode per unit blackbody surface is q s ( v , T ) = 2 q ( v , T ) / λ 2 = 2 k B T c 2 v 2 u e u − 1 = 2 h c 2 v 3 e u − 1 . To find the maximum frequency on the energy distribution we ask that d q s ( v , T ) d v = 0 or,

d q s ( v , T ) d v = ( 2 h c 2 ) d d v v 3 e u − 1 = 2 h c 2 ( e u − 1 ) ( 3 v 2 − v 2 u e u e u − 1 ) = 0

Namely,

u e u e u − 1 = 3 (22)

The function u e u e u − 1 has poles; nevertheless, we can calculate it in the limits. In the classical limit u ≪ 1 , so then we obtain that e u max ≈ 3 and therefore h v max k B T = ln 3 .

In the quantum limit u ≫ 1 e u ≈ e u − 1 , therefore, h v max k B T = 3 .

Substitute v max λ max = c and obtain,

λ max = b T (23)

In the classical limit b c l a s s = h c k B ln 3 and in the quantum limit b q u n t = h c 3 k B .

Usually, Wein’s law applied to Maxwell-Boltzmann distribution, which is in the quantum regime. Therefore, Equation (23) converges very fast to the area where b is constant.

Temperature is considered a macroscopic property. However, a microscopic entity which its energy

can be measured and its entropy can be calculated from Clausius entropy definition s = q T , has a temperature. One should ask, what is the physical meaning of it?

In the classical approximation, the answer is straightforward: energy flows from the hot to the cold. I.e. two identical classical pendulums, one hot with high amplitude and one cold with low amplitude, if they linked together, their amplitudes will eventually equate. The laser beam is a classical oscillator as n ≫ 1 , its beam temperature is the maximum temperature that it can heat a materialistic object. Consider a laser beam having power P and coherence length δ x , therefore, its coherence time is δ t = δ x / c . The

energy per mode is the energy within the coherence time. Therefore, q = P δ t = P δ x c . The temperature of the beam is,

T c l a s s = q k B = P δ x c k B = n h v k B (24)

I.e. HeNe laser having a power of 1 mW and coherence length 30 cm yields a temperature of about 10^{10} K and the energy of a mode is about 10^{−13} J (1 mW∙0.3/c). This is a very high temperature, but not infinite as it is usually assumed. Nevertheless, for any practical reason, this temperature is equivalent to infinite temperature. The reason that the HeNe laser beam does not fit for heating is its negligible energy per mode 10^{−13} J. Namely, to heat 1 gram of water in one degree Kelvin one needs about 4 × 10^{13} modes, which means, of about 4 hours of heating.

Digital communication is done by a transmission of a file that is a sequence of classical oscillators from Bob to Alice. Therefore, the entropy of the file which reflects its amount of data in it is not a function of the temperature of the bits. Nevertheless, we want that it will be well above the ambient’s temperature such that the thermal noise will not add noise bits the file. Therefore the power namely, the temperature of our antenna has practical importance.

The general expression of the temperature of EM radiation is given by Planck expression of Equation (12) for the occupation number n, that can be rewritten as,

T = h v k B ln ( 1 + 1 n ) (25)

In

It is seen that when n approach to zero both the entropy and the temperature are approaching zero. Entropy is heat divided by the temperature, therefore we might expect that zero temperature will yield infinite entropy, like in a zero temperature bath. The explanation to this “abnormality”, namely that

lim n → 0 q T = 0 , is that the energy reduces exponentially to zero, while the temperature decay linearly. From Equation (17) we obtain that in the quantum limit q = n h v = h v e h v k B T − 1 ≈ h v e − h v k B T . The temperature of Equation (25), when we multiply by n both nominator and the denominator, yields T = n h v k B ln ( 1 + 1 n ) n . Therefore in the quantum limit lim n → 0 T = lim n → 0 n h v k B the temperature converges to zero linearly with n while

the heat converges exponentially so that both temperature and entropy approach to zero when n approach to zero.

How to calculate this temperature? Suppose that we want to measure the temperature of a single photon, a single spatial mode source. We need a sensitive detector that can detect a single photon. Suppose that we measure on the average, P single photons over a time period Δ t . The detector is equipped with a filter that passes the frequency v, and its bandwidth is δ v . The number of modes N = 2 π Δ t δ v . The

occupation number is n = B N . Photons with a narrower bandwidth δ v can also pass the filter, therefore, N ≤ 2 π Δ t δ v and n ≥ B 2 π Δ t δ v .

Can we transmit data with single-photon files? Suppose that we have a file composed of HeNe single-photons bits. We ask that its temperature will be at least ten times higher than that of the ambient, to

reduce thermal noise, namely T = 2.28 × 10 4 ln ( 1 + 1 n ) ≥ 10 T a m b i e n t . At room temperature T = 300 K , thus we obtain ln ( 1 + 1 n ) ≥ 7.6 or n ≥ 1 2000 , that is an acceptable limit for communication.

EM radiation is heat. In thermodynamics, heat is characterized by two independent quantities, energy and entropy. EM radiation has many more characteristics namely; frequency, spectral distribution, temporal distribution, temporal coherence, spatial coherence, beam profile, polarization, and direction. How all these properties are connected to that of heat in general? Planck’s law connects between the numbers of photons in a mode (occupation number) and their frequency in a thermal bath (blackbody) of a given temperature. Here we use the spectral distribution, temporal distribution, temporal coherence, spatial

coherence, beam profile, and polarization to calculate the number of modes. From the EM beam’s power, we calculate the number of photons. This enables us to calculate for EM radiation the entropy and the energy of a mode that is consistent with Planck’s law.

It is shown that the EM radiation statistics vary smoothly from the classical limit in which all the radiation modes in equilibrium have an energy of k_{B}T like in ideal gas oscillators and the quantum limit in which most of the modes are empty and some have very high energy as compared to k_{B}T. Using Planck law, we review in a consistent way the calculation of all the thermodynamical properties, namely classical energy distribution—Zipf law; quantum energy distribution—Maxwell Boltzmann law; Energy flux—Stephan-Boltzmann law, and maximum frequency—Wein displacement law.

The entropy of a single-mode as given by Equation (18) namely s ( v , T ) = q ( v , T ) T = k B u exp ( u ) − 1 where u = h v k B T , monotonically increases from zero at very low occupation number in which u is very

large namely, h v ≫ k B T to k B when the occupation number is going to infinity namely h v ≪ k B T and u → 0 . From Equation (25) it is seen that for low occupation number n → 0 temperature increases monotonically from zero in low occupation number to infinity in a very large occupation number.

The author declares no conflicts of interest regarding the publication of this paper.

δ v : spectral width

δ t : coherence time

δ k x : wavenumber’s width

δ x : spatial coherence

ϵ : energy

h: Planck constant

v: frequency

λ : wavelength

c : speed of light

T: temperature

V: volume

N: radiation modes

N_{l}: longitudinal modes

N_{s}: Spatial modes

Ω : Solid angle

a: beam waist

k_{B}: Boltzmann constant

Q: Heat

q: heat per mode

Δ t : measuring time

S: entropy

s: entropy per mode

n: occupation number – the number of photons in a mode

B: total number of photons

W: number of microstates

β : Lagrange multiplier

ε ( n ) : normalized distribution of frequency

H_{N}: Harmonic number

Z_{P}: canonical partition function

u: the energy of a photon relative to the average energy.

I: energy flux

v max : The frequency of maximum emission of a blackbody.

λ max : The wavelength of maximum emission of a blackbody

b: Stephan Boltzmann constant

P: power